title: $\beta$-CuI: a Dirac semimetal without surface Fermi arcs

abstract: Anomalous surface states with Fermi arcs are commonly considered to be a fingerprint of Dirac semimetals (DSMs). In contrast to Weyl semimetals, however, Fermi arcs of DSMs are not topologically protected. Using first-principles calculations, we predict that $\beta$-CuI is a peculiar DSM whose surface states form closed Fermi pockets instead of Fermi arcs. In such a fermiological Dirac semimetal, the deformation mechanism from Fermi arcs to Fermi pockets stems from a large cubic term preserving all crystal symmetries, and the small energy difference between the surface and bulk Dirac points. The cubic term in $\beta$-CuI, usually negligible in prototypical DSMs, becomes relevant because of the particular crystal structure. As such, we establish a concrete material example manifesting the lack of topological protection for surface Fermi arcs in DSMs


title: Towards a theoretical determination of the geographical probability distribution of meteoroid impacts on Earth

abstract: Tunguska and Chelyabinsk impact events occurred inside a geographical area of only 3.4\% of the Earth's surface. Although two events hardly constitute a statistically significant demonstration of a geographical pattern of impacts, their spatial coincidence is at least tantalizing. To understand if this concurrence reflects an underlying geographical and/or temporal pattern, we must aim at predicting the spatio-temporal distribution of meteoroid impacts on Earth. For this purpose we designed, implemented and tested a novel numerical technique, the "Gravitational Ray Tracing" (GRT) designed to compute the relative impact probability (RIP) on the surface of any planet. GRT is inspired by the so-called ray-casting techniques used to render realistic images of complex 3D scenes. In this paper we describe the method and the results of testing it at the time of large impact events. Our findings suggest a non-trivial pattern of impact probabilities at any given time on Earth. Locations at $60-90\deg$ from the apex are more prone to impacts, especially at midnight. Counterintuitively, sites close to apex direction have the lowest RIP, while in the antapex RIP are slightly larger than average. We present here preliminary maps of RIP at the time of Tunguska and Chelyabinsk events and found no evidence of a spatial or temporal pattern, suggesting that their coincidence was fortuitous. We apply the GRT method to compute theoretical RIP at the location and time of 394 large fireballs. Although the predicted spatio-temporal impact distribution matches marginally the observed events, we successfully predict their impact speed distribution.


title: A compressed classical description of quantum states

abstract: We show how to approximately represent a quantum state using the square root of the usual amount of classical memory. The classical representation of an $n$-qubit state $\psi$ consists of its inner products with $O(\sqrt{2^n})$ stabilizer states. A quantum state initially specified by its $2^n$ entries in the computational basis can be compressed to this form in time $O(2^n \mathrm{poly}(n))$, and, subsequently, the compressed description can be used to additively approximate the expectation value of an arbitrary observable. Our compression scheme directly gives a new protocol for the vector in subspace problem with randomized one-way communication complexity that matches (up to polylogarithmic factors) the optimal upper bound, due to Raz. We obtain an exponential improvement over Raz's protocol in terms of computational efficiency.


title: Heegaard Floer homology, degree-one maps and splicing knot complements

abstract: Let $K$ denote a knot inside the homology sphere $Y$ and $K'$ denote a knot inside a homology sphere $L$-space. Let $X=Y(K,K')$ denote the 3-manifold obtained by splicing the complements of $K$ and $K'$. We show that $\text{rank}(\widehat{HF}(X)) \ge \text{rank}(\widehat{HF}(Y))$.


title: Entanglement between a photonic time-bin qubit and a collective atomic spin excitation

abstract: Entanglement between light and matter combines the advantage of long distance transmission of photonic qubits with the storage and processing capabilities of atomic qubits. To distribute photonic states efficiently over long distances several schemes to encode qubits have been investigated -- time-bin encoding being particularly promising due to its robustness against decoherence in optical fibers. Here, we demonstrate the generation of entanglement between a photonic time-bin qubit and a single collective atomic spin excitation (spin-wave) in a cold atomic ensemble, followed by the mapping of the atomic qubit onto another photonic qubit. A magnetic field that induces a periodic dephasing and rephasing of the atomic excitation ensures the temporal distinguishability of the two time-bins and plays a central role in the entanglement generation. To analyse the generated quantum state, we use largely imbalanced Mach-Zehnder interferometers to perform projective measurements in different qubit bases and verify the entanglement by violating a CHSH Bell inequality.


title: On an extremal property of Jordan algebras of Clifford type

abstract: If $V$ is a finite-dimensional unital commutative (maybe nonassociative) algebra carrying an associative positive definite bilinear form then there exist a nonzero idempotent $c\ne e$ ($e$ being the algebra unit) of the shortest possible length $|c|^2$. In particular, $|c|^2\le \frac12|e|^2$. We prove that the equality holds exactly when $V$ is a Jordan algebra of Clifford type.


title: Bayesian Estimation of Gaussian Graphical Models with Predictive Covariance Selection

abstract: Gaussian graphical models are used for determining conditional relationships between variables. This is accomplished by identifying off-diagonal elements in the inverse-covariance matrix that are non-zero. When the ratio of variables (p) to observations (n) approaches one, the maximum likelihood estimator of the covariance matrix becomes unstable and requires shrinkage estimation. Whereas several classical (frequentist) methods have been introduced to address this issue, fully Bayesian methods remain relatively uncommon in practice and methodological literatures. Here we introduce a Bayesian method for estimating sparse matrices, in which conditional relationships are determined with projection predictive selection. With this method, that uses Kullback-Leibler divergence and cross-validation for neighborhood selection, we reconstruct the inverse-covariance matrix in both low and high-dimensional settings. Through simulation and applied examples, we characterized performance compared to several Bayesian methods and the graphical lasso, in addition to TIGER that similarly estimates the inverse-covariance matrix with regression. Our results demonstrate that projection predictive selection not only has superior performance compared to selecting the most probable model and Bayesian model averaging, particularly for high-dimensional data, but also compared to the the Bayesian and classical glasso methods. Further, we show that estimating the inverse-covariance matrix with multiple regression is often more accurate, with respect to various loss functions, and efficient than direct estimation. In low-dimensional settings, we demonstrate that projection predictive selection also provides competitive performance. We have implemented the projection predictive method for covariance selection in the R package GGMprojpred


title: A Markov Process Approach to the asymptotic Theory of abstract Cauchy Problems driven by Poisson Processes

abstract: In this paper, we employ Markov process theory to prove asymptotic results for a class of stochastic processes which arise as solutions of a stochastic evolution inclusion and are given by the representation formula \begin{align*} \mathbb{X}_{x}(t)=\sum \limits_{m=0}\limits^{\infty}T((t-\alpha_{m})_{+})(x_{x,m})1\hspace{-0,9ex}1_{[\alpha_{m},\alpha_{m+1})}(t), \end{align*} where $(T(t))_{t \geq 0}$ is a (nonlinear) time-continuous, contractive semigroup acting on a separable Banach space $(V,||\cdot||_{V})$, $(\alpha_{m})_{m \in \mathbb{N}}$ is the sequence of arrival times of a homogeneous Poisson process, $x$ is a $V$-valued random variable and $(x_{x,m})_{m \in \mathbb{N}}$ is a recursively defined sequence of $V$-valued random variables, fulfilling $x_{x,0}=x$. It will be demonstrated that $\mathbb{X}_{x}$ is, under some distributional assumptions on the involved random variables, a time-continuous Markov process and that it obeys, under polynomial decay conditions on $T$, the strong law of large numbers (SLLN) and, if the decay rate is sufficiently fast, also the central limit theorem (CLT). Finally, we consider two examples: A nonlinear ordinary differential equation and the (weighted) $p$-Laplacian evolution equation for $p \in (2,\infty)$.


title: CutLang: A Particle Physics Analysis Description Language and Runtime Interpreter

abstract: This note introduces CutLang, a domain specific language that aims to provide a clear, human readable way to define analyses in high energy particle physics (HEP) along with an interpretation framework of that language. A proof of principle (PoP) implementation of the CutLang interpreter, achieved using C++ as a layer over the CERN data analysis framework ROOT, is presently available. This PoP implementation permits writing HEP analyses in an unobfuscated manner, as a set of commands in human readable text files, which are interpreted by the framework at runtime. We describe the main features of CutLang and illustrate its usage with two analysis examples. Initial experience with CutLang has shown that a just-in-time interpretation of a human readable HEP specific language is a practical alternative to analysis writing using compiled languages such as C++.


title: Resurrecting the partially isotropic Haldane-Shastry model

abstract: We present a new and simpler expression for the Hamiltonian of the partially isotropic (XXZ-like) version of the Haldane-Shastry model, which was derived by D. Uglov over two decades ago in an apparently little-known preprint. While resembling the pairwise long-range form of the Haldane-Shastry model our formula accounts for the multi-spin interactions obtained by Uglov. Our expression is physically meaningful, makes hermiticity manifest, and is computationally more efficient. We discuss the model's properties, including its limits and (ordinary and quantum-affine) symmetries. In particular we introduce the appropriate notions of translational invariance and momentum. We review the model's exact spectrum found by Uglov for finite spin-chain length, which parallels the isotropic case up to level splitting due to the anisotropy. We also extend the partially isotropic model to higher rank, with $SU(n)$ 'spins', for which the spectrum is determined by $\mathfrak{sl}_n$-motifs.

categories: cond-mat.str-el hep-th math-ph math.MP math.QA nlin.SI


title: Weak Mixing in Switched Systems

abstract: Given a switched system, we introduce weakly mixing sets of type 1, 2 and Xiong chaotic sets of type 1, 2 with respect to a given set and show that they are equivalent respectively.


title: A Complete Analytical Study on the Dynamics of Simple Chaotic Systems

abstract: We report in this paper a complete analytical study on the bifurcations and chaotic phenomena observed in certain second-order, non-autonomous, dissipative chaotic systems. One-parameter bifurcation diagrams obtained from the analytical solutions proving the numerically observed chaotic phenomena such as antimonotonicity, period-doubling sequences, Feignbaum remerging have been presented. Further, the analytical solutions are used to obtain the basins of attraction, phase-portraits and Poincare maps for different chaotic systems. Experimentally observed chaotic attractors in some of the systems are presented to confirm the analytical results. The bifurcations and chaotic phenomena studied through explicit analytical solutions is reported in the literature for the first time.


title: In-network Neural Networks

abstract: We present N2Net, a system that implements binary neural networks using commodity switching chips deployed in network switches and routers. Our system shows that these devices can run simple neural network models, whose input is encoded in the network packets' header, at packet processing speeds (billions of packets per second). Furthermore, our experience highlights that switching chips could support even more complex models, provided that some minor and cheap modifications to the chip's design are applied. We believe N2Net provides an interesting building block for future end-to-end networked systems.


title: Homogeneous deformations of toric pairs

abstract: We extend the Altmann--Mavlyutov construction of homogeneous deformations of affine toric varieties to the case of toric pairs $(X, \partial X)$, where $X$ is an affine or projective toric variety and $\partial X$ is its toric boundary. As an application, we generalise a result due to Ilten to the case of Fano toric pairs.


title: Deep Convolutional Neural Networks for Eigenvalue Problems in Mechanics

abstract: We show that deep convolutional neural networks (CNN) can massively outperform traditional densely-connected neural networks (both deep or shallow) in predicting eigenvalue problems in mechanics. In this sense, we strike out in a new direction in mechanics computations with strongly predictive NNs whose success depends not only on architectures being deep, but also being fundamentally different from the widely-used to date. We consider a model problem: predicting the eigenvalues of 1-D and 2-D phononic crystals. For the 1-D case, the optimal CNN architecture reaches $98\%$ accuracy level on unseen data when trained with just 20,000 samples, compared to $85\%$ accuracy even with $100,000$ samples for the typical network of choice in mechanics research. We show that, with relatively high data-efficiency, CNNs have the capability to generalize well and automatically learn deep symmetry operations, easily extending to higher dimensions and our 2D case. Most importantly, we show how CNNs can naturally represent mechanical material tensors, with its convolution kernels serving as local receptive fields, which is a natural representation of mechanical response. Strategies proposed are applicable to other mechanics' problems and may, in the future, be used to sidestep cumbersome algorithms with purely data-driven approaches based upon modern deep architectures.


title: Eliminating the effect of rating bias on reputation systems

abstract: The ongoing rapid development of the e-commercial and interest-base websites make it more pressing to evaluate objects' accurate quality before recommendation by employing an effective reputation system. The objects' quality are often calculated based on their historical information, such as selected records or rating scores, to help visitors to make decisions before watching, reading or buying. Usually high quality products obtain a higher average ratings than low quality products regardless of rating biases or errors. However many empirical cases demonstrate that consumers may be misled by rating scores added by unreliable users or deliberate tampering. In this case, users' reputation, i.e., the ability to rating trustily and precisely, make a big difference during the evaluating process. Thus, one of the main challenges in designing reputation systems is eliminating the effects of users' rating bias on the evaluation results. To give an objective evaluation of each user's reputation and uncover an object's intrinsic quality, we propose an iterative balance (IB) method to correct users' rating biases. Experiments on two online video-provided Web sites, namely MovieLens and Netflix datasets, show that the IB method is a highly self-consistent and robust algorithm and it can accurately quantify movies' actual quality and users' stability of rating. Compared with existing methods, the IB method has higher ability to find the "dark horses", i.e., not so popular yet good movies, in the Academy Awards.

categories: physics.soc-ph cs.CY cs.IR cs.SI econ.TH


title: Simulations of Coulomb systems with slab geometry using an efficient 3d Ewald summation method

abstract: We present a new approach to efficiently simulate electrolytes confined between infinite charged walls using a 3d Ewald summation method. The optimal performance is achieved by separating the electrostatic potential produced by the charged walls from the electrostatic potential of electrolyte. The electric field produced by 3d periodic images of the plates is constant, with the field produced by the transverse images of the charged plates canceling out. We show that under suitable renormalization, the non-neutral electrolyte confined between charged plates can be simulated using 3d Ewald summation with a correction that accounts for the conditional convergence of the resulting lattice sum. The new algorithm is at least an order of magnitude more rapid than the usual simulation methods for the slab geometry and can be further sped up by adopting Particle-Particle Particle-Mesh (P 3 M ) approach.


title: Probing the BSM physics with CMB precision cosmology: an application to supersymmetry

abstract: The cosmic history before the BBN is highly determined by the physics that operates beyond the Standard Model (BSM) of particle physics and it is poorly constrained observationally. Ongoing and future precision measurements of the CMB observables can provide us with significant information about the pre-BBN era and hence possibly test the cosmological predictions of different BSM scenarios. Supersymmetry is a particularly motivated BSM theory and it is often the case that different superymmetry breaking schemes require different cosmic histories with specific reheating temperatures or low entropy production in order to be cosmologically viable. In this paper we quantify the effects of the possible alternative cosmic histories on the $n_s$ and $r$ CMB observables assuming a generic non-thermal stage after cosmic inflation. We analyze TeV and especially multi-TeV supersymmetry breaking schemes assuming the neutralino and gravitino dark matter scenarios. We complement our analysis considering the Starobinsky $R^2$ inflation model to exemplify the improved CMB predictions that a unified description of the early universe cosmic evolution yields. Our analysis underlines the importance of the CMB precision measurements that can be viewed, to some extend, as complementary to the laboratory experimental searches for supersymmetry or other BSM theories.


title: A fast and robust numerical scheme for solving models of charge carrier transport and ion vacancy motion in perovskite solar cells

abstract: Drift-diffusion models that account for the motion of both electronic and ionic charges are important tools for explaining the hysteretic behaviour and guiding the development of metal halide perovskite solar cells. Furnishing numerical solutions to such models for realistic operating conditions is challenging owing to the extreme values of some of the parameters. In particular, those characterising (i) the short Debye lengths (giving rise to rapid changes in the solutions across narrow layers), (ii) the relatively large potential differences across devices and (iii) the disparity in timescales between the motion of the electronic and ionic species give rise to significant stiffness. We present a finite difference scheme with an adaptive time step that is posed on a non-uniform staggered grid that provides second order accuracy in the mesh spacing. The method is able to cope with the stiffness of the system for realistic parameters values whilst providing high accuracy and maintaining modest computational costs. For example, a transient sweep of a current-voltage curve can be computed in only a few minutes on a standard desktop computer.


title: Characterization of Time Series Via R\'enyi Complexity-Entropy Curves

abstract: One of the most useful tools for distinguishing between chaotic and stochastic time series is the so-called complexity-entropy causality plane. This diagram involves two complexity measures: the Shannon entropy and the statistical complexity. Recently, this idea has been generalized by considering the Tsallis monoparametric generalization of the Shannon entropy, yielding complexity-entropy curves. These curves have proven to enhance the discrimination among different time series related to stochastic and chaotic processes of numerical and experimental nature. Here we further explore these complexity-entropy curves in the context of the R\'enyi entropy, which is another monoparametric generalization of the Shannon entropy. By combining the R\'enyi entropy with the proper generalization of the statistical complexity, we associate a parametric curve (the R\'enyi complexity-entropy curve) with a given time series. We explore this approach in a series of numerical and experimental applications, demonstrating the usefulness of this new technique for time series analysis. We show that the R\'enyi complexity-entropy curves enable the differentiation among time series of chaotic, stochastic, and periodic nature. In particular, time series of stochastic nature are associated with curves displaying positive curvature in a neighborhood of their initial points, whereas curves related to chaotic phenomena have a negative curvature; finally, periodic time series are represented by vertical straight lines.


title: Avoiding apparent signaling in Bell tests for quantitative applications

abstract: Bell tests have become a powerful tool for quantifying security, randomness, entanglement, and many other properties, as well as for investigating fundamental physical limits. In all these cases, the specific experimental value of the Bell parameter is important as it leads to a quantitative conclusion. However, most experimental implementations aiming for high values of the Bell parameter suffer from the defect of showing signaling. This signaling can be attributed to systematic errors occurring due to weaknesses in the experimental designs. Here we point out the importance, for quantitative applications, to identify and address this problem. We present a set of experiments with polarization-entangled photons in which we point out common sources of systematic errors and demonstrate approaches to avoid them. This allows us to establish a reliable estimate for the Bell parameter.


title: Effective sup-norm bounds on average for cusp forms of even weight

abstract: Let $\Gamma\subset\mathrm{PSL}_{2}(\mathbb{R})$ be a Fuchsian subgroup of the first kind acting on the upper half-plane $\mathbb{H}$. Consider the $d_{2k}$-dimensional space of cusp forms $\mathcal{S}_{2k}^{\Gamma}$ of weight $2k$ for $\Gamma$, and let $\{f_{1},\ldots,f_{d_{2k}}\}$ be an orthonormal basis of $\mathcal{S}_{2k}^{\Gamma}$ with respect to the Petersson inner product. In this paper we will give effective upper and lower bounds for the supremum of the quantity $S_{2k}^{\Gamma}(z):=\sum_{j=1}^{d_{2k}}\vert f_{j}(z)\vert^{2}\,\mathrm{Im}(z)^{2k}$ as $z$ ranges through $\mathbb{H}$.


title: StreetGen : In base city scale procedural generation of streets: road network, road surface and street objects

abstract: Streets are large, diverse, and used for several (and possibly conflicting) transport modalities as well as social and cultural activities. Proper planning is essential and requires data. Manually fabricating data that represent streets (street reconstruction) is error-prone and time consuming. Automatising street reconstruction is a challenge because of the diversity, size, and scale of the details (few centimetres for cornerstone) required. The state-of-the-art focuses on roads (no context, no urban features) and is strongly determined by each application (simulation, visualisation, planning). We propose a unified framework that works on real Geographic Information System (GIS) data and uses a strong, yet simple modelling hypothesis when possible to robustly model streets at the city level or street level. Our method produces a coherent street-network model containing topological traffic information, road surface and street objects. We demonstrate the robustness and genericity of our method by reconstructing the entire city of Paris streets and exploring other similar reconstruction (airport driveway).


title: Flux cost functions and optimal metabolic states

abstract: The metabolic fluxes in cells follow physical, biochemical, and economic principles. Some flux balance analysis (FBA) methods trade flux benefit against flux cost. However, if flux cost functions are linear and meant to describe underlying enzyme costs, this entails that enzyme efficiencies are constant and ignores the interplay between fluxes, metabolite concentrations and enzyme levels in cells. Here I introduce realistic flux cost functions that describe an "overhead cost", namely the minimum enzyme and metabolite cost associated with the fluxes in a kinetic model. These flux cost functions have general mathematical properties. Enzymatic flux cost functions, which represent enzyme costs, scale proportionally with the flux profile and are concave on the flux polytope. Kinetic flux cost functions represent the sum of enzyme and metabolite costs. If two flux profiles are superimposed, their different demands for metabolite concentrations cause an extra compromise cost, which makes flux cost functions strictly concave in almost all cases. When fluxes change their direction, the enzymatic cost jumps abruptly. Here I propose two methods for flux modelling: Flux Cost Minimisation, a nonlinear variant of FBA with flux minimisation, and Flux Benefit Optimisation, a nonlinear variant of FBA with molecular crowding. The optimal flux profiles, at a given flux benefit, are vertices of the flux polytope. Linear approximations of enzymatic flux cost can be used in FBA. In contrast to flux costs chosen ad hoc, these functions reflect the enzyme kinetics and extracellular concentrations in realistic kinetic models. Based on enzymatic flux costs, we can describe the cell growth rate as a convex function on the flux polytope and derive growth-optimal metabolic states and statistical distributions for the fluxes in cell populations.


title: The mechanical influence of single fibre inclusions in discrete fibre networks

abstract: Semiflexible biopolymer networks are commonly found in biological systems, from the cytoskeleton of cells to the extracellular matrix. Such networks often naturally occur as composites, in which various components interact to generate rich mechanical behaviours. In this work we examine the mechanics of composites formed when a single fibre inclusion is placed within discrete fibre networks of two distinct architectures. In particular, we computationally and theoretically investigate the mechanics of composites formed when an inclusion is introduced to Voronoi- and Mikado--type networks within the nonaffine regime. On subjecting these single inclusion composites to small shear deformations, we observe different behaviours dependent on the choice of network geometry. This divergence in mechanical responses is interpreted as a consequence of architecture-dependent differences in the nonaffine displacement field. Scaling laws for the increase in network energy due to inclusions of different lengths, orientations, and elastic properties are proposed and computationally verified. We present theoretical predictions for critical values of the inclusion bending and stretching stiffnesses, above which no further increases in network energy are observed for a given shear deformation. These predictions are supported by extensive computational evidence. We expect the architectural differences identified in this work to be pertinent to theoretical investigations into the complex behaviour of biopolymer networks, and to experimental work, where complex mechanics has been observed in composite fibre networks.


title: Forecasting Solar Flares Using Magnetogram-based Predictors and Machine Learning

abstract: We propose a forecasting approach for solar flares based on data from Solar Cycle 24, taken by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) mission. In particular, we use the Space-weather HMI Active Region Patches (SHARP) product that facilitates cut-out magnetograms of solar active regions (AR) in the Sun in near-real-time (NRT), taken over a five-year interval (2012 - 2016). Our approach utilizes a set of thirteen predictors, which are not included in the SHARP metadata, extracted from line-of-sight and vector photospheric magnetograms. We exploit several Machine Learning (ML) and Conventional Statistics techniques to predict flares of peak magnitude >M1 and >C1, within a 24 h forecast window. The ML methods used are multi-layer perceptrons (MLP), support vector machines (SVM) and random forests (RF). We conclude that random forests could be the prediction technique of choice for our sample, with the second best method being multi-layer perceptrons, subject to an entropy objective function. A Monte Carlo simulation showed that the best performing method gives accuracy ACC=0.93(0.00), true skill statistic TSS=0.74(0.02) and Heidke skill score HSS=0.49(0.01) for >M1 flare prediction with probability threshold 15% and ACC=0.84(0.00), TSS=0.60(0.01) and HSS=0.59(0.01) for >C1 flare prediction with probability threshold 35%.


title: Fast Generation of Covariance Matrices for Weak Lensing

abstract: Upcoming weak lensing surveys will probe large fractions of the sky with unprecedented accuracy. To infer cosmological constraints, a large ensemble of survey simulations are required to accurately model cosmological observables and their covariances. We develop a parallelized multi-lens-plane pipeline called UFalcon, designed to generate full-sky weak lensing maps from lightcones within a minimal runtime. It makes use of L-PICOLA, an approximate numerical code, which provides a fast and accurate alternative to cosmological $N$-Body simulations. The UFalcon maps are constructed by nesting 2 simulations covering a redshift-range from $z=0.1$ to $1.5$ without replicating the simulation volume. We compute the convergence and projected overdensity maps for L-PICOLA in the lightcone or snapshot mode. The generation of such a map, including the L-PICOLA simulation, takes about 3 hours walltime on 220 cores. We use the maps to calculate the spherical harmonic power spectra, which we compare to theoretical predictions and to UFalcon results generated using the full $N$-Body code GADGET-2. We then compute the covariance matrix of the full-sky spherical harmonic power spectra using 150 UFalcon maps based on L-PICOLA in lightcone mode. We consider the PDF, the higher-order moments and the variance of the smoothed field variance to quantify the accuracy of the covariance matrix, which we find to be a few percent for scales $\ell \sim 10^2$ to $10^3$. We test the impact of this level of accuracy on cosmological constraints using an optimistic survey configuration, and find that the final results are robust to this level of uncertainty. The speed and accuracy of our developed pipeline provides a basis to also include further important features such as masking, varying noise and will allow us to compute covariance matrices for models beyond $\Lambda$CDM. [abridged]


title: TernausNet: U-Net with VGG11 Encoder Pre-Trained on ImageNet for Image Segmentation

abstract: Pixel-wise image segmentation is demanding task in computer vision. Classical U-Net architectures composed of encoders and decoders are very popular for segmentation of medical images, satellite images etc. Typically, neural network initialized with weights from a network pre-trained on a large data set like ImageNet shows better performance than those trained from scratch on a small dataset. In some practical applications, particularly in medicine and traffic safety, the accuracy of the models is of utmost importance. In this paper, we demonstrate how the U-Net type architecture can be improved by the use of the pre-trained encoder. Our code and corresponding pre-trained weights are publicly available at https://github.com/ternaus/TernausNet. We compare three weight initialization schemes: LeCun uniform, the encoder with weights from VGG11 and full network trained on the Carvana dataset. This network architecture was a part of the winning solution (1st out of 735) in the Kaggle: Carvana Image Masking Challenge.


title: The simultaneous discharge of liquid and grains from a silo

abstract: The flow rate of water through an orifice at the bottom of a container depends on the hydrostatic pressure whereas for a dry granular material is nearly constant. But what happens during the simultaneous discharge of grains and liquid from a silo? By measuring the flow rate as a function of time, we found that: (i) different regimes appear, going from constant flow rate dominated by the effective fluid viscosity to a hydrostatic-like discharge depending on the aperture and grain size, (ii) the mixed material is always discharged faster than dry grains but slower than liquid, (iii) for the mixture, the liquid level drops faster than the grains level but they are always linearly proportional to one another, and (iv) a sudden growth in the flow rate happens during the transition from a biphasic to a single phase discharge. These results are associated to the competition between the decrease of hydrostatic pressure above the granular bed and the hydrodynamic resistance. A model combining the Kozeny-Carman, Bernoulli and mass conservation equations is proposed and the numerical results are in good agreement with experiments.


title: Supersimple omega-categorical theories and pregeometries

abstract: We prove that if $T$ is an $\omega$-categorical supersimple theory with nontrivial dependence (given by forking), then there is a nontrivial regular 1-type over a finite set of reals which is realized by real elements; hence forking induces a nontrivial pregeometry on the solution set of this type and the pregeometry is definable (using only finitely many parameters). The assumption about $\omega$-categoricity is necessary. This result is used to prove the following: If $V$ is a finite relational vocabulary with maximal arity 3 and $T$ is a supersimple $V$-theory with elimination of quantifiers, then $T$ has trivial dependence and finite SU-rank. This immediately gives the following strengthening of a previous result of the author: if $\mathcal{M}$ is a ternary simple homogeneous structure with only finitely many constraints, then $Th(\mathcal{M})$ has trivial dependence and finite SU-rank.


title: On Gaussian random matrices coupled to the discrete Laplacian

abstract: We study operators obtained by coupling an $n \times n$ random matrix from one of the Gaussian ensembles to the discrete Laplacian. We find the joint distribution of the eigenvalues and resonances of such operators. This is one of the possible mathematical models for quantum scattering in a complex physical system with one semi-infinite lead attached.


title: No exceptional precision of exceptional point sensors

abstract: Recently, sensors with resonances at exceptional points (EPs) have been suggested to have a vastly improved sensitivity due to the extraordinary scaling of the complex frequency splitting of the $n$ initially degenerate modes with the $n$-th root of the perturbation. We show here that the resulting quantum-limited signal to noise at EPs is proportional to the perturbation, and comparable to other sensors, thus providing the same precision. The complex frequency splitting close to EPs is therefore not suited to estimate the precision of EP sensors. The underlying reason of this counter-intuitive result is that the mode fields, described by the eigenvectors, are equal for all modes at the EP, and are strongly changing with the perturbation.


title: On the equidistribution of some Hodge loci

abstract: We prove the equidistribution of the Hodge locus for certain non-isotrivial, polarized variations of Hodge structure of weight $2$ with $h^{2,0}=1$ over complex, quasi-projective curves. Given some norm condition, we also give an asymptotic on the growth of the Hodge locus. In particular, this implies the equidistribution of elliptic fibrations in quasi-polarized, non-isotrivial families of $K3$ surfaces.


title: Part 1: Training Sets & ASG Transforms

abstract: In this paper, I discuss a method to tackle the issues arising from the small data-sets available to data-scientists when building price predictive algorithms that use monthly/quarterly macro-financial indicators. I approach this by training separate classifiers on the equivalent dataset from a range of countries. Using these classifiers, a three level meta learning algorithm (MLA) is developed. I develop a transform, ASG, to create a country agnostic proxy for the macro-financial indicators. Using these proposed methods, I investigate the degree to which a predictive algorithm for the US 5Y bond price, predominantly using macro-financial indicators, can outperform an identical algorithm which only uses statistics deriving from previous price. This was an undergraduate project, subsequently the research was not exhaustive.


title: The Lipman-Zariski conjecture in low genus

abstract: We prove the Lipman-Zariski conjecture for complex surface singularities of genus one, and also for those of genus two whose link is not a rational homology sphere. As an application, we characterize complex $2$-tori as the only normal compact complex surfaces whose smooth locus has trivial tangent bundle. We also deduce that all complex-projective surfaces with locally free and generically nef tangent sheaf are smooth, and we classify them.


title: Experiments Indicate Regolith is Looser in the Lunar Polar Regions than at the Lunar Landing Sites

abstract: Since the Apollo program or earlier it has been widely believed that the lunar regolith was compacted through vibrations including nearby impact events, thermal stress release in the regolith, deep moon quakes, and shallow moon quakes. Experiments have shown that vibrations both compact and re-loosen regolith as a function of depth in the lunar soil column and amplitude of the vibrational acceleration. Experiments have also identified another process that is extremely effective at compacting regolith: the expansion and contraction of individual regolith grains due to thermal cycling in the upper part of the regolith where the diurnal thermal wave exists. Remote sensing data sets from the Moon suggest that the soil is less compacted in regions where there is less thermal cycling, including infrared emissions measured by the Diviner radiometer on the Lunar Reconnaissance Orbiter (LRO). Here, we performed additional experiments in thermal cycling simulated lunar regolith and confirm that it is an effective compaction mechanism and may explain the remote sensing data. This creates a consistent picture that the soil really is looser in the upper layers in polar regions, which may be a challenge for rovers that must drive in the looser soil.


title: Discussions on non-probabilistic convex modelling for uncertain problems

abstract: Non-probabilistic convex model utilizes a convex set to quantify the uncertainty domain of uncertain-but-bounded parameters, which is very effective for structural uncertainty analysis with limited or poor-quality experimental data. To overcome the complexity and diversity of the formulations of current convex models, in this paper, a unified framework for construction of the non-probabilistic convex models is proposed. By introducing the correlation analysis technique, the mathematical expression of a convex model can be conveniently formulated once the correlation matrix of the uncertain parameters is created. More importantly, from the theoretic analysis level, an evaluation criterion for convex modelling methods is proposed, which can be regarded as a test standard for validity verification of subsequent newly proposed convex modelling methods. And from the practical application level, two model assessment indexes are proposed, by which the adaptabilities of different convex models to a specific uncertain problem with given experimental samples can be estimated. Four numerical examples are investigated to demonstrate the effectiveness of the present study.


title: Content Placement in Cache-Enabled Sub-6 GHz and Millimeter-Wave Multi-antenna Dense Small Cell Networks

abstract: This paper studies the performance of cache-enabled dense small cell networks consisting of multi-antenna sub-6 GHz and millimeter-wave base stations. Different from the existing works which only consider a single antenna at each base station, the optimal content placement is unknown when the base stations have multiple antennas. We first derive the successful content delivery probability by accounting for the key channel features at sub-6 GHz and mmWave frequencies. The maximization of the successful content delivery probability is a challenging problem. To tackle it, we first propose a constrained cross-entropy algorithm which achieves the near-optimal solution with moderate complexity. We then develop another simple yet effective heuristic probabilistic content placement scheme, termed two-stair algorithm, which strikes a balance between caching the most popular contents and achieving content diversity. Numerical results demonstrate the superior performance of the constrained cross-entropy method and that the two-stair algorithm yields significantly better performance than only caching the most popular contents. The comparisons between the sub-6 GHz and mmWave systems reveal an interesting tradeoff between caching capacity and density for the mmWave system to achieve similar performance as the sub-6 GHz system.


title: Experience-driven Networking: A Deep Reinforcement Learning based Approach

abstract: Modern communication networks have become very complicated and highly dynamic, which makes them hard to model, predict and control. In this paper, we develop a novel experience-driven approach that can learn to well control a communication network from its own experience rather than an accurate mathematical model, just as a human learns a new skill (such as driving, swimming, etc). Specifically, we, for the first time, propose to leverage emerging Deep Reinforcement Learning (DRL) for enabling model-free control in communication networks; and present a novel and highly effective DRL-based control framework, DRL-TE, for a fundamental networking problem: Traffic Engineering (TE). The proposed framework maximizes a widely-used utility function by jointly learning network environment and its dynamics, and making decisions under the guidance of powerful Deep Neural Networks (DNNs). We propose two new techniques, TE-aware exploration and actor-critic-based prioritized experience replay, to optimize the general DRL framework particularly for TE. To validate and evaluate the proposed framework, we implemented it in ns-3, and tested it comprehensively with both representative and randomly generated network topologies. Extensive packet-level simulation results show that 1) compared to several widely-used baseline methods, DRL-TE significantly reduces end-to-end delay and consistently improves the network utility, while offering better or comparable throughput; 2) DRL-TE is robust to network changes; and 3) DRL-TE consistently outperforms a state-ofthe-art DRL method (for continuous control), Deep Deterministic Policy Gradient (DDPG), which, however, does not offer satisfying performance.


title: On Grothendieck's section conjecture for orbicurves

abstract: As already noted by Niels Borne and Michel Emsalem, there is a natural generalization of the section conjecture for proper orbicurves. Combined with the reformulation by Niels Borne and Angelo Vistoli of the conjecture in terms of the \'etale fundamental gerbe, this suggests an even stronger conjecture for orbicurves, asking an equivalence of categories instead of a mere bijection. We prove that the three versions of the conjecture are in fact equivalent, and that "injectivity" (i.e. full faithfulness) holds in the case of orbicurves. As a byproduct, we obtain a new proof of the fact that the section conjecture for proper curves implies the section conjecture for open curves.


title: Evaluating the role of risk networks on risk identification, classification and emergence

abstract: Modern society heavily relies on strongly connected, socio-technical systems. As a result, distinct risks threatening the operation of individual systems can no longer be treated in isolation. Consequently, risk experts are actively seeking for ways to relax the risk independence assumption that undermines typical risk management models. Prominent work has advocated the use of risk networks as a way forward. Yet, the inevitable biases introduced during the generation of these survey-based risk networks limit our ability to examine their topology, and in turn challenge the utility of the very notion of a risk network. To alleviate these concerns, we proposed an alternative methodology for generating weighted risk networks. We subsequently applied this methodology to an empirical dataset of financial data. This paper reports our findings on the study of the topology of the resulting risk network. We observed a modular topology, and reasoned on its use as a robust risk classification framework. Using these modules, we highlight a tendency of specialization during the risk identification process, with some firms being solely focused on a subset of the available risk classes. Finally, we considered the independent and systemic impact of some risks and attributed possible mismatches to their emerging nature.


title: CryptoRuble: From Russia with Love

abstract: We discuss Russia's underlying motives for issuing its government-backed cryptocurrency, CryptoRuble, and the implications thereof and of other likely-soon-forthcoming government-issued cryptocurrencies to some stakeholders (populace, governments, economy, finance, etc.), existing decentralized cryptocurrencies (such as Bitcoin and Ethereum), as well as the future of the world monetary system (the role of the U.S. therein and a necessity for the U.S. to issue CryptoDollar), including a future algorithmic universal world currency that may also emerge. We further provide a comprehensive list of references on cryptocurrencies.


title: Rasch Analysis of the Mathematics Self Concept Questionnaire

abstract: TIMSS & PIRLS International Study Center is a research center at Boston College that conducts a series of assessments in a number of countries to measure trends in mathematics and science achievement at the fourth and eighth grades. In general, TIMSS assessments include achievement tests as well as questionnaires for student, parent, teacher, school, and curricular. There are 63 participating countries and 14 benchmarking participants in TIMSS 2011, including 608,641 students, 49,429 teachers, 19,612 school principals, and the National Research Coordinators of each country. TIMSS data are valuable for researchers and analysts from all over the world, especially for those from participating countries, to conduct related studies to improve education.


title: Heights in families of abelian varieties and the Geometric Bogomolov Conjecture

abstract: On an abelian scheme over a smooth curve over $\overline{\mathbb Q}$ a symmetric relatively ample line bundle defines a fiberwise N\'eon-Tate height. If the base curve is inside a projective space, we also have a height on its $\overline{\mathbb Q}$-points that serves as a measure of each fiber, an abelian variety. Silverman proved an asymptotic equality between these two heights on a curve in the abelian scheme. In this paper we prove an inequality between these heights on a subvariety of any dimension of the abelian scheme. As an application we prove the Geometric Bogomolov Conjecture for the function field of a curve defined over $\overline{\mathbb Q}$. Using Moriwaki's height we sketch how to extend our result when the base field of the curve has characteristic 0.


title: Cavity solitons in a microring dimer with gain and loss

abstract: We address a pair of vertically coupled microring resonators with gain and loss pumped by a single-frequency field. Coupling between microrings results in a twofold splitting of the single microring resonance that increases when gain and losses decrease and that gives rise to two different cavity soliton (CS) families. We show that the existence regions of CSs are tunable and that both CS families can be stable in the presence of an imbalance between gain and losses in the two microrings. These findings enable experimental realization of frequency combs in configurations with active microrings and contribute towards the realization of compact multisoliton comb sources.


title: M-STAR: A Modular, Evidence-based Software Trustworthiness Framework

abstract: Despite years of intensive research in the field of software vulnerabilities discovery, exploits are becoming ever more common. Consequently, it is more necessary than ever to choose software configurations that minimize systems' exposure surface to these threats. In order to support users in assessing the security risks induced by their software configurations and in making informed decisions, we introduce M-STAR, a Modular Software Trustworthiness ARchitecture and framework for probabilistically assessing the trustworthiness of software systems, based on evidence, such as their vulnerability history and source code properties. Integral to M-STAR is a software trustworthiness model, consistent with the concept of computational trust. Computational trust models are rooted in Bayesian probability and Dempster-Shafer Belief theory, offering mathematical soundness and expressiveness to our framework. To evaluate our framework, we instantiate M-STAR for Debian Linux packages, and investigate real-world deployment scenarios. In our experiments with real-world data, M-STAR could assess the relative trustworthiness of complete software configurations with an error of less than 10%. Due to its modular design, our proposed framework is agile, as it can incorporate future advances in the field of code analysis and vulnerability prediction. Our results point out that M-STAR can be a valuable tool for system administrators, regular users and developers, helping them assess and manage risks associated with their software configurations.


title: Cosmological constraints from low-redshift data

abstract: In this paper we summarise the constraints that low-redshift data --such as supernovae Ia (SN Ia), baryon acoustic oscillations (BAO) and cosmic chronometers (CC)-- are able to set on the concordance model and its extensions, as well as on inhomogeneous but isotropic models. We provide a broad overlook into these cosmological scenarios and several aspects of data analysis. In particular, we review a number of systematic issues of SN Ia analysis that include magnitude correction techniques, selection bias and their influence on the inferred cosmological constraints. Furthermore, we examine the isotropic and anisotropic components of the BAO data and their individual relevance for cosmological model-fitting. We extend the discussion presented in earlier works regarding the inferred dynamics of cosmic expansion and its present rate from the low-redshift data. Specifically, we discuss the cosmological constraints on the accelerated expansion and related model-selections. In addition, we extensively talk about the Hubble constant problem, then focus on the low-redshift data constraint on $H_0$ that is based on CC. Finally, we present the way in which this result compares the two of high-redshift $H_0$ estimate and local (redshift zero) measurements that are in tension.


title: A model of dynamic stability of H3K9me3 heterochromatin to explain the resistance to reprogramming of differentiated cells

abstract: Despite their dynamic nature, certain chromatin marks must be maintained over the long term. This is particulary true for histone 3 lysine 9 (H3K9) trimethylation, that is involved in the maintenance of healthy differentiated cellular states by preventing inappropriate gene expression, and has been recently identified as the most efficient barrier to cellular reprogramming in nuclear transfer experiments. We propose that the capacity of the enzymes SUV39H1/2 to rebind to a minor fraction of their products, either directly or via HP1, contributes to the solidity of this mark through (i) a positive feedback involved in its establishment by the mutual enforcement of H3K9me3 and SUV39H1/2 and then (ii) a negative feedback sufficient to strongly stabilize H3K9me3 heterochromatin in post-mitotic cells by generating local enzyme concentrations capable of counteracting transient bursts of demethylation. This model does not require direct molecular interactions with adjacent nucleosomes and is favoured by a series of additional mechanisms including (i) the protection of chromatin-bound SUV39H1/2 from the turnovers of soluble proteins, which can explain the uncoupling between the cellular contents in SUV39H1 mRNA and protein; (ii) the cooperative dependence on the local density of the H3K9me3 of HP1-dependent heterochomatin condensation and, dispensably (iii) restricted enzyme exchanges with chromocenters confining the reactive bursts of SUV39H1/2 in heterochromatin. This mechanism illustrates how seemingly static epigenetic states can be firmly maintained by dynamic and reversible modifications.


title: Polylogarithmic equilibrium treatment of molecular aggregation and critical concentrations

abstract: A full equilibrium treatment of molecular aggregation is presented for prototypes of 1D and 3D aggregates, with and without nucleation. By skipping complex kinetic parameters like aggregate size-dependent diffusion, the equilibrium treatment allows to predict directly time-independent quantities such as critical concentrations. The relationships between the macroscopic equilibrium constants for the different paths are first established by statistical corrections and so as to comply with the detailed balance constraints imposed by nucleation, and the composition of the mixture resulting from homogeneous aggregation is then analyzed using the polylogarithm function. Several critical concentrations are distinguished: the residual monomer concentation in equilibrium (RMC) and the critical nucleation concentration (CNC), that is the threshold concentration of total subunits necessary for initiating aggregation. When increasing the concentration of total subunits, the RMC converges more strongly to its asymptotic value, the equilibrium constant of depolymerization, for 3D aggregates and in case of nucleation. The CNC moderately depends on the number of subunits in the nucleus, but sharply increases with the difference between the equilibrium constants of polymerization and nucleation. As the RMC and CNC can be numerically but not analytically determined, ansatz equations connecting them to thermodynamic parameters are proposed.


title: The Asymptotic Capacity of Private Search

abstract: The private search problem is introduced, where a dataset comprised of $L$ i.i.d. records is replicated across $N$ non-colluding servers, each record takes values uniformly from an alphabet of size $K$, and a user wishes to search for all records that match a privately chosen value, without revealing any information about the chosen value to any individual server. The capacity of private search is the maximum number of bits of desired information that can be retrieved per bit of download. The asymptotic (large $K$) capacity of private search is shown to be $1-1/N$, even as the scope of private search is further generalized to allow approximate (OR) search over a number of realizations that grows with $K$. The results are based on the asymptotic behavior of a new converse bound for private information retrieval with arbitrarily dependent messages.


title: Single crystal growth of Cu4(OH)6BrF and universal behavior in quantum spin liquid candidates synthetic barlowite and herbertsmithite

abstract: Synthetic barlowite, Cu4(OH)6BrF, has emerged as a new quantum spin liquid (QSL) host, containing kagome layers of S=1/2 Cu2+ ions separated by interlayer Cu2+ ions. Similar to synthetic herbertsmithite, ZnCu3(OH)6Cl2, it has been reported that Zn2+ substitution for the interlayer Cu2+ induces a QSL ground state. Here we report a scalable synthesis of single crystals of Cu4(OH)6BrF. Through x-ray, neutron, and electron diffraction measurements coupled with magic angle spinning 19F and 1H NMR spectroscopy, we resolve the previously reported positional disorder of the interlayer Cu2+ ions and find that the structure is best described in the orthorhombic space group, Cmcm, with lattice parameters a = 6.665(13) A, b = 11.521(2) A, c = 9.256(18) A and an ordered arrangement of interlayer Cu2+ ions. Infrared spectroscopy measurements of the O-H and F-H stretching frequencies demonstrate that the orthorhombic symmetry persists upon substitution of Zn2+ for Cu2+. Specific heat and magnetic susceptibility measurements of Zn-substituted barlowite, ZnxCu4-x(OH)6BrF, reveal striking similarities with the behavior of ZnxCu4-x(OH)6Cl2. These parallels imply universal behavior of copper kagome lattices even in the presence of small symmetry breaking distortions. Thus synthetic barlowite demonstrates universality of the physics of synthetic Cu2+ kagome minerals and furthers the development of real QSL states.


title: The macroeconomics determinants of default of the borrowers: The case of Moroccan bank

abstract: This article aims to explore an empirical approach to analyze the macroeconomicsdeterminants of default of borrowers. For this purpose, we have measured the impact of the adverse economic conditions on the degradation of the credit portfolio quality.In our paper, we have shed more light on the question of the aggravation of default rate. For this, we have undertaken econometric modeling of the default rate distribution of a Moroccan bank while we inspired from some studies carried out. Our findings demonstrate that the decline in the economic situation has a positive impact on default of borrowers. Hence, the bank also has responsibility for monitoring the adverse economic conditions.


title: Propagation velocity of slip front and emergence of macroscopic static friction in the system with vanishing local static friction

abstract: We investigate the propagation of the slip front in the elastic body on the rigid substrate. We first obtain the slip profile and the slip front velocity of the steady state by employing the local friction law with the quadratic form of the slip velocity and with vanishing static friction stress. The macroscopic static friction stress emerges spontaneously, which is expressed in terms of the parameter emerging in the friction law. For the model with viscosity, the macroscopic static friction stress again emerges spontaneously. The analytical treatment gives estimations for two slip front propagation velocities. They corresponds to two different boundary conditions, and one of them describes the framework employed here. Linear Marginal Stability Hypothesis based on the linearized equation of motion shows that two slip front propagation velocities exist in this system, both of which coincide with the analytical solutions noted above. These imply that the linearized friction law dominantly governs the slip front propagation behavior. Seismological implications are also given based on the analytical and numerical results.


title: Ranking Data with Continuous Labels through Oriented Recursive Partitions

abstract: We formulate a supervised learning problem, referred to as continuous ranking, where a continuous real-valued label Y is assigned to an observable r.v. X taking its values in a feature space $\mathcal{X}$ and the goal is to order all possible observations x in $\mathcal{X}$ by means of a scoring function $s:\mathcal{X}\rightarrow \mathbb{R}$ so that s(X) and Y tend to increase or decrease together with highest probability. This problem generalizes bi/multi-partite ranking to a certain extent and the task of finding optimal scoring functions s(x) can be naturally cast as optimization of a dedicated functional criterion, called the IROC curve here, or as maximization of the Kendall ${\tau}$ related to the pair (s(X), Y ). From the theoretical side, we describe the optimal elements of this problem and provide statistical guarantees for empirical Kendall ${\tau}$ maximization under appropriate conditions for the class of scoring function candidates. We also propose a recursive statistical learning algorithm tailored to empirical IROC curve optimization and producing a piecewise constant scoring function that is fully described by an oriented binary tree. Preliminary numerical experiments highlight the difference in nature between regression and continuous ranking and provide strong empirical evidence of the performance of empirical optimizers of the criteria proposed.


title: Radiative Transfer: Asymptotic Solutions of the Kinetic Equation of Radiation Propagation, $n$th Order Asymptotic Approximation and Improved Boundary Conditions

abstract: In the article, new asymptotic approximation of the $n$th order is obtained and proposed to be used in calculations of radiation propagation without scattering in optically thick media; the asymptotic approximation is much simpler and more precise than the known diffusion approximation. The rigorous derivation of the diffusion approximation equation and the equation of the radiation heat conduction approximation is obtained from the constructed asymptotic solution of the kinetic equation of radiation propagation in optically thick media. It is shown, that for optically thick media the asymptotic solution of the kinetic equation of radiation propagation without scattering is asymptotic expansion of the exact integral solution of that kinetic equation. Improved boundary conditions, which are essential for practical application in calculations of radiation propagation, are derived (for inner boundaries and outer boundaries with vacuum).


title: Product of three octonions

abstract: This paper is devoted to octonions that are the eight-dimensional hypercomplex numbers characterized by multiplicative non-associativity. The decomposition of the product of three octonions with the conjugated central factor into the sum of mutually orthogonal anticommutator, commutator and associator, is introduced in an obvious way by commuting of factors and alternating the multiplication order. The commutator is regarded as a generalization of the cross product to the case of three arguments both for quaternions and for octonions. It is verified that the resulting additive decomposition is equivalent to the known solution derived and presented by S. Okubo in a cumbersome form.


title: Local Coordinate Spaces: a proposed unification of manifolds and fiber bundles, and associated machinery

abstract: This paper presents a unified view of manifolds and fiber bundles, which, while superficially different, have strong parallels. It introduces the notions of an m-atlas and of a local coordinate space, and shows that special cases are equivalent to fiber bundles and manifolds. Along the way it defines some convenient notation, defines categories of atlases, and constructs potentially useful functors.


title: Metastable transitions in inertial Langevin systems: what can be different from the overdamped case?

abstract: Metastable transitions in Langevin dynamics can exhibit rich behaviors that are markedly different from its overdamped limit. In addition to local alterations of the transition path geometry, more fundamental global changes may exist. For instance, when the dissipation is weak, heteroclinic connections that exist in the overdamped limit do not necessarily have a counterpart in the Langevin system, potentially leading to different transition rates. Furthermore, when the friction coefficient depends on the velocity, the overdamped limit no longer exists, but it is still possible to efficiently find instantons. The approach we employed for these discoveries was based on (i) a simple rewriting of the Freidlin-Wentzell action in terms of time-reversed dynamics, and (ii) an adaptation of the string method, which was originally designed for gradient systems, to this specific non-gradient system.

categories: physics.comp-ph cond-mat.stat-mech math.DS math.NA math.PR


title: On the number of maximal paths in directed last-passage percolation

abstract: We show that the number of maximal paths in directed last-passage percolation on the hypercubic lattice ${\mathbb Z}^d$ $(d\geq2)$ in which weights take finitely many values is typically exponentially large.


title: A radio frequency study of the accreting millisecond X-ray pulsar, IGR J16597$-$3704, in the globular cluster NGC 6256

abstract: We present Karl G. Jansky Very Large Array radio frequency observations of the new accreting millisecond X-ray pulsar (AMXP), IGR J16597$-$3704, located in the globular cluster NGC 6256. With these data, we detect a radio counterpart to IGR J16597$-$3704, and determine an improved source position. Pairing our radio observations with quasi-simultaneous Swift/XRT X-ray observations, we place IGR J16597$-$3704 on the radio -- X-ray luminosity plane, where we find that IGR J16597$-$3704 is one of the more radio-quiet neutron star low-mass X-ray binaries known to date. We discuss the mechanisms that may govern radio luminosity (and in turn jet production and evolution) in AMXPs. Further, we use our derived radio position to search for a counterpart in archival Hubble Space Telescope and Chandra X-ray Observatory data, and estimate an upper limit on the X-ray luminosity of IGR J16597$-$3704 during quiescence.


title: Rationality of Darmon points over genus fields of non-maximal orders

abstract: Stark-Heegner points, also known as Darmon points, were introduced by H. Darmon as certain local points on rational elliptic curves, conjecturally defined over abelian extensions of real quadratic fields. The rationality conjecture for these points is only known in the unramified case, namely, when these points are specializations of global points defined over the strict Hilbert class field $H^+_F$ of the real quadratic field $F$ and twisted by (unramified) quadratic characters of $Gal(H_c^+/F)$. We extend these results to the situation of ramified quadratic characters; more precisely, we show that Darmon points of conductor $c\geq 1$ twisted by quadratic characters of $G_c^+=Gal(H_c^+/F)$, where $H_c^+$ is the strict ring class field of $F$ of conductor $c$, come from rational points on the elliptic curve defined over $H_c^+$.


title: A New Type of Extreme-mass-ratio Inspirals Produced by Tidal Capture of Binary Black Holes

abstract: Extreme-mass-ratio inspiral (EMRI) is an important gravitational-wave (GW) source and it normally consists of one stellar-mass black hole (BH) whirling closely around a supermassive black hole (SMBH). In this Letter, we demonstrate that the small body, in fact, could be a BH binary (BHB). Previous numerical scatting experiments have shown that SMBHs can tidally capture BHBs to bound orbits. Here we investigate the subsequent long-term evolution. We find that those BHBs with a semi-major axis of $a\lesssim5\times10^{-3}$ AU can be captured to tightly-bound orbits such that they will successfully inspiral towards the central SMBHs without being scattered away by stellar relaxation processes. We estimate that these binary-EMRIs (b-EMRIs) could constitute at most $10\%$ of the EMRI population. Moreover, we show that when the eccentricity of a b-EMRI drops to about $0.85$, the two stellar BHs will quickly merge due to the tidal perturbation by the SMBH. The high-frequency ($\sim10^2$ Hz) GWs generated during the coalescence coincide with the low-frequency ($\sim10^{-3}$ Hz) waves from the b-EMRI, making this system an ideal target for future multi-band GW observations.


title: The Propitious Role of Solar Energetic Particles in the Origin of Life

abstract: We carry out 3-D numerical simulations to assess the penetration and bombardment effects of Solar Energetic Particles (SEPs), i.e. high-energy particle bursts during large flares and superflares, on ancient and current Mars. We demonstrate that the deposition of SEPs is non-uniform at the planetary surface, and that the corresponding energy flux is lower than other sources postulated to have influenced the origin of life. Nevertheless, SEPs may have been capable of facilitating the synthesis of a wide range of vital organic molecules (e.g. nucleobases and amino acids). Owing to the relatively high efficiency of these pathways, the overall yields might be comparable to (or even exceed) the values predicted for some conventional sources such as electrical discharges and exogenous delivery by meteorites. We also suggest that SEPs could have played a role in enabling the initiation of lightning. A notable corollary of our work is that SEPs may constitute an important mechanism for prebiotic synthesis on exoplanets around M-dwarfs, thereby mitigating the deficiency of biologically active ultraviolet radiation on these planets. Although there are several uncertainties associated with (exo)planetary environments and prebiotic chemical pathways, our study illustrates that SEPs represent a potentially important factor in understanding the origin of life.

categories: astro-ph.EP astro-ph.HE astro-ph.SR physics.bio-ph physics.space-ph


title: Superpositions of the cosmological constant allow for singularity resolution and unitary evolution in quantum cosmology

abstract: A novel approach to quantization is shown to allow for superpositions of the cosmological constant in isotropic and homogeneous mini-superspace models. Generic solutions featuring such superpositions display unitary evolution and resolution of the classical singularity. Physically well-motivated cosmological solutions are constructed. These particular solutions exhibit characteristic features of a cosmic bounce including universal phenomenology that can be rendered insensitive to Planck-scale physics in a natural manner.


title: The Uniform Electron Gas at Warm Dense Matter Conditions

abstract: We review the uniform electron gas (UEG) at finite temperature and over a broad density range relevant for warm dense matter (WDM) applications. We provide an overview of different simulation techniques, focusing on recent developments in the dielectric formalism and quantum Monte Carlo (QMC). Our primary focus is on two novel QMC methods: Permutation blocking path integral MC (PB-PIMC) and configuration PIMC (CPIMC). In fact, a combination of PB-PIMC and CPIMC has allowed for a highly accurate description of the warm dense UEG over a broad density-temperature range. We are able to effectively avoid the notorious fermion sign problem, without invoking uncontrolled approximations such as the fixed node approximation. Furthermore, a new finite-size correction scheme is presented that makes it possible to treat the UEG in the thermodynamic limit without loss of accuracy. In addition, we in detail discuss the construction of a parametrization of the exchange-correlation free energy that provides a complete description of the UEG and is of crucial importance as input for the simulation of real WDM applications. Further, we test previous theories, including restricted PIMC, finite-temperature Green functions, the classical mapping by Perrot and Dharma-wardana, and various dielectric methods such as the random phase approximation, or the Singwi-Tosi-Land-Sj\"olander, Vashishta-Singwi and the recent Tanaka scheme for the local field correction. Thus, for the first time, thorough benchmarks of important approximation schemes regarding various quantities such as different energies, in particular the exchange-correlation free energy, and the static structure factor, are possible. Finally, we outline a way how to rigorously extend our QMC studies to the inhomogeneous electron gas and present first ab initio data for the static density response and for the static local field correction.


title: Machine learning action parameters in lattice quantum chromodynamics

abstract: Numerical lattice quantum chromodynamics studies of the strong interaction are important in many aspects of particle and nuclear physics. Such studies require significant computing resources to undertake. A number of proposed methods promise improved efficiency of lattice calculations, and access to regions of parameter space that are currently computationally intractable, via multi-scale action-matching approaches that necessitate parametric regression of generated lattice datasets. The applicability of machine learning to this regression task is investigated, with deep neural networks found to provide an efficient solution even in cases where approaches such as principal component analysis fail. The high information content and complex symmetries inherent in lattice QCD datasets require custom neural network layers to be introduced and present opportunities for further development.


title: Flavor-Safe Light Squarks in Higgs-Anomaly Mediation

abstract: We consider a simple setup with light squarks which is free from the gravitino and SUSY flavor problems. In our setup, a SUSY breaking sector is sequestered from the matter and gauge sectors, and it only couples to the Higgs sector directly with $\mathcal{O}(100)\,$TeV gravitino. Resulting mass spectra of sfermions are split: the first and second generation sfermions are light as $\mathcal{O}(1)\,$TeV while the third generation sfermions are heavy as $\mathcal{O}(10)\,$TeV. The light squarks of $\mathcal{O}(1)\,$TeV can be searched at the (high-luminosity) LHC and future collider experiments. Our scenario can naturally avoid too large flavor-changing neutral currents and it is consistent with the $\epsilon_K$ constraint. Moreover, there are regions explaining the muon $g-2$ anomaly and bottom-tau/top-bottom-tau Yukawa coupling unification simultaneously.


title: Covariant formulation of scalar-torsion gravity

abstract: We consider a generalized teleparallel theory of gravitation, where the action contains an arbitrary function of the torsion scalar and a scalar field, $f(T,\phi)$, thus encompassing the cases of $f(T)$ gravity and nonminimally coupled scalar field as subclasses. The action is manifestly Lorentz invariant when besides the tetrad one allows for flat but nontrivial spin connection. We derive the field equations and demonstrate how the antisymmetric part of the tetrad equations is automatically satisfied when the spin connection equation holds. The spin connection equation is a vital part of the covariant formulation, since it determines the spin connection compatible with a given tetrad. We discuss how the spin connection equation can be solved in general, and provide the cosmological and spherically symmetric examples. Finally we generalize the theory to an arbitrary number of scalar fields.


title: Faster gaze prediction with dense networks and Fisher pruning

abstract: Predicting human fixations from images has recently seen large improvements by leveraging deep representations which were pretrained for object recognition. However, as we show in this paper, these networks are highly overparameterized for the task of fixation prediction. We first present a simple yet principled greedy pruning method which we call Fisher pruning. Through a combination of knowledge distillation and Fisher pruning, we obtain much more runtime-efficient architectures for saliency prediction, achieving a 10x speedup for the same AUC performance as a state of the art network on the CAT2000 dataset. Speeding up single-image gaze prediction is important for many real-world applications, but it is also a crucial step in the development of video saliency models, where the amount of data to be processed is substantially larger.


title: Soft-gluon and Coulomb corrections to hadronic top-quark pair production beyond NNLO

abstract: We construct a resummation at partial next-to-next-to-next-to-leading logarithmic accuracy for hadronic top-quark pair production near partonic threshold, including simultaneously soft-gluon and Coulomb corrections, and use this result to obtain approximate next-to-next-to-next-to-leading order predictions for the total top-quark pair-production cross section at the LHC. We generalize a required one-loop potential in non-relativistic QCD to the colour-octet case and estimate the remaining unknown two-loop potentials and three-loop anomalous dimensions. We obtain a moderate correction of $1.5\%$ relative to the next-to-next-to-leading order prediction and observe a reduction of the perturbative uncertainty below $\pm 5\%$.


title: Bouncing Unitary Cosmology I: Mini-Superspace General Solution

abstract: We offer a new proposal for cosmic singularity resolution based upon a quantum cosmology with a unitary bounce. This proposal is illustrated via a novel quantization of a mini-superspace model in which there can be superpositions of the cosmological constant. This possibility leads to a finite, bouncing unitary cosmology. Whereas the usual Wheeler--DeWitt cosmology generically displays pathological behaviour in terms of non-finite expectation values and non-unitary dynamics, the finiteness and unitarity of our model are formally guaranteed. For classically singular models with a massless scalar field and cosmological constant, we show that well-behaved quantum observables can be constructed and generic solutions to the universal Schr\"{o}dinger equation are singularity-free. Generic solutions of our model displays novel features including: i) superpositions of values of the cosmological constant; ii) universal effective physics due to non-trivial self-adjoint extensions of the Hamiltonian; and iii) bound `Efimov universe' states for negative cosmological constant. The last feature provides a new platform for quantum simulation of the early universe. A companion paper provides detailed interpretation and analysis of particular cosmological solutions that display a cosmic bounce due to quantum gravitational effects, a well-defined FLRW limit far from the bounce, and a semi-classical turnaround point in the dynamics of the scalar field which resembles an effective inflationary epoch.


title: Exploring Hidden Acoustic Spin Underwater

abstract: Investigating wave propagation in fluid enables a variety of important applications in underwater communications, object detections and unmanned robot control. Conventionally, momentum and spin reveal fundamental physical properties about propagating waves. Yet, vast previous research focused on the orbital angular momentum of acoustics without thinking about the existence possibility of spin due to the longitudinal wave nature. Here, we show that underwater acoustic wave processes the non-trivial spin angular momentum intrinsically, which is associated with its special spin-orbital coupling relation for longitudinal waves. Furthermore, we demonstrate that this intrinsic spin, although unobservable in plane wave form, can be detected by four approaches: wave interference, Gaussian exponential decay form, boundary evanescent wave, and waveguide mode. We further show that the strong spin-orbital coupling can be exploited to achieve unidirectional excitation and backscattering immune transport. We hope the present results can improve the geometric and topological understanding about underwater acoustic wave and pave the way on the spin-related underwater applications.


title: Pathwise Convergence of the Hard Spheres Kac Process

abstract: We derive two estimates for the deviation of the $N$-particle, hard-spheres Kac process from the corresponding Boltzmann equation, measured in expected Wasserstein distance. Particular care is paid to the long-time properties of our estimates, exploiting the stability properties of the limiting Boltzmann equation at the level of realisations of the interacting particle system. As a consequence, we obtain an estimate for the propagation of chaos, uniformly in time and with polynomial rates, as soon as the initial data has a $k^\mathrm{th}$ moment, $k>2$. Our approach is similar to Kac's proposal of relating the long-time behaviour of the particle system to that of the limit equation. Along the way, we prove a new estimate for the continuity of the Boltzmann flow measured in Wasserstein distance.


title: A Missing Proof of The Gibbard-Satterthwaite Theorem

abstract: A short and direct proof of the Gibbard-Satterthwaite theorem \`{a} la Amartya Sen's proof of Arrow's impossibility theorem is given.


title: First-principles quantum transport modeling of spin-transfer and spin-orbit torques in magnetic multilayers

abstract: We review a unified approach for computing: (i) spin-transfer torque in magnetic trilayers like spin-valves and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces; and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled-material, where injected charge current flows parallel to the interface. Our approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green's functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. We provide illustrative examples by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin-valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations; as well as spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide.


title: An unconditionally stable semi-implicit CutFEM for an interaction problem between an elastic membrane and an incompressible fluid

abstract: In this paper we introduce a finite element method for the Stokes equations with a massless immersed membrane. This membrane applies normal and tangential forces affecting the velocity and pressure of the fluid. Additionally, the points representing this membrane move with the local fluid velocity. We design and implement a high-accuracy cut finite element method (CutFEM) which enables the use of a structured mesh that is not aligned with the immersed membrane and then we formulate a time discretization that yields an unconditionally energy stable scheme. We prove that the stability is not restricted by the parameter choices that constrained previous finite element immersed boundary methods and illustrate the theoretical results with numerical simulations.


title: Shortest Path and Maximum Flow Problems Under Service Function Chaining Constraints

abstract: With the advent of Network Function Virtualization (NFV), Physical Network Functions (PNFs) are gradually being replaced by Virtual Network Functions (VNFs) that are hosted on general purpose servers. Depending on the call flows for specific services, the packets need to pass through an ordered set of network functions (physical or virtual) called Service Function Chains (SFC) before reaching the destination. Conceivably for the next few years during this transition, these networks would have a mix of PNFs and VNFs, which brings an interesting mix of network problems that are studied in this paper: (1) How to find an SFC-constrained shortest path between any pair of nodes? (2) What is the achievable SFC-constrained maximum flow? (3) How to place the VNFs such that the cost (the number of nodes to be virtualized) is minimized, while the maximum flow of the original network can still be achieved even under the SFC constraint? In this work, we will try to address such emerging questions. First, for the SFC-constrained shortest path problem, we propose a transformation of the network graph to minimize the computational complexity of subsequent applications of any shortest path algorithm. Second, we formulate the SFC-constrained maximum flow problem as a fractional multicommodity flow problem, and develop a combinatorial algorithm for a special case of practical interest. Third, we prove that the VNFs placement problem is NP-hard and present an alternative Integer Linear Programming (ILP) formulation. Finally, we conduct simulations to elucidate our theoretical results.


title: Towards a computational model of social norms

abstract: We describe a computational model of social norms based on identifying values that a certain culture finds desirable such as dignity, generosity and politeness. The model quantifies these values in the form of Culture-Sanctioned Social Metrics (CSSMs) and treats social norms as the requirement to maximize these metrics from the perspective of the self, peers and public. This model can be used to create realistic social simulations, to explain or predict human behavior in specific scenarios, or as a component of robots or agents that need to interact with humans in specific social-cultural settings. We validate the model by using it to represent a complex deception scenario and showing that it can yield non-trivial insights such as the explanation of apparently irrational human behavior.


title: Real Time Impact Control on Charging Energy Storage For Shipboard Power Systems

abstract: Medium voltage direct-current based integrated power system is projected as one of the solutions for powering the all-electric ship. It faces significant challenges for accurately energizing advanced loads, especially the pulsed power load, which can be rail gun, high power radar, and other state of art equipment. Energy storage based on supercapacitors is proposed as a technique for buffering the direct impact of pulsed power load on the power systems. However, the high magnitude of charging current of the energy storage can pose as a disturbance to both distribution and generation systems. This paper presents a fast switching device based real time control system that can achieve a desired balance between maintaining the required power quality and fast charging the energy storage in required time. Test results are shown to verify the performance of the proposed control algorithm.


title: An effect-theoretic reconstruction of quantum theory

abstract: An often used model for quantum theory is to associate to every physical system a C*-algebra. From a physical point of view it is unclear why operator algebras would form a good description of nature. In this paper, we find a set of physically meaningful assumptions such that any physical theory satisfying these assumptions must embed into the category of finite-dimensional C*-algebras. These assumptions were originally introduced in the setting of effectus theory, a categorical logical framework generalizing classical and quantum logic. As these assumptions have a physical interpretation, this motivates the usage of operator algebras as a model for quantum theory. In contrast to other reconstructions of quantum theory, we do not start with the framework of generalized probabilistic theories and instead use effect theories where no convex structure and no tensor product needs to be present. The lack of this structure in effectus theory has led to a different notion of pure maps. A map in an effectus is pure when it is a composition of a compression and a filter. These maps satisfy particular universal properties and respectively correspond to `forgetting' and `measuring' the validity of an effect. We define a pure effect theory (PET) to be an effect theory where the pure maps form a dagger-category and filters and compressions are adjoint. We show that any convex finite-dimensional PET must embed into the category of Euclidean Jordan algebras. Moreover, if the PET also has monoidal structure, then we show that it must embed into either the category of real or complex C*-algebras, which completes our reconstruction.


title: Symmetrization, factorization and arithmetic of quasi-Banach function spaces

abstract: We investigate relations between symmetrizations of quasi-Banach function spaces and constructions such as Calderon-Lozanovskii spaces, pointwise product spaces and pointwise multipliers. We show that under reasonable assumptions the symmetrization commutes with these operations. We determine also the spaces of pointwise multipliers between Lorentz spaces and Cesaro spaces. Developed methods may be regarded as an arithmetic of quasi-Banach function spaces and proofs of Theorems 3, 4 and 6 give a kind of tutorial for these methods. Finally, the above results will be used in proofs of some factorization results.


title: Interactive in-base street model edit: how common GIS software and a database can serve as a custom Graphical User Interface

abstract: Our modern world produces an increasing quantity of data, and especially geospatial data, with advance of sensing technologies, and growing complexity and organisation of vector data. Tools are needed to efficiently create and edit those vector geospatial data. Procedural generation has been a tool of choice to generate strongly organised data, yet it may be hard to control. Because those data may be involved to take consequence-full real life decisions, user interactions are required to check data and edit it. The classical process to do so would be to build an adhoc Graphical User Interface (GUI) tool adapted for the model and method being used. This task is difficult, takes a large amount of resources, and is very specific to one model, making it hard to share and re-use. Besides, many common generic GUI already exists to edit vector data, each having its specialities. We propose a change of paradigm; instead of building a specific tool for one task, we use common GIS software as GUIs, and deport the specific interactions from the software to within the database. In this paradigm, GIS software simply modify geometry and attributes of database layers, and those changes are used by the database to perform automated task. This new paradigm has many advantages. The first one is genericity. With in-base interaction, any GIS software can be used to perform edition, whatever the software is a Desktop sofware or a web application. The second is concurrency and coherency. Because interaction is in-base, use of database features allows seamless multi-user work, and can guarantee that the data is in a coherent state. Last we propose tools to facilitate multi-user edits, both during the edit phase (each user knows what areas are edited by other users), and before and after edit (planning of edit, analyse of edited areas).


title: Invariant random subgroups of groups acting on rooted trees

abstract: We investigate invariant random subgroups in groups acting on rooted trees. Let $\mathrm{Alt}_f(T)$ be the group of finitary even automorphisms of the $d$-ary rooted tree $T$. We prove that a nontrivial ergodic IRS of $\mathrm{Alt}_f(T)$ that acts without fixed points on the boundary of $T$ contains a level stabilizer, in particular it is the random conjugate of a finite index subgroup. Applying the technique to branch groups we prove that an ergodic IRS in a finitary regular branch group contains the derived subgroup of a generalized rigid level stabilizer. We also prove that every weakly branch group has continuum many distinct atomless ergodic IRS's. This extends a result of Benli, Grigorchuk and Nagnibeda who exhibit a group of intermediate growth with this property.


title: Measuring, Understanding, and Classifying News Media Sympathy on Twitter after Crisis Events

abstract: This paper investigates bias in coverage between Western and Arab media on Twitter after the November 2015 Beirut and Paris terror attacks. Using two Twitter datasets covering each attack, we investigate how Western and Arab media differed in coverage bias, sympathy bias, and resulting information propagation. We crowdsourced sympathy and sentiment labels for 2,390 tweets across four languages (English, Arabic, French, German), built a regression model to characterize sympathy, and thereafter trained a deep convolutional neural network to predict sympathy. Key findings show: (a) both events were disproportionately covered (b) Western media exhibited less sympathy, where each media coverage was more sympathetic towards the country affected in their respective region (c) Sympathy predictions supported ground truth analysis that Western media was less sympathetic than Arab media (d) Sympathetic tweets do not spread any further. We discuss our results in light of global news flow, Twitter affordances, and public perception impact.


title: Model-free inference of direct network interactions from nonlinear collective dynamics

abstract: The topology of interactions in network dynamical systems fundamentally underlies their function. Accelerating technological progress creates massively available data about collective nonlinear dynamics in physical, biological, and technological systems. Detecting direct interaction patterns from those dynamics still constitutes a major open problem. In particular, current nonlinear dynamics approaches mostly require to know a priori a model of the (often high dimensional) system dynamics. Here we develop a model-independent framework for inferring direct interactions solely from recording the nonlinear collective dynamics generated. Introducing an explicit dependency matrix in combination with a block-orthogonal regression algorithm, the approach works reliably across many dynamical regimes, including transient dynamics toward steady states, periodic and non-periodic dynamics, and chaos. Together with its capabilities to reveal network (two point) as well as hypernetwork (e.g., three point) interactions, this framework may thus open up nonlinear dynamics options of inferring direct interaction patterns across systems where no model is known.

categories: physics.soc-ph nlin.CD physics.comp-ph physics.data-an q-bio.QM


title: Resolving discrete pulsar spin-down states with current and future instrumentation

abstract: An understanding of pulsar timing noise offers the potential to improve the timing precision of a large number of pulsars as well as facilitating our understanding of pulsar magnetospheres. For some sources, timing noise is attributable to a pulsar switching between two different spin-down rates $(\dot{\nu})$. Such transitions may be common but difficult to resolve using current techniques. In this work, we use simulations of $\dot{\nu}$-variable pulsars to investigate the likelihood of resolving individual $\dot{\nu}$ transitions. We inject step-changes in the value of $\dot{\nu}$ with a wide range of amplitudes and switching timescales. We then attempt to redetect these transitions using standard pulsar timing techniques. The pulse arrival-time precision and the observing cadence are varied. Limits on $\dot{\nu}$ detectability based on the effects such transitions have on the timing residuals are derived. With the typical cadences and timing precision of current timing programs, we find we are insensitive to a large region of $\Delta \dot{\nu}$ parameter space which encompasses small, short timescale switches. We find, where the rotation and emission states are correlated, that using changes to the pulse shape to estimate $\dot{\nu}$ transition epochs, can improve detectability in certain scenarios. The effects of cadence on $\Delta \dot{\nu}$ detectability are discussed and we make comparisons with a known population of intermittent and mode-switching pulsars. We conclude that for short timescale, small switches, cadence should not be compromised when new generations of ultra-sensitive radio telescopes are online.


title: Cosmology and Accelerator Tests of Strongly Interacting Dark Matter

abstract: A natural possibility for dark matter is that it is composed of the stable pions of a QCD-like hidden sector. Existing literature largely assumes that pion self-interactions alone control the early universe cosmology. We point out that processes involving vector mesons typically dominate the physics of dark matter freeze-out and significantly widen the viable mass range for these models. The vector mesons also give rise to striking signals at accelerators. For example, in most of the cosmologically favored parameter space, the vector mesons are naturally long-lived and produce Standard Model particles in their decays. Electron and proton beam fixed-target experiments such as HPS, SeaQuest, and LDMX can exploit these signals to explore much of the viable parameter space. We also comment on dark matter decay inherent in a large class of previously considered models and explain how to ensure dark matter stability.