title: Critical phases in the raise and peel model

abstract: The raise and peel model (RPM) is a nonlocal stochastic model describing the space and time fluctuations of an evolving one dimensional interface. Its relevant parameter $u$ is the ratio between the rates of local adsorption and nonlocal desorption processes (avalanches) processes. The model at $u=1$ give us the first example of a conformally invariant stochastic model. For small values $u<u_0$ the model is known to be noncritical, while for $u>u_0$ it is critical. By calculating the structure function of the height profiles in the reciprocal space we confirm with good precision that indeed $u_0=1$. We establish that at the conformal invariant point $u=1$ the RPM has a roughness transition with dynamical and roughness critical exponents $z=1$ and $\alpha=0$, respectively. For $u>1$ the model is critical with an $u$-dependent dynamical critical exponent $z(u)$ that tends towards zero as $u\to \infty$. However at $1/u=0$ the RPM is exactly mapped into the totally asymmetric exclusion problem (TASEP). This last model is known to be noncritical (critical) for open (periodic) boundary conditions. Our studies indicate that the RPM as $u \to \infty$, due to its nonlocal dynamics processes, has the same large-distance physics no matter what boundary condition we chose. For $u>1$, our analysis show that differently from previous predictions, the region is composed by two distinct critical phases. For $u\leq u < u_c\approx 40$ the height profiles are rough ($\alpha = \alpha(u) >0$), and for $u>u_c$ the height profiles are flat at large distances ($\alpha = \alpha(u) <0$). We also observed that in both critical phases ($u>1$) the RPM at short length scales, has an effective behavior in the Kardar-Parisi-Zhang (KPZ) critical universality class, that is not the true behavior of the system at large length scales.


title: On the role of multiple scales in metapopulations of public good producers

abstract: Multiple scales in metapopulations can give rise to paradoxical behaviour: in a conceptual model for a public goods game, the species associated with a fitness cost due to the public good production can be stabilised in the well-mixed limit due to the mere existence of these scales. The scales in this model involve a length scale corresponding to separate patches, coupled by mobility, and separate time scales for reproduction and interaction with a local environment. Contrary to the well-mixed high mobility limit, we find that for low mobilities, the interaction rate progressively stabilises this species due to stochastic effects, and that the formation of spatial patterns is not crucial for this stabilisation.


title: Exponential equilibration of genetic circuits using entropy methods

abstract: We analyse a continuum model for genetic circuits based on a partial integro-differential equation initially proposed in Friedman, Cai \& Xie (2006) as an approximation of a chemical master equation. We use entropy methods to show exponentially fast convergence to equilibrium for this model with explicit bounds. The asymptotic equilibration for the multidimensional case of more than one gene is also obtained under suitable assumptions on the equilibrium stationary states. The asymptotic equilibration property for networks involving one and more than one gene is investigated via numerical simulations.


title: Monogamy of particle statistics in tripartite systems simulating bosons and fermions

abstract: In the quantum world correlations can take form of entanglement which is known to be monogamous. In this work we show that another type of correlations, indistinguishability, is also restricted by some form of monogamy. Namely, if particles A and B simulate bosons, then A and C cannot perfectly imitate fermions. Our main result consists in demonstrating to what extent it is possible.


title: Semi Concurrent vector fields in Finsler geometry

abstract: In the present paper, we introduce and investigate the notion of a semi concurrent vector field on a Finsler manifold. We show that some special Finsler manifolds admitting such vector fields turn out to be Riemannian. We prove that Tachibana's characterization of Finsler manifolds admitting a concurrent vector field leads to Riemannain metrics. We give an answer to the question raised in \cite{DWF}: "Is any n-dimensional Finsler manifold $(M,F)$, admitting a non-constant smooth function $f$ on $M$ such that $\frac{\partial f}{\partial x^i}\frac{\partial g^{ij}}{\partial y^k}=0$, a Riemannian manifold?". Various examples for conic Finsler and Riemannian spaces that admit semi-concurrent vector field are presented. Finally, we conjectured that there is no regular Finsler non-Riemannian metric that admits a semi-concurrent vector field. In other words, a Finsler metric admitting a semi-concurrent vector field is necessarily either Riemannian or conic Finslerian.


title: Nonrelativistic Limit of Dirac Theory From Effective Field Theory

abstract: In this work we analyze the low energy nonrelativistic limit of Dirac theory in the framework of effective field theory. By integrating out the high energy modes of Dirac field, given in terms of a combination of the two-components Weyl spinors, we obtain a low energy effective action for the remaining components, whose equation of motion can then be compared to the Pauli-Schr\"odinger equation after demanding normalization of the wave function. We then discuss the relevance of the terms in the effective action in the context of an anisotropic dimensional analysis which is suitable for nonrelativistic theories.


title: Management of the Correlations of Ultracold Bosons in Triple Wells

abstract: Ultracold interacting atoms are an excellent tool to study correlation functions of many-body systems that are generally eluding detection and manipulation. Herein, we investigate the ground state of bosons in a tilted triple-well potential and characterize the many-body state by the eigenvalues of its reduced one-body density matrix and Glauber correlation functions. We unveil how the interplay between the interaction strength and the tilt can be used to control the number of correlated wells as well as the fragmentation, i.e. the number of macroscopic eigenvalues of the reduced one-body density matrix.


title: Exploiting ecology in drug pulse sequences in favour of population reduction

abstract: A deterministic population dynamics model involving birth and death for a two-species system, comprising a wild-type and more resistant species competing via logistic growth, is subjected to two distinct stress environments designed to mimic those that would typically be induced by temporal variation in the concentration of a drug (antibiotic or chemotherapeutic) as it permeates through the population and is progressively degraded. Different treatment regimes, involving single or periodical doses, are evaluated in terms of the minimal population size (a measure of the extinction probability), and the population composition (a measure of the selection pressure for resistance or tolerance during the treatment). We show that there exist timescales over which the low-stress regime is as effective as the high-stress regime, due to the competition between the two species. For multiple periodic treatments, competition can ensure that the minimal population size is attained during the first pulse when the high-stress regime is short, which implies that a single short pulse can be more effective than a more protracted regime. Our results suggest that when the duration of the high-stress environment is restricted, a treatment with one or multiple shorter pulses can produce better outcomes than a single long treatment. If ecological competition is to be exploited for treatments, it is crucial to determine these timescales, and estimate for the minimal population threshold that suffices for extinction. These parameters can be quantified by experiment.


title: Unique Quasi-Stationary Distribution, with a possibly stabilizing extinction

abstract: We establish sufficient conditions for exponential convergence to a unique quasi-stationary distribution in the total variation norm. These conditions also ensure the existence and exponential ergodicity of the Q-process, the process conditionned upon never being absorbed. The technique relies on a coupling procedure that is related to Harris recurrence (for Markov Chains). It applies to general continuous-time and continuous-space Markov processes. The main novelty is that we modulate each coupling step depending both on a final horizon of time (for survival) and on the initial distribution. By this way, we could notably include in the convergence a dependency on the initial condition. As an illustration, we consider a continuous-time birth-death process with catastrophes and a diffusion process describing a (localized) population adapting to its environment.


title: Gundy-Varopoulos martingale transforms and their projection operators on manifolds and vector bundles

abstract: This paper proves the $L^p$ boundedness of generalized first order Riesz transforms obtained as conditional expectations of martingale transforms \`a la Gundy-Varopoulos for quite general diffusions on manifolds and vector bundles. Several specific examples and applications are presented: Lie groups of compact type, the Heisenberg group, SU(2), and Riesz transforms on forms and spinors.


title: Gauge-invariant formulation of time-dependent configuration interaction singles method

abstract: We propose a gauge-invariant formulation of the channel orbital-based time-dependent configuration interaction singles (TDCIS) method [Phys. Rev. A 74, 043420 (2006)], one of the powerful ab initio methods to investigate electron dynamics in atoms and molecules subject to an external laser field. In the present formulation, we derive the equations of motion (EOMs) in the velocity gauge using gauge-transformed orbitals, not fixed orbitals, that are equivalent to the conventional EOMs in the length gauge using fixed orbitals. The new velocity-gauge EOMs avoid the use of the length-gauge dipole operator, which diverges at large distance, and allows to exploit computational advantages of the velocity-gauge treatment over the length-gauge one, e.g, a faster convergence in simulations with intense and long-wavelength lasers, and the feasibility of exterior complex scaling as an absorbing boundary. The reformulated TDCIS method is applied to an exactly solvable model of one-dimensional helium atom in an intense laser field to numerically demonstrate the gauge invariance. We also discuss the consistent method for evaluating the time derivative of an observable, relevant e.g, in simulating high-harmonic generation.


title: Analysis of stochastic bifurcations with phase portraits

abstract: We propose a method to obtain phase portraits for stochastic systems. Starting from the Fokker-Planck equation, we separate the dynamics into a convective and a diffusive part. We show that stable and unstable fixed points of the convective field correspond to maxima and minima of the stationary probability distribution if the probability current vanishes at these points. Stochastic phase portraits, which are vector plots of the convective field, therefore indicate the extrema of the stationary distribution and can be used to identify stochastic bifurcations that change the number and stability of these extrema. We show that limit cycles in stochastic phase portraits can indicate ridges of the probability distribution, and we identify a novel type of stochastic bifurcations, where the probability maximum moves to the edge of the system through a gap between the two nullclines of the convective field.


title: Structure, magnetic and transport properties of epitaxial thin films of equiatomic CoFeMnGe quaternary Heusler alloy

abstract: Future spintronics requires the realization of thin film of half-metallic ferromagnets having high Curie temperature and 100\% spin polarization at the Fermi level for potential spintronics applications. In this paper, we report the epitaxial thin films growth of half-metallic CoFeMnGe Heusler alloy on MgO (001) substrate using pulsed laser deposition system, along with the study of structural, magnetic and transport properties. The magnetic property measurements of the thin film suggest a soft ferromagnetic state at room temperature with an in-plane magnetic anisotropy and a Curie temperature well above the room temperature. Anisotropic magnetoresistance (AMR) ratio and temperature dependent electrical resistivity measurements of the thin film indicate the compound to be half-metallic in nature and therefore suitable for the fabrications of spintronics devices.


title: Cosmology and Dark Matter

abstract: This lecture course covers cosmology from the particle physicist perspective. Therefore, the emphasis will be on the evidence for the new physics in cosmological and astrophysical data together with minimal theoretical frameworks needed to understand and appreciate the evidence. I review the case for non-baryonic dark matter and describe popular models which incorporate it. In parallel, the story of dark energy will be developed, which includes accelerated expansion of the Universe today, the Universe origin in the Big Bang, and support for the Inflationary theory in CMBR data.


title: Breaking the current density threshold in spin-orbit-torque magnetic random access memory

abstract: Spin-orbit-torque magnetic random access memory (SOT-MRAM) is a promising technology for the next generation of data storage devices. The main bottleneck of this technology is the high reversal current density threshold. This outstanding problem of SOT-MRAM is now solved by using a current density of constant magnitude and varying flow direction that reduces the reversal current density threshold by a factor of more than the Gilbert damping coefficient. The Euler-Lagrange equation for the fastest magnetization reversal path and the optimal current pulse are derived for an arbitrary magnetic cell. The theoretical limit of minimal reversal current density and current density for a GHz switching rate of the new reversal strategy for CoFeB/Ta SOT-MRAMs are respectively of the order of $10^5$ A/cm$^2$ and $10^6$ A/cm$^2$ far below $10^7$ A/cm$^2$ and $10^8$ A/cm$^2$ in the conventional strategy. Furthermore, no external magnetic field is needed for a deterministic reversal in the new strategy.


title: Constraining Gluonic Quartic Gauge Coupling Operators with $gg \to \gamma \gamma$

abstract: Gluon-gluon to photon-photon scattering $gg \to \gamma\gamma$ offers to the LHC experiments a uniquely powerful probe of dimension-8 operators in the Standard Model Effective Field Theory (SMEFT) that are quadratic in both the electromagnetic and gluonic field-strength tensors, such as would appear in the Born-Infeld extension of the Standard Model (SM). We use 13-TeV ATLAS data on the production of isolated photon pairs to set lower limits on the scales of dimension-8 operators $M \gtrsim 1$~TeV, and discuss the prospective sensitivities of possible future hadron colliders.


title: Parton distributions and $\cos 2\phi_h$ asymmetry induced by anomalous photon-quark coupling

abstract: In the spectator models of the nucleon with scalar and axial-vector diquarks, we show that the effect of Pauli coupling in photon-quark vertex to the parton distribution functions (PDFs) of nucleon and azimuthal asymmetry in the unpolarized semi-inclusive deep inelastic scattering (SIDIS). This anomalous coupling gives obvious contribution to unpolarized and polarized PDFs, and also leads to a $\cos 2\phi_h$ azimuthal asymmetry proportional to the squared Pauli form factor, due to the helicity flip of the struck quark. After determining the model parameters by fitting PDFs to the global fits, this new distribution for $\cos 2\phi_h$ asymmetry is given numerically. In the framework of transverse momentum dependence (TMD), we find that it is positive and of a few percent in the kinematical regime of HERMES and COMPASS collaborations, in the same order of magnitude with Cahn effect.


title: Target selection for a small low-thrust mission to near-Earth asteroids

abstract: The preliminary mission design of spacecraft missions to asteroids often involves, in the early phases, the selection of candidate target asteroids. The final result of such an analysis is a list of asteroids, ranked with respect to the necessary propellant to be used, that the spacecraft could potentially reach. In this paper we investigate the sensitivity of the produced asteroids rank to the employed trajectory model in the specific case of a small low-thrust propelled spacecraft beginning its journey from the Sun-Earth $L_2$ Lagrangian point and heading to a rendezvous with some near-Earth asteroid. We consider five increasingly complex trajectory models: impulsive, Lambert, nuclear electric propulsion, nuclear electric propulsion including the Earth's gravity, solar electric propulsion including the Earth's gravity and we study the final correlation between the obtained target rankings. We find that the use of a low-thrust trajectory model is of great importance for target selection, since the use of chemical propulsion surrogates leads to favouring less attractive options 19\% of times, a percentage that drops to 8\% already using a simple nuclear electric propulsion model that neglects the Earth's gravity effects and thrust dependence on the solar distance. We also find that for the study case considered, a small interplanetary CubeSat named M-ARGO, the inclusion of the Earth's gravity in the considered dynamics does not affect the target selection significantly.


title: Algorithmic transformation of multi-loop Feynman integrals to a canonical basis

abstract: The evaluation of multi-loop Feynman integrals is one of the main challenges in the computation of precise theoretical predictions for the cross sections measured at the LHC. In recent years, the method of differential equations has proven to be a powerful tool for the computation of Feynman integrals. It has been observed that the differential equation of Feynman integrals can in many instances be transformed into a so-called canonical form, which significantly simplifies its integration in terms of iterated integrals. The main result of this thesis is an algorithm to compute rational transformations of differential equations of Feynman integrals into a canonical form. Apart from requiring the existence of such a rational transformation, the algorithm needs no further assumptions about the differential equation. In particular, it is applicable to problems depending on multiple kinematic variables and also allows for a rational dependence on the dimensional regulator. First, the transformation law is expanded in the dimensional regulator to derive differential equations for the coefficients of the transformation. Using an ansatz in terms of rational functions, these differential equations are then solved to determine the transformation. This thesis also presents an implementation of the algorithm in the Mathematica package CANONICA, which is the first publicly available program to compute transformations to a canonical form for differential equations depending on multiple variables. The main functionality and its usage are illustrated with some simple examples. Furthermore, the package is applied to state-of-the-art integral topologies appearing in recent multi-loop calculations. These topologies depend on up to three variables and include previously unknown topologies contributing to higher-order corrections to the cross section of single top-quark production at the LHC.


title: A group-theoretical interpretation of the word problem for free idempotent generated semigroups

abstract: The set of idempotents of any semigroup carries the structure of a biordered set, which contains a great deal of information concerning the idempotent generated subsemigroup of the semigroup in question. This leads to the construction of a free idempotent generated semigroup $\mathsf{IG}(\mathcal{E})$ - the `free-est' semigroup with a given biordered set $\mathcal{E}$ of idempotents. We show that when $\mathcal{E}$ is finite, the word problem for $\mathsf{IG}(\mathcal{E})$ is equivalent to a family of constraint satisfaction problems involving rational subsets of direct products of pairs of maximal subgroups of $\mathsf{IG}(\mathcal{E})$. As an application, we obtain decidability of the word problem for an important class of examples. Also, we prove that for finite $\mathcal{E}$, $\mathsf{IG}(\mathcal{E})$ is always a weakly abundant semigroup satisfying the congruence condition.


title: Size and host-medium effects on topologically protected surface states in bi-anisotropic 3D optical waveguides

abstract: We study the optical properties of bi-anisotropic optical waveguides with nontrivial topological structure in wavevector space, placed in an ordinary dielectric matrix. We derive an exact analytical description of the eigenmodes of the systems with arbitrary parameters that allows us to investigate topologically protected surface states (TPSS) in details. In particular, we find that the TPSS on the waveguides would disappear (1) if their radius is smaller than a critical radius due to the dimensional quantization of azimuthal wavenumber, and also (2) if the permittivity of the host-medium exceeds a critical value. Interestingly, we also find that the TPSS in the waveguides have negative refraction for some geometries. We have found a TPSS phase diagram that will pave the way for development of the topological waveguides for optical interconnects and devices.


title: Revisiting the Inverted Indices for Billion-Scale Approximate Nearest Neighbors

abstract: This work addresses the problem of billion-scale nearest neighbor search. The state-of-the-art retrieval systems for billion-scale databases are currently based on the inverted multi-index, the recently proposed generalization of the inverted index structure. The multi-index provides a very fine-grained partition of the feature space that allows extracting concise and accurate short-lists of candidates for the search queries. In this paper, we argue that the potential of the simple inverted index was not fully exploited in previous works and advocate its usage both for the highly-entangled deep descriptors and relatively disentangled SIFT descriptors. We introduce a new retrieval system that is based on the inverted index and outperforms the multi-index by a large margin for the same memory consumption and construction complexity. For example, our system achieves the state-of-the-art recall rates several times faster on the dataset of one billion deep descriptors compared to the efficient implementation of the inverted multi-index from the FAISS library.


title: On the Generalizability of Linear and Non-Linear Region of Interest-Based Multivariate Regression Models for fMRI Data

abstract: In contrast to conventional, univariate analysis, various types of multivariate analysis have been applied to functional magnetic resonance imaging (fMRI) data. In this paper, we compare two contemporary approaches for multivariate regression on task-based fMRI data: linear regression with ridge regularization and non-linear symbolic regression using genetic programming. The data for this project is representative of a contemporary fMRI experimental design for visual stimuli. Linear and non-linear models were generated for 10 subjects, with another 4 withheld for validation. Model quality is evaluated by comparing $R$ scores (Pearson product-moment correlation) in various contexts, including single run self-fit, within-subject generalization, and between-subject generalization. Propensity for modelling strategies to overfit is estimated using a separate resting state scan. Results suggest that neither method is objectively or inherently better than the other.


title: Interview validation of the Physics Lab Inventory of Critical thinking (PLIC)

abstract: Although an important goal of introductory physics labs is to train students in scientific reasoning and critical thinking, currently there are no standard tests in physics designed to assess such skills. We are in the process of developing and validating the Physics Lab Inventory of Critical thinking (PLIC), an assessment to probe students' critical thinking abilities in physics lab courses. The instrument asks students to critique a set of experimental methods and data and use them to evaluate a particular physical model (the period of a mass on a spring). Currently, we are validating the closed-response survey through interviews with students and present the results of 12 such interviews here. We describe a trend that has emerged from these interviews, with students' reasoning falling into three main patterns of behavior: selecting all options, cuing on keywords, and critically analyzing. We have found ways to shift students to the last and more desirable behavior. We discuss ways in which these findings are likely relevant to the design of other concept inventories.


title: Impact of Land Use on the DOM Composition in Different Seasons in a Subtropical River Flowing through the Region Undergoing Rapid Urbanization

abstract: The dissolved organic matter (DOM) composition in river ecosystems could reflect the human impacts on the river ecosystem, and plays an important role in the carbon cycling process. We collected water and phytoplankton samples at 107 sites in the Dongjiang River in two seasons to assess the impact of the sub-catchments land use structure on the DOM composition. The results showed that (1) the forested sub-catchments had higher humic-like C1 (16.45%) and C2 (25.04%) and lower protein-like C3 (22.57%) and C4 (35.95%) than urbanized and mixed forest-agriculture sub-catchments, while the urbanized sub-catchments showed an inverse trend (4.54%, 15.51%, 33.97% and 45.98%, respectively). (2) The significant variation in the proportion of C1 and C4 between the dry and rainy seasons was recorded in both the forested and the mixed forest-agriculture sub-catchments (p<0.01), but only C4 showed an obvious seasonal variation in the urbanized sub-catchments (p<0.01). While the DOM composition was mainly related to the proportion of urbanized land and forested land year-round (p<0.01), it had stronger correlation with forested land in the dry season and urbanized land in the rainy season. (3) No significant correlation between the DOM composition and the agricultural land proportion was found in either season (p>0.05). Our findings indicated that the DOM composition was strongly dependent on the land use structure of the sub-catchments and varied seasonally, but the seasonal variation pattern could be disturbed by human activities in the extensively urbanized catchments. Notably, the higher C4 proportion compared with those in temperate rivers implied a larger amount of CO2 released from the subtropical rivers into the atmosphere when considering bioavailability.


title: The linearization problem of a binary quadratic problem and its applications

abstract: We provide several applications of the linearization problem of a binary quadratic problem. We propose a new lower bounding strategy, called the linearization-based scheme, that is based on a simple certificate for a quadratic function to be non-negative on the feasible set. Each linearization-based bound requires a set of linearizable matrices as an input. We prove that the Generalized Gilmore-Lawler bounding scheme for binary quadratic problems provides linearization-based bounds. Moreover, we show that the bound obtained from the first level reformulation linearization technique is also a type of linearization-based bound, which enables us to provide a comparison among mentioned bounds. However, the strongest linearization-based bound is the one that uses the full characterization of the set of linearizable matrices. Finally, we present a polynomial-time algorithm for the linearization problem of the quadratic shortest path problem on directed acyclic graphs. Our algorithm gives a complete characterization of the set of linearizable matrices for the quadratic shortest path problem.


title: MRI Tumor Segmentation with Densely Connected 3D CNN

abstract: Glioma is one of the most common and aggressive types of primary brain tumors. The accurate segmentation of subcortical brain structures is crucial to the study of gliomas in that it helps the monitoring of the progression of gliomas and aids the evaluation of treatment outcomes. However, the large amount of required human labor makes it difficult to obtain the manually segmented Magnetic Resonance Imaging (MRI) data, limiting the use of precise quantitative measurements in the clinical practice. In this work, we try to address this problem by developing a 3D Convolutional Neural Network~(3D CNN) based model to automatically segment gliomas. The major difficulty of our segmentation model comes with the fact that the location, structure, and shape of gliomas vary significantly among different patients. In order to accurately classify each voxel, our model captures multi-scale contextual information by extracting features from two scales of receptive fields. To fully exploit the tumor structure, we propose a novel architecture that hierarchically segments different lesion regions of the necrotic and non-enhancing tumor~(NCR/NET), peritumoral edema~(ED) and GD-enhancing tumor~(ET). Additionally, we utilize densely connected convolutional blocks to further boost the performance. We train our model with a patch-wise training schema to mitigate the class imbalance problem. The proposed method is validated on the BraTS 2017 dataset and it achieves Dice scores of 0.72, 0.83 and 0.81 for the complete tumor, tumor core and enhancing tumor, respectively. These results are comparable to the reported state-of-the-art results, and our method is better than existing 3D-based methods in terms of compactness, time and space efficiency.


title: An elementary conjecture which implies the Goldbach conjecture

abstract: Let $p_{1}$, ..., $p_{k}$ be the first $k$ odd primes in succession. Let $n$ be an even integer such that $n > p_{k}$. We conjecture that if none of $n - p_{1}$, ..., $n - p_{k}$ are prime, then at least one of them has a prime factor which is greater than or equal to $p_{k}$. In this brief note, we observe that Goldbach's conjecture follows from this conjecture.


title: Maintenance of diversity in a hierarchical host-parasite model with balancing selection and reinfection

abstract: Inspired by DNA data of the human cytomegalovirus we propose a model of a two-type parasite population distributed over its hosts. The parasite is capable to persist in its host till the host dies, and to reinfect other hosts. To maintain type diversity within a host, balancing selection is assumed. For a suitable parameter regime we show that in the limit of large host and parasite populations the host state frequencies follow a dynamical system with a globally stable equilibrium, guaranteeing that both types are maintained in the parasite population for a long time on the host time scale.


title: From Elasticity to Electromagnetism: Beyond the Mirror

abstract: The first purpose of this short but striking paper is to revisit Elasticity (EL) and Electromagnetism (EM) by comparing the structure of these two theories and examining with details their well known couplings, in particular piezoelectricity and photoelasticity. Despite the strange Helmholtz and Mach-Lippmann analogies existing between them, no classical technique may provide a common setting. However, unexpected arguments discovered independently by the brothers E. and F. Cosserat in 1909 for EL and by H. Weyl in 1918 for EM are leading to construct a new differential sequence called Spencer sequence in the framework of the formal theory of Lie pseudogroups and to introduce it for the conformal group of space-time with 15 parameters. Then, all the previous explicit couplings can be deduced abstractly and one must just go to a laboratory in order to know about the coupling constants on which they are depending, like in the Hooke or Minkowski constitutive relations existing respectively in EL or EM separately. We finally provide a new combined experimental and theoretical proof of the fact that any 1-form with value in the second order jets (elations) of the conformal group of space-time can be uniquely decomposed into the direct sum of the Ricci tensor R and the electromagnetic field F. This result questions the mathematical foundations of both General Relativity (GR) and Gauge Theory (GT). In particular, the Einstein operator (6 terms) must be thus replaced by the formal adjoint of the Ricci operator (4 terms only) in the study of gravitational waves.


title: Canonicality of Makanin-Razborov Diagrams - Counterexample

abstract: Sets of solutions to finite systems of equations in a free group, are equivalent to sets of homomorphisms from a fixed f.p. group into a free group. The latter can be encoded in a diagram, the construction of which is valid also for f.g. groups. The diagram is known to be canonical for a fixed f.g. group with a fixed generating set. In this paper we prove that the construction depends on the chosen generating set of the given f.g. group.


title: CryptoRec: Privacy-preserving Recommendation as a Service

abstract: Recommender systems rely on large datasets of historical data and entail serious privacy risks. A server offering Recommendation as a Service to a client might leak more information than necessary regarding its recommendation model and dataset. At the same time, the disclosure of the client's preferences to the server is also a matter of concern. Devising privacy-preserving protocols using general cryptographic primitives (e.g., secure multi-party computation or homomorphic encryption), is a typical approach to overcome privacy concerns, but in conjunction with state-of-the-art recommender systems often yields far-from-practical solutions. In this paper, we tackle this problem from the direction of constructing crypto-friendly machine learning algorithms. In particular, we propose CryptoRec, a secure two-party computation protocol for Recommendation as a Service, which encompasses a novel recommender system. This model possesses two interesting properties: (1) It models user-item interactions in an item-only latent feature space in which personalized user representations are automatically captured by an aggregation of pre-learned item features. This means that a server with a pre-trained model can provide recommendations for a client whose data is not in its training set. Nevertheless, re-training the model with the client's data still improves accuracy. (2) It only uses addition and multiplication operations, making the model straightforwardly compatible with homomorphic encryption schemes. We demonstrate the efficiency and accuracy of CryptoRec on three real-world datasets. CryptoRec allows a server with thousands of items to privately answer a prediction query within a few seconds on a single PC, while its prediction accuracy is still competitive with state-of-the-art recommender systems computing over clear data.


title: The Binary $\mathfrak{aff}(n|1)$-Invariant Differential Operators On Weighted Densities On The Superspace $\mathbb{R}^{1|n}$ And $\mathfrak{aff}(n|1)$-Relative Cohomology

abstract: We consider the $\mathfrak{aff}(n|1)-$module structure on the spaces of differential bilinear operators acting on the superspaces of weighted densities. We classify $\mathfrak{aff}(n|1)-$invariant binary differential operators acting on the spaces of weighted densities. This result allows us to compute the first $\mathfrak{aff}(n|1)-$relative differential cohomology of $\mathcal{K}(n)$ with coefficients in the superspace of linear differential operators acting on the superspaces of weighted densities.


title: Evolutionary Computation plus Dynamic Programming for the Bi-Objective Travelling Thief Problem

abstract: This research proposes a novel indicator-based hybrid evolutionary approach that combines approximate and exact algorithms. We apply it to a new bi-criteria formulation of the travelling thief problem, which is known to the Evolutionary Computation community as a benchmark multi-component optimisation problem that interconnects two classical NP-hard problems: the travelling salesman problem and the 0-1 knapsack problem. Our approach employs the exact dynamic programming algorithm for the underlying Packing-While-Travelling (PWT) problem as a subroutine within a bi-objective evolutionary algorithm. This design takes advantage of the data extracted from Pareto fronts generated by the dynamic program to achieve better solutions. Furthermore, we develop a number of novel indicators and selection mechanisms to strengthen synergy of the two algorithmic components of our approach. The results of computational experiments show that the approach is capable to outperform the state-of-the-art results for the single-objective case of the problem.


title: Mixed-state localization operators: Cohen's class and trace class operators

abstract: We study mixed-state localization operators from the perspective of Werner's operator convolutions which allows us to extend known results from the rank-one case to trace class operators. The idea of localizing a signal to a domain in phase space is approached from various directions such as bounds on the spreading function, probability densities associated to mixed-state localization operators, positive operator valued measures, positive correspondence rules and variants of Tauberian theorems for operator translates. Our results include a rigorous treatment of multiwindow-STFT filters and a characterization of mixed-state localization operators as positive correspondence rules. Furthermore, we provide a description of the Cohen class in terms of Werner's convolution of operators and deduce consequences on positive Cohen class distributions, an uncertainty principle, uniqueness and phase retrieval for general elements of Cohen's class.


title: Stochastic Deconvolutional Neural Network Ensemble Training on Generative Pseudo-Adversarial Networks

abstract: The training of Generative Adversarial Networks is a difficult task mainly due to the nature of the networks. One such issue is when the generator and discriminator start oscillating, rather than converging to a fixed point. Another case can be when one agent becomes more adept than the other which results in the decrease of the other agent's ability to learn, reducing the learning capacity of the system as a whole. Additionally, there exists the problem of Mode Collapse which involves the generators output collapsing to a single sample or a small set of similar samples. To train GANs a careful selection of the architecture that is used along with a variety of other methods to improve training. Even when applying these methods there is low stability of training in relation to the parameters that are chosen. Stochastic ensembling is suggested as a method for improving the stability while training GANs.


title: On the formation of neutron stars via accretion-induced collapse in binaries

abstract: We investigate evolutionary pathways leading to neutron star formation through the collapse of oxygen-neon white dwarf (ONe WD) stars in interacting binaries. We consider (1) non-dynamical mass transfer where an ONe WD approaches the Chandrasekhar mass leading to accretion-induced collapse (AIC) and (2) dynamical timescale merger-induced collapse (MIC) between an ONe WD and another WD. We present rates, delay times, and progenitor properties for two different treatments of common envelope evolution. We show that AIC neutron stars are formed via many different channels and the most dominant channel depends on the adopted common envelope physics. Most AIC and MIC neutron stars are born shortly after star formation, though some have delay times >10 Gyr. The shortest delay time (25-50 Myr) AIC neutron stars have stripped-envelope, compact, helium-burning star donors, though many prompt AIC neutron stars form via wind-accretion from an asymptotic giant branch star. The longest delay time AIC neutron stars, which may be observed as young milli-second pulsars among globular clusters, have a red giant or main sequence donor at the time of NS formation and will eventually evolve into NS + helium WD binaries. We discuss AIC & MIC binaries as potential gravitational wave sources for LISA. Neutron stars created via AIC undergo a LMXB phase, offering an electromagnetic counterpart for those shortest orbital period sources that LISA could identify. The formation of neutron stars from interacting WDs in binaries is likely to be a key mechanism for the production of LIGO/Virgo gravitational wave sources (NS-NS and BH-NS mergers) in globular clusters.


title: Pixel-Level Alignment of Facial Images for High Accuracy Recognition Using Ensemble of Patches

abstract: The variation of pose, illumination and expression makes face recognition still a challenging problem. As a pre-processing in holistic approaches, faces are usually aligned by eyes. The proposed method tries to perform a pixel alignment rather than eye-alignment by mapping the geometry of faces to a reference face while keeping their own textures. The proposed geometry alignment not only creates a meaningful correspondence among every pixel of all faces, but also removes expression and pose variations effectively. The geometry alignment is performed pixel-wise, i.e., every pixel of the face is corresponded to a pixel of the reference face. In the proposed method, the information of intensity and geometry of faces are separated properly, trained by separate classifiers, and finally fused together to recognize human faces. Experimental results show a great improvement using the proposed method in comparison to eye-aligned recognition. For instance, at the false acceptance rate of 0.001, the recognition rates are respectively improved by 24% and 33% in Yale and AT&T datasets. In LFW dataset, which is a challenging big dataset, improvement is 20% at FAR of 0.1.


title: Approximate Spacetime for Neutron Stars

abstract: An approximate realistic metric representing the spacetime of neutron stars is obtained by perturbing the Kerr metric. This metric has five parameters, namely the mass, spin or angular momentum, mass quadrupole, spin octupole and mass hexadecapole. Moreover, a version of the Hartle-Thorne metric containing these parameters is constructed by means of a series transformation between these spacetimes and solving the Einstein field equations. The form of the Pappas metric in Schwarzschild spherical coordinates is found. The three relativistic multipole structures are compared.


title: Tropicalized quartics and canonical embeddings for tropical curves of genus 3

abstract: Brodsky, Joswig, Morrison and Sturmfels showed that not all abstract tropical curves of genus $3$ can be realized as a tropicalization of a quartic in the euclidean plane. In this article, we focus on the interior of the maximal cones in the moduli space and classify all curves which can be realized as a faithful tropicalization in a tropical plane. Reflecting the algebro-geometric world, we show that these are exactly those which are not realizably hyperelliptic. Our approach is constructive: For any not realizably hyperelliptic curve, we explicitly construct a realizable model of the tropical plane and a faithfully tropicalized quartic in it. These constructions rely on modifications resp. tropical refinements. Conversely, we prove that any realizably hyperelliptic curve cannot be embedded in such a fashion. For that, we rely on the theory of tropical divisors and embeddings from linear systems, and recent advances in the realizability of sections of the tropical canonical divisor.


title: Transverse Single Spin Asymmetry in $J/\psi$ Production

abstract: We estimate transverse single spin asymmetry (TSSA) in electroproduction of $J/\psi$ for J-Lab and EIC energies. We present estimates of TSSAs in $J/\psi$ production within generalized parton model (GPM) using recent parametrizations of gluon Sivers function (GSF) and compare the results obtained using color singlet model (CSM) with those obtained using color evaporation model (CEM) of quarkonium production.


title: Heavy flavour energy loss from AdS/CFT: A novel diffusion coefficient

abstract: Two AdS/CFT based energy loss models are used to compute the suppression and azimuthal correlations of heavy quarks in heavy ion collisions. The model with a velocity independent diffusion coefficient is in good agreement with B and D meson data up to high $p_T$. The partonic azimuthal correlations we calculate exhibit an order of magnitude difference in low momentum correlations to pQCD calculations [arXiv:1305.3823]. We thus propose heavy flavour momentum correlations as a distinguishing observable of weakly- and strongly-coupled energy loss mechanisms.


title: Holomorphic polygons and the bordered Heegaard Floer homology of link complements

abstract: We describe the construction of an $\mathcal{A}_\infty$ multi-module in terms of counts of holomorphic polygons in a series Heegaard multi-diagrams. We show that this is quasi-isomorphic to the type-A bordered-sutured invariant of a link complement with a view to calculating, in the sequel, these invariants in terms of the link Floer homology of the corresponding link.


title: What do we talk about when we speak of cosmological redshift?

abstract: From the first observations made by Slipher, our understanding and interpretation of the cosmological redshift was evolving until reaching the current consensus, through the expanding universe and the emergence of modern physical cosmology within the framework of General Relativity. The redshift is one of the most basic concepts of astronomy, and is one of the few observational parameters that can be measured directly. To refer to the temporal evolution of objects or cosmic structures in the universe, we often do so indistinctly through cosmic time or cosmological redshift. But repeatedly this connection ends up generating confusion not only among popular science communicators but also within the professional astronomical community. In this article, we will make a pedagogical approach to the link between cosmic time and cosmological redshift, and we will also clarify several common misunderstandings around this relation.


title: Second-order QCD effects in Higgs boson production through vector boson fusion

abstract: We compute the factorising second-order QCD corrections to the electroweak production of a Higgs boson through vector boson fusion. Our calculation is fully differential in the kinematics of the Higgs boson and of the final state jets, and uses the antenna subtraction method to handle infrared singular configurations in the different parton-level contributions. Our results allow us to reassess the impact of the next-to-leading order (NLO) QCD corrections to electroweak Higgs-plus-three-jet production and of the next-to-next-to-leading order (NNLO) QCD corrections to electroweak Higgs-plus-two-jet production. The NNLO corrections are found to be limited in magnitude to around $\pm 5\%$ and are uniform in several of the kinematical variables, displaying a kinematical dependence only in the transverse momenta and rapidity separation of the two tagging jets.


title: Coupling between Spin and Charge Order Driven by Magnetic Field in Triangular Ising System LuFe2O4+{\delta}

abstract: We present a study of the magnetic-field effect on spin correlations in the charge ordered triangular Ising system LuFe2O4+{\delta} through single crystal neutron diffraction. In the absence of a magnetic field, the strong diffuse neutron scattering observed below the Neel temperature (TN = 240 K) indicates that LuFe2O4+{\delta} shows short-range, two-dimensional (2D) correlations in the FeO5 triangular layers, characterized by the development of a magnetic scattering rod along the 1/3 1/3 L direction, persisting down to 5 K. We also found that on top of the 2D correlations, a long range ferromagnetic component associated with the propagation vector k1 = 0 sets in at around 240 K. On the other hand, an external magnetic field applied along the c-axis effectively favours a three-dimensional (3D) spin correlation between the FeO5 bilayers evidenced by the increase of the intensity of satellite reflections with propagation vector k2 = (1/3, 1/3, 3/2). This magnetic modulation is identical to the charge ordered superstructure, highlighting the field-promoted coupling between the spin and charge degrees of freedom. Formation of the 3D spin correlations suppresses both the rod-type diffuse scattering and the k1 component. Simple symmetry-based arguments provide a natural explanation of the observed phenomenon and put forward a possible charge redistribution in the applied magnetic field.


title: Field extensions, Derivations, and Matroids over Skew Hyperfields

abstract: We show that a field extension $K\subseteq L$ in positive characteristic $p$ and elements $x_e\in L$ for $e\in E$ gives rise to a matroid $M^\sigma$ on ground set $E$ with coefficients in a certain skew hyperfield $L^\sigma$. This skew hyperfield $L^\sigma$ is defined in terms of $L$ and its Frobenius action $\sigma:x\mapsto x^p$. The matroid underlying $M^\sigma$ describes the algebraic dependencies over $K$ among the $x_e\in L$ , and $M^\sigma$ itself comprises, for each $m\in \mathbb{Z}^E$, the space of $K$-derivations of $K\left(x_e^{p^{m_e}}: e\in E\right)$. The theory of matroid representation over hyperfields was developed by Baker and Bowler for commutative hyperfields. We partially extend their theory to skew hyperfields. To prove the duality theorems we need, we use a new axiom scheme in terms of quasi-Pl\"ucker coordinates.


title: Analytical Noise Bias Correction for Weak Lensing Shear Analysis with ERA

abstract: Highly precise weak lensing shear measurement is required for statistical weak gravitational lensing analysis such as cosmic shear measurement to achieve severe constraint on the cosmological parameters. For this purpose, the accurate shape measurement of background galaxies is absolutely important in which any systematic error in the measurement should be carefully corrected. One of the main systematic error comes from photon noise which is Poisson noise of flux from the atmosphere(noise bias). We investigate how the photon noise makes a systematic error in shear measurement within the framework of ERA method we developed in earlier papers and gives a practical correction method. The method is tested by simulations with real galaxy images with various conditions and it is confirmed that it can correct $80 \sim 90\%$ of the noise bias except for galaxies with very low signal to noise ratio.


title: Representation Theory of $\mathbb{Z}_2^{*n}$

abstract: We study the representations of the group $\mathbb{Z}_2^{*n}$, the free product of $\mathbb{Z}_2$ with itself $n$-times. We use the action of $B_n = S_2 \wr S_n $ as algebra automorphisms on the group algebra $\mathbb{C}(\mathbb{Z}_2^{*n})$ to find the components that contain simple representations and to study smoothness of their GIT-quotients. In particular, all the possible local quiver settings are studied for the component containing the standard $n$-dimensional representation of $S_{n+1}$.


title: Latitudinal structure and dynamic of the photospheric magnetic field

abstract: Analysis of the structure and dynamics of the magnetic field of the Sun is fundamental for understanding of the origin of solar activity and variability as well as for the study of solar-terrestrial relations. Observations of the large scale magnetic field in the photosphere taken at the Wilcox Solar Observatory from 1976 up to 2007 have been analysed to deduce its latitudinal and longitudinal structures, its differential rotation, and their variability in time. This paper is dedicated to the analysis and dynamics of the latitudinal structure of the solar magnetic field over three solar cycles 21, 22, 23. The main results discussed in this paper are the following: the large scale latitudinal structure is antisymmetric and composed of four zones with boundaries located at the equator, -25 and + 25 degrees, stable over 10-11 years with a time delay of about 5-6 years in near-equatorial zones. The variability and North-South asymmetry of polarity waves running from the equator to the poles with 2-3 - year period was studied in detail.


title: A Novel Supergeometric Generalization of Grassmannians

abstract: A new generalization of Grassmannians to supergeometry, different from the well known supergrassmannian, is introduced. These are constructed by gluing a finite number of copies of a \nu\- domain, i.e. a superdomain with an odd involution, say \nu\, on their structure sheaf considered as a sheaf of C^\infty_{R^m}-modules.


title: Some Properties of Fibonacci-Sum Set-Graphs

abstract: In this paper we study some properties of Fibonacci-sum set-graphs. The aforesaid graphs are an extension of the notion of Fibonacci-sum graphs to the notion of set-graphs. The colouring of Fibonacci-sum graphs is also discussed. A number of challenging research problems are posed in the conclusion.


title: A Posteriori Error Estimates for Non-Stationary Non-Linear Convection-Diffusion Equations

abstract: Motivated by stochastic convection-diffusion problems we derive a posteriori error estimates for non-stationary non-linear convection-diffusion equations acting as a deterministic paradigm. The problem considered here neither fits into the standard linear framework due to its non-linearity nor into the standard non-linear framework due to the lacking differentiability of the non-linearity. Particular attention is paid to the interplay of the various parameters controlling the relative sizes of diffusion, convection, reaction, and non-linearity (noise).


title: Markov spectrum near Freiman's isolated points in $M\setminus L$

abstract: Freiman proved in 1968 that the Lagrange and Markov spectra do not coincide by exhibiting a countable infinite collection $\mathcal{F}$ of isolated points of the Markov spectrum which do not belong the Lagrange spectrum. In this paper, we describe the structure of the elements of the Markov spectrum in the largest interval $(c_{\infty}, C_{\infty})$ containing $\mathcal{F}$ and avoiding the Lagrange spectrum. In particular, we compute the smallest known element $f$ of $M\setminus L$, and we show that the Hausdorff dimension of the portion of the Markov spectrum between $c_{\infty}$ and $C_{\infty}$ is $> 0.2628$.


title: Explaining a changeover from normal to super diffusion in time-dependent billiards

abstract: The changeover from normal to super diffusion in time dependent billiards is explained analytically. The unlimited energy growth for an ensemble of bouncing particles in time dependent billiards is obtained by means of a two dimensional mapping of the first and second moments of the velocity distribution function. We prove that for low initial velocities the mean velocity of the ensemble grows with exponent ~1/2 of the number of collisions with the border, therefore exhibiting normal diffusion. Eventually, this regime changes to a faster growth characterized by an exponent ~1 corresponding to super diffusion. For larger initial velocities, the temporary symmetry in the diffusion of velocities explains an initial plateau of the average velocity.


title: Ordinary $GL_2(F)$-representations in characteristic two via affine Deligne-Lusztig constructions

abstract: The group $\GL_2$ over a local field with (residue) characteristic $2$ possesses much more smooth supercuspidal $\ell$-adic representations, than over a local field of residue characteristic $> 2$. One way to construct these representations is via the theory of types of Bushnell-Kutzko. We construct many of them in the cohomology of certain extended affine Deligne-Lusztig varieties attached to $\GL_2$ and wildly ramified maximal tori in it. Then we compare our construction with the type-theoretic one. The corresponding extended affine Deligne-Lusztig varieties were introduced in a preceding article. Also in the present case they turn out to be zero-dimensional.


title: Asymptotically Locally Euclidean/Kaluza-Klein Stationary Vacuum Black Holes in 5 Dimensions

abstract: We produce new examples, both explicit and analytical, of bi-axisymmetric stationary vacuum black holes in 5 dimensions. A novel feature of these solutions is that they are asymptotically locally Euclidean in which spatial cross-sections at infinity have lens space $L(p,q)$ topology, or asymptotically Kaluza-Klein so that spatial cross-sections at infinity are topologically $S^1\times S^2$. These are nondegenerate black holes of cohomogeneity 2, with any number of horizon components, where the horizon cross-section topology is any one of the three admissible types: $S^3$, $S^1\times S^2$, or $L(p,q)$. Uniqueness of these solutions is also established. Our method is to solve the relevant harmonic map problem with prescribed singularities, having target symmetric space $SL(3,\mathbb{R})/SO(3)$. In addition, we analyze the possibility of conical singularities and find a large family for which geometric regularity is guaranteed.


title: Group theoretical formulation of free fall and projectile motion

abstract: In this work we formulate the group theoretical description of free fall and projectile motion. We show that the kinematic equations for constant acceleration form a one parameter group acting on a phase space. We define the group elements $\phi_t$ by their action on the points in the phase space. We also generalize this approach to projectile motion. We evaluate the group orbits regarding their relations to the physical orbits of particles and unphysical solutions. We note that the group theoretical formulation does not apply to more general cases involving a time dependent acceleration. This method improves our understanding of the constant acceleration problem with its global approach. It is especially beneficial for students who want to pursue a career in theoretical physics.


title: Vortex Solution of the Gravitational Field Equation of a Twisted Skyrme Strings

abstract: We construct non-linear sigma model plus Skyrme term (Skyrme model) with a twist in the gravitational field. To simplify the solution, first we examine non-linear sigma model without Skyrme term, in particular with a twist, which comprises a vortex solution with an added dependence on a twist term $mkz$, where $z$ is the vertical coordinate. We find that vortex solution for non-linear sigma model with a twist is similar with vortex solution without a twist. The work is still progress.


title: Upscaling Singular Sources in Weighted Sobolev Spaces by Sub-Grid Corrections

abstract: In this paper, we develop a numerical multiscale method to solve elliptic boundary value problems with heterogeneous diffusion coefficients and with singular source terms. When the diffusion coefficient is heterogeneous, this adds to the computational costs, and this difficulty is compounded by a singular source term. For singular source terms, the solution does not belong to the Sobolev space $H^1$, but to the space $W^{1,p}$ for some $p<2$. Hence, the problem may be reformulated in a distance-weighted Sobolev space. Using this formulation, we develop a method to upscale the multiscale coefficient near the singular sources by incorporating corrections into the coarse-grid. Using a sub-grid correction method, we correct the basis functions in a distance-weighted Sobolev space and show that these corrections can be truncated to design a computationally efficient scheme with optimal convergence rates. Due to the nature of the formulation in weighted spaces, the variational form must be posed on the cross product of complementary spaces. Thus, two such sub-grid corrections must be computed, one for each multiscale space of the cross product. A key ingredient of this method is the use of quasi-interpolation operators to construct the fine scale spaces. Therefore, we develop a weighted projective quasi-interpolation that can be used for a general class of Muckenhoupt weight functions. We verify the optimal convergence of the method in some numerical examples with singular point sources and line fractures, and with oscillatory and heterogeneous diffusion coefficients.


title: Rotation of the photospheric magnetic field through solar cycles 21, 22, 23

abstract: Rotation of the large scale solar magnetic field has a great importance for the understanding of solar dynamic, for the search of longitidinal structure and for the study of solar-terrestrial relations. 30-year long observations taken at the Wilcox Solar Observatory (USA) in 21-23 cycles have been analysed carefully to deduce magnetic field rotation rate at different latitudes in both hemispheres and its variability in time. The WSO data appear to indicate that additionally to the differential rotation along the latitudes there are running waves of fast rotation of the magnetic field. These torsional waves are running from the poles to the equator with a period of 11 years. The rotation of the magnetic field (RMF) is almost rigid at latitudes above 55 degrees in both hemispheres. The rotation rate in the sub-polar regions is slower when the magnetic field is strong there (during minima of solar activity), and faster when the magnetic field changes polarity (during maxima of solar activity).


title: Systematic Analysis of Low/Hard State RXTE Spectra of GX 339-4 to Constrain the Geometry of the System

abstract: One of the popular models for the low/hard state of Black Hole Binaries is that the standard accretion disk is truncated and the hot inner region produces via Comptonization, the hard X-ray flux. This is supported by the value of the high energy photon index, which is often found to be small $\sim$ 1.7 ($<$ 2) implying that the hot medium is seed photons starved. On the other hand, the suggestive presence of a broad relativistic Fe line during the hard state would suggest that the accretion disk is not truncated but extends all the way to the inner most stable circle orbit. In such a case, it is a puzzle why the hot medium would remain photon starved. The broad Fe line should be accompanied by a broad smeared reflection hump at $\sim$ 30 keV and it may be that this additional component makes the spectrum hard and the intrinsic photon index is larger, i.e. $>$ 2. This would mean that the medium is not photon deficient, reconciling the presence of a broad Fe line in the observed hard state. To test this hypothesis, we have analyzed the RXTE observations of GX 339-4 from the four outbursts during 2002-2011 and identify the observations when the system was in the hard state and showed a broad Fe line. We have then attempted to fit these observations with models, which include smeared reflection to understand whether the intrinsic photon index can indeed be large. We find that, while for some observations the inclusion of reflection does increase the photon index, there are hard state observations with broad Fe line that have photon indices less than 2.


title: Second order differentiation formula on $RCD^*(K,N)$ spaces

abstract: Aim of this paper is to prove the second order differentiation formula for $H^{2,2}$ functions along geodesics in $RCD^*(K,N)$ spaces with $N < \infty$. This formula is new even in the context of Alexandrov spaces, where second order differentiation is typically related to semiconvexity. We establish this result by showing that $W_2$-geodesics can be approximated up to second order, in a sense which we shall make precise, by entropic interpolation. In turn this is achieved by proving new, even in the smooth setting, estimates concerning entropic interpolations which we believe are interesting on their own. In particular we obtain: - equiboundedness of the densities along the entropic interpolations, - local equi-Lipschitz continuity of the Schr\"odinger potentials, - a uniform weighted $L^2$ control of the Hessian of such potentials. Finally, the techniques adopted in this paper can be used to show that in the $RCD$ setting the viscous solution of the Hamilton-Jacobi equation can be obtained via a vanishing viscosity method, in accordance with the smooth case. With respect to a previous version, where the space was assumed to be compact, in this paper the second order differentiation formula is proved in full generality.


title: Improved Time of Arrival measurement model for non-convex optimization with noisy data

abstract: The quadratic system provided by the Time of Arrival technique can be solved analytical or by optimization algorithms. In real environments the measurements are always corrupted by noise. This measurement noise effects the analytical solution more than non-linear optimization algorithms. On the other hand it is also true that local optimization tends to find the local minimum, instead of the global minimum. This article presents an approach how this risk can be significantly reduced in noisy environments. The main idea of our approach is to transform the local minimum to a saddle point, by increasing the number of dimensions.


title: Phantom fluid supporting traversable wormholes in alternative gravity with extra material terms

abstract: Wormholes are tunnels connecting different regions in space-time. They were obtained originally as a solution for Einstein's General Relativity theory and according to this theory they need to be filled by an exotic kind of anisotropic matter. In the present sense, by "exotic matter" we mean matter that does not satisfy the energy conditions. In this article we propose the modelling of wormholes within an alternative gravity theory that proposes an extra material (rather than geometrical) term in its gravitational action. Our solutions are obtained from well-known particular cases of the wormhole metric potentials, named redshift and shape functions, and yield the wormholes to be filled by a phantom fluid, that is, a fluid with equation of state parameter $\omega<-1$. In possession of the solutions for the wormhole material content, we also apply the energy conditions to them. The features of those are carefully discussed.


title: Three Random Intercepts of a Segment

abstract: We construct random triangles via uniform sampling of certain families of lines in the plane. Two examples are given. The word "uniform" turns out to be vague; two competing models are examined. Everything we write is well-known to experts. Which model is more appropriate? Our hope is to engage a larger audience in answering this question.


title: Mixtures of Factor Analyzers with Fundamental Skew Symmetric Distributions

abstract: Mixtures of factor analyzers (MFA) provide a powerful tool for modelling high-dimensional datasets. In recent years, several generalizations of MFA have been developed where the normality assumption of the factors and/or of the errors was relaxed to allow for skewness in the data. However, due to the form of the adopted component densities, the distribution of the factors/errors in most of these models is typically limited to modelling skewness oncentrated in a single direction. Here, we introduce a more flexible finite mixture of factor analyzers based on the class of scale mixtures of canonical fundamental skew normal (SMCFUSN) distributions. This very general class of skew distributions can capture various types of skewness and asymmetry in the data. In particular, the proposed mixture model of SMCFUSN factor analyzers(SMCFUSNFA) can simultaneously accommodate multiple directions of skewness. As such, it encapsulates many commonly used models as special and/or limiting cases, such as models of some versions of skew normal and skew t-factor analyzers, and skew hyperbolic factor analyzers. For illustration, we focus on the t-distribution member of the class of SMCFUSN distributions, leading to mixtures of canonical fundamental skew t-factor analyzers (CFUSTFA). Parameter estimation can be carried out by maximum likelihood via an EM-type algorithm. The usefulness and potential of the proposed model are demonstrated using two real datasets.


title: Efficient Learning of Bounded-Treewidth Bayesian Networks from Complete and Incomplete Data Sets

abstract: Learning a Bayesian networks with bounded treewidth is important for reducing the complexity of the inferences. We present a novel anytime algorithm (k-MAX) method for this task, which scales up to thousands of variables. Through extensive experiments we show that it consistently yields higher-scoring structures than its competitors on complete data sets. We then consider the problem of structure learning from incomplete data sets. This can be addressed by structural EM, which however is computationally very demanding. We thus adopt the novel k-MAX algorithm in the maximization step of structural EM, obtaining an efficient computation of the expected sufficient statistics. We test the resulting structural EM method on the task of imputing missing data, comparing it against the state-of-the-art approach based on random forests. Our approach achieves the same imputation accuracy of the competitors, but in about one tenth of the time. Furthermore we show that it has worst-case complexity linear in the input size, and that it is easily parallelizable.


title: A complete framework for linear filtering of bivariate signals

abstract: A complete framework for the linear time-invariant (LTI) filtering theory of bivariate signals is proposed based on a tailored quaternion Fourier transform. This framework features a direct description of LTI filters in terms of their eigenproperties enabling compact calculus and physically interpretable filtering relations in the frequency domain. The design of filters exhibiting fondamental properties of polarization optics (birefringence, diattenuation) is straightforward. It yields an efficient spectral synthesis method and new insights on Wiener filtering for bivariate signals with prescribed frequency-dependent polarization properties. This generic framework facilitates original descriptions of bivariate signals in two components with specific geometric or statistical properties. Numerical experiments support our theoretical analysis and illustrate the relevance of the approach on synthetic data.


title: On the gaps between consecutive primes

abstract: Let $p_n$ denote the $n$-th prime. For any $m\geq 1$, there exist infinitely many $n$ such that $p_{n}-p_{n-m}\leq C_m$ for some large constant $C_m>0$, and $$p_{n+1}-p_n\geq \frac{c_m\log n\log\log n\log\log\log\log n}{\log\log\log n}, $$ for some small constant $c_m>0$. Furthermore, we also obtain a related result concerning the least primes in arithmetic progressions.


title: Local Self-Energies for V and Pd Emergent from a Non-Local LDA+FLEX Implementation

abstract: In the spirit of recently developed LDA+U and LDA+DMFT methods we implement a combination of density functional theory in its local density approximation (LDA) with a $k$- and $\omega -$dependent self-energy found from diagrammatic fluctuational exchange (FLEX) approximation. The active Hilbert space here is described by the correlated subset of electrons which allows to tremendously reduce the sizes of matrices needed to represent charge and spin susceptibilities. The method is perturbative in nature but accounts for both bubble and ladder diagrams and accumulates the physics of momentum resolved spin fluctuations missing in such popular approach as GW. As an application, we study correlation effects on band structures in V and Pd. The d-electron self-energies emergent from this calculation are found to be remarkably k-independent. However, when we compare our calculated electronic mass enhancements against LDA+DMFT, we find that for a long standing problem of spin fluctuations in Pd, LDA+FLEX delivers a better agreement with experiment, although this conclusion depends on a particular value of Hubbard $U$ used in the simulation. We also discuss outcomes of recently proposed combinations of k-dependent FLEX with DMFT.


title: Forceless Sadowsky strips are spherical

abstract: We show that thin rectangular ribbons, defined as energy-minimising configurations of the Sadowsky functional for narrow developable elastic strips, have a propensity to form spherical shapes in the sense that forceless solutions lie on a sphere. This has implications for ribbonlike objects in (bio)polymer physics and nanoscience that cannot be described by the classical wormlike chain model. A wider class of functionals with this property is identified.


title: Radiating solitary waves in coupled Boussinesq equations

abstract: In this paper we are consider radiating solitary wave solutions of coupled regularised Boussinesq equations. This type of solution consists of a leading solitary wave with a small-amplitude co-propagating oscillatory tail, and emerges from a pure solitary wave solution of a symmetric reduction of the full system. We construct an asymptotic solution, where the leading order approximation in both components is obtained as a particular solution of the regularised Boussinesq equations in the symmetric case. At the next order, the system uncouples into two linear non-homogeneous ordinary differential equations with variable coefficients, one correcting the localised part of the solution, which we find analytically, and the other describing the co-propagating oscillatory tail. This latter equation is a fourth-order ordinary differential equation and is solved approximately by two different methods, each exploiting the assumption that the leading solitary wave has a small amplitude, and thus enabling an explicit estimate for the amplitude of the oscillating tail. These estimates are compared with corresponding numerical simulations.


title: High-level python abstractions for optimal checkpointing in inversion problems

abstract: Inversion and PDE-constrained optimization problems often rely on solving the adjoint problem to calculate the gradient of the objec- tive function. This requires storing large amounts of intermediate data, setting a limit to the largest problem that might be solved with a given amount of memory available. Checkpointing is an approach that can reduce the amount of memory required by redoing parts of the computation instead of storing intermediate results. The Revolve checkpointing algorithm o ers an optimal schedule that trades computational cost for smaller memory footprints. Integrat- ing Revolve into a modern python HPC code and combining it with code generation is not straightforward. We present an API that makes checkpointing accessible from a DSL-based code generation environment along with some initial performance gures with a focus on seismic applications.


title: Out-of-Band Radiation from Antenna Arrays Clarified

abstract: Non-linearities in radio-frequency (RF) transceiver hardware, particularly in power amplifiers, cause distortion in-band and out-of-band. Contrary to claims made in recent literature, in a multiple-antenna system this distortion is correlated across the antennas in the array. A significant implication of this fact is that out-of-band emissions caused by non-linearities are beamformed, in some cases into the same direction as the useful signal.


title: Order continuous operators on pre-Riesz spaces and embeddings

abstract: We investigate properties of order continuous operators on pre-Riesz spaces with respect to the embedding of the range space into a vector lattice cover or, in particular, into its Dedekind completion. We show that order continuity is preserved under this embedding for positive operators, but not in general. For the vector lattice $\ell_0^\infty$ of eventually constant sequences, we consider the pre-Riesz space of regular operators on $\ell_0^\infty$ and show that making the range space Dedekind complete does not provide a vector lattice cover of the pre-Riesz space. A similar counterexample is obtained for the directed part of the space of order continuous operators on $\ell_0^\infty$.


title: Common solar wind drivers behind magnetic storm-magnetospheric substorm dependency

abstract: The dynamical relationship between magnetic storms and magnetospheric substorms presents one of the most controversial problems of contemporary geospace research. Here, we tackle this issue by applying a causal inference approach to two corresponding indices in conjunction with several relevant solar wind variables. We demonstrate that the vertical component of the interplanetary magnetic field is the strongest and common driver of both, storms and substorms, and explains their the previously reported association. These results hold during both solar maximum and minimum phases and suggest that, at least based on the analyzed indices, there is no statistical evidence for a direct or indirect dependency between substorms and storms. A physical mechanism by which substorms drive storms or vice versa is, therefore, unlikely.


title: On the Stability of Independence Polynomials

abstract: The independence polynomial of a graph is the generating polynomial for the number of independent sets of each size, and its roots are called {\em independence roots}. We investigate the stability of such polynomials, that is, conditions under which the roots lie in the left half-plane (all of the real roots of independence polynomial are negative and hence lie in this half-plane). We show stability for all independence polynomials of graphs with independence number at most three, but for larger independence number we show that the independence polynomials can have roots arbitrarily far to the right. We provide families of graphs whose independence polynomials are stable and ones that are not, utilizing various graph operations.


title: Virtual photon polarization and dilepton anisotropy in relativistic nucleus-nucleus collisions

abstract: The polarization of virtual photons produced in relativistic nucleus-nucleus collisions provides information on the conditions in the emitting medium. In a hydrodynamic framework, the resulting angular anisotropy of the dilepton final state depends on the flow as well as on the transverse momentum and invariant mass of the photon. We illustrate these effects in dilepton production from quark-antiquark annihilation in the QGP phase and $\pi^+\pi^-$ annihilation in the hadronic phase for a static medium in global equilibrium and for a longitudinally expanding system.


title: Click Spam Prevention Model for On-Line Advertisement

abstract: This paper shows a vulnerability of the pay-per-click accounting of Google Ads and proposes a statistical tradeoff-based approach to manage this vulnerability. The result of this paper is a model to calculate the overhead cost per click necessary to protect the subscribers and a simple algorithm to implement this protection. Simulations validate the correctness of the model and the economical applicability.


title: Weighted Morrey spaces related to Schrodinger operators with potentials satisfying a reverse Holder inequality and fractional integrals

abstract: Let $\mathcal L=-\Delta+V$ be a Schr\"odinger operator on $\mathbb R^d$, $d\geq3$, where $\Delta$ is the Laplacian operator on $\mathbb R^d$ and the nonnegative potential $V$ belongs to the reverse H\"older class $RH_s$ for $s\geq d/2$. For given $0<\alpha<d$, the fractional integrals associated to the Schr\"odinger operator $\mathcal L$ is defined by $\mathcal I_{\alpha}={\mathcal L}^{-{\alpha}/2}$.Suppose that $b$ is a locally integrable function on $\mathbb R^d$, the commutator generated by $b$ and $\mathcal I_{\alpha}$ is defined by $[b,\mathcal I_{\alpha}]f(x)=b(x)\cdot \mathcal I_{\alpha}f(x)-\mathcal I_{\alpha}(bf)(x)$. In this paper, we first introduce some kinds of weighted Morrey spaces related to certain nonnegative potentials belonging to the reverse H\"older class $RH_s$ for $s\geq d/2$. Then we will establish the boundedness properties of the fractional integrals $\mathcal I_{\alpha}$ on these new spaces. Furthermore, weighted strong-type estimate for the corresponding commutator $[b,\mathcal I_{\alpha}]$ in the framework of Morrey spaces is also obtained. The classes of weights, the classes of symbol functions as well as weighted Morrey spaces discussed in this paper are larger than $A_{p,q}$, $\mathrm{BMO}(\mathbb R^d)$ and $L^{p,\kappa}(\mu,\nu)$ corresponding to the classical case (that is $V\equiv0$).


title: The paradigm of complexity. Contributions for hypertext's formal approaches

abstract: This article argues for a return to formal approaches of hypertext and builds on the paradigm of complexity to develop the idea of "hypermediator website". A hypermediator website is an intermediate device between a digitalization of book culture and a "real" hypertext writing. If our thinking on the hypermediator website joined the hypertext's notions and the databases, it differs by the relationship reader-device no longer based on information search query but using the visualization of the information.


title: Power-Controlled Hamiltonian Systems: Application to Electrical Systems with Constant Power Loads

abstract: We study a type of port-Hamiltonian system, in which the controller or disturbance is not applied to the flow variables, but to the systems power, a scenario that appears in many practical applications. A suitable framework is provided to model these systems and to investigate their shifted passivity properties, based on which, a stability analysis is carried out. The applicability of the results is illustrated with the important problem of stability analysis of electrical circuits with constant power loads.


title: Dark states in spin-polarized transport through triple quantum dot molecules

abstract: We study the spin-polarized transport through a triple quantum dot molecule weakly coupled to ferromagnetic leads. The analysis is performed by means of the real-time diagrammatic technique including up to the second order of perturbation expansion with respect to the tunnel coupling. The emphasis is put on the impact of dark states on spin-resolved transport characteristics. It is shown that the interplay of coherent population trapping and cotunneling processes results in a highly nontrivial behavior of the tunnel magnetoresistance, which can take negative values. Moreover, a super-Poissonian shot noise is found in transport regimes where the current is blocked by the formation of dark states, which can be additionally enhanced by spin-dependence of tunneling processes, depending on magnetic configuration of the device. The mechanisms leading to those effects are thoroughly discussed.


title: BROJA-2PID: A robust estimator for bivariate partial information decomposition

abstract: Makkeh, Theis, and Vicente found in [8] that Cone Programming model is the most robust to compute the Bertschinger et al. partial information decompostion (BROJA PID) measure [1]. We developed a production-quality robust software that computes the BROJA PID measure based on the Cone Programming model. In this paper, we prove the important property of strong duality for the Cone Program and prove an equivalence between the Cone Program and the original Convex problem. Then describe in detail our software and how to use it.\newline\indent


title: The field of quantum $GL(N,\mathbb{C})$ in the C$^*$-algebraic setting

abstract: Given a unital $*$-algebra $\mathscr{A}$ together with a suitable positive filtration of its set of irreducible bounded representations, one can construct a C$^*$-algebra $A_0$ with a dense two-sided ideal $A_c$ such that $\mathscr{A}$ maps into the multiplier algebra of $A_c$. When the filtration is induced from a central element in $\mathscr{A}$, we say that $\mathscr{A}$ is an s$^*$-algebra. We also introduce the notion of $\mathscr{R}$-algebra relative to a commutative s$^*$-algebra $\mathscr{R}$, and of Hopf $\mathscr{R}$-algebra. We formulate conditions such that the completion of a Hopf $\mathscr{R}$-algebra gives rise to a continuous field of Hopf C$^*$-algebras over the spectrum of $R_0$. We apply the general theory to the case of quantum $GL(N,\mathbb{C})$ as constructed from the FRT-formalism.


title: Invariant states on noncommutative tori

abstract: For any number $h$ such that $\hbar:=h/2\pi$ is irrational and any skew-symmetric, non-degenerate bilinear form $\sigma:\mathbb{Z}^{2g}\times \mathbb{Z}^{2g} \to \mathbb{Z}$, let be $\mathcal{A}^h_{g,\sigma}$ be the twisted group $*$-algebra $\mathbb{C}[\mathbb{Z}^{2g}]$ and consider the ergodic group of $*$-automorphisms of $\mathcal{A}^h_{g,\sigma}$ induced by the action of the symplectic group Sp$(\mathbb{Z}^{2g},\sigma)$. We show that the only Sp$(\mathbb{Z}^{2g},\sigma)$-invariant state on $\mathcal{A}^h_{g,\sigma}$ is the trace state $\tau$.


title: SCH-GAN: Semi-supervised Cross-modal Hashing by Generative Adversarial Network

abstract: Cross-modal hashing aims to map heterogeneous multimedia data into a common Hamming space, which can realize fast and flexible retrieval across different modalities. Supervised cross-modal hashing methods have achieved considerable progress by incorporating semantic side information. However, they mainly have two limitations: (1) Heavily rely on large-scale labeled cross-modal training data which are labor intensive and hard to obtain. (2) Ignore the rich information contained in the large amount of unlabeled data across different modalities, especially the margin examples that are easily to be incorrectly retrieved, which can help to model the correlations. To address these problems, in this paper we propose a novel Semi-supervised Cross-Modal Hashing approach by Generative Adversarial Network (SCH-GAN). We aim to take advantage of GAN's ability for modeling data distributions to promote cross-modal hashing learning in an adversarial way. The main contributions can be summarized as follows: (1) We propose a novel generative adversarial network for cross-modal hashing. In our proposed SCH-GAN, the generative model tries to select margin examples of one modality from unlabeled data when giving a query of another modality. While the discriminative model tries to distinguish the selected examples and true positive examples of the query. These two models play a minimax game so that the generative model can promote the hashing performance of discriminative model. (2) We propose a reinforcement learning based algorithm to drive the training of proposed SCH-GAN. The generative model takes the correlation score predicted by discriminative model as a reward, and tries to select the examples close to the margin to promote discriminative model by maximizing the margin between positive and negative data. Experiments on 3 widely-used datasets verify the effectiveness of our proposed approach.


title: Pressure-driven dynamics of liquid plugs in rectangular microchannels: influence of the transition between quasi-static and dynamic film deposition regimes

abstract: In this paper, we study experimentally and theoretically the dynamics of liquid plugs in rectangular microchannels for both unidirectional and cyclic pressure forcing. In both cases, it is shown that the transition between quasi-static and dynamic film deposition behind the liquid plug leads to a dramatic acceleration of the plug, rapidly leading to its rupture. This behaviour proper to channels with sharp corners is recovered from a reduced dimension model based on previous theoretical and numerical developments. In addition, it is shown for cyclic periodic forcing that the plug undergoes stable periodic oscillations if it remains in the quasi-static film deposition regime during the first cycle, while otherwise it accelerates cyclically and ruptures. The transition between these two regimes occurs at a pressure-dependent critical initial length, whose value can be predicted theoretically.


title: Conformal invariance and the Lundgren-Monin-Novikov equations for vorticity fields in 2D turbulence: Refuting a recent claim

abstract: The recent claim by Grebenev et al. [J. Phys. A: Math. Theor. 50, 435502 (2017)] that the inviscid 2D Lundgren-Monin-Novikov (LMN) equations on a zero vorticity characteristic naturally would reveal local conformal invariance when only analyzing these by means of a classical Lie-group symmetry approach, is invalid and will be refuted in the present comment. To note is that within this comment the (possible) existence of conformal invariance in 2D turbulence is not questioned, only the conclusion as is given in Grebenev et al. (2017) and their approach how this invariance was derived is what is being criticized and refuted herein. In fact, the algebraic derivation for conformal invariance of the 2D LMN vorticity equations in Grebenev et al. (2017) is flawed. A key constraint of the LMN equations has been wrongly transformed. Providing the correct transformation instead will lead to a breaking of the proclaimed conformal group. The corrected version of Grebenev et al. (2017) just leads to a globally constant scaling in the fields and not to a local one as claimed. In consequence, since in Grebenev et al. (2017) only the first equation within the infinite and unclosed LMN chain is considered, also different Lie-group infinitesimals for the one- and two-point probability density functions (PDFs) will result from this correction, replacing thus the misleading ones proposed.

categories: physics.flu-dyn cond-mat.stat-mech math-ph math.MP


title: The law of a point process of Brownian excursions in a domain is determined by the law of its trace

abstract: We show the result that is stated in the title of the paper, which has consequences about decomposition of Brownian loop-soup clusters in two dimensions.


title: $T$-odd correlations in polarized top quark decays in the sequential decay $t(\uparrow) \to X_b+W^+(\to \ell^+ + \nu_\ell)$ and in the quasi three-body decay $t(\uparrow) \to X_b+ \ell^+ + \nu_\ell$

abstract: We identify the $T$-odd structure functions that appear in the description of polarized top quark decays in the sequential decay $t(\uparrow) \to X_b+W^+(\to \ell^+ + \nu_\ell)$ (two structure functions) and the quasi-three-body decay $t(\uparrow) \to X_b+ \ell^+ + \nu_\ell$ (one structure function). A convenient measure of the magnitude of the $T$-odd structure functions is the contribution of the imaginary part Im $g_R$ of the right-chiral tensor coupling $g_R$ to the $T$-odd structure functions which we work out. Contrary to the case of QCD the NLO electroweak corrections to polarized top quark decays admit of absorptive one-loop vertex contributions. We analytically calculate the imaginary parts of the relevant four electroweak one-loop triangle vertex diagrams and determine their contributions to the $T$-odd helicity structure functions that appear in the description of polarized top quark decays.


title: Up to topological concordance links are strongly quasipositive

abstract: We generalize an algorithm of Rudolph to establish that every link is topologically concordant to a strongly quasipositive link.


title: Viscous effects on the dynamical evolution of QCD matter during the first-order confinement phase transition in heavy-ion collisions

abstract: We investigate viscous effects on the dynamical evolution of QCD matter during the first-order phase transition, which may happen in heavy-ion collisions. We first obtain the first-order phase transition line in the QCD phase diagram under the Gibbs condition by using the MIT bag model and the hadron resonance gas model for the equation of state of partons and hadrons. The viscous pressure, which corresponds to the friction in the energy balance, is then derived from the energy and net baryon number conservation during the phase transition. We find that the viscous pressure relates to the thermodynamic change of the two-phase state and thus affects the timescale of the phase transition. Numerical results are presented for demonstrations.


title: The merger of two compact stars: a tool for dense matter nuclear physics

abstract: We discuss the different signals, in gravitational and electromagnetic waves, emitted during the merger of two compact stars. We will focus in particular on the possible contraints that those signals can provide on the equation of state of dense matter. Indeed, the stiffness of the equation of state and the particle composition of the merging compact stars, strongly affect e.g. the life time of the post-merger remnant and its gravitational wave signal, the emission of the short gamma-ray-burst, the amount of ejected mass and the related kilonova. The first detection of gravitational waves from the merger of two compact stars in August 2017, GW170817, and the subsequent detections of its electromagnetic counterparts, GRB170817A and AT2017gfo, is the first example of the era of "multi-messenger astronomy": we discuss what we have learned from this detection on the equation of state of compact stars and we provide a tentative interpretation of this event, within the two families scenario, as due to the merger of a hadronic star with a quark star.


title: Defining the principles for the motion of blood through arteries

abstract: This is an annotated translation from Latin of 'Principia pro motu sanguinis per arterias determinando' in which Euler develops the first known work on the mechanics of flows in elastic tubes, intended to the first contest of the Dijon Academy of 1742. In this work, Euler applied the principles of mass conservation and momentum conservation to the one dimensional flow of an incompressible fluid through an elastic tube driven by a piston pump. These would allow the analytical treatment of circulatory physiology and hemodynamics, which still undergird the most advanced numerical methods in use today for blood flow analysis in arterial networks.


title: Privacy preserving clustering with constraints

abstract: The $k$-center problem is a classical combinatorial optimization problem which asks to find $k$ centers such that the maximum distance of any input point in a set $P$ to its assigned center is minimized. The problem allows for elegant $2$-approximations. However, the situation becomes significantly more difficult when constraints are added to the problem. We raise the question whether general methods can be derived to turn an approximation algorithm for a clustering problem with some constraints into an approximation algorithm that respects one constraint more. Our constraint of choice is privacy: Here, we are asked to only open a center when at least $\ell$ clients will be assigned to it. We show how to combine privacy with several other constraints.


title: Spectral Learning of Binomial HMMs for DNA Methylation Data

abstract: We consider learning parameters of Binomial Hidden Markov Models, which may be used to model DNA methylation data. The standard algorithm for the problem is EM, which is computationally expensive for sequences of the scale of the mammalian genome. Recently developed spectral algorithms can learn parameters of latent variable models via tensor decomposition, and are highly efficient for large data. However, these methods have only been applied to categorial HMMs, and the main challenge is how to extend them to Binomial HMMs while still retaining computational efficiency. We address this challenge by introducing a new feature-map based approach that exploits specific properties of Binomial HMMs. We provide theoretical performance guarantees for our algorithm and evaluate it on real DNA methylation data.


title: Coarse median algebras: The intrinsic geometry of coarse median spaces and their intervals

abstract: This paper establishes a new combinatorial framework for the study of coarse median spaces, bridging the worlds of asymptotic geometry, algebra and combinatorics. We introduce a simple and entirely algebraic notion of coarse median algebra which simultaneously generalises the concepts of bounded geometry coarse median spaces and classical discrete median algebras. We study the coarse median universe from the perspective of intervals, with a particular focus on cardinality as a proxy for distance. In particular we prove that the metric on a quasi-geodesic coarse median space of bounded geometry can be constructed up to quasi-isometry using only the coarse median operator. Finally we develop a concept of rank for coarse median algebras in terms of the geometry of intervals and show that the notion of finite rank coarse median algebra provides a natural higher dimensional analogue of Gromov's concept of $\delta$-hyperbolicity.


title: Cadre Modeling: Simultaneously Discovering Subpopulations and Predictive Models

abstract: We consider the problem in regression analysis of identifying subpopulations that exhibit different patterns of response, where each subpopulation requires a different underlying model. Unlike statistical cohorts, these subpopulations are not known a priori; thus, we refer to them as cadres. When the cadres and their associated models are interpretable, modeling leads to insights about the subpopulations and their associations with the regression target. We introduce a discriminative model that simultaneously learns cadre assignment and target-prediction rules. Sparsity-inducing priors are placed on the model parameters, under which independent feature selection is performed for both the cadre assignment and target-prediction processes. We learn models using adaptive step size stochastic gradient descent, and we assess cadre quality with bootstrapped sample analysis. We present simulated results showing that, when the true clustering rule does not depend on the entire set of features, our method significantly outperforms methods that learn subpopulation-discovery and target-prediction rules separately. In a materials-by-design case study, our model provides state-of-the-art prediction of polymer glass transition temperature. Importantly, the method identifies cadres of polymers that respond differently to structural perturbations, thus providing design insight for targeting or avoiding specific transition temperature ranges. It identifies chemically meaningful cadres, each with interpretable models. Further experimental results show that cadre methods have generalization that is competitive with linear and nonlinear regression models and can identify robust subpopulations.