title: Phase structure of pure SU(3) lattice gauge theory in 5 dimensions

abstract: We investigate the nonperturbative phase structure of five-dimensional SU(3) pure Yang-Mills theory on the lattice. We perform numerical simulations using the Wilson plaquette gauge action on an anisotropic lattice with a four-dimensional lattice spacing (a4) and with an independent value in the fifth dimension (a5). We investigate both cases of a4 > a5 and a4 < a5. The Polyakov loops in the fourth and the fifth directions are observed, and we find that there are four possible phases for the anisotropic five-dimensional quenched QCD theory on the lattice. We determine the critical values of the lattice bare coupling and the anisotropy parameter for each phase transition. Furthermore, we find that the two center domains where the Polyakov loop have locally different charges of the center symmetry appears in single configuration in the specific region of the lattice parameters.


title: Riesz transform on manifolds with quadratic curvature decay

abstract: We investigate the $L^p$-boundness of the Riesz transform on Riemannian manifolds whose Ricci curvature has quadratic decay. Two criteria for the $L^p$-unboundness of the Riesz transform are given. We recover known results about manifolds that are Euclidean or conical at infinity.


title: Quantum(-like) decision making: on validity of the Aumann theorem

abstract: Through set-theoretic formalization of the notion of common knowledge, Aumann proved that if two agents have the common priors, and their posteriors for a given event are common knowledge, then their posteriors must be equal. In this paper we investigate the problem of validity of this theorem in the framework of quantum(-like) decision making.


title: The young binary HD 102077: Orbit, spectral type, kinematics, and moving group membership

abstract: The K-type binary star HD 102077 was proposed as a candidate member of the TW Hydrae Association (TWA) which is a young (5-15 Myr) moving group in close proximity (~50 pc) to the solar system. The aim of this work is to verify this hypothesis by different means. We first combine diffraction-limited observations from the ESO NTT 3.5m telescope in SDSS-i' and -z' passbands and ESO 3.6m telescope in H-band with literature data to obtain a new, amended orbit fit, estimate the spectral types of both components, and reanalyse the Hipparcos parallax and proper motion taking the orbital motion into account. Moreover, we use two high-resolution spectra of HD 102077 obtained with the fibre-fed optical echelle spectrograph FEROS at the MPG/ESO 2.2m telescope to determine the radial velocity and the lithium equivalent width of the system. The trajectory of HD 102077 is well constrained and we derive a total system mass of $2.6 \pm 0.8\,$ M$_{\odot}$ and a semi-major axis of $14.9 \pm 1.6\,$AU. From the i'-z' colours we infer an integrated spectral type of K2V, and individual spectral types of K0 +/- 1 and K5 +/- 1. The radial velocity corrected for the orbital motion of the system is $17.6 \pm 2\,$km/s. Even though the parallax determination from the Hipparcos data is not influenced by the orbital motion, the proper motion changes to $\mu_\alpha*\cos(\delta) = -137.84 \pm 1.26\,$ mas/yr and $\mu_\delta = -33.53 \pm 1.45 \,$mas/yr. With the resultant space motion, the probability of HD 102077 being a member of TWA is less than 1%. Furthermore, the lithium equivalent width of $200 \pm 4\,$m\AA $\,$ is consistent with an age between 30 Myr and 120 Myr and thus older than the predicted age of TWA. In conclusion, HD 102077's age, galactic space motion, and position do not fit TWA or any other young moving group.


title: Exponential decay properties of a mathematical model for a certain fluid-structure interaction

abstract: In this work, we derive a result of exponential stability for a coupled system of partial differential equations (PDEs) which governs a certain fluid-structure interaction. In particular, a three-dimensional Stokes flow interacts across a boundary interface with a two-dimensional mechanical plate equation. In the case that the PDE plate component is rotational inertia-free, one will have that solutions of this fluid-structure PDE system exhibit an exponential rate of decay. By way of proving this decay, an estimate is obtained for the resolvent of the associated semigroup generator, an estimate which is uniform for frequency domain values along the imaginary axis. Subsequently, we proceed to discuss relevant point control and boundary control scenarios for this fluid-structure PDE model, with an ultimate view to optimal control studies on both finite and infinite horizon. (Because of said exponential stability result, optimal control of the PDE on time interval $(0,\infty)$ becomes a reasonable problem for contemplation.)


title: Estimating the aggregate/intraaggregate mass ratio of a shrinking soil

abstract: A recently introduced parameter, the ratio of an aggregate solid mass to a solid mass of an intraaggregate matrix (K ratio) is connected with the mean thickness of a deformable, but non-shrinking surface layer of aggregates and is a fundamental property of aggregated soils that essentially influences their shrinkage. The objective of this work is to suggest and validate an approach to estimating the K ratio at any soil clay content through characteristics of soil texture and structure. We derive an equation that reflects the interrelation between the K ratio and soil texture and structure. The K ratio can be estimated as the solution of the equation and is determined by the mean size of soil solids and the maximum size of soil aggregates in the oven-dried state, independently of a measured shrinkage curve. To validate the approach we use available data for eight soils.


title: Critical temperature of site-diluted spin-1/2 systems with long-range ferromagnetic interactions

abstract: In the paper the Pair Approximation (PA) method for studies of the site-diluted spin-1/2 systems of arbitrary dimensionality with the long-range ferromagnetic interactions is adopted. The method allows to take into account arbitrary anisotropy of the interactions in the spin space, so it is not limited to purely Ising couplings. Within this approach, the Gibbs free energy is obtained, which allows to derive all the further interesting thermodynamic properties. In particular, we obtain an equation for the critical temperature of the second-order phase transitions for the model in question. In the study we focus our attention on the systems with ferromagnetic interactions decaying with the distance according to the power law $J(r)\propto r^{-n}$. We discuss the dependence of the critical temperature on the concentration of magnetic component and the index $n$ for selected one-, two- and three-dimensional lattices.We confirm the absence of the critical concentration for a diluted magnet with infinite interaction range. In the regime of the low concentrations of magnetic component, we find a non-linear increase of the critical temperature with the concentration in the form of $T_{c}\propto p^{n/d}$, depending on the system dimensionality $d$ and the index $n$, whereas $n > d$.


title: Directional superradiant emission from statistically independent incoherent non-classical and classical sources

abstract: Superradiance is one of the outstanding problems in quantum optics since Dicke introduced the concept of enhanced directional spontaneous emission by an ensemble of identical two-level atoms. The effect is based on correlated collective Dicke states which turn out to be highly entangled. Here we show that enhanced directional emission of spontaneous radiation can be produced also with statistically independent incoherent sources via the measurement of higher order correlation functions of the emitted radiation. Our analysis is applicable to a wide variety of quantum systems like trapped atoms, ions, quantum dots or NV-centers, and is also valid for statistically independent incoherent classical emitters. This is experimentally confirmed with up to eight independent thermal light sources.


title: Direct comparison of optical lattice clocks with an intercontinental baseline of 9 000 km

abstract: We have demonstrated a direct frequency comparison between two $^{87}{\rm Sr}$ lattice clocks operated in intercontinentally separated laboratories in real time. Two-way satellite time and frequency transfer technique based on the carrier phase was employed for a direct comparison with a baseline of 9 000 km between Japan and Germany. A clock comparison was achieved for 83 640 s resulting in a fractional difference of $(1.1\pm1.6) \times 10^{-15}$, where the statistical part is the biggest contribution to the uncertainty. This measurement directly confirms the agreement of the two optical clocks on an intercontinental scale.


title: Ultrahigh energy particle collisions near the black hole horizon in the strong magnetic field

abstract: We consider collision between two charged (or charged and neutral) particles near the black hole horizon in the strong magnetic field B. It is shown that there exists a strip near the horizon within which collision of any two such particles leads to ultrahigh energy in the centre of mass frame. The results apply to generic (not necessarily vacuum) black holes.


title: Controlling the dynamics of a bidimensional gel above and below its percolation transition

abstract: The morphology and the microscopic internal dynamics of a bidimensional gel formed by spontaneous aggregation of gold nanoparticles confined at the water surface are investigated by a suite of techniques, including grazing-incidence x-ray photon correlation spectroscopy (GI-XPCS). The range of concentrations studied spans across the percolation transition for the formation of the gel. The dynamical features observed by GI-XPCS are interpreted in view of the results of microscopical imaging; an intrinsic link between the mechanical modulus and internal dynamics is demonstrated for all the concentrations. Our work presents, to the best of our knowledge, the first example of a transition from stretched to compressed correlation function actively controlled by quasistatically varying the relevant thermodynamic variable. Moreover, by applying a model proposed time ago by Duri and Cipelletti [A. Duri and L. Cipelletti, Europhys. Lett. 76, 972 (2006)] we are able to build a novel master curve for the shape parameter, whose scaling factor allows us to quantify a 'long time displacement length'. This characteristic length is shown to converge, as the concentration is increased, to the 'short time localization length' determined by pseudo Debye-Waller analysis of the initial contrast. Finally, the intrinsic dynamics of the system are then compared with that induced by means of a delicate mechanical perturbation applied to the interface.


title: Sigma Partitioning: Complexity and Random Graphs

abstract: A $\textit{sigma partitioning}$ of a graph $G$ is a partition of the vertices into sets $P_1, \ldots, P_k$ such that for every two adjacent vertices $u$ and $v$ there is an index $i$ such that $u$ and $v$ have different numbers of neighbors in $P_i$. The $\textit{ sigma number}$ of a graph $G$, denoted by $\sigma(G)$, is the minimum number $k$ such that $ G $ has a sigma partitioning $P_1, \ldots, P_k$. Also, a $\textit{ lucky labeling}$ of a graph $G$ is a function $ \ell :V(G) \rightarrow \mathbb{N}$, such that for every two adjacent vertices $ v $ and $ u$ of $ G $, $ \sum_{w \sim v}\ell(w)\neq \sum_{w \sim u}\ell(w) $ ($ x \sim y $ means that $ x $ and $y$ are adjacent). The $\textit{ lucky number}$ of $ G $, denoted by $\eta(G)$, is the minimum number $k $ such that $ G $ has a lucky labeling $ \ell :V(G) \rightarrow \mathbb{N}_k$. It was conjectured in [Inform. Process. Lett., 112(4):109--112, 2012] that it is $ \mathbf{NP} $-complete to decide whether $ \eta(G)=2$ for a given 3-regular graph $G$. In this work, we prove this conjecture. Among other results, we give an upper bound of five for the sigma number of a uniformly random graph.


title: Distribution of electrons and holes in cuprate superconductors as determined from $^{17}$O and $^{63}$Cu nuclear magnetic resonance

abstract: The distribution of electrons and holes in the CuO$_2$ plane of the high-temperature superconducting cuprates is determined with nuclear magnetic resonance through the quadrupole splittings of $^{17}$O and $^{63}$Cu. Based on new data for single crystals of electron-doped Pr$_{2-x}$Ce$_x$CuO$_4$(x=0, 0.05, 0.10, 0.15) as well as Nd$_{2-x}$Ce$_x$CuO$_4$ (x=0, 0.13) the changes in hole contents $n_d$ of Cu 3d$(x^2-y^2)$ and $n_p$ of O 2$p_\sigma$ orbitals are determined and they account for the stoichiometrically doped charges, similar to hole-doped \lsco. It emerges that while $n_d+2n_p=1$ in all parent materials as expected, $n_d$ and $n_p$ vary substantially between different groups of materials. Doping holes increases predominantly $n_p$, but also $n_d$. To the contrary, doping electrons predominantly decreases $n_d$ and only slightly $n_p$. However, $n_p$ for the electron doped systems is higher than that in hole doped La$_{1.85}$Sr$_{0.15}$CuO$_4$. Cuprates with the highest maximum $T_{\rm c}$s appear to have a comparably low $n_d$ while, at the same time, $n_p$ is very high. The rather high oxygen hole content of the Pr$_2$CuO$_4$ and Nd$_2$CuO$_4$ with the low $n_d$ seems to make them ideal candidates for hole doping to obtain the highest $T_{\rm c}$.


title: Spectral Sparse Representation for Clustering: Evolved from PCA, K-means, Laplacian Eigenmap, and Ratio Cut

abstract: Dimensionality reduction, cluster analysis, and sparse representation are basic components in machine learning. However, their relationships have not yet been fully investigated. In this paper, we find that the spectral graph theory underlies a series of these elementary methods and can unify them into a complete framework. The methods include PCA, K-means, Laplacian eigenmap (LE), ratio cut (Rcut), and a new sparse representation method developed by us, called spectral sparse representation (SSR). Further, extended relations to conventional over-complete sparse representations (e.g., method of optimal directions, KSVD), manifold learning (e.g., kernel PCA, multidimensional scaling, Isomap, locally linear embedding), and subspace clustering (e.g., sparse subspace clustering, low-rank representation) are incorporated. We show that, under an ideal condition from the spectral graph theory, PCA, K-means, LE, and Rcut are unified together. And when the condition is relaxed, the unification evolves to SSR, which lies in the intermediate between PCA/LE and K-mean/Rcut. An efficient algorithm, NSCrt, is developed to solve the sparse codes of SSR. SSR combines merits of both sides: its sparse codes reduce dimensionality of data meanwhile revealing cluster structure. For its inherent relation to cluster analysis, the codes of SSR can be directly used for clustering. Scut, a clustering approach derived from SSR reaches the state-of-the-art performance in the spectral clustering family. The one-shot solution obtained by Scut is comparable to the optimal result of K-means that are run many times. Experiments on various data sets demonstrate the properties and strengths of SSR, NSCrt, and Scut.


title: Brackets with $(\tau,\sigma)$-derivations and $(p,q)$-deformations of Witt and Virasoro algebras

abstract: The aim of this paper is to study some brackets defined on $(\tau,\sigma)$-derivations satisfying quasi-Lie identities. Moreover, we provide examples of $(p,q)$-deformations of Witt and Virasoro algebras as well as $\mathfrak{sl}(2)$ algebra. These constructions generalize the results obtained by Hartwig, Larsson and Silvestrov on $\sigma$-derivations, arising in connection with discretizations and deformations of algebras of vector fields.


title: On Hardy q-inequalities

abstract: Some q-analysis variants of Hardy type inequalities of the form \int_0^b (x^{\alpha-1} \int_0^x t^{-\alpha} f(t) d_qt)^p d_qx \leq C \int_0^b f^p(t) d_qt with sharp constant C are proved and discussed. A similar result with the Riemann-Liouville operator involved is also proved. Finally, it is pointed out that by using these techniques we can also obtain some new discrete Hardy and Copson type inequalities in the classical case.


title: Search for top quark decays $t\rightarrow qH$ with $H\to\gamma\gamma$ using the ATLAS detector

abstract: A search is performed for flavour-changing neutral currents in the decay of a top quark to an up-type ($c, u$) quark and a Higgs boson, where the Higgs boson decays to two photons. The proton-proton collision data set used corresponds to 4.7 fb$^{-1}$ at $\sqrt{s}$ = 7 \tev\ and 20.3 fb$^{-1}$ at $\sqrt{s}$ = 8 \tev\ collected by the ATLAS experiment at the LHC. Top quark pair events are searched for in which one top quark decays to $qH$ and the other decays to $bW$. Both the hadronic and the leptonic decay modes of the $W$ boson are used. No significant signal is observed and an upper limit is set on the $t\rightarrow qH$ branching ratio of 0.79% at the 95% confidence level. The corresponding limit on the $tqH$ coupling combination $\sqrt{\lambda_{tcH}^{2} + \lambda_{tuH}^{2}}$ is 0.17.


title: On the origin of bursts in blue compact dwarf galaxies: clues from kinematics and stellar populations

abstract: Blue compact dwarf galaxies (BCDs) form stars at, for their sizes, extraordinarily high rates. In this paper, we study what triggers this starburst and what is the fate of the galaxy once its gas fuel is exhausted. We select four BCDs with smooth outer regions, indicating them as possible progenitors of dwarf elliptical galaxies. We have obtained photometric and spectroscopic data with the FORS and ISAAC instruments on the VLT. We analyse their infra-red spectra using a full spectrum fitting technique which yields the kinematics of their stars and ionized gas together with their stellar population characteristics. We find that the_stellar_ velocity to velocity dispersion ratio (v/sigma) of our BCDs is of the order of 1.5, similar to that of dwarf elliptical galaxies. Thus, those objects do not require significant (if any) loss of angular momentum to fade into early type dwarfs. This finding is in discordance with previous studies, which however compared the stellar kinematics of dwarf elliptical galaxies with the gaseous kinematics of star forming dwarfs. The stellar velocity fields of our objects are very disturbed and the star-formation regions are often kinematically decoupled from the rest of the galaxy. These regions can be more or less metal rich with respect to the galactic body, and sometimes they are long lived. These characteristics prevent us from pinpointing a unique trigger of the star formation, even within the same galaxy. Gas impacts, mergers, and in-spiraling gas clumps are all possible star-formation ignitors for our targets.


title: Asymptotic Properties of Minimum S-Divergence Estimator for Discrete Models

abstract: Robust inference based on the minimization of statistical divergences has proved to be a useful alternative to the classical techniques based on maximum likelihood and related methods. Recently Ghosh et al. (2013) proposed a general class of divergence measures, namely the S-Divergence Family and discussed its usefulness in robust parametric estimation through some numerical illustrations. In this present paper, we develop the asymptotic properties of the proposed minimum S-Divergence estimators under discrete models.


title: On Consistent Hypothesis Testing

abstract: This paper explores conditions of existence of different types of consistent tests. New links of these types of consistency are also established. The existence of discernible (strong consistent) tests follows from the existence of pointwise consistent tests. If there are consistent tests, then the set of alternatives can be represented as countable union of nested subsets such that there are uniformly consistent tests for the hypothesis and each of this subset of alternatives. Implementing these results we explore both sufficient conditions and necessary conditions for existence of consistent, pointwise consistent, uniformly consistent and discernible tests for hypothesis testing on a probability measure of independent sample, on a mean measure of Poisson process, on a solution of linear ill-posed problems in Gaussian noise, on a solution of deconvolution problem and for the problem of signal detection in Gaussian white noise. In the last three cases the necessary conditions and sufficient conditions coincide.


title: Construction of the B88 exchange-energy functional in two dimensions

abstract: We construct a generalized-gradient approximation for the exchange-energy density of finite two-dimensional systems. Guided by non-empirical principles, we include the proper small-gradient limit and the proper tail for the exchange-hole potential. The observed performance is superior to that of the two-dimensional local-density approximation, which underlines the usefulness of the approach in practical applications.


title: Analytic and geometric properties of generic Ricci solitons

abstract: The aim of this paper is to prove some classification results for generic shrinking Ricci solitons. In particular, we show that every three dimensional generic shrinking Ricci soliton is given by quotients of either $\mathds{S}^3$, $\erre\times\mathds{S}^2$ or $\erre^3$, under some very weak conditions on the vector field $X$ generating the soliton structure. In doing so we introduce analytical tools that could be useful in other settings; for instance we prove that the Omori-Yau maximum principle holds for the $X$-Laplacian on every generic Ricci soliton, without any assumption on $X$.


title: QCD sum rules on the complex Borel plane

abstract: Borel transformed QCD sum rules conventionally use a real valued parameter (the Borel mass) for specifying the exponential weight over which hadronic spectral functions are averaged. In this paper, it is shown that the Borel mass can be generalized to have complex values and that new classes of sum rules can be derived from the resulting averages over the spectral functions. The real and imaginary parts of these novel sum rules turn out to have damped oscillating kernels and potentially contain a larger amount of information on the hadronic spectrum than the real valued QCD sum rules. As a first practical test, we have formulated the complex Borel sum rules for the phi meson channel and have analyzed them using the maximum entropy method, by which we can extract the most probable spectral function from the sum rules without strong assumptions on its functional form. As a result, it is demonstrated that, compared to earlier studies, the complex valued sum rules allow us to extract the spectral function with a significantly improved resolution and thus to study more detailed structures of the hadronic spectrum than previously possible.


title: From Galois to Hopf Galois: theory and practice

abstract: Hopf Galois theory expands the classical Galois theory by considering the Galois property in terms of the action of the group algebra k[G] on K/k and then replacing it by the action of a Hopf algebra. We review the case of separable extensions where the Hopf Galois property admits a group-theoretical formulation suitable for counting and classifying, and also to perform explicit computations and explicit descriptions of all the ingredients involved in a Hopf Galois structure. At the end we give just a glimpse of how this theory is used in the context of Galois module theory for wildly ramified extensions.


title: Muon anomalous magnetic moment and positron excess at AMS-02 in a gauged horizontal symmetric model

abstract: We studied an extension of the standard model with a fourth generation of fermions to explain the discrepancy in the muon $(g-2)$ and explain the positron excess seen in the AMS-02 experiment. We introduce a gauged $SU(2)_{HV}$ horizontal symmetry between the muon and the 4th generation lepton families. The 4th generation right-handed neutrino is identified as the dark matter with mass $\sim 700$ GeV. The dark matter annihilates only to $(\mu^+ \mu^-)$ and $(\nu^c_\mu ~\nu_\mu)$ states via $SU(2)_{HV}$ gauge boson. The $SU(2)_{HV}$ gauge boson with mass $\sim 1.4$ TeV gives an adequate contribution to the $(g-2)$ of muon and fulfill the experimental constraint from BNL measurement. The higgs production constraints from 4th generation fermions is evaded by extending the higgs sector.


title: Microscopic description of $^7$Li in the $^{7}\text{Li}+^{12}\text{C}$ and $^7\text{Li}+^{28}\text{Si}$ elastic scattering at high energies

abstract: We employ a microscopic continuum-discretized coupled-channels reaction framework (MCDCC) to study the elastic angular distribution of the $^7$Li$=\alpha+t$ nucleus colliding with $^{12}$C and $^{28}$Si targets at $E_{\text{Lab}}$=350 MeV. In this framework, the $^7$Li projectile is described in a microscopic cluster model and impinges on non-composite targets. The diagonal and coupling potentials are constructed from nucleon-target interactions and $^7$Li microscopic wave functions. We obtain a fair description of the experimental data, in the whole angular range studied, when continuum channels are included. The inelastic and breakup angular distributions on the lightest target are also investigated. In addition, we compute $^{7}$Li$+^{12}$C MCDCC elastic cross sections at energies much higher than the Coulomb barrier and we use them as reference calculations to test the validity of multichannel eikonal cross sections.


title: Verification for Timed Automata extended with Unbounded Discrete Data Structures

abstract: We study decidability of verification problems for timed automata extended with unbounded discrete data structures. More detailed, we extend timed automata with a pushdown stack. In this way, we obtain a strong model that may for instance be used to model real-time programs with procedure calls. It is long known that the reachability problem for this model is decidable. The goal of this paper is to identify subclasses of timed pushdown automata for which the language inclusion problem and related problems are decidable.


title: Estimation of Multivariate Location and Covariance using the S -Hellinger Distance

abstract: This paper describes a generalization of the Hellinger distance which we call the S -Hellinger distance; this general family connects the Hellinger distance smoothly with the $L_2$-divergence by a tuning parameter $\alpha$ and is indeed a subfamily of the S -Divergence family of Ghosh et al. (2013 a, b). We use this general divergence in the context of estimating the location and covariances under (continuous) multivariate models and show that the proposed minimum S -Hellinger distance estimator is affine equivariant, asymptotically consistent and have high breakdown point under suitable conditions. We also illustrate its performance through an extensive simulation study which show that the proposed estimators give more robust estimator than the minimum Hellinger distance estimator for the location and correlation parameters under different types of contamination with the contamination proportion being as high as 20%.


title: A Pattern based methodology for evolution management in business process reuse

abstract: Today, there are Process-Aware Information Systems (PAIS) with a set of business process models which vary over time to meet the new requirements. In a competitive environment, the key challenge of enterprises is to reduce the cost and time of process design and application development. For this purpose, research on reuse in business process management have introduced the concept of configurable process models which attempts to manage business process variability, by integrating a set of process variants in a single model. In this context, many research works were interested in creating and elaborating configurable process models. However, this has become insufficient since the configurable process model should itself evolve to add new variations. In turn, this requires a comprehensive support for managing the evolution of configurable process models. In this paper, we present a complete pattern based methodology for managing the evolution of configurable process models in terms of activities, data and resources. Our objective is to propose a process patterns system for guiding designers in modeling and evolving configurable process models. Furthermore, our process patterns system will be used for an automated support so as to manage the evolution of configurable process models.


title: Negative-U properties for substitutional Au in Si

abstract: The isolated substitutional gold impurity in bulk silicon is studied in detail using electronic structure calculations based on density-functional theory. The defect system is found to be a non-spin-polarized negative-U centre, thus providing a simple solution to the long-standing debate over the electron paramagnetic resonance signal for gold in silicon. There is an excellent agreement (within 0.03 eV) between the well-established experimental donor and acceptor levels and the predicted stable charge state transition levels, allowing for the unambiguous assignment of the two experimental levels to the (1+/1-) and (1-/3-) transitions, respectively, in contrast to previously held assumptions about the system.


title: A sharp refinement of a result of Alon, Ben-Shimon and Krivelevich on bipartite graph vertex sequences

abstract: We give a sharp refinement of a result of Alon, Ben-Shimon and Krivelevich. This gives a sufficient condition for a finite sequence of positive integers to be the vertex degree list of both parts of a bipartite graph. The condition depends only on the length of the sequence and its largest and smallest elements.


title: Non-linear Kalman filters for calibration in radio interferometry

abstract: We present a new calibration scheme based on a non-linear version of Kalman filter that aims at estimating the physical terms appearing in the Radio Interferometry Measurement Equation (RIME). We enrich the filter's structure with a tunable data representation model, together with an augmented measurement model for regularization. We show using simulations that it can properly estimate the physical effects appearing in the RIME. We found that this approach is particularly useful in the most extreme cases such as when ionospheric and clock effects are simultaneously present. Combined with the ability to provide prior knowledge on the expected structure of the physical instrumental effects (expected physical state and dynamics), we obtain a fairly cheap algorithm that we believe to be robust, especially in low signal-to-noise regime. Potentially the use of filters and other similar methods can represent an improvement for calibration in radio interferometry, under the condition that the effects corrupting visibilities are understood and analytically stable. Recursive algorithms are particularly well adapted for pre-calibration and sky model estimate in a streaming way. This may be useful for the SKA-type instruments that produce huge amounts of data that have to be calibrated before being averaged.


title: "TNOs are Cool": A survey of the trans-Neptunian region X. Analysis of classical Kuiper belt objects from Herschel and Spitzer observations

abstract: The classical Kuiper belt contains objects both from a low-inclination, presumably primordial, distribution and from a high-inclination dynamically excited population. Based on a sample of classical TNOs with observations at thermal wavelengths we determine radiometric sizes, geometric albedos and thermal beaming factors as well as study sample properties of dynamically hot and cold classicals. Observations near the thermal peak of TNOs using infra-red space telescopes are combined with optical magnitudes using the radiometric technique with near-Earth asteroid thermal model (NEATM). We have determined three-band flux densities from Herschel/PACS observations at 70.0, 100.0 and 160.0 $\mu$m and Spitzer/MIPS at 23.68 and 71.42 $\mu$m when available. We have analysed 18 classical TNOs with previously unpublished data and re-analysed previously published targets with updated data reduction to determine their sizes and geometric albedos as well as beaming factors when data quality allows. We have combined these samples with classical TNOs with radiometric results in the literature for the analysis of sample properties of a total of 44 objects. We find a median geometric albedo for cold classical TNOs of 0.14 and for dynamically hot classical TNOs, excluding the Haumea family and dwarf planets, 0.085. We have determined the bulk densities of Borasisi-Pabu (2.1 g/cm^3), Varda-Ilmare (1.25 g/cm^3) and 2001 QC298 (1.14 g/cm^3) as well as updated previous density estimates of four targets. We have determined the slope parameter of the debiased cumulative size distribution of dynamically hot classical TNOs as q=2.3 +- 0.1 in the diameter range 100<D<500 km. For dynamically cold classical TNOs we determine q=5.1 +- 1.1 in the diameter range 160<D<280 km as the cold classical TNOs have a smaller maximum size.


title: How much can we learn about the physics of inflation?

abstract: The recent BICEP2 measurement of B-modes in the polarization of the cosmic microwave background suggests that inflation was driven by a field at an energy scale of $2\times 10^{16}$ GeV. I explore the potential of upcoming CMB polarization experiments to further constrain the physics underlying inflation. If the signal is confirmed, then two sets of experiments covering larger area will shed light on inflation. Low resolution measurements can pin down the tensor to scalar ratio at the percent level, thereby distinguishing models from one another. A high angular resolution experiment will be necessary to measure the tilt of the tensor spectrum, testing the consistency relation that relates the tilt to the amplitude.


title: Strangeness enhancement at the hadronic chemical freeze-out

abstract: The chemical freeze-out of hadrons created in the high energy nuclear collisions is studied within the realistic version of the hadron resonance gas model. The chemical non-equilibrium of strange particles is accounted via the usual $\gamma_{s}$ factor which gives us an opportunity to perform a high quality fit with $\chi^2/dof \simeq 63.5/55 \simeq 1.15$ of the hadronic multiplicity ratios measured from the low AGS to the highest RHIC energies. In contrast to previous findings, at low energies we observe the strangeness enhancement instead of a suppression. In addition, the performed $\gamma_{s}$ fit allows us to achieve the highest quality of the Strangeness Horn description with $\chi^2/dof=3.3/14$. For the first time the top point of the Strangeness Horn is perfectly reproduced, which makes our theoretical horn as sharp as an experimental one. However, the $\gamma_{s}$ fit approach does not sizably improve the description of the multi-strange baryons and antibaryons. Therefore, an apparent deviation of multi-strange baryons and antibaryons from chemical equilibrium requires further explanation.


title: Dynamical systems and forward-backward algorithms associated with the sum of a convex subdifferential and a monotone cocoercive operator

abstract: In a Hilbert framework, we introduce continuous and discrete dynamical systems which aim at solving inclusions governed by structured monotone operators $A=\partial\Phi+B$, where $\partial\Phi$ is the subdifferential of a convex lower semicontinuous function $\Phi$, and $B$ is a monotone cocoercive operator. We first consider the extension to this setting of the regularized Newton dynamic with two potentials. Then, we revisit some related dynamical systems, namely the semigroup of contractions generated by $A$, and the continuous gradient projection dynamic. By a Lyapunov analysis, we show the convergence properties of the orbits of these systems. The time discretization of these dynamics gives various forward-backward splitting methods (some new) for solving structured monotone inclusions involving non-potential terms. The convergence of these algorithms is obtained under classical step size limitation. Perspectives are given in the field of numerical splitting methods for optimization, and multi-criteria decision processes.


title: Extremality conditions and regularity of solutions to optimal partition problems involving Laplacian eigenvalues

abstract: Let $\Omega\subset \mathbb{R}^N$ be an open bounded domain and $m\in \mathbb{N}$. Given $k_1,\ldots,k_m\in \mathbb{N}$, we consider a wide class of optimal partition problems involving Dirichlet eigenvalues of elliptic operators, of the following form \[ \inf\left\{F(\lambda_{k_1}(\omega_1),\ldots, \lambda_{k_m}(\omega_m)):\ (\omega_1,\ldots, \omega_m)\in \mathcal{P}_m(\Omega)\right\}, \] where $\lambda_{k_i}(\omega_i)$ denotes the $k_i$--th eigenvalue of $(-\Delta,H^1_0(\omega_i))$ counting multiplicities, and $\mathcal{P}_m(\Omega)$ is the set of all open partitions of $\Omega$, namely \[ \mathcal{P}_m(\Omega)=\left\{(\omega_1,\ldots,\omega_m):\ \omega_i\subset \Omega \text{ open},\ \omega_i\cap \omega_j=\emptyset\ \forall i\neq j\right\}. \] While existence of a quasi-open optimal partition $(\omega_1,\ldots, \omega_m)$ follows from a general result by Bucur, Buttazzo and Henrot [Adv. Math. Sci. Appl. 8, 1998], the aim of this paper is to associate with such minimal partitions and their eigenfunctions some suitable extremality conditions and to exploit them, proving as well the Lipschitz continuity of some eigenfunctions, and regularity of the partition in the sense that the free boundary $\cup_{i=1}^m \partial \omega_i\cap \Omega$ is, up to a residual set, locally a $C^{1,\alpha}$ hypersurface. This last result extend the ones in the paper by Caffarelli and Lin [J. Sci. Comput. 31, 2007] to the case of higher eigenvalues.


title: Large-Scale Distribution of Arrival Directions of Cosmic Rays Detected at the Pierre Auger Observatory Above $10~$PeV

abstract: Searches for large-scale anisotropies in the distribution of arrival directions of cosmic rays detected above $\simeq 10 $PeV at the Pierre Auger Observatory are presented. Although no significant deviation from isotropy is revealed at present, some of the measurements suggest that future data will provide hints for large-scale anisotropies over a wide energy range. Those anisotropies would have amplitudes which are too small to be significantly observed within the current statistics. Assuming that the cosmic ray anisotropy is dominated by dipole and quadrupole moments in the EeV-energy range, some consequences of the present upper limits on their amplitudes are presented.


title: Cost Effective Rumor Containment in Social Networks

abstract: The spread of rumors through social media and online social networks can not only disrupt the daily lives of citizens but also result in loss of life and property. A rumor spreads when individuals, who are unable decide the authenticity of the information, mistake the rumor as genuine information and pass it on to their acquaintances. We propose a solution where a set of individuals (based on their degree) in the social network are trained and provided resources to help them distinguish a rumor from genuine information. By formulating an optimization problem we calculate the optimum set of individuals, who must undergo training, and the quality of training that minimizes the expected training cost and ensures an upper bound on the size of the rumor outbreak. Our primary contribution is that although the optimization problem turns out to be non convex, we show that the problem is equivalent to solving a set of linear programs. This result also allows us to solve the problem of minimizing the size of rumor outbreak for a given cost budget. The optimum solution displays an interesting pattern which can be implemented as a heuristic. These results can prove to be very useful for social planners and law enforcement agencies for preventing dangerous rumors and misinformation epidemics.


title: On the quantisation of gravity by embedding spacetime in a higher dimensional space

abstract: Certain difficulties of quantum gravity can be avoided if we embed the spacetime $V_4$ into a higher dimensional space $V_N$; then our spacetime is merely a 4-surface in $V_N$.What remains is conceptually not so difficult: just to quantise this 4-surface. Our formal procedure generalises our version of Stueckelberg's proper time method of worldline quantisation. We write the equations of $V_4$ in the covariant canonical form starting from a model Lagrangian which contains the classical Einstein gravity as a particular case. Then we perform quantisation in the Schr\"odinger picture by using the concepts of a phase functional and wave functional. As a result we obtain the uncertainty relations which imply that an observer is `aware' either of a particular spacetime surface and has no information about other spacetime surfaces (which represent alternative histories); or conversely, he loses information about a particular $V_4$ whilst he obtains some information about other spacetimes (and histories). Equivalently, one cannot measure to an arbitrary precision both the metric on $V_4$ and matter distribution on various alternative spacetime surfaces. We show how this special case in the `coordinate' representations can be generalised to an arbitrary vector in an abstract Hilbert space.


title: Comparison of spectral invariants in Lagrangian and Hamiltonian Floer theory

abstract: We compare spectral invariants in periodical orbits and Lagrangian Floer homology case, for closed symplectic manifold $P$ and its closed Lagrangian submanifolds $L$, when $\omega|_{\pi_2(P,L)}=0$, and $\mu|_{\pi_2(P,L)}=0$. From this result, we derive a corollary considering comparison of Hofer's distance in periodic orbits and Lagrangian case. We also define a product $HF_*(H)\otimes HF_*(L,\phi^1_H(L))\to HF_*(L,\phi^1_H(L))$ and prove subadditivity of invariants with respect to this product.


title: Stabilizing dual-energy X-ray computed tomography reconstructions using patch-based regularization

abstract: Recent years have seen growing interest in exploiting dual- and multi-energy measurements in computed tomography (CT) in order to characterize material properties as well as object shape. Material characterization is performed by decomposing the scene into constitutive basis functions, such as Compton scatter and photoelectric absorption functions. While well motivated physically, the joint recovery of the spatial distribution of photoelectric and Compton properties is severely complicated by the fact that the data are several orders of magnitude more sensitive to Compton scatter coefficients than to photoelectric absorption, so small errors in Compton estimates can create large artifacts in the photoelectric estimate. To address these issues, we propose a model-based iterative approach which uses patch-based regularization terms to stabilize inversion of photoelectric coefficients, and solve the resulting problem though use of computationally attractive Alternating Direction Method of Multipliers (ADMM) solution techniques. Using simulations and experimental data acquired on a commercial scanner, we demonstrate that the proposed processing can lead to more stable material property estimates which should aid materials characterization in future dual- and multi-energy CT systems.


title: Broad 1st Order Phase Transitions: phase coexistence, unusual hysteresis and glasslike kinetic arrest

abstract: We discuss the phenomenology of phase transitions studied with two control variables. Such measurements have become routine with magnetic field and temperature being varied for 1st order magnetic transitions with an ease not conceivable with pressure and temperature. Similar ease may be possible with electric field as the second control variable, as in current studies on multiferroics. We develop this phenomenology for broad 1st order transitions that occur over a range of the control variable, as is common in substitutional alloys with inherent quenched disorder, including the possibility of a glasslike arrest of the transition occurring within this range of either control variable. The interrupted transition then results in the two competing phases coexisting to the lowest temperature. The new experimental protocol with acronym CHUF, where the sample is cooled and heated under different values of the second control variable, allows a detailed study of glasslike arrest of kinetics of the phase transformation. The many experimental studies on 1st order magnetic transitions in a large number of materials, that provide the cause for this phenomenology, shall be discussed subsequently.


title: Twisting and Bending Stress in DNA Minicircles

abstract: The interplay between bending of the molecule axis and appearance of disruptions in circular DNA molecules, with $\sim 100$ base pairs, is addressed. Three minicircles with different radii and almost equal content of AT and GC pairs are investigated. The DNA sequences are modeled by a mesoscopic Hamiltonian which describes the essential interactions in the helix at the level of the base pair and incorporates twisting and bending degrees of freedom. Helix unwinding and bubble formation patterns are consistently computed by a path integral method that sums over a large number of molecule configurations compatible with the model potential. The path ensembles are determined, as a function of temperature, by minimizing the free energy of the system. Fluctuational openings appear along the helix to release the stress due to the bending of the molecule backbone. In agreement with the experimental findings, base pair disruptions are found with larger probability in the smallest minicircle of \textit{66-bps} whose bending angle is $\sim 6^{o} $. For this minicircle, a sizeable untwisting is obtained with the helical repeat showing a step-like increase at $\tau =\,315K$. The method can be generalized to determine the bubble probability profiles of open ends linear sequences.

categories: cond-mat.soft cond-mat.stat-mech q-bio.BM q-bio.QM


title: Evolution of transverse flow and effective temperatures in the parton phase from a multi-phase transport model

abstract: I study the space-time evolution of transverse flow and effective temperatures in the dense parton phase with the string melting version of a multi-phase transport model. Parameters of the model are first constrained to reproduce the bulk data on the rapidity density, $p_{\rm T}$ spectrum and elliptic flow at low $p_{\rm T}$ for central and mid-central Au+Au collisions at $200A$ GeV and Pb+Pb collisions at $2760A$ GeV. I then calculate the transverse flow and effective temperatures in volume cells within mid-spacetime-rapidity $|\eta|<1/2$. I find that the effective temperatures extracted from different variables, which are all evaluated in the rest frame of a volume cell, can be very different; this indicates that the parton system in the model is not in full chemical or thermal equilibrium locally, even after averaging over many events. In particular, the effective temperatures extracted from the parton energy density or number density are often quite different than those extracted from the parton mean $p_{\rm T}$ or mean energy. For these collisions in general, effective temperatures extracted from the parton energy density or number density are higher than those extracted from the parton mean $p_{\rm T}$ in the inner part of the overlap volume, while the opposite occurs in the outer part of the overlap volume. I argue that this indicates that the dense parton matter in the inner part of the overlap volume is over-populated; I also find that all cells with energy density above 1 GeV/fm$^3$ are over-populated after a couple of fm/$c$.


title: Do the Fix Ingredients Already Exist? An Empirical Inquiry into the Redundancy Assumptions of Program Repair Approaches

abstract: Much initial research on automatic program repair has focused on experimental results to probe their potential to find patches and reduce development effort. Relatively less effort has been put into understanding the hows and whys of such approaches. For example, a critical assumption of the GenProg technique is that certain bugs can be fixed by copying and re-arranging existing code. In other words, GenProg assumes that the fix ingredients already exist elsewhere in the code. In this paper, we formalize these assumptions around the concept of ''temporal redundancy''. A temporally redundant commit is only composed of what has already existed in previous commits. Our experiments show that a large proportion of commits that add existing code are temporally redundant. This validates the fundamental redundancy assumption of GenProg.


title: On Progressive Filtration Expansions with a Process; Applications to Insider Trading

abstract: In this paper we study progressive filtration expansions with c\`adl\`ag processes. Using results from the theory of the weak convergence of $\sigma$-fields, we first establish a semimartingale convergence theorem. Then we apply it in a filtration expansion with a process setting and provide sufficient conditions for a semimartingale of the base filtration to remain a semimartingale in the expanded filtration. Applications to the expansion of a Brownian filtration are given. The paper concludes with applications to models of insider trading in financial mathematics.


title: A characterization of sub-game perfect Nash equilibria for SDEs of mean field type

abstract: We study a class of dynamic decision problems of mean field type with time inconsistent cost functionals, and derive a stochastic maximum principle to characterize subgame perfect Nash equilibrium points. Subsequently, this approach is extended to a mean field game to construct decentralized strategies and obtain an estimate of their performance.


title: The Clusters AgeS Experiment (CASE). VI. Analysis of two detached eclipsing binaries in the Globular Cluster M55

abstract: We present an analysis of the detached eclipsing binaries V44 and V54 belonging to the globular cluster M55. For V54 we obtain the following absolute parameters: $M_p=0.726\pm 0.015\,M_\odot$, $R_p=1.006\pm 0.009\,R_\odot$, $L_p=1.38\pm 0.07\,L_\odot$ for the primary, and $M_s=0.555\pm 0.008\,M_\odot$, $R_s=0.528\pm 0.005\,R_\odot$, $L_s=0.16\pm0.01\,L_\odot$ for the secondary. The age and apparent distance modulus of V54 are estimated at 13.3 - 14.7 Gyr and $13.94\pm0.05$ mag, respectively. This age is substantially larger than ages we have derived from the analysis of binary systems in 47 Tuc and M4. The secondary of V44 is so weak in the optical domain that only mass function and relative parameters are obtained for the components of this system. However, there is a good chance that the velocity curve of the secondary could be derived from near-IR spectra. As the primary of V44 is more evolved than that of V54, such data would impose much tighter limits on the age and distance of M55.


title: Distributive Gruppen endlicher Ordnung [Finite Distributive Groups]

abstract: This is a translation. I have added translations for (possibly) outdated definitions in an appendix at the end. In this paper, we define distributive groups and show some properties of them. We then concern ourselves with the homogeinity of distributive groups, before showing how they can be generated from both associative and distributive groups. After that, we look at subgroups and define an index theorem for distributive groups before describing the structure of distributive groups. Finally, we present two addenda with several theorems that we proved while trying to prove that $\left|(A.p).(A.q)\right|=\left|A\right|$.


title: The birth of black holes: neutron star collapse times, gamma-ray bursts and fast radio bursts

abstract: Recent observations of short gamma-ray bursts (SGRBs) suggest that binary neutron star (NS) mergers can create highly magnetised, millisecond NSs. Sharp cut-offs in X-ray afterglow plateaus of some SGRBs hint at the gravitational collapse of these remnant NSs to black holes. The collapse of such `supramassive' NSs also describes the blitzar model, a leading candidate for the progenitors of fast radio bursts (FRBs). The observation of an FRB associated with an SGRB would provide compelling evidence for the blitzar model and the binary NS merger scenario of SGRBs, and lead to interesting constraints on the NS equation of state. We predict the collapse times of supramassive NSs created in binary NS mergers, finding that such stars collapse $\sim10\,{\rm s}$ -- $4.4\times10^{4}\,{\rm s}$ (95% confidence) after the merger. This directly impacts observations targeting NS remnants of binary NS mergers, providing the optimal window for high time resolution radio and X-ray follow-up of SGRBs and gravitational wave bursts.


title: Analysis of self-overlap reveals trade-offs in plankton swimming trajectories

abstract: Movement is a fundamental behaviour of organisms that brings about beneficial encounters with resources and mates, but at the same time exposes the organism to dangerous encounters with predators. The movement patterns adopted by organisms should reflect a balance between these contrasting processes. This trade-off can be hypothesized as being evident in the behaviour of plankton, which inhabit a dilute 3D environment with few refuges or orienting landmarks. We present an analysis of the swimming path geometries based on a volumetric Monte Carlo sampling approach, which is particularly adept at revealing such trade-offs by measuring the self-overlap of the trajectories. Application of this method to experimentally measured trajectories reveals that swimming patterns in copepods are shaped to efficiently explore volumes at small scales, while achieving a large overlap at larger scales. Regularities in the observed trajectories make the transition between these two regimes always sharper than in randomized trajectories or as predicted by random walk theory. Thus real trajectories present a stronger separation between exploration for food and exposure to predators. The specific scale and features of this transition depend on species, gender, and local environmental conditions, pointing at adaptation to state and stage dependent evolutionary trade-offs.


title: Simpson's Paradox and Collapsibility

abstract: Simpson's paradox and collapsibility are two closely related concepts in the context of data analysis. While the knowledge about the occurrence of Simpson's paradox helps a statistician to draw correct and meaningful conclusions, the concept of collapsibility deals with dimension-reduction aspects, when Simpson's paradox does not occur. We discuss in this paper in some detail the nature and the genesis of Simpson's paradox with respect to well-known examples and also various concepts of collapsiblity. The main aim is to bring out the close connections between these two phenomena, especially with regard to the analysis of contingency tables, regression models and a certain measure of association or a dependence function. There is a vast literature on these topics and so we focus only on certain aspects, recent developments and some important results in the above-mentioned areas.


title: Problem Complexity in Parallel Problem Solving

abstract: Recent works examine the relationship between the communication structure and the performance of a group in a problem solving task. Some conclude that inefficient communication networks with long paths outperform efficient networks on the long run. Others find no influence of the network topology on group performance. We contribute to this discussion by examining the role of problem complexity. In particular, we study whether and how the complexity of the problem at hand moderates the influence of the communication network on group performance. Results obtained from multi-agent modelling suggest that problem complexity indeed has an influence. We observe an influence of the network only for problems of moderate difficulty. For easier or harder problems, the influence of network topology becomes weaker or irrelevant, which offers a possible explanation for inconsistencies in the literature.


title: On the Computational Complexity of Vertex Integrity and Component Order Connectivity

abstract: The Weighted Vertex Integrity (wVI) problem takes as input an $n$-vertex graph $G$, a weight function $w:V(G)\to\mathbb{N}$, and an integer $p$. The task is to decide if there exists a set $X\subseteq V(G)$ such that the weight of $X$ plus the weight of a heaviest component of $G-X$ is at most $p$. Among other results, we prove that: (1) wVI is NP-complete on co-comparability graphs, even if each vertex has weight $1$; (2) wVI can be solved in $O(p^{p+1}n)$ time; (3) wVI admits a kernel with at most $p^3$ vertices. Result (1) refutes a conjecture by Ray and Deogun and answers an open question by Ray et al. It also complements a result by Kratsch et al., stating that the unweighted version of the problem can be solved in polynomial time on co-comparability graphs of bounded dimension, provided that an intersection model of the input graph is given as part of the input. An instance of the Weighted Component Order Connectivity (wCOC) problem consists of an $n$-vertex graph $G$, a weight function $w:V(G)\to \mathbb{N}$, and two integers $k$ and $l$, and the task is to decide if there exists a set $X\subseteq V(G)$ such that the weight of $X$ is at most $k$ and the weight of a heaviest component of $G-X$ is at most $l$. In some sense, the wCOC problem can be seen as a refined version of the wVI problem. We prove, among other results, that: (4) wCOC can be solved in $O(\min\{k,l\}\cdot n^3)$ time on interval graphs, while the unweighted version can be solved in $O(n^2)$ time on this graph class; (5) wCOC is W[1]-hard on split graphs when parameterized by $k$ or by $l$; (6) wCOC can be solved in $2^{O(k\log l)} n$ time; (7) wCOC admits a kernel with at most $kl(k+l)+k$ vertices. We also show that result (6) is essentially tight by proving that wCOC cannot be solved in $2^{o(k \log l)}n^{O(1)}$ time, unless the ETH fails.


title: Variable speed branching Brownian motion 1. Extremal processes in the weak correlation regime

abstract: We prove the convergence of the extremal processes for variable speed branching Brownian motions where the "speed functions", that describe the time-inhomogeneous variance, lie strictly below their concave hull and satisfy a certain weak regularity condition. These limiting objects are universal in the sense that they only depend on the slope of the speed function at $0$ and the final time $t$. The proof is based on previous results for two-speed BBM obtained in a recent paper of ours and uses Gaussian comparison arguments to extend these to the general case.


title: Range Expansion of Heterogeneous Populations

abstract: Risk spreading in bacterial populations is generally regarded as a strategy to maximize survival. Here, we study its role during range expansion of a genetically diverse population where growth and motility are two alternative traits. We find that during the initial expansion phase fast growing cells do have a selective advantage. By contrast, asymptotically, generalists balancing motility and reproduction are evolutionarily most successful. These findings are rationalized by a set of coupled Fisher equations complemented by stochastic simulations.


title: A Tight Bound on the Projective Dimension of Four Quadrics

abstract: Motivated by Stillman's question, we show that the projective dimension of an ideal generated by four quadric forms in a polynomial ring is at most 6; moreover, this bound is tight. We achieve this bound, in part, by giving a characterization of the low degree generators of ideals primary to height three primes of multiplicities one and two.


title: Herschel/PACS observations of the 69 $\mu m$ band of crystalline olivine around evolved stars

abstract: We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $\mu$m. This wavelength range covers the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{2-2x}\text{Fe}_{(2x)}\text{SiO}_{4}$). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, $\dot M \ge 10^{-5}$ M$_\odot$/yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. The goal of this study is to exploit the spectral properties of the 69 $\mu$m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. We fit the 69 $\mu$m band and use its width and position to probe the composition and temperature of the crystalline olivine. For 27 sources in the sample, we detected the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $\mu$m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 $\mu$m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 $\mu$m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. continued...


title: Magnetic instability in AdS/CFT : Schwinger effect and Euler-Heisenberg Lagrangian of Supersymmetric QCD

abstract: To reveal the Schwinger effect for quarks, i.e., pair creation process of quarks and antiquarks, we derive the vacuum decay rate at strong coupling using AdS/CFT correspondence. Magnetic fields, in addition to the electric field responsible for the pair creation, causes prominent effects on the rate, and is important also in experiments such as RHIC/LHC heavy ion collisions. In this paper, through the gravity dual we obtain the full Euler-Heisenberg Lagrangian of N=2 supersymmetric QCD and study the Schwinger mechanism with not only a constant electric field but also a constant magnetic field as external fields. We determine the quark mass and temperature dependence of the Lagrangian. In sharp contrast with the zero magnetic field case, we find that the imaginary part, and thus the vacuum decay rate, diverges in the massless zero-temperature limit. This may be related to a strong instability of the QCD vacuum in strong magnetic fields. The real part of the Lagrangian serves as a generating function for non-linear electro-magnetic responses, and is found such that the Cotton-Mouton effect vanishes. Interestingly, our results of the Schwinger / Cotton-Mouton effects coincide precisely with those of N=2 supersymmetric QED.


title: A derivative-free optimization method to compute scalar perturbation of AdS black holes

abstract: In this work we describe an interesting application of a simple derivative-free optimization method to extract the quasinormal modes (QNM's) of a massive scalar field propagating in a 4-dimensional Schwarzschild anti-de Sitter black hole (Sch-AdS$_4$). In this approach, the problem to find the QNM's is reduced to minimize a real valued function of two variables and does not require any information about derivatives. In fact, our strategy requires only evaluations of the objective function to search global minimizers of the optimization problem. Firstly, numerical experiments were performed to find the QNM's of a massless scalar field propagating in intermediate and large Sch-AdS$_4$ black holes. The performance of this optimization algorithm was compared with other numerical methods used in previous works. Our results showed to be in good agreement with those obtained previously. Finally, the massive scalar field case and its QNM's were also obtained and discussed.


title: Multi-wavelength study of the MGRO J2019+37

abstract: MGRO J2019+37 within the Cygnus region is a bright and extended source revealed by Milagro at 12-35 TeV. This source is almost as bright as Crab Nebula in northern sky, while it is not confirmed by ARGO-YBJ around TeV. Up to now, no obvious counterpart at low energy wavelengths has been found. Hence, MGRO J2019+37 becomes mysterious and its VHE gamma-ray emission mechanism is attractive. In this paper, a brief summary of the multi-wavelength observations from Radio to gamma-ray is presented. All the available data from the XMM-Newton and INTEGRAL at X-ray, and the Fermi-LAT at gamma-ray bands were used to get constraint on its emission flux at low energy wavelengths. Then, its possible counterparts and the VHE emission mechanism are discussed.


title: Brain modularity controls the critical behavior of spontaneous activity

abstract: The human brain exhibits a complex structure made of scale-free highly connected modules loosely interconnected by weaker links to form a small-world network. These features appear in healthy patients whereas neurological diseases often modify this structure. An important open question concerns the role of brain modularity in sustaining the critical behaviour of spontaneous activity. Here we analyse the neuronal activity of a model, successful in reproducing on non-modular networks the scaling behaviour observed in experimental data, on a modular network implementing the main statistical features measured in human brain. We show that on a modular network, regardless the strength of the synaptic connections or the modular size and number, activity is never fully scale-free. Neuronal avalanches can invade different modules which results in an activity depression, hindering further avalanche propagation. Critical behaviour is solely recovered if inter-module connections are added, modifying the modular into a more random structure.


title: Feasibility of single-photon cross-phase modulation using metastable xenon in a high finesse cavity

abstract: Cross-phase modulation at the single-photon level has a wide variety of fundamental applications in quantum optics including the generation of macroscopic entangled states. Here we describe a practical method for producing a weak cross-phase modulation at the single-photon level using metastable xenon in a high finesse cavity. We estimate the achievable phase shift and give a brief update on the experimental progress towards its realization. A single-photon cross-phase modulation of approximately 20 milliradians is predicted by both a straightforward perturbation theory calculation and a numerical matrix diagonalization method.


title: A Note on the Sagnac Effect and Current Terrestrial Experiments

abstract: We focus on the Sagnac effect for light beams in order to evaluate if the higher order relativistic corrections of kinematic origin could be relevant for actual terrestrial experiments. Moreover, we discuss to what extent the analogy with the Aharonov-Bohm effect holds true in a fully relativistic framework. We show that the analogy with the Aharonov-Bohm is not true in general, but is recovered in a suitable low order approximation, and that even though the Sagnac effect is influenced by both the position of the interferometer in the rotating frame and its extension, these effects are negligible for current terrestrial experiments.


title: Behavioral and Network Origins of Wealth Inequality: Insights from a Virtual World

abstract: Almost universally, wealth is not distributed uniformly within societies or economies. Even though wealth data have been collected in various forms for centuries, the origins for the observed wealth-disparity and social inequality are not yet fully understood. Especially the impact and connections of human behavior on wealth could so far not be inferred from data. Here we study wealth data from the virtual economy of the massive multiplayer online game (MMOG) Pardus. This data not only contains every player's wealth at every point in time, but also all actions of every player over a timespan of almost a decade. We find that wealth distributions in the virtual world are very similar to those in western countries. In particular we find an approximate exponential for low wealth and a power-law tail. The Gini index is found to be $g=0.65$, which is close to the indices of many Western countries. We find that wealth-increase rates depend on the time when players entered the game. Players that entered the game early on tend to have remarkably higher wealth-increase rates than those who joined later. Studying the players' positions within their social networks, we find that the local position in the trade network is most relevant for wealth. Wealthy people have high in- and out-degree in the trade network, relatively low nearest-neighbor degree and a low clustering coefficient. Wealthy players have many mutual friendships and are socially well respected by others, but spend more time on business than on socializing. We find that players that are not organized within social groups with at least three members are significantly poorer on average. We observe that high `political' status and high wealth go hand in hand. Wealthy players have few personal enemies, but show animosity towards players that behave as public enemies.


title: Equilibrium interfacial free energies and Turnbull coefficient for bcc crystallizing colloidal charged sphere suspensions

abstract: We extend previous analysis of data for the melt-nucleus interfacial free energy, $\gamma$, gained from optical experiments on suspensions of charged colloidal spheres, which crystallize with body centred cubic (bcc) crystal structures. Compiling data from five pure species with different polydispersities and one binary mixture, we find the equilibrium melt-crystal interfacial energy to be considerably larger than the hard sphere reference value. Both this quantity and the entropy of freezing decrease with increasing polydispersity. Moreover, we give a first experimental determination of the Turnbull coefficient for a bcc crystallizing material. The observed value $C_{T, bcc} \approx 0.3$ agrees well with theoretical expectations for bcc systems with short to medium ranged interactions.


title: Neutrino electromagnetic interactions: a window to new physics

abstract: We review the theory and phenomenology of neutrino electromagnetic interactions, which give us powerful tools to probe the physics beyond the Standard Model. After a derivation of the general structure of the electromagnetic interactions of Dirac and Majorana neutrinos in the one-photon approximation, we discuss the effects of neutrino electromagnetic interactions in terrestrial experiments and in astrophysical environments. We present the experimental bounds on neutrino electromagnetic properties and we confront them with the predictions of theories beyond the Standard Model.


title: On variants of $H$-measures and compensated compactness

abstract: We introduce new variant of $H$-measures defined on spectra of general algebra of test symbols and derive the localization properties of such $H$-measures. Applications for the compensated compactness theory are given. In particular, we present new compensated compactness results for quadratic functionals in the case of general pseudo-differential constraints. The case of inhomogeneous second order differential constraints is also studied.


title: How to find metal-rich asteroids

abstract: The metal content of asteroids is of great interest, not only for theories of their origins and the evolution of the solar system but, in the case of near-Earth objects (NEOs), also for impact mitigation planning and endeavors in the field of planetary resources. However, since the reflection spectra of metallic asteroids are largely featureless, it is difficult to identify them and relatively few are known. We show how data from the Wide-field Infrared Survey Explorer (WISE)/NEOWISE thermal-infrared survey and similar surveys, fitted with a simple thermal model, can reveal objects likely to be metal rich. We provide a list of candidate metal-rich NEOs. Our results imply that future infrared surveys with the appropriate instrumentation could discover many more metal-rich asteroids, providing valuable data for assessment of the impact hazard and the potential of NEOs as reservoirs of vital materials for future interplanetary space activities and, eventually perhaps, for use on Earth.


title: Finding Shortest Paths between Graph Colourings

abstract: The $k$-colouring reconfiguration problem asks whether, for a given graph $G$, two proper $k$-colourings $\alpha$ and $\beta$ of $G$, and a positive integer $\ell$, there exists a sequence of at most $\ell+1$ proper $k$-colourings of $G$ which starts with $\alpha$ and ends with $\beta$ and where successive colourings in the sequence differ on exactly one vertex of $G$. We give a complete picture of the parameterized complexity of the $k$-colouring reconfiguration problem for each fixed $k$ when parameterized by $\ell$. First we show that the $k$-colouring reconfiguration problem is polynomial-time solvable for $k=3$, settling an open problem of Cereceda, van den Heuvel and Johnson. Then, for all $k \geq 4$, we show that the $k$-colouring reconfiguration problem, when parameterized by $\ell$, is fixed-parameter tractable (addressing a question of Mouawad, Nishimura, Raman, Simjour and Suzuki) but that it has no polynomial kernel unless the polynomial hierarchy collapses.


title: Updating Formulas and Algorithms for Computing Entropy and Gini Index from Time-Changing Data Streams

abstract: Despite growing interest in data stream mining the most successful incremental learners, such as VFDT, still use periodic recomputation to update attribute information gains and Gini indices. This note provides simple incremental formulas and algorithms for computing entropy and Gini index from time-changing data streams.


title: Scoping Studies for NBI Launch Geometries on DEMO

abstract: Scans of Neutral Beam Injection (NBI) tangency radii and elevation on two possible DEMO scenarios have been performed with two beam energies, 1.5MeV and 1.0MeV, in order to determine the most favourable options for Neutral Beam Current Drive (NBCD) efficiency. In addition, a method using a genetic algorithm has been used to seek optimised solutions of NBI source locations and powers to synthesize a target total plasma driven-current profile. It is found that certain beam trajectories may be proscribed by limitations on shinethrough onto the vessel wall. This may affect the ability of NBCD to extend the duration of a pulse in a scenario where it must complement the induced plasma current. Operating at the lower beam energy reduces the restrictions due to shinethrough and is attractive for technical reasons, but in the scenarios examined here this results in a spatial broadening of the NBCD profile, which may make it more challenging to achieve desired total driven-current profiles.


title: On Modeling and Measuring the Temperature of the z~5 IGM

abstract: The temperature of the low-density intergalactic medium (IGM) at high redshift is sensitive to the timing and nature of hydrogen and HeII reionization, and can be measured from Lyman-alpha forest absorption spectra. Since the memory of intergalactic gas to heating during reionization gradually fades, measurements as close as possible to reionization are desirable. In addition, measuring the IGM temperature at sufficiently high redshifts should help to isolate the effects of hydrogen reionization since HeII reionization starts later, at lower redshift. Motivated by this, we model the IGM temperature at z>5 using semi-numeric models of patchy reionization. We construct mock Lyman-alpha forest spectra from these models and consider their observable implications. We find that the small-scale structure in the Lyman-alpha forest is sensitive to the temperature of the IGM even at redshifts where the average absorption in the forest is as high as 90%. We forecast the accuracy at which the z~5 IGM temperature can be measured using existing samples of high resolution quasar spectra, and find that interesting constraints are possible. For example, an early reionization model in which reionization ends at z~10 should be distinguishable -- at high statistical significance -- from a lower redshift model where reionization completes at z~6. We discuss improvements to our modeling that may be required to robustly interpret future measurements.


title: Submodularity of Energy Related Controllability Metrics

abstract: The quantification of controllability and observability has recently received new interest in the context of large, complex networks of dynamical systems. A fundamental but computationally difficult problem is the placement or selection of actuators and sensors that optimize real-valued controllability and observability metrics of the network. We show that several classes of energy related metrics associated with the controllability Gramian in linear dynamical systems have a strong structural property, called submodularity. This property allows for an approximation guarantee by using a simple greedy heuristic for their maximization. The results are illustrated for randomly generated systems and for placement of power electronic actuators in a model of the European power grid.


title: Spin waves in the anisotropic fcc kagome antiferromagnet

abstract: Spin wave calculations demonstrate that the macroscopic continuous degeneracy associated with the two-dimensional kagome Heisenberg spin lattice persists in the case of the stacked fcc structure giving rise to zero energy modes in three dimensions. The addition of an effective local cubic anisotropy is shown to remove this continuous degeneracy and introduce a gap in the spectrum as well as modify the inelastic scattering function $S({\bf q},\omega)$. This scenario supports earlier Monte Carlo simulations which indicate that the phase transition to long range $q=0$ magnetic order is driven to be discontinuous by critical fluctuations associated with the large degeneracy in the absence of anisotropy, but becomes continuous with the addition of anisotropy. The results are relevant to Ir-Mn alloys which are widely used in the magnetic storage industry in thin film form as the antiferromagnetic pinning layer in spin-valve structures.


title: Singular stochastic PDEs

abstract: We present a series of recent results on the well-posedness of very singular parabolic stochastic partial differential equations. These equations are such that the question of what it even means to be a solution is highly non-trivial. This problem can be addressed within the framework of the recently developed theory of "regularity structures", which allows to describe candidate solutions locally by a "jet", but where the usual Taylor polynomials are replaced by a sequence of custom-built objects. In order to illustrate the theory, we focus on the particular example of the Kardar-Parisi-Zhang equation, a popular model for interface propagation.


title: Intermittency and Alignment in Strong RMHD Turbulence

abstract: We develop an analytic model of intermittent, three-dimensional, strong, reduced magnetohydrodynamic turbulence with zero cross helicity. We take the fluctuation amplitudes to have a log-Poisson distribution and incorporate into the model a new phenomenology of scale-dependent dynamic alignment. The log-Poisson distribution in our model is characterized by two parameters. To calculate these parameters, we make use of three assumptions: that the energy cascade rate is independent of scale within the inertial range, that the most intense coherent structures at scale $\lambda$ are sheet-like with a volume filling factor proportional to $\lambda$, and that most of the cascade power arises from interactions between exceptionally intense fluctuations and much weaker fluctuations. We then compute the scalings of the power spectrum, the kurtosis, higher-order structure functions, and three different average alignment angles. These scalings appear to be consistent with numerous results from direct numerical simulations.


title: Continuum limit of total variation on point clouds

abstract: We consider point clouds obtained as random samples of a measure on a Euclidean domain. A graph representing the point cloud is obtained by assigning weights to edges based on the distance between the points they connect. Our goal is to develop mathematical tools needed to study the consistency, as the number of available data points increases, of graph-based machine learning algorithms for tasks such as clustering. In particular, we study when is the cut capacity, and more generally total variation, on these graphs a good approximation of the perimeter (total variation) in the continuum setting. We address this question in the setting of $\Gamma$-convergence. We obtain almost optimal conditions on the scaling, as number of points increases, of the size of the neighborhood over which the points are connected by an edge for the $\Gamma$-convergence to hold. Taking the limit is enabled by a transportation based metric which allows to suitably compare functionals defined on different point clouds.


title: The Logarithmic Contributions to the O(\alpha_s^3) Asymptotic Massive Wilson Coefficients and Operator Matrix Elements in Deeply Inelastic Scattering

abstract: We calculate the logarithmic contributions to the massive Wilson coefficients for deep-inelastic scattering in the asymptotic region $Q^2 \gg m^2$ to 3-loop order in the fixed-flavor number scheme and present the corresponding expressions for the massive operator matrix elements needed in the variable flavor number scheme. Explicit expressions are given both in Mellin-$N$ space and $z$-space.


title: Universal Duality in Luttinger Liquid Coupled to Generic Environment

abstract: We study a Luttinger Liquid (LL) coupled to a generic environment consisting of bosonic modes with arbitrary density-density and current-current interactions. The LL can be either in the conducting phase and perturbed by a weak scatterer or in the insulating phase and perturbed by a weak link. The environment modes can also be scattered by the imperfection in the system with arbitrary transmission and reflection amplitudes. We present a general method of calculating correlation functions under the presence of the environment and prove the duality of exponents describing the scaling of the weak scatterer and of the weak link. This duality holds true for a broad class of models and is sensitive to neither interaction nor environmental modes details, thus it shows up as the universal property. It ensures that the environment cannot generate new stable fixed points of the RG flow. Thus, the LL always flows toward either conducting or insulating phase. Phases are separated by a sharp boundary which is shifted by the influence of the environment. Our results are relevant, for example, for low-energy transport in (i) an interacting quantum wire or a carbon nanotube where the electrons are coupled to the acoustic phonons scattered by the lattice defect; (ii) a mixture of interacting fermionic and bosonic cold atoms where the bosonic modes are scattered due to an abrupt local change of the interaction, (iii) mesoscopic electric circuits.


title: Immunophenotypes of Acute Myeloid Leukemia From Flow Cytometry Data Using Templates

abstract: Motivation: We investigate whether a template-based classification pipeline could be used to identify immunophenotypes in (and thereby classify) a heterogeneous disease with many subtypes. The disease we consider here is Acute Myeloid Leukemia, which is heterogeneous at the morphologic, cytogenetic and molecular levels, with several known subtypes. The prognosis and treatment for AML depends on the subtype. Results: We apply flowMatch, an algorithmic pipeline for flow cytometry data created in earlier work, to compute templates succinctly summarizing classes of AML and healthy samples. We develop a scoring function that accounts for features of the AML data such as heterogeneity to identify immunophenotypes corresponding to various AML subtypes, including APL. All of the AML samples in the test set are classified correctly with high confidence. Availability: flowMatch is available at www.bioconductor.org/packages/devel/bioc/html/flowMatch.html; programs specific to immunophenotyping AML are at www.cs.purdue.edu/homes/aazad/software.html.


title: Coexistence of fast photodarkening and slow photobleaching in Ge19As21Se60 thin films

abstract: We experimentally demonstrate the coexistence of two opposite photo-effects, viz. fast photodarkening (PD) and slow photobleaching (PB) in Ge19As21Se60 thin films, when illuminated with a laser of wavelength 671nm, PD appears to begin instantaneously upon light illumination and saturates in tens of seconds. By comparison, PB is a slower process that starts only after PD has saturated. Although we could observe the coexistence of PD/PB even at moderate, one order of magnitude lower intensity of 0.2 W/cm2, the kinetics of transformation is significantly slowed down. However, both PD and PB follow stretched exponetial dependence on time. Modeling of overall change as a linear sum of two contributions suggests that the changes in As and Ge parts of glass network respond to light indepndent of each other.


title: Multipacting saturation in parallel plate and micro-pulse electron gun

abstract: A novel parallel plate model is proposed that divided the electron cloud into three parts at saturation, and it is studied in detail using both an analytical approach and PIC (Particle In Cell) code simulations. As one part of the electron cloud, ribbons modes are suggested by tracking the trajectory of individual particle, and the aim of this mode form is to simplify the progress of multipacting effect in the parallel plate so as to be eliminated by optimizing RF parameters. The micro-pulse electron gun (MPG) has demonstrated the potential to address the need for high peak and average current electron beams, hence studying the multipacting in MPG is essential. On the basis of multipacting studying in the parallel plate, it is clear that increasing the cavity voltage is of interest in yielding high quality beams in the gun.


title: Weak locking capacity of quantum channels can be much larger than private capacity

abstract: We show that it is possible for the so-called weak locking capacity of a quantum channel [Guha et al., PRX 4:011016, 2014] to be much larger than its private capacity. Both reflect different ways of capturing the notion of reliable communication via a quantum system while leaking almost no information to an eavesdropper; the difference is that the latter imposes an intrinsically quantum security criterion whereas the former requires only a weaker, classical condition. The channels for which this separation is most straightforward to establish are the complementary channels of classical-quantum (cq-)channels, and hence a subclass of Hadamard channels. We also prove that certain symmetric channels (related to photon number splitting) have positive weak locking capacity in the presence of a vanishingly small pre-shared secret, whereas their private capacity is zero. These findings are powerful illustrations of the difference between two apparently natural notions of privacy in quantum systems, relevant also to quantum key distribution (QKD): the older, naive one based on accessible information, contrasting with the new, composable one embracing the quantum nature of the eavesdropper's information. Assuming an additivity conjecture for constrained minimum output Renyi entropies, the techniques of the first part demonstrate a single-letter formula for the weak locking capacity of complements to cq-channels, coinciding with a general upper bound of Guha et al. for these channels. Furthermore, still assuming this additivity conjecture, this upper bound is given an operational interpretation for general channels as the maximum weak locking capacity of the channel activated by a suitable noiseless channel.


title: Rotating Drops with Helicoidal Symmetry

abstract: See http://youtu.be/Mf4IE8gWcJs for a YouTube video showing part of the results in this paper. We consider helicoidal immersions in the Euclidean space whose axis of symmetry is the z-axis that are solutions of the equation 2 H=\Lambda_0-a 1/2 R^2 where H is the mean curvature of the surface, R is the distance form the point in the surface to the z-axis and a is a real number. We refer to these surfaces as helicoidal rotating drops. We prove the existence of properly immersed solutions that contain the z-axis. We also show the existence of several families of embedded examples. We describe the set of possible solutions and we show that most of these solutions are not properly immerse and are dense in the region bounded by two concentric cylinders. We show that all properly immersed solutions, besides being invariant under a one parameter helicoidal group, they are invariant under a cyclic group of rotations of the variables x and y. The second variation of energy for the volume constrained problem with Dirichlet boundary conditions is also studied.


title: Secret-Sharing Matroids need not be Algebraic

abstract: We combine some known results and techniques with new ones to show that there exists a non-algebraic, multi-linear matroid. This answers an open question by Matus (Discrete Mathematics 1999), and an open question by Pendavingh and van Zwam (Advances in Applied Mathematics 2013). The proof is constructive and the matroid is explicitly given.


title: Lipkin method of particle-number restoration to higher orders

abstract: Background: On the mean-field level, pairing correlations are incorporated through the Bogoliubov-Valatin transformation, whereupon the particle degrees of freedom are replaced by quasiparticles. This approach leads to a spontaneous breaking of the particle-number symmetry and mixing of states with different particle numbers. In order to restore the particle number, various methods have been employed, which are based on projection approaches before or after variation. Approximate variation-after-projection (VAP) schemes, utilizing the Lipkin method, have mostly been used within the Lipkin-Nogami prescription. Purpose: Without recurring to the Lipkin-Nogami prescription, and using instead states rotated in the gauge space, we derive the Lipkin method of particle-number restoration up to sixth order and we test the convergence and accuracy of the obtained expansion. Methods: We perform self-consistent calculations using the higher-order Lipkin method to restore the particle-number symmetry in the framework of superfluid nuclear energy-density functional theory. We also apply the Lipkin method to a schematic exactly solvable two-level pairing model. Results: Calculations performed in open-shell tin and lead isotopes show that the Lipkin method converges at fourth order and satisfactorily reproduces the VAP ground-state energies and energy kernels. Near closed shells, the higher-order Lipkin method cannot be applied because of a non-analytic kink in the ground-state energies in function of the particle number. Conclusions: In open-shell nuclei, the higher-order Lipkin method provides a good approximation to the exact VAP energies. The method is computationally inexpensive, making it particularly suitable, for example, for future optimizations of the nuclear energy-density functionals and simultaneous restoration of different symmetries.


title: Stem Cell Transplantation As A Dynamical System: Are Clinical Outcomes Deterministic?

abstract: Outcomes in stem cell transplantation (SCT) are modeled using probability theory. However the clinical course following SCT appears to demonstrate many characteristics of dynamical systems, especially when outcomes are considered in the context of immune reconstitution. Dynamical systems tend to evolve over time according to mathematically determined rules. Characteristically, the future states of the system are predicated on the states preceding them, and there is sensitivity to initial conditions. In SCT, the interaction between donor T cells and the recipient may be considered as such a system in which, graft source, conditioning and early immunosuppression profoundly influence immune reconstitution over time. This eventually determines clinical outcomes, either the emergence of tolerance or the development of graft versus host disease. In this paper parallels between SCT and dynamical systems are explored and a conceptual framework for developing mathematical models to understand disparate transplant outcomes is proposed.


title: Time evolution of cascade decay

abstract: We study non-perturbatively the time evolution of cascade decay for generic fields $\pi \rightarrow \phi_1\phi_2\rightarrow \phi_2\chi_1\chi_2$ and obtain the time dependence of amplitudes and populations for the resonant and final states. We analyze in detail the different time scales and the manifestation of unitary time evolution in the dynamics of production and decay of resonant intermediate and final states. The probability of occupation (population) "flows" as a function of time from the initial to the final states. When the decay width of the parent particle $\Gamma_\pi$ is much larger than that of the intermediate resonant state $\Gamma_{\phi_1}$ there is a "bottleneck" in the flow, the population of resonant states builds up to a maximum at $t^* = \ln[\Gamma_\pi/\Gamma_{\phi_1}]/(\Gamma_\pi-\Gamma_{\phi_1})$ nearly saturating unitarity and decays to the final state on the longer time scale $1/\Gamma_{\phi_1}$. As a consequence of the wide separation of time scales in this case the cascade decay can be interpreted as evolving sequentially $\pi \rightarrow \phi_1\phi_2; ~ \phi_1\phi_2\rightarrow \phi_2\chi_1\chi_2$. In the opposite limit the population of resonances ($\phi_1$) does not build up substantially and the cascade decay proceeds almost directly from the initial parent to the final state without resulting in a large amplitude of the resonant state. An alternative but equivalent non-perturbative method useful in cosmology is presented. Possible phenomenological implications for heavy sterile neutrinos as resonant states and consequences of quantum entanglement and correlations in the final state are discussed.


title: A Framework for Hybrid Systems with Denial-of-Service Security Attack

abstract: Hybrid systems are integrations of discrete computation and continuous physical evolution. The physical components of such systems introduce safety requirements, the achievement of which asks for the correct monitoring and control from the discrete controllers. However, due to denial-of-service security attack, the expected information from the controllers is not received and as a consequence the physical systems may fail to behave as expected. This paper proposes a formal framework for expressing denial-of-service security attack in hybrid systems. As a virtue, a physical system is able to plan for reasonable behavior in case the ideal control fails due to unreliable communication, in such a way that the safety of the system upon denial-of-service is still guaranteed. In the context of the modeling language, we develop an inference system for verifying safety of hybrid systems, without putting any assumptions on how the environments behave. Based on the inference system, we implement an interactive theorem prover and have applied it to check an example taken from train control system.


title: Differential equations for singular values of products of Ginibre random matrices

abstract: It was proved by Akemann, Ipsen and Kieburg that squared singular values of products of $M$ complex Ginibre random matrices form a determinantal point process whose correlation kernel is expressible in terms of Meijer's $G$-functions. Kuijlaars and Zhang recently showed that at the edge of the spectrum, this correlation kernel has a remarkable scaling limit $K_M(x,y)$ which can be understood as a generalization of the classical Bessel kernel of Random Matrix Theory. In this paper we investigate the Fredholm determinant of the operator with the kernel $K_M(x,y)\chi_J(y)$, where $J$ is a disjoint union of intervals, $J=\cup_j(a_{2j-1},a_{2j})$, and $\chi_J$ is the characteristic function of the set $J$. This Fredholm determinant is equal to the probability that $J$ contains no particles of the limiting determinantal point process defined by $K_M(x,y)$ (the gap probability). We derive Hamiltonian differential equations associated with the corresponding Fredholm determinant, and relate them with the monodromy preserving deformation equations of the Jimbo, Miwa, Mori, Ueno and Sato theory. In the special case $J=(0,s)$ we give a formula for the gap probability in terms of a solution of a system of non-linear ordinary differential equations.


title: Dimensional Analysis and the Time Required to Urinate

abstract: According to the recently discovered 'Law of Urination', mammals, ranging in size from mice to elephants, take, on the average, 21s to urinate. We attempt to gain insights into the physical processes responsible for this uniformity using simple dimensional analysis. We assume that the biological apparatus for urination in mammals simply scales with linear size, and consider the scenarios where the driving force is gravity or elasticity, and where the response is dominated by inertia or viscosity. We ask how the time required for urination depends on the length scale, and find that for the time to be independent of body size, the dominant driving force must be elasticity, and the dominant response viscosity. Our note demonstrates that dimensional analysis can indeed readily give insights into complex physical and biological processes.