title: Equilibration properties of small quantum systems: further examples

abstract: It has been proposed to investigate the equilibration properties of a small isolated quantum system by means of the matrix of asymptotic transition probabilities in some preferential basis. The trace $T$ of this matrix measures the degree of equilibration of the system prepared in a typical state of the preferential basis. This quantity may vary between unity (ideal equilibration) and the dimension $N$ of the Hilbert space (no equilibration at all). Here we analyze several examples of simple systems where the behavior of $T$ can be investigated by analytical means. We first study the statistics of $T$ when the Hamiltonian governing the dynamics is random and drawn from a distribution invariant under the group U$(N)$ or O$(N)$. We then investigate a quantum spin $S$ in a tilted magnetic field making an arbitrary angle with the preferred quantization axis, as well as a tight-binding particle on a finite electrified chain. The last two cases provide examples of the interesting situation where varying a system parameter -- such as the tilt angle or the electric field -- through some scaling regime induces a continuous transition from good to bad equilibration properties.


title: Spacings Around An Order Statistic

abstract: We determine the joint limiting distribution of adjacent spacings around a central, intermediate, or an extreme order statistic $X_{k:n}$ of a random sample of size $n$ from a continuous distribution $F$. For central and intermediate cases, normalized spacings in the left and right neighborhoods are asymptotically i.i.d. exponential random variables. The associated independent Poisson arrival processes are independent of $X_{k:n}$. For an extreme $X_{k:n}$, the asymptotic independence property of spacings fails for $F$ in the domain of attraction of Fr\'{e}chet and Weibull ($\alpha \neq 1$) distributions. This work also provides additional insight into the limiting distribution for the number of observations around $X_{k:n}$ for all three cases.


title: Efficient Large-scale Approximate Nearest Neighbor Search on the GPU

abstract: We present a new approach for efficient approximate nearest neighbor (ANN) search in high dimensional spaces, extending the idea of Product Quantization. We propose a two-level product and vector quantization tree that reduces the number of vector comparisons required during tree traversal. Our approach also includes a novel highly parallelizable re-ranking method for candidate vectors by efficiently reusing already computed intermediate values. Due to its small memory footprint during traversal, the method lends itself to an efficient, parallel GPU implementation. This Product Quantization Tree (PQT) approach significantly outperforms recent state of the art methods for high dimensional nearest neighbor queries on standard reference datasets. Ours is the first work that demonstrates GPU performance superior to CPU performance on high dimensional, large scale ANN problems in time-critical real-world applications, like loop-closing in videos.


title: Routes Obey Hierarchy in Complex Networks

abstract: Various hypotheses exist about the paths used for communication between the nodes of complex networks. Most studies simply suppose that communication goes via shortest paths, while others have more explicit assumptions about how routing (alternatively navigation or search) works or should work in real networks. However, these assumptions are rarely checked against real data. Here we directly analyze the structure of operational paths using real measurements. For this purpose we use existing and newly created datasets having both the topology of the network and a sufficient number of empirically-determined paths over it. Such datasets are processed for air transportation networks, the human brain, the Internet and the fit-fat-cat word ladder game. Our results suggest that from the great number of possible paths, nature seems to pick according to some simple rules, which we will refer to as routing policies. First we confirm, that the preference of short paths is an inevitable policy element, however the observed stretch of the paths suggests that there are other policies at work simultaneously. We identify two additional policies common in our networks: the "conform hierarchy", meaning that the paths should obey the structural hierarchy of the network, and the "prefer downstream" policy which promotes avoiding the network core if possible. Building upon these simple policies, we propose a synthetic routing policy which can recover the basic statistical properties of the operational paths in networks. Our results can be helpful in estimating the reaction of complex systems for stress coming from the outside more accurately than the shortest path assumption permits.


title: Why Do Cascade Sizes Follow a Power-Law?

abstract: We introduce random directed acyclic graph and use it to model the information diffusion network. Subsequently, we analyze the cascade generation model (CGM) introduced by Leskovec et al. [19]. Until now only empirical studies of this model were done. In this paper, we present the first theoretical proof that the sizes of cascades generated by the CGM follow the power-law distribution, which is consistent with multiple empirical analysis of the large social networks. We compared the assumptions of our model with the Twitter social network and tested the goodness of approximation.


title: A complete characterization on the robust isolated calmness of the nuclear norm regularized convex optimization problems

abstract: In this paper, we provide a complete characterization on the robust isolated calmness of the Karush-Kuhn-Tucker (KKT) solution mapping for convex constrained optimization problems regularized by the nuclear norm function. This study is motivated by the recent work in [8], where the authors show that under the Robinson constraint qualification at a local optimal solution, the KKT solution mapping for a wide class of conic programming problem is robustly isolated calm if and only if both the second order sufficient condition (SOSC) and the strict Robinson constraint qualification (SRCQ) are satisfied. Based on the variational properties of the nuclear norm function and its conjugate, we establish the equivalence between the primal/dual SOSC and the dual/primal SRCQ. The derived results lead to several equivalent characterizations of the robust isolated calmness of the KKT solution mapping and add insights to the existing literature on the stability of the nuclear norm regularized convex optimization problems.


title: Classical and semi-classical energy conditions

abstract: The standard energy conditions of classical general relativity are (mostly) linear in the stress-energy tensor, and have clear physical interpretations in terms of geodesic focussing, but suffer the significant drawback that they are often violated by semi-classical quantum effects. In contrast, it is possible to develop non-standard energy conditions that are intrinsically non-linear in the stress-energy tensor, and which exhibit much better well-controlled behaviour when semi-classical quantum effects are introduced, at the cost of a less direct applicability to geodesic focussing. In this article we will first review the standard energy conditions and their various limitations. (Including the connection to the Hawking--Ellis type I, II, III, and IV classification of stress-energy tensors). We shall then turn to the averaged, nonlinear, and semi-classical energy conditions, and see how much can be done once semi-classical quantum effects are included.


title: Comment: Energy Spectrum of a Graphene Quantum Dot in a Perpendicular Magnetic Field

abstract: In a recent comment (arXiv:1607.06081), Falaye et al. claim that there are certain flaws in our publication (Phys. Rev. B, 78, 195427 (2008)). We point out that our results, in particular the analytic derivation of the energy spectrum of a circular graphene quantum dot exposed to a perpendicular magnetic field, are correct and equivalent to the result of Falaye et al. A misleading notation error is corrected.


title: On a variable step size modification of Hines' method in computational neuroscience

abstract: For simulating large networks of neurons Hines proposed a method which uses extensively the structure of the arising systems of ordinary differential equations in order to obtain an efficient implementation. The original method requires constant step sizes and produces the solution on a staggered grid. In the present paper a one-step modification of this method is introduced and analyzed with respect to their stability properties. The new method allows for step size control. Local error estimators are constructed. The method has been implemented in matlab and tested using simple Hodgkin-Huxley type models. Comparisons with standard state-of-the-art solvers are provided.


title: Quantum corpuscular corrections to the Newtonian potential

abstract: We study an effective quantum description of the static gravitational potential for spherically symmetric systems up to the first post-Newtonian order. We start by obtaining a Lagrangian for the gravitational potential coupled to a static matter source from the weak field expansion of the Einstein-Hilbert action. By analysing a few classical solutions of the resulting field equation, we show that our construction leads to the expected post-Newtonian expressions. Next, we show that one can reproduce the classical Newtonian results very accurately by employing a coherent quantum state and modifications to include the first post-Newtonian corrections are considered. Our findings establish a connection between the corpuscular model of black holes and post-Newtonian gravity, and set the stage for further investigations of these quantum models.


title: Revisiting the Abraham-Minkowski Dilemma

abstract: Here, the Abraham-Minkowski controversy on the correct definition of the light momentum in a macroscopic medium is revisited with the purpose to highlight that an effective medium formalism necessarily restricts the available information on the internal state of a system, and that this is ultimately the reason why the dilemma has no universal solution. Despite these difficulties, it is demonstrated that in the limit of no material absorption and under steady-state conditions, the time-averaged light (kinetic) momentum may be unambiguously determined by the Abraham result, both for bodies at rest and for circulatory flows of matter. The implications of these findings are discussed in the context of quantum optics of moving media, and we examine in detail the fundamental role of the Minkowski momentum in such a context.


title: Broadening of the thermal component of the prompt GRB emission due to rapid temperature evolution

abstract: The observations of the prompt emission of gamma ray bursts (GRB) by GLAST Burst Monitor (GBM), on board Fermi Gamma-ray Space Telescope, suggest the presence of a significant thermal spectral component, whose origin is not well understood. Recently, it has been shown that for long duration GRBs, the spectral width as defined as the logarithm of the ratio of the energies at which the spectrum falls to half its peak value, lie in the range of 0.84-1.3 with a median value of 1.07. Thus, while most of the GRB spectra are found to be too narrow to be explained by synchrotron emission from an electron distribution, they are also significantly broader than a blackbody spectrum whose width should be 0.54. Here, we consider the possibility that an intrinsic thermal spectrum from a fire-ball like model, may be observed to be broadened if the system undergoes a rapid temperature evolution. We construct a toy-model to show that for bursts with durations in the range 5-70 seconds, the widths of their 1 second time-averaged spectra can be at the most $\lesssim 0.557$. Thus, while rapid temperature variation can broaden the detected spectral shape, the observed median value of $\sim 1.07$ requires that there must be significant sub-photospheric emission and/or an anisotropic explosion to explain the broadening for most GRB spectra.


title: A Mean-field Stochastic Control Problem with Partial Observations

abstract: In this paper we are interested in a new type of {\it mean-field}, non-Markovian stochastic control problems with partial observations. More precisely, we assume that the coefficients of the controlled dynamics depend not only on the paths of the state, but also on the conditional law of the state, given the observation to date. Our problem is strongly motivated by the recent study of the mean field games and the related McKean-Vlasov stochastic control problem, but with added aspects of path-dependence and partial observation. We shall first investigate the well-posedness of the state-observation dynamics, with combined reference probability measure arguments in nonlinear filtering theory and the Schauder fixed point theorem. We then study the stochastic control problem with a partially observable system in which the conditional law appears nonlinearly in both the coefficients of the system and cost function. As a consequence the control problem is intrinsically "time-inconsistent", and we prove that the Pontryagin Stochastic Maximum Principle holds in this case and characterize the adjoint equations, which turn out to be a new form of mean-field type BSDEs.


title: Variational Problems for Foppl-von Karman plates

abstract: Some variational problems for a Foppl-von Karman plate subject to general equilibrated loads are studied. The existence of global minimizers is proved under the assumption that the out-of-plane displacement fulfils homogeneous Dirichlet condition on the whole boundary while the in-plane displacement fulfils nonhomogeneous Neumann condition. If the Dirichlet condition is prescribed only on a subset of the boundary, then the energy may be unbounded from below over the set of admissible configurations, as shown by several explicit conterexamples: in these cases the analysis of critical points is addressed through an asymptotic development of the energy functional in a neighborhood of the flat configuration. By a $\Gamma$-convergence approach we show that critical points of the Foppl-von Karman energy can be strongly approximated by uniform Palais-Smale sequences of suitable functionals: this property leads to identify relevant features for critical points of approximating functionals, e.g. buckled configurations of the plate. Eventually we perform further analysis as the plate thickness tends to 0, by assuming that the plate is prestressed and the energy functional depends only on the transverse displacement around the given prestressed state: by this approach, first we identify suitable exponents of plate thickness for load scaling, then we show explicit asymptotic oscillating minimizers as a mechanism to relax compressive states in an annular plate.


title: On-chip linear and nonlinear control of single molecules coupled to a nanoguide

abstract: While experiments with one or two quantum emitters have become routine in various laboratories, scalable platforms for efficient optical coupling of many quantum systems remain elusive. To address this issue, we report on chip-based systems made of one-dimensional subwavelength dielectric waveguides (nanoguides) and polycyclic aromatic hydrocarbon molecules. After discussing the design and fabrication requirements, we present data on coherent linear and nonlinear spectroscopy of single molecules coupled to a nanoguide mode. Our results show that external microelectrodes as well as optical beams can be used to switch the propagation of light in a nanoguide via the Stark effect and a nonlinear optical process, respectively. The presented nanoguide architecture paves the way for the investigation of many-body phenomena and polaritonic states and can be readily extended to more complex geometries for the realization of quantum integrated photonic circuits.


title: On the dark radiation problem in the axiverse

abstract: String scenarios generically predict that we live in a so called axiverse: the Universe with about a hundred of light axion species which are decoupled from the Standard Model particles. However, the axions can couple to the inflaton which leads to their production after inflation. Then, these axions remain in the expanding Universe contributing to the dark radiation component, which is severely bounded from present cosmological data. We place a general constraint on the axion production rate and apply it to several variants of reasonable inflaton-to-axion couplings. The limit merely constrains the number of ultralight axions and the relative strength of inflaton-to-axion coupling. It is valid in both large and small field inflationary models irrespectively of the axion energy scales and masses. Thus, the limit is complementary to those associated with the Universe overclosure and axion isocurvature fluctuations. In particular, a hundred of axions is forbidden if inflaton universally couples to all the fields at reheating. In the case of gravitational sector being responsible for the reheating of the Universe (which is a natural option in all inflationary models with modified gravity), the axion production can be efficient. We find that in the Starobinsky $R^2$-inflation even a single axion (e.g. the standard QCD-axion) is in tension with the Planck data, making the model inconsistent with the axiverse. The general conclusion is that an inflation with inefficient reheating mechanism and low reheating temperature may be in tension with the presence of light scalars.


title: Integrable and non-integrable structures in Einstein - Maxwell equations with Abelian isometry group $\mathcal{G}_2$

abstract: The classes of electrovacuum Einstein - Maxwell fields (with a cosmological constant), which metrics admit an Abelian two-dimensional isometry group $\mathcal{G}_2$ with non-null orbits and electromagnetic fields possess the same symmetry, are considered. For these fields we describe the structures of so called "nondynamical degrees of freedom" which presence as well as the presence of a cosmological constant change (in a strikingly similar ways) the dynamical equations and destroy their known integrable structures. Modifications of the known reduced forms of Einstein - Maxwell equations -- the Ernst equations and self-dual Kinnersley equations in the presence of non-dynamical degrees of freedom are found and subclasses of fields with non-dynamical degrees of freedom are considered for : (I) vacuum metrics with cosmological constant, (II) vacuum space-times with isometry groups $\mathcal{G}_2$ which orbits do not admit the orthogonal 2-surfaces (none-orthogonally-transitive isometry groups) and (III) electrovacuum fields with more general structures of electromagnetic fields than in the known integrable cases. For these classes, in the case of diagonal metrics, the field equations can be reduced to the only nonlinear equation of the fourth order for one real function $\alpha(x^1,x^2)$ which characterises the element of area on the orbits of the isometry group $\mathcal{G}_2$ . Simple examples of solutions for each of these classes are presented. It is pointed out that if for some two-dimensional reduction of Einstein's field equations in four or higher dimensions, the function $\alpha(x^1,x^2)$ possess a "harmonic" structure, instead of being (together with other field variables) a solution of some nonlinear equations, this can be an indication of possible complete integrability of these reduced dynamical equations for the fields with vanishing of all non-dynamical degrees of freedom.


title: Radio Galaxies with the Cherenkov Telescope Array

abstract: Misaligned AGN (MAGNs), i.e., radio-loud AGNs with the jet not pointing directly towards us, represent a new class of GeV emitters revealed by the Fermi space telescope. Although they comprise only a small fraction of the high-energy sources, MAGNs are extremely interesting objects offering a different perspective to study high-energy processes with respect to blazars. The aim of this work is to evaluate the impact of the new-generation Cherenkov Telescope Array (CTA) on the MAGN class and propose possible observational strategies to optimize their detection.


title: Implementing Parrondo's paradox with two coin quantum walks

abstract: Parrondo's paradox is ubiquitous in games, ratchets and random walks.The apparent paradox, devised by J.~M.~R.~Parrondo, that two losing games $A$ and $B$ can produce an winning outcome has been adapted in many physical and biological systems to explain their working. However, proposals on demonstrating Parrondo's paradox using quantum walks failed {for large number of steps}. In this work, we show that instead of a single coin if we consider a two coin initial state which may or may not be entangled, we can observe a genuine Parrondo's paradox with quantum walks. Further we focus on reasons for this and pin down the asymmetry in initial two-coin state or asymmetry in shift operator, either of which are necessary for observing a genuine Parrondo's paradox. We extend our work to a 3-coin initial state too with similar results. The implications of our work for observing quantum ratchet like behavior using quantum walks is also discussed.


title: Study of corrections to the eikonal approximation

abstract: For the last decades, multiple international facilities have developed Radioactive-Ion Beams (RIB) to measure reaction processes including exotic nuclei. These measurements coupled with an accurate theoretical model of the reaction enable us to infer information about the structure of these nuclei. The partial-wave expansion provides a precise description of two-body collisions but has a large computational cost. To cope with this issue, the eikonal approximation is a powerful tool as it reduces the computational time and still describes the quantum effects observed in reaction observables. However, its range of validity is restricted to high energy and to forward scattering angles. In this work, we analyse the extension of the eikonal approximation to lower energies and larger angles through the implementation of two corrections. These aim to improve the treatment of the nuclear and Coulomb interactions within the eikonal model. The first correction is based on an expansion of the $T$-matrix while the second relies on a semi-classical approach. They permit to better account for the deflection of the projectile by the target, which is neglected in the standard eikonal model. The gain in accuracy of each correction is evaluated through the analyses of angular cross sections computed with the standard eikonal model, its corrections and the partial-wave expansion. These analyses have been performed for tightly bound projectiles ($^{10}\mathrm{Be}$) from intermediate energies (67 MeV/nucleon) down to energies of interest of future RIB facilities such as HIE-ISOLDE and ReA12 at MSU (10 MeV/nucleon).


title: Decay angular distributions of $K^*$ and $D^*$ vector mesons in pion-nucleon scattering

abstract: The production mechanisms of open strangeness ($K^*$) and open charm ($D^*$) vector mesons in $\pi^- p$ scattering, namely, $\pi^- +p \to K^{*0} + \Lambda$ and $\pi^- + p \to D^{*-} + \Lambda_c^+$, are investigated within the modified quark-gluon string model. In order to identify the major reaction mechanisms, we consider the subsequent decays of the produced vector mesons into two pseudoscalar mesons, i.e., $K^* \to K + \pi$ and $D^* \to D + \pi$. We found that the decay distributions and density matrix elements are sensitive to the production mechanisms and can be used to disentangle the vector trajectory and pseudoscalar trajectory exchange models. Our results for $K^*$ production are compared with the currently available experimental data and the predictions for $D^*$ production processes are presented as well. Our predictions can be tested at the present or planned experimental facilities.


title: Broadening of the Berezinskii-Kosterlitz-Thouless transition by correlated disorder

abstract: The Berezinskii-Kosterlitz-Thouless (BKT) transition in two-dimensional superconductors is usually expected to be protected against disorder. However, its typical signatures in real system, like e.g. the superfluid-density jump, are often at odd with this expectation. Here we show that the disorder-induced granularity of the superconducting state modifies the nucleation mechanism for vortex-antivortex pairs. This leads to a considerable smearing of the universal superfluid-density jump as compared to the paradigmatic clean case, in agreement with experimental observations.


title: The importance of stain normalization in colorectal tissue classification with convolutional networks

abstract: The development of reliable imaging biomarkers for the analysis of colorectal cancer (CRC) in hematoxylin and eosin (H&E) stained histopathology images requires an accurate and reproducible classification of the main tissue components in the image. In this paper, we propose a system for CRC tissue classification based on convolutional networks (ConvNets). We investigate the importance of stain normalization in tissue classification of CRC tissue samples in H&E-stained images. Furthermore, we report the performance of ConvNets on a cohort of rectal cancer samples and on an independent publicly available dataset of colorectal H&E images.


title: Robust and adaptive search

abstract: Binary search finds a given element in a sorted array with an optimal number of $\log n$ queries. However, binary search fails even when the array is only slightly disordered or access to its elements is subject to errors. We study the worst-case query complexity of search algorithms that are robust to imprecise queries and that adapt to perturbations of the order of the elements. We give (almost) tight results for various parameters that quantify query errors and that measure array disorder. In particular, we exhibit settings where query complexities of $\log n + ck$, $(1+\varepsilon)\log n + ck$, and $\sqrt{cnk}+o(nk)$ are best-possible for parameter value $k$, any $\varepsilon>0$, and constant $c$.


title: Qualitative robustness for bootstrap approximations

abstract: An important property of statistical estimators is qualitative robustness, that is small changes in the distribution of the data only result in small chances of the distribution of the estimator. Moreover, in practice, the distribution of the data is commonly unknown, therefore bootstrap approximations can be used to approximate the distribution of the estimator. Hence qualitative robustness of the statistical estimator under the bootstrap approximation is a desirable property. Currently most theoretical investigations on qualitative robustness assume independent and identically distributed pairs of random variables. However, in practice this assumption is not fulfilled. Therefore, we examine the qualitative robustness of bootstrap approximations for non-i.i.d. random variables, for example $\alpha$-mixing and weakly dependent processes. In the i.i.d. case qualitative robustness is ensured via the continuity of the statistical operator, representing the estimator, see Hampel (1971) and Cuevas and Romo (1993). We show, that qualitative robustness of the bootstrap approximation is still ensured under the assumption that the statistical operator is continuous and under an additional assumption on the stochastic process. In particular, we require a convergence condition of the empirical measure of the underlying process, the so called Varadarajan property.


title: On the efficient computation of a generalized Jacobian of the projector over the Birkhoff polytope

abstract: We derive an explicit formula, as well as an efficient procedure, for constructing a generalized Jacobian for the projector of a given square matrix onto the Birkhoff polytope, i.e., the set of doubly stochastic matrices. To guarantee the high efficiency of our procedure, a semismooth Newton method for solving the dual of the projection problem is proposed and efficiently implemented. Extensive numerical experiments are presented to demonstrate the merits and effectiveness of our method by comparing its performance against other powerful solvers such as the commercial software Gurobi and the academic code PPROJ [{\sc Hager and Zhang}, SIAM Journal on Optimization, 26 (2016), pp.~1773--1798]. In particular, our algorithm is able to solve the projection problem with over one billion variables and nonnegative constraints to a very high accuracy in less than 15 minutes on a modest desktop computer. More importantly, based on our efficient computation of the projections and their generalized Jacobians, we can design a highly efficient augmented Lagrangian method (ALM) for solving a class of convex quadratic programming (QP) problems constrained by the Birkhoff polytope. The resulted ALM is demonstrated to be much more efficient than Gurobi in solving a collection of QP problems arising from the relaxation of quadratic assignment problems.


title: Spin-resolved dynamical conductance of correlated large-spin magnetic molecules

abstract: The finite-frequency transport properties of a large-spin molecule attached to ferromagnetic contacts are studied theoretically in the Kondo regime. The focus is on the behavior of the dynamical conductance in the linear response regime, which is determined by using the numerical renormalization group method. It is shown that the dynamical conductance depends greatly on the magnetic configuration of the device and intrinsic parameters of the molecule. In particular, conductance exhibits characteristic features for frequencies corresponding to the dipolar and quadrupolar exchange fields resulting from the presence of spin-dependent tunneling. Moreover, a dynamical spin accumulation in the molecule, associated with the off-diagonal-in-spin component of the conductance, is predicted. This spin accumulation becomes enhanced with increasing the spin polarization of the leads, and it results in a nonmonotonic dependence of the conductance on frequency, with local maxima occurring for characteristic energy scales.


title: Equation of State of Neutron Stars with Junction Conditions in the Starobinsky Model

abstract: We study the Starobinsky or $R^2$ model of $f(R)=R+\alpha R^2$ for neutron stars with the structure equations represented by the coupled differential equations and the \emph{polytropic} type of the matter equation of state. The junction conditions of $f(R)$ gravity are used as the boundary conditions to match the Schwarschild solution at the surface of the star. Based on these the conditions, we demonstrate that the coupled differential equations can be solved \emph{directly}. In particular, from the dimensionless equation of state $\bar{\rho} = \bar{k}\, \bar{p}^{\,\gamma}$ with $\bar{k}\sim5.0$ and $\gamma\sim0.75$ and the constraint of $\alpha\lesssim {1.47722}\times 10^{7}\, \text{m}^2$, we obtain the \emph{minimal} mass of the NS to be around 1.44 $M_{\odot}$. In addition, if $\bar{k}$ is larger than 5.0, the mass and radius of the NS would be smaller.


title: Isometries in spaces of K\"ahler potentials

abstract: The space of K\"ahler potentials in a compact K\"ahler manifold, endowed with Mabuchi's metric, is an infinite dimensional Riemannian manifold. We characterize local isometries between spaces of K\"ahler potentials, and prove existence and uniqueness for such isometries.


title: Fission dynamics of $^{252}$Cf

abstract: The time-dependent generator coordinate method with the gaussian overlap approximation (TDGCM+GOA) formalism is applied to describe the fission of $^{252}$Cf. We perform analysis of fission from the initial states laying in the energetic range from the ground state to the state located 4 MeV above the fission barrier. The fission fragment mass distributions, obtained for different parity, energy of levels and types of mixed states, are calculated and compared with experimental data. The impact of the total time of wave packet propagation on the final results is studied as well. The weak dependence of obtained mass yields on the initial conditions is shown.


title: An Efficient Method for online Detection of Polychronous Patterns in Spiking Neural Network

abstract: Polychronous neural groups are effective structures for the recognition of precise spike-timing patterns but the detection method is an inefficient multi-stage brute force process that works off-line on pre-recorded simulation data. This work presents a new model of polychronous patterns that can capture precise sequences of spikes directly in the neural simulation. In this scheme, each neuron is assigned a randomized code that is used to tag the post-synaptic neurons whenever a spike is transmitted. This creates a polychronous code that preserves the order of pre-synaptic activity and can be registered in a hash table when the post-synaptic neuron spikes. A polychronous code is a sub-component of a polychronous group that will occur, along with others, when the group is active. We demonstrate the representational and pattern recognition ability of polychronous codes on a direction selective visual task involving moving bars that is typical of a computation performed by simple cells in the cortex. The computational efficiency of the proposed algorithm far exceeds existing polychronous group detection methods and is well suited for online detection.


title: Global stability and $H$-theorem in lattice models with non-conservative interactions

abstract: In kinetic theory, a system is usually described by its one-particle distribution function $f(\mathbf{r},\mathbf{v},t)$, such that $f(\mathbf{r},\mathbf{v},t)d\mathbf{r} d\mathbf{v}$ is the fraction of particles with positions and velocities in the intervals $(\mathbf{r}, \mathbf{r}+d\mathbf{r})$ and $(\mathbf{v}, \mathbf{v}+d\mathbf{v})$, respectively. Therein, global stability and the possible existence of an associated Lyapunov function or $H$-theorem are open problems when non-conservative interactions are present, as in granular fluids. Here, we address this issue in the framework of a lattice model for granular-like velocity fields. For a quite general driving mechanism, including both boundary and bulk driving, we show that the steady state reached by the system in the long time limit is globally stable. This is done by proving analytically that a certain $H$-functional is non-increasing in the long time limit. Moreover, for two specific energy injection mechanisms, we are able to demonstrate that the proposed $H$-functional is non-increasing for all times. Also, we put forward a proof that clearly illustrates why the "classical" Boltzmann functional $H_{B}[f]=\int\! d\mathbf{r} \, d\mathbf{v} f(\mathbf{r},\mathbf{v},t) \ln f(\mathbf{r},\mathbf{v},t)$ is inadequate for systems with non-conservative interactions. Not only is this done for the simplified kinetic description that holds in the lattice models analysed here but also for a general kinetic equation, like Boltzmann's or Enskog's.

categories: cond-mat.stat-mech cond-mat.soft math-ph math.MP


title: SurvivalNet: Predicting patient survival from diffusion weighted magnetic resonance images using cascaded fully convolutional and 3D convolutional neural networks

abstract: Automatic non-invasive assessment of hepatocellular carcinoma (HCC) malignancy has the potential to substantially enhance tumor treatment strategies for HCC patients. In this work we present a novel framework to automatically characterize the malignancy of HCC lesions from DWI images. We predict HCC malignancy in two steps: As a first step we automatically segment HCC tumor lesions using cascaded fully convolutional neural networks (CFCN). A 3D neural network (SurvivalNet) then predicts the HCC lesions' malignancy from the HCC tumor segmentation. We formulate this task as a classification problem with classes being "low risk" and "high risk" represented by longer or shorter survival times than the median survival. We evaluated our method on DWI of 31 HCC patients. Our proposed framework achieves an end-to-end accuracy of 65% with a Dice score for the automatic lesion segmentation of 69% and an accuracy of 68% for tumor malignancy classification based on expert annotations. We compared the SurvivalNet to classical handcrafted features such as Histogram and Haralick and show experimentally that SurvivalNet outperforms the handcrafted features in HCC malignancy classification. End-to-end assessment of tumor malignancy based on our proposed fully automatic framework corresponds to assessment based on expert annotations with high significance (p>0.95).


title: Terahertz detection using mechanical resonators based on 2D materials

abstract: We have investigated a THz detection scheme based on mixing of electrical signals in a voltage-dependent capacitance made out of suspended graphene. We have analyzed both coherent and incoherent detection regimes and compared their performance with the state of the art. Using a high-amplitude local oscillator, we anticipate potential for quantum limited detection in the coherent mode. The sensitivity stems from the extraordinary mechanical and electrical properties of atomically thin graphene or graphene-related 2D materials.


title: Moments of volumes of lower-dimensional random simplices are not monotone

abstract: In a $d$-dimensional convex body $K$, for $n \leq d+1$, random points $X_0, \dots, X_{n-1}$ are chosen according to the uniform distribution in $K$. Their convex hull is a random $(n-1)$-simplex with probability $1$. We denote its $(n-1)$-dimensional volume by $V_{K[n]}$. The $k$-th moment of the $(n-1)$-dimensional volume of a random $(n-1)$-simplex is monotone under set inclusion, if $K \subseteq L$ implies that the $k$-th moment of $V_{K[n]}$ is not larger than that of $V_{L[n]}$. Extending work of Rademacher [On the monotonicity of the expected volume of a random simplex. Mathematika 58 (2012), 77--91] and Reichenwallner and Reitzner [On the monotonicity of the moments of volumes of random simplices. Mathematika 62 (2016), 949--958], it is shown that for $n \leq d$, the moments of $V_{K[n]}$ are not monotone under set inclusion.


title: Relation between regional uncertainty spillovers in the global banking system

abstract: We report on time-varying network connectedness within three banking systems: North America, the EU, and ASEAN. The original method by Diebold and Yilmaz is improved by using exponentially weighted daily returns and ridge regularization on vector autoregression (VAR) and forecast error variance decomposition (FEVD). We compute the total network connectedness for each of the three banking systems, which quantifies regional uncertainty. Results over rolling windows of 300 days during the period between 2005 and 2015 reveal changing uncertainty patterns which are similar across regions, with common peaks associated with identifiable exogenous events. Lead-lag relationships among changes of total network connectedness of the three systems, quantified by transfer entropy, reveal that uncertainties in the three regional systems are significantly causally related, with the North American system having the largest influence on EU and ASEAN.


title: On the Comparison of Context-Free Grammars

abstract: In this paper we consider the problem of context-free grammars comparison from the analysis point of view. We show that the problem can be reduced to numerical solution of systems of nonlinear matrix equations. The approach presented here forms a basis for probabilistic comparison algorithms oriented to automatic assessment of of student's answers in computer science.


title: Cartesian products of directed graphs with loops

abstract: We show that every nontrivial finite or infinite connected directed graph with loops and at least one vertex without a loop is uniquely representable as a Cartesian or weak Cartesian product of prime graphs. For finite graphs the factorization can be computed in linear time and space.


title: Cutting Planes for Families Implying Frankl's Conjecture

abstract: We find previously unknown families of sets which ensure Frankl's conjecture holds for all families that contain them using an algorithmic framework. The conjecture states that for any nonempty union-closed (UC) family there exists an element of the ground set in at least half the sets of the considered UC family. Poonen's Theorem characterizes the existence of weights which determine whether a given UC family implies the conjecture for all UC families which contain it. We design a cutting-plane method that computes the explicit weights which imply the existence conditions of Poonen's Theorem. This method enables us to answer several open questions regarding structural properties of UC families, including the construction of a counterexample to a conjecture of Morris from 2006.


title: Periodic and q-potent elements in rings

abstract: This paper, we consider some properties of rings via q-potent and periodic elements. In this paper we give some results of rings in which every element is a sum of an idempotent and a q-potent that commute; periodic rings and k-potent elements of algebras.


title: Kinetic formulation of a 2 $\times$ 2 hyperbolic system arising in gas chromatography

abstract: The PSA system commonly used in the context of gas-solid chromatography is reformulated as a single kinetic equation using an additional kinetic variable. A kinetic numerical scheme is built from this new formulation and its behavior is tested on solving the Riemann problem in different configurations leading to single or composite waves.


title: Quantum critical phase with infinite projected entangled paired states

abstract: A classification of SU(2)-invariant Projected Entangled Paired States (PEPS) on the square lattice, based on a unique site tensor, has been recently introduced by Mambrini et al.~\cite{Mambrini2016}. It is not clear whether such SU(2)-invariant PEPS can either i) exhibit long-range magnetic order (like in the N\'eel phase) or ii) describe a genuine quantum critical point (QCP) or quantum critical phase (QCPh) separating two ordered phases. Here, we identify a specific family of SU(2)-invariant PEPS of the classification which provides excellent variational energies for the $J_1-J_2$ frustrated Heisenberg model, especially at $J_2=0.5$, corresponding to the approximate location of the QCP or QCPh separating the N\'eel phase from a dimerized phase. The PEPS are build from virtual states belonging to the $\frac{1}{2}^{\otimes N} \oplus 0$ SU(2)-representation, i.e. with $N$ "colors" of virtual \hbox{spin-$\frac{1}{2}$}. Using a full update infinite-PEPS approach directly in the thermodynamic limit, based on the Corner Transfer Matrix renormalization algorithm supplemented by a Conjugate Gradient optimization scheme, we provide evidence of i) the absence of magnetic order and of ii) diverging correlation lengths (i.e. showing no sign of saturation with increasing environment dimension) in both the singlet and triplet channels, when the number of colors $N\ge 3$. We argue that such a PEPS gives a qualitative description of the QCP or QCPh of the $J_1-J_2$ model.


title: Refined Vertex Sparsifiers of Planar Graphs

abstract: We study the following version of cut sparsification. Given a large edge-weighted network $G$ with $k$ terminal vertices, compress it into a smaller network $H$ with the same terminals, such that every minimum terminal cut in $H$ approximates the corresponding one in $G$, up to a factor $q\geq 1$ that is called the quality. (The case $q=1$ is known also as a mimicking network). We provide new insights about the structure of minimum terminal cuts, leading to new results for cut sparsifiers of planar graphs. Our first contribution identifies a subset of the minimum terminal cuts, which we call elementary, that generates all the others. Consequently, $H$ is a cut sparsifier if and only if it preserves all the elementary terminal cuts (up to this factor $q$). This structural characterization lead to improved bounds on the size of $H$. For example, it improve the bound of mimicking-network size for planar graphs into a near-optimal one. Our second and main contribution is to refine the known bounds in terms of $\gamma=\gamma(G)$, which is defined as the minimum number of faces that are incident to all the terminals in a planar graph $G$. We prove that the number of elementary terminal cuts is $O((2k/\gamma)^{2\gamma})$ (compared to $O(2^k)$ terminal cuts), and furthermore obtain a mimicking-network of size $O(\gamma 2^{2\gamma} k^4)$, which is near-optimal as a function of $\gamma$. In the analysis we break the elementary terminal cuts into fragments, and count them carefully. Our third contribution is a duality between cut sparsification and distance sparsification for certain planar graphs, when the sparsifier $H$ is required to be a minor of $G$. This duality connects problems that were previously studied separately, implying new results, new proofs of known results, and equivalences between open gaps.


title: Interplay between social influence and competitive strategical games in multiplex networks

abstract: We present a model that takes into account the coupling between evolutionary game dynamics and social influence. Importantly, social influence and game dynamics take place in different domains, which we model as different layers of a multiplex network. We show that the coupling between these dynamical processes can lead to cooperation in scenarios where the pure game dynamics predicts defection. In addition, we show that the structure of the network layers and the relation between them can further increase cooperation. Remarkably, if the layers are related in a certain way, the system can reach a polarized metastable state.These findings could explain the prevalence of polarization observed in many social dilemmas.


title: Giant exciton Fano resonance in quasi-one-dimensional Ta$_2$NiSe$_5$

abstract: We report the complex dielectric function of the quasi-one-dimensional chalcogenide Ta$_2$NiSe$_5$, which exhibits a structural phase transition that has been attributed to exciton condensation below $T_c = 326$ K, and of the isostructural Ta$_2$NiS$_5$ which does not exhibit such a transition. Using spectroscopic ellipsometry, we have detected exciton doublets with pronounced Fano lineshapes in both the compounds. The exciton Fano resonances in Ta$_2$NiSe$_5$ display an order of magnitude higher intensity than those in Ta$_2$NiS$_5$. In conjunction with prior theoretical work by E. Rashba, we attribute this observation to the giant oscillator strength of spatially extended exciton-phonon bound states in Ta$_2$NiSe$_5$. The formation of exciton-phonon complexes in Ta$_2$NiS$_5$ and Ta$_2$NiSe$_5$ is confirmed by the pronounced temperature dependence of sharp interband transitions in the optical spectra, whose peak energies and widths scale with the thermal population of optical phonon modes. The description of the optically excited states in terms of strongly overlapping exciton complexes is in good agreement with the hypothesis of an EI ground state.


title: Simulating quantum light propagation through atomic ensembles using matrix product states

abstract: A powerful method to interface quantum light with matter is to propagate the light through an ensemble of atoms. Recently, a number of such interfaces have emerged, most prominently Rydberg ensembles, that enable strong nonlinear interactions between propagating photons. A largely open problem is whether these systems produce exotic many-body states of light and developing new tools to study propagation in the large photon number limit is highly desirable. Here, we provide a method based on a "spin model" that maps quasi one-dimensional (1D) light propagation to the dynamics of an open 1D interacting spin system, where all photon correlations are obtained from those of the spins. The spin dynamics in turn are numerically solved using the toolbox of matrix product states. We apply this formalism to investigate vacuum induced transparency, wherein the different photon number components of a pulse propagate with number-dependent group velocity and separate at output.


title: Northwest Africa 5958: a weakly altered CM-related ungrouped chondrite, not a CI3

abstract: Northwest Africa (NWA) 5958 is a carbonaceous chondrite found in Morocco in 2009. Preliminary chemical and isotopic data leading to its initial classification as C3.0 ungrouped have prompted us to conduct a multi-technique study of this meteorite and present a general description here. The petrography and chemistry of NWA 5958 is most similar to a CM chondrite, with a low degree of aqueous alteration, apparently under oxidizing conditions, and evidence of a second, limited alteration episode manifested by alteration fronts. The oxygen isotopic composition, with $\Delta^{17}$O = -4.3 $\permil$, is more 16O-rich than all CM chondrites, indicating, along with other compositional arguments, a separate parent body of origin. We suggest that NWA 5958 be reclassified as an ungrouped carbonaceous chondrite related to the CM group.


title: Hydrogen in diffuse molecular clouds in the Milky Way: Atomic column densities and molecular fraction along prominent lines of sight

abstract: Recent submillimeter and far-infrared wavelength observations of absorption in the rotational ground-state lines of various simple molecules against distant Galactic continuum sources have opened the possibility of studying the chemistry of diffuse molecular clouds throughout the Milky Way. In order to calculate abundances, the column densities of molecular and atomic hydrogen, HI, must be known. We aim at determining the atomic hydrogen column densities for diffuse clouds located on the sight lines toward a sample of prominent high-mass star-forming regions that were intensely studied with the HIFI instrument onboard Herschel. Based on Jansky Very Large Array data, we employ the 21 cm HI absorption-line technique to construct profiles of the HI opacity versus radial velocity toward our target sources. These profiles are combined with lower resolution archival data of extended HI emission to calculate the HI column densities of the individual clouds along the sight lines. We employ Bayesian inference to estimate the uncertainties of the derived quantities. Our study delivers reliable estimates of the atomic hydrogen column density for a large number of diffuse molecular clouds at various Galactocentric distances. Together with column densities of molecular hydrogen derived from its surrogates observed with HIFI, the measurements can be used to characterize the clouds and investigate the dependence of their chemistry on the molecular fraction, for example.


title: I Ate This: A Photo-based Food Journaling System with Expert Feedback

abstract: What we eat is one of the most frequent and important health decisions we make in daily life, yet it remains notoriously difficult to capture and understand. Effective food journaling is thus a grand challenge in personal health informatics. In this paper we describe a system for food journaling called I Ate This, which is inspired by the Remote Food Photography Method (RFPM). I Ate This is simple: you use a smartphone app to take a photo and give a very basic description of any food or beverage you are about to consume. Later, a qualified dietitian will evaluate your photo, giving you feedback on how you did and where you can improve. The aim of I Ate This is to provide a convenient, visual and reliable way to help users learn from their eating habits and nudge them towards better choices each and every day. Ultimately, this incremental approach can lead to long-term behaviour change. Our goal is to bring RFPM to a wider audience, through APIs that can be incorporated into other apps.


title: Reflection Separation Using Guided Annotation

abstract: Photographs taken through a glass surface often contain an approximately linear superposition of reflected and transmitted layers. Decomposing an image into these layers is generally an ill-posed task and the use of an additional image prior and user provided cues is presently necessary in order to obtain good results. Current annotation approaches rely on a strong sparsity assumption. For images with significant texture this assumption does not typically hold, thus rendering the annotation process unviable. In this paper we show that using a Gaussian Mixture Model patch prior, the correct local decomposition can almost always be found as one of 100 likely modes of the posterior. Thus, the user need only choose one of these modes in a sparse set of patches and the decomposition may then be completed automatically. We demonstrate the performance of our method using synthesized and real reflection images.


title: Optimal control for perfect state transfer in linear quantum memory

abstract: A quantum memory is a system that enables transfer, storage, and retrieval of optical quantum states by ON/OFF switching of the control signal in each stages of the memory. In particular, it is known that, for perfect transfer of a single-photon state, appropriate shaping of the input pulse is required. However, in general, such a desirable pulse shape has a complicated form, which would be hard to generate in practice. In this paper, for a wide class of linear quantum memory systems, we develop a method that reduces the complexity of the input pulse shape of a single-photon while maintaining the perfect state transfer. The key idea is twofold; (i) the control signal is allowed to vary continuously in time to introduce an additional degree of freedom, and then (ii) an optimal control problem is formulated to design a simple-formed input pulse and the corresponding control signal. Numerical simulations are conducted for Lambda-type atomic media and a networked atomic ensembles, to show the effectiveness of the proposed method.


title: A Statistical Learning Approach to Modal Regression

abstract: This paper studies the nonparametric modal regression problem systematically from a statistical learning view. Originally motivated by pursuing a theoretical understanding of the maximum correntropy criterion based regression (MCCR), our study reveals that MCCR with a tending-to-zero scale parameter is essentially modal regression. We show that nonparametric modal regression problem can be approached via the classical empirical risk minimization. Some efforts are then made to develop a framework for analyzing and implementing modal regression. For instance, the modal regression function is described, the modal regression risk is defined explicitly and its \textit{Bayes} rule is characterized; for the sake of computational tractability, the surrogate modal regression risk, which is termed as the generalization risk in our study, is introduced. On the theoretical side, the excess modal regression risk, the excess generalization risk, the function estimation error, and the relations among the above three quantities are studied rigorously. It turns out that under mild conditions, function estimation consistency and convergence may be pursued in modal regression as in vanilla regression protocols, such as mean regression, median regression, and quantile regression. However, it outperforms these regression models in terms of robustness as shown in our study from a re-descending M-estimation view. This coincides with and in return explains the merits of MCCR on robustness. On the practical side, the implementation issues of modal regression including the computational algorithm and the tuning parameters selection are discussed. Numerical assessments on modal regression are also conducted to verify our findings empirically.


title: Gevrey estimates for one dimensional parabolic invariant manifolds of non-hyperbolic fixed points

abstract: We study the Gevrey character of a natural parameterization of one dimensional invariant manifolds associated to a parabolic direction of fixed points of analytic maps, that is, a direction associated with an eigenvalue equal to $1$. We show that, under general hypotheses, these invariant manifolds are Gevrey with type related to some explicit constants. We provide examples of the optimality of our results as well as some applications to celestial mechanics, namely, the Sitnikov problem and the restricted planar three body problem.


title: Latent Variable Dialogue Models and their Diversity

abstract: We present a dialogue generation model that directly captures the variability in possible responses to a given input, which reduces the `boring output' issue of deterministic dialogue models. Experiments show that our model generates more diverse outputs than baseline models, and also generates more consistently acceptable output than sampling from a deterministic encoder-decoder model.


title: Markov $L_2$ inequality with the Gegenbauer weight

abstract: For the Gegenbauer weight function $w_{\lambda}(t)=(1-t^2)^{\lambda-1/2}$, $\lambda>-1/2$, we denote by $\Vert\cdot\Vert_{w_{\lambda}}$ the associated $L_2$-norm, $$ \Vert f\Vert_{w_{\lambda}}:=\Big(\int_{-1}^{1}w_{\lambda}(t)f^2(t)\,dt\Big)^{1/2}. $$ We study the Markov inequality $$ \Vert p^{\prime}\Vert_{w_{\lambda}}\leq c_{n}(\lambda)\,\Vert p\Vert_{w_{\lambda}},\qquad p\in \mathcal{P}_n, $$ where $\mathcal{P}_n$ is the class of algebraic polynomials of degree not exceeding $n$. Upper and lower bounds for the best Markov constant $c_{n}(\lambda)$ are obtained, which are valid for all $n\in \mathbb{N}$ and $\lambda>-\frac{1}{2}$.


title: Virtual quandle for links in lens spaces

abstract: We construct a virtual quandle for links in lens spaces $L(p,q)$, with $q=1$. This invariant has two valuable advantages over an ordinary fundamental quandle for links in lens spaces: the virtual quandle is an essential invariant and the presentation of the virtual quandle can be easily written from the band diagram of a link.


title: Evolution of magnetism in cerium doped LaCo2P2 crystals: a magnetic phase diagram

abstract: ThCr2Si2-type phosphide ACo2P2 (A=Rare earth elements) has the same structure as iron arsenides, but their magnetic behaviors are quite distinct. In this paper, we for the first time grew a series of La1-xCexCo2P2 single crystals (x=0.0 to1.0), and made structural and magnetic characterizations. This allows us to carry out a careful investigation on the evolution of magnetism with cerium content and build a magnetic phase diagram. We found that the introduction of cerium induces a rapid decrease of c-axis and a change from ferromagnetic (FM) to antiferromagnetic (AFM) states. By employing first-principles band-structure calculations, we identify the formation of P-P bonding with the shortening of c-axis, which effectively drives an increase of AFM interaction and eventually leads to AFM ordering in the high doping region. The present study may shed light on the interplay between the structural collapsing and electronic/magnetic properties in 122 iron pnictides.


title: Does warm debris dust stem from asteroid belts?

abstract: Many debris discs reveal a two-component structure, with a cold outer and a warm inner component. While the former are likely massive analogues of the Kuiper belt, the origin of the latter is still a matter of debate. In this work we investigate whether the warm dust may be a signature of asteroid belt analogues. In the scenario tested here the current two-belt architecture stems from an originally extended protoplanetary disc, in which planets have opened a gap separating it into the outer and inner discs which, after the gas dispersal, experience a steady-state collisional decay. This idea is explored with an analytic collisional evolution model for a sample of 225 debris discs from a Spitzer/IRS catalogue that are likely to possess a two-component structure. We find that the vast majority of systems (220 out of 225, or 98%) are compatible with this scenario. For their progenitors, original protoplanetary discs, we find an average surface density slope of $-0.93\pm0.06$ and an average initial mass of $\left(3.3^{+0.4}_{-0.3}\right)\times 10^{-3}$ solar masses, both of which are in agreement with the values inferred from submillimetre surveys. However, dust production by short-period comets and - more rarely - inward transport from the outer belts may be viable, and not mutually excluding, alternatives to the asteroid belt scenario. The remaining five discs (2% of the sample: HIP 11486, HIP 23497, HIP 57971, HIP 85790, HIP 89770) harbour inner components that appear inconsistent with dust production in an "asteroid belt." Warm dust in these systems must either be replenished from cometary sources or represent an aftermath of a recent rare event, such as a major collision or planetary system instability.


title: Practical, Linear-time, Fully Distributed Algorithms for Irregular Gather and Scatter

abstract: We present new, simple, fully distributed, practical algorithms with linear time communication cost for irregular gather and scatter operations in which processors contribute or consume possibly different amounts of data. In a linear cost transmission model with start-up latency $\alpha$ and cost per unit $\beta$, the new algorithms take time $3|{\log_2 p}|\alpha+\beta \sum_{i\neq r}m_i$ where $p$ is the number of processors, $m_i$ the amount of data for processor $i, 0\leq i<p$, and processor $r, 0\leq r<p$ a root processor determined by the algorithm. For a fixed, externally given root processor $r$, there is an additive penalty of at most $\beta(M_{d'}-m_{r_{d'}}-\sum_{0\leq j<d'}M_j)$ time steps where each $M_j$ is the total amount of data in a tree of $2^j$ different processors with roots $r_j$ as constructed by the algorithm. The worst-case penalty is less than $\beta \sum_{i\neq r}m_i$ time steps. The algorithms have attractive properties for implementing the operations for MPI (the Message-Passing Interface). Standard algorithms using fixed trees take time either $|{\log_2 p}|(\alpha+\beta \sum_{i\neq r} m_i)$ in the worst case, or $\sum_{i\neq r}(\alpha+\beta m_i)$. We have used the new algorithms to give prototype implementations for the MPI_Gatherv and MPI_Scatterv collectives of MPI, and present benchmark results from a small and a medium-large InfiniBand cluster. In order to structure the experimental evaluation we formulate new performance guidelines for irregular collectives that can be used to assess the performance in relation to the corresponding regular collectives. We show that the new algorithms can fulfill these performance expectations with a large margin, and that standard implementations do not.


title: Effects of magnetic fields on quark-antiquark interactions

abstract: We discuss some recent results obtained in the study of strong quark-antiquark interactions in the presence of intense external magnetic fields by means of lattice QCD simulations. We confirm previous findings and show that both at zero and finite temperature the external field induces anisotropies in the static quark potential. An in-depth study suggests that the effects are essentially due to the variation of the string tension whose angular dependence can be nicely parametrized by the first allowed term in a Fourier expansion. In the confined phase at high temperature, we observe that the suppression of the string tension is enhanced as the strength of the external field increases. Our results support the idea that the loss of confining properties is the dominant effect related to the decrease of $T_c$ as a function of $B$.


title: Interaction of a Laguerre-Gaussian beam with Rydberg atoms

abstract: Transfer mechanism of orbital angular moment(OAM) of light to trapped ground-state atoms under paraxial approximation is well known. Here we show how optical OAM of a Laguerre-Gaussian(LG) beam under paraxial approximation can be transferred to trapped Rydberg atoms. Optical OAM is shown to be transferable to a Rydberg electronic state in dipole transition. The Gaussian part of the profile of the LG beam, which is generally neglected , is found to have an important effect on the OAM transfer to the Rydberg atoms. Numerical calculations are calculated based on this theory for Rubidium Rydberg atoms trapped in a harmonic potential. Our results exhibit the mixing of final states of different parities.


title: Automatic Liver and Tumor Segmentation of CT and MRI Volumes using Cascaded Fully Convolutional Neural Networks

abstract: Automatic segmentation of the liver and hepatic lesions is an important step towards deriving quantitative biomarkers for accurate clinical diagnosis and computer-aided decision support systems. This paper presents a method to automatically segment liver and lesions in CT and MRI abdomen images using cascaded fully convolutional neural networks (CFCNs) enabling the segmentation of a large-scale medical trial or quantitative image analysis. We train and cascade two FCNs for a combined segmentation of the liver and its lesions. In the first step, we train a FCN to segment the liver as ROI input for a second FCN. The second FCN solely segments lesions within the predicted liver ROIs of step 1. CFCN models were trained on an abdominal CT dataset comprising 100 hepatic tumor volumes. Validations on further datasets show that CFCN-based semantic liver and lesion segmentation achieves Dice scores over 94% for liver with computation times below 100s per volume. We further experimentally demonstrate the robustness of the proposed method on an 38 MRI liver tumor volumes and the public 3DIRCAD dataset.


title: Regularization by noise in one-dimensional continuity equation

abstract: A linear stochastic continuity equation with non-regular coefficients is considered. We prove existence and uniqueness of strong solution, in the probabilistic sense, to the Cauchy problem when the vector field has low regularity, in which the classical DiPerna-Lions-Ambrosio theory of uniqueness of distributional solutions does not apply. We solve partially the open problem that is the case when the vector-field has random dependence. In addition, we prove a stability result for the solutions.


title: Generalising rate heterogeneity across sites in statistical phylogenetics

abstract: Phylogenetics uses alignments of molecular sequence data to learn about evolutionary trees relating species. Along branches, sequence evolution is modelled using a continuous-time Markov process characterised by an instantaneous rate matrix. Early models assumed the same rate matrix governed substitutions at all sites of the alignment, ignoring variation in evolutionary pressures. Substantial improvements in phylogenetic inference and model fit were achieved by augmenting these models with multiplicative random effects that describe the result of variation in selective constraints and allow sites to evolve at different rates which linearly scale a baseline rate matrix. Motivated by this pioneering work, we consider an extension using a quadratic, rather than linear, transformation. The resulting models allow for variation in the selective coefficients of different types of point mutation at a site in addition to variation in selective constraints. We derive properties of the extended models. For certain non-stationary processes, the extension gives a model that allows variation in sequence composition both across sites and taxa. We adopt a Bayesian approach, describe an MCMC algorithm for posterior inference and provide software. Our quadratic models are applied to alignments spanning the tree of life and compared with site-homogeneous and linear models.


title: Asymptotic Freedom in the BV Formalism

abstract: We define the beta-function of a perturbative quantum field theory in the mathematical framework introduced by Costello -- combining perturbative renormalization and the BV formalism -- as the cohomology class of a certain element in the obstruction-deformation complex. We show that the one-loop beta-function is a well-defined element of the local deformation complex for translation-invariant and classically scale-invariant theories, and furthermore that it is locally constant as a function on the space of classical interactions and computable as a rescaling anomaly, or as the logarithmic one-loop counterterm. We compute the one-loop beta-function in first-order Yang--Mills theory, recovering the famous asymptotic freedom for Yang--Mills in a mathematical context.


title: Topological states at the (001) surface of SrTiO3

abstract: Defect-free SrTiO3 (STO) is a band insulator but Angle Resolved Photoemission Spectroscopy (ARPES) experiments have demonstrated the existence of a nanometer thin two-dimensional electron liquid (2DEG) at the (001) oriented surface of this compound. The bulk is a trivial insulator, but our theoretical study reveals that the parity of electronic wavefunctions in this 2DEG is inverted in the vicinity of special points in reciprocal space where the low-energy dispersion consists of four gapped Dirac cones with a tilted and anisotropic shape. This gives rise to linearly dispersing topological edge states at the one-dimensional boundary. We propose to probe these modes by measuring the Josephson radiation from gapless bound Andreev states in STO based junctions, as it is predicted that they display distinctive signatures of topology.


title: On square functions with independent increments and Sobolev spaces on the line

abstract: We prove a characterization of some $L^p$-Sobolev spaces involving the quadratic symmetrization of the Calder\'on commutator kernel, which is related to a square function with differences of difference quotients. An endpoint weak type estimate is established for functions in homogeneous Hardy-Sobolev spaces $\dot H^1_\alpha$. We also use a local version of this square function to characterize pointwise differentiability for functions in the Zygmund class.


title: Mean-field theory for confinement transitions and magnetization plateaux in spin ice

abstract: We study phase transitions in classical spin ice at nonzero magnetization, by introducing a mean-field theory designed to capture the interplay between confinement and topological constraints. The method is applied to a model of spin ice in an applied magnetic field along the [100] crystallographic direction and yields a phase diagram containing the Coulomb phase as well as a set of magnetization plateaux. We argue that the lobe structure of the phase diagram, strongly reminiscent of the Bose-Hubbard model, is generic to Coulomb spin liquids.


title: On the Spectral Efficiency and Fairness in Full-Duplex Cellular Networks

abstract: To increase the spectral efficiency of wireless networks without requiring full-duplex capability of user devices, a potential solution is the recently proposed three-node full-duplex mode. To realize this potential, networks employing three-node full-duplex transmissions must deal with self-interference and user-to-user interference, which can be managed by frequency channel and power allocation techniques. Whereas previous works investigated either spectral efficient or fair mechanisms, a scheme that balances these two metrics among users is investigated in this paper. This balancing scheme is based on a new solution method of the multi-objective optimization problem to maximize the weighted sum of the per-user spectral efficiency and the minimum spectral efficiency among users. The mixed integer non-linear nature of this problem is dealt by Lagrangian duality. Based on the proposed solution approach, a low-complexity centralized algorithm is developed, which relies on large scale fading measurements that can be advantageously implemented at the base station. Numerical results indicate that the proposed algorithm increases the spectral efficiency and fairness among users without the need of weighting the spectral efficiency. An important conclusion is that managing user-to-user interference by resource assignment and power control is crucial for ensuring spectral efficient and fair operation of full-duplex networks.


title: Berezin-toeplitz quantization and complex weyl quantization of the torus t${}^2$

abstract: In this paper, we give a correspondence between the Berezin-Toeplitz and the complex Weyl quantizations of the torus $ \mathbb{T}^2$. To achieve this, we use the correspondence between the Berezin-Toeplitz and the complex Weyl quantizations of the complex plane and a relation between the Berezin-Toeplitz quantization of a periodic symbol on the real phase space $\mathbb{R}^2$ and the Berezin-Toeplitz quantization of a symbol on the torus $ \mathbb{T}^2 $.

categories: math-ph math.AP math.FA math.MP math.OA math.SG


title: Permanency of the age-structured population model on several temporally variable patches

abstract: We consider a system of nonlinear partial differential equations that describes an age-structured population inhabiting several temporally varying patches. We prove existence and uniqueness of solution and analyze its large-time behavior in cases when the environment is constant and when it changes periodically. A pivotal assumption is that individuals can disperse and that each patch can be reached from every other patch, directly or through several intermediary patches. We introduce the net reproductive operator and characteristic equations for time-independent and periodical models and prove that permanency is defined by the net reproductive rate for the whole system. If the net reproductive rate is less or equal to one, extinction on all patches is imminent. Otherwise, permanency on all patches is guaranteed. The proof is based on a new approach to analysis of large-time stability.


title: An experimental study of low-velocity impacts into granular material in reduced gravity

abstract: In order to improve our understanding of landing on small bodies and of asteroid evolution, we use our novel drop tower facility to perform low-velocity (2-40 cm s^-1), shallow impact experiments of a 10 cm diameter aluminum sphere into quartz sand in low effective gravities (~0.2-1 m s^-2). Using in situ accelerometers, we measure the acceleration profile during the impacts and determine the peak accelerations, collision durations and maximum penetration depth. We find that the penetration depth scales linearly with the collision velocity but is independent of the effective gravity for the experimental range tested, and that the collision duration is independent of both the effective gravity and the collision velocity. No rebounds are observed in any of the experiments. Our low-gravity experimental results indicate that the transition from the quasi-static regime to the inertial regime occurs for impact energies two orders of magnitude smaller than in similar impact experiments under terrestrial gravity. The lower energy regime change may be due to the increased hydrodynamic drag of the surface material in our experiments, but may also support the notion that the quasi-static regime reduces as the effective gravity becomes lower.


title: Representation Variety of Surface Groups

abstract: We give an exact formula for the dimension of the variety of homomorphisms from $S_g$ to $\mathit{any}$ semisimple real algebraic group, where $S_g$ is a surface group of genus $g \geq 2$.


title: Wages of wins: could an amateur make money from match outcome predictions?

abstract: Evaluating the accuracies of models for match outcome predictions is nice and well but in the end the real proof is in the money to be made by betting. To evaluate the question whether the models developed by us could be used easily to make money via sports betting, we evaluate three cases: NCAAB post-season, NBA season, and NFL season, and find that it is possible yet not without its pitfalls. In particular, we illustrate that high accuracy does not automatically equal high pay-out, by looking at the type of match-ups that are predicted correctly by different models.


title: Generating Adversarial Malware Examples for Black-Box Attacks Based on GAN

abstract: Machine learning has been used to detect new malware in recent years, while malware authors have strong motivation to attack such algorithms. Malware authors usually have no access to the detailed structures and parameters of the machine learning models used by malware detection systems, and therefore they can only perform black-box attacks. This paper proposes a generative adversarial network (GAN) based algorithm named MalGAN to generate adversarial malware examples, which are able to bypass black-box machine learning based detection models. MalGAN uses a substitute detector to fit the black-box malware detection system. A generative network is trained to minimize the generated adversarial examples' malicious probabilities predicted by the substitute detector. The superiority of MalGAN over traditional gradient based adversarial example generation algorithms is that MalGAN is able to decrease the detection rate to nearly zero and make the retraining based defensive method against adversarial examples hard to work.


title: Tensor Network Models of Multiboundary Wormholes

abstract: We study the entanglement structure of states dual to multiboundary wormhole geometries using tensor network models. Perfect and random tensor networks tiling the hyperbolic plane have been shown to provide good models of the entanglement structure in holography. We extend this by quotienting the plane by discrete isometries to obtain models of the multiboundary states. We show that there are networks where the entanglement structure is purely bipartite, extending results obtained in the large temperature limit. We analyse the entanglement structure in a range of examples.


title: Fano's inequality for random variables

abstract: We extend Fano's inequality, which controls the average probability of events in terms of the average of some $f$--divergences, to work with arbitrary events (not necessarily forming a partition) and even with arbitrary $[0,1]$--valued random variables, possibly in continuously infinite number. We provide two applications of these extensions, in which the consideration of random variables is particularly handy: we offer new and elegant proofs for existing lower bounds, on Bayesian posterior concentration (minimax or distribution-dependent) rates and on the regret in non-stochastic sequential learning.


title: The Guiding Influence of Stanley Mandelstam, from S-Matrix Theory to String Theory

abstract: The guiding influence of some of Stanley Mandelstam's key contributions to the development of theoretical high energy physics is discussed, from the motivation for the study of the analytic properties of the scattering matrix through to dual resonance models and their evolution into string theory.


title: Bounds for multivariate residues and for the polynomials in the elimination theorem

abstract: We present several upper bounds for the height of global residues of rational forms on an affine variety. As a consequence, we deduce upper bounds for the height of the coefficients in the Bergman-Weil trace formula. We also present upper bounds for the degree and the height of the polynomials in the elimination theorem on an affine variety. This is an arithmetic analogue of Jelonek's effective elimination theorem, that plays a crucial role in the proof of our bounds for the height of global residues.


title: Mass Gap without Confinement

abstract: We revisit a one-parameter family of three-dimensional gauge theories with known supergravity duals. We show that three infrared behaviors are possible. For generic values of the parameter, the theories exhibit a mass gap but no confinement, meaning no linear quark-antiquark potential; for one limiting value of the parameter the theory flows to an infrared fixed point; and for another limiting value it exhibits both a mass gap and confinement. Theories close to these limiting values exhibit quasi-conformal and quasi-confining dynamics, respectively. Eleven-dimensional supergravity provides a simple, geometric explanation of these features.


title: Rigidity of square-tiled interval exchange transformations

abstract: We look at interval exchange transformations defined as first return maps on the set of diagonals of a flow of direction $\theta$ on a square-tiled surface: using a combinatorial approach, we show that, when the surface has at least one true singularity both the flow and the interval exchange are rigid if and only if tan $\theta$ has bounded partial quotients. Moreover, if all vertices of the squares are singularities of the flat metric, and tan $\theta$ has bounded partial quotients, the square-tiled interval exchange transformation T is not of rank one. Finally, for another class of surfaces, those defined by the unfolding of billiards in Veech triangles, we build an uncountable set of rigid directional flows and an uncountable set of rigid interval exchange transformations.


title: Light enhancement by quasi-bound states in the continuum in dielectric arrays

abstract: The article reports on light enhancement by structural resonances in linear periodic arrays of identical dielectric elements. As the basic elements both spheres and rods with circular cross section have been considered. In either case, it has been demonstrated that high-$Q$ structural resonant modes originated from bound states in the continuum enable near-field amplitude enhancement by factor of $10$--$25$ in the red-to-near infrared range in lossy silicon. The asymptotic behavior of the $Q$-factor with the number of elements in the array is explained theoretically by analyzing quasi-bound states propagation bands.


title: Finite-size corrections in the random assignment problem

abstract: We analytically derive, in the context of the replica formalism, the first finite size corrections to the average optimal cost in the random assignment problem for a quite generic distribution law for the costs. We show that, when moving from a power-law distribution to a $\Gamma$ distribution, the leading correction changes both in sign and in its scaling properties. We also examine the behavior of the corrections when approaching a $\delta$-function distribution. By using a numerical solution of the saddle-point equations, we provide predictions that are confirmed by numerical simulations.


title: Approximate and exact solutions of intertwining equations through random spanning forests

abstract: For different reversible Markov kernels on finite state spaces, we look for families of probability measures for which the time evolution almost remains in their convex hull. Motivated by signal processing problems and metastability studies we are interested in the case when the size of such families is smaller than the size of the state space, and we want such distributions to be with small overlap among them. To this aim we introduce a squeezing function to measure the common overlap of such families, and we use random forests to build random approximate solutions of the associated intertwining equations for which we can bound from above the expected values of both squeezing and total variation errors. We also explain how to modify some of these approximate solutions into exact solutions by using those eigenvalues of the associated Laplacian with the largest absolute values.


title: An Extended Framework for Marginalized Domain Adaptation

abstract: We propose an extended framework for marginalized domain adaptation, aimed at addressing unsupervised, supervised and semi-supervised scenarios. We argue that the denoising principle should be extended to explicitly promote domain-invariant features as well as help the classification task. Therefore we propose to jointly learn the data auto-encoders and the target classifiers. First, in order to make the denoised features domain-invariant, we propose a domain regularization that may be either a domain prediction loss or a maximum mean discrepancy between the source and target data. The noise marginalization in this case is reduced to solving the linear matrix system $AX=B$ which has a closed-form solution. Second, in order to help the classification, we include a class regularization term. Adding this component reduces the learning problem to solving a Sylvester linear matrix equation $AX+BX=C$, for which an efficient iterative procedure exists as well. We did an extensive study to assess how these regularization terms improve the baseline performance in the three domain adaptation scenarios and present experimental results on two image and one text benchmark datasets, conventionally used for validating domain adaptation methods. We report our findings and comparison with state-of-the-art methods.


title: Homoclinic tangencies and singular hyperbolicity for three-dimensional vector fields

abstract: We prove that any vector field on a three-dimensional compact manifold can be approximated in the C1-topology by one which is singular hyperbolic or by one which exhibits a homoclinic tangency associated to a regular hyperbolic periodic orbit. This answers a conjecture by Palis. During the proof we obtain several other results with independent interest: a compactification of the rescaled sectional Poincar\'e flow and a generalization of Ma\~n\'e-Pujals-Sambarino theorem for three-dimensional C2 vector fields with singularities.


title: On energy-critical half-wave maps into $\mathbb{S}^2$

abstract: We consider the energy-critical half-wave maps equation $$\partial_t \mathbf{u} + \mathbf{u} \wedge |\nabla| \mathbf{u} = 0$$ for $\mathbf{u} : [0,T) \times \mathbb{R} \to \mathbb{S}^2$. We give a complete classification of all traveling solitary waves with finite energy. The proof is based on a geometric characterization of these solutions as minimal surfaces with (not necessarily free) boundary on $\mathbb{S}^2$. In particular, we discover an explicit Lorentz boost symmetry, which is implemented by the conformal M\"obius group on the target $\mathbb{S}^2$ applied to half-harmonic maps from $\mathbb{R}$ to $\mathbb{S}^2$. Complementing our classification result, we carry out a detailed analysis of the linearized operator $L$ around half-harmonic maps $\mathbf{Q}$ with arbitrary degree $m \geq 1$. Here we explicitly determine the nullspace including the zero-energy resonances; in particular, we prove the nondegeneracy of $\mathbf{Q}$. Moreover, we give a full description of the spectrum of $L$ by finding all its $L^2$-eigenvalues and proving their simplicity. Furthermore, we prove a coercivity estimate for $L$ and we rule out embedded eigenvalues inside the essential spectrum. Our spectral analysis is based on a reformulation in terms of certain Jacobi operators (tridiagonal infinite matrices) obtained from a conformal transformation of the spectral problem posed on $\mathbb{R}$ to the unit circle $\mathbb{S}$. Finally, we construct a unitary map which can be seen as a gauge transform tailored for a future stability and blowup analysis close to half-harmonic maps. Our spectral results also have potential applications to the half-harmonic map heat flow, which is the parabolic counterpart of the half-wave maps equation.


title: Quantitative statistical stability and speed of convergence to equilibrium for partially hyperbolic skew products

abstract: We consider a general relation between fixed point stability of suitably perturbed transfer operators and convergence to equilibrium (a notion which is strictly related to decay of correlations). We apply this relation to deterministic perturbations of a class of (piecewise) partially hyperbolic skew products whose behavior on the preserved fibration is dominated by the expansion of the base map. In particular we apply the results to power law mixing toral extensions. It turns out that in this case, the dependence of the physical measure on small deterministic perturbations, in a suitable anisotropic metric is at least Holder continuous, with an exponent which is explicitly estimated depending on the arithmetical properties of the system. We show explicit examples of toral extensions having actually Holder stability and non differentiable dependence of the physical measure on perturbations.


title: Refinement-based Specification and Security Analysis of Separation Kernels

abstract: Assurance of information-flow security by formal methods is mandated in security certification of separation kernels. As an industrial standard for improving safety, ARINC 653 has been complied with by mainstream separation kernels. Due to the new trend of integrating safe and secure functionalities into one separation kernel, security analysis of ARINC 653 as well as a formal specification with security proofs are thus significant for the development and certification of ARINC 653 compliant Separation Kernels (ARINC SKs). This paper presents a specification development and security analysis method for ARINC SKs based on refinement. We propose a generic security model and a stepwise refinement framework. Two levels of functional specification are developed by the refinement. A major part of separation kernel requirements in ARINC 653 are modeled, such as kernel initialization, two-level scheduling, partition and process management, and inter-partition communication. The formal specification and its security proofs are carried out in the Isabelle/HOL theorem prover. We have reviewed the source code of one industrial and two open-source ARINC SK implementations, i.e. VxWorks 653, XtratuM, and POK, in accordance with the formal specification. During the verification and code review, six security flaws, which can cause information leakage, are found in the ARINC 653 standard and the implementations.


title: A New Route for the Determination of Protein Structure in Physiological Environment

abstract: Revealing the structure of complex biological macromolecules, such as proteins, is an essential step for understanding the chemical mechanisms that determine the diversity of their functions. Synchrotron based x-ray crystallography and cryo-electron microscopy have made major contributions in determining thousands of protein structures even from micro-sized crystals. They suffer from some limitations that have not been overcome, such as radiation damage, the natural inability to crystallize of a number of proteins and experimental conditions for structure determination that are incompatible with the physiological environment. Today the ultrashort and ultra-bright pulses of X-ray free-electron lasers (XFELs) have made attainable the dream to determine protein structure before radiation damage starts to destroy the samples. However, the signal-to-noise ratio remains a great challenge to obtain usable diffraction patterns from a single protein molecule. We describe here a new methodology that should overcome the signal and protein crystallization limits. Using a multidisciplinary approach, we propose to create a two dimensional protein array with defined orientation attached on a self-assembled-monolayer. We develop a literature-based, flexible toolbox capable of assembling different proteins on a functionalized surface while keeping them under physiological conditions during the experiment, using a water-confining graphene cover.


title: X-ray survey of the North-America and Pelican star-forming complex (NGC7000/IC5070)

abstract: We present the first extensive X-ray study of the North-America and Pelican star-forming region (NGC7000/IC5070), with the aim of finding and characterizing its young stellar population. X-ray data from Chandra (four pointings) and XMM-Newton (seven pointings) were reduced and source detection was performed on each image. We complement the X-ray data with optical and near-IR data from the IPHAS, UKIDSS, and 2MASS catalogs, and with other published optical and Spitzer IR data. More than 700 X-ray sources are detected, the majority of which have an optical/NIR counterpart. This allows to identify young stars in different stages of formation. Less than 30% of X-ray sources are identified with a previously known young star. We argue that most X-ray sources with an optical/NIR counterpart, except perhaps for a few tens at near-zero reddening, are likely candidate members of the star-forming region, on the basis of both their optical/NIR magnitudes and colors, and of X-ray properties like spectrum hardness or flux variations. They are characterized by a wide range of extinction, and sometimes near-IR excesses, both of which prevent derivation of accurate stellar parameters. The optical color-magnitude diagram suggests ages between 1-10~Myrs. The X-ray members have a very complex spatial distribution with some degree of subclustering, qualitatively similar to that of previously known members. The distribution of X-ray sources relative to IR-excess objects found with Spitzer is sometimes suggestive of sequential star formation, especially near the Gulf of Mexico region, probably triggered by the O5 star illuminating the whole region. Around this latter star no enhancement in the young star density is found, in agreement with previous results. We also determine the local optical/IR reddening law, and compute an updated reddening map of the entire region.


title: 'Viral' Turing Machines, Computation from Noise and Combinatorial Hierarchies

abstract: The interactive computation paradigm is reviewed and a particular example is extended to form the stochastic analog of a computational process via a transcription of a minimal Turing Machine into an equivalent asynchronous Cellular Automaton with an exponential waiting times distribution of effective transitions. Furthermore, a special toolbox for analytic derivation of recursive relations of important statistical and other quantities is introduced in the form of an Inductive Combinatorial Hierarchy.


title: Optimal design of experiments by combining coarse and fine measurements

abstract: In many contexts it is extremely costly to perform enough high quality experimental measurements to accurately parameterize a predictive quantitative model. However, it is often much easier to carry out large numbers of experiments that indicate whether each sample is above or below a given threshold. Can many such categorical or "coarse" measurements be combined with a much smaller number of high resolution or "fine" measurements to yield accurate models? Here, we demonstrate an intuitive strategy, inspired by statistical physics, wherein the coarse measurements are used to identify the salient features of the data, while the fine measurements determine the relative importance of these features. A linear model is inferred from the fine measurements, augmented by a quadratic term that captures the correlation structure of the coarse data. We illustrate our strategy by considering the problems of predicting the antimalarial potency and aqueous solubility of small organic molecules from their 2D molecular structure.

categories: physics.data-an cond-mat.soft cond-mat.stat-mech physics.chem-ph


title: Reply to: "TOV or OV? The Whole Story"

abstract: We argued previously that the well-known equation for hydrostatic equilibrium in a static spherically symmetric spacetime supported by an isotropic perfect fluid should be called the Oppenheimer-Volkoff (OV) equation, rather than the Tolman-Oppenheimer-Volkoff (TOV) equation, a later ArXiv manuscript has disagreed. Here we reply to that comment, reaffirming the original argument.


title: Exploring anomalous $hb\bar b$ and $h b\bar b\gamma$ couplings in the context of the LHC and an $e^+e^-$ collider

abstract: In the light of the 125 GeV Higgs ($h$) discovery at the Large Hadron Collider (LHC), one of the primary goals of the LHC and possible future colliders is to understand its interactions more precisely. Here we have studied the $h$-$b$-$\bar b$-$\gamma$ effective interaction terms arising out of gauge invariant dimension six operators in a model independent setting, as a potential source of new physics. Their role in some detectable final states have been compared with those coming from anomalous $h$-$b$-$\bar b$ interactions. We have considered the bounds coming from the existing collider and other low energy experimental data in order to derive constraints on the potential new physics couplings and predict possible collider signals for the two different new physics scenarios in the context of 14 TeV LHC and and a future $e^+e^-$ machine. We conclude that the anomalous $h$-$b$-$\bar b$-$\gamma$ coupling can be probed at the LHC at 14 TeV at the 3$\sigma$ level with an integrated luminosity of $\sim 2000~{\rm fb}^{-1}$, which an $e^+e^-$ collider can probe at the 3$\sigma$ level with $\sim 12(7)~{\rm fb}^{-1}$ at $\sqrt{s}=250(500)~{\rm GeV}$. It is also found that anomalous $h$-$b$-$\bar b$ interactions, subject to the existing LHC constraints, can not compete with the rates driven by $h$-$b$-$\bar b$-$\gamma$ effective interactions.


title: Nonclassical photon number distribution in a superconducting cavity under a squeezed drive

abstract: A superconducting qubit in the strong dispersive regime of a circuit quantum electrodynamics system is a powerful probe for microwave photons in a cavity mode. In this regime, a qubit spectrum is split into multiple peaks, with each peak corresponding to an individual photon number in the cavity (discrete ac Stark shift). Here, we measure the qubit spectrum in the cavity that is driven continuously with a squeezed vacuum field generated by a Josephson parametric amplifier. By fitting the qubit spectrum with a model which takes into account the finite qubit excitation power, the photon number distribution, which is dissimilar from the apparent peak area ratio in the spectrum, is determined. The photon number distribution shows the even-odd photon number oscillation and quantitatively fulfills Klyshko's criterion for the nonclassicality.