title: Structured Shrinkage Priors

abstract: In many regression settings the unknown coefficients may have some known structure, for instance they may be ordered in space or correspond to a vectorized matrix or tensor. At the same time, the unknown coefficients may be sparse, with many nearly or exactly equal to zero. However, many commonly used priors and corresponding penalties for coefficients do not encourage simultaneously structured and sparse estimates. In this paper we develop structured shrinkage priors that generalize multivariate normal, Laplace, exponential power and normal-gamma priors. These priors allow the regression coefficients to be correlated a priori without sacrificing elementwise sparsity or shrinkage. The primary challenges in working with these structured shrinkage priors are computational, as the corresponding penalties are intractable integrals and the full conditional distributions that are needed to approximate the posterior mode or simulate from the posterior distribution may be non-standard. We overcome these issues using a flexible elliptical slice sampling procedure, and demonstrate that these priors can be used to introduce structure while preserving sparsity.


title: A framework for synchronizing a team of aerial robots in communication-limited environments

abstract: This paper addresses a synchronization problem that arises when a team of aerial robots (ARs) need to communicate while performing assigned tasks in a cooperative scenario. Each robot has a limited communication range and flies within a previously assigned closed trajectory. When two robots are close enough, a communication link may be established, allowing the robots to exchange information. The goal is to schedule the flights such that the entire system can be synchronized for maximum information exchange, that is, every pair of neighbors always visit the feasible communication link at the same time. We propose an algorithm for scheduling a team of robots in this scenario and propose a robust framework in which the synchronization of a large team of robots is assured. The approach allows us to design a fault-tolerant system that can be used for multiple tasks such as surveillance, area exploration, searching for targets in a hazardous environment, and assembly and structure construction, to name a few.


title: Realism and causality I: Pilot wave and retrocausal models as possible facilitators

abstract: Of all basic principles of classical physics, realism should arguably be the last to be given up when seeking a better interpretation of quantum mechanics. We examine the de Broglie-Bohm pilot wave theory as a well developed example of a realistic theory. We present three challenges to a naive reading of pilot-wave theory, each based on a system of several entangled particles. With the help of a coarse graining of pilot wave theory into a discrete system, we show how these challenges can be answered. However this comes with a cost. In the description of individual systems, particles appear to scatter off empty branches of the wave function as if they were particles, and conversely travel through particles as if they were waves. More generally, the "particles" of pilot wave theory are led by the guidance equation to move in ways no classical particle would, involving apparent violations of the principles of inertia and momentum conservation.We next argue that the aforementioned cost can be avoided within a retrocausal model. In the proposed version of the pilot wave theory, the particle is guided by a combination of advanced and retarded waves. The resulting account for quantum physics seems to have greater heuristic power, it demands less damage to intuition, and moreover provides some general hints regarding spacetime and causality. This is the first of two papers. In the second [1] we show that, in the context of an explicit model, retrocausality, with respect to an effective, emergent spacetime metric, can coexist with a strict irreversibility of causal processes.


title: SUSY Signals from QCD Production at the Upgraded LHC

abstract: Weak-scale supersymmetry remains to be one of the best-motivated theories of physics beyond the Standard Model. We evaluate the sensitivities of the High Luminosity (HL) and High Energy (HE) upgrades of the LHC to gluinos and stops, decaying through the simplified topologies $\tilde{g} \to q \bar{q} \chi^0$, $\tilde{g} \to t \bar{t} \chi^0$ and $\tilde{t} \to t \tilde{\chi}^0$. Our HL-LHC analyses improve on existing experimental projections by optimizing the acceptance of kinematic variables. The HE-LHC studies represent the first 27 TeV analyses. We find that the HL-(HE-)LHC with 3 ab$^{-1}$ (15 ab$^{-1}$) of integrated luminosity will be sensitive to the masses of gluinos and stops at 3.2 (5.7) TeV and 1.5 (2.7) TeV, respectively, decaying to massless neutralinos.


title: Monge-Amp\`ere of Pac-Man

abstract: We show that the Monge-Amp\`ere density of the extremal function $V_P$ for a non-convex Pac-Man set $P\subset {\bf R}^2$ tends to a finite limit as we approach the vertex $p$ of $P$ linearly but with a value that may vary with the line. On the other hand, along a tangential approach to $p$ the Monge-Amp\`ere density becomes unbounded. This partially mimics the behavior of the Monge-Amp\`ere density of the union of two quarter disks set $S$ of Sigurdsson and Snaebjarnarson. We also recover their formula for $V_S$ by elementary methods.


title: Graphene thermal break-down induced by anharmonic bending mode

abstract: The abrupt loss of mechanical stability of two-dimensional graphene-type crystals at a certain transition temperature is described. At this temperature, the graphene state with practically zero-speed bending sound and developed bending fluctuations becomes energetically favorable. Such phenomenon, akin to melting, is naturally caused by the anharmonicity of crystal oscillations. In order to circumvent the known difficulties associated with taking into account the anharmonic effects, we propose an original pseudo-harmonic approximation, within which we determine the free energy of the anharmonic crystal and find a numerical characteristic for the intensity of bending vibrations at transition temperature. This characteristic is similar to the empiric Lindemann criterion for the melting phenomenon. At the same time, in contrast to the conventional Lindemann criterion, the found characteristic is explicitly expressed through the graphene bending moduli of the second, third, and fourth orders.


title: From Time-Domain Data to Low-Dimensional Structured Models

abstract: We present a framework for constructing a structured realization of a linear time-invariant dynamical system solely from a discrete sampling of an input and output trajectory of the system. We estimate the transfer function of the original model at selected frequencies using a modification of the empirical transfer function estimation that was recently presented in [Peherstorfer, Gugercin, Willcox, SIAM J. Sci. Comput., 39(5):2152-2178, 2017]. Our realization interpolates the transfer function estimates and can be seen as a generalization of the Loewner framework to structured systems. We demonstrate the presented framework by means of a delay example.


title: A Study on Graph-Structured Recurrent Neural Networks and Sparsification with Application to Epidemic Forecasting

abstract: We study epidemic forecasting on real-world health data by a graph-structured recurrent neural network (GSRNN). We achieve state-of-the-art forecasting accuracy on the benchmark CDC dataset. To improve model efficiency, we sparsify the network weights via transformed-$\ell_1$ penalty and maintain prediction accuracy at the same level with 70% of the network weights being zero.


title: Quark Gluon Plasma Formation in Proton-Proton Collisions Using PYTHIA

abstract: The aim of this study is to investigate the possibility of the Quark Gluon Plasma (QGP) formation in proton-proton (p-p) collisions. The p-p collisions at center of mass energies $\sqrt{s_{_{NN}}} = 200$ $GeV$ and $\sqrt{s_{_{NN}}} = 13$ $TeV$ were simulated using PYTHIA. These are the available energies of the current collider experiments; the Relativistic Heavy Ion Collider (RHIC), and the Large Hadron Collider (LHC). To search for the QGP Phase, the near-side yield to the away-side yield ratios were calculated for high transverse momentum particles. These ratios were then compared for the low and high multiplicity events as a function of transverse momentum. At the LHC energy in the high multiplicity events, the away-side yields show suppression compared with the near-side yields. This suppression might indicate the possibility of forming the QGP at the high-energy and high-multiplicity events in p-p collisions.


title: An extension of the fewest switches surface hopping algorithm to complex Hamiltonians and photophysics in magnetic fields: Berry's phase and "magnetic" forces

abstract: We present a preliminary extension of the fewest switches surface hopping (FSSH) algorithm to the case of complex Hamiltonians as appropriate for modeling the dynamics of photoexcited molecules in magnetic fields. We make ansatze for the direction of momentum rescaling and we account for Berry's phase effects through "magnetic" forces as applicable in the adiabatic limit. Because Berry's phase is a nonlocal, topological characteristic of a set of entangled potential energy surfaces, we find that Tully's local FSSH algorithm can only partially capture the correct physics.


title: Probabilistic Neural Architecture Search

abstract: In neural architecture search (NAS), the space of neural network architectures is automatically explored to maximize predictive accuracy for a given task. Despite the success of recent approaches, most existing methods cannot be directly applied to large scale problems because of their prohibitive computational complexity or high memory usage. In this work, we propose a Probabilistic approach to neural ARchitecture SEarCh (PARSEC) that drastically reduces memory requirements while maintaining state-of-the-art computational complexity, making it possible to directly search over more complex architectures and larger datasets. Our approach only requires as much memory as is needed to train a single architecture from our search space. This is due to a memory-efficient sampling procedure wherein we learn a probability distribution over high-performing neural network architectures. Importantly, this framework enables us to transfer the distribution of architectures learnt on smaller problems to larger ones, further reducing the computational cost. We showcase the advantages of our approach in applications to CIFAR-10 and ImageNet, where our approach outperforms methods with double its computational cost and matches the performance of methods with costs that are three orders of magnitude larger.


title: Structure theorems for actions of homeomorphism groups

abstract: We give general classification and structure theorems for actions of groups of homeomorphisms and diffeomorphisms on manifolds, reminiscent of classical results for actions of (locally) compact groups. This gives a negative answer to Ghys' "extension problem" for diffeomorphisms of manifolds with boundary, as well as a classification of all homomorphisma $\mathrm{Homeo}_0(M) \to \mathrm{Homeo}_0(N)$ when dim(M) = dim(N) (and related results for diffeomorphisms), and a complete classification of actions of $\mathrm{Homeo}_0(S^1)$ on surfaces. This resolves many problems in a program initiated by Ghys, and gives definitive answers to conjectures of Militon and Hurtado and a question of Rubin.


title: Unitarity of Singh-Hagen model in $D$ dimensions

abstract: The particle content of the Singh-Hagen model ($SH$) in $D$ dimensions is revisited. We suggest a complete set of spin-projection operators acting on totally symmetric rank-3 fields. We give a general expression for the propagator and determine the coefficients of the $SH$ model confirming previous results of the literature. Adding totally symmetric source terms we provide an unitarity analysis in $D$ dimensions.


title: Overdetermined Systems of Equations on Toric, Spherical, and Other Algebraic Varieties

abstract: Let $E_1,\ldots,E_k$ be a collection of linear series on an algebraic variety $X$ over $\mathbb{C}$. That is, $E_i\subset H^0(X, \mathcal{L}_i)$ is a finite dimensional subspace of the space of regular sections of line bundles $ \mathcal{L}_i$. Such a collection is called overdetermined if the generic system \[ s_1 = \ldots = s_k = 0, \] with $s_i\in E_i$ does not have any roots on $X$. In this paper we study solvable systems which are given by an overdetermined collection of linear series. Generalizing the notion of a resultant hypersurface we define a consistency variety $R\subset \prod_{i=1}^k E_i$ as the closure of the set of all systems which have at least one common root and study general properties of zero sets $Z_{\bf s}$ of a generic consistent system ${\bf s}\in R$. Then, in the case of equivariant linear series on spherical homogeneous spaces we provide a strategy for computing discrete invariants of such generic non-empty set $Z_{\bf s}$. For equivariant linear series on the torus $(\mathbb{C}^*)^n$ this strategy provides explicit calculations and generalizes the theory of Newton polyhedra.


title: Design of Face Centered Cubic Co81.8Si9.1B9.1 with High Magnetocrystalline Anisotropy

abstract: Despite the composition close to glassy forming alloys, face centered cubic (FCC) Co81.8Si9.1B9.1, designed based on Co9B atomic cluster (polyhedral), are synthesized as singlephase ribbons successfully. These ribbons, with grain sizes of ca. 92 nm, show supreme ductility and strong orientation along (111), which couples with shape anisotropy leading to high magnetocrystalline anisotropy comparable to Co rich Co-Pt nanoscale thin films, with a coercivity of 430 Oe and squareness of 0.82 at room temperature. The stability and magnetic behaviors of the phase are discussed based on experimental electronic structure. This work not only develops low cost Co-based materials for hard magnetic applications, but also extends the atomic cluster model developed for amorphous alloys into the design of new crystalline materials.


title: Loop interactions and their representations in Fock space

abstract: Given a weighted graph, we show how to define two types of natural interactions which correspond to local interactions between two Fock spaces. The first type of interaction involves loop ensembles and spanning trees. The second type of interaction involves loop holonomies and random connections.


title: Active Suppression of Microphonics Detuning in high $Q_L$ Cavities

abstract: Operation of Superconducting Radio Frequency (SRF) cavities with high loaded quality factors is becoming increasingly preferred for applications which involve low beam loading including Energy Recovery Linacs (ERL). Vibration induced microphonics detuning poses a major operational bottleneck in these low bandwidth systems, adversely affecting field stability. Besides passive measures of mitigating the vibration sources, modern SRF cavities are also attached to fast tuners incorporating piezo-electric actuators. We demonstrate the narrow band active noise control algorithm for realizing active resonance control and propose a modification based on the Least Mean Square approach to adaptively tune the control parameters and study it's stability and performance. We discuss our experience of using passive mitigation techniques while commissioning the Main Linac Cryomodule of the Cornell-BNL ERL Test Accelerator (CBETA) and report a net reduction in peak detuning by more than a factor of 2 in its unstiffened cavities. Finally, we demonstrate stable performance of our resonance control system with consistent reduction of peak microphonics detuning by almost a factor of 2 on multiple cavities.


title: Helical Miura Origami

abstract: We characterize the phase-space of all Helical Miura Origami. These structures are obtained by taking a partially folded Miura parallelogram as the unit cell, applying a generic helical or rod group to the cell, and characterizing all the parameters that lead to a globally compatible origami structure. When such compatibility is achieved, the result is cylindrical-type origami that can be manufactured from a suitably designed flat tessellation and "rolled-up" by a rigidly foldable motion into a cylinder. We find that the closed Helical Miura Origami are generically rigid to deformations that preserve cylindrical symmetry, but multistable. We are inspired by the ways atomic structures deform [1] to develop two broad strategies for reconfigurability: motion by slip, which involves relaxing the closure condition; and motion by phase transformation, which exploits multistability. Taken together, these results provide a comprehensive description of the phase-space of cylindrical origami, as well as quantitative design guidance for their use as actuators or metamaterials that exploit twist, axial extension, radial expansion, and symmetry.


title: Adsorbate-enhanced field-emission from single-walled carbon nanotubes: a comparative first-principles study

abstract: We report a comparative ab-initio study of different types of adsorbates and their adsorption mechanism on the field-emission performance of single-walled carbon nanotubes by analyzing electrical properties and transport characteristics and considering the thermal stability of the adsorbed structure. Adsorbates were found to reduce the work function by up to 1.3 eV, enhance tunneling near the carbon nanotube tip, and increase the field-emission current by as much as two orders of magnitude. A significant localization of the electron cloud was also observed near the adsorbates under a high applied electric field .


title: Force sensing in hybrid Bose-Einstein condensate optomechanics based on parametric amplification

abstract: In this paper, the scheme of a force sensor is proposed which has been composed of a hybrid optomechanical cavity containing an interacting cigar-shaped Bose-Einstein condensate (BEC) where the \textit{s}-wave scattering frequency of the BEC atoms as well as the spring coefficient of the cavity moving end-mirror (the mechanical oscillator) are parametrically modulated. It is shown that in the red-detuned regime and under the so-called impedance-matching condition, the mechanical response of the system to the input signal is enhanced substantially, which leads to the amplification of the weak input signal while the added noises of measurement (backaction noises) can be suppressed and lowered much below the standard quantum limit (SQL). In this way, such a hybrid system operates as an ultra-sensitive force sensor which can amplify the input signal and simultaneously suppress the added noises by controlling the amplitudes of modulation and the system cooperativities. The advantage of the presented nonlinear hybrid system accompanied with the mechanical and atomic modulations in comparison to the bare optomechanical cavities is the enhancement of signal amplification as well as the extension of amplification bandwidth.


title: Constraints on nuclear parton distributions from dijet photoproduction at the LHC

abstract: Using QCD calculations of the cross section of inclusive dijet photoproduction in Pb-Pb ultraperipheral collisions in the LHC kinematics as pseudo-data, we study the effect of including these data using the Bayesian reweighting technique on nCTEQ15, nCTEQ15np, and EPPS16 nuclear parton distribution functions (nPDFs). We find that, depending on the assumed error of the pseudo-data, it leads to a significant reduction of the nPDF uncertainties at small values of the momentum fraction $x_A$. Taking the error to be 5\%, the uncertainty of nCTEQ15 and nCTEQ15np nPDFs reduces approximately by a factor of two at $x_A=10^{-3}$. At the same time, the reweighting effect on EPPS16 nPDFs is much smaller due to the higher value of the tolerance and a more flexible parametrization form.


title: The H$_x$Li$_{1-x}$TaWO$_6.n$H$_2$O trirutile structure characterization through Raman and IR spectra comparison: The non-center symmetric case

abstract: We have performed the LiTaWO$_6$ vibrational studies for both Raman and infrared actives modes. Although we do not have crystalline samples, the study was conducted qualitatively, firstly comparing the sample with similar compounds in the literature, and secondly comparing both Raman and infrared spectra of the LiTaWO$_6$ compound. Although the procedure is limited, we can assign 8 of the 10 Raman-IR active modes expected by Group Theory.


title: A nonlinear hyperelasticity model for single layer blue phosphorus based on ab-initio calculations

abstract: A new hyperelastic membrane material model is proposed for single layer blue phosphorus ($\beta\text{-P}$), also known as blue phosphorene. The model is fully nonlinear and captures the anisotropy of $\beta\text{-P}$ at large strains. The material model is calibrated from density functional theory (DFT) calculations considering a set of elementary deformation states. Those are pure dilatation and uniaxial stretching along the armchair and zigzag directions. The material model is compared and validated with additional DFT results and existing DFT results from the literature, and the comparison shows good agreement. The new material model can be directly used within computational shell formulations that are for example based on rotation-free isogeometric finite elements. This is demonstrated by simulations of the indentation and vibration of single layer blue phosphorus sheets at micrometer scales. The elasticity constants at small deformations are also reported.


title: Boundary Orbits: 1 Static Spacetimes

abstract: The study of circular orbits in spacetime is of astrophysical importance. The identification and classification of circular orbits in both static and stationary spacetimes remains an active area of interest. Even in the simplest static spherically symmetric case, it is well known that the introduction of a cosmological constant in vacuum leads to the study of quartic polynomials in order to locate \textit{boundary orbits}, those that straddle between stable and unstable orbits. These orbits are often referred to as `marginally stable orbits' or `indifferently stable orbits'. A comprehensive study of texts offers little clarification as to the stability or instability of these boundary orbits. Here we argue that the direct use of second order perturbation theory immediately shows that these boundary orbits are unstable in the perturbative sense. Our study here includes the two-particle Curzon-Chazy solution.


title: Well-posed Cauchy formulation for Einstein-\ae ther theory

abstract: We study the well-posedness of the initial value (Cauchy) problem of vacuum Einstein-aether theory. The latter is a Lorentz-violating gravitational theory consisting of General Relativity with a dynamical timelike 'aether' vector field, which selects a 'preferred time' direction at each spacetime event. The Einstein-aether action is quadratic in the aether, and thus yields second order field equations for the metric and the aether. However, the well-posedness of the Cauchy problem is not easy to prove away from the simple case of perturbations over flat space. This is particularly problematic because well-posedness is a necessary requirement to ensure stability of numerical evolutions of the initial value problem. Here, we employ a first-order formulation of Einstein-aether theory in terms of projections on a tetrad frame. We show that under suitable conditions on the coupling constants of the theory, the resulting evolution equations can be cast into strongly or even symmetric hyperbolic form, and therefore they define a well-posed Cauchy problem.


title: Neural-Network Approach to Dissipative Quantum Many-Body Dynamics

abstract: In experimentally realistic situations, quantum systems are never perfectly isolated and the coupling to their environment needs to be taken into account. Often, the effect of the environment can be well approximated by a Markovian master equation. However, solving this master equation for quantum many-body systems, becomes exceedingly hard due to the high dimension of the Hilbert space. Here we present an approach to the effective simulation of the dynamics of open quantum many-body systems based on machine learning techniques. We represent the mixed many-body quantum states with neural networks in the form of restricted Boltzmann machines and derive a variational Monte-Carlo algorithm for their time evolution and stationary states. We document the accuracy of the approach with numerical examples for a dissipative spin lattice system.


title: Fast Multipole Method For 3-D Helmholtz Equation In Layered Media

abstract: In this paper, a fast multipole method (FMM) is proposed to compute long-range interactions of wave sources embedded in 3-D layered media. The layered media Green's function for the Helmholtz equation, which satisfies the transmission conditions at material interfaces, is decomposed into a free space component and four types of reaction field components arising from wave reflections and transmissions through the layered media. The proposed algorithm is a combination of the classic FMM for the free space component and FMMs specifically designed for the four types reaction components, made possible by new multipole expansions (MEs) and local expansions (LEs) as well as the multipole-to-local translation (M2L) operators for the reaction field components. { Moreover, equivalent polarization source can be defined for each reaction component based on the convergence analysis of its ME. The FMMs for the reaction components, implemented with the target particles and equivalent polarization sources, are found to be much more efficient than the classic FMM for the free space component due to the fact that the equivalent polarization sources and the target particles are always separated by a material interface.} As a result, the FMM algorithm developed for layered media has a similar computational cost as that for the free space. Numerical results validate the fast convergence of the MEs and the $O(N)$ complexity of the FMM for interactions of low-frequency wave sources in 3-D layered media.


title: Counting lines on projective surfaces

abstract: We prove a bound on the number of lines on a smooth degree-d surface in three-dimensional projective space for $d \geq 3$. This bound improves a bound due to Segre and renders some of his arguments rigorous. It is the best known bound for $d \geq 6$.


title: Efficient Continuous Multi-Query Processing over Graph Streams

abstract: Graphs are ubiquitous and ever-present data structures that have a wide range of applications involving social networks, knowledge bases and biological interactions. The evolution of a graph in such scenarios can yield important insights about the nature and activities of the underlying network, which can then be utilized for applications such as news dissemination, network monitoring, and content curation. Capturing the continuous evolution of a graph can be achieved by long-standing sub-graph queries. Although, for many applications this can only be achieved by a set of queries, state-of-the-art approaches focus on a single query scenario. In this paper, we therefore introduce the notion of continuous multi-query processing over graph streams and discuss its application to a number of use cases. To this end, we designed and developed a novel algorithmic solution for efficient multi-query evaluation against a stream of graph updates and experimentally demonstrated its applicability. Our results against two baseline approaches using real-world, as well as synthetic datasets, confirm a two orders of magnitude improvement of the proposed solution.


title: A Low-overhead Kernel Object Monitoring Approach for Virtual Machine Introspection

abstract: Monitoring kernel object modification of virtual machine is widely used by virtual-machine-introspection-based security monitors to protect virtual machines in cloud computing, such as monitoring dentry objects to intercept file operations, etc. However, most of the current virtual machine monitors, such as KVM and Xen, only support page-level monitoring, because the Intel EPT technology can only monitor page privilege. If the out-of-virtual-machine security tools want to monitor some kernel objects, they need to intercept the operation of the whole memory page. Since there are some other objects stored in the monitored pages, the modification of them will also trigger the monitor. Therefore, page-level memory monitor usually introduces overhead to related kernel services of the target virtual machine. In this paper, we propose a low-overhead kernel object monitoring approach to reduce the overhead caused by page-level monitor. The core idea is to migrate the target kernel objects to a protected memory area and then to monitor the corresponding new memory pages. Since the new pages only contain the kernel objects to be monitored, other kernel objects will not trigger our monitor. Therefore, our monitor will not introduce runtime overhead to the related kernel service. The experimental results show that our system can monitor target kernel objects effectively only with very low overhead.


title: Local loop opening in untangled ring polymer melts: A detailed "Feynman test" of models for the large scale structure

abstract: The conformational statistics of ring polymers in melts or dense solutions is strongly affected by their quenched microscopic topological state. The effect is particularly strong for untangled (i.e. non-concatenated and unknotted) rings, which are known to crumple and segregate. Here we study these systems using a computationally efficient multi-scale approach, where we combine massive simulations on the fiber level with the explicit construction of untangled ring melt configurations based on theoretical ideas for their large scale structure. We find (i) that topological constraints may be neglected on scales below the standard entanglement length, $L_e$, (ii) that rings with a size $1 \le L_r/L_e \le 30$ exhibit nearly ideal lattice tree behavior characterized by primitive paths which are randomly branched on the entanglement scale, and (iii) that larger rings are compact with gyration radii $\langle R_g^2(L_r) \rangle \propto L_r^{2/3}$. The detailed comparison between equilibrated and constructed ensembles allows us to perform a `Feynman test' of our understanding of untangled rings: can we convert ideas for the large-scale ring structure into algorithms for constructing (nearly) equilibrated ring melt samples? We show that most structural observables are quantitatively reproduced by two different construction schemes: hierarchical crumpling and ring melts derived from the analogy to interacting branched polymers. However, the latter fail the `Feynman test' with respect to the magnetic radius, $R_m$, which we have defined based on an analogy to magnetostatics. While $R_m$ is expected to vanish for double-folded structures, the observed values of $\langle R_m^2(L_r) \rangle \propto \langle R_g^2(L_r) \rangle$ provide a simple and computationally convenient measure of the presence of a non-negligible amount of local loop opening in crumpled rings.


title: The source number counts at high energies: SWIFT vs. NuSTAR

abstract: The hard X-ray sky at energies above 10 keV, has been extensively explored by the Swift and the NuSTAR missions in the 14-195 keV and the 3-24 keV bands respectively. The major population of the hard X-ray detected sources are AGN. We compare the number count distributions derived from the swift 105-month catalogue with those from the serendipitous NuSTAR source catalogue. In the 14-195 keV band, the Swift/BAT counts follow a Euclidean slope with alpha=1.51+/-0.10 (90% confidence level). The NuSTAR counts in the 8-24 keV band present a steeper slope with alpha = 1.71+/- 0.2$, because of an upturn at fluxes below ~2x10^{-13} erg/cm2/s. The same upturn is observed in the soft (3-8 keV) NuSTAR number counts, which in overall also present a steep slope with alpha = 1.82+/- 0.15. Only the bright part of the NuSTAR 3-8 keV counts agrees with the chandra number counts in the 2-10 keV band while the fainter part (<7x10^{-13} erg/cm2/s ) of the soft NuSTAR counts is in marked disagreement with the chandra counts. Next, we compare the derived number counts in the different bands using our X-ray AGN population synthesis models. The comparison between the Chandra and the Swift/BAT number counts shows a very good agreement for the 'standard' AGN spectrum with a power-law slope Gamma=1.9, a high energy cut-off at ~130 keV and a 2-10 keV reflection component of 3%. On the other hand, using the above standard AGN spectral model, only the bright part of the NuSTAR 8-24 keV and 3-8 keV number counts, agree with the model predictions. Then it is most likely that the disagreement between the X-ray number counts in the different bands is because of the faint NuSTAR number counts. We discuss various possibilities for the origin of this disagreement.


title: On the sensitivity of the diffusion MRI signal to brain activity in response to a motor cortex paradigm

abstract: Diffusion functional MRI (dfMRI) is a promising technique to map functional activations by acquiring diffusion-weighed spin-echo images. In previous studies, dfMRI showed higher spatial accuracy at activation mapping compared to classic functional MRI approaches. However, it remains unclear whether dfMRI measures result from changes in the intra-/extracellular environment, perfusion and/or T2 values. We designed an acquisition/quantification scheme to disentangle such effects in the motor cortex during a finger tapping paradigm. dfMRI was acquired at specific diffusion weightings to selectively suppress perfusion and free-water diffusion, then times series of the apparent diffusion coefficient (ADC-fMRI) and of the perfusion signal fraction (IVIM-fMRI) were derived. ADC-fMRI provided ADC estimates sensitive to changes in perfusion and free-water volume, but not to T2/T2* values. With IVIM-fMRI we isolated the perfusion contribution to ADC, while suppressing T2 effects. Compared to conventional gradient-echo BOLD fMRI, activation maps obtained with dfMRI and ADC-fMRI had smaller clusters, and the spatial overlap between the three techniques was below 50%. Increases of perfusion fractions were observed during task in both dfMRI and ADC-fMRI activations. Perfusion effects were more prominent with ADC-fMRI than with dfMRI but were significant in less than 25% of activation ROIs. Taken together, our results suggest that the sensitivity to task of dfMRI derives from a decrease of hindered diffusion and an increase of the pseudo-diffusion signal fraction, leading to different, more confined spatial activation patterns compared to classic functional MRI.


title: Equimultiplicity of families of map germs from $\mathbb{C}^2$ to $\mathbb{C}^3$

abstract: In 1971, Zariski proposed some questions in Theory of Singularities. One of such problems is the so-called, nowadays, Zariski's multiplicity conjecture. In this work, we consider the version of this conjecture for families. We answer positively Zariski's multiplicity conjecture for a special class of non-isolated singularities.


title: Discovery of Strongly-lensed Gravitational Waves - Implications for the LSST Observing Strategy

abstract: LSST's wide-field of view and sensitivity will revolutionize studies of the transient sky by finding extraordinary numbers of new transients every night. The recent discovery of a kilonova counterpart to LIGO/Virgo's first detection of gravitational waves (GWs) from a double neutron star (NS-NS) merger also creates an exciting opportunity for LSST to offer a Target of Opportunity (ToO) mode of observing. We have been exploring the possibility of detecting strongly lensed GWs, that would enable new tests of GR, extend multi-messenger astronomy out to $z\gtrsim1$, and deliver a new class of sub-millisecond precision time-delay constraints on lens mass distributions. We forecast that the rate of detection of lensed NS-NS mergers in the 2020s will be $\sim0.1$ per Earth year, that the typical source will be at $z\simeq2$, and that the multiply-imaged kilonova counterpart will have a magnitude of ${\rm AB}\simeq25.4$ in $g/r/i$-band filters - i.e. fainter than the sensitivity of a single LSST WFD visit. We therefore advocate (1) creating a flexible and efficient Target of Opportunity programme within the LSST observing strategy that is capable of discovering sources fainter than single-visit depth, and (2) surveying the entire observable extragalactic sky as rapidly as possible in the WFD survey. The latter will enable a very broad range of early science that relies on wide survey area for detection of large samples of objects and/or maximizing the fraction of sky over which reference imaging is available. For example, it will enable prompt discovery of a uniform and all-sky sample of galaxy/group/cluster-scale lenses that will underpin LSST strong-lensing science. This white paper complements submissions from DESC, SLSC, and TVSSC, that discuss kilonova, GW, and strong lensing.


title: Strong Lensing considerations for the LSST observing strategy

abstract: Strong gravitational lensing enables a wide range of science: probing cosmography; testing dark matter models; understanding galaxy evolution; and magnifying the faint, small and distant Universe. However to date exploiting strong lensing as a tool for these numerous cosmological and astrophysical applications has been severely hampered by limited sample sized. LSST will drive studies of strongly lensed galaxies, galaxy groups and galaxy clusters into the statistical age. Time variable lensing events, e.g. measuring cosmological time delays from strongly lensed supernovae and quasars, place the strongest constraints on LSST's observing strategy and have been considered in the DESC observing strategy white papers. Here we focus on aspects of `static' lens discovery that will be affected by the observing strategy. In summary, we advocate (1) ensuring comparable (sub-arcsecond) seeing in the g-band as in r and i to facilitate discovery of gravitational lenses, and (2) initially surveying the entire observable extragalactic sky as rapidly as possible to enable early science spanning a broad range of static and transient interests.


title: Electroweak Precision Tests of the Standard Model after the Discovery of the Higgs Boson

abstract: The global fit of the Standard Model predictions to electroweak precision data, which has been routinely performed in the past decades by several groups, led to the prediction of the top quark and the Higgs boson masses before their respective discoveries. With the measurement of the Higgs boson mass at the Large Hadron Collider (LHC) in 2012 by the ATLAS and CMS collaborations, the last free parameter of the Standard Model of particle physics has been fixed, and the global electroweak fit can be used to test the full internal consistency of the electroweak sector of the Standard Model and constrain models beyond. In this article, we review the current state-of-the-art theoretical calculations, as well as the precision measurements performed at the LHC, and interpret them within the context of the global electroweak fit. Special focus is drawn in the impact of the Higgs boson mass on the fit.


title: Dust unveils the formation of a mini-Neptune planet in a protoplanetary ring

abstract: Rings and radial gaps are ubiquitous in protoplanetary disks, yet their possible connection to planet formation is currently subject to intense debates. In principle, giant planet formation leads to wide gaps which separate the gas and dust mass reservoir in the outer disk, while lower mass planets lead to shallow gaps which are manifested mainly on the dust component. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the star HD169142, host to a prominent disk with deep wide gaps that sever the disk into inner and outer regions. The new ALMA high resolution images allow for the outer ring to be resolved as three narrow rings. The HD169142 disk thus hosts both the wide gaps trait of transition disks and a narrow ring system similar to those observed in the TW Hya and HL Tau systems. The mass reservoir beyond a deep gap can thus host ring systems. The observed rings are narrow in radial extent (width/radius of 1.5/57.3, 1.8/64.2 and 3.4/76.0, in au) and have asymmetric mutual separations: the first and middle ring are separated by 7 au while the middle and outermost ring are distanced by ~12 au. Using hydrodynamical modeling we found that a simple explanation, involving a single migrating low mass planet (10 M$_\oplus$), entirely accounts for such an apparently complex phenomenon. Inward migration of the planet naturally explains the ring's asymmetric mutual separation. The isolation of HD169142's outer rings thus allows a proof of concept to interpret the detailed architecture of the outer region of protoplanetary disks with low mass planet formation of mini-Neptune's size, i.e. as in the protosolar nebula.


title: The Origin of Kepler-419b: A Path to Tidal Migration Via Four-body Secular Interactions

abstract: We test the high-eccentricity tidal migration scenario for Kepler-419b, a member of the eccentric warm Jupiter class of planets whose origin is debated. Kepler-419 hosts two known planets (b,c). However, in its current configuration, planet c cannot excite the eccentricity of planet b enough to undergo high-eccentricity tidal migration. We investigate whether the presence of an undiscovered fourth body could explain the orbit of Kepler-419b. We explore the parameter space of this potential third giant planet using a suite of N-body simulations with a range of initial conditions. From the results of these simulations, coupled with observational constraints, we can rule out this mechanism for much of the parameter space of initial object d conditions. However, for a small range of parameters (masses between 0.5 and 7 $m_{\rm{Jup}}$, semi-major axes between 4 and 7.5 AU, eccentricities between 0.18 and 0.35, and mutual inclinations near 0$^{\circ}$) an undiscovered object d could periodically excite the eccentricity of Kepler-419b without destabilizing the system over 1 Gyr while producing currently undetectable radial velocity and transit timing variation signals.


title: A Robust Riemann Solver for Multiple Hydro-Elastoplastic Solid Mediums

abstract: We propose a robust approximate solver for the hydro-elastoplastic solid material, a general constitutive law extensively applied in explosion and high speed impact dynamics, and provide a natural transformation between the fluid and solid in the case of phase transitions. The hydrostatic components of the solid is described by a family of general Mie-Gr\"uneisen equation of state (EOS), while the deviatoric component includes the elastic phase, linearly hardened plastic phase and fluid phase. The approximate solver provides the interface stress and normal velocity by an iterative method. The well-posedness and convergence of our solver are proved with mild assumptions on the equations of state. The proposed solver is applied in computing the numerical flux at the phase interface for our compressible multi-medium flow simulation on Eulerian girds. Several numerical examples, including Riemann problems, shock-bubble interactions, implosions and high speed impact applications, are presented to validate the approximate solver.


title: On the CVP for the root lattices via folding with deep ReLU neural networks

abstract: Point lattices and their decoding via neural networks are considered in this paper. Lattice decoding in Rn, known as the closest vector problem (CVP), becomes a classification problem in the fundamental parallelotope with a piecewise linear function defining the boundary. Theoretical results are obtained by studying root lattices. We show how the number of pieces in the boundary function reduces dramatically with folding, from exponential to linear. This translates into a two-layer ReLU network requiring a number of neurons growing exponentially in n to solve the CVP, whereas this complexity becomes polynomial in n for a deep ReLU network.


title: An algorithmic approach to the multiple impact of a disk in a corner

abstract: We present the algorithmic procedure determining the impulsive behavior of a rigid disk having a single or possibly multiple frictionless impact with two walls forming a corner. The algorithmic procedure represents an application of the general theory of multiple impacts as presented in \cite{Pasquero2016Multiple} for the ideal case. In the first part, two theoretical algorithms are presented for the cases of ideal impact and Newtonian frictionless impact with global dissipation index. The termination analysis of the algorithms differentiates the two cases: in the ideal case, we show that the algorithm always terminates and the disk exits from the corner after a finite number of steps independently of the initial impact velocity of the disk and the angle formed by the walls; in the non--ideal case, although is not proved that the disk exits from the corner in a finite number of steps, we show that its velocity decreases to zero and the termination of the algorithm can be fixed through an "almost at rest" condition. In the second part, we present a numerical version of both the theoretical algorithms that is more robust than the theoretical ones with respect to noisy initial data and floating point arithmetic computation. Moreover, we list and analyze the outputs of the numerical algorithm in several cases.


title: Probabilistic Generative Deep Learning for Molecular Design

abstract: Probabilistic generative deep learning for molecular design involves the discovery and design of new molecules and analysis of their structure, properties and activities by probabilistic generative models using the deep learning approach. It leverages the existing huge databases and publications of experimental results, and quantum-mechanical calculations, to learn and explore molecular structure, properties and activities. We discuss the major components of probabilistic generative deep learning for molecular design, which include molecular structure, molecular representations, deep generative models, molecular latent representations and latent space, molecular structure-property and structure-activity relationships, molecular similarity and molecular design. We highlight significant recent work using or applicable to this new approach.


title: Seismic characterization of multiple fracture sets at Rulison Field, Colorado

abstract: Conventional fracture-characterization methods assume the presence of a single set of aligned, vertical cracks in the subsurface. We relax this assumption and demonstrate the feasibility of seismic characterization of multiple fracture sets. Our technique relies on recent numerical findings, indicating that multiple, differently oriented, possibly intersecting planar cracks embedded in an otherwise isotropic host rock result in a nearly orthorhombic (or orthotropic) effective medium. Here, the governing parameters of crack-induced orthotropy are estimated from 3D, wide-azimuth, multicomponent seismic reflection data acquired over the tight-gas Rulison Field in Colorado, USA. We translate strong azimuthal variations of the normal-moveout velocities into interval crack densities, fracture orientations, type of fluid infill, and velocities of the P- and S-waves in an unfractured rock. Our inversion procedure identifies a set of cracks aligned in approximately WNW-ESE direction in the western part of the study area and multiple, likely intersecting fractures in its eastern part. We validate both our underlying theoretical model and the obtained estimates by two independent measurements: (i) the estimated fluid-infill parameter indicates dry cracks as expected for the gas-producing sandstones at Rulison; and (ii) the obtained crack orientations are supported by measurements in a well. As a by-product of fracture characterization, we build an anisotropic velocity model of the Rulison reservoir that, we believe, is the first orthorhombic velocity field constructed from surface seismic data.


title: A fiber-integrated single photon source emitting at telecom wavelengths

abstract: Fiber-coupled single photon sources are essential components of photonics-based quantum information processors. Most fiber-coupled single photon sources require careful alignment between fibers and quantum emitters. In this work, we present an alignment-free fiber-integrated single photon source based on an InAs/InP quantum dot emitting at telecom wavelengths. We designed a nanobeam containing the quantum dots attached to a fiber taper. The adiabatic tapered coupler of the nanobeam enables efficient light coupling to the fiber taper. Using a tungsten probe in a focused ion beam system, we transferred the nanobeam to the fiber taper. The observed fiber-coupled single photon emission occurs with a brightness of 1.5% and purity of 86%. This device provides a building block for fiber-optic quantum circuits that have various applications, such as quantum communication and distributed quantum computing.


title: Correlated Insulating and Superconducting States in Twisted Bilayer Graphene Below the Magic Angle

abstract: The emergence of flat bands and correlated behaviors in 'magic angle' twisted bilayer graphene (tBLG) has sparked tremendous interest, though many aspects of the system are under intense debate. Here we report observation of both superconductivity and the Mott-like insulating state in a tBLG device with a twist angle of approximately 0.93, which is smaller than the magic angle by 15%. At an electron concentration of +/-5 electrons per moire unit cell, we observe a narrow resistance peak with an activation energy gap of approximately 0.1 meV, indicating the existence of an additional correlated insulating state. This is consistent with theory predicting the presence of a high-energy band with an energetically flat dispersion. At a doping of +/-12 electrons per moire unit cell we observe a resistance peak due to the presence of Dirac points in the spectrum. Our results reveal that the magic range of tBLG is in fact larger than what is previously expected, and provide a wealth of new information to help decipher the strongly correlated phenomena observed in tBLG.


title: The Cost of Monitoring Alone

abstract: We compare the succinctness of two monitoring systems for properties of infinite traces, namely parallel and regular monitors. Although a parallel monitor can be turned into an equivalent regular monitor, the cost of this transformation is a double-exponential blowup in the syntactic size of the monitors, and a triple-exponential blowup when the goal is a deterministic monitor. We show that these bounds are tight and that they also hold for translations between corresponding fragments of Hennessy-Milner logic with recursion over infinite traces.


title: GROJ1750-27: a neutron star far behind the Galactic Center switching into the propeller regime

abstract: We report on analysis of properties of the X-ray binary pulsar GROJ1750-27 based on X-ray (Chandra, Swift, and Fermi/GBM), and near-infrared (VVV and UKIDSS surveys) observations. An accurate position of the source is determined for the first time and used to identify its infrared counterpart. Based on the VVV data we investigate the spectral energy distribution (SED) of the companion, taking into account a non-standard absorption law in the source direction. A comparison of this SED with those of known Be/X-ray binaries and early type stars has allowed us to estimate a lower distance limit to the source at $>12$ kpc. An analysis of the observed spin-up torque during a giant outburst in 2015 provides an independent distance estimate of $14-22$ kpc, and also allows to estimate the magnetic field on the surface of the neutron star at $B\simeq(3.5-4.5)\times10^{12}$ G. The latter value is in agreement with the possible transition to the propeller regime, a strong hint for which was revealed by Swift/XRT and Chandra. We conclude that GROJ1750-27 is located far behind the Galactic Center, which makes it one of the furthest Galactic X-ray binaries known.


title: $L^1$-spaces of vector measures with vector density

abstract: Let $F$ be a function with values in a Banach space. When $F$ is locally (Pettis or Bochner) integrable with respect to a locally determined positive measure, a vector measure $\nu_F$ with density $F$ defined on a $\delta$-ring is obtained. We present the existing connection between the spaces $L^1_w(\nu_F)$, $L^1(\nu_F)$ and $L^1(|\nu_F|)$ and the spaces of Dunford, Pettis or Bochner integrable functions.


title: An instance where the major and minor arc integrals meet

abstract: We apply the circle method to obtain an asymptotic formula for the number of integral points on a certain sliced cubic hypersurface related to the Segre cubic. Unusually, the major and minor arc integrals in this application are both positive and of the same order of magnitude.


title: Multiple Systems Estimation for Sparse Capture Data: Inferential Challenges when there are Non-Overlapping Lists

abstract: Multiple systems estimation strategies have recently been applied to quantify hard-to-reach populations, particularly when estimating the number of victims of human trafficking and modern slavery. In such contexts, it is not uncommon to see sparse or even no overlap between some of the lists on which the estimates are based. These create difficulties in model fitting and selection, and we develop inference procedures to address these challenges. The approach is based on Poisson log-linear regression modeling. Issues investigated in detail include taking proper account of data sparsity in the estimation procedure, as well as the existence and identifiability of maximum likelihood estimates. A stepwise method for choosing the most suitable parameters is developed, together with a bootstrap approach to finding confidence intervals for the total population size. We apply the strategy to two empirical data sets of trafficking in US regions, and find that the approach results in stable, reasonable estimates. An accompanying R software implementation has been made publicly available.


title: The Role of Energy Minimization in Algebraic Multigrid Interpolation

abstract: Algebraic multigrid (AMG) methods are powerful solvers with linear or near-linear computational complexity for certain classes of linear systems, Ax=b. Broadening the scope of problems that AMG can effectively solve requires the development of improved interpolation operators. Such development is often based on AMG convergence theory. However, convergence theory in AMG tends to have a disconnect with AMG in practice due to the practical constraints of (i) maintaining matrix sparsity in transfer and coarse-grid operators, and (ii) retaining linear complexity in the setup and solve phase. This paper presents a review of fundamental results in AMG convergence theory, followed by a discussion on how these results can be used to motivate interpolation operators in practice. A general weighted energy minimization functional is then proposed to form interpolation operators, and a novel `diagonal' preconditioner for Sylvester- or Lyapunov-type equations developed simultaneously. Although results based on the weighted energy minimization typically underperform compared to a fully constrained energy minimization, numerical results provide new insight into the role of energy minimization and constraint vectors in AMG interpolation.


title: Unsupervised machine learning in atomistic simulations, between predictions and understanding

abstract: Automated analyses of the outcome of a simulation have been an important part of atomistic modeling since the early days, addressing the need of linking the behavior of individual atoms and the collective properties that are usually the final quantity of interest. Methods such as clustering and dimensionality reduction have been used to provide a simplified, coarse-grained representation of the structure and dynamics of complex systems, from proteins to nanoparticles. In recent years, the rise of machine learning has led to an even more widespread use of these algorithms in atomistic modeling, and to consider different classification and inference techniques as part of a coherent toolbox of data-driven approaches. This perspective briefly reviews some of the unsupervised machine-learning methods -- that are geared towards classification and coarse-graining of molecular simulations -- seen in relation to the fundamental mathematical concepts that underlie all machine-learning techniques. It discusses the importance of using concise yet complete representations of atomic structures as the starting point of the analyses, and highlights the risk of introducing preconceived biases when using machine learning to rationalize and understand structure-property relations. Supervised machine-learning techniques, that explicitly attempt to predict the properties of a material given its structure, are less susceptible to such biases. Current developments in the field suggest that using these two classes of approaches side-by-side and in a fully integrated mode, while keeping in mind the relations between the data analysis framework and the fundamental physical principles, will be key for realizing the full potential of machine learning to help understanding the behavior of complex molecules and materials.


title: Observing Microscopic Transitions from Macroscopic Bursts: Instability-Mediated Resetting in the Incoherent Regime of the $D$-dimensional Generalized Kuramoto Model

abstract: We consider a recently introduced $D$-dimensional generalized Kuramoto model for many $(N\gg 1)$ interacting agents in which the agent states are $D$-dimensional unit vectors. It was previously shown that, for even $D$, similar to the original Kuramoto model ($D=2$), there exists a continuous dynamical phase transition from incoherence to coherence of the time asymptotic attracting state as the coupling parameter $K$ increases through a critical value $K_c^{(+)}>0$. We consider this transition from the point of view of the stability of an incoherent state, i.e., where the $N\to\infty$ distribution function is time-independent and the macroscopic order parameter is zero. In contrast with $D=2$, for even $D>2$ there is an infinity of possible incoherent equilibria, each of which becomes unstable with increasing $K$ at a different point $K=K_c$. We show that there are incoherent equilibria for all $K_c$ within the range $(K_c^{(+)}/2)\leq K_c \leq K_c^{(+)}$. How can the possible instability of incoherent states arising at $K=K_c<K_c^{(+)}$ be reconciled with the previous finding that, at large time, the state is always incoherent unless $K>K_c^{(+)}$? We find, for a given incoherent equilibrium, that, if $K$ is rapidly increased from $K<K_c$ to $K_c<K<K_c^{(+)}$, due to the instability, a short, macroscopic burst of coherence is observed, which initially grows exponentially, but then reaches a maximum, past which it decays back into incoherence. After this decay, we observe that the equilibrium has been reset to one whose $K_c$ value exceeds that of the increased $K$. Thus this process, which we call `Instability-Mediated Resetting,' leads to an increase in the effective $K_c$ with continuously increasing $K$, until the equilibrium has been effectively set to one for which for which $K_c\approx K_c^{(+)}$, leading to a unique critical point of the time asymptotic state.


title: Ultrastrong time-dependent light-matter interactions are gauge-relative

abstract: Time-dependent light-matter interactions are a widespread means by which to describe controllable experimental operations. They can be viewed as an approximation in which a third system - the control system - is treated as external within the Hamiltonian. We demonstrate that this results in non-equivalence between gauges. We provide a physical example in which each different non-equivalent model coincides with a gauge-invariant description applied in a different experimental situation. The qualitative final-time predictions obtained from these models, including entanglement and photon number, depend on the gauge within which the time-dependent coupling assumption is made. This occurs whenever the interaction switching is sufficiently strong and non-adiabatic even if the coupling vanishes at the preparation and measurement stages of the protocol, at which times the subsystems are unique and experimentally addressable.


title: Fitting rectangles under vulnerability curves: optimal water flow through plants

abstract: We study an optimization problem for a model of steady state water transport through plants that maximizes water flow subject to the constraints on hydraulic conductance due to vulnerability to embolism (air blockage of conduits). The model has an elementary geometric interpretation, and exhibits bottleneck behavior where one of the plant segments limits the overall optimal flow, sometimes in a counterintuitive way. The results show good agreement with experimental measurements and provide support for the hypothesis that leaves serve as a safety buffer protecting stems against excessive embolism.


title: Quantum Language Processing

abstract: We present a representation for linguistic structure that we call a Fock-space representation, which allows us to embed problems in language processing into small quantum devices. We further develop a formalism for understanding both classical as well as quantum linguistic problems and phrase them both as a Harmony optimization problem that can be solved on a quantum computer which we show is related to classifying vectors using quantum Boltzmann machines. We further provide a new training method for learning quantum Harmony operators that describe a language. This also provides a new algorithm for training quantum Boltzmann machines that requires no approximations and works in the presence of hidden units. We additionally show that quantum language processing is BQP-complete, meaning that it is polynomially equivalent to the circuit model of quantum computing which implies that quantum language models are richer than classical models unless BPP=BQP. It also implies that, under certain circumstances, quantum Boltzmann machines are more expressive than classical Boltzmann machines. Finally, we examine the performance of our approach numerically for the case of classical harmonic grammars and find that the method is capable of rapidly parsing even non-trivial grammars. This suggests that the work may have value as a quantum inspired algorithm beyond its initial motivation as a new quantum algorithm.


title: Prediction of Glass Elasticity from Free Energy Density of Topological Constraints

abstract: Despite the critical importance of the elastic properties of modern materials, there is not a singular model that can predict the modulus to an accuracy needed for industrial glass design. To address this problem, we propose an approach to calculate the elastic modulus based on the free energy density of topological constraints in the glass-forming network. Our approach shows quantitatively accurate agreement with glasses across a variety of compositional families. Moreover, using temperature-dependent constraint theory, the temperature dependence of the modulus can also be predicted. Our approach is general and theoretically can be applied to any network glass.


title: SoK: Transparent Dishonesty: front-running attacks on Blockchain

abstract: We consider front-running to be a course of action where an entity benefits from prior access to privileged market information about upcoming transactions and trades. Front-running has been an issue in financial instrument markets since the 1970s. With the advent of the blockchain technology, front-running has resurfaced in new forms we explore here, instigated by blockchains decentralized and transparent nature. In this paper, we draw from a scattered body of knowledge and instances of front-running across the top 25 most active decentral applications (DApps) deployed on Ethereum blockchain. Additionally, we carry out a detailed analysis of Status.im initial coin offering (ICO) and show evidence of abnormal miners behavior indicative of front-running token purchases. Finally, we map the proposed solutions to front-running into useful categories.


title: Non-thermal emission from the interaction of magnetized exoplanets with the wind of their host star

abstract: We study the non-thermal emission from the interaction between magnetized Jupiter-like exoplanets and the wind from their host star. The supersonic motion of planets through the wind forms a bow shock that accelerates electrons which produces non-thermal radiation across a broad wavelength range. We discuss three wind mass loss rates: $\dot{M}_{\rm w}\sim10^{-14}$, $10^{-9}$, $10^{-6}\,M_{\odot}\,\rm yr^{-1}$ corresponding to solar-type, T Tauri and massive O/B type stars, respectively. We find that the expected radio synchrotron emission from a Jupiter-like planet is detectable by the Jansky Very Large Array and the Square Kilometer Array at ~1-10 GHz out to a distance ~ 100 pc, whereas the infrared emission is detectable by the James Webb Space Telescope out to a similar distance. Inverse Compton scattering of the stellar radiation results in X-ray emission detectable by Chandra X-ray Observatory out to ~ 150 pc. Finally, we apply our model to the upper limit constraints on V380 Tau, the first star-hot Jupiter system observed in radio wavelength. Our bow shock model provides constraints on the magnetic field, the interplanetary medium and the non-thermal emission efficiency in V380 Tau.


title: Space-Efficient Data Structures for Lattices

abstract: A lattice is a partially-ordered set in which every pair of elements has a unique meet (greatest lower bound) and join (least upper bound). We present new data structures for lattices that are simple, efficient, and nearly optimal in terms of space complexity. Our first data structure can answer partial order queries in constant time and find the meet or join of two elements in $O(n^{3/4})$ time, where $n$ is the number of elements in the lattice. It occupies $O(n^{3/2}\log n)$ bits of space, which is only a $\Theta(\log n)$ factor from the $\Theta(n^{3/2})$-bit lower bound for storing lattices. The preprocessing time is $O(n^2)$. This structure admits a simple space-time tradeoff so that, for any $c \in [\frac{1}{2}, 1]$, the data structure supports meet and join queries in $O(n^{1-c/2})$ time, occupies $O(n^{1+c}\log n)$ bits of space, and can be constructed in $O(n^2 + n^{1+3c/2})$ time. Our second data structure uses $O(n^{3/2}\log n)$ bits of space and supports meet and join in $O(d \frac{\log n}{\log d})$ time, where $d$ is the maximum degree of any element in the transitive reduction graph of the lattice. This structure is much faster for lattices with low-degree elements. This paper also identifies an error in a long-standing solution to the problem of representing lattices. We discuss the issue with this previous work.


title: Some Interesting Features of Memristor CNN

abstract: In this paper, we introduce some interesting features of a memristor CNN (Cellular Neural Network). We first show that there is the similarity between the dynamics of memristors and neurons. That is, some kind of flux-controlled memristors can not respond to the sinusoidal voltage source quickly, namely, they can not switch `on' rapidly. Furthermore, these memristors have refractory period after switch `on', which means that it can not respond to further sinusoidal inputs until the flux is decreased. We next show that the memristor-coupled two-cell CNN can exhibit chaotic behavior. In this system, the memristors switch `off' and `on' at irregular intervals, and the two cells are connected when either or both of the memristors switches `on'. We then propose the modified CNN model, which can hold a binary output image, even if all cells are disconnected and no signal is supplied to the cell after a certain point of time. However, the modified CNN requires power to maintain the output image, that is, it is volatile. We next propose a new memristor CNN model. It can also hold a binary output state (image), even if all cells are disconnected, and no signal is supplied to the cell, by memristor's switching behavior. Furthermore, even if we turn off the power of the system during the computation, it can resume from the previous average output state, since the memristor CNN has functions of both short-term (volatile) memory and long-term (non-volatile) memory. The above suspend and resume feature are useful when we want to save the current state, and continue work later from the previous state. Finally, we show that the memristor CNN can exhibit interesting two-dimensional waves, if an inductor is connected to each memristor CNN cell.


title: Nonlinear Depolarization of Light in Optical Communication Fiber

abstract: We report the experimental observation of a novel transmission phenomenon in optical long-haul communication systems. Un-polarized ASE depolarizes via nonlinear fiber interactions a cw laser light during their co-propagation which leads to small but measurable ultra-fast polarization state fluctuations at the fiber output. We provide a phenomenological approach and a theory that qualitatively corroborates our experimental results. One of our major findings suggests that the applicability of the often used Manakov equation needs to be scrutinized for highly accurate studies of nonlinear polarization state evolutions in noisy environments. The described phenomenon leads to a qualitatively different microscopic understanding of nonlinear light propagation in fiber and can contribute towards an explanation for todays commonly perceived gap between simulated and experimentally obtained system performance in optical data transmission.


title: Quantized \Delta S=2 Excitation Spectra by Confinement in an S=1 Spin Chain

abstract: We calculate the dynamical spin-structure factor of the $S=1$ Ising spin chain with negative single-ion anisotropy in magnetic fields using the infinite time-evolving-block-decimation algorithm. We show that when a transverse magnetic field is applied, both the $\Delta S=2$ excitation continuum and one-magnon mode appear in the low-lying excitation. When a longitudinal magnetic field is further applied, the excitation continuum changes into quantized excitation spectra. The quantized $\Delta S=2$ excitation spectra originate from the confinement of two domain walls, each of which carries $\Delta S=1$. The quantized excitation energies are explained by the negative zeros in the Airy function.


title: GrapAL: Connecting the Dots in Scientific Literature

abstract: We introduce GrapAL (Graph database of Academic Literature), a versatile tool for exploring and investigating a knowledge base of scientific literature, that was semi-automatically constructed using NLP methods. GrapAL satisfies a variety of use cases and information needs requested by researchers. At the core of GrapAL is a Neo4j graph database with an intuitive schema and a simple query language. In this paper, we describe the basic elements of GrapAL, how to use it, and several use cases such as finding experts on a given topic for peer reviewing, discovering indirect connections between biomedical entities and computing citation-based metrics. We open source the demo code to help other researchers develop applications that build on GrapAL.


title: Accelerating dynamical peakons and their behaviour

abstract: A wide class of nonlinear dispersive wave equations are shown to possess a novel type of peakon solution in which the amplitude and speed of the peakon are time-dependent. These novel dynamical peakons exhibit a wide variety of different behaviours for their amplitude, speed, and acceleration, including an oscillatory amplitude and constant speed which describes a peakon breather. Examples are presented of families of nonlinear dispersive wave equations that illust rate various interesting behaviours, such as asymptotic travelling-wave peakons, dissipating/anti-dissipating peakons, direction-reversing peakons, runaway and blow up peakons, among others.


title: Computability logic: Giving Caesar what belongs to Caesar

abstract: The present article is a brief informal survey of computability logic --- the game-semantically conceived formal theory of computational resources and tasks. This relatively young nonclassical logic is a conservative extension of classical first order logic but is much more expressive than the latter, yielding a wide range of new potential application areas. In a reasonable (even if not strict) sense the same holds for intuitionistic and linear logics, which allows us to say that CoL reconciles and unifies the three traditions of logical thought (and beyond) on the basis of its natural and "universal" game semantics. A comprehensive online survey of the subject can be found at http://www.csc.villanova.edu/~japaridz/CL/ .


title: Estimation of Gaussian directed acyclic graphs using partial ordering information with an application to dairy cattle data

abstract: Estimating a directed acyclic graph (DAG) from observational data represents a canonical learning problem and has generated a lot of interest in recent years. Research has focused mostly on the following two cases: when no information regarding the ordering of the nodes in the DAG is available, and when a domain-specific complete ordering of the nodes is available. In this paper, motivated by a recent application in dairy science, we develop a method for DAG estimation for the middle scenario, where partition based partial ordering of the nodes is known based on domain specific knowledge.We develop an efficient algorithm that solves the posited problem, coined Partition-DAG. Through extensive simulations using the DREAM3 Yeast data, we illustrate that Partition-DAG effectively incorporates the partial ordering information to improve both speed and accuracy. We then illustrate the usefulness of Partition-DAG by applying it to recently collected dairy cattle data, and inferring relationships between various variables involved in dairy agroecosystems.


title: Global solutions to the supercooled Stefan problem with blow-ups: regularity and uniqueness

abstract: We consider the supercooled Stefan problem, which captures the freezing of a supercooled liquid, in one space dimension. A probabilistic reformulation of the problem allows to define global solutions, even in the presence of blow-ups of the freezing rate. We provide a complete description of such solutions, by relating the temperature distribution in the liquid to the regularity of the ice growth process. The latter is shown to transition between (i) continuous differentiability, (ii) Holder continuity, and (iii) discontinuity. In particular, in the second regime we rediscover the square root behavior of the growth process pointed out by Stefan in his seminal paper [Ste89] from 1889 for the ordinary Stefan problem. In our second main theorem, we establish the uniqueness of the global solutions, a first result of this kind in the context of growth processes with singular self-excitation when blow-ups are present.


title: Dualities of Corner Configurations and Supersymmetric Indices

abstract: We compute supersymmetric indices which count local operators at certain half-BPS interfaces and quarter-BPS junctions of interfaces in four-dimensional $\mathcal{N}=4$ Super Yang-Mills theory. We use the indices as very stringent tests of a variety of string theory-inspired conjectures about the action of S-duality on such defects.


title: Toward Ergonomic Risk Prediction via Segmentation of Indoor Object Manipulation Actions Using Spatiotemporal Convolutional Networks

abstract: Automated real-time prediction of the ergonomic risks of manipulating objects is a key unsolved challenge in developing effective human-robot collaboration systems for logistics and manufacturing applications. We present a foundational paradigm to address this challenge by formulating the problem as one of action segmentation from RGB-D camera videos. Spatial features are first learned using a deep convolutional model from the video frames, which are then fed sequentially to temporal convolutional networks to semantically segment the frames into a hierarchy of actions, which are either ergonomically safe, require monitoring, or need immediate attention. For performance evaluation, in addition to an open-source kitchen dataset, we collected a new dataset comprising twenty individuals picking up and placing objects of varying weights to and from cabinet and table locations at various heights. Results show very high (87-94)\% F1 overlap scores among the ground truth and predicted frame labels for videos lasting over two minutes and consisting of a large number of actions.


title: Multi-Objective Policy Generation for Multi-Robot Systems Using Riemannian Motion Policies

abstract: In many applications, multi-robot systems are required to achieve multiple objectives. For these multi-objective tasks, it is oftentimes hard to design a single control policy that fulfills all the objectives simultaneously. In this paper, we focus on multi-objective tasks that can be decomposed into a set of simple subtasks. Controllers for these subtasks are individually designed and then combined into a control policy for the entire team. One significant feature of our work is that the subtask controllers are designed along with their underlying manifolds. When a controller is combined with other controllers, their associated manifolds are also taken into account. This formulation yields a policy generation framework for multi-robot systems that can combine controllers for a variety of objectives while implicitly handling the interaction among robots and subtasks. To describe controllers on manifolds, we adopt Riemannian Motion Policies (RMPs), and propose a collection of RMPs for common multi-robot subtasks. Centralized and decentralized algorithms are designed to combine these RMPs into a final control policy. Theoretical analysis shows that the system under the control policy is stable. Moreover, we prove that many existing multi-robot controllers can be closely approximated by the framework. The proposed algorithms are validated through both simulated tasks and robotic implementations.


title: Spectre is here to stay: An analysis of side-channels and speculative execution

abstract: The recent discovery of the Spectre and Meltdown attacks represents a watershed moment not just for the field of Computer Security, but also of Programming Languages. This paper explores speculative side-channel attacks and their implications for programming languages. These attacks leak information through micro-architectural side-channels which we show are not mere bugs, but in fact lie at the foundation of optimization. We identify three open problems, (1) finding side-channels, (2) understanding speculative vulnerabilities, and (3) mitigating them. For (1) we introduce a mathematical meta-model that clarifies the source of side-channels in simulations and CPUs. For (2) we introduce an architectural model with speculative semantics to study recently-discovered vulnerabilities. For (3) we explore and evaluate software mitigations and prove one correct for this model. Our analysis is informed by extensive offensive research and defensive implementation work for V8, the production JavaScript virtual machine in Chrome. Straightforward extensions to model real hardware suggest these vulnerabilities present formidable challenges for effective, efficient mitigation. As a result of our work, we now believe that speculative vulnerabilities on today's hardware defeat all language-enforced confidentiality with no known comprehensive software mitigations, as we have discovered that untrusted code can construct a universal read gadget to read all memory in the same address space through side-channels. In the face of this reality, we have shifted the security model of the Chrome web browser and V8 to process isolation.


title: Multi-task learning with compressible features for Collaborative Intelligence

abstract: A promising way to deploy Artificial Intelligence (AI)-based services on mobile devices is to run a part of the AI model (a deep neural network) on the mobile itself, and the rest in the cloud. This is sometimes referred to as collaborative intelligence. In this framework, intermediate features from the deep network need to be transmitted to the cloud for further processing. We study the case where such features are used for multiple purposes in the cloud (multi-tasking) and where they need to be compressible in order to allow efficient transmission to the cloud. To this end, we introduce a new loss function that encourages feature compressibility while improving system performance on multiple tasks. Experimental results show that with the compression-friendly loss, one can achieve around 20% bitrate reduction without sacrificing the performance on several vision-related tasks.


title: The Many-to-Many Mapping Between the Concordance Correlation Coefficient and the Mean Square Error

abstract: We derive the mapping between two of the most pervasive utility functions, the mean square error ($MSE$) and the concordance correlation coefficient (CCC, $\rho_c$). Despite its drawbacks, $MSE$ is one of the most popular performance metrics (and a loss function); along with lately $\rho_c$ in many of the sequence prediction challenges. Despite the ever-growing simultaneous usage, e.g., inter-rater agreement, assay validation, a mapping between the two metrics is missing, till date. While minimisation of $L_p$ norm of the errors or of its positive powers (e.g., $MSE$) is aimed at $\rho_c$ maximisation, we reason the often-witnessed ineffectiveness of this popular loss function with graphical illustrations. The discovered formula uncovers not only the counterintuitive revelation that `$MSE_1<MSE_2$' does not imply `$\rho_{c_1}>\rho_{c_2}$', but also provides the precise range for the $\rho_c$ metric for a given $MSE$. We discover the conditions for $\rho_c$ optimisation for a given $MSE$; and as a logical next step, for a given set of errors. We generalise and discover the conditions for any given $L_p$ norm, for an even p. We present newly discovered, albeit apparent, mathematical paradoxes. The study inspires and anticipates a growing use of $\rho_c$-inspired loss functions e.g., $\left|\frac{MSE}{\sigma_{XY}}\right|$, replacing the traditional $L_p$-norm loss functions in multivariate regressions.


title: Data Correlation-Aware Resource Management in Wireless Virtual Reality (VR): An Echo State Transfer Learning Approach

abstract: In this paper, the problem of wireless virtual reality (VR) resource management is investigated for a wireless VR network in which VR contents are sent by a cloud to cellular small base stations (SBSs). The SBSs will collect tracking data from the VR users, over the uplink, in order to generate the VR content and transmit it to the end-users using downlink cellular links. For this model, the data requested or transmitted by the users can exhibit correlation, since the VR users may engage in the same immersive virtual environment with different locations and orientations. As such, the proposed resource management framework can factor in such spatial data correlation, to better manage uplink and downlink traffic. This potential spatial data correlation can be factored into the resource allocation problem to reduce the traffic load in both uplink and downlink. In the downlink, the cloud can transmit 360 contents or specific visible contents that are extracted from the original 360 contents to the users according to the users' data correlation to reduce the backhaul traffic load. For uplink, each SBS can associate with the users that have similar tracking information so as to reduce the tracking data size. This data correlation-aware resource management problem is formulated as an optimization problem whose goal is to maximize the users' successful transmission probability, defined as the probability that the content transmission delay of each user satisfies an instantaneous VR delay target. To solve this problem, an echo state networks (ESNs) based transfer learning is introduced. By smartly transferring information on the SBS's utility, the proposed transfer-based ESN algorithm can quickly cope with changes in the wireless networking environment.


title: On reconstruction in the inverse conductivity problem with one measurement

abstract: We consider an inverse problem for electrically conductive material occupying a domain $\Omega$ in $\Bbb R^2$. Let $\gamma$ be the conductivity of $\Omega$, and $D$ a subdomain of $\Omega$. We assume that $\gamma$ is a positive constant $k$ on $D$, $k\not=1$ and is $1$ on $\Omega\setminus D$; both $D$ and $k$ are unknown. The problem is to find a reconstruction formula of $D$ from the Cauchy data on $\partial\Omega$ of a non-constant solution $u$ of the equation $\nabla\cdot\gamma\nabla u=0$ in $\Omega$. We prove that if $D$ is known to be a convex polygon such that $\text{diam}\,D<\text{dist}\,(D,\partial\Omega)$, there are two formulae for calculating the support function of $D$ from the Cauchy data.


title: Manipulating Soft Tissues by Deep Reinforcement Learning for Autonomous Robotic Surgery

abstract: In robotic surgery, pattern cutting through a deformable material is a challenging research field. The cutting procedure requires a robot to concurrently manipulate a scissor and a gripper to cut through a predefined contour trajectory on the deformable sheet. The gripper ensures the cutting accuracy by nailing a point on the sheet and continuously tensioning the pinch point to different directions while the scissor is in action. The goal is to find a pinch point and a corresponding tensioning policy to minimize damage to the material and increase cutting accuracy measured by the symmetric difference between the predefined contour and the cut contour. Previous study considers finding one fixed pinch point during the course of cutting, which is inaccurate and unsafe when the contour trajectory is complex. In this paper, we examine the soft tissue cutting task by using multiple pinch points, which imitates human operations while cutting. This approach, however, does not require the use of a multi-gripper robot. We use a deep reinforcement learning algorithm to find an optimal tensioning policy of a pinch point. Simulation results show that the multi-point approach outperforms the state-of-the-art method in soft pattern cutting task with respect to both accuracy and reliability.


title: A Covariance-Based Hybrid Channel Feedback in FDD Massive MIMO Systems

abstract: In this paper, a novel covariance-based channel feedback mechanism is investigated for frequency division duplexing (FDD) massive multi-input multi-output (MIMO) systems. The concept capitalizes on the notion of user statistical separability which was hinted in several prior works in the massive antenna regime but not fully exploited so far. We here propose a hybrid statistical-instantaneous feedback mechanism where the users are separated into two classes of feedback design based on their channel covariance. Under the hybrid framework, each user either operates on a statistical feedback mode or quantized instantaneous channel feedback mode depending on their so-called statistical isolability. The key challenge lies in the design of a covariance-aware classification algorithm which can handle the complex mutual interactions between all users. The classification is derived from rate bound principles. A suitable precoding method is also devised under the mixed statistical and instantaneous feedback model. Simulations are performed to validate our analytical results and illustrate the sum rate advantages of the proposed feedback scheme under a global feedback overhead constraint.


title: Experimental demonstration of the correction of coupled transverse dynamics aberration in an rf photoinjector

abstract: The production of electron bunches with low transverse emittance approaches the thermal emittance of the photocathode as various aberrations are corrected. Recently, the coupled transverse dynamics aberration was theoretically identified as a significant source of emittance growth and a corrector magnet was proposed for its elimination [D.H. Dowell, F. Zhou, and J. Schmerge, PRAB 21, 010101 (2018)]. This aberration arises when the beam acquires an asymmetric distribution that is then rotated with respect to the transverse reference axis thus introducing a correlation in the vertical and horizontal planes. The asymmetry is introduced by a weak quadrupole field in the rf gun or emittance compensation solenoid and the rotation is caused by the solenoid. This Letter presents an experimental study of the coupled transverse dynamics aberration in an rf photoinjector and demonstrates its elimination by a quadrupole corrector consisting of a normal and a skew quadrupole. The experimental results agree well with theoretical predictions and numerical simulations. The study also demonstrates the emittance of a low charge beam can be preserved during transportation at its thermal value, which was 1.05 mm mrad/mm, for the cesium telluride photocathode and 248 nm UV laser used.


title: On reconstruction from a partial knowledge of the Neumann-to-Dirichlet operator

abstract: We give formulae that yield an information about the location of an unknown polygonal inclusion having unknown constant conductivity inside a known conductive material having known constant conductivity from a partial knowledge of the Neumann -to-Dirichlet operator.


title: Liouville theorems on the upper half space

abstract: In this paper we shall establish some Liouville theorems for solutions bounded from below to certain linear elliptic equations on the upper half space. In particular, we show that for $a \in (0, 1)$ constants are the only $C^1$ up to the boundary positive solutions to $div(x_n^a \nabla u)=0$ on the upper half space.


title: A Comparison of Photometric Redshift Techniques for Large Radio Surveys

abstract: Future radio surveys will generate catalogues of tens of millions of radio sources, for which redshift estimates will be essential to achieve many of the science goals. However, spectroscopic data will be available for only a small fraction of these sources, and in most cases even the optical and infrared photometry will be of limited quality. Furthermore, radio sources tend to be at higher redshift than most optical sources and so a significant fraction of radio sources hosts differ from those for which most photometric redshift templates are designed. We therefore need to develop new techniques for estimating the redshifts of radio sources. As a starting point in this process, we evaluate a number of machine-learning techniques for estimating redshift, together with a conventional template-fitting technique. We pay special attention to how the performance is affected by the incompleteness of the training sample and by sparseness of the parameter space or by limited availability of ancillary multi-wavelength data. As expected, we find that the quality of the photometric-redshift degrades as the quality of the photometry decreases, but that even with the limited quality of photometry available for all sky-surveys, useful redshift information is available for the majority of sources, particularly at low redshift. We find that a template-fitting technique performs best with high-quality and almost complete multi-band photometry, especially if radio sources that are also X-ray emitting are treated separately. When we reduced the quality of photometry to match that available for the EMU all-sky radio survey, the quality of the template-fitting degraded and became comparable to some of the machine learning methods. Machine learning techniques currently perform better at low redshift than at high redshift, because of incompleteness of the currently available training data at high redshifts.


title: OPENMENDEL: A Cooperative Programming Project for Statistical Genetics

abstract: Statistical methods for genomewide association studies (GWAS) continue to improve. However, the increasing volume and variety of genetic and genomic data make computational speed and ease of data manipulation mandatory in future software. In our view, a collaborative effort of statistical geneticists is required to develop open source software targeted to genetic epidemiology. Our attempt to meet this need is called the OPENMENDELproject (https://openmendel.github.io). It aims to (1) enable interactive and reproducible analyses with informative intermediate results, (2) scale to big data analytics, (3) embrace parallel and distributed computing, (4) adapt to rapid hardware evolution, (5) allow cloud computing, (6) allow integration of varied genetic data types, and (7) foster easy communication between clinicians, geneticists, statisticians, and computer scientists. This article reviews and makes recommendations to the genetic epidemiology community in the context of the OPENMENDEL project.


title: On global solutions to a viscous compressible two-fluid model with unconstrained transition to single-phase flow in three dimensions

abstract: We consider the Dirichlet problem for a compressible two-fluid model in three dimensions, and obtain the global existence of weak solution with large initial data and independent adiabatic constants \Gamma,\gamma>=9/5. The pressure functions are of two components solving the continuity equations. Two typical cases for the pressure are considered, which are motivated by the compressible two-fluid model with possibly unequal velocities [3] and by a limiting system from the Vlasov-Fokker-Planck/compressible Navier-Stokes system [27] (see also some other relevant models like compressible MHD system for two-dimensional case [24] and compressible Oldroyd-B model with stress diffusion [1]). The lack of enough regularity for the two densities turns out some essential difficulties in the two-component pressure compared with the single-phase model, i.e., compressible Navier-Stokes equations. In this paper, the global existence theory does not require any domination conditions for the initial densities, which implies that transition to each single-phase flow is allowed.


title: Extraction formulae for an inverse boundary value problem for the equation $\nabla\cdot(\sigma-i\omega\epsilon)\nabla u=0$

abstract: We consider an inverse boundary value problem for the equation $\nabla\cdot(\sigma-i\omega\epsilon)\nabla u=0$ in a given bounded domain $\Omega$ at a fixed $\omega>0$. $\sigma$ and $\epsilon$ denote the conductivity and permittivity of the material forming $\Omega$, respectively. We give some formulae for extracting information about the location of the discontinuity surface of $(\sigma,\epsilon)$ from the Dirichlet-to-Neumann map. In order to obtain results we make use of two methods. The first is the enclosure method which is based on a new role of the exponentially growing solutions of the equation for the background material. The second is a generalization of the enclosure method based on a new role of Mittag-Leffler's function.


title: Potentially Information-theoretic Secure Y00 Quantum Stream Cipher with Limited Key Lengths beyond One-Time Pad

abstract: The previous work showed that the Y00 protocol could stay secure with the eavesdropper's guessing probability on the secret keys being strictly less than one under an unlimitedly long known-plaintext attack with quantum memory. However, an assumption that at least a fast correlation attack is completely disabled by irregular mapping. The present study shows that the Y00 protocol can be information-theoretic secure under any quantum-computational crypto-analyses if the Y00 system is well designed. The Y00 protocol directly encrypts messages with short secret keys expanded into pseudo-random running keys unlike One-Time Pad. However, it may offer information-theoretic security beyond the Shannon limit of cryptography.


title: An application of parallel cut elimination in multiplicative linear logic to the Taylor expansion of proof nets

abstract: We examine some combinatorial properties of parallel cut elimination in multiplicative linear logic (MLL) proof nets. We show that, provided we impose a constraint on some paths, we can bound the size of all the nets satisfying this constraint and reducing to a fixed resultant net. This result gives a sufficient condition for an infinite weighted sum of nets to reduce into another sum of nets, while keeping coefficients finite. We moreover show that our constraints are stable under reduction. Our approach is motivated by the quantitative semantics of linear logic: many models have been proposed, whose structure reflect the Taylor expansion of multiplicative exponential linear logic (MELL) proof nets into infinite sums of differential nets. In order to simulate one cut elimination step in MELL, it is necessary to reduce an arbitrary number of cuts in the differential nets of its Taylor expansion. It turns out our results apply to differential nets, because their cut elimination is essentially multiplicative. We moreover show that the set of differential nets that occur in the Taylor expansion of an MELL net automatically satisfies our constraints. Interestingly, our nets are untyped: we only rely on the sequentiality of linear logic nets and the dynamics of cut elimination. The paths on which we impose bounds are the switching paths involved in the Danos--Regnier criterion for sequentiality. In order to accommodate multiplicative units and weakenings, our nets come equipped with jumps: each weakening node is connected to some other node. Our constraint can then be summed up as a bound on both the length of switching paths, and the number of weakenings that jump to a common node.