title: Universal electric current of interacting resonant-level models with asymmetric interactions: An extension of the Landauer formula

abstract: We study the electron transport in open quantum-dot systems described by the interacting resonant-level models with Coulomb interactions. We consider the situation in which the quantum dot is connected to the left and right leads asymmetrically. We exactly construct many-electron scattering eigenstates for the two-lead system, where two-body bound states appear as a consequence of one-body resonances and the Coulomb interactions. By using an extension of the Landauer formula, we calculate the average electric current for the system under bias voltages in the first order of the interaction parameters. Through a renormalization-group technique, we arrive at the universal electric current, where we observe the suppression of the electric current for large bias voltages, i.e., negative differential conductance. We find that the suppressed electric current is restored by the asymmetry of the system parameters.


title: A Resilient Quantum Secret Sharing Scheme

abstract: A resilient secret sharing scheme is supposed to generate the secret correctly even after some shares are damaged. In this paper, we show how quantum error correcting codes can be exploited to design a resilient quantum secret sharing scheme, where a quantum state is shared among more than one parties.


title: Precise measurement of the 222Rn half-life: a probe to monitor the stability of radioactivity

abstract: We give the results of a study on the 222Rn decay we performed in the Gran Sasso Laboratory (LNGS) by detecting the gamma rays from the radon progeny. The motivation was to monitor the stability of radioactivity measuring several times per year the half-life of a short lifetime (days) source instead of measuring over a long period the activity of a long lifetime (tens or hundreds of years) source. In particular, we give a possible reason of the large periodical fluctuations in the count rate of the gamma rays due to radon inside a closed canister which has been described in literature and which has been attributed to a possible influence of a component in the solar irradiation affecting the nuclear decay rates. We then provide the result of four half-life measurements we performed underground at LNGS in the period from May 2014 to January 2015 with radon diffused into olive oil. Briefly, we did not measure any change of the 222Rn half-life with a 8*10^-5 precision. Finally, we provide the most precise value for the 222Rn half-life: 3.82146(16){stat}(4){syst} days.


title: Gibbs-Ringing Artifact Removal Based on Local Subvoxel-shifts

abstract: Gibbs-ringing is a well known artifact which manifests itself as spurious oscillations in the vicinity of sharp image transients, e.g. at tissue boundaries. The origin can be seen in the truncation of k-space during MRI data-acquisition. Consequently, correction techniques like Gegenbauer reconstruction or extrapolation methods aim at recovering these missing data. Here, we present a simple and robust method which exploits a different view on the Gibbs-phenomena. The truncation in k-space can be interpreted as a convolution with a sinc-function in image space. Hence, the severity of the artifacts depends on how the sinc-function is sampled. We propose to re-interpolate the image based on local, subvoxel shifts to sample the ringing pattern at the zero-crossings of the oscillating sinc-function. With this, the artifact can effectively and robustly be removed with a minimal amount of smoothing.


title: The Tate spectrum of the higher real $K$-theories at height $n=p-1$

abstract: Let $E_n$ be Morava $E$-theory and let $G \subset G_n$ be a finite subgroup of $G_n$, the extended Morava stabilizer group. Let $E_{n}^{tG}$ be the Tate spectrum, defined as the cofiber of the norm map $N:(E_n)_{hG} \to E_n^{hG}$. We use the Tate spectral sequence to calculate $\pi_*E_{p-1}^{tG}$ for $G$ a maximal finite $p$-subgroup, and $p$ an odd prime. We show that $E_{p-1}^{tG} \simeq \ast$, so that the norm map gives equivalence between homotopy fixed point and homotopy orbit spectra. Our methods also give a calculation of $\pi_*E_{p-1}^{hG}$, which is a folklore calculation of Hopkins and Miller.


title: The population of white dwarf binaries with hot subdwarf companions

abstract: Hot subdwarfs (sdBs) are core helium-burning stars, which lost almost their entire hydrogen envelope in the red-giant phase. Since a high fraction of those stars are in close binary systems, common envelope ejection is an important formation channel. We identified a total population of 51 close sdB+WD binaries based on time-resolved spectroscopy and multi-band photometry, derive the WD mass distribution and constrain the future evolution of these systems. Most WDs in those binaries have masses significantly below the average mass of single WDs and a high fraction of them might therefore have helium cores. We found 12 systems that will merge in less than a Hubble time and evolve to become either massive C/O WDs, AM\,CVn systems, RCrB stars or even explode as supernovae type Ia.


title: Sufficient Covariate, Propensity Variable and Doubly Robust Estimation

abstract: Statistical causal inference from observational studies often requires adjustment for a possibly multi-dimensional variable, where dimension reduction is crucial. The propensity score, first introduced by Rosenbaum and Rubin, is a popular approach to such reduction. We address causal inference within Dawid's decision-theoretic framework, where it is essential to pay attention to sufficient covariates and their properties. We examine the role of a propensity variable in a normal linear model. We investigate both population-based and sample-based linear regressions, with adjustments for a multivariate covariate and for a propensity variable. In addition, we study the augmented inverse probability weighted estimator, involving a combination of a response model and a propensity model. In a linear regression with homoscedasticity, a propensity variable is proved to provide the same estimated causal effect as multivariate adjustment. An estimated propensity variable may, but need not, yield better precision than the true propensity variable. The augmented inverse probability weighted estimator is doubly robust and can improve precision if the propensity model is correctly specified.


title: Carter Subgroups, Amalgams, Simple Groups, and the Zp*-theorem

abstract: We consider an amalgam of groups constructed from fusion systems for different odd primes p and q. This amalgam contains a self-normalizing cyclic subgroup of order pq and isolated elements of order p and q.


title: Recovering Dirac systems with singularities in interior points

abstract: We study the non-selfadjoint Dirac system on a finite interval having non-integrable regular singularities in interior points with additional matching conditions at these points. Properties of spectral characteristics are established, and the inverse spectral problem is investigated. We provide a constructive procedure for the solution of the inverse problem, and prove its uniqueness. Moreover, necessary and sufficient conditions for the global solvability of this nonlinear inverse problem are obtained.


title: Flow Patterns around Dark Matter Halos: the Link between Halo Dynamical Properties and Large Scale Tidal Field

abstract: We study how halo intrinsic dynamical properties are linked to their formation processes for halos in two mass ranges, $10^{12}-10^{12.5}h^{-1}{\rm M_\odot}$ and $\ge 10^{13}h^{-1}{\rm M_\odot}$, and how both are correlated with the large scale tidal field within which the halos reside at present. Halo merger trees obtained from cosmological $N$-body simulations are used to identify infall halos that are about to merge with their hosts. We find that the tangential component of the infall velocity increases significantly with the strength of the local tidal field, but no strong correlation is found for the radial component. These results can be used to explain how the internal velocity anisotropy and spin of halos depend on environment. The position vectors and velocities of infall halos are aligned with the principal axes of the local tidal field, and the alignment depends on the strength of the tidal field. Opposite accretion patterns are found in weak and strong tidal fields, in the sense that in a weak field the accretion flow is dominated by radial motion within the local structure, while a large tangential component is present in a strong field. These findings can be used to understand the strong alignments we find between the principal axes of the internal velocity ellipsoids of halos and the local tidal field, and their dependence on the strength of tidal field. They also explain why halo spin increases with the strength of local tidal field, but only in weak tidal fields does the spin-tidal field alignment follow the prediction of the tidal torque theory. We discuss how our results may be used to understand the spins of disk galaxies and velocity structures of elliptical galaxies and their correlations with large-scale structure.


title: A realizability-preserving high-order kinetic scheme using WENO reconstruction for entropy-based moment closures of linear kinetic equations in slab geometry

abstract: We develop a high-order kinetic scheme for entropy-based moment models of a one-dimensional linear kinetic equation in slab geometry. High-order spatial reconstructions are achieved using the weighted essentially non-oscillatory (WENO) method, and for time integration we use multi-step Runge-Kutta methods which are strong stability preserving and whose stages and steps can be written as convex combinations of forward Euler steps. We show that the moment vectors stay in the realizable set using these time integrators along with a maximum principle-based kinetic-level limiter, which simultaneously dampens spurious oscillations in the numerical solutions. We present numerical results both on a manufactured solution, where we perform convergence tests showing our scheme converges of the expected order up to the numerical noise from the numerical optimization, as well as on two standard benchmark problems, where we show some of the advantages of high-order solutions and the role of the key parameter in the limiter.


title: Light axial vector mesons

abstract: Inspired by the abundant experimental observation of axial-vector states, we study whether the observed axial-vector states can be categorized into the conventional axial-vector meson family. In this paper we carry out an analysis based on the mass spectra and two-body Okubo-Zweig-Iizuka-allowed decays. Besides testing the possible axial-vector meson assignments, we also predict abundant information for their decays and the properties of some missing axial-vector mesons, which are valuable for further experimental exploration of the observed and predicted axial-vector mesons.


title: Inverse Spectral Problems for Sturm-Liouville Operators on Hedgehog-type Graphs with General Matching Conditions

abstract: Boundary value problems on hedgehog-type graphs for Sturm-Liouville differential operators with general matching conditions are studied. We investigate inverse spectral problems of recovering the coefficients of the differential equation from the spectral data. For this inverse problem we prove a uniqueness theorem and provide a procedure for constructing its solution.


title: Confidence intervals for AB-test

abstract: AB-testing is a very popular technique in web companies since it makes it possible to accurately predict the impact of a modification with the simplicity of a random split across users. One of the critical aspects of an AB-test is its duration and it is important to reliably compute confidence intervals associated with the metric of interest to know when to stop the test. In this paper, we define a clean mathematical framework to model the AB-test process. We then propose three algorithms based on bootstrapping and on the central limit theorem to compute reliable confidence intervals which extend to other metrics than the common probabilities of success. They apply to both absolute and relative increments of the most used comparison metrics, including the number of occurrences of a particular event and a click-through rate implying a ratio.


title: Probability distribution functions of gas in M31 and M51

abstract: We present probability distribution functions (PDFs) of the surface densities of ionized and neutral gas in the nearby spiral galaxies M31 and M51, as well as of dust emission and extinction Av in M31. The PDFs are close to lognormal and those for HI and Av in M31 are nearly identical. However, the PDFs for H2 are wider than the HI PDFs and the M51 PDFs have larger dispersions than those for M31. We use a simple model to determine how the PDFs are changed by variations in the line-of-sight (LOS) pathlength L through the gas, telescope resolution and the volume filling factor of the gas, f_v. In each of these cases the dispersion sigma of the lognormal PDF depends on the variable with a negative power law. We also derive PDFs of mean LOS volume densities of gas components in M31 and M51. Combining these with the volume density PDFs for different components of the ISM in the Milky Way (MW), we find that sigma decreases with increasing length L with an exponent of -0.76 +/- 0.06, which is steeper than expected. We show that the difference is due to variations in f_v. As f_v is similar in M31, M51 and the MW, the density structure in the gas in these galaxies must be similar. Finally, we demonstrate that an increase in f_v with increasing distance to the Galactic plane explains the decrease in sigma with latitude of the PDFs of emission measure and FUV emission observed for the MW.


title: Matter-wave interferometry with composite quantum objects

abstract: We discuss modern developments in quantum optics with organic molecules, clusters and nanoparticles -- in particular recent realizations of near-field matter-wave interferometry. A unified theoretical description in phase space allows us to describe quantum interferometry in position space and in the time domain on an equal footing. In order to establish matter-wave interferometers as a universal tool, which can accept and address a variety of nanoparticles, we elaborate on new quantum optical elements, such as diffraction gratings made of matter and light, as well as their absorptive and dispersive interaction with complex materials. We present Talbot-Lau interferometry (TLI), the Kapitza-Dirac-Talbot-Lau interferometer (KDTLI) and interferometry with pulsed optical gratings (OTIMA) as the most advanced devices to study the quantum wave nature of composite matter. These experiments define the current mass and complexity record in interferometric explorations of quantum macroscopicity and they open new avenues to quantum assisted metrology with applications in physical chemistry and biomolecular physics.


title: On Quasi-periodic Differential Pencils with Jump Conditions Inside the Interval

abstract: Non-self-adjoint second-order differential pencils on a finite interval with non-separated quasi-periodic boundary conditions and jump conditions are studied. We establish properties of spectral characteristics and investigate the inverse spectral problem of recovering the operator from its spectral data. For this inverse problem we prove the corresponding uniqueness theorem and provide an algorithm for constructing its solution.


title: Stability of magnetized strange quark matter in the MIT bag model with the density dependent bag pressure

abstract: The stability of magnetized strange quark matter (MSQM) is studied in the MIT bag model with the density dependent bag pressure. In the consistent thermodynamic description of MSQM, the quark chemical potentials, the total thermodynamic potential and the anisotropic pressure acquire the corresponding additional term proportional to the density derivative of the bag pressure. The model parameter space is determined, for which MSQM is absolutely stable, i.e., its energy per baryon is less than that of the most stable $^{56}$Fe nucleus under the zero external pressure and vanishing temperature. It is shown that there exists the magnetic field strength $H_{u\,max}$ at which the upper bound $B_\infty^u$ on the asymptotic bag pressure $B_\infty\equiv B(\varrho_B\gg \varrho_0$) ($\varrho_0$ being the nuclear saturation density) from the absolute stability window vanishes. The value of this field, \hbox{$H_{u\,max}\sim$$(1$--$3)\cdot10^{18}$}~G, represents the upper bound on the magnetic field strength, which can be reached in a strongly magnetized strange quark star. It is clarified how the absolute stability window and upper bound on the magnetic field strength are affected by varying the parameters in the Gaussian parametrization for the density dependence of the bag pressure.


title: Polyhedral Omega: A New Algorithm for Solving Linear Diophantine Systems

abstract: Polyhedral Omega is a new algorithm for solving linear Diophantine systems (LDS), i.e., for computing a multivariate rational function representation of the set of all non-negative integer solutions to a system of linear equations and inequalities. Polyhedral Omega combines methods from partition analysis with methods from polyhedral geometry. In particular, we combine MacMahon's iterative approach based on the Omega operator and explicit formulas for its evaluation with geometric tools such as Brion decompositions and Barvinok's short rational function representations. In this way, we connect two recent branches of research that have so far remained separate, unified by the concept of symbolic cones which we introduce. The resulting LDS solver Polyhedral Omega is significantly faster than previous solvers based on partition analysis and it is competitive with state-of-the-art LDS solvers based on geometric methods. Most importantly, this synthesis of ideas makes Polyhedral Omega the simplest algorithm for solving linear Diophantine systems available to date. Moreover, we provide an illustrated geometric interpretation of partition analysis, with the aim of making ideas from both areas accessible to readers from a wide range of backgrounds.


title: Near Optimal Subdivision Algorithms for Real Root Isolation

abstract: We describe a subroutine that improves the running time of any subdivision algorithm for real root isolation. The subroutine first detects clusters of roots using a result of Ostrowski, and then uses Newton iteration to converge to them. Near a cluster, we switch to subdivision, and proceed recursively. The subroutine has the advantage that it is independent of the predicates used to terminate the subdivision. This gives us an alternative and simpler approach to recent developments of Sagraloff (2012) and Sagraloff-Mehlhorn (2013), assuming exact arithmetic. The subdivision tree size of our algorithm using predicates based on Descartes's rule of signs is bounded by $O(n\log n)$, which is better by $O(n\log L)$ compared to known results. Our analysis differs in two key aspects. First, we use the general technique of continuous amortization from Burr-Krahmer-Yap (2009), and second, we use the geometry of clusters of roots instead of the Davenport-Mahler bound. The analysis naturally extends to other predicates.


title: Enumeration and structure of inhomogeneous graphs

abstract: We analyze a general model of weighted graphs, introduced by de Panafieu and Ravelomanana (2014) and similar to the "inhomogeneous graph model" of S\"oderberg (2002). Each vertex receives a "type" among a set of $q$ possibilities as well as a "weight" corresponding to this type, and each edge is weighted according to the types of the vertices it links. The weight of the graph is then the product of the weights of its vertices and edges. We investigate the sum of the weights of such graphs and prove that when the number of edges is small, almost all of them contain no component with more than one cycle. Those results allow us to give a new proof in a more general setting of a theorem of Wright (1961) on the enumeration of properly colored graphs. We also discuss applications related to social networks.


title: Images of the lines under the MS transformations and the Concept of Velocity in the DSR theories

abstract: The effect of the Maguejo-Smolin (MS) transformations on a straight line in the energy- momentum space will be studied. We will interpret the slope of this line as velocity $dE/dp$, which can leads to addition rule for the velocities in the MS doubly special relativity (DSR) case. Relation between two expressions $dE/dp$ and $p/E$ for velocity in the momentum space will be investigated more and the energy dependency of the velocities in the DSR theories is related to the geometrical properties of the lines under DSR transformations. The images of two parallel lines under the MS transformations will be studied and we will compute crossing point of these lines under the MS transformations in the energy-momentum space. The linear-fractional transformations don't keep parallelism. The crossing point is on a line in the energy-momentum space with a constant momentum $E_p/c$.


title: Mean-field solution of the Blume-Capel model under a random crystal field

abstract: In this work we investigate the Blume-Capel model with infinite-range ferromagnetic interactions and under the influence of a quenched disorder - a random crystal field. For a suitable choice of the random crystal field the model displays a wealth of multicritical behavior, continuous and first-order transition lines, as well as re-entrant behavior. The resulting phase diagrams show a variety of topologies as a function of the disorder parameter $\textit{p}$. A comparison with recent results on the Blume-Capel model in random crystal field is discussed.


title: Foreign Exchange Market Microstructure and the WM/Reuters 4pm Fix

abstract: A market fix serves as a benchmark for foreign exchange (FX) execution, and is employed by many institutional investors to establish an exact reference at which execution takes place. The currently most popular FX fix is the World Market Reuters (WM/R) 4pm fix. Execution at the WM/R 4pm fix is a service offered by FX brokers (normally banks), who deliver execution at the fix provided they obtain the trade order until a certain time prior to 4pm. In this paper, we study the market microstructure around 4pm. We demonstrate that market dynamics can be distinguished from other times during the day through increased volatility and size of movements. Our findings question the aggregate benefit to the client base of using the 4pm fix in its current form.


title: A family of stabilizer codes for $D({\mathbb Z}_2)$ anyons and Majorana modes

abstract: We study and generalize the class of qubit topological stabilizer codes that arise in the Abelian phase of the honeycomb lattice model. The resulting family of codes, which we call `matching codes' realize the same anyon model as the surface codes, and so may be similarly used in proposals for quantum computation. We show that these codes are particularly well suited to engineering twist defects that behave as Majorana modes. A proof of principle system that demonstrates the braiding properties of the Majoranas is discussed that requires only three qubits.


title: On the regularity set and angular integrability for the Navier--Stokes equation

abstract: We investigate the size of the regular set for suitable weak solutions of the Navier--Stokes equation, in the sense of Caffarelli--Kohn--Nirenberg. We consider initial data in weighted Lebesgue spaces with mixed radial-angular integrability, and we prove that the regular set increases if the data have higher angular integrability, invading the whole half space $\{t>0\}$ in an appropriate limit. In particular, we obtain that if the $L^{2}$ norm with weight $|x|^{-\frac12}$ of the data tends to 0, the regular set invades $\{t>0\}$.


title: Parametrizations in scalar-tensor theories of gravity and the limit of general relativity

abstract: We consider a general scalar-tensor theory of gravity and review briefly different forms it can be presented (different conformal frames and scalar field parametrizations). We investigate the conditions under which its field equations and the parametrized post-Newtonian parameters coincide with those of general relativity. We demonstrate that these so-called limits of general relativity are independent of the parametrization of the scalar field, although the transformation between scalar fields may be singular at the corresponding value of the scalar field. In particular, the limit of general relativity can equivalently be determined and investigated in the commonly used Jordan and Einstein frames.


title: One- and two-dimensional solitons supported by singular modulation of quadratic nonlinearity

abstract: We introduce a model of one- and two-dimensional (1D and 2D) optical media with the $\chi ^{(2)}$ nonlinearity whose local strength is subject to cusp-shaped spatial modulation, $\chi ^{(2)}\sim r^{-\alpha }$, with $\alpha >0$, which can be induced by spatially nonuniform poling. Using analytical and numerical methods, we demonstrate that this setting supports 1D and 2D fundamental solitons, at $\alpha <1$ and $\alpha <2$, respectively. The 1D solitons have a small instability region, while the 2D solitons have a stability region at $\alpha <0.5$ and are unstable at $\alpha >0.5$. 2D solitary vortices are found too. They are unstable, splitting into a set of fragments, which eventually merge into a single fundamental soliton pinned to the cusp. Spontaneous symmetry breaking of solitons is studied in the 1D system with a symmetric pair of the cusp-modulation peaks.


title: Solvation of the Fluorine Containing Anions and Their Lithium Salts in Propylene Carbonate and Dimethoxyethane

abstract: Electrolyte solutions based on the propylene carbonate (PC)-dimethoxyethane (DME) mixtures are of significant importance and urgency due to emergence of lithium-ion batteries. Solvation and coordination of the lithium cation in these systems have been recently attended in detail. However, analogous information concerning anions (tetrafluoroborate, hexafluorophosphate) is still missed. This work reports PM7-MD simulations (electronic-structure level of description) to include finite-temperature effects on the anion solvation regularities in the PC-DME mixture. The reported result evidences that the anions appear weakly solvated. This observation is linked to the absence of suitable coordination sites in the solvent molecules. In the concentrated electrolyte solutions, both BF4 and PF6 prefer to exist as neutral ion pairs (LiBF4, LiPF6).


title: Ultrafast helicity control of surface currents in topological insulators with near-unity fidelity

abstract: In recent years, a class of solid state materials, called three-dimensional topological insulators, has emerged. In the bulk, a topological insulator behaves like an ordinary insulator with a band gap. At the surface, conducting gapless states exist showing remarkable properties such as helical Dirac dispersion and suppression of backscattering of spin-polarized charge carriers. The characterization and control of the surface states via transport experiments is often hindered by residual bulk contributions yet at cryogenic temperatures. Here, we show that surface currents in Bi2Se3 can be controlled by circularly polarized light on a picosecond time scale with a fidelity near unity even at room temperature. We re-veal the temporal separation of such ultrafast helicity-dependent surface currents from photo-induced thermoelectric and drift currents in the bulk. Our results uncover the functionality of ultrafast optoelectronic devices based on surface currents in topological insulators.


title: Renormalization group invariant of lepton Yukawa couplings

abstract: By using quark Yukawa matrices only, we can construct renormalization invariants that are exact at the one-loop level in the standard model. One of them $I^q$ is accidentally consistent with unity, even though quark masses are strongly hierarchical. We calculate a lepton version of the invariant $I^l$ for Dirac and Majorana neutrino cases and find that $I^l$ can also be close to unity. For the Dirac neutrino and inverted hierarchy case, if the lightest neutrino mass is 3.0 meV to 8.8 meV, an equality $I^q=I^l$ can be satisfied. These invariants are not changed even if new particles couple to the standard model particles, as long as those couplings are generation independent.


title: Stability and Ramsey numbers for cycles and wheels

abstract: We study the structure of red-blue edge colorings of complete graphs, with no copies of the $n$-cycle $C_n$ in red, and no copies of the $n$-wheel $W_n = C_n \ast K_1$ in blue, for an odd integer $n$. Our first main result is that in any such coloring, deleting at most two vertices we obtain a vertex-partition of $G$ into three sets such that the edges inside the partition classes are red, and edges between partition classes are blue. As a second result, we obtain bounds for the Ramsey numbers of $r(C_{2k+1},W_{2j})$ for $k < j$ integers, which asymptotically confirm the values of $4j+1$, as it were conjectured by Zhang et al.


title: Phase diagram and optimal switching induced by spin Hall effect in a perpendicular magnetic layer

abstract: In a ferromagnet/heavy-metal bilayer device with strong spin Hall effect an in-plane current excites magnetic dynamics through spin torque. We analyze bilayers with perpendicular magnetization and calculate three-dimensional phase diagrams describing switching by external magnetic field at a fixed current. We then concentrate on the case of a field applied in the plane formed by the film normal and the current direction. Here we analytically study the evolution of both the conventional "up"/"down" magnetic equilibria and the additional equilibria created by the spin torque. Expressions for the stability regions of all equilibria are derived, and the nature of switching at each critical boundary is discussed. The qualitative picture obtained this way predicts complex hysteresis patterns that should occur in bilayers. By analyzing the phase portraits of the system we show that when the spin torque induced equilibrium exists, switching between "up" and "down" states proceeds through it as an intermediate state. Using numeric simulations we analyze the switching time and compare it to that of a conventional spin torque device with collinear magnetizations of the polarizer and the free layer.


title: Human diffusion and city influence

abstract: Cities are characterized by concentrating population, economic activity and services. However, not all cities are equal and a natural hierarchy at local, regional or global scales spontaneously emerges. In this work, we introduce a method to quantify city influence using geolocated tweets to characterize human mobility. Rome and Paris appear consistently as the cities attracting most diverse visitors. The ratio between locals and non-local visitors turns out to be fundamental for a city to truly be global. Focusing only on urban residents' mobility flows, a city to city network can be constructed. This network allows us to analyze centrality measures at different scales. New York and London play a predominant role at the global scale, while urban rankings suffer substantial changes if the focus is set at a regional level.

categories: physics.soc-ph cs.CY cs.SI physics.data-an stat.AP


title: Probability, Logic and Objectivity - The concept of probability of Carl Stumpf

abstract: We will show that Carl Stumpf's interpretation of the concept of probability is best understood as that of an objective Bayesian. First we analyse Stumpf's work in relation to that of his contemporary Johannes von Kries, and after that we will discuss various ways in which Stumpf's probability-concept has been construed. By showing that the construals of Stumpf's account by Hans Reichenbach, Richard von Mises and Andreas Kamlah are unfair - and at some points incorrect - we uncover the aspects that are essential to Stumpf's probability interpretation.


title: On the automorphism group of a binary $q$-analog of the Fano plane

abstract: The smallest set of admissible parameters of a $q$-analog of a Steiner system is $S_2[2,3,7]$. The existence of such a Steiner system -- known as a binary $q$-analog of the Fano plane -- is still open. In this article, the automorphism group of a putative binary $q$-analog of the Fano plane is investigated by a combination of theoretical and computational methods. As a conclusion, it is either rigid or its automorphism group is cyclic of order $2$, $3$ or $4$. Up to conjugacy in $\operatorname{GL}(7,2)$, there remains a single possible group of order $2$ and $4$, respectively, and two possible groups of order $3$. For the automorphisms of order $2$, we give a more general result which is valid for any binary $q$-Steiner triple system.


title: Bright sink-type localized states in exciton-polariton condensates

abstract: The family of one-dimensional localized solutions to dissipative nonlinear equations includes a variety of objects such as sources, sinks, shocks (kinks), and pulses. These states are in general accompanied by nontrivial density currents, which are not necessarily related to the movement of the object itself. We investigate the existence and physical properties of sink-type solutions in nonresonantly pumped exciton-polariton condensates modeled by an open-dissipative Gross-Pitaevskii equation. While sinks possess density profiles similar to bright solitons, they are qualitatively different objects as they exist in the case of repulsive interactions and represent a heteroclinic solution. We show that sinks can be created in realistic systems with appropriately designed pumping profiles. We also conclude that in two-dimensional configurations, due to the proliferation of vortices, sinks do not appear.


title: Frequency correlations in reflection from random media

abstract: We present a theoretical study of frequency correlations of light backscattered from a random scattering medium. This statistical quantity provides insight into the dynamics of multiple scattering processes accessible both, in theoretical and experimental investigations. For frequency correlations between field amplitudes, we derive a simple expression in terms of the path length distribution of the underlying backscattering processes. In a second step, we apply this relation to describe frequency correlations between intensities in the regime of weak disorder. Since, with increasing disorder strength, an unexplained breakdown of the angular structure of the frequency correlation function has recently been reported in experimental studies, we explore extensions of our model to the regime of stronger disorder. In particular, we show that closed scattering trajectories tend to suppress the angular dependence of the frequency correlation function.


title: Free energies of molecular clusters determined by guided mechanical disassembly

abstract: The excess free energy of a molecular cluster is a key quantity in models of the nucleation of droplets from a metastable vapour phase; it is often viewed as the free energy arising from the presence of an interface between the two phases. We show how this quantity can be extracted from simulations of the mechanical disassembly of a cluster using guide particles in molecular dynamics. We disassemble clusters ranging in size from 5 to 27 argon-like Lennard-Jones atoms, thermalised at 60 K, and obtain excess free energies, by means of the Jarzynski equality, that are consistent with previous studies. We only simulate the cluster of interest, in contrast to approaches that require a series of comparisons to be made between clusters differing in size by one molecule. We discuss the advantages and disadvantages of the scheme and how it might be applied to more complex systems.


title: Polaron Model of the Formation of Hydrated Electron States

abstract: A computer simulation of the formation of photoexcited electrons in water is performed within the framework of a dynamic model. The obtained results are discussed in comparison with experimental data and theoretical estimates.


title: Death ligand concentration and the membrane proximal signaling module regulate the type 1/ type 2 choice in apoptotic death signaling

abstract: Apoptotic death pathways are frequently activated by death ligand induction and subsequent activation of the membrane proximal signaling module. Death receptors cluster upon binding to death ligands, leading to formation of a membrane proximal death-inducing-signaling-complex (DISC). In this membrane proximal signalosome, initiator caspases (caspase 8) are processed resulting in activation of both type 1 and type 2 pathways of apoptosis signaling. How the type 1/type 2 choice is made is an important question in the systems biology of apoptosis signaling. In this study, we utilize a Monte Carlo based in silico approach to elucidate the role of membrane proximal signaling module in the type 1/type 2 choice of apoptosis signaling. Our results provide crucial mechanistic insights into the formation of DISC signalosome and caspase 8 activation. Increased concentration of death ligands was shown to correlate with increased type 1 activation. We also study the caspase 6 mediated system level feedback activation of apoptosis signaling and its role in the type 1/type 2 choice. Our results clarify the basis of cell-to-cell stochastic variability in apoptosis activation and ramifications of this issue is further discussed in the context of therapies for cancer and neurodegenerative disorders.


title: Revisiting perturbations in extended quasidilaton massive gravity

abstract: In this work we study the theory of extended quasidilaton massive gravity together with the presence of matter fields. After discussing the homogeneous and isotropic fully dynamical background equations, which governs the exact expansion history of the universe, we consider small cosmological perturbations around these general FLRW solutions. The stability of tensor, vector and scalar perturbations on top of these general background solutions give rise to slightly different constraints on the parameters of the theory than those obtained in the approximative assumption of the late-time asymptotic form of the expansion history, which does not correspond to our current epoch. This opens up the possibility of stable FLRW solutions to be compared with current data on cosmic expansion with the restricted parameter space based on theoretical ground.


title: Continuum Model of the One-Dimensional Holstein Bipolaron in DNA

abstract: The work is devoted to obtaining the 1D bipolaron functional by Holstein method for the continuum approximation. We analyzed the effect of electron correlations associated with the direct dependence of the wave function of the electron system from the electron-electron distance on the binding energy of the bipolaron. Numerical calculations of the bipolaron binding energy were performed according to the parameters of the system.


title: Exploring the Universe with Very High Energy Neutrinos

abstract: With the discovery of a high-energy neutrino flux in the 0.1 PeV to PeV range from beyond the Earth's atmosphere with the IceCube detector, neutrino astronomy has achieved a major breakthrough in the exploration of the high-energy universe. One of the main goals is the identification and investigation of the still mysterious sources of the cosmic rays which are observed at Earth with energies up to several $10^5$ PeV. In addition to being smoking-gun evidence for the presence of cosmic rays in a specific object, neutrinos escape even dense environments and can reach us from distant places in the universe, thereby providing us with a unique tool to explore cosmic accelerators. This article summarizes our knowledge about the observed astrophysical neutrino flux and current status of the search for individual cosmic neutrino sources. At the end, it gives an overview of plans for future neutrino telescope projects.


title: Renewal Approach to the Analysis of the Asynchronous State for Coupled Noisy Oscillators

abstract: We develop a framework in which the activity of nonlinear pulse-coupled oscillators is posed within the renewal theory. In this approach, the evolution of inter-event density allows for a self-consistent calculation that determines the asynchronous state and its stability. This framework, can readily be extended to the analysis of systems with more state variables. To exhibit this, we study a nonlinear pulse-coupled system, where couplings are dynamic and activity dependent. We investigate stability of this system and we show it undergoes a super-critical Hopf bifurcation to collective synchronization.


title: Locality-aware parallel block-sparse matrix-matrix multiplication using the Chunks and Tasks programming model

abstract: We present a method for parallel block-sparse matrix-matrix multiplication on distributed memory clusters. By using a quadtree matrix representation, data locality is exploited without prior information about the matrix sparsity pattern. A distributed quadtree matrix representation is straightforward to implement due to our recent development of the Chunks and Tasks programming model [Parallel Comput. 40, 328 (2014)]. The quadtree representation combined with the Chunks and Tasks model leads to favorable weak and strong scaling of the communication cost with the number of processes, as shown both theoretically and in numerical experiments. Matrices are represented by sparse quadtrees of chunk objects. The leaves in the hierarchy are block-sparse submatrices. Sparsity is dynamically detected by the matrix library and may occur at any level in the hierarchy and/or within the submatrix leaves. In case graphics processing units (GPUs) are available, both CPUs and GPUs are used for leaf-level multiplication work, thus making use of the full computing capacity of each node. The performance is evaluated for matrices with different sparsity structures, including examples from electronic structure calculations. Compared to methods that do not exploit data locality, our locality-aware approach reduces communication significantly, achieving essentially constant communication per node in weak scaling tests.


title: Heisenberg model and Rigged Configurations

abstract: We show a correspondence of all the solutions of the spin-1/2 isotropic Heisenberg model for N=12 to the rigged configurations based on the comparison of the set of Takahashi quantum numbers in lexicographical order with the set of riggings of the rigged configurations in co-lexicographical order.

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


title: Differential topology interacts with isoparametric foliations

abstract: In this note, we discuss the interactions between differential topology and isoparametric foliations, surveying some recent progress and open problems.


title: The Composite Twin Higgs scenario

abstract: Based on an explicit model, we propose and discuss the generic features of a possible implementation of the Twin Higgs program in the context of composite Higgs models. We find that the Twin Higgs quadratic divergence cancellation argument can be uplifted to a genuine protection of the Higgs potential, based on symmetries and selection rules, but only under certain conditions which are not fulfilled in some of the existing models. We also find that a viable scenario, not plagued by a massless Twin Photon, can be obtained by not gauging the Twin Hypercharge and taking this as the only source of Twin Symmetry breaking at a very high scale.


title: Lagrangian transport through surfaces in volume-preserving flows

abstract: Advective transport of scalar quantities through surfaces is of fundamental importance in many scientific applications. From the Eulerian perspective of the surface it can be quantified by the well-known integral of the flux density. The recent development of highly accurate semi-Lagrangian methods for solving scalar conservation laws and of Lagrangian approaches to coherent structures in turbulent (geophysical) fluid flows necessitate a new approach to transport from the (Lagrangian) material perspective. We present a Lagrangian framework for calculating transport of conserved quantities through a given surface in $n$-dimensional, fully aperiodic, volume-preserving flows. Our approach does not involve any dynamical assumptions on the surface or its boundary.

categories: nlin.CD math.CA math.DS nlin.SI physics.flu-dyn


title: Inverse diffraction for the Atmospheric Imaging Assembly in the Solar Dynamics Observatory

abstract: The Atmospheric Imaging Assembly in the Solar Dynamics Observatory provides full Sun images every 1 seconds in each of 7 Extreme Ultraviolet passbands. However, for a significant amount of these images, saturation affects their most intense core, preventing scientists from a full exploitation of their physical meaning. In this paper we describe a mathematical and automatic procedure for the recovery of information in the primary saturation region based on a correlation/inversion analysis of the diffraction pattern associated to the telescope observations. Further, we suggest an interpolation-based method for determining the image background that allows the recovery of information also in the region of secondary saturation (blooming).


title: Random walks on binary strings applied to the somatic hypermutation of B-cells

abstract: Within the germinal center in follicles, B-cells proliferate, mutate and differentiate, while being submitted to a powerful selection~: a micro-evolutionary mechanism at the heart of adaptive immunity. A new foreign pathogen is confronted to our immune system, the mutation mechanism that allows B-cells to adapt to it is called {\em somatic hypermutation}~: a programmed process of mutation affecting B-cell receptors at extremely high rate. By considering random walks on graphs, we introduce and analyze a simplified mathematical model in order to understand this extremely efficient learning process. The structure of the graph reflects the choice of the mutation rule. We focus on the impact of this choice on typical time-scales of the graphs' exploration. We derive explicit formulas to evaluate the expected hitting time to cover a given Hamming distance on the graphs under consideration. This characterizes the efficiency of these processes in driving antibody affinity maturation. In a further step we present a biologically more involved model and discuss its numerical outputs within our mathematical framework. We provide as well limitations and possible extensions of our approach.


title: On the $\ell$-adic Fourier transform and the determinant of the middle convolution

abstract: We study the relation of the middle convolution to the $\ell$-adic Fourier transformation in the \'etale context. Using Katz' work and Laumon's theory of local Fourier transformations we obtain a detailed description of the local monodromy and the determinant of Katz' middle convolution functor $\MC_\chi$ in the tame case. The theory of local $\epsilon$-constants then implies that the property of an \'etale sheaf of having an at most quadratic determinant is often preserved under $\MC_\chi$ if $\chi$ is quadratic.


title: Multiwave imaging in an enclosure with variable wave speed

abstract: In this paper we consider the mathematical model of thermo- and photo-acoustic tomography for the recovery of the initial condition of a wave field from knowledge of its boundary values. Unlike the free-space setting, we consider the wave problem in a region enclosed by a surface where an impedance boundary condition is imposed. This condition models the presence of physical boundaries such as interfaces or acoustic mirrors which reflect some of the wave energy back into the enclosed domain. By recognizing that the inverse problem is equivalent to a statement of boundary observability, we use control operators to prove the unique and stable recovery of the initial wave profile from knowledge of boundary measurements. Since our proof is constructive, we explicitly derive a solvable equation for the unknown initial condition. This equation can be solved numerically using the conjugate gradient method. We also propose an alternative approach based on the stabilization of waves. This leads to an exponentially and uniformly convergent Neumann series reconstruction when the impedance coefficient is not identically zero. In both cases, if well-known geometrical conditions are satisfied, our approaches are naturally suited for variable wave speed and for measurements on a subset of the boundary.


title: A Quantum Phase Transition in the Cosmic Ray Energy Distribution

abstract: We here argue that the "knee" of the cosmic ray energy distribution at $E_c \sim 1$ PeV represents a second order phase transition of cosmic proportions. The discontinuity of the heat capacity per cosmic ray particle is given by $\Delta c=0.450196\ k_B$. However the idea of a deeper critical point singularity cannot be ruled out by present accuracy in neither theory nor experiment. The quantum phase transition consists of cosmic rays dominated by bosons for the low temperature phase E<E_c and dominated by fermions for high temperature phase $E > E_c$. The low temperature phase arises from those nuclei described by the usual and conventional collective boson models of nuclear physics. The high temperature phase is dominated by protons. The transition energy $E_c$ may be estimated in terms of the photo-disintegration of nuclei.


title: The Theoretical Power Law Exponent for Electron and Positron Cosmic Rays: A Comment on the Recent Letter of the AMS Collaboration

abstract: In a recent letter, the AMS collaboration reported the detailed and extensive data concerning the distribution in energy of electron and positron cosmic rays. A central result of the experimental work resides in the energy regime $30\ {\rm GeV}< E < 1\ {\rm TeV}$ wherein the power law exponent of the energy distribution is measured to be $\alpha ({\rm experiment})=3.17$. In virtue of the Fermi statistics obeyed by electrons and positrons, a theoretical value was predicted as $\alpha ({\rm theory})=3.151374$ in very good agreement with experimental data. The consequences of this agreement between theory and experiment concerning the sources of cosmic ray electrons and positrons are briefly explored.


title: Anderson localization of Bogoliubov excitations on quasi-1D strips

abstract: Anderson localization of Bogoliubov excitations is studied for disordered lattice Bose gases in planar quasi-one-dimensional geometries. The inverse localization length is computed as function of energy by a numerical transfer-matrix scheme, for strips of different widths. These results are described accurately by analytical formulas based on a weak-disorder expansion of backscattering mean free paths.


title: Quasiseparable Hessenberg reduction of real diagonal plus low rank matrices and applications

abstract: We present a novel algorithm to perform the Hessenberg reduction of an $n\times n$ matrix $A$ of the form $A = D + UV^*$ where $D$ is diagonal with real entries and $U$ and $V$ are $n\times k$ matrices with $k\le n$. The algorithm has a cost of $O(n^2k)$ arithmetic operations and is based on the quasiseparable matrix technology. Applications are shown to solving polynomial eigenvalue problems and some numerical experiments are reported in order to analyze the stability of the approach


title: Shock Tube Design for High Intensity Blast Waves for Laboratory Testing of Armor and Combat Materiel

abstract: Shock tubes create simulated blast waves which can be directed and measured to study blast wave effects under laboratory conditions. It is desirable to increase available peak pressure from ~1 MPa to ~5 MPa to simulate closer blast sources and facilitate development and testing of personal and vehicle armors. Three methods were investigated to increase peak simulated blast pressure produced by an oxy-acetylene driven shock tube while maintaining suitability for laboratory studies. The first method is the addition of a Shchelkin spiral priming section which works by increasing the turbulent flow of the deflagration wave, thus increasing its speed and pressure. This approach increased the average peak pressure from 1.17 MPa to 5.33 MPa while maintaining a relevant pressure-time curve (Friedlander waveform). The second method is a bottleneck between the driving and driven sections. Coupling a 79 mm diameter driving section to a 53 mm driven section increased the peak pressure from 1.17 MPa to 2.25 MPa. Using a 103 mm driving section increased peak pressure to 2.64 MPa. The third method, adding solid fuel to the driving section with the oxy-acetylene, resulted in a peak pressure increase to 1.70 MPa.


title: Valued Workflow Satisfiability Problem

abstract: A workflow is a collection of steps that must be executed in some specific order to achieve an objective. A computerised workflow management system may enforce authorisation policies and constraints, thereby restricting which users can perform particular steps in a workflow. The existence of policies and constraints may mean that a workflow is unsatisfiable, in the sense that it is impossible to find an authorised user for each step in the workflow and satisfy all constraints. In this paper, we consider the problem of finding the "least bad" assignment of users to workflow steps by assigning a weight to each policy and constraint violation. To this end, we introduce a framework for associating costs with the violation of workflow policies and constraints and define the \emph{valued workflow satisfiability problem} (Valued WSP), whose solution is an assignment of steps to users of minimum cost. We establish the computational complexity of Valued WSP with user-independent constraints and show that it is fixed-parameter tractable. We then describe an algorithm for solving Valued WSP with user-independent constraints and evaluate its performance, comparing it to that of an off-the-shelf mixed integer programming package.


title: Parsimonious Tensor Response Regression

abstract: Aiming at abundant scientific and engineering data with not only high dimensionality but also complex structure, we study the regression problem with a multidimensional array (tensor) response and a vector predictor. Applications include, among others, comparing tensor images across groups after adjusting for additional covariates, which is of central interest in neuroimaging analysis. We propose parsimonious tensor response regression adopting a generalized sparsity principle. It models all voxels of the tensor response jointly, while accounting for the inherent structural information among the voxels. It effectively reduces the number of free parameters, leading to feasible computation and improved interpretation. We achieve model estimation through a nascent technique called the envelope method, which identifies the immaterial information and focuses the estimation based upon the material information in the tensor response. We demonstrate that the resulting estimator is asymptotically efficient, and it enjoys a competitive finite sample performance. We also illustrate the new method on two real neuroimaging studies.


title: Observations of the extreme runaway HD271791: nucleosynthesis in a core-collapse supernova

abstract: Some young, massive stars can be found in the Galactic halo. As star formation is unlikely to occur in the halo, they must have been formed in the disk and been ejected shortly afterwards. One explanation is a supernova in a tight binary system. The companion is ejected and becomes a runaway star. HD\,271791 is the kinematically most extreme runaway star known (Galactic restframe velocity $725 \pm 195\, \rm km\,s^{-1}$, which is even larger than the Galactic escape velocity). Moreover, an analysis of the optical spectrum showed an enhancement of the $\alpha$-process elements. This indicates the capture of supernova ejecta, and therefore an origin in a core-collapse supernova. As such high space velocities are not reached by the runaway stars in classical binary supernova ejection scenarios, a very massive but compact primary, probably of Wolf-Rayet type is required. HD\,271791 is therefore a perfect candidate for studying nucleosynthesis in a supernova of probably type Ibc. The goal of this project is to determine the abundances of a large number of elements from the $\alpha$-process, the iron group, and heavier elements by a quantitative analysis of the optical and UV spectral range. Detailed line-formation calculations are employed that account for deviations from local thermodynamic equilibrium (non-LTE). We intend to verify whether core-collapse supernova are a site of r-process element production. Here, we state the current status of the project.


title: A method to discriminate solar and antisolar differential rotation in high-precision light curves

abstract: Surface differential rotation (DR) is one major ingredient of the magnetic field generation process in the Sun and likely in other stars. The term solar-like differential rotation describes the observation that solar equatorial regions rotate faster than polar ones. The opposite effect of polar regions rotating faster than equatorial ones (termed as antisolar DR) has only been observed in a few stars, although there is evidence from theoretical dynamo models. We present a new method to detect the sign of DR (i.e. solar-like or antisolar DR) by analyzing long-term high-precision light curves with the Lomb-Scargle periodogram.We compute the Lomb-Scargle periodogram and identify a set of significant periods $P_k$, which we associate with active regions located at different latitudes on the the stellar surface. If detectable, the first harmonics ($P_k'$) of these periods were identified to compute their peak-height-ratios $r_k:=h(P_k')/h(P_k)$. Spots rotating at lower latitudes generate less sine-shaped light curves, which requires additional power in the harmonics, and results in larger ratios $r_k$. Comparing different ratios $r_k$ and the associated periods $P_k$ yields information about the spot latitudes, and reveals the sign of DR. We tested our method on different sets of synthetic light curves all exhibiting solar-like DR. The number of cases where our method detects antisolar DR is the false-positive rate of our method. Depending on the set of light curves, the noise level, the required minimum peak separation, and the presence or absence of spot evolution, our method fails to detect the correct sign in at most 20%. We applied our method to 50 Kepler G stars and found 21-34 stars with solar-like DR and 5-10 stars with antisolar DR, depending on the minimum peak separation.


title: Modeling Superionic Behavior of Plutonium Dioxide

abstract: The Bredig transition to the superionic phase indicated with the lambda-peak in Cp was highly expected for PuO2 as other actinide dioxides. However, least-square fit and local smoothing techniques applied to the experimental enthalpy data of plutonium dioxide in 80's could not detect a lambda-peak in specific heat that might be due to too scattered and insufficient experimental data. Therefore, this issue has not been yet put beyond the doubts. In the current article, a superionic model of PuO2 is developed with partially ionic model of a rigid ion potential. Thermophysical properties were calculated in constant pressure-temperature ensemble using molecular dynamics simulation. The Bredig transition with vicinity of a lambda-peak in specific heat was a successfully observed for the model system at about 2100K. Moreover, the experimental enthalpy change was well reproduced before and after the estimated transition temperature.


title: Some Cosmological Consequences of Weyl Invariance

abstract: We examine some Weyl invariant cosmological models in the framework of generalized dilaton gravity, in which the action is made of a set of $N$ conformally coupled scalar fields. It will be shown that when the FRW ansatz for the spacetime metric is assumed, the Ward identity for conformal invariance guarantees that the gravitational equations hold whenever the scalar fields EM do so. It follows that any scale factor can solve the theory provided a non-trivial profile for a dilaton field. In particular, accelerated expansion is a natural solution to the full set of equations.


title: Propagation of chaos, Wasserstein gradient flows and toric Kahler-Einstein metrics

abstract: Motivated by a probabilistic approach to Kahler-Einstein metrics we consider a general non-equilibrium statistical mechanics model in Euclidean space consisting of the stochastic gradient flow of a given (possibly singular) quasi-convex N-particle interaction energy. We show that a deterministic "macroscopic" evolution equation emerges in the large N-limit of many particles. This is a strengthening of previous results which required a uniform two-sided bound on the Hessian of the interaction energy. The proof uses the theory of weak gradient flows on the Wasserstein space. Applied to the setting of permanental point processes at "negative temperature" the corresponding limiting evolution equation yields a drift-diffusion equation, coupled to the Monge-Ampere operator, whose static solutions correspond to toric Kahler-Einstein metrics. This drift-diffusion equation is the gradient flow on the Wasserstein space of probability measures of the K-energy functional in Kahler geometry and it can be seen as a fully non-linear version of various extensively studied dissipative evolution equations and conservations laws, including the Keller-Segel equation and Burger's equation. We also obtain a real probabilistic (and tropical) analog of the complex geometric Yau-Tian-Donaldson conjecture in this setting. In a companion paper applications to singular pair interactions are given.


title: Control and generation of localized pulses in passively mode-locked semiconductor lasers

abstract: We show experimentally and theoretically that localized pulses can be generated from an electrically biased $200\,\mu$m multi-transverse mode Vertical-Cavity Surface-Emitting Laser. The device is passively mode-locked using optical feedback from a distant Resonant Saturable Absorber Mirror and it is operated below threshold. We observe multistability between the off solution and a large variety of pulsating solutions with different number and arrangements of pulses per round-trip, thus indicating that the mode-locked pulses are localized, i.e. mutually independent. We show that a modulation of the bias current allows controlling the number of the pulses travelling within the cavity, thus suggesting that our system can be operated as an arbitrary pattern generator of 10 ps pulses and 1 W peak power.


title: TRiPoD (Temporal Relationalism incorporating Principles of Dynamics)

abstract: Temporal Relationalism is that there is no time for the universe as a whole at the primary level. Time emerges rather at a secondary level; one compelling idea for this is Mach's: that time is to be abstracted from change. Temporal Relationalism leads to, and better explains, the well-known Frozen Formalism Problem encountered in GR and other background-independent theories at the quantum level. Abstraction from change is then a type of emergent time resolution of this. Moreover, the Frozen Formalism Problem is but one of the many Problem of Time facets, which are notoriously interconnected. The current article concerns modifications of physical formalism which ensure that once Temporal Relationalism is resolved, it stays incorporated. At the classical level, this involves modifying much of the Principles of Dynamics. I first introduce the anti-Routhian to complete the Legendre square of Lagrangian, Hamiltonian and Routhian. I next pass from velocities $\dot{Q}\mbox{}^{A}$ to changes d$Q^{A}$. Then Lagrangians are supplanted by Jacobi arc elements, Euler--Lagrange equations by Jacobi--Mach ones, and momentum requires redefining but actions remain unchanged. A differential (d) version of the Hamiltonian is required, giving rise to a variant of the Dirac approach based on a d-almost Hamiltonian subcase of the d-anti Routhian. On the other hand, the forms of the constraints themselves, and of Hamilton--Jacobi theory, remain unaltered.


title: On short time existence of Lagrangian mean curvature flow

abstract: We consider a short time existence problem motivated by a conjecture of Joyce. Specifically we prove that given any compact Lagrangian $L\subset \mathbb{C}^n$ with a finite number of singularities, each asymptotic to a pair of non-area-minimising, transversally intersecting Lagrangian planes, there is a smooth Lagrangian mean curvature flow existing for some positive time, that attains $L$ as $t \searrow 0$ as varifolds, and smoothly locally away from the singularities.


title: Nearby Low-Mass Hypervelocity Stars

abstract: Hypervelocity stars are those that have speeds exceeding the escape speed and are hence unbound from the Milky Way. We investigate a sample of low-mass hypervelocity candidates obtained using data from the high-precision SDSS Stripe 82 catalogue, which we have combined with spectroscopy from the 200-inch Hale Telescope at Palomar Observatory. We find four good candidates, but without metallicities it is difficult to pin-down their distances and therefore total velocities. Our best candidate has a significant likelihood that it is escaping the Milky Way for a wide-range of metallicities.


title: Molecular Dynamics Simulations of Solutions at Constant Chemical Potential

abstract: Molecular Dynamics studies of chemical processes in solution are of great value in a wide spectrum of applications, which range from nano-technology to pharmaceutical chemistry. However, these calculations are affected by severe finite-size effects, such as the solution being depleted as the chemical process proceeds, which influence the outcome of the simulations. To overcome these limitations, one must allow the system to exchange molecules with a macroscopic reservoir, thus sampling a Grand-Canonical ensemble. Despite the fact that different remedies have been proposed, this still represents a key challenge in molecular simulations. In the present work we propose the Constant Chemical Potential Molecular Dynamics (C$\mu$MD) method, which introduces an external force that controls the environment of the chemical process of interest. This external force, drawing molecules from a finite reservoir, maintains the chemical potential constant in the region where the process takes place. We have applied the C$\mu$MD method to the paradigmatic case of urea crystallization in aqueous solution. As a result, we have been able to study crystal growth dynamics under constant supersaturation conditions, and to extract growth rates and free-energy barriers.


title: Continued fractions and parametric geometry of numbers

abstract: Recently, W. M. Schmidt and L. Summerer developed a new theory called Parametric Geometry of Numbers which approximates the behaviour of the successive minima of a family of convex bodies in $\mathbb{R}^{n}$ related to the problem of simultaneous rational approximation to given real numbers. In the case of one number, we show that the qualitative behaviour of the minima reflects the continued fraction expansion of the smallest distance from this number to an integer.


title: Electronic transition from graphite to graphene via controlled movement of the top layer with scanning tunneling microscopy

abstract: A series of measurements using a technique called electrostatic-manipulation scanning tunneling microscopy (EM-STM) were performed on a highly-oriented pyrolytic graphite surface. The electrostatic interaction between the STM tip and the sample can be tuned to produce both reversible and irreversible large-scale movement of the graphite surface. Under this influence, atomic-resolution STM images reveal that a continuous electronic transition from triangular symmetry, where only alternate atoms are imaged, to hexagonal symmetry can be systematically controlled. Density functional theory (DFT) calculations reveal that this transition can be related to vertical displacements of the top layer of graphite relative to the bulk. Evidence for horizontal shifts in the top layer of graphite is also presented. Excellent agreement is found between experimental STM images and those simulated using DFT.


title: Deconfining temperatures in SO(N) and SU(N) gauge theories

abstract: We present our current results for the deconfining temperatures in SO(N) gauge theories in 2+1 dimensions. SO(2N) theories may help us to understand QCD at finite chemical potential since there is a large-N orbifold equivalence between SO(2N) QCD-like theories and SU(N) QCD, and SO(2N) theories do not have the sign problem present in QCD. We show that the deconfining temperatures in these two theories match at the large-N limit. We also present results for SO(2N+1) gauge theories and compare results for SO(6) with SU(4) gauge theories, which have the same Lie algebras but different centres.


title: An Action Principle for Relativistic MHD

abstract: A covariant action principle for ideal relativistic magnetohydrodynamics (MHD) in terms of natural Eulerian field variables is given. This is done by generalizing the covariant Poisson bracket theory of Marsden et al., which uses a noncanonical bracket to effect constrained variations of an action functional. Various implications and extensions of this action principle are also discussed. Two significant by-products of this formalism are the introduction of a new divergence-free 4-vector variable for the magnetic field, and a new Lie-dragged form for the theory.


title: On the symplectic realization of Poisson-Nijenhuis manifolds

abstract: We consider the problem of the symplectic realization of a Poisson-Nijenhuis manifold. By applying a new technique developed by M. Crainic and I. Marcut for the study of the above problem in the case of a Poisson manifold, we establish the existence, under a condition, of a nondegenerate Poisson-Nijenhuis structure on an open neighborhood of the zero-section of the cotangent bundle of the manifold, which symplectizes the initial structure. Additionally, we present some examples.


title: Post measurement bipartite entanglement entropy in conformal field theories

abstract: We derive exact formulas for bipartite von Neumann entanglement entropy after partial projective local measurement in $1+1$ dimensional conformal field theories with periodic and open boundary conditions. After defining the set up we will check numerically the validity of our results in the case of Klein-Gordon field theory (coupled harmonic oscillators) and spin-$1/2$ XX chain in a magnetic field. The agreement between analytical results and the numerical calculations is very good. We also find a lower bound for localizable entanglement in coupled harmonic oscillators.


title: How does a flexible chain of active particles swell?

abstract: We study the swelling of a flexible linear chain composed of active particles by analytical theory and computer simulation. Three different situations are considered: a free chain, a chain confined to an external harmonic trap, and a chain dragged at one end. First we consider an ideal chain with harmonic springs and no excluded volume between the monomers. The Rouse model of polymers is generalized to the case of self-propelled monomers and solved analytically. The swelling, as characterized by the spatial extension of the chain, scales with the monomer number defining a Flory exponent $\nu$ which is $\nu =1/2, 0, 1$ in the three different situations. As a result, we find that activity does not change the Flory exponent but affects the prefactor of the scaling law. This can be quantitatively understood by mapping the system onto an equilibrium chain with a higher effective temperature such that the chain swells under an increase of the self-propulsion strength. We then use computer simulations to study the effect of self-avoidance on active polymer swelling. In the three different situations, the Flory exponent is now $\nu = 3/4, 1/4, 1$ and again unchanged under self-propulsion. However, the chain extension behaves non-monotonic in the self-propulsion strength.


title: Quantification of Entanglement Entropies for Doubly Excited States in Helium

abstract: In this work, we study the quantum entanglement for doubly excited resonance states in helium by using highly correlated Hylleraas type functions to represent such states of the two-electron system. The doubly-excited resonance states are determined by calculation of density of resonance states under the framework of the stabilization method. The spatial (electron-electron orbital) entanglement measures for the low-lying doubly excited 2s2, 2s3s, and 2p2 1Se states are carried out. Once a resonance state wave function is obtained, the linear entropy and von Neumann entropy for such a state are quantified using the Schmidt-Slater decomposition method. To check the consistence, linear entropy is also determined by solving analytically the needed four-electron (12-dimensional) integrals.

categories: quant-ph physics.atm-clus physics.atom-ph physics.chem-ph


title: States generated in the $K$-multi-$\rho$ interactions

abstract: In the present work we use three-body interaction formalism to investigate the $K$-multi-$\rho$ interactions. First, we reproduced the resonances $f_2(1270)$ and $K_1(1270)$ in the $\rho \rho$ and $\rho K$ two-body interactions respectively, as the clusters of the fixed center approximation. Then, we study the three-body $K -\rho \rho (f_2)$ and $\rho - \rho K (K_1)$ interactions with the fixed center approximation of the Faddeev equations. Furthermore, we extrapolate the formalism to study the four-body, five-body and six-body systems, containing one $K$ meson and multi-$\rho$ mesons. In our research, without introducing any free parameters, we generate the $K_2 (1770)$ state in the three-body interaction with the mass of 1707 MeV and a width about 113 MeV, which are consistent with the experiments. We also find a clear resonant structure in our results of the five-body interaction, with a mass 2505 MeV and a width about 32 MeV or more, which is associated to the $K_4 (2500)$ state, where we obtain consistent results with the experimental findings. Furthermore, we predict some new states in the other many body interactions, $K_3 (2080)$, $K_5 (2670)$ (isospin $I=1/2$), and $K_4 (2640)$ (isospin $I=3/2$), with uncertainties.


title: The size of the giant component in random hypergraphs

abstract: The phase transition in the size of the giant component in random graphs is one of the most well-studied phenomena in random graph theory. For hypergraphs, there are many possible generalisations of the notion of a component, and for all but the simplest example, the phase transition phenomenon was first proved by Cooley, Kang and Person. In this paper we build on this and determine the asymptotic size of the unique giant component.


title: Time and spatial parity operations with trapped ions

abstract: We propose a physical implementation of time and spatial parity transformations, as well as Galilean boosts, in a trapped-ion quantum simulator. By embedding the simulated model into an enlarged simulating Hilbert space, these fundamental symmetry operations can be fully realized and measured with ion traps. We illustrate our proposal with analytical and numerical techniques of prototypical examples with state-of-the-art trapped-ion platforms. These results pave the way for the realization of time and spatial parity transformations in other models and quantum platforms.


title: Death by Dynamics: Planetoid-Induced Explosions on White Dwarfs

abstract: At intervals as short as ten thousand years, each white dwarf (WD) passes within a solar radius of a planetoid, i.e., a comet, asteroid, or planet. Gravitational tidal forces tear the planetoid apart; its metal-rich debris falls onto the WD, enriching the atmosphere. A third of WDs exhibit atmospheric "pollution". For roughly every hundred planetoid disruptions, a planetoid collides with a WD. We simulate a small number of collisions, in which "death-by-dynamics" refers to the fate of the planetoid. We also compute the energies and likely durations of a broad sample of collision events, and identify detection strategies at optical and X-ray wavelengths. Collisions with the most massive planetoids can be detected in external galaxies. Some may trigger nuclear burning. If one in $\sim 10^7-10^8$ of WD-planetoid collisions creates the conditions needed for a Type Ia supernova (SN~Ia), "death-by-dynamics" would also refer to the fate of the WD, and could provide a novel channel for the production of SN~Ia. We consider the circumstances under which the rate of SNe~Ia can be increased by interactions with planetoids.


title: New G2 holonomy cones and exotic nearly Kaehler structures on the 6-sphere and the product of a pair of 3-spheres

abstract: There is a rich theory of so-called (strict) nearly Kaehler manifolds, almost-Hermitian manifolds generalising the famous almost complex structure on the 6-sphere induced by octonionic multiplication. Nearly Kaehler 6-manifolds play a distinguished role both in the general structure theory and also because of their connection with singular spaces with holonomy group the compact exceptional Lie group G2: the metric cone over a Riemannian 6-manifold M has holonomy contained in G2 if and only if M is a nearly Kaehler 6-manifold. A central problem in the field has been the absence of any complete inhomogeneous examples. We prove the existence of the first complete inhomogeneous nearly Kaehler 6-manifolds by proving the existence of at least one cohomogeneity one nearly Kaehler structure on the 6-sphere and on the product of a pair of 3-spheres. We conjecture that these are the only simply connected (inhomogeneous) cohomogeneity one nearly Kaehler structures in six dimensions.


title: Improved spectral stability in spin transfer nano-oscillators: single vortex versus coupled vortices dynamics

abstract: We perform a comparative study of spin transfer induced excitation of the gyrotropic motion of a vortex core with either uniform or vortex spin polarizers. The microwave output voltage associated with the vortex dynamics, detected in both cases, displays a strong reduction of phase fluctuations in the case of the vortex polarizer, with a decrease of the peak linewidth by one order of magnitude down to 200kHz at zero field. A thorough study of rf emission features for the different accessible vortex configurations shows that this improvement is related to the excitation of coupled vortex dynamics by spin transfer torques.


title: Local Single Ring Theorem

abstract: The Single Ring Theorem, by Guionnet, Krishnapur and Zeitouni, describes the empirical eigenvalue distribution of a large generic matrix with prescribed singular values, i.e. an $N\times N$ matrix of the form $A=UTV$, with $U, V$ some independent Haar-distributed unitary matrices and $T$ a deterministic matrix whose singular values are the ones prescribed. In this text, we give a local version of this result, proving that it remains true at the microscopic scale $(\log N)^{-1/4}$. On our way to prove it, we prove a matrix subordination result for singular values of sums of non Hermitian matrices, as Kargin did for Hermitian matrices. This allows to prove a local law for the singular values of the sum of two non Hermitian matrices and a delocalization result for singular vectors.


title: Non-equilibrium dynamics of the Dicke model for mesoscopic aggregates: signatures of superradiance

abstract: In the dissipative quantum dynamics of a mesoscopic aggregate of excited two level systems (atoms) coupled to a single resonance mode of a cavity, two physical phenomena associated with superradiance appear. A pronounced emission peak on short time scales is related to the known superradiant burst of excited atoms cooperatively radiating into free space. It is followed by relaxation to a stationary state of the composite system such that for sufficiently large atom-field coupling a strongly correlated state emerges. The crossover to this state can be interpreted as the precursor of the transition from normal to superradiant phase in the thermodynamic limit of the Dicke Hamiltonian. Motivated by recent experimental activities, these features are investigated in detail for a mesoscopic number of atoms and a cavity embedded in a dissipative medium described by a damped Dicke model. We identify observables and characteristics of the quantum dynamics on shorter time scales which allow to clearly distinguish weakly correlated from strongly correlated Dicke physics.


title: Observations on the Partial Breaking of $N=2$ Rigid Supersymmetry

abstract: We study the partial breaking of $N=2$ rigid supersymmetry for a generic rigid special geometry of $n$ abelian vector multiplets in the presence of Fayet-Iliopoulos terms induced by the Hyper-K\"ahler momentum map. By exhibiting the symplectic structure of the problem we give invariant conditions for the breaking to occur, which rely on a quartic invariant of the Fayet-Iliopoulos charges as well as on a modification of the $N=2$ rigid symmetry algebra by a vector central charge.


title: ANTARES Constrains a Blazar Origin of Two IceCube PeV Neutrino Events

abstract: The source(s) of the neutrino excess reported by the IceCube Collaboration is unknown. The TANAMI Collaboration recently reported on the multiwavelength emission of six bright, variable blazars which are positionally coincident with two of the most energetic IceCube events. Objects like these are prime candidates to be the source of the highest-energy cosmic rays, and thus of associated neutrino emission. We present an analysis of neutrino emission from the six blazars using observations with the ANTARES neutrino telescope.The standard methods of the ANTARES candidate list search are applied to six years of data to search for an excess of muons --- and hence their neutrino progenitors --- from the directions of the six blazars described by the TANAMI Collaboration, and which are possibly associated with two IceCube events. Monte Carlo simulations of the detector response to both signal and background particle fluxes are used to estimate the sensitivity of this analysis for different possible source neutrino spectra. A maximum-likelihood approach, using the reconstructed energies and arrival directions of through-going muons, is used to identify events with properties consistent with a blazar origin.Both blazars predicted to be the most neutrino-bright in the TANAMI sample (1653$-$329 and 1714$-$336) have a signal flux fitted by the likelihood analysis corresponding to approximately one event. This observation is consistent with the blazar-origin hypothesis of the IceCube event IC14 for a broad range of blazar spectra, although an atmospheric origin cannot be excluded. No ANTARES events are observed from any of the other four blazars, including the three associated with IceCube event IC20. This excludes at a 90\% confidence level the possibility that this event was produced by these blazars unless the neutrino spectrum is flatter than $-2.4$.


title: A Proximal Bregman Projection Approach to Continuous Max-Flow Problems Using Entropic Distances

abstract: One issue limiting the adaption of large-scale multi-region segmentation is the sometimes prohibitive memory requirements. This is especially troubling considering advances in massively parallel computing and commercial graphics processing units because of their already limited memory compared to the current random access memory used in more traditional computation. To address this issue in the field of continuous max-flow segmentation, we have developed a \textit{pseudo-flow} framework using the theory of Bregman proximal projections and entropic distances which implicitly represents flow variables between labels and designated source and sink nodes. This reduces the memory requirements for max-flow segmentation by approximately 20\% for Potts models and approximately 30\% for hierarchical max-flow (HMF) and directed acyclic graph max-flow (DAGMF) models. This represents a great improvement in the state-of-the-art in max-flow segmentation, allowing for much larger problems to be addressed and accelerated using commercially available graphics processing hardware.


title: A sharp quantitative version of Alexandrov's theorem via the method of moving planes

abstract: We prove the following quantitative version of the celebrated Soap Bubble Theorem of Alexandrov. Let $S$ be a $C^2$ closed embedded hypersurface of $\mathbb{R}^{n+1}$, $n\geq1$, and denote by $osc(H)$ the oscillation of its mean curvature. We prove that there exists a positive $\varepsilon$, depending on $n$ and upper bounds on the area and the $C^2$-regularity of $S$, such that if $osc(H) \leq \varepsilon$ then there exist two concentric balls $B_{r_i}$ and $B_{r_e}$ such that $S \subset \overline{B}_{r_e} \setminus B_{r_i}$ and $r_e -r_i \leq C \, osc(H)$, with $C$ depending only on $n$ and upper bounds on the surface area of $S$ and the $C^2$ regularity of $S$. Our approach is based on a quantitative study of the method of moving planes and the quantitative estimate on $r_e-r_i$ we obtain is optimal. As a consequence of this theorem, we also prove that if $osc(H)$ is small then $S$ is diffeomorphic to a sphere and give a quantitative bound which implies that $S$ is $C^1$-close to a sphere.


title: An injectivity radius estimate in terms of metric sphere

abstract: In this paper we prove that if a point $p$ in a complete Riemannian manifold is not a cut point of any point whose distance to $p$ is $r$, then the injectivity radius of $p$ is strictly large than $r$. As a corollary we give a positive answer to a problem raised by Z. Sun and J. Wan.


title: Electronic Medication Prescribing Support System for Diagnosing Tropical Diseases

abstract: This paper presents the development of an e-prescription system for diagnosing tropical diseases.Results after testing the developed system by medical experts indicated that the e-prescription systems is more efficient and less susceptible to common errors associated with the conventional handwritten medical prescription and can also go a long way to help to improve patients health outcome in the health industry especially in the tropics.


title: Improved gradient flow for step scaling function and scale setting

abstract: The gradient flow renormalized coupling offers a simple and relatively inexpensive way to calculate the step scaling function and the lattice scale, but both applications can be hindered by large lattice artifacts. Recently we introduced an empirical non-perturbative improvement that can reduce, even remove $\mathcal{O}(a^2)$ lattice artifacts. The method is easy to implement and can be applied to any lattice gauge theory of interest both in step scaling studies and for scale setting. In this talk I will briefly review this improvement method and discuss its application for determining the discrete $\beta$ function of the 8 and 12 flavor SU(3) systems and for improved scale setting in 2+1+1 flavor QCD


title: Thermodynamics of the Bose-Hubbard model in a Bogoliubov+U theory

abstract: We derive the Bogoliubov+U formalism to study the thermodynamical properties of the Bose-Hubbard model. The framework can be viewed as the zero-frequency limit of bosonic dynamical mean-field theory (B-DMFT), but equally well as an extension of the mean-field decoupling approximation in which pair creation and annihilation of depleted particles is taken into account. The self-energy on the impurity site is treated variationally, minimizing the grand potential. The theory containing just three parameters that are determined self-consistently reproduces the T=0 phase diagrams of the three-dimensional and two-dimensional Bose-Hubbard model with an accuracy of 1% or better. The superfluid to normal transition at finite temperature is also reproduced well and only slightly less accurately than in B-DMFT.


title: Estimating viral infection parameters using Markov Chain Monte Carlo simulations

abstract: Given a mathematical model quantifying the viral infection of pandemic influenza H1N1pdm09-H275 wild type (WT) and H1N1pdm09-H275Y mutant (MUT) strains, we describe a simple method of estimating the model's constant parameters using Monte Carlo methods. Monte Carlo parameter estimation methods present certain advantages over the bootstrapping methods previously used in such studies: the result comprises actual parameter distributions (posteriors) that can be used to compare different viral strains; the recovered parameter distributions offer an exact method to compute credible intervals (similar to the frequentist 95% parametric confidence intervals (CI)), that, in turn, using a suitable analysis statistic, will be narrower than the ones obtained from bootstrapping; given an appropriate computational parallelization, Monte Carlo methods are also faster and less computationally intensive than bootstrapping. We fit Gaussian distributions to the parameter posterior distributions and use a two-sided Kolmogorov-Smirnoff test to compare the two strains from a parametric point of view; our example result shows that the two strains are 94% different. Furthermore, based on the obtained parameter values, we estimate the reproductive number R0 for each strain and show that the infectivity of the mutant strain is larger than the wild type strain.


title: On the analytic torsion of hyperbolic manifolds of finite volume

abstract: In this paper we study the analytic torsion for a complete oriented hyperbolic manifold of finite volume. This requires the definition of a regularized trace of heat operators. We use the Selberg trace formula to study the asymptotic behavior of the regularized trace for small time. The main result of the paper is a new approach to deal with the weighted orbital integrals on the geometric side of the trace formula.


title: Alignment and resolution studies of a MARS scanner

abstract: The MARS scanner is designed for the x-ray spectroscopic study of samples with the aid of computer tomography methods. Computer tomography allows the reconstruction of slices of an investigated sample using a set of shadow projections obtained for different angles. Projections in the MARS scanner are produced using a cone x-ray beam geometry. Correct reconstruction in this scheme requires precise knowledge of several geometrical parameters of a tomograph, such as displacement of a rotation axis, x-ray source position with respect to a camera, and camera inclinations. Use of inaccurate parameters leads to a poor sample reconstruction. Non-ideal positioning of camera, x-ray source and cylindrical rotating frame (gantry) itself on which these parts are located, leads to the need for tomograph alignment. In this note we describe the alignment procedure that was used to get different geometrical corrections for the reconstruction. Also, several different estimations of the final spatial resolution for reconstructed images are presented.


title: Longitudinal wave function control in single quantum dots with an applied magnetic field

abstract: Controlling single-particle wave functions in single semiconductor quantum dots is in demand to implement solid-state quantum information processing and spintronics. Normally, particle wave functions can be tuned transversely by an perpendicular magnetic field. We report a longitudinal wave function control in single quantum dots with a magnetic field. For a pure InAs quantum dot with a shape of pyramid or truncated pyramid, the hole wave function always occupies the base because of the less confinement at base, which induces a permanent dipole oriented from base to apex. With applying magnetic field along the base-apex direction, the hole wave function shrinks in the base plane. Because of the linear changing of the confinement for hole wave function from base to apex, the center of effective mass moves up during shrinking process. Due to the uniform confine potential for electrons, the center of effective mass of electrons does not move much, which results in a permanent dipole moment change and an inverted electron-hole alignment along the magnetic field direction. Manipulating the wave function longitudinally not only provides an alternative way to control the charge distribution with magnetic field but also a new method to tune electron-hole interaction in single quantum dots.


title: Kinetic Monte Carlo study of the type1/type 2 choice in apoptosis elucidates selective killing of cancer cells under death ligand induction

abstract: Death ligand mediated apoptotic activation is a mode of programmed cell death that is widely used in cellular and physiological situations. Interest in studying death ligand induced apoptosis has increased due to the promising role of recombinant soluble forms of death ligands (mainly recombinant TRAIL) in anti-cancer therapy. A clear elucidation of how death ligands activate the type 1 and type 2 apoptotic pathways in healthy and cancer cells may help develop better chemotherapeutic strategies. In this work, we use kinetic Monte Carlo simulations to address the problem of type 1/ type 2 choice in death ligand mediated apoptosis of cancer cells. Our study provides insights into the activation of membrane proximal death module that results from complex interplay between death and decoy receptors. Relative abundance of death and decoy receptors was shown to be a key parameter for activation of the initiator caspases in the membrane module. Increased concentration of death ligands frequently increased the type 1 activation fraction in cancer cells, and, in certain cases changes the signaling phenotype from type 2 to type 1. Results of this study also indicate that inherent differences between cancer and healthy cells, such as in the membrane module, may allow robust activation of cancer cell apoptosis by death ligand induction. At the same time, large cell-to-cell variability through the type 2 pathway was shown to provide protection for healthy cells. Such elucidation of selective activation of apoptosis in cancer cells addresses a key question in cancer biology and cancer therapy.


title: Contact Geometry of the Pontryagin Maximum Principle

abstract: This paper gives a brief contact-geometric account of the Pontryagin maximum principle. We show that key notions in the Pontryagin maximum principle---such as the separating hyperplanes, costate, necessary condition, and normal/abnormal minimizers---have natural contact-geometric interpretations. We then exploit the contact-geometric formulation to give a simple derivation of the transversality condition for optimal control with terminal cost.


title: Busting Up Binaries: Encounters Between Compact Binaries and a Supermassive Black Hole

abstract: Given the stellar density near the galactic center, close encounters between compact object binaries and the supermassive black hole are a plausible occurrence. We present results from a numerical study of close to 13 million such encounters. Consistent with previous studies, we corroborate that, for binary systems tidally disrupted by the black hole, the component of the binary remaining bound to the hole has eccentricity ~ 0.97 and circularizes dramatically by the time it enters the classical LISA band. Our results also show that the population of surviving binaries merits attention. These binary systems experience perturbations to their internal orbital parameters with potentially interesting observational consequences. We investigated the regions of parameter space for survival and estimated the distribution of orbital parameters post-encounter. We found that surviving binaries harden and their eccentricity increases, thus accelerating their merger due gravitational radiation emission and increasing the predicted merger rates by up to 1%.


title: Predicting Results of the Research Excellence Framework using Departmental h-Index -- Revisited

abstract: We revisit our recent study [Predicting results of the Research Excellence Framework using departmental h-index, Scientometrics, 2014, 1-16; arXiv:1411.1996] in which we attempted to predict outcomes of the UK's Research Excellence Framework (REF~2014) using the so-called departmental $h$-index. Here we report that our predictions failed to anticipate with any accuracy either overall REF outcomes or movements of individual institutions in the rankings relative to their positions in the previous Research Assessment Exercise (RAE~2008).


title: Testing for Structural Breaks via Ordinal Pattern Dependence

abstract: We propose new concepts in order to analyze and model the dependence structure between two time series. Our methods rely exclusively on the order structure of the data points. Hence, the methods are stable under monotone transformations of the time series and robust against small perturbations or measurement errors. Ordinal pattern dependence can be characterized by four parameters. We propose estimators for these parameters, and we calculate their asymptotic distributions. Furthermore, we derive a test for structural breaks within the dependence structure. All results are supplemented by simulation studies and empirical examples. For three consecutive data points attaining different values, there are six possibilities how their values can be ordered. These possibilities are called ordinal patterns. Our first idea is simply to count the number of coincidences of patterns in both time series, and to compare this with the expected number in the case of independence. If we detect a lot of coincident patterns, this means that the up-and-down behavior is similar. Hence, our concept can be seen as a way to measure non-linear `correlation'. We show in the last section, how to generalize the concept in order to capture various other kinds of dependence.