title: Can Data Generated by Connected Vehicles Enhance Safety? A proactive approach to intersection safety management

abstract: Traditionally, evaluation of intersection safety has been largely reactive, based on historical crash frequency data. However, the emerging data from Connected and Automated Vehicles (CAVs) can complement historical data and help in proactively identify intersections which have high levels of variability in instantaneous driving behaviors prior to the occurrence of crashes. Based on data from Safety Pilot Model Deployment in Ann Arbor, Michigan, this study developed a unique database that integrates intersection crash and inventory data with more than 65 million real-world Basic Safety Messages logged by 3,000 connected vehicles, providing a more complete picture of operations and safety performance of intersections. As a proactive safety measure and a leading indicator of safety, this study introduces location-based volatility (LBV), which quantifies variability in instantaneous driving decisions at intersections. LBV represents the driving performance of connected vehicle drivers traveling through a specific intersection. As such, by using coefficient of variation as a standardized measure of relative dispersion, LBVs are calculated for 116 intersections in Ann Arbor. To quantify relationships between intersection-specific volatilities and crash frequencies, rigorous fixed- and random-parameter Poisson regression models are estimated. While controlling for exposure related factors, the results provide evidence of statistically significant (5% level) positive association between intersection-specific volatility and crash frequencies for signalized intersections. The implications of the findings for proactive intersection safety management are discussed in the paper.


title: An Improved Algorithm for E-Generalization

abstract: E-generalization computes common generalizations of given ground terms w.r.t. a given equational background theory E. In 2005 [arXiv:1403.8118], we had presented a computation approach based on standard regular tree grammar algorithms, and a Prolog prototype implementation. In this report, we present algorithmic improvements, prove them correct and complete, and give some details of an efficiency-oriented implementation in C that allows us to handle problems larger by several orders of magnitude.


title: $L^p$ Christoffel-Minkowski problem: the case $1< p<k+1$

abstract: We consider a fully nonlinear partial differential equation associated to the intermediate $L^p$ Christoffel-Minkowski problem in the case $1<p<k+1$. We establish the existence of convex body with prescribed $k$-th even $p$-area measure on $\mathbb S^n$, under an appropriate assumption on the prescribed function. We construct examples to indicate certain geometric condition on the prescribed function is needed for the existence of smooth strictly convex body. We also obtain $C^{1,1}$ regularity estimates for admissible solutions of the equation when $ p\ge \frac{k+1}2$.


title: An Optimal Strategy for Accurate Bulge-to-disk Decomposition of Disk Galaxies

abstract: The development of two-dimensional (2D) bulge-to-disk decomposition techniques has shown their advantages over traditional one-dimensional (1D) techniques, especially for galaxies with non-axisymmetric features. However, the full potential of 2D techniques has yet to be fully exploited. Secondary morphological features in nearby disk galaxies, such as bars, lenses, rings, disk breaks, and spiral arms, are seldom accounted for in 2D image decompositions, even though some image-fitting codes, such as GALFIT, are capable of handling them. We present detailed, 2D multi-model and multi-component decomposition of high-quality $R$-band images of a representative sample of nearby disk galaxies selected from the Carnegie-Irvine Galaxy Survey, using the latest version of GALFIT. The sample consists of five barred and five unbarred galaxies, spanning Hubble types from S0 to Sc. Traditional 1D decomposition is also presented for comparison. In detailed case studies of the 10 galaxies, we successfully model the secondary morphological features. Through a comparison of best-fit parameters obtained from different input surface brightness models, we identify morphological features that significantly impact bulge measurements. We show that nuclear and inner lenses/rings and disk breaks must be properly taken into account to obtain accurate bulge parameters, whereas outer lenses/rings and spiral arms have a negligible effect. We provide an optimal strategy to measure bulge parameters of typical disk galaxies, as well as prescriptions to estimate realistic uncertainties of them, which will benefit subsequent decomposition of a larger galaxy sample.


title: Koml\'os-Major-Tusn\'ady approximations to increments of uniform empirical processes

abstract: The well-known Koml\'os-Major-Tusn\'ady inequalities [Z. Wahrsch. Verw. Gebiete 32 (1975) 111-131; Z. Wahrsch. Verw. Gebiete 34 (1976) 33-58] provide sharp inequalities to partial sums of iid standard exponential random variables by a sequence of standard Brownian motions. In this paper, we employ these results to establish Gaussian approximations to weighted increments of uniform empirical and quantile processes. This approach provides rates to the approximations which, among others, have direct applications to statistics of extreme values for randomly censored data.


title: Recovery of the singularities of a potential from backscattering data in general dimension

abstract: We prove that in dimension $n\ge 2$ the main singularities of a complex potential $q$ having a certain a priori regularity are contained in the Born approximation $q_{B}$ constructed from backscattering data. This is archived using a new explicit formula for the multiple dispersion operators in the Fourier transform side. We also show that ${q-q_{B}}$ can be up to one derivative more regular than $q$ in the Sobolev scale. On the other hand, we construct counterexamples showing that in general it is not possible to have more than one derivative gain, sometimes even strictly less, depending on the a priori regularity of $q$.


title: Efimov states near a Feshbach resonance and the limits of van der Waals universality at finite background scattering length

abstract: We calculate the spectrum of three-body Efimov bound states near a Feshbach resonance within a model which accounts both for the finite range of interactions and the presence of background scattering. The latter may be due to direct interactions in an open channel or a second overlapping Feshbach resonance. It is found that background scattering gives rise to substantial changes in the trimer spectrum as a function of the detuning away from a Feshbach resonance, in particular in the regime where the background channel supports Efimov states on its own. Compared to the situation with negligible background scattering, the regime where van der Waals universality applies is shifted to larger values of the resonance strength if the background scattering length is positive. For negative background scattering lengths, in turn, van der Waals universality extends to even small values of the resonance strength parameter, consistent with experimental results on Efimov states in $^{39}$K. Within a simple model, we show that short-range three-body forces do not affect van der Waals universality significantly. Repulsive three-body forces may, however, explain the observed variation between around $-8$ and $-10$ of the ratio between the scattering length where the first Efimov trimer appears and the van der Waals length.


title: Functional equations in algebra

abstract: We study flat deformations of quotients of a polynomial algebra in a class of graded commutative associative algebras. Functional equations and their solutions in terms of theta functions play important role in these studies. An analog of this theory in a fermionic case is also briefly discussed.


title: Sushi Dish - Object detection and classification from real images

abstract: In conveyor belt sushi restaurants, billing is a burdened job because one has to manually count the number of dishes and identify the color of them to calculate the price. In a busy situation, there can be a mistake that customers are overcharged or under-charged. To deal with this problem, we developed a method that automatically identifies the color of dishes and calculate the total price using real images. Our method consists of ellipse fitting and convol-utional neural network. It achieves ellipse detection precision 85% and recall 96% and classification accuracy 92%.


title: A short note on the joint entropy of n/2-wise independence

abstract: In this note, we prove a tight lower bound on the joint entropy of $n$ unbiased Bernoulli random variables which are $n/2$-wise independent. For general $k$-wise independence, we give new lower bounds by adapting Navon and Samorodnitsky's Fourier proof of the `LP bound' on error correcting codes. This counts as partial progress on a problem asked by Gavinsky and Pudl\'ak.


title: Compressed Sensing MRI Reconstruction using a Generative Adversarial Network with a Cyclic Loss

abstract: Compressed Sensing MRI (CS-MRI) has provided theoretical foundations upon which the time-consuming MRI acquisition process can be accelerated. However, it primarily relies on iterative numerical solvers which still hinders their adaptation in time-critical applications. In addition, recent advances in deep neural networks have shown their potential in computer vision and image processing, but their adaptation to MRI reconstruction is still in an early stage. In this paper, we propose a novel deep learning-based generative adversarial model, RefineGAN, for fast and accurate CS-MRI reconstruction. The proposed model is a variant of fully-residual convolutional autoencoder and generative adversarial networks (GANs), specifically designed for CS-MRI formulation; it employs deeper generator and discriminator networks with cyclic data consistency loss for faithful interpolation in the given under-sampled k-space data. In addition, our solution leverages a chained network to further enhance the reconstruction quality. RefineGAN is fast and accurate -- the reconstruction process is extremely rapid, as low as tens of milliseconds for reconstruction of a 256x256 image, because it is one-way deployment on a feed-forward network, and the image quality is superior even for extremely low sampling rate (as low as 10%) due to the data-driven nature of the method. We demonstrate that RefineGAN outperforms the state-of-the-art CS-MRI methods by a large margin in terms of both running time and image quality via evaluation using several open-source MRI databases.


title: Rapid green biosynthesis of Silvermontmorillonite nanocomposite using water extract of Satureja hortensis L and evaluation of the antibacterial properties

abstract: In this work, facile, rapid and green biosynthesis and characterization of silvermontmorillonite nanocomposite is reported at room temperature. Silver nanoparticles were prepared by water extract of Satureja hortensis as the reducing agent and MMT as interlamellar space for controlling the size of Ag-NPs. MMT was suspended in the aqueous AgNO3 solution, and after the absorption of silver ions, Ag+ was reduced to Ag0 using water extract of Satureja hortensis. Evaluations of the antibacterial properties are also reported. The nanocomposite was characterized by ultraviolet-visible spectroscopy, powder Xray diffraction, scanning electron microscopy and transmission electron microscopy . TEM study showed the formations of nanocomposite using water extract of Satureja hortensis L in the 4.88-26.70 nm range and an average particles size of 15.79 nm. In addition, XRD study showed that the particles have a face centered cubic structure. The nanocomposite showed antibacterial properties against Gram positive and negative bacteria.


title: Spectral Triples for the Variants of the Sierpinski Gasket

abstract: Fractal geometry is the study of sets which exhibit the same pattern at multiple scales. Developing tools to study these sets is of great interest. One step towards developing some of these tools is recognizing the duality between topological spaces and commutative $C^\ast$-algebras. When one lifts the commutativity axiom, one gets what are called noncommutative spaces and the study of noncommutative geometry. The tools built to study noncommutative spaces can in fact be used to study fractal sets. In what follows we will use the spectral triples of noncommutative geometry to describe various notions from fractal geometry. We focus on the fractal sets known as the harmonic Sierpinski gasket and the stretched Sierpinski gasket, and show that the spectral triples constructed by Christensen, Ivan, and Lapidus in 2008 and Lapidus and Sarhad in 2015, can recover the standard self-affine measure in the case of the harmonic Sierpinski gasket and the Hausdorff dimension, geodesic metric, and Hausdorff measure in the case of the stretched Sierpinski gasket.


title: Improved measurement of $^8$B solar neutrinos with 1.5 kt y of Borexino exposure

abstract: We report on an improved measurement of the $^8$B solar neutrino interaction rate with the Borexino experiment at the Laboratori Nazionali del Gran Sasso. Neutrinos are detected via their elastic scattering on electrons in a large volume of liquid scintillator. The measured rate of scattered electrons above 3 MeV of energy is $0.223\substack{+0.015 \\ -0.016}\,(stat)\,\substack{+0.006 \\ -0.006}\,(syst)$ cpd/100 t, which corresponds to an observed solar neutrino flux assuming no neutrino flavor conversion of $\Phi\substack{\rm ES \\ ^8\rm B}=2.57\substack{+0.17 \\ -0.18}(stat)\substack{+0.07\\ -0.07}(syst)\times$10$^6$ cm$^{-2}\,$s$^{-1}$. This measurement exploits the active volume of the detector in almost its entirety for the first time, and takes advantage of a reduced radioactive background following the 2011 scintillator purification campaign and of novel analysis tools providing a more precise modeling of the background. Additionally, we set a new limit on the interaction rate of solar $hep$ neutrinos, searched via their elastic scattering on electrons as well as their neutral current-mediated inelastic scattering on carbon, $^{12}$C($\nu,\nu'$)$^{12}$C* ($E_{\gamma}$= 15.1 MeV).


title: On the Continuity of the Topological Entropy of Non-autonomous Dynamical Systems

abstract: Let $M$ be a compact Riemannian manifold. The set $\text{F}^{r}(M)$ consisting of sequences $(f_{i})_{i\in\mathbb{Z}}$ of $C^{r}$-diffeomorphisms on $M$ can be endowed with the compact topology or with the strong topology. A notion of topological entropy is given for these sequences. I will prove this entropy is discontinuous at each sequence if we consider the compact topology on $\text{F}^{r}(M)$. On the other hand, if $ r\geq 1$ and we consider the strong topology on $\text{F}^{r}(M)$, this entropy is a continuous map.


title: A Method for Preparation and Readout of Polyatomic Molecules in Single Quantum States

abstract: Polyatomic molecular ions contain many desirable attributes of a useful quantum system, including rich internal degrees of freedom and highly controllable coupling to the environment. To date, the vast majority of state-specific experimental work on molecular ions has concentrated on diatomic species. The ability to prepare and readout polyatomic molecules in single quantum states would enable diverse experimental avenues not available with diatomics, including new applications in precision measurement, sensitive chemical and chiral analysis at the single molecule level, and precise studies of Hz-level molecular tunneling dynamics. While cooling the motional state of a polyatomic ion via sympathetic cooling with a laser cooled atomic ion is straightforward, coupling this motional state to the internal state of the molecule has proven challenging. Here we propose a new method for readout and projective measurement of the internal state of a trapped polyatomic ion. The method exploits the rich manifold of technically accessible rotational states in the molecule to realize robust state-preparation and readout with far less stringent engineering than quantum logic methods recently demonstrated on diatomic molecules. The method can be applied to any reasonably small (less than about 10 atoms) polyatomic ion with an anisotropic polarizability.


title: Boundary Control of a Nonhomogeneous Flexible Wing with Bounded Input Disturbances

abstract: This note deals with the boundary control problem of a nonhomogeneous flexible wing evolving under unsteady aerodynamic loads. The wing is actuated at its tip by flaps and is modeled by a distributed parameter system consisting of two coupled partial differential equations. Based on the proposed boundary control law, the well-posedness of the underlying Cauchy problem is first investigated by resorting to the semigroup theory. Then, a Lyapunov-based approach is employed to assess the stability of the closed-loop system in the presence of bounded input disturbances.


title: Meromorphic continuation of Selberg zeta functions with twists having non-expanding cusp monodromy

abstract: We initiate the study of Selberg zeta functions $Z_{\Gamma,\chi}$ for geometrically finite Fuchsian groups $\Gamma$ and finite-dimensional representations $\chi$ with non-expanding cusp monodromy. We show that for all choices of $(\Gamma,\chi)$, the Selberg zeta function $Z_{\Gamma,\chi}$ converges on some half-plane in $\mathbb{C}$. In addition, under the assumption that $\Gamma$ admits a strict transfer operator approach, we show that $Z_{\Gamma,\chi}$ extends meromorphically to all of $\mathbb{C}$.


title: Eisenstein series twisted with non-expanding cusp monodromies

abstract: Let $\Gamma$ be a geometrically finite Fuchsian group and suppose that $\chi\colon\Gamma\to\mathrm{GL}(V)$ is a finite-dimensional representation with non-expanding cusp monodromy. We show that the parabolic Eisenstein series for $\Gamma$ with twist $\chi$ converges on some half-plane. Further, we develop Fourier-type expansions for these Eisenstein series.


title: Nanoscale chiral valley-photon interface through optical spin-orbit coupling

abstract: The emergence of two-dimensional transition metal chalcogenide materials has sparked an intense activity in valleytronics since their valley information can be directly encoded and detected by using the spin angular momentum of light. For their practical applications such as on-chip valley logic gates and chip-to-chip valley transport, the encoding and processing of valley pseudospin using light should be extended to an integrated, on-chip nanophotonic system. Here, we successfully demonstrate, at room temperature, the valley-dependent directional coupling of light using a plasmonic nanowire-WS2 layers system. Our calculations show that the local transverse spin angular momentum of the mode of the plasmonic nanowire provides robust optical spin-path locking of up to 91 %. Experimentally we demonstrate that valley pseudospin in WS2 is coupled with optical spin of the same handedness and exhibits a high directional coupling efficiency up to 90 % to the plasmonic guided mode. The result opens up new avenues of controlling, detecting and processing valley and spin information with precise optical control at the nanoscale.


title: A Nova Outburst Powered by Shocks

abstract: Classical novae are runaway thermonuclear burning events on the surfaces of accreting white dwarfs in close binary star systems, sometimes appearing as new naked-eye sources in the night sky. The standard model of novae predicts that their optical luminosity derives from energy released near the hot white dwarf which is reprocessed through the ejected material. Recent studies with the Fermi Large Area Telescope have shown that many classical novae are accompanied by gigaelectronvolt gamma-ray emission. This emission likely originates from strong shocks, providing new insights into the properties of nova outflows and allowing them to be used as laboratories to study the unknown efficiency of particle acceleration in shocks. Here we report gamma-ray and optical observations of the Milky Way nova ASASSN-16ma, which is among the brightest novae ever detected in gamma-rays. The gamma-ray and optical light curves show a remarkable correlation, implying that the majority of the optical light comes from reprocessed emission from shocks rather than the white dwarf. The ratio of gamma-ray to optical flux in ASASSN-16ma directly constrains the acceleration efficiency of non-thermal particles to be ~0.005, favouring hadronic models for the gamma-ray emission. The need to accelerate particles up to energies exceeding 100 gigaelectronvolts provides compelling evidence for magnetic field amplification in the shocks.


title: Moduli spaces of low dimensional Lie superalgebras

abstract: In this paper, we study moduli spaces of low dimensional complex Lie superalgebras. We discover a similar pattern for the structure of these moduli spaces as we observed for ordinary Lie algebras, namely, that there is a stratification of the moduli space by projective orbifolds. The moduli spaces consist of some families as well as some singleton elements. The different strata are linked by jump deformations, which gives a uniques manner of decomposing the moduli space which is consistent with deformation theory.


title: Number of hidden states needed to physically implement a given conditional distribution

abstract: We consider the problem of how to construct a physical process over a finite state space $X$ that applies some desired conditional distribution $P$ to initial states to produce final states. This problem arises often in the thermodynamics of computation and nonequilibrium statistical physics more generally (e.g., when designing processes to implement some desired computation, feedback controller, or Maxwell demon). It was previously known that some conditional distributions cannot be implemented using any master equation that involves just the states in $X$. However, here we show that any conditional distribution $P$ can in fact be implemented---if additional "hidden" states not in $X$ are available. Moreover, we show that it is always possible to implement $P$ in a thermodynamically reversible manner. We then investigate a novel cost of the physical resources needed to implement a given distribution $P$: the minimal number of hidden states needed to do so. We calculate this cost exactly for the special case where $P$ represents a single-valued function, and provide an upper bound for the general case, in terms of the nonnegative rank of $P$. These results show that having access to one extra binary degree of freedom, thus doubling the total number of states, is sufficient to implement any $P$ with a master equation in a thermodynamically reversible way, if there are no constraints on the allowed form of the master equation. (Such constraints can greatly increase the minimal needed number of hidden states.) Our results also imply that for certain $P$ that can be implemented without hidden states, having hidden states permits an implementation that generates less heat.

categories: cond-mat.stat-mech cs.ET physics.bio-ph q-bio.QM


title: BSM Landscape

abstract: BSM landscape of motivated and unmotivated theories is overviewed with an emphasis on new developments guided by naturalness principle.


title: Hofstadter-Hubbard model with opposite magnetic fields: Bardeen-Cooper-Schrieffer pairing and superfluidity in the nearly flat butterfly bands

abstract: Despite the multi-band spectrum of the widely-known Hofstadter butterfly, it turns out that the pairing correlations of the time-reversal-symmetric Hofstadter-Hubbard model are well-described by a single order parameter that is uniform in real space. By exploiting a BCS mean-field theory for the nearly-flat butterfly-bands regime of low magnetic-flux limit, here we reveal a number of unusual superfluid properties both in the ground state and at finite temperatures. Our thorough analysis includes but is not limited to the order parameter, condensate and superfluid fractions, and the critical BCS and BKT transition temperatures.


title: Femtosecond formation dynamics of the spin Seebeck effect revealed by terahertz spectroscopy

abstract: Understanding the transfer of spin angular momentum is essential in modern magnetism research. A model case is the generation of magnons in magnetic insulators by heating an adjacent metal film. Here, we reveal the initial steps of this spin Seebeck effect with <27fs time resolution using terahertz spectroscopy on bilayers of ferrimagnetic yttrium-iron garnet and platinum. Upon exciting the metal with an infrared laser pulse, a spin Seebeck current $j_\textrm{s}$ arises on the same ~100fs time scale on which the metal electrons thermalize. This observation highlights that efficient spin transfer critically relies on carrier multiplication and is driven by conduction electrons scattering off the metal-insulator interface. Analytical modeling shows that the electrons' dynamics are almost instantaneously imprinted onto $j_\textrm{s}$ because their spins have a correlation time of only ~4fs and deflect the ferrimagnetic moments without inertia. Applications in material characterization, interface probing, spin-noise spectroscopy and terahertz spin pumping emerge.


title: The growth of Betti numbers and approximation theorems

abstract: These short lecture notes provide a brief introduction to the field of homology growth. They are composed out of two lectures, which I have given at the Borel seminar 2017 in Les Diablerets. We give a proof of L\"uck's approximation theorem, discuss generalizations and mention some related open problems. Then we discuss the growth of mod-$p$ Betti numbers, where many problems remain open. We take a closer look at the special case of $p$-adic analytic towers and discuss an approximation theorem due to Bergeron-Linnell-L\"uck-Sauer and Calegari-Emerton.


title: Understanding the Logical and Semantic Structure of Large Documents

abstract: Current language understanding approaches focus on small documents, such as newswire articles, blog posts, product reviews and discussion forum entries. Understanding and extracting information from large documents like legal briefs, proposals, technical manuals and research articles is still a challenging task. We describe a framework that can analyze a large document and help people to know where a particular information is in that document. We aim to automatically identify and classify semantic sections of documents and assign consistent and human-understandable labels to similar sections across documents. A key contribution of our research is modeling the logical and semantic structure of an electronic document. We apply machine learning techniques, including deep learning, in our prototype system. We also make available a dataset of information about a collection of scholarly articles from the arXiv eprints collection that includes a wide range of metadata for each article, including a table of contents, section labels, section summarizations and more. We hope that this dataset will be a useful resource for the machine learning and NLP communities in information retrieval, content-based question answering and language modeling.


title: Timing Observations of Diffusions

abstract: This paper addresses a problem in experimental design: We consider It\^o diffusions specified by some $\theta \in \mathbb{R}$ and assume that we are allowed to observe their sample paths only $n$ times before a terminal time $\tau < \infty$. We propose a policy for timing these observations to optimally estimate $\theta$. Our policy is adaptive (meaning it leverages earlier observations), and it maximizes the expected Fisher information for $\theta$ carried by the observations. In numerical studies, this design reduces the variation of estimated parameters by as much as 75% relative to observations spaced uniformly in time. The policy depends on the value of the parameter being estimated, so we also discuss strategies for incorporating Bayesian priors over $\theta$.


title: Modulated magnetic structure of Fe3PO7 as seen by 57Fe M\"ossbauer spectroscopy

abstract: The paper reports new results of the 57Fe M\"ossbauer measurements on Fe3PO4O3 powder sample recorded at various temperatures including the point of magnetic phase transition TN ~ 163K. The spectra measured above TN consist of quadrupole doublet with high quadrupole splitting of D300K ~ 1.10 mm/s, emphasizing that Fe3+ ions are located in crystal positions with a strong electric field gradient (EFG). In order to predict the sign and orientation of the main components of the EFG tensor we calculated monopole lattice contributions to the EFG. In the temperature range T < TN, the experimental spectra were fitted assuming that the electric hyperfine interactions are modulated when the Fe3+ spin (S) rotates with respect to the EFG axis and emergence of spatial anisotropy of the hyperfine field Hhf = S\~AI at 57Fe nuclei. These data were analyzed to estimate the components of the anisotropic hyperfine coupling tensor (\~A). The large anharmonicity parameter, m ~ 0.94, of the spiral spin structure results from easy-axis anisotropy in the plane of the iron spin rotation. The temperature evolution of the hyperfine field Hhf(T) was described by Bean-Rodbell model that takes into account that the exchange magnetic interactions are strong function of the lattice spacing. The obtained M\"ossbauer data are in qualitative agreement with previous neutron diffraction data for a modulated helical magnetic structure in strongly frustrated Fe3PO4O3.


title: Collision Energy Dependence of Moments of Net-Kaon Multiplicity Distributions at RHIC

abstract: Fluctuations of conserved quantities such as baryon number, charge, and strangeness are sensitive to the correlation length of the hot and dense matter created in relativistic heavy-ion collisions and can be used to search for the QCD critical point. We report the first measurements of the moments of net-kaon multiplicity distributions in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV. The collision centrality and energy dependence of the mean ($M$), variance ($\sigma^2$), skewness ($S$), and kurtosis ($\kappa$) for net-kaon multiplicity distributions as well as the ratio $\sigma^2/M$ and the products $S\sigma$ and $\kappa\sigma^2$ are presented. Comparisons are made with Poisson and negative binomial baseline calculations as well as with UrQMD, a transport model (UrQMD) that does not include effects from the QCD critical point. Within current uncertainties, the net-kaon cumulant ratios appear to be monotonic as a function of collision energy.


title: On the Viscous Camassa-Holm Equations with Fractional Diffusion

abstract: We study Cauchy problem of a class of viscous Camassa-Holm equations (or Lagrangian averaged Navier-Stokes equations) with fractional diffusion in both smooth bounded domains and in the whole space in two and three dimensions. Order of the fractional diffusion is assumed to be $2s$ with $s\in [n/4,1)$, which seems to be sharp for the validity of the main results of the paper; here $n=2,3$ is the dimension of space. We prove global well-posedness in $C_{[0,+\infty)}(D(A))\cap L^2_{[0,+\infty),loc}(D(A^{1+s/2}))$ whenever the initial data $u_0\in D(A)$, where $A$ is the Stokes operator. We also prove that such global solutions gain regularity instantaneously after the initial time. A bound on a higher-order spatial norm is also obtained.


title: Spin controlled atom-ion inelastic collisions

abstract: The control of the ultracold collisions between neutral atoms is an extensive and successful field of study. The tools developed allow for ultracold chemical reactions to be managed using magnetic fields, light fields and spin-state manipulation of the colliding particles among other methods. The control of chemical reactions in ultracold atom-ion collisions is a young and growing field of research. Recently, the collision energy and the ion electronic state were used to control atom-ion interactions. Here, we demonstrate spin-controlled atom-ion inelastic processes. In our experiment, both spin-exchange and charge-exchange reactions are controlled in an ultracold Rb-Sr$^+$ mixture by the atomic spin state. We prepare a cloud of atoms in a single hyperfine spin-state. Spin-exchange collisions between atoms and ion subsequently polarize the ion spin. Electron transfer is only allowed for (RbSr)$^+$ colliding in the singlet manifold. Initializing the atoms in various spin states affects the overlap of the collision wavefunction with the singlet molecular manifold and therefore also the reaction rate. We experimentally show that by preparing the atoms in different spin states one can vary the charge-exchange rate in agreement with theoretical predictions.


title: Estimation of interventional effects of features on prediction

abstract: The interpretability of prediction mechanisms with respect to the underlying prediction problem is often unclear. While several studies have focused on developing prediction models with meaningful parameters, the causal relationships between the predictors and the actual prediction have not been considered. Here, we connect the underlying causal structure of a data generation process and the causal structure of a prediction mechanism. To achieve this, we propose a framework that identifies the feature with the greatest causal influence on the prediction and estimates the necessary causal intervention of a feature such that a desired prediction is obtained. The general concept of the framework has no restrictions regarding data linearity; however, we focus on an implementation for linear data here. The framework applicability is evaluated using artificial data and demonstrated using real-world data.


title: Parity Games and Automata for Game Logic (Extended Version)

abstract: Parikh's game logic is a PDL-like fixpoint logic interpreted on monotone neighbourhood frames that represent the strategic power of players in determined two-player games. Game logic translates into a fragment of the monotone $\mu$-calculus, which in turn is expressively equivalent to monotone modal automata. Parity games and automata are important tools for dealing with the combinatorial complexity of nested fixpoints in modal fixpoint logics, such as the modal $\mu$-calculus. In this paper, we (1) discuss the semantics a of game logic over neighbourhood structures in terms of parity games, and (2) use these games to obtain an automata-theoretic characterisation of the fragment of the monotone $\mu$-calculus that corresponds to game logic. Our proof makes extensive use of structures that we call syntax graphs that combine the ease-of-use of syntax trees of formulas with the flexibility and succinctness of automata. They are essentially a graph-based view of the alternating tree automata that were introduced by Wilke in the study of modal $\mu$-calculus.


title: On the powers of the descent set statistic

abstract: We study the sum of the $r$th powers of the descent set statistic and how many small prime factors occur in these numbers. Our results depend upon the base $p$ expansion of $n$ and $r$.


title: Directional Cell Search Delay Analysis for Cellular Networks with Static Users

abstract: Cell search is the process for a user to detect its neighboring base stations (BSs) and make a cell selection decision. Due to the importance of beamforming gain in millimeter wave (mmWave) and massive MIMO cellular networks, the directional cell search delay performance is investigated. A cellular network with fixed BS and user locations is considered, so that strong temporal correlations exist for the SINR experienced at each BS and user. For Poisson cellular networks with Rayleigh fading channels, a closed-form expression for the spatially averaged mean cell search delay of all users is derived. This mean cell search delay for a noise-limited network (e.g., mmWave network) is proved to be infinite whenever the non-line-of-sight (NLOS) path loss exponent is larger than 2. For interference-limited networks, a phase transition for the mean cell search delay is shown to exist in terms of the number of BS antennas/beams $M$: the mean cell search delay is infinite when $M$ is smaller than a threshold and finite otherwise. Beam-sweeping is also demonstrated to be effective in decreasing the cell search delay, especially for the cell edge users.


title: A deep Parkes HI survey of the Sculptor group and filament: HI mass function and environment

abstract: We present the results of a deep survey of the nearby Sculptor group and the associated Sculptor filament taken with the Parkes 64-m radio telescope in the 21-cm emission line of neutral hydrogen. We detect 31 HI sources in the Sculptor group/filament, eight of which are new HI detections. We derive a slope of the HI mass function along the Sculptor filament of $\alpha = -1.10^{+0.20}_{-0.11}$, which is significantly flatter than the global mass function and consistent with the flat slopes previously found in other low-density group environments. Some physical process, such as star formation, photoionisation or ram-pressure stripping, must therefore be responsible for removing neutral gas predominantly from low-mass galaxies. All of our HI detections have a confirmed or tentative optical counterpart and are likely associated with luminous rather than 'dark' galaxies. Despite a column density sensitivity of about $4 \times 10^{17}~\mathrm{cm}^{-2}$, we do not find any traces of extragalactic gas or tidal streams, suggesting that the Sculptor filament is, at the current time, a relatively quiescent environment that has not seen any recent major interactions or mergers.


title: Non-Uniform Wavelet Sampling for RF Analog-to-Information Conversion

abstract: Feature extraction, such as spectral occupancy, interferer energy and type, or direction-of-arrival, from wideband radio-frequency~(RF) signals finds use in a growing number of applications as it enhances RF transceivers with cognitive abilities and enables parameter tuning of traditional RF chains. In power and cost limited applications, e.g., for sensor nodes in the Internet of Things, wideband RF feature extraction with conventional, Nyquist-rate analog-to-digital converters is infeasible. However, the structure of many RF features (such as signal sparsity) enables the use of compressive sensing (CS) techniques that acquire such signals at sub-Nyquist rates. While such CS-based analog-to-information (A2I) converters have the potential to enable low-cost and energy-efficient wideband RF sensing, they suffer from a variety of real-world limitations, such as noise folding, low sensitivity, aliasing, and limited flexibility. This paper proposes a novel CS-based A2I architecture called non-uniform wavelet sampling (NUWS). Our solution extracts a carefully-selected subset of wavelet coefficients directly in the RF domain, which mitigates the main issues of existing A2I converter architectures. For multi-band RF signals, we propose a specialized variant called non-uniform wavelet bandpass sampling (NUWBS), which further improves sensitivity and reduces hardware complexity by leveraging the multi-band signal structure. We use simulations to demonstrate that NUWBS approaches the theoretical performance limits of $\ell_1$-norm-based sparse signal recovery. We investigate hardware-design aspects and show ASIC measurement results for the wavelet generation stage, which highlight the efficacy of NUWBS for a broad range of RF feature extraction tasks in cost- and power-limited applications.


title: TARN: A SDN-based Traffic Analysis Resistant Network Architecture

abstract: Destination IP prefix-based routing protocols are core to Internet routing today. Internet autonomous systems (AS) possess fixed IP prefixes, while packets carry the intended destination AS's prefix in their headers, in clear text. As a result, network communications can be easily identified using IP addresses and become targets of a wide variety of attacks, such as DNS/IP filtering, distributed Denial-of-Service (DDoS) attacks, man-in-the-middle (MITM) attacks, etc. In this work, we explore an alternative network architecture that fundamentally removes such vulnerabilities by disassociating the relationship between IP prefixes and destination networks, and by allowing any end-to-end communication session to have dynamic, short-lived, and pseudo-random IP addresses drawn from a range of IP prefixes rather than one. The concept is seemingly impossible to realize in todays Internet. We demonstrate how this is doable today with three different strategies using software defined networking (SDN), and how this can be done at scale to transform the Internet addressing and routing paradigms with the novel concept of a distributed software defined Internet exchange (SDX). The solution works with both IPv4 and IPv6, whereas the latter provides higher degrees of IP addressing freedom. Prototypes based on OpenvSwitches (OVS) have been implemented for experimentation across the PEERING BGP testbed. The SDX solution not only provides a technically sustainable pathway towards large-scale traffic analysis resistant network (TARN) support, it also unveils a new business model for customer driven, customizable and trustable end-to-end network services.


title: Charged Particle Diffusion in a Magnetic Dipole Trap

abstract: When particles are magnetized, a diffusion process is influenced by the ambient magnetic field. While the entropy increases, the constancy of the magnetic moment puts a constraint. Here, we compare the E-cross-B diffusion caused by random fluctuations of the electric field in two different systems, the Penning-Malmberg trap and the magnetic dipole trap. A Fokker-Planck equation is derived by applying the ergodic ansatz on the invariant measure of the system. In the dipole magnetic field particles diffuse inward and accumulate in the higher magnetic field region, while, in a homogeneous magnetic field, particles diffuse out from the confinement region. The properties of analogous transport in a more general class of magnetic fields are also briefly discussed.


title: Novel gas target for Laser Wakefield Accelerators

abstract: A novel gas target for interactions between high power lasers and gaseous medium, especially for laser wakefield accelerators, has been designed, manufactured and characterized. The gas target has been designed to provide a uniform density profile along the central gas cell axis by combining a gas cell and slit nozzle. The gas density can be tuned from 10^16 atoms/cm^3 to 10^19 atoms/cm^3 and the gas target length can be varied from 0 to 10 cm, both changes can be made simultaneously while keeping the flat-top gas profile. The gas distributions inside the gas cell have been measured using interferometry and validated using computational fluid dynamics.


title: Strong Coupling of Alfv\'en and Fast Modes in Compressible Relativistic Magnetohydrodynamic Turbulence in Magnetically-Dominated Plasmas

abstract: In this paper, we report our detailed analysis of the new strong-coupling regime between Alfv\'en and fast modes in Poynting-dominated plasma turbulence, reported in our previous work Takamoto & Lazarian (2016), which is an important effect for many relativistic plasma phenomena, and calls for new theories of Poynting-dominated MHD turbulence. We performed numerical simulations of relativistic MHD turbulence in isothermal plasmas, and analyzed the ratio of fast to Alfv\'en mode energy. We found that the increase of the fast mode with the background $\sigma$-parameter can be observed even in isothermal plasma, showing that such a phenomena is universal in trans-Alfv\'enic turbulence in Poynting-dominated plasmas. To study the detailed energy conversion process, we also performed a series of simulations of decaying turbulence injecting pure Alfv\'en, fast, and slow modes, respectively, and investigated the development of the mode conversion from the each mode. We also found that the mode conversion between Alfv\'en and fast modes is nearly insensitive to the background temperature. Finally, we report a result of a simulation with initially fast mode dominated turbulence. It developed into a temporally strong-coupling regime, which is a strong evidence for the existence of our suggesting strong-coupling regime of fast and Alfv\'en modes. Our result suggests that the strong turbulence in Poynting-dominated plasma is very different from that in the non-relativistic plasma. It will also give an important guidance to studies of particle acceleration and non-thermal photon emission from Poynting-dominated plasma.

categories: astro-ph.HE astro-ph.SR physics.flu-dyn physics.plasm-ph


title: Machine learning methods for histopathological image analysis

abstract: Abundant accumulation of digital histopathological images has led to the increased demand for their analysis, such as computer-aided diagnosis using machine learning techniques. However, digital pathological images and related tasks have some issues to be considered. In this mini-review, we introduce the application of digital pathological image analysis using machine learning algorithms, address some problems specific to such analysis, and propose possible solutions.


title: The Topological Period-Index Problem over 8-Complexes, I

abstract: We study the Postnikov tower of the classifying space of a compact Lie group P(n,mn), which gives obstructions to lifting a topological Brauer class of period $n$ to a PU_{mn}-torsor, where the base space is a CW complex of dimension 8. Combined with the study of a twisted version of Atiyah-Hirzebruch spectral sequence, this solves the topological period-index problem for CW complexes of dimension 8.


title: Tropicalization of 1-tacnodal curves on toric surfaces

abstract: A degeneration of a singular curve on a toric surface, called a tropicalization, was constructed by E. Shustin. He classified the degeneration of 1-cuspidal curves using polyhedral complexes called tropical curves. In this paper, we define a tropical version of a 1-tacnodal curve, that is, a curve having exactly one singular point whose topological type is $A_3$, and by applying the tropicalization method, we classify tropical curves which correspond to 1-tacnodal curves.


title: The combinatorics of the colliding bullets problem

abstract: The finite colliding bullets problem is the following simple problem: consider a gun, whose barrel remains in a fixed direction; let $(V_i)_{1\le i\le n}$ be an i.i.d.\ family of random variables with uniform distribution on $[0,1]$; shoot $n$ bullets one after another at times $1,2,\dots, n$, where the $i$th bullet has speed $V_i$. When two bullets collide, they both annihilate. We give the distribution of the number of surviving bullets, and in some generalisation of this model. While the distribution is relatively simple (and we found a number of bold claims online), our proof is surprisingly intricate and mixes combinatorial and geometric arguments; we argue that any rigorous argument must very likely be rather elaborate.


title: Two-dimensional magneto-optical trap as a source for cold strontium atoms

abstract: We report on the realization of a transversely loaded two-dimensional magneto-optical trap serving as a source for cold strontium atoms. We analyze the dependence of the source's properties on various parameters, in particular the intensity of a pushing beam accelerating the atoms out of the source. An atomic flux exceeding $10^9\,\mathrm{atoms/s}$ at a rather moderate oven temperature of $500\,^\circ\mathrm{C}$ is achieved. The longitudinal velocity of the atomic beam can be tuned over several tens of m/s by adjusting the power of the pushing laser beam. The beam divergence is around $60$ mrad, determined by the transverse velocity distribution of the cold atoms. The slow atom source is used to load a three-dimensional magneto-optical trap realizing loading rates up to $10^9\,\mathrm{atoms/s}$ without indication of saturation of the loading rate for increasing oven temperature. The compact setup avoids undesired effects found in alternative sources like, e.g., Zeeman slowers, such as vacuum contamination and black-body radiation due to the hot strontium oven.


title: Residue Integrals and Waring's Formulas for a Class of Systems of Transcendental Equations in $\mathbb{C}^n$

abstract: The present article is focused on the study of a special class of systems of nonlinear transcendental equations for which classical algebraic and symbolic methods are inapplicable. For the purpose of the study of such systems, we develop a method for computing residue integrals with integration over certain cycles. We describe conditions under which the mentioned residue integrals coincide with power sums of the inverses to the roots of a system of equations (i.e., multidimensional Waring's formulas). As an application of the suggested method, we consider a problem of finding sums of multi-variable number series.


title: Graphs determined by their $A_{\alpha}$-spectra

abstract: Let $G$ be a graph with $n$ vertices, and let $A(G)$ and $D(G)$ denote respectively the adjacency matrix and the degree matrix of $G$. Define $$ A_{\alpha}(G)=\alpha D(G)+(1-\alpha)A(G) $$ for any real $\alpha\in [0,1]$. The collection of eigenvalues of $A_{\alpha}(G)$ together with multiplicities are called the \emph{$A_{\alpha}$-spectrum} of $G$. A graph $G$ is said to be \emph{determined by its $A_{\alpha}$-spectrum} if all graphs having the same $A_{\alpha}$-spectrum as $G$ are isomorphic to $G$. We first prove that some graphs are determined by its $A_{\alpha}$-spectrum for $0\leq\alpha<1$, including the complete graph $K_m$, the star $K_{1,n-1}$, the path $P_n$, the union of cycles and the complement of the union of cycles, the union of $K_2$ and $K_1$ and the complement of the union of $K_2$ and $K_1$, and the complement of $P_n$. Setting $\alpha=0$ or $\frac{1}{2}$, those graphs are determined by $A$- or $Q$-spectra. Secondly, when $G$ is regular, we show that $G$ is determined by its $A_{\alpha}$-spectrum if and only if the join $G\vee K_m$ is determined by its $A_{\alpha}$-spectrum for $\frac{1}{2}<\alpha<1$. Furthermore, we also show that the join $K_m\vee P_n$ is determined by its $A_{\alpha}$-spectrum for $\frac{1}{2}<\alpha<1$. In the end, we pose some related open problems for future study.


title: Cell contraction induces long-ranged stress stiffening in the extracellular matrix

abstract: Animal cells in tissues are supported by biopolymer matrices, which typically exhibit highly nonlinear mechanical properties. While the linear elasticity of the matrix can significantly impact cell mechanics and functionality, it remains largely unknown how cells, in turn, affect the nonlinear mechanics of their surrounding matrix. Here we show that living contractile cells are able to generate a massive stiffness gradient in three distinct 3D extracellular matrix model systems: collagen, fibrin, and Matrigel. We decipher this remarkable behavior by introducing Nonlinear Stress Inference Microscopy (NSIM), a novel technique to infer stress fields in a 3D matrix from nonlinear microrheology measurement with optical tweezers. Using NSIM and simulations, we reveal a long-ranged propagation of cell-generated stresses resulting from local filament buckling. This slow decay of stress gives rise to the large spatial extent of the observed cell-induced matrix stiffness gradient, which could form a mechanism for mechanical communication between cells.

categories: physics.bio-ph cond-mat.soft q-bio.QM q-bio.TO


title: The GALAH survey: Chemical Tagging of Star Clusters and New Members in the Pleiades

abstract: The technique of chemical tagging uses the elemental abundances of stellar atmospheres to `reconstruct' chemically homogeneous star clusters that have long since dispersed. The GALAH spectroscopic survey --which aims to observe one million stars using the Anglo-Australian Telescope -- allows us to measure up to 30 elements or dimensions in the stellar chemical abundance space, many of which are not independent. How to find clustering reliably in a noisy high-dimensional space is a difficult problem that remains largely unsolved. Here we explore t-distributed stochastic neighbour embedding (t-SNE) -- which identifies an optimal mapping of a high-dimensional space into fewer dimensions -- whilst conserving the original clustering information. Typically, the projection is made to a 2D space to aid recognition of clusters by eye. We show that this method is a reliable tool for chemical tagging because it can: (i) resolve clustering in chemical space alone, (ii) recover known open and globular clusters with high efficiency and low contamination, and (iii) relate field stars to known clusters. t-SNE also provides a useful visualization of a high-dimensional space. We demonstrate the method on a dataset of 13 abundances measured in the spectra of 187,000 stars by the GALAH survey. We recover 7 of the 9 observed clusters (6 globular and 3 open clusters) in chemical space with minimal contamination from field stars and low numbers of outliers. With chemical tagging, we also identify two Pleiades supercluster members (which we confirm kinematically), one as far as 6$^\circ$ -- one tidal radius away from the cluster centre.


title: Plum pudding random medium model of biological tissue toward remote microscopy from spectroscopic light scattering

abstract: Biological tissue has a complex structure and exhibits rich spectroscopic behavior. There is \emph{no} tissue model up to now able to account for the observed spectroscopy of tissue light scattering and its anisotropy. Here we present, \emph{for the first time}, a plum pudding random medium (PPRM) model for biological tissue which succinctly describes tissue as a superposition of distinctive scattering structures (plum) embedded inside a fractal continuous medium of background refractive index fluctuation (pudding). PPRM faithfully reproduces the wavelength dependence of tissue light scattering and attributes the "anomalous" trend in the anisotropy to the plum and the powerlaw dependence of the reduced scattering coefficient to the fractal scattering pudding. Most importantly, PPRM opens up a novel venue of quantifying the tissue architecture and microscopic structures on average from macroscopic probing of the bulk with scattered light alone without tissue excision. We demonstrate this potential by visualizing the fine microscopic structural alterations in breast tissue (adipose, glandular, fibrocystic, fibroadenoma, and ductal carcinoma) deduced from noncontact spectroscopic measurement.


title: Counting Unlabelled Chord Diagrams of Maximal Genus

abstract: Maximal chord diagrams up to all isomorphisms are enumerated. The enumerating formula is based on a bijection between rooted one-vertex one-face maps on locally orientable surfaces and a certain class of symmetric chord diagrams. This result extends the one of Cori and Marcus regarding maximal chord diagrams enumerated up to rotations.


title: An Extension of the Non-Inferior Set Estimation Algorithm for Many Objectives

abstract: This work proposes a novel multi-objective optimization approach that globally finds a representative non-inferior set of solutions, also known as Pareto-optimal solutions, by automatically formulating and solving a sequence of weighted sum method scalarization problems. The approach is called MONISE (Many-Objective NISE) because it represents an extension of the well-known non-inferior set estimation (NISE) algorithm, which was originally conceived to deal with two-dimensional objective spaces. The proposal is endowed with the following characteristics: (1) uses a mixed-integer linear programming formulation to operate in two or more dimensions, thus properly supporting many (i.e., three or more) objectives; (2) relies on an external algorithm to solve the weighted sum method scalarization problem to optimality; and (3) creates a faithful representation of the Pareto frontier in the case of convex problems, and a useful approximation of it in the non-convex case. Moreover, when dealing specifically with two objectives, some additional properties are portrayed for the estimated non-inferior set. Experimental results validate the proposal and indicate that MONISE is competitive, in convex and non-convex (combinatorial) problems, both in terms of computational cost and the overall quality of the non-inferior set, measured by the acquired hypervolume.


title: Note on the backwards uniqueness of mean curvature flow

abstract: In this note, we will show a backwards uniqueness theorem of the mean curvature flow with bounded second fundamental form in arbitrary codimension.


title: Non-rigid image registration using fully convolutional networks with deep self-supervision

abstract: We propose a novel non-rigid image registration algorithm that is built upon fully convolutional networks (FCNs) to optimize and learn spatial transformations between pairs of images to be registered. Different from most existing deep learning based image registration methods that learn spatial transformations from training data with known corresponding spatial transformations, our method directly estimates spatial transformations between pairs of images by maximizing an image-wise similarity metric between fixed and deformed moving images, similar to conventional image registration algorithms. At the same time, our method also learns FCNs for encoding the spatial transformations at the same spatial resolution of images to be registered, rather than learning coarse-grained spatial transformation information. The image registration is implemented in a multi-resolution image registration framework to jointly optimize and learn spatial transformations and FCNs at different resolutions with deep self-supervision through typical feedforward and backpropagation computation. Since our method simultaneously optimizes and learns spatial transformations for the image registration, our method can be directly used to register a pair of images, and the registration of a set of images is also a training procedure for FCNs so that the trained FCNs can be directly adopted to register new images by feedforward computation of the learned FCNs without any optimization. The proposed method has been evaluated for registering 3D structural brain magnetic resonance (MR) images and obtained better performance than state-of-the-art image registration algorithms.


title: Ultra-Low Energy Calibration of LUX Detector using $^{127}$Xe Electron Capture

abstract: We report an absolute calibration of the ionization yields($\textit{Q$_y$})$ and fluctuations for electronic recoil events in liquid xenon at discrete energies between 186 eV and 33.2 keV. The average electric field applied across the liquid xenon target is 180 V/cm. The data are obtained using low energy $^{127}$Xe electron capture decay events from the 95.0-day first run from LUX (WS2013) in search of Weakly Interacting Massive Particles (WIMPs). The sequence of gamma-ray and X-ray cascades associated with $^{127}$I de-excitations produces clearly identified 2-vertex events in the LUX detector. We observe the K- (binding energy, 33.2 keV), L- (5.2 keV), M- (1.1 keV), and N- (186 eV) shell cascade events and verify that the relative ratio of observed events for each shell agrees with calculations. The N-shell cascade analysis includes single extracted electron (SE) events and represents the lowest-energy electronic recoil $\textit{in situ}$ measurements that have been explored in liquid xenon.


title: Some theoretical results on tensor elliptical distribution

abstract: The multilinear normal distribution is a widely used tool in tensor analysis of magnetic resonance imaging (MRI). Diffusion tensor MRI provides a statistical estimate of a symmetric 2nd-order diffusion tensor, for each voxel within an imaging volume. In this article, tensor elliptical (TE) distribution is introduced as an extension to the multilinear normal (MLN) distribution. Some properties including the characteristic function and distribution of affine transformations are given. An integral representation connecting densities of TE and MLN distributions is exhibited that is used in deriving the expectation of any measurable function of a TE variate.


title: Global analysis of Skyrme forces with higher-order density dependence

abstract: The density dependent term in Skyrme forces is essential, which simulates three-body and many-body correlations beyond the low-momentum two-body interaction. We speculate that a single density term may be insufficient and a higher-order density dependent term is added. The present work investigates the influences of higher-order density dependencies based on extended UNEDF0 and SkM* forces. The global descriptions of nuclear masses and charge radii have been presented. Consequently the extended UNEDF0 force gives a global rms error on binding energies of 1.29 MeV. The influences on fission barriers and equation of state have also been investigated. The perspectives to improve Skyrme forces have also been discussed, including global center-of-mass corrections and Lipkin-Nogami pairing corrections.


title: An Implicit Gas-Kinetic Scheme for Turbulent Flow on Unstructured Hybrid Mesh

abstract: In this study, an implicit scheme for the gas-kinetic scheme (GKS) on the unstructured hybrid mesh is proposed. The Spalart-Allmaras (SA) one equation turbulence model is incorporated into the implicit gas-kinetic scheme (IGKS) to predict the effects of turbulence. The implicit macroscopic governing equations are constructed and solved by the matrix-free lower-upper symmetric-Gauss-Seidel (LU-SGS) method. To reduce the number of cells and computational cost, the hybrid mesh is applied. A modified non-manifold hybrid mesh data(NHMD) is used for both unstructured hybrid mesh and uniform grid. Numerical investigations are performed on different 2D laminar and turbulent flows. The convergence property and the computational efficiency of the present IGKS method are investigated. Much better performance are obtained compared with the standard explicit gas-kinetic scheme. Also, our numerical results are found to be in good agreement with experiment data and other numerical solutions, demonstrating the good applicability and high efficiency of the present IGKS for the simulations of laminar and turbulent flows.


title: Anisotropic Radial Layout for Visualizing Centrality and Structure in Graphs

abstract: This paper presents a novel method for layout of undirected graphs, where nodes (vertices) are constrained to lie on a set of nested, simple, closed curves. Such a layout is useful to simultaneously display the structural centrality and vertex distance information for graphs in many domains, including social networks. Closed curves are a more general constraint than the previously proposed circles, and afford our method more flexibility to preserve vertex relationships compared to existing radial layout methods. The proposed approach modifies the multidimensional scaling (MDS) stress to include the estimation of a vertex depth or centrality field as well as a term that penalizes discord between structural centrality of vertices and their alignment with this carefully estimated field. We also propose a visualization strategy for the proposed layout and demonstrate its effectiveness using three social network datasets.


title: Approximation of stable law in Wasserstein-1 distance by Stein's method

abstract: Let $n \in \mathbb N$, let $\zeta_{n,1},...,\zeta_{n,n}$ be a sequence of independent random variables with $\mathbb E \zeta_{n,i}=0$ and $\mathbb E |\zeta_{n,i}|<\infty$ for each $i$, and let $\mu$ be an $\alpha$-stable distribution having characteristic function $e^{-|\lambda|^{\alpha}}$ with $\alpha\in (1,2)$. Denote $S_{n}=\zeta_{n,1}+...+\zeta_{n,n}$ and its distribution by $\mathcal L(S_n)$, we bound the Wasserstein distance of $\mathcal L(S_{n})$ and $\mu$ essentially by an $L^{1}$ discrepancy between two kernels, this bound can be interpreted as a generalization of the Stein discrepancy (in $L^{2}$ sense) introduced by Ledoux, Nourdin and Peccati. More precisely, we prove the following inequality: \begin{equation} \begin{split} d_W\left(\mathcal L (S_n), \mu\right) \ \le C \left[\sum_{i=1}^n\int_{-N}^N \left|\frac{\mathcal K_\alpha(t,N)}n -\frac{ K_i(t,N)}{\alpha}\right| d t \ +\ \mathcal R_{N,n}\right], \end{split} \end{equation} where $d_{W}$ is the Wasserstein distance of probability measures, $\mathcal K_\alpha(t,N)$ is the kernel of a decomposition of the fractional Laplacian $\Delta^{\frac \alpha2}$, $ K_i(t,N)$ is a kernel introduced by Chen, Goldstein and Shao with a truncation which can be interpreted as an $L^1$ Stein kernel, and $\mathcal R_{N,n}$ is a small remainder. The integral term $$\sum_{i=1}^n\int_{-N}^N \left|\frac{\mathcal K_\alpha(t,N)}n -\frac{ K_i(t,N)}{\alpha}\right| d t$$ can be interpreted as an $L^{1}$ Stein discrepancy. As an application, we prove a general theorem of stable law convergence rate when $\zeta_{n,i}$ are i.i.d. and the distribution falls in the normal domain of attraction of $\mu$. We also study four examples with comparing our convergence rates and those known for these examples, among which the distribution in the second example is not in the normal domain of attraction of $\mu$.


title: Combining Spitzer parallax and Keck II adaptive optics imaging to measure the mass of a solar-like star orbited by a cold gaseous planet discovered by microlensing

abstract: To obtain accurate mass measurements for cold planets discovered by microlensing, it is usually necessary to combine light curve modeling with at least two lens mass-distance relations. Often, a constraint on the Einstein ring radius measurement is obtained from the caustic crossing time: This is supplemented by secondary constraints such as precise parallax measurements and/or measures of the lens luminosity using high angular resolution observations. We resolved the source+lens star from sub-arcsecond blends in H band using adaptive optics (AO) observations with NIRC2 mounted on Keck II telescope. We identify additional flux, coincident with the source to within 160 mas. We estimate the potential contributions to this blended light (chance-aligned star, additional companion to the lens or to the source) and find that 85 % of of the NIR flux is due to the lens star at H_L=16.63 +- 0.06 and K_L=16.46 +- 0.06. We combined the parallax constraint and the AO constraint to derive the physical parameters of the system. The lensing system is composed of a mid-late type G main sequence star of M_L=0.89 +- 0.05 Mo located at D_L = 3.6 +- 0.3 kpc in the Galactic disk. Taking the mass ratio and projected separation from the original study leads to a planet of M_p= 0.64 +- 0.044 M_Jupiter at 3.48 +- 0.22 AU. Excellent parallax measurement from simultaneous ground-space observations have been obtained on the microlensing event OGLE-2014-BLG-0124, but it is only when they are combined with ~ 30 min of Keck II AO observations that the physical parameters of the host star are well measured.


title: An Ore-type Condition for Large $k$-factor and Disjoint Perfect Matchings

abstract: Win [\emph{J. Graph Theory} {\bf 6}(1982), 489--492] conjectured that a graph $G$ on $n$ vertices contains $k$ disjoint perfect matchings, if the degree sum of any two nonadjacent vertices is at least $n+k-2$, where $n$ is even and $n\geq k+2$. In this paper, we prove that Win's conjecture is true for $k\geq n/2$, where $n$ is sufficiently large. To show this result, we prove a theorem on $k$-factor in a graph under some Ore-type condition. Our main tools include Tutte's $k$-factor theorem, the Karush-Kuhn-Tucker theorem on convex optimization, and the solution to the longstanding 1-factor decomposition conjecture.


title: Charmed Baryon Weak Decays with SU(3) Flavor Symmetry

abstract: We study the semileptonic and non-leptonic charmed baryon decays with $SU(3)$ flavor symmetry, where the charmed baryons can be ${\bf B}_{c}=(\Xi_c^0,\Xi_c^+,\Lambda_c^+)$, ${\bf B}'_{c}=(\Sigma_c^{(++,+,0)},\Xi_{c}^{\prime(+,0)},\Omega_c^0)$, ${\bf B}_{cc}=(\Xi_{cc}^{++},\Xi_{cc}^+,\Omega_{cc}^+)$, or ${\bf B}_{ccc}=\Omega^{++}_{ccc}$. With ${\bf B}_n^{(\prime)}$ denoted as the baryon octet (decuplet), we find that the ${\bf B}_{c}\to {\bf B}'_n\ell^+\nu_\ell$ decays are forbidden, while the $\Omega_c^0\to \Omega^-\ell^+\nu_\ell$, $\Omega_{cc}^+\to\Omega_c^0\ell^+\nu_\ell$, and $\Omega_{ccc}^{++}\to \Omega_{cc}^+\ell^+\nu_\ell$ decays are the only existing Cabibbo-allowed modes for ${\bf B}'_{c}\to {\bf B}'_n\ell^+\nu_\ell$, ${\bf B}_{cc}\to {\bf B}'_c\ell^+\nu_\ell$, and ${\bf B}_{ccc}\to {\bf B}_{cc}^{(\prime)}\ell^+\nu_\ell$, respectively. We predict the rarely studied ${\bf B}_{c}\to {\bf B}_n^{(\prime)}M$ decays, such as ${\cal B}(\Xi_c^0\to\Lambda^0\bar K^0,\,\Xi_c^+\to\Xi^0\pi^+)=(8.3\pm 0.9,8.0\pm 4.1)\times 10^{-3}$ and ${\cal B}(\Lambda_c^+\to \Delta^{++}\pi^-,\,\Xi_c^0\to\Omega^- K^+)=(5.5\pm 1.3,4.8\pm 0.5)\times 10^{-3}$. For the observation, the doubly and triply charmed baryon decays of $\Omega_{cc}^{+}\to \Xi_c^+\bar K^0$, $\Xi_{cc}^{++}\to (\Xi_c^+\pi^+$, $\Sigma_c^{++}\bar K^0)$, and $\Omega_{ccc}^{++}\to (\Xi_{cc}^{++}\bar K^0,\Omega_{cc}^+\pi^+,\Xi_c^+ D^+)$ are the favored Cabibbo-allowed decays, which are accessible to the BESIII and LHCb experiments.


title: Large time behavior of solutions of the heat equation with inverse square potential

abstract: Let $L:=-\Delta+V$ be a nonnegative Schr\"odinger operator on $L^2({\bf R}^N)$, where $N\ge 2$ and $V$ is a radially symmetric inverse square potential. In this paper we assume either $L$ is subcritical or null-critical and we establish a method for obtaining the precise description of the large time behavior of $e^{-tL}\varphi$, where $\varphi\in L^2({\bf R}^N,e^{|x|^2/4}\,dx)$.


title: Fixed function and its application to Medical Science

abstract: In the present paper, the concept of contraction has been extended in a refined manner by introducing $\mathfrak{D}$-Contraction defined on a family $\mathfrak{F}$ of bounded functions. Also, a new notion of fixed function has been introduced for a metric space. Some fixed function theorems along with illustrative examples have also been given to verify the effectiveness of our results. In addition, an application to medical science has also been presented. This application is based on best approximation of treatment plan for tumor patients getting intensity modulated radiation therapy(IMRT). In this technique, a proper DDC matrix truncation has been used that significantly improves accuracy of results. In 2013, Z. Tian \textit{et al.} presented a fluence map optimization(FMO) model for dose calculation by splitting the DDC matrix into two components on the basis of a threshold intensity value. Following this concept, a sequence of functions can be constructed through the presented results which contains different dose distributions corresponding to different patients and finally it converges to a fixed function. The fixed function obtained in this application represents the suitable doses of a number of tumor patients at the same time. A nice explanation has also been given to check the authenticity of the application and the results.


title: Critical P\'olya urn

abstract: We propose a variant model of P{\'o}lya urn process, where the dynamics consist of two competing elements namely, suppression of growth and enhancement of dormant character. Here the level of such features are controlled by an internal parameter in a random manner, so as the system self-organizes and characteristic observables exhibit scale invariance suggesting the existence of criticality. Varying the internal control parameter, one can explain interesting universality classes for avalanche-type events. We also discuss different versions of the model. It is pointed that such an underlying mechanism has earlier been noted to operate in a neural network.


title: Phase Diagrams of Three-Dimensional Anderson and Quantum Percolation Models using Deep Three-Dimensional Convolutional Neural Network

abstract: The three-dimensional Anderson model is a well-studied model of disordered electron systems that shows the delocalization--localization transition. As in our previous papers on two- and three-dimensional (2D, 3D) quantum phase transitions [J. Phys. Soc. Jpn. {\bf 85}, 123706 (2016), {\bf 86}, 044708 (2017)], we used an image recognition algorithm based on a multilayered convolutional neural network. However, in contrast to previous papers in which 2D image recognition was used, we applied 3D image recognition to analyze entire 3D wave functions. We show that a full phase diagram of the disorder-energy plane is obtained once the 3D convolutional neural network has been trained at the band center. We further demonstrate that the full phase diagram for 3D quantum bond and site percolations can be drawn by training the 3D Anderson model at the band center.


title: From Review to Rating: Exploring Dependency Measures for Text Classification

abstract: Various text analysis techniques exist, which attempt to uncover unstructured information from text. In this work, we explore using statistical dependence measures for textual classification, representing text as word vectors. Student satisfaction scores on a 3-point scale and their free text comments written about university subjects are used as the dataset. We have compared two textual representations: a frequency word representation and term frequency relationship to word vectors, and found that word vectors provide a greater accuracy. However, these word vectors have a large number of features which aggravates the burden of computational complexity. Thus, we explored using a non-linear dependency measure for feature selection by maximizing the dependence between the text reviews and corresponding scores. Our quantitative and qualitative analysis on a student satisfaction dataset shows that our approach achieves comparable accuracy to the full feature vector, while being an order of magnitude faster in testing. These text analysis and feature reduction techniques can be used for other textual data applications such as sentiment analysis.


title: How vortices and shocks provide for a flux loop in two-dimensional compressible turbulence

abstract: Large-scale turbulence in fluid layers and other quasi-two-dimensional compressible systems consists of planar vortices and waves. Separately, wave turbulence usually produces a direct energy cascade, while solenoidal planar turbulence transports energy to large scales by an inverse cascade. Here, we consider turbulence at finite Mach numbers when the interaction between acoustic waves and vortices is substantial. We employ solenoidal pumping at intermediate scales and show how both direct and inverse energy cascades are formed starting from the pumping scale. We show that there is an inverse cascade of kinetic energy up to a scale $\ell$, where a typical velocity reaches the speed of sound; this creates shock waves, which provide for a compensating direct cascade. When the system size is less than $\ell$, the steady state contains a system-size pair of long-living condensate vortices connected by a system of shocks. Thus turbulence in fluid layers processes energy via a loop: Most energy first goes to large scales via vortices and is then transported by waves to small-scale dissipation.

categories: physics.flu-dyn astro-ph.GA nlin.CD physics.ao-ph


title: Switching of the topologically trivial and non-trivial quantum phase transitions in compressed 1T-TiTe2: Experiments and Theory

abstract: We report the structural, vibrational and electrical transport properties up to 16 GPa of the 1T-TiTe2, a prominent layered 2D system, which is predicted to show a series of topologically trivial - nontrivial transitions under hydrostatic compression. We clearly show signatures of two iso-structural transition at 2 GPa and 4 GPa obtained from the minima in c/a ratio concomitant with the phonon linewidth anomalies of Eg and A1g modes at around the same pressures, providing strong indication of unusual electron-phonon coupling associated to these transitions. Resistivity presents nonlinear behavior over similar pressure ranges providing a strong indication of the electronic origin of these pressure driven isostructural transitions. Our data thus provide clear evidences of topological changes at A and L point of the Brillouin zone predicted to be present in the compressed 1T-TiTe2. Between 4 GPa and 8 GPa, the c/a ratio shows a plateau suggesting a transformation from an anisotropic 2D layer to a quasi 3D crystal network. First principles calculations suggest that the 2D to quasi 3D evolution without any structural phase transitions is mainly due to the increased interlayer Te-Te interactions (bridging) via the charge density overlap. In addition to the pressure dependent isostructural phase transitions, our data also evidences the occurrence of a first order structural phase transition from the trigonal (P-3m1) phase at higher pressures. We estimate the start of this structural phase transition to be 8 GPa and the symmetric of the new high-pressure phase to be monoclinic (C2/m).


title: Revealing correlation effect of Co 3d electrons in La$_3$Co$_4$Sn$_{13}$ and Ce$_3$Co$_4$Sn$_{13}$ by infrared spectroscopy study

abstract: We report resistivity, specific heat and optical spectroscopy measurements on single crystal samples of La$_3$Co$_4$Sn$_{13}$ and Ce$_3$Co$_4$Sn$_{13}$. We observed clear temperature-induced spectral weight suppression below 4000 \cm for both compounds in conductivity spectra $\sigma_1(\omega)$, indicating the progressive formation of gap-like features with decreasing temperature. The suppressed spectral weight transfers mostly to the higher energy region. The observation reflects the presence of correlation effect in the compounds. We attribute the correlation effect to the Co 3d electrons.


title: An unknotting index for virtual knots

abstract: In this paper we introduce the notion of an unknotting index for virtual knots. We give some examples of computation by using writhe invariants, and discuss a relationship between the unknotting index and the virtual knot module. In particular, we show that for any non-negative integer $n$ there exists a virtual knot whose unknotting index is $(1,n)$.


title: Microstructures and Dynamics of Tetraalkylphosphonium Chloride Ionic Liquids

abstract: Atomistic simulations have been performed to investigate the effect of aliphatic chain length in tetraalkylphosphonium cations on liquid morphologies, microscopic ionic structures and dynamical properties of tetraalkylphosphonium chloride ionic liquids. The liquid morphologies are characterized by sponge-like interpenetrating polar and apolar networks in ionic liquids consisting of tetraalkylphosphonium cations with short aliphatic chains. The lengthening aliphatic chains in tetraalkylphosphonium cations leads to polar domains consisting of chloride anions and central polar groups in cations being partially or totally segregated in ionic liquid matrices due to a progressive expansion of apolar domains in between. Prominent polarity alternation peaks and adjacency correlation peaks are observed at low and high $q$ range in total X-ray scattering structural functions, respectively, and their peak positions gradually shift to lower q values with lengthening aliphatic chains in tetraalkylphosphonium cations. The charge alternation peaks registered in intermediate q range exhibit complicated tendencies due to the complete cancellations of peaks and anti-peaks in partial structural functions for ionic subcomponents. The particular microstructures and liquid morphologies in tetraalkylphosphonium chloride ionic liquids intrinsically contribute to distinct dynamics characterized by translational diffusion coefficients, van Hove correlation functions, and non-Gaussian parameters for ionic species in heterogeneous ionic environment. The increase of aliphatic chain length in tetraalkylphosphonium cations leads to concomitant shift of van Hove correlation functions and non-Gaussian parameters to larger radial distances and longer timescales, respectively, indicating the enhanced translational dynamical heterogeneities of tetraalkylphosphonium cations and the corresponding chloride anions.


title: Hidden-charm and bottom meson-baryon molecules coupled with five-quark states

abstract: In this paper, we investigate the hidden-charm pentaquarks as $\bar{D}^{(\ast)}\Lambda_{\rm c}$ and $\bar{D}^{(\ast)}\Sigma^{(\ast)}_{\rm c}$ molecules coupled to the five-quark states. Furthermore, we extend our calculations to the hidden-bottom sector. The coupling to the five-quark states is treated as the short range potential, where the relative strength for the meson-baryon channels is determined by the structure of the five-quark states. We found that resonant and/or bound states appear in both the charm and bottom sectors. The five-quark state potential turned out to be attractive and, for this reason, it plays an important role to produce these states. In the charm sector, we need the five-quark potential in addition to the pion exchange potential in producing bound and resonant states, whereas, in the bottom sector, the pion exchange interaction is strong enough to produce states. Thus, from this investigation, it emerges that the hidden-bottom pentaquarks are more likely to form than their hidden-charm counterparts; for this reason, we suggest that the experimentalists should look for states in the bottom sector.


title: On the distribution of divisors of monic polynomials over function fields

abstract: This paper deals with function field analogues of famous theorems of Laudau which counted the number of integers which have $t$ prime factors and R. Hall which researched the distribution of divisors of integers in residue classes.\;We extend the Selberg-Delange method to handle the following problems.\;The number of monic polynomials with degree $n$ have $t$ irreducible factors;\;The number of monic polynomials with degree $n$ in some residue classes have $t$ irreducible factors and the residue classes distribution of divisors of monic polynomial.\;


title: The electromagnetic waves propagation in unmagnetized plasma media using parallelized finite-difference time-domain method

abstract: The finite-difference time-domain (FDTD) method has been commonly utilized to simulate the electromagnetic (EM) waves propagation in the plasma media. However, the FDTD method may bring about extra run-time on concerning computationally large and complicated EM problems. Fortunately, the FDTD method is easy to parallelize. Besides, GPU has been widely used for parallel computing due to its unique SPMD (Single Program Multiple Data) architecture. In this paper, we represent the parallel Runge-Kutta exponential time differencing scheme FDTD (RKETD) method for the unmagnetized plasma implemented on GPU. The detailed flowchart of parallel RKETD-FDTD method is described. The accuracy and acceleration performance of the proposed parallel RKETD-FDTD method implemented on GPU are substantiated by calculating the reflection and transmission coefficients for one-dimensional unmagnetized plasma slab. The results indicate that the numerical precision of the parallel RKETD-FDTD scheme is consistent with that of the code implemented on CPU. The computation efficiency is greatly improved compared with merely CPU-based serial RKETD-FDTD method. Moreover, the comparisons of the performance of CUDA-based GPU parallel program, OpenMP (Open Multi-Processing)-based CPU parallel program, and single-CPU serial program on the same host computer are done. Compared with the serial program, both parallel programs get good results, while GPU-based parallel program gains better result.


title: An efficient algorithm to compute the genus of discrete surfaces and applications to turbulent flows

abstract: A simple and efficient algorithm to numerically compute the genus of surfaces of three-dimensional objects using the Euler characteristic formula is presented. The algorithm applies to objects obtained by thresholding a scalar field in a structured-collocated grid, and does not require any triangulation of the data. This makes the algorithm fast, memory-efficient and suitable for large datasets. Applications to the characterization of complex surfaces in turbulent flows are presented to illustrate the method.


title: Multiscale analysis of the topological invariants in the logarithmic region of turbulent channels at $Re_{\tau}=932$

abstract: The invariants of the velocity gradient tensor, $R$ and $Q$, and their enstrophy and strain components are studied in the logarithmic layer of an incompressible turbulent channel flow. The velocities are filtered in the three spatial directions and the results analyzed at different scales. We show that the $R$--$Q$ plane does not capture the changes undergone by the flow as the filter width increases, and that the enstrophy/enstrophy-production and strain/strain-production planes represent better choices. We also show that the conditional mean trajectories may differ significantly from the instantaneous behavior of the flow since they are the result of an averaging process where the mean is 3-5 times smaller than the corresponding standard deviation. The orbital periods in the $R$--$Q$ plane are shown to be independent of the intensity of the events, and of the same order of magnitude than those in the enstrophy/enstrophy-production and strain/strain-production planes. Our final goal is to test whether the dynamics of the flow are self-similar in the inertial range, and the answer turns out to be no. The mean shear is found to be responsible for the absence of self-similarity and progressively controls the dynamics of the eddies observed as the filter width increases. However, a self-similar behavior emerges when the calculations are repeated for the fluctuating velocity gradient tensor. Finally, the turbulent cascade in terms of vortex stretching is considered by computing the alignment of the vorticity at a given scale with the strain at a different one. These results generally support a non-negligible role of the phenomenological energy-cascade model formulated in terms of vortex stretching.


title: Making Bidirected Graphs Strongly Connected

abstract: We consider problems to make a given bidirected graph strongly connected with minimum cardinality of additional signs or additional arcs. For the former problem, we show the minimum number of additional signs and give a linear-time algorithm for finding an optimal solution. For the latter problem, we give a linear-time algorithm for finding a feasible solution whose size is equal to the obvious lower bound or more than that by one.


title: Lower bounds on the lifting degree of single-edge and multiple-edge QC-LDPC codes by difference matrices

abstract: In this paper, we define two matrices named as "difference matrices", denoted by $D$ and $DD$ which significantly contribute to achieve regular single-edge QC-LDPC codes with the shortest length and the certain girth as well as regular and irregular multiple-edge QC-LDPC codes. Making use of these matrices, we obtain necessary and sufficient conditions to have single-edge $(m,n)$-regular QC-LDPC codes with girth 6, 10 and 12. Additionally, for girth 6, we achieve all non-isomorphic codes with the minimum lifting degree, $N$, for $m=4$ and $5\leq n\leq 11$, and present an exponent matrix for each minimum distance. For girth 10, we provide a lower bound on the lifting degree which is tighter than the existing bound. More important, for an exponent matrix whose first row and first column are all-zero, we demonstrate that the non-existence of 8-cycles proves the non-existence of 6-cycles related to the first row of the exponent matrix too. All non-isomorphic QC-LDPC codes with girth 10 and $n=5,6$ whose numbers are more than those presented in the literature are provided. For $n=7,8$ we decrease the lifting degrees from 159 and 219 to 145 and 211, repectively. For girth 12, a lower bound on the lifting degree is achieved. For multiple-edge category, for the first time a lower bound on the lifting degree for both regular and irregular QC-LDPC codes with girth 6 is achieved. We also demonstrate that the achieved lower bounds on multiple-edge $(4,n)$-regular QC-LDPC codes with girth 6 are tight and the resultant codes have shorter length compared to their counterparts in single-edge codes. Additionally, difference matrices help to reduce the conditions of considering 6-cycles from seven states to five states. We obtain multiple-edge $(4,n)$-regular QC-LDPC codes with girth 8 and $n=4,6,8$ with the shortest length.


title: Analysing Mutual Exclusion using Process Algebra with Signals

abstract: In contrast to common belief, the Calculus of Communicating Systems (CCS) and similar process algebras lack the expressive power to accurately capture mutual exclusion protocols without enriching the language with fairness assumptions. Adding a fairness assumption to implement a mutual exclusion protocol seems counter-intuitive. We employ a signalling operator, which can be combined with CCS, or other process calculi, and show that this minimal extension is expressive enough to model mutual exclusion: we confirm the correctness of Peterson's mutual exclusion algorithm for two processes, as well as Lamport's bakery algorithm, under reasonable assumptions on the underlying memory model. The correctness of Peterson's algorithm for more than two processes requires stronger, less realistic assumptions on the underlying memory model.


title: Bisimulation and Hennessy-Milner Logic for Generalized Synchronization Trees

abstract: In this work, we develop a generalization of Hennessy-Milner Logic (HML) for Generalized Synchronization Trees (GSTs) that we call Generalized Hennessy Milner Logic (GHML). Importantly, this logic suggests a strong relationship between (weak) bisimulation for GSTs and ordinary bisimulation for Synchronization Trees (STs). We demonstrate that this relationship can be used to define the GST analog for image-finiteness of STs. Furthermore, we demonstrate that certain maximal Hennessy-Milner classes of STs have counterparts in maximal Hennessy-Milner classes of GSTs with respect to GST weak bisimulation. We also exhibit some interesting characteristics of these maximal Hennessy-Milner classes of GSTs.


title: Reversing Imperative Parallel Programs

abstract: We propose an approach and a subsequent extension for reversing imperative programs. Firstly, we produce both an augmented version and a corresponding inverted version of the original program. Augmentation saves reversal information into an auxiliary data store, maintaining segregation between this and the program state, while never altering the data store in any other way than that of the original program. Inversion uses this information to revert the final program state to the state as it was before execution. We prove that augmentation and inversion work as intended, and illustrate our approach with several examples. We also suggest a modification to our first approach to support non-communicating parallelism. Execution interleaving introduces a number of challenges, each of which our extended approach considers. We define annotation and redefine inversion to use a sequence of statement identifiers, making the interleaving order deterministic in reverse.


title: Using Session Types for Reasoning About Boundedness in the Pi-Calculus

abstract: The classes of depth-bounded and name-bounded processes are fragments of the pi-calculus for which some of the decision problems that are undecidable for the full calculus become decidable. P is depth-bounded at level k if every reduction sequence for P contains successor processes with at most k active nested restrictions. P is name-bounded at level k if every reduction sequence for P contains successor processes with at most k active bound names. Membership of these classes of processes is undecidable. In this paper we use binary session types to decise two type systems that give a sound characterization of the properties: If a process is well-typed in our first system, it is depth-bounded. If a process is well-typed in our second, more restrictive type system, it will also be name-bounded.


title: Distributive Laws for Monotone Specifications

abstract: Turi and Plotkin introduced an elegant approach to structural operational semantics based on universal coalgebra, parametric in the type of syntax and the type of behaviour. Their framework includes abstract GSOS, a categorical generalisation of the classical GSOS rule format, as well as its categorical dual, coGSOS. Both formats are well behaved, in the sense that each specification has a unique model on which behavioural equivalence is a congruence. Unfortunately, the combination of the two formats does not feature these desirable properties. We show that monotone specifications - that disallow negative premises - do induce a canonical distributive law of a monad over a comonad, and therefore a unique, compositional interpretation.


title: Hypothesis Testing based Intrinsic Evaluation of Word Embeddings

abstract: We introduce the cross-match test - an exact, distribution free, high-dimensional hypothesis test as an intrinsic evaluation metric for word embeddings. We show that cross-match is an effective means of measuring distributional similarity between different vector representations and of evaluating the statistical significance of different vector embedding models. Additionally, we find that cross-match can be used to provide a quantitative measure of linguistic similarity for selecting bridge languages for machine translation. We demonstrate that the results of the hypothesis test align with our expectations and note that the framework of two sample hypothesis testing is not limited to word embeddings and can be extended to all vector representations.


title: Extended Phase Graph formalism for systems with Magnetization Transfer and Chemical Exchange

abstract: An Extended Phase Graph framework for modelling systems with exchange or magnetization transfer (MT) is proposed. The framework, referred to as EPG-X, models coupled two-compartment systems by describing each compartment with separate phase graphs that exchange during evolution periods. There are two variants: EPG-X(BM) for systems governed by the Bloch-McConnell equations; and EPG-X(MT) for the pulsed MT formalism. For the MT case the "bound" protons have no transverse components so their phase graph consists only longitudinal states. EPG-X was used to model steady-state gradient echo imaging, MT effects in multislice Turbo Spin Echo imaging, multiecho CPMG for multicomponent T2 relaxometry and transient variable flip angle gradient echo imaging of the type used for MR Fingerprinting. Experimental data were also collected for the final case. Steady-state predictions from EPG-X closely match directly derived steady-state solutions which differ substantially from classic "single pool" EPG predictions. EPG-X(MT) predicts similar MT related levels of signal attenuation in white matter as have been reported elsewhere in the literature. Modelling of CPMG echo trains with EPG-X(BM) suggests that exchange processes can lead to an underestimate of the fraction of short T2 species. Modelling of transient gradient echo sequences with EPG-X(MT) suggests that measurable MT effects result from variable saturation of bound protons, particularly after inversion pulses. In conclusion, EPG-X can be used for modelling of the transient signal response of systems exhibiting chemical exchange or MT. This may be particularly beneficial for relaxometry approaches that rely on characterising transient rather than steady-state sequences.


title: Code Staging in GNU Guix

abstract: GNU Guix is a " functional " package manager that builds upon earlier work on Nix. Guix implements high-level abstractions such as packages and operating system services as domain-specic languages (DSLs) embedded in Scheme. It also implements build actions and operating system orchestration in Scheme. This leads to a multi-tier programming environment where embedded code snippets are staged for eventual execution. This paper presents G-expressions or " gexps ", the staging mechanism we devised for Guix. We explain our journey from traditional Lisp S-expressions to G-expressions, which augment the former with contextual information and ensure hygienic code staging. We discuss the implementation of gexps and report on our experience using them in a variety of operating system use cases-from package build processes to system services. Gexps provide a novel way to cover many aspects of OS connguration in a single, multi-tier language, while facilitating code reuse and code sharing .


title: How does the cosmic web impact assembly bias?

abstract: The mass, accretion rate and formation time of dark matter haloes near proto-filaments (identified as saddle points of the potential) are analytically predicted using a conditional version of the excursion set approach in its so-called "upcrossing" approximation. The model predicts that at fixed mass, mass accretion rate and formation time vary with orientation and distance from the saddle, demonstrating that assembly bias is indeed influenced by the tides imposed by the cosmic web. Starved, early forming haloes of smaller mass lie preferentially along the main axis of filaments, while more massive and younger haloes are found closer to the nodes. Distinct gradients for distinct tracers such as typical mass and accretion rate occur because the saddle condition is anisotropic, and because the statistics of these observables depend on both the conditional means and their covariances. The theory is extended to other critical points of the potential field. The response of the mass function to variations of the matter density field (the so-called large scale bias) is computed, and its trend with accretion rate is shown to invert along the filament. The signature of this model should correspond at low redshift to an excess of reddened galactic hosts at fixed mass along preferred directions, as recently reported in spectroscopic and photometric surveys and in hydrodynamical simulations. The anisotropy of the cosmic web emerges therefore as a significant ingredient to describe jointly the dynamics and physics of galaxies, e.g. in the context of intrinsic alignments or morphological diversity.