title: Lattice study of the Boer-Mulders transverse momentum distribution in the pion

abstract: The Boer-Mulders transverse momentum-dependent parton distribution (TMD) characterizes polarized quark transverse momentum in an unpolarized hadron. Techniques previously developed for lattice calculations of nucleon TMDs are applied to the pion. These techniques are based on the evaluation of matrix elements of quark bilocal operators containing a staple-shaped Wilson connection. Results for the Boer-Mulders transverse momentum shift in the pion, obtained at a pion mass of $m_{\pi} = 518\, \mbox{MeV} $, are presented and compared to corresponding results in the nucleon.


title: Extraction of 3D field maps of magnetic multipoles from 2D surface measurements with applications to the optics calculations of the large-acceptance superconducting fragment separator BigRIPS

abstract: The fringing fields of magnets with large apertures and short lengths greatly affect ion-optical calculations. In particular, for a high magnetic field where the iron core becomes saturated, the effective lengths and shapes of the field distribution must be considered because they change with the excitation current. Precise measurement of the three-dimensional magnetic fields and the correct application of parameters in the ion-optical calculations are necessary. First we present a practical numerical method of extracting full 3D magnetic field maps of magnetic multipoles from 2D field measurements of the surface of a cylinder. Using this novel method we extracted the distributions along the beam axis for the coefficient of the first-order quadrupole component, which is the leading term of the quadrupole components in the multipole expansion of magnetic fields and proportional to the distance from the axis. Higher order components of the 3D magnetic field can be extracted from the leading term via recursion relations. The measurements were done for many excitation current values for the large-aperture superconducting triplet quadrupole magnets (STQs) in the BigRIPS fragment separator at the RIKEN Nishina Center RI Beam Factory. These distributions were parameterized using the Enge functions to fit the fringe field shapes at all excitation current values, so that unmeasured values are interpolated. The extracted distributions depend only on the position along the beam axis, and thus the measured three-dimensional field can easily be parameterized for ion-optical calculations. We implemented these parameters in the ion-optical calculation code COSY INFINITY and realized a first-order calculation that incorporates the effect of large and varying fringe fields more accurately. We applied the calculation to determine the excitation current settings of the STQs...


title: Exotic open-flavor $bc\bar{q}\bar{q}$, $bc\bar{s}\bar{s}$ and $qc\bar{q}\bar{b}$, $sc\bar{s}\bar{b}$ tetraquark states

abstract: We study the exotic $bc\bar{q}\bar{q}$, $bc\bar{s}\bar{s}$ and $qc\bar{q}\bar{b}$, $sc\bar{s}\bar{b}$ systems by constructing the corresponding tetraquark currents with $J^P=0^+$ and $1^+$. After investigating the two-point correlation functions and the spectral densities, we perform QCD sum rule analysis and extract the masses of these open-flavor tetraquark states. Our results indicate that the masses of both the scalar and axial vector tetraquark states are about $7.1-7.2$ GeV for the $bc\bar{q}\bar{q}$ system and $7.2-7.3$ GeV for the $bc\bar{s}\bar{s}$ system. For the $qc\bar{q}\bar{b}$ tetraquark states with $J^P=0^+$ and $1^+$, their masses are extracted to be around $7.1$ GeV. The masses for the scalar and axial vector $sc\bar{s}\bar{b}$ states are $7.1$ GeV and $6.9-7.1$ GeV, respectively. The tetraquark states $qc\bar{q}\bar{b}$ and $sc\bar{s}\bar{b}$ lie below the thresholds of $D^{(\ast)}B^{(\ast)}$ and $D_s^{(\ast)}B_s^{(\ast)}$ respectively, but they can decay into $B_c$ plus a light meson. However, the tetraquark states $bc\bar{q}\bar{q}$ and $bc\bar{s}\bar{s}$ lie below the $D^{(\ast)}\bar B^{(\ast)}$ and $D_s^{(\ast)}\bar B_s^{(\ast)}$ thresholds, suggesting dominantly weak decay mechanisms.


title: Complex box bounds for real maps

abstract: In this paper we prove complex bounds, also referred to as a priori bounds, for real analytic (and even C3) interval maps. This means that we associate to such a map a complex box mapping (which provides a kind of Markov structure), together with universal geometric bounds on the shape of the domains. Such bounds show that the first return maps to these domains are well-controlled, and consequently form one of the corner stones in many recent results on one-dimensional dynamics: renormalisation theory, rigidity results, density of hyperbolicity, local connectivity.


title: Smoothed and Iterated Bootstrap Confidence Regions for Parameter Vectors

abstract: The construction of confidence regions for parameter vectors is a difficult problem in the nonparametric setting, particularly when the sample size is not large. The bootstrap has shown promise in solving this problem, but empirical evidence often indicates that some bootstrap methods have difficulty in maintaining the correct coverage probability, while other methods may be unstable, often resulting in very large confidence regions. One way to improve the performance of a bootstrap confidence region is to restrict the shape of the region in such a way that the error term of an expansion is as small an order as possible. To some extent, this can be achieved by using the bootstrap to construct an ellipsoidal confidence region. This paper studies the effect of using the smoothed and iterated bootstrap methods to construct an ellipsoidal confidence region for a parameter vector. The smoothed estimate is based on a multivariate kernel density estimator. This paper establishes a bandwidth matrix for the smoothed bootstrap procedure that reduces the asymptotic coverage error of the bootstrap percentile method ellipsoidal confidence region. We also provide an analytical adjustment to the nominal level to reduce the computational cost of the iterated bootstrap method. Simulations demonstrate that the methods can be successfully applied in practice.


title: A note on singularity of fibers of singular sets

abstract: A general theorem on fibers of singular sets is presented.


title: Unexpected links reflect the noise in networks

abstract: Gene covariation networks are commonly used to study biological processes. The inference of gene covariation networks from observational data can be challenging, especially considering the large number of players involved and the small number of biological replicates available for analysis. We propose a new statistical method for estimating the number of erroneous edges in reconstructed networks that strongly enhances commonly used inference approaches. This method is based on a special relationship between sign of correlation (positive/negative) and directionality (up/down) of gene regulation, and allows for the identification and removal of approximately half of all erroneous edges. Using the mathematical model of Bayesian networks and positive correlation inequalities we establish a mathematical foundation for our method. Analyzing existing biological datasets, we find a strong correlation between the results of our method and false discovery rate (FDR). Furthermore, simulation analysis demonstrates that our method provides a more accurate estimate of network error than FDR.


title: On the Optimum Energy Efficiency for Flat-fading Channels with Rate-dependent Circuit Power

abstract: This paper investigates the optimum energy efficiency (EE) and the corresponding spectral efficiency (SE) for a communication link operating over a flat-fading channel. The EE is evaluated by the total energy consumption for transmitting per message bit. Three channel cases are considered, namely static channel with channel state information available at transmitter (CSIT), fast-varying (FV) channel with channel distribution information available at transmitter (CDIT), and FV channel with CSIT. A general circuit power model is considered. For all the three channel cases, the tradeoff between the EE and SE is studied. It is shown that the EE improves strictly as the SE increases from 0 to the optimum SE, and then strictly degrades as the SE increases beyond the optimum SE. The impact of {\kappa}, {\rho} and other system parameters on the optimum EE and corresponding SE is investigated to obtain insight.Some of the important and interesting results for all the channel cases include: (1) when {\kappa} increases the SE corresponding to the optimum EE should keep unchanged if {\phi}(R) = R, but reduced if {\phi}(R) is strictly convex of R; (2) when the rate-independent circuit power {\rho} increases, the SE corresponding to the optimum EE has to be increased. A polynomial-complexity algorithm is developed with the bisection method to find the optimum SE. The insight is corroborated and the optimum EE for the three cases are compared by simulation results.


title: Matter-wave interference with particles selected from a molecular library with masses exceeding 10000 amu

abstract: The quantum superposition principle, a key distinction between quantum physics and classical mechanics, is often perceived as a philosophical challenge to our concepts of reality, locality or space-time since it contrasts our intuitive expectations with experimental observations on isolated quantum systems. While we are used to associating the notion of localization with massive bodies, quantum physics teaches us that every individual object is associated with a wave function that may eventually delocalize by far more than the body's own extension. Numerous experiments have verified this concept at the microscopic scale but intuition wavers when it comes to delocalization experiments with complex objects. While quantum science is the uncontested ideal of a physics theory, one may ask if the superposition principle can persist on all complexity scales. This motivates matter-wave diffraction and interference studies with large compounds in a three-grating interferometer configuration which also necessitates the preparation of high-mass nanoparticle beams at low velocities. Here we demonstrate how synthetic chemistry allows us to prepare libraries of fluorous porphyrins which can be tailored to exhibit high mass, good thermal stability and relatively low polarizability, which allows us to form slow thermal beams of these high-mass compounds, which can be detected in electron ionization mass spectrometry. We present successful superposition experiments with selected species from these molecular libraries in a quantum interferometer, which utilizes the diffraction of matter waves at an optical phase grating. We observe high-contrast quantum fringe patterns with molecules exceeding a mass of 10 000 amu and 810 atoms in a single particle.


title: Absolutely calibrated radio polarimetry of the inner Galaxy at 2.3 GHz and 4.8 GHz

abstract: We present high sensitivity and absolutely calibrated images of diffuse radio polarisation at a resolution of about 10 arcmin covering the range 10 degr <l<34 degr and |b|<5 degr at 2.3 GHz from the S-band Parkes All Sky Survey and at 4.8 GHz from the Sino-German 6 cm polarisation survey of the Galactic plane. Strong depolarisation near the Galactic plane is seen at 2.3 GHz, which correlates with strong Halpha emission. We ascribe the depolarisation to spatial Faraday rotation measure fluctuations of about 65 rad/m2 on scales smaller than 6-9 pc. We argue that most (about 90%) of the polarised emission seen at 4.8 GHz originates from a distance of 3-4 kpc in the Scutum arm and that the random magnetic field dominates the regular field there. A branch extending from the North Polar Spur towards lower latitudes can be identified from the polarisation image at 4.8 GHz but only partly from the polarised image at 2.3 GHz, implying the branch is at a distance larger than 2-3 kpc. We show that comparison of structure functions of complex polarised intensity with those of polarised intensity can indicate whether the observed polarised structures are intrinsic or caused by Faraday screens. The probability distribution function of gradients from the polarisation images at 2.3 GHz indicates the turbulence in the warm ionised medium is transonic.


title: Entanglement entropy on the fuzzy sphere

abstract: We obtain entanglement entropy on the noncommutative (fuzzy) two-sphere. To define a subregion with a well defined boundary in this geometry, we use the symbol map between elements of the noncommutative algebra and functions on the sphere. We find that entanglement entropy is not proportional to the length of the region's boundary. Rather, in agreement with holographic predictions, it is extensive for regions whose area is a small (but fixed) fraction of the total area of the sphere. This is true even in the limit of small noncommutativity. We also find that entanglement entropy grows linearly with N, where N is the size of the irreducible representation of SU(2) used to define the fuzzy sphere.


title: Orbital control of effective dimensionality: from spin-orbital fractionalization to confinement in the anisotropic ladder system CaCu2O3

abstract: Fractionalization of an electronic quasiparticle into spin, charge and orbital parts is a fundamental and characteristic property of interacting electrons in one dimension. However, real materials are never strictly one-dimensional and the fractionalization phenomena are hard to observe. Here we studied the spin and orbital excitations of the anisotropic ladder material CaCu2O3, whose electronic structure is not one-dimensional. Combining high-resolution resonant inelastic x-ray scattering experiments with theoretical model calculations we show that: (i) spin-orbital fractionalization occurs in CaCu2O3 along the leg direction x through the xz orbital channel as in a 1D system; and (ii) no fractionalization is observed for the xy orbital, which extends in both leg and rung direction, contrary to a 1D system. We conclude that the directional character of the orbital hopping can select different degrees of dimensionality. Using additional model calculations, we show that spin-orbital separation is generally far more robust than the spin-charge separation. This is not only due to the already mentioned selection realized by the orbital hopping, but also due to the fact that spinons are faster than the orbitons.


title: Quantum Imaging of High-Dimensional Hilbert Spaces with Radon Transform

abstract: High-dimensional Hilbert spaces possess large information encoding and transmission capabilities. Characterizing exactly the real potential of high-dimensional entangled systems is a cornerstone of tomography and quantum imaging. The accuracy of the measurement apparatus and devices used in quantum imaging is physically limited, which allows no further improvements to be made. To extend the possibilities, we introduce a post-processing method for quantum imaging that is based on the Radon transform and the projection-slice theorem. The proposed solution leads to an enhanced precision and a deeper parameterization of the information conveying capabilities of high-dimensional Hilbert spaces. We demonstrate the method for the analysis of high-dimensional position-momentum photonic entanglement. We show that the entropic separability bound in terms of standard deviations is violated considerably more strongly in comparison to the standard setting and current data processing. The results indicate that the possibilities of the quantum imaging of high-dimensional Hilbert spaces can be extended by applying appropriate calculations in the post-processing phase.


title: Lifshitz and Schrodinger Vacua, Superstar Resolution in Gauged Maximal Supergravities

abstract: We consider the subset of gauged maximal supergravities that consists of the SO(n+1) gauge fields A^{ij} and the scalar deformation T^{ij} of the S^n in the spherical reduction of M-theory or type IIB. We focus on the Abelian Cartan subgroup and the diagonal entries of T^{ij}. The resulting theories can be viewed as the STU models with additional hyperscalars. We find that the theories with only one or two such vectors can be generalized naturally to arbitrary dimensions. The same is true for the D=4 or 5 Einstein-Maxwell theory with such a hyperscalar. The gauge fields become massive, determined by stationary points of the hyperscalars a la the analogous Abelian Higgs mechanism. We obtain classes of Lifshitz and Schrodinger vacua in these theories. The scaling exponent z turns out to be rather restricted, taking fractional or irrational numbers. Tweaking the theories by relaxing the mass parameter or making a small change of the superpotential, we find that solutions with z=2 can emerge. In a different application, we find that the resolution of superstar singularity in the STU models by using bubbling-AdS solitons can be generalized to arbitrary dimensions in our theories. In particular, we obtain the smooth AdS solitons that can be viewed as the resolution of the Reissner-Nordstrom superstars in general dimensions.


title: Thermal machines beyond the weak coupling regime

abstract: How much work can be extracted from a heat bath using a thermal machine? The study of this question has a very long tradition in statistical physics in the weak-coupling limit, applied to macroscopic systems. However, the assumption that thermal heat baths remain uncorrelated with physical systems at hand is less reasonable on the nano-scale and in the quantum setting. In this work, we establish a framework of work extraction in the presence of quantum correlations. We show in a mathematically rigorous and quantitative fashion that quantum correlations and entanglement emerge as a limitation to work extraction compared to what would be allowed by the second law of thermodynamics. At the heart of the approach are operations that capture naturally non-equilibrium dynamics encountered when putting physical systems into contact with each other. We discuss various limits that relate to known results and put our work into context of approaches to finite-time quantum thermodynamics.


title: Billiards on Graphs

abstract: We prove some partial results on the periodicity of billiard systems on graphs. The results specialize to the case of $n$ billiards with equal mass on the unit interval or circle traveling at the same speed.


title: Identification and Separation of Radioactive Isotope Beams by the BigRIPS Separator at the RIKEN RI Beam Factory

abstract: We have developed a method for achieving excellent resolving power in in-flight particle identification of radioactive isotope (RI) beams at the BigRIPS fragment separator at the RIKEN Nishina Center RI Beam Factory (RIBF). In the BigRIPS separator, RI beams are identified by their atomic number Z and mass-to-charge ratio A/Q which are deduced from the measurements of time of flight (TOF), magnetic rigidity (Brho) and energy loss (delta-E), and delivered as tagged RI beams to a variety of experiments including secondary reaction measurements. High A/Q resolution is an essential requirement for this scheme, because the charge state Q of RI beams has to be identified at RIBF energies such as 200-300 MeV/nucleon. By precisely determining the Brho and TOF values, we have achieved relative A/Q resolution as good as 0.034% (root-mean-square value). The achieved A/Q resolution is high enough to clearly identify the charge state Q in the Z versus A/Q particle identification plot, where fully-stripped and hydrogen-like peaks are very closely located. The precise Brho determination is achieved by refined particle trajectory reconstruction, while a slew correction is performed to precisely determine the TOF value. Furthermore background events are thoroughly removed to improve reliability of the particle identification. In the present paper we present the details of the particle identification scheme in the BigRIPS separator. The isotope separation in the BigRIPS separator is also briefly introduced.


title: Self-quartic interaction for a scalar field in an extended DFR noncommutative spacetime

abstract: The framework Doplicher-Fredenhagen-Roberts (DFR) of a noncommutative (NC) space-time is considered as a alternative approach to study the NC space-time of the early Universe. In this formalism, the parameter of noncommutative $\theta^{\mu\nu}$ is promoted to a coordinate of the space-time, and consequently, we are describing a field theory in a space-time with extra-dimension. Consequently, there is a canonical momentum associated to this new coordinate in which the effects of a new physics can emerge in the propagation of the fields along the extra-dimension. The Fourier space of this framework is automatically extended by the addition of new momenta components. The main concept that we would like to emphasize from the outset is that the formalism demonstrated here will not be constructed introducing a NC parameter in the system, as usual. It will be generated naturally from an already NC space. When the components of the new momentum are zero, the DFR approach is reduced to the usual NC case, in which $\theta^{\mu\nu}$ is a antisymmetric constant matrix. We study a scalar field action with self-quartic interaction $\phi^{4}\star$ defined in the DFR NC spacetime, obtaining the Feynman rules in the Fourier space for the scalar propagator and vertex of the model. With these rules we are able to build out the radiative corrections to one loop order for the model propagator. The influence of the NC scale, as well as the propagation of the field in the extra-dimension, are analyzed in the ultraviolet divergences scenario. We investigate the actual possibility if this $\theta^{\mu\nu}$ conjugate momentum has the property of healing the mixing IR/UV divergences that emerges in this recently new NC spacetime quantum field theory.


title: Stochastically Symplectic Maps and Their Applications to Navier-Stokes Equation

abstract: Poincare's invariance principle for Hamiltonian flows implies Kelvin's principle for solution to Incompressible Euler Equation. Iyer-Constantin Circulation Theorem offers a stochastic analog of Kelvin's principle for Navier-Stokes Equation. Weakly symplectic diffusions are defined to produce stochastically symplectic flows in a systematic way. With the aid of symplectic diffusions, we produce a family of martigales associated with solutions to Navier-Stokes Equation that in turn can be used to prove Iyer-Constantin Circulation Theorem. We also review some basic facts in symplectic and contact geometry and their applications to Euler Equation.


title: Thermalization of field driven quantum systems

abstract: There is much interest in how quantum systems thermalize after a sudden change, because unitary evolution should preclude thermalization. The eigenstate thermalization hypothesis resolves this because all observables for quantum states in a small energy window have essentially the same value; it is violated for integrable systems due to the infinite number of conserved quantities. Here, we show that when a system is driven by a DC electric field there are five generic behaviors: (i) monotonic or (ii) oscillatory approach to an infinite-temperature steady state; (iii) monotonic or (iv) oscillatory approach to a nonthermal steady state; or (v) evolution to an oscillatory state. Examining the Hubbard model (which thermalizes under a quench) and the Falicov-Kimball model (which does not), we find both exhibit scenarios (i-iv), while only Hubbard shows scenario (v). This shows richer behavior than in interaction quenches and integrability in the absence of a field plays no role.


title: New Evidence for Efficient Collisionless Heating of Electrons at the Reverse Shock of a Young Supernova Remnant

abstract: Although collisionless shocks are ubiquitous in astrophysics, certain key aspects of them are not well understood. In particular, the process known as collisionless electron heating, whereby electrons are rapidly energized at the shock front, is one of the main open issues in shock physics. Here we present the first clear evidence for efficient collisionless electron heating at the reverse shock of Tycho's supernova remnant (SNR), revealed by Fe-K diagnostics using high-quality X-ray data obtained by the Suzaku satellite. We detect K-beta (3p->1s) fluorescence emission from low-ionization Fe ejecta excited by energetic thermal electrons at the reverse shock front, which peaks at a smaller radius than Fe K-alpha (2p->1s) emission dominated by a relatively highly-ionized component. Comparison with our hydrodynamical simulations implies instantaneous electron heating to a temperature 1000 times higher than expected from Coulomb collisions alone. The unique environment of the reverse shock, which is propagating with a high Mach number into rarefied ejecta with a low magnetic field strength, puts strong constraints on the physical mechanism responsible for this heating, and favors a cross-shock potential created by charge deflection at the shock front. Our sensitive observation also reveals that the reverse shock radius of this SNR is about 10% smaller than the previous measurement using the Fe K-alpha morphology from the Chandra observations. Since strong Fe K-beta fluorescence is expected only from low-ionization plasma where Fe ions still have many 3p electrons, this feature is key to diagnosing the plasma state and distribution of the immediate postshock ejecta in a young SNR.


title: The Nature of Millisecond Pulsars with Helium White Dwarf Companions

abstract: We examine the growing data set of binary millisecond pulsars that are thought to have a helium white dwarf companion. These systems are believed to form when a low- to intermediate-mass companion to a neutron star fills its Roche lobe between central hydrogen exhaustion and core helium ignition. We confirm that our own stellar models reproduce a well-defined period-companion mass relation irrespective of the details of the mass transfer process. With magnetic braking this relation extends to periods of less than 1d for a 1Msun giant donor. With this and the measured binary mass functions we calculate the orbital inclination of each system for a given pulsar mass. We expect these inclinations to be randomly oriented in space. If the masses of the pulsars were typically 1.35Msun then there would appear to be a distinct dearth of high-inclination systems. However if the pulsar masses are more typically 1.55 to 1.65Msun then the distribution of inclinations is indeed indistinguishable from random. If it were as much as 1.75Msun then there would appear to be an excess of high-inclination systems. Thus with the available data we can argue that the neutron star masses in binary millisecond pulsars recycled by mass transfer from a red giant typically lie around 1.6Msun and that there is no preferred inclination at which these systems are observed. Hence there is reason to believe that pulsar beams are either sufficiently broad or show no preferred direction relative to the pulsar's spin axis which is aligned with the binary orbit. This is contrary to some previous claims, based on a subset of the data available today, that there might be a tendency for the pulsar beams to be perpendicular to their spin.


title: Point-process deconvolution of fMRI reveals effective connectivity alterations in chronic pain patients

abstract: It is now recognized that important information can be extracted from the brain spontaneous activity, as exposed by recent analysis using a repertoire of computational methods. In this context a novel method, based on a blind deconvolution technique, is used in this paper to analyze potential changes due to chronic pain in the brain pain matrix's effective connectivity. The approach is able to deconvolve the hemodynamic response function to spontaneous neural events, i.e., in the absence of explicit onset timings, and to evaluate information transfer between two regions as a joint probability of the occurrence of such spontaneous events. The method revealed that the chronic pain patients exhibit important changes in the Insula's effective connectivity which can be relevant to understand the overall impact of chronic pain on brain function.


title: Love can be Tough to Measure

abstract: The waveform phase for a neutron star binary can be split into point-particle terms and finite-size terms (characterized by the Love number) that account for equation of state effects. The latter first enter at 5 post-Newtonian (PN) order (i.e. proportional to the tenth power of the orbital velocity), but the former are only known completely to 3.5 PN order, with higher order terms only known to leading-order in the mass-ratio. We here find that not including point-particle terms at 4PN order to leading- and first-order in the mass ratio in the template model can severely deteriorate our ability to measure the equation of state. This problem can be solved if one uses numerical waveforms once their own systematic errors are under control.


title: Cohomological resolutions for anomalous Lie constraints

abstract: It is shown that the BRST resolution of the spaces of physical states of the systems with anomalies can be consistently defined. The appropriate anomalous complexes are obtained by canonical restrictions of the ghost extended spaces to the kernel of anomaly operator without any modifications of the "matter" sector. The cohomologies of the anomalous complex for the case of constarints constituting a centrally extended simple Lie algebra of compact type are calculated and analyzed in details within the framework of Hodge - deRham - Kaehler theory: the vanishing theorem of the relative cohomologies is proved and the absolute cohomologies are reconstructed.


title: A SIS reaction-diffusion-advection model in a low-risk and high-risk domain

abstract: A simplified SIS reaction-diffusion-advection model is proposed and investigated to understand the impact of spatial heterogeneity of environment and advection on the persistence and eradication of an infectious disease. The free boundary is introduced to model the contact transmission at the spreading front of the disease. The behavior of positive solutions to a reaction-diffusion-advection system are discussed. The basic reproduction number $R_0^F(t)$ associated with the diseases in the spatial setting is introduced for this diffusive SIS model with the free boundary, we prove that fast diffusion, small expanding rate and small initial infected domain are benefit for the control of the spatial spread of the disease. Sufficient conditions for the disease to be eradicated or to spread are also given, our result shows that the disease will spread to the whole area if there exists a $t_0\geq 0$ such that $R_{0}^F(t_0)\geq 1$, that is, if the spreading domain is high-risk at some time, the disease will continue to spread till the whole area is infected; while if $R_{0}^F(0)<1$, the disease may be vanishing or keep spreading depends on the expanding rate and the initial number of the infective individuals. The spreading speeds are also given when spreading happens, and numerical simulations are also given to illustrate the impacts of the advection and the expanding rate on the spreading fronts.


title: Higgs Working Group Report of the Snowmass 2013 Community Planning Study

abstract: This report summarizes the work of the Energy Frontier Higgs Boson working group of the 2013 Community Summer Study (Snowmass). We identify the key elements of a precision Higgs physics program and document the physics potential of future experimental facilities as elucidated during the Snowmass study. We study Higgs couplings to gauge boson and fermion pairs, double Higgs production for the Higgs self-coupling, its quantum numbers and $CP$-mixing in Higgs couplings, the Higgs mass and total width, and prospects for direct searches for additional Higgs bosons in extensions of the Standard Model. Our report includes projections of measurement capabilities from detailed studies of the Compact Linear Collider (CLIC), a Gamma-Gamma Collider, the International Linear Collider (ILC), the Large Hadron Collider High-Luminosity Upgrade (HL-LHC), Very Large Hadron Colliders up to 100 TeV (VLHC), a Muon Collider, and a Triple-Large Electron Positron Collider (TLEP).


title: A diagrammatic approach to Kronecker squares

abstract: In this paper we apply a method of Robinson and Taulbee for computing Kronecker coefficients together with other ingredients and show that the multiplicity of each component in a Kronecker square can be obtained from an evaluation of a certain polynomial, which depends only on the component and is computed combinatorially. This polynomial has as many variables as the set of isomorphism classes of connected skew diagrams of size at most the depth of the component. We present two applications. The first is a contribution to Saxl conjecture, which asserts that the Kronecker square of the staircase partition, contains every irreducible character of the symmetric group as a component. We prove that for any partition there is a piecewise polynomial function in one real variable such that for all k, such that the multiplicity of this partition in the Kronecker square of the staircase partition of size k is given by the evaluation of the polynomial function in k. The second application is a proof of a new stability property for Kronecker coefficients.


title: Improving frequency selection of driven pulses using derivative-based transition suppression

abstract: Many techniques in quantum control rely on frequency separation as a means for suppressing unwanted couplings. In its simplest form, the mechanism relies on the low bandwidth of control pulses of long duration. Here we perform a higher-order quantum-mechanical treatment that allows for higher precision and shorter times. In particular, we identify three kinds of off-resonant effects: i) simultaneous unwanted driven couplings (e.g. due to drive crosstalk), ii) additional (initially undriven) transitions such as those in an infinite ladder system, and iii) sideband frequencies of the driving waveform such as we find in corrections to the rotating wave approximation. With a framework that is applicable to all three cases, in addition to the known adiabatic error responsible for a shift of the energy levels we typically see in the spectroscopy of such systems, we derive error terms in a controlled expansion corresponding to higher order adiabatic effects and diabatic excitations. We show, by also expanding the driving waveform in a basis of different order derivatives of a trial function (typically a Gaussian) these different error terms can be corrected for in a systematic way hence strongly improving quantum control of systems with dense spectra.


title: Incentive Mechanisms for Mobile Crowd Sensing: Current States and Challenges of Work

abstract: Mobile crowd sensing (MCS) is a new paradigm which leverages the ubiquity of sensor-equipped mobile devices such as smartphones, music players, and in-vehicle sensors at the edge of the Internet, to collect data. The new paradigm will fuel the evolution of the Internet of Things to three changes as follows: First, the terminal devices at the edge of the Internet change from PCs to mobile phones. Second, the interactive mode extends from the virtual space to the real physical world. Thirdly, the forwarding manner of sensing data are undergoing the transition from the priori to the opportunistic. To better meet the demands of MCS applications at a societal scale, incentive mechanisms are indispensable. In this paper, we will first overview three categories of MCS applications, and then propose a new architecture for MCS applications. Based on the architecture, we discuss various research challenges about incentive mechanism designs for MCS applications, followed by potential future work discussions. Finally, we present potential future works.


title: The Pseudospectrum of the Zakharov-Shabat System

abstract: We study the pseudospectrum of the non-selfadjoint Zakharov-Shabat system in the semiclassical regime. The pseudospectrum may be defined as the union of the spectra of perturbations of the Zakharov-Shabat system, thus it is relevant to the numerical computation of true eigenvalues.


title: On spontaneous breaking of conformal symmetry by probe flavour D-branes

abstract: We explore the possibilities of breaking conformal symmetry spontaneously by introducing flavour branes into conformal holographic backgrounds in the probe limit. A prototype model of such a mechanism is based on placing D7 - anti-D7 configuration in the Klebanov-Witten conifold based model. In this paper we generalize this model. We conjecture on the required topology of the backgrounds and the corresponding probe brane embeddings. We identify several models that obey these requirements and admit spontaneous breaking of conformal invariance. These include type IIB conifold based examples, dual to defect field theories based on the conifold, and type IIA constructions based on the ABJM model. We identify the dilaton, the corresponding Goldstone boson, discuss its effective action and address the "a-term". We briefly discuss the relevance of these models to the pseudo dilaton.


title: Electron acceleration in a nonrelativistic shock with very high Alfv\'en Mach number

abstract: Electron acceleration associated with various plasma kinetic instabilities in a nonrelativistic, very-high-Alfv\'en Mach-number ($M_A \sim 45$) shock is revealed by means of a two-dimensional fully kinetic PIC simulation. Electromagnetic (ion Weibel) and electrostatic (ion-acoustic and Buneman) instabilities are strongly activated at the same time in different regions of the two-dimensional shock structure. Relativistic electrons are quickly produced predominantly by the shock surfing mechanism with the Buneman instability at the leading edge of the foot. The energy spectrum has a high-energy tail exceeding the upstream ion kinetic energy accompanying the main thermal population. This gives a favorable condition for the ion acoustic instability at the shock front, which in turn results in additional energization. The large-amplitude ion Weibel instability generates current sheets in the foot, implying another dissipation mechanism via magnetic reconnection in a three-dimensional shock structure in the very-high-$M_A$ regime.


title: A Non-Centrosymmetric Superconductor with a Bulk 3D Dirac Cone Gapped by Strong Spin Orbit Coupling

abstract: Layered, non-centrosymmetric, heavy element PbTaSe2 is found to be superconducting. We report its electronic properties accompanied by electronic structure calculations. Specific heat, electrical resistivity and magnetic susceptibility measurements indicate that PbTaSe2 is a moderately coupled, type-II BCS superconductor (Tc = 3.72 K, Ginzburg-Landau parameter Kappa = 14) with an electronphonon coupling constant of Lambda_ep = 0.74. Electronic structure calculations reveal a single bulk 3D Dirac cone at the K point of the Brillouin Zone derived exclusively from its hexagonal Pb layer; it is similar to the feature found in graphene except there is a 0.8 eV gap opened by spin-orbit coupling. The combination of large spin-orbit coupling and lack of inversion symmetry also results in large Rashba splitting on the order of tenths of eV.


title: On the inverses of some classes of permutations of finite fields

abstract: We study the compositional inverses of some general classes of permutation polynomials over finite fields. We show that we can write these inverses in terms of the inverses of two other polynomials bijecting subspaces of the finite field, where one of these is a linearized polynomial. In some cases we are able to explicitly obtain these inverses, thus obtaining the compositional inverse of the permutation in question. In addition we show how to compute a linearized polynomial inducing the inverse map over subspaces on which a prescribed linearized polynomial induces a bijection. We also obtain the explicit compositional inverses of two classes of permutation polynomials generalizing those whose compositional inverses were recently obtained in [22] and [24], respectively.


title: Orbital magnetic moment and coercivity of SiO$_{2}$-coated FePt nanoparticles studied by x-ray magnetic circular dichroism

abstract: We have investigated the spin and orbital magnetic moments of Fe in FePt nanoparticles in the $L$1$_{0}$-ordered phase coated with SiO$_{2}$ by x-ray absorption spectroscopy (XAS) and x-ray magnetic circular dichroism (XMCD) measurements at the Fe $L_{\rm 2,3}$ absorption edges. Using XMCD sum rules, we evaluated the ratio of the orbital magnetic moment ($M_{\rm orb}$) to the spin magnetic moment ($M_{\rm spin}$) of Fe to be $M_{\rm orb}/M_{\rm spin}$ = 0.08. This $M_{\rm orb}/M_{\rm spin}$ value is comparable to the value (0.09) obtained for FePt nanoparticles prepared by gas phase condensation, and is larger than the values ($\sim$0.05) obtained for FePt thin films, indicating a high degree of $L$1$_{0}$ order. The hysteretic behavior of the FePt component of the magnetization was measured by XMCD. The magnetic coercivity ($H_{\rm c}$) was found to be as large as 1.8 T at room temperature, $\sim$3 times larger than the thin film value and $\sim$50 times larger than that of the gas phase condensed nanoparticles. The hysteresis curve is well explained by the Stoner-Wohlfarth model for non-interacting single-domain nanoparticles with the $H_{\rm c}$ distributed from 1 T to 5 T.


title: Tensor product weight representations of the Neveu-Schwarz algebra

abstract: In this paper, the tensor product of highest weight modules with intermediate series modules over the Neveu-Schwarz algebra is studied. The weight spaces of such tensor products are all infinitely dimensional if the highest weight module is nontrivial. We find that all such tensor products are indecomposable. We give the necessary and sufficient conditions for these tensor product modules to be irreducible by using shifting technique established for the Virasoro case in [13]. The necessary and sufficient conditions for any two such tensor products to be isomorphic are also determined.


title: Scaling of entanglement entropy in the (branching) multi-scale entanglement renormalization ansatz

abstract: We investigate the scaling of entanglement entropy in both the multi-scale entanglement renormalization ansatz (MERA) and in its generalization, the branching MERA. We provide analytical upper bounds for this scaling, which take the general form of a boundary law with various types of multiplicative corrections, including power-law corrections all the way to a bulk law. For several cases of interest, we also provide numerical results that indicate that these upper bounds are saturated to leading order. In particular we establish that, by a suitable choice of holographic tree, the branching MERA can reproduce the logarithmic multiplicative correction of the boundary law observed in Fermi liquids and spin-Bose metals in $D\geq 2$ dimensions.


title: Lipschitz-Volume rigidity on limit spaces with Ricci curvature bounded from below

abstract: We prove a Lipschitz-Volume rigidity theorem for the non-collapsed Gromov-Hausdorff limits of manifolds with Ricci curvature bounded from below. This is a counterpart of the Lipschitz-Volume rigidity in Alexandrov geometry.


title: Capacity and Delay-Throughput Tradeoff in ICMNs with Poisson Meeting Process

abstract: Intermittently connected mobile networks (ICMNs) serve as an important network model for many critical applications. This paper focuses on a continuous ICMN model where the pair-wise meeting process between network nodes follows a homogeneous and independent Poisson process. This ICMN model is known to serve as good approximations to a class of important ICMNs with mobility models like random waypoint and random direction, so it is widely adopted in the performance study of ICMNs. This paper studies the throughput capacity and delay-throughput tradeoff in the considered ICMNs with Poisson meeting process. For the concerned ICMN, we first derive an exact expression of its throughput capacity based on the pairwise meeting rate therein and analyze the expected end-to-end packet delay under a routing algorithm that can achieve the throughput capacity. We then explore the inherent tradeoff between delay and throughput and establish a necessary condition for such tradeoff that holds under any routing algorithm in the ICMN. To illustrate the applicability of the theoretical results, case studies are further conducted for the random waypoint and random direction mobility models. Finally, simulation and numerical results are provided to verify the efficiency of our theoretical capacity/delay results and to illustrate our findings.


title: BCS-like superconductivity in the noncentrosymmetric compounds Nb$_x$Re$_{1-x}$

abstract: We present research on the superconducting properties of Nb$_{x}$Re$_{1-x}$ ($x$ = 0.13-0.38) obtained by measuring the electrical resistivity $\rho(T)$, magnetic susceptibility $\chi(T)$, specific heat $C_P(T)$, and London penetration depth $\Delta\lambda(T)$. It is found that the superconducting transition temperature $T_c$ decreases monotonically with an increase of $x$. The upper critical field $B_{c2}(T)$ for various $x$ can be nicely scaled by its corresponding $T_c$. The electronic specific heat $C_e(T)/T$, penetration depth $\Delta\lambda(T)$, and superfluid density $\rho_{s}(T)$ demonstrate exponential behavior at low temperatures and can be well fitted by a one-gap BCS model. The residual Sommerfeld coefficient $\gamma_0(B)$ in the superconducting state follows a linear field dependence. All these properties suggest an \emph{s}-wave BCS-type of superconductivity with a very large $B_{c2}(0)$ for Nb$_{x}$Re$_{1-x}$ (0.13 $\leq x \leq$ 0.38).


title: Can we observationally test the weak cosmic censorship conjecture?

abstract: In general relativity, gravitational collapse of matter fields ends with the formation of a spacetime singularity, where the matter density becomes infinite and standard physics breaks down. According to the weak cosmic censorship conjecture, singularities produced in the gravitational collapse cannot be seen by distant observers and must be hidden within black holes. The validity of this conjecture is still controversial and at present we cannot exclude that naked singularities can be created in our Universe from regular initial data. In this paper, we study the radiation emitted by a collapsing cloud of dust and check whether it is possible to distinguish the birth of a black hole from the one of a naked singularity. In our simple dust model, we find that the properties of the radiation emitted in the two scenarios is qualitatively similar. That suggests that observational tests of the cosmic censorship conjecture may be very difficult, even in principle.


title: Event-based Bipartite Consensus on Signed Networks

abstract: In this paper, a bipartite consensus problem for a multi-agent system is formulated firstly. Then an event-based interaction rule is proposed for the multi-agent system with antagonistic interactions. The bipartite consensus stability is analyzed on the basis of spectral properties of the signed Laplacian matrix associated with multi-agent networks. Some simulation results are presented to illustrate the bipartite consensus with the proposed interaction rule.


title: Stanley-Wilf limits are typically exponential

abstract: For a permutation $\pi$, let $S_{n}(\pi)$ be the number of permutations on $n$ letters avoiding $\pi$. Marcus and Tardos proved the celebrated Stanley-Wilf conjecture that $L(\pi)= \lim_{n \to \infty} S_n(\pi)^{1/n}$ exists and is finite. Backed by numerical evidence, it has been conjectured by many researchers over the years that $L(\pi)=\Theta(k^2)$ for every permutation $\pi$ on $k$ letters. We disprove this conjecture, showing that $L(\pi)=2^{k^{\Theta(1)}}$ for almost all permutations $\pi$ on $k$ letters.


title: Oblique propagation of dust ion-acoustic solitary waves in a magnetized dusty pair-ion plasma

abstract: We investigate the propagation characteristics of electrostatic waves in a magnetized pair-ion plasma with immobile charged dusts. It is shown that obliquely propagating (OP) low-frequency (in comparison with the negative-ion cyclotron frequency) long-wavelength "slow" and "fast" modes can propagate, respectively, as dust ion-acoustic (DIA) and dust ion-cyclotron (DIC)-like waves. The properties of these modes are studied with the effects of obliqueness of propagation $(\theta)$, the static magnetic field, the ratios of the negative to positive ion masses $(m)$ and temperatures $(T)$ as well as the dust to negative-ion number density ratio $(\delta)$. Using the standard reductive perturbation technique, we derive a Korteweg-de Vries (KdV) equation which governs the evolution of small-amplitude OP DIA waves. It is found that the KdV equation admits only rarefactive solitons in plasmas with $m$ well below its critical value $m_c~(\gg1)$ which typically depends on $T$ and $\delta$. It is shown that the nonlinear coefficient of the KdV equation vanishes at $m=m_c$, i.e., for plasmas with much heavier negative ions, and the evolution of the DIA waves is then described by a modified KdV (mKdV) equation. The latter is shown to have only compressive soliton solution. The properties of both the KdV and mKdV solitons are studied with the system parameters as above, and possible applications of our results to laboratory and space plasmas are briefly discussed.


title: Interplay of electron correlation and covalency in FeTe(0.6)Se(0.4)

abstract: We investigate the electronic structure of FeTe(0.6)Se(0.4) employing high resolution photoemission spectroscopy and ab initio band structure calculations. Fe 2p core level and the valence band spectra exhibit signature of strong electron correlation in the electronic structure. The electronic states near the Fermi level reduces in intensity with the decrease in temperature in conformity with the insulating transport observed near 300 K. The observation of an insulator to metal transition around 150 K in the transport properties may be related to the spectral lineshape change in the vicinity of the Fermi level observed in this study. The spectral features near Fermi level exhibit significant p orbital character due to the correlation induced Fe d spectral weight transfer. The experimental spectra reveal dominant temperature dependence of the spectral functions possessing large p-character. These results demonstrate significant renormalization of the character of the conduction electrons due to electron correlation and emphasizes the importance of ligand states in the superconductivity of these materials.


title: A Labeling Approach to Incremental Cycle Detection

abstract: In the \emph{incremental cycle detection} problem arcs are added to a directed acyclic graph and the algorithm has to report if the new arc closes a cycle. One seeks to minimize the total time to process the entire sequence of arc insertions, or until a cycle appears. In a recent breakthrough, Bender, Fineman, Gilbert and Tarjan \cite{BeFiGiTa11} presented two different algorithms, with time complexity $O(n^2 \log n)$ and $O(m \cdot \min \{m^{1/2}, n^{2/3} \})$, respectively. In this paper we introduce a new technique for incremental cycle detection that allows us to obtain both bounds (up to a logarithmic factor). Furthermore, our approach seems more amiable for distributed implementation.


title: Novel K Rattling: A New Route to Thermoelectric Materials?

abstract: We have performed ab initio molecular dynamics (MD) simulations to study the alkali-metal dynamics in the Al-doped (KAl0.33W1.67O6 and RbAl0.33W1.67O6) and undoped (KW2O6 and RbW2O6) defect pyrochlore tungstates. The K atoms exhibit novel rattling dynamics in both the doped and undoped tungstates while the Rb atoms do not. The KAl0.33W1.67O6 experimental thermal conductivity curve shows an unusual depression between ~ 50 K and ~ 250 K, coinciding with two crossovers in the K dynamics: the first at ~ 50 K, from oscillatory to diffusive, and the second at ~ 250 K, from diffusive back to oscillatory. We found that the low-temperature crossover is a result of the system transitioning below the activation energy of the diffusive dynamics whereas the high-temperature crossover is driven by a complex reconstruction of the local potential around the K atoms due to the cage dynamics. This leads to a hardening of the K potential with increasing temperature. This unusual reconstruction of the potential may have important implications for the interpretation of finite-temperature dynamics based on zero-temperature potentials in similar materials. The key result is that the novel K rattling, involving local diffusion, leads to a significant reduction in the thermal conductivity. We suggest that this may open a new route in the phonon engineering of cage compounds for thermoelectric materials where the rattlers are specifically selected to reduce the lattice thermal conductivity by the mechanism of local diffusion.


title: Bioinformatics Knowledge Transmission (training, learning, and teaching): overview and flexible comparison of computer based training approaches

abstract: The merger of computer science, mathematics, and life sciences has brought about the discipline known as bioinformatics. However, the transmission (e.g. training, learning, and teaching) of this knowledge becomes an important issue. Many tools have been developed to help the bioinformatics community with that transmission challenge. When selecting the best of these tools, called here BKTMS (Bioinformatics Knowledge Transmission Management Systems), there may be confusion. What makes a good BKTMS? How can we make this choice efficiently? These questions remain unanswered for many users (e.g. learner, teacher and student, trainer and trainee, administrator). This paper provides a critical review of 32 existing BKTMS and a flexible comparison. This review and evaluation will be used to gain insight into the tools, systems, and capabilities that will be added to or excluded from a new proposed model for the next generation of BKTMS, involving multidisciplinary, web semantic tools (e.g. web services, workflow) and standards like LOM, or SCORM.


title: Laser damage helps the eavesdropper in quantum cryptography

abstract: We propose a class of attacks on quantum key distribution (QKD) systems where an eavesdropper actively engineers new loopholes by using damaging laser illumination to permanently change properties of system components. This can turn a perfect QKD system into a completely insecure system. A proof-of-principle experiment performed on an avalanche photodiode-based detector shows that laser damage can be used to create loopholes. After about 1 W illumination, the detectors' dark count rate reduces 2 to 5 times, permanently improving single-photon counting performance. After about 1.5 W, the detectors switch permanently into the linear photodetection mode and become completely insecure for QKD applications.


title: Local Fourier Analysis of Multigrid Methods with Polynomial Smoothers and Aggressive coarsening

abstract: We focus on the study of multigrid methods with aggressive coarsening and polynomial smoothers for the solution of the linear systems corresponding to finite difference/element discretizations of the Laplace equation. Using local Fourier analysis we determine automatically the optimal values for the parameters involved in defining the polynomial smoothers and achieve fast convergence of cycles with aggressive coarsening. We also present numerical tests supporting the theoretical results and the heuristic ideas. The methods we introduce are highly parallelizable and efficient multigrid algorithms on structured and semi-structured grids in two and three spatial dimensions.


title: Intrinsic VHE Gamma-ray spectra of Blazars as a probe for Extragalactic Background Light

abstract: Very high energy (VHE) $\gamma$-rays above 10$'$s of GeV energy, emitted from distant blazars, are attenuated by photons from the extragalactic background light (EBL). Unfortunately, neither the EBL nor the intrinsic blazar spectrum is accurately known to derive one quantity from the other. In this work we use a homogeneous one zone model involving synchrotron, synchrotron self Compton (SSC) and external Compton (EC) emission mechanisms to estimate the intrinsic VHE spectra of blazars. The model is applied on three VHE blazars, namely PKS2155-304, RGB J0710+591 and 3C279, for which simultaneous multi-wavelength data are available from various observations. The predicted values of the intrinsic VHE fluxes are then compared with the observations by imaging atmospheric Cherenkov telescopes to determine the optical depth of VHE $\gamma$-rays. On comparing these optical depth values with those predicted by four different EBL models, we observe a somewhat pronounced systematic deviation for PKS2155-304 and 3C279 at higher energies, especially for the EBL model proposed by Finke et al.(2010). We attribute this deviation to be an outcome of either the failure of the extrapolation of blazar SED to VHE energies and/or due to various assumptions buried in the EBL models.


title: Density-based and transport-based core-periphery structures in networks

abstract: Networks often possess mesoscale structures, and studying them can yield insights into both structure and function. It is most common to study community structure, but numerous other types of mesoscale structures also exist. In this paper, we examine core-periphery structures based on both density and transport. In such structures, core network components are well-connected both among themselves and to peripheral components, which are not well-connected to anything. We examine core-periphery structures in a wide range of examples of transportation, social, and financial networks---including road networks in large urban areas, a rabbit warren, a dolphin social network, a European interbank network, and a migration network between counties in the United States. We illustrate that a recently developed transport-based notion of node coreness is very useful for characterizing transportation networks. We also generalize this notion to examine core versus peripheral edges, and we show that the resulting diagnostic is also useful for transportation networks. To examine the properties of transportation networks further, we develop a family of generative models of roadlike networks. We illustrate the effect of the dimensionality of the embedding space on transportation networks, and we demonstrate that the correlations between different measures of coreness can be very different for different types of networks.

categories: physics.soc-ph cond-mat.dis-nn cs.SI q-bio.OT


title: Provable Security of Networks

abstract: We propose a definition of {\it security} and a definition of {\it robustness} of networks against the cascading failure models of deliberate attacks and random errors respectively, and investigate the principles of the security and robustness of networks. We propose a {\it security model} such that networks constructed by the model are provably secure against any attacks of small sizes under the cascading failure models, and simultaneously follow a power law, and have the small world property with a navigating algorithm of time complex $O(\log n)$. It is shown that for any network $G$ constructed from the security model, $G$ satisfies some remarkable topological properties, including: (i) the {\it small community phenomenon}, that is, $G$ is rich in communities of the form $X$ of size poly logarithmic in $\log n$ with conductance bounded by $O(\frac{1}{|X|^{\beta}})$ for some constant $\beta$, (ii) small diameter property, with diameter $O(\log n)$ allowing a navigation by a $O(\log n)$ time algorithm to find a path for arbitrarily given two nodes, and (iii) power law distribution, and satisfies some probabilistic and combinatorial principles, including the {\it degree priority theorem}, and {\it infection-inclusion theorem}. By using these principles, we show that a network $G$ constructed from the security model is secure for any attacks of small scales under both the uniform threshold and random threshold cascading failure models. Our security theorems show that networks constructed from the security model are provably secure against any attacks of small sizes, for which natural selections of {\it homophyly, randomness} and {\it preferential attachment} are the underlying mechanisms.


title: Gas clumping in galaxy clusters

abstract: The reconstruction of galaxy cluster's gas density profiles is usually performed by assuming spherical symmetry and averaging the observed X-ray emission in circular annuli. In the case of a very inhomogeneous and asymmetric gas distribution, this method has been shown to return biased results in numerical simulations because of the $n^2$ dependence of the X-ray emissivity. We propose a method to recover the true density profiles in the presence of inhomogeneities, based on the derivation of the azimuthal median of the surface brightness in concentric annuli. We demonstrate the performance of this method with numerical simulations, and apply it to a sample of 31 galaxy clusters in the redshift range 0.04-0.2 observed with ROSAT/PSPC. The clumping factors recovered by comparing the mean and the median are mild and show a slight trend of increasing bias with radius. For $R<R_{500}$, we measure a clumping factor $\sqrt{C}<1.1$, which indicates that the thermodynamic properties and hydrostatic masses measured in this radial range are only mildly affected by this effect. Comparing our results with three sets of hydrodynamical numerical simulations, we found that non-radiative simulations significantly overestimate the level of inhomogeneities in the ICM, while the runs including cooling, star formation, and AGN feedback reproduce the observed trends closely. Our results indicate that most of the accretion of X-ray emitting gas is taking place in the diffuse, large-scale accretion patterns rather than in compact structures.


title: Harnack's inequality and Green functions on locally finite graphs

abstract: In this paper we study the gradient estimate for positive solutions of Schrodinger equations on locally finite graph. Then we derive Harnack's inequality for positive solutions of the Schrodinger equations. We also set up some results about Green functions of the Laplacian equation on locally finite graph. Interesting properties of Schrodinger equation are derived.

categories: math.DG cs.IT math.AP math.CA math.IT math.MG


title: Harnack inequalities for $1$-$d$ stochastic Klein-Gordon type equations

abstract: By the coupling method, we establish the Harnack inequalities, derivative formula and Driver's integration by parts formula for the stochastic Klein-Gordon type equations in the interval. We provide a detailed discussion about the nonlinear term. Some applications are given.


title: Cloud computing security using encryption technique

abstract: Cloud Computing has been envisioned as the next generation architecture of IT Enterprise. The Cloud computing concept offers dynamically scalable resources provisioned as a service over the Internet. Economic benefits are the main driver for the Cloud, since it promises the reduction of capital expenditure and operational expenditure. In order for this to become reality, however, there are still some challenges to be solved. Most important among these are security and trust issues,since the users data has to be released to the Cloud and thus leaves the protection sphere of the data owner.In contrast to traditional solutions, where the IT services are under proper physical,logical and personnel controls, Cloud Computing moves the application software and databases to the large data centers, where the management of the data and services may not be fully trustworthy. This unique attribute, however, poses many new security challenges which have not been well understood. Security is to save data from danger and vulnerability. There are so many dangers and vulnerabilities to be handled. Various security issues and some of their solution are explained and are concentrating mainly on public cloud security issues and their solutions. Data should always be encrypted when stored(using separate symmetric encryption keys)and transmitted. If this is implemented appropriately, even if another tenant can access the data, all that will appear is gibberish. So a method is proposed such that we are encrypting the whole data along with the cryptographic key.


title: Surface bulk differences in a conventional superconductor, ZrB(12)

abstract: We studied the electronic structure of a conventional superconductor, ZrB(12) using high resolution x-photoemission spectroscopy and single crystalline samples. Experimental results with different bulk sensitivity reveals boron deficiency and different valence states of Zr at the surface relative to the bulk. Signature of a satellite features is observed in the Zr core level spectra corresponding to the bulk of the material suggesting importance of electron correlation among the conduction electrons in the bulk while the surface appears to be uncorrelated. These results provide an insight in fabricating devices based on such superconductors.


title: Impurity Screening and Surface Acoustic Wave Absorption in a Dipolar Exciton Condensate at Finite Temperatures

abstract: We describe the behavior of a repulsively interacting Bose-Einstein condensate of indirect dipolar exciton gas in a double quantum well (QW) system under external static or dynamic electric fields at finite temperatures. Electrostatic perturbation is considered to be created by an impurity atom or shot-range defect of QW fluctuation. The screening of this defect potential by an exciton condensate is studied. We find asymptotic spatial dependence of the screened potential and analyse its dependence on the temperature and exciton concentration. It is shown that the asymptotic of the screened potential has a steep power law dependence in contrast to the well known results of electron gas. This peculiarity reflects the bosonic nature of the exciton condensate. The behavior of exciton condensate under external alternative field created by a surface acoustic wave (SAW) is examined in detail. We focus our attention on the dependence of SAW absorption coefficient on temperature and exciton concentration. We found that at zero temperatures Landau damping does not contribute to the SAW absorption, but the Belyaev mechanism produces unusual behavior of SAW absorption coefficient on exciton concentration: if the exciton concentration exceeds some critical value, the SAW absorption vanishes. At finite temperatures Landau damping comes into action and results in washing out the sharp absorption behavior. Such unusual SAW absorption properties can be used for experimental evidence of the exciton condensation. This method is also applicable to the experimental testing of both dark and bright exciton condensates, that is impossible to do with the optical luminescence technique.


title: Sedimentation of granular columns in the viscous and weakly inertial regimes

abstract: We investigate the dynamics of granular columns of point particles that interact via long-range hydrodynamic interactions and fall under the action of gravity. We investigate the influence of inertia using the Green's function for the Oseen equation. The initial conditions (density and aspect ratio) are systematically varied. Our results suggest that universal self-similar laws may be sufficient to characterize the temporal and structural evolution of the granular columns. A characteristic time above which an instability is triggered (which may enable the formation of clusters) is also retrieved and discussed.


title: Tunable and Growing Network Generation Model with Community Structures

abstract: Recent years have seen a growing interest in the modeling and simulation of social networks to understand several social phenomena. Two important classes of networks, small world and scale free networks have gained a lot of research interest. Another important characteristic of social networks is the presence of community structures. Many social processes such as information diffusion and disease epidemics depend on the presence of community structures making it an important property for network generation models to be incorporated. In this paper, we present a tunable and growing network generation model with small world and scale free properties as well as the presence of community structures. The major contribution of this model is that the communities thus created satisfy three important structural properties: connectivity within each community follows power-law, communities have high clustering coefficient and hierarchical community structures are present in the networks generated using the proposed model. Furthermore, the model is highly robust and capable of producing networks with a number of different topological characteristics varying clustering coefficient and inter-cluster edges. Our simulation results show that the model produces small world and scale free networks along with the presence of communities depicting real world societies and social networks.


title: Linear Convergence on Positively Homogeneous Functions of a Comparison Based Step-Size Adaptive Randomized Search: the (1+1) ES with Generalized One-fifth Success Rule

abstract: In the context of unconstraint numerical optimization, this paper investigates the global linear convergence of a simple probabilistic derivative-free optimization algorithm (DFO). The algorithm samples a candidate solution from a standard multivariate normal distribution scaled by a step-size and centered in the current solution. This solution is accepted if it has a better objective function value than the current one. Crucial to the algorithm is the adaptation of the step-size that is done in order to maintain a certain probability of success. The algorithm, already proposed in the 60's, is a generalization of the well-known Rechenberg's $(1+1)$ Evolution Strategy (ES) with one-fifth success rule which was also proposed by Devroye under the name compound random search or by Schumer and Steiglitz under the name step-size adaptive random search. In addition to be derivative-free, the algorithm is function-value-free: it exploits the objective function only through comparisons. It belongs to the class of comparison-based step-size adaptive randomized search (CB-SARS). For the convergence analysis, we follow the methodology developed in a companion paper for investigating linear convergence of CB-SARS: by exploiting invariance properties of the algorithm, we turn the study of global linear convergence on scaling-invariant functions into the study of the stability of an underlying normalized Markov chain (MC). We hence prove global linear convergence by studying the stability (irreducibility, recurrence, positivity, geometric ergodicity) of the normalized MC associated to the $(1+1)$-ES. More precisely, we prove that starting from any initial solution and any step-size, linear convergence with probability one and in expectation occurs. Our proof holds on unimodal functions that are the composite of strictly increasing functions by positively homogeneous functions with degree $\alpha$ (assumed also to be continuously differentiable). This function class includes composite of norm functions but also non-quasi convex functions. Because of the composition by a strictly increasing function, it includes non continuous functions. We find that a sufficient condition for global linear convergence is the step-size increase on linear functions, a condition typically satisfied for standard parameter choices. While introduced more than 40 years ago, we provide here the first proof of global linear convergence for the $(1+1)$-ES with generalized one-fifth success rule and the first proof of linear convergence for a CB-SARS on such a class of functions that includes non-quasi convex and non-continuous functions. Our proof also holds on functions where linear convergence of some CB-SARS was previously proven, namely convex-quadratic functions (including the well-know sphere function).


title: Weak Minkowski Spaces

abstract: We define the notion of weak Minkowski metric and prove some basic properties of such metrics. We also highlight some of the important analogies between Minkowski geometry and the Funk and Hilbert geometries.


title: Absolute Algebra III - The saturated spectrum

abstract: Let B1 denote the set {0,1} with the usual operations except that $1+1=1$, in other words, the smallest characteristic 1 semifield . We compare two possible analogues of the notion of prime ideal for B1--algebras. We then consider the relations between these notions and Deitmar's theory of F1--schemes.


title: Prime and primary ideals in characteristic one

abstract: We study zero divisors and minimal prime ideals in semirings of characteristic one. Thereafter we find a counterexample to the most obvious version of primary decomposition, but are able to establish a weaker version. Lastly, we study Evans'condition in this context.


title: On the commuting probability and supersolvability of finite groups

abstract: For a finite group $G$, let $d(G)$ denote the probability that a randomly chosen pair of elements of $G$ commute. We prove that if $d(G)>1/s$ for some integer $s>1$ and $G$ splits over an abelian normal nontrivial subgroup $N$, then $G$ has a nontrivial conjugacy class inside $N$ of size at most $s-1$. We also extend two results of Barry, MacHale, and N\'{\i} Sh\'{e} on the commuting probability in connection with supersolvability of finite groups. In particular, we prove that if $d(G)>5/16$ then either $G$ is supersolvable, or $G$ isoclinic to $A_4$, or $G/\Center(G)$ is isoclinic to $A_4$.


title: Spectral evolution in a Shastry-Sutherland lattice, HoB(4)

abstract: We studied the electronic structure of a Shastry-Sutherland lattice system, HoB4 employing high resolution photoemission spectroscopy and ab initio band structure calculations. The surface and bulk borons exhibit subtle differences, and loss of boron compared to the stoichiometric bulk. However, the surface and bulk conduction bands near Fermi level are found to be similar. Evolution of the electronic structure with temperature is found to be similar to that observed in a typical charge-disordered system. A sharp dip is observed at the Fermi level in the low temperature spectra revealing signature of antiferromagnetic gap. Asymmetric spectral weight transfer with temperature manifests particle-hole asymmetry that may be related to the exotic properties of these systems.


title: An Existential Proof of the Conjecture on Packing Anchored Rectangles

abstract: Let $P_{n}$ be a set of $n$ points, including the origin, in the unit square $U = [0,1]^2$. We consider the problem of constructing $n$ axis-parallel and mutually disjoint rectangles inside $U$ such that the bottom-left corner of each rectangle coincides with a point in $P_{n}$ and the total area covered by the rectangles is maximized \cite{ibmpuzzle}, \cite{Winkler2007}, \cite{Winkler2010a}, \cite{Winkler2010b}. The longstanding conjecture has been that at least half of $U$ can be covered when such rectangles are properly placed. In this paper, we give an existential proof of the conjecture.


title: Higgs decays to $\gamma$ and invisible particles in the standard model

abstract: Using the Higgs boson mass $m_h=125$ GeV, the radiative Higgs decays $h\rightarrow\gamma \nu_l\bar\nu_l$ with $\nu_l = \nu_e,\,\nu_\mu$ and $\nu_\tau$ are analyzed in the standard model. Our calculation shows that the inclusive width of these processes, i.e., the sum of $\Gamma(h\to\gamma\nu_l\bar\nu_l)$ for ${\nu_l=\nu_e,\nu_\mu,\nu_\tau}$, is $1.41$ keV, which is about $15\%$ of $\Gamma(h\to\gamma\gamma)$. Therefore, the observation of these channels in the future precise experiments may provide us some useful information on the Higgs physics both in the standard model and in its possible extensions.


title: Hard Lefschetz Property for Hamiltonian torus actions on 6-dimensional GKM manifolds

abstract: In this paper, we study the hard Lefschetz property of a symplectic manifold which admits a Hamiltonian torus action. More precisely, let $(M,\omega)$ be a 6-dimensional compact symplectic manifold with a Hamiltonian $T^2$-action. We will show that if the moment map image of $M$ is a GKM-graph and if the graph is index-increasing, then $(M,\omega)$ satisfies the hard Lefschetz property.


title: Flows on Bidirected Graphs

abstract: The study of nowhere-zero flows began with a key observation of Tutte that in planar graphs, nowhere-zero k-flows are dual to k-colourings (in the form of k-tensions). Tutte conjectured that every graph without a cut-edge has a nowhere-zero 5-flow. Seymour proved that every such graph has a nowhere-zero 6-flow. For a graph embedded in an orientable surface of higher genus, flows are not dual to colourings, but to local-tensions. By Seymour's theorem, every graph on an orientable surface without the obvious obstruction has a nowhere-zero 6-local-tension. Bouchet conjectured that the same should hold true on non-orientable surfaces. Equivalently, Bouchet conjectured that every bidirected graph with a nowhere-zero $\mathbb{Z}$-flow has a nowhere-zero 6-flow. Our main result establishes that every such graph has a nowhere-zero 12-flow.


title: Explicit round fold maps on some fundamental manifolds

abstract: Stable fold maps are fundamental tools in a generalization of the theory of Morse functions on smooth manifolds and its application to studies of topological properties of smooth manifolds. Round fold maps were introduced as stable fold maps with singular value sets, defined as the set consisting of all the singular values, of concentric spheres by the author in 2013; for example, some special generic maps on spheres are regarded as round fold maps whose singular value sets are connected. Algebraic invariants such as homology and homotopy groups of manifolds admitting round fold maps and more precisely, the homeomorphism and diffeomorphism types of manifolds admitting such maps having appropriate differential topological structures were studied. Moreover, explicit round fold maps into the Eucidean space of dimension larger than $1$ are constructed on some fundamental manifolds such as manifolds having the structures of bundles over the standard sphere of dimension equal to the Euclidean space whose fibers are closed smooth manifolds and manifolds of dimension not smaller than twice the dimension of the Euclidean space represented as the connected sum of manifolds having the structures of bundles over the standard sphere of dimension equal to the Euclidean space whose fibers are diffeomorphic to standard spheres. In this paper, we construct new explicit fold maps on some fundamental manifolds including the manifolds before by using extended methods of ones used in the constructions before.


title: All Linear-Time Congruences for Familiar Operators

abstract: The detailed behaviour of a system is often represented as a labelled transition system (LTS) and the abstract behaviour as a stuttering-insensitive semantic congruence. Numerous congruences have been presented in the literature. On the other hand, there have not been many results proving the absence of more congruences. This publication fully analyses the linear-time (in a well-defined sense) region with respect to action prefix, hiding, relational renaming, and parallel composition. It contains 40 congruences. They are built from the alphabet, two kinds of traces, two kinds of divergence traces, five kinds of failures, and four kinds of infinite traces. In the case of finite LTSs, infinite traces lose their role and the number of congruences drops to 20. The publication concentrates on the hardest and most novel part of the result, that is, proving the absence of more congruences.


title: Enhanced resonant force between two entangled identical atoms in a photonic crystal

abstract: We consider the resonant interaction energy and force between two identical atoms, one in an excited state and the other in the ground state, placed inside a photonic crystal. The atoms, having the same orientation of their dipole moment, are supposed prepared in their symmetrical state and interact with the quantum electromagnetic field. We consider two specific models of photonic crystals: a one-dimensional model and an isotropic model. We show that in both cases the resonant interatomic force can be strongly enhanced by the presence of the photonic crystal, as a consequence of the modified dispersion relation and density of states, in particular if the transition frequency of the atoms is close to the edge of a photonic gap. Differences between the two models considered of photonic crystal are discussed in detail, as well as comparison with the analogous system of two impurity atoms in a quantum semiconductor wire. A numerical estimate of the effect in a realistic situation is also discussed.


title: Photonic processes in Born-Infeld theory

abstract: We study the processes of photon-photon scattering and photon splitting in a magnetic field in Born-Infeld theory. In both cases we combine the terms from the tree-level Born-Infeld Lagrangian with the usual one-loop QED contributions, where those are approximated by the Euler-Heisenberg Lagrangian, including also the interference terms. For photon-photon scattering we obtain the total cross section in the low-energy approximation. For photon splitting we compute the total absorption coefficient in the hexagon (weak field) approximation, and also show that, due to the non-birefringence property of Born-Infeld theory, the selection rules found by Adler for the QED case continue to hold in this more general setting. We discuss the bounds on the free parameter of Born-Infeld theory that may be obtained from this type of processes.


title: Verification by stochastic Perron's method in stochastic exit time control problems

abstract: We apply the Stochastic Perron method, created by Bayraktar and S\^irbu, to a stochastic exit time control problem. Our main assumption is the validity of the Strong Comparison Result for the related Hamilton-Jacobi-Bellman (HJB) equation. Without relying on Bellman's optimality principle we prove that inside the domain the value function is continuous and coincides with a viscosity solution of the Dirichlet boundary value problem for the HJB equation.


title: Wetting Boundary Conditions in Phase-Field-Based Simulation of Binary Fluids: Some Comparative Studies and New Development

abstract: We studied several wetting boundary conditions (WBCs) in the simulation of binary fluids based on phase-field theory. Five WBCs, three belonging to the surface energy (SE) formulation using the linear, cubic and sine functions (denoted as LinSE, CubSE and SinSE), the fourth using a geometric formulation (Geom), and the fifth using a characteristic interpolation (CI), were compared with each other through the study of several problems: (1) the static contact angle of a drop; (2) a Poiseuille flow-driven liquid column; (3) a wettability gradient (WG)-driven liquid column; (4) drop dewetting. It was found that while all WBCs can predict the static contact angle fairly accurately, they may affect the simulation outcomes of dynamic problems differently, depending on the driving mechanism. For the flow-driven problem, to use different WBCs had almost no effect on the flow characteristics over a large scale. But for other capillarity-driven problems, the WBC had some noticeable effects. For the WG-driven liquid column, Geom gave the most consistent prediction between the drop velocity and dynamic contact angles, and LinSE delivered the poorest prediction in this aspect. Except for Geom, the dynamic contact angle differed from the prescribed (static) one when other WBCs were used. For drop dewetting, Geom led to the most violent drop motion whereas CubSE caused the weakest motion; the initial contact line velocity was also found to be dependent on the WBC. For several problems, CubSE and SinSE gave almost the same results, and those by Geom and CI were close as well, possibly due to similar consideration in their design. Besides various comparisons, a new implementation that may be used for all WBCs was proposed to mimic the wall energy relaxation and control the degree of slip. This new procedure made it possible to allow the simulations to match experimental measurements well.


title: Nilpotent and Abelian Hall Subgroups in Finite Groups

abstract: We give a characterization of the finite groups having nilpotent or abelian Hall $\pi$-subgroups which can easily be verified from the character table.


title: Combinatorics of Link Diagrams and Volume

abstract: We show that the volumes of certain hyperbolic A-adequate links can be bounded (above and) below in terms of two diagrammatic quantities: the twist number and the number of certain alternating tangles in an A-adequate diagram. We then restrict our attention to plat closures of certain braids, a rich family of links whose volumes can be bounded in terms of the twist number alone. Furthermore, in the absence of special tangles, our volume bounds can be expressed in terms of a single stable coefficient of the colored Jones polynomial. Consequently, we are able to provide a new collection of links that satisfy a Coarse Volume Conjecture.


title: "Cracking-the-whip" effect stretches driven polymers

abstract: Cracking the whip accelerates the tail of a chain to hit the air loudly and clearly. We proved that the similar acceleration effect causes coil deformation of driven chain-like polymers. We first preformed Monte Carlo simulations of a single driven polymer coil to demonstrate its deformation in company with faster or slower deviations of velocities. We then performed parallel Brownian Dynamics simulations to demonstrate that the coil deformation can be caused by the so-called "cracking-the-whip" effect due to non-synchronous biased Brownian motions of monomers inherited in Monte Carlo simulations. Since such non-synchronous motions represent random perturbations in the environmentally dependent potential energy landscape or mobility, reflecting heterogeneous dynamics of polymers in the liquid phase, our observations bring new insights into the non-linear dynamics of driven chain-like polymers.


title: Holographic Approach to Nonequilibrium Dynamics of Moving Mirrors Coupled to Quantum Critical Theories

abstract: We employ the holographic method to study fluctuations and dissipation of an $n$-dimensional moving mirror coupled to quantum critical theories in $d$ spacetime dimensions. The bulk counterpart of the mirror with perfect reflection is a D$(n+1)$ brane in the Lifshitz geometry of $d+1$ dimensions. The motion of the mirror can be realized from the dynamics of the brane at the boundary of the bulk. The excited modes of the brane in the bulk render the mirror undergoing Brownian motion. For small displacement of the mirror, we derive the analytical results of the correlation functions and response functions. The dynamics of the mirror due to small fluctuations around the brane vacuum state in the bulk is found supraohmic so that after initial growth, the velocity fluctuations approach a saturated value at late time with a power-law behavior. On the contrary, in the Lifshitz black hole background, the mirror in thermal fluctuations shows that its relaxation dynamics becomes ohmic, and the saturation of velocity fluctuations is reached exponentially in time. Finally a comparison is made with the result of a moving mirror driven by free fields.


title: Cycling Dynamo in a Young Sun: Grand Minima and Equatorward Propagation

abstract: We assess the global-scale dynamo action achieved in a simulation of a sun-like star rotating at three times the solar rate. The 3-D MHD Anelastic Spherical Harmonic code, using slope-limited diffusion, is employed to capture convection and dynamo processes in such a young sun. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. We find that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals every 6.2 years. The magnetic energy also rises and falls with a regular period, with two magnetic energy cycles required to complete a polarity cycle. These magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the polarity reversals are present due to the spatial separation of the equatorial and polar dynamos. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation and to a near-surface shear layer. This simulation also enters a grand minimum lasting roughly 20 years, after which the dynamo recovers its regular polarity cycles.


title: An efficient distributed learning algorithm based on effective local functional approximations

abstract: Scalable machine learning over big data is an important problem that is receiving a lot of attention in recent years. On popular distributed environments such as Hadoop running on a cluster of commodity machines, communication costs are substantial and algorithms need to be designed suitably considering those costs. In this paper we give a novel approach to the distributed training of linear classifiers (involving smooth losses and L2 regularization) that is designed to reduce the total communication costs. At each iteration, the nodes minimize locally formed approximate objective functions; then the resulting minimizers are combined to form a descent direction to move. Our approach gives a lot of freedom in the formation of the approximate objective function as well as in the choice of methods to solve them. The method is shown to have $O(log(1/\epsilon))$ time convergence. The method can be viewed as an iterative parameter mixing method. A special instantiation yields a parallel stochastic gradient descent method with strong convergence. When communication times between nodes are large, our method is much faster than the Terascale method (Agarwal et al., 2011), which is a state of the art distributed solver based on the statistical query model (Chuet al., 2006) that computes function and gradient values in a distributed fashion. We also evaluate against other recent distributed methods and demonstrate superior performance of our method.


title: Valence double parton distributions of the nucleon in a simple model

abstract: Valence double parton distribution functions of the nucleon are evaluated in the framework of a simple model, where the conservation of the longitudinal momentum is taken into account. The leading-order DGLAP QCD evolution from the low quark-model scale to higher renormalization scales is carried out via the Mellin moments of the distributions. Results of the valence quark correlation function show that in general the double distributions cannot be approximated as a product of the single-particle distributions.


title: SMaSH: A Benchmarking Toolkit for Human Genome Variant Calling

abstract: Motivation: Computational methods are essential to extract actionable information from raw sequencing data, and to thus fulfill the promise of next-generation sequencing technology. Unfortunately, computational tools developed to call variants from human sequencing data disagree on many of their predictions, and current methods to evaluate accuracy and computational performance are ad-hoc and incomplete. Agreement on benchmarking variant calling methods would stimulate development of genomic processing tools and facilitate communication among researchers. Results: We propose SMaSH, a benchmarking methodology for evaluating human genome variant calling algorithms. We generate synthetic datasets, organize and interpret a wide range of existing benchmarking data for real genomes, and propose a set of accuracy and computational performance metrics for evaluating variant calling methods on this benchmarking data. Moreover, we illustrate the utility of SMaSH to evaluate the performance of some leading single nucleotide polymorphism (SNP), indel, and structural variant calling algorithms. Availability: We provide free and open access online to the SMaSH toolkit, along with detailed documentation, at smash.cs.berkeley.edu.


title: Multiple positive solutions for a nonlinear three-point integral boundary-value problem

abstract: We investigate the existence of positive solutions to the nonlinear second-order three-point integral boundary value problem \begin{equation*} \label{eq-1} \begin{gathered} {u^{\prime \prime}}(t)+f(t, u(t))=0,\ 0<t<T, \\ u(0)={\beta}u(\eta),\ u(T)={\alpha}\int_{0}^{\eta}u(s)ds, \end{gathered} \end{equation*} where $0<{\eta}<T$, $0<{\alpha}< \frac{2T}{{\eta}^{2}}$, $0<{\beta}<\frac{2T-\alpha\eta^{2}}{\alpha\eta^{2}-2\eta+2T}$ are given constants. We establish the existence of at least three positive solutions by using the Leggett-Williams fixed-point theorem.


title: Recurrence statistics for the space of Interval Exchange maps and the Teichm\"uller flow on the space of translation surfaces

abstract: In this note we show that the transfer operator of a Rauzy-Veech-Zorich renormalization map acting on a space of quasi-H\"older functions is quasicompact and derive certain statistical recurrence properties for this map and its associated Teichm\"uller flow. We establish Borel-Cantelli lemmas, Extreme Value statistics and return time statistics for the map and flow. Previous results have established quasicompactness in H\"older or analytic function spaces, for example the work of M. Pollicott and T. Morita. The quasi-H\"older function space is particularly useful for investigating return time statistics. In particular we establish the shrinking target property for nested balls in the setting of Teichm\"uller flow. Our point of view, approach and terminology derives from the work of M. Pollicott augmented by that of M. Viana.


title: Flux tubes and coherence length in the SU(3) vacuum

abstract: An estimate of the London penetration and coherence lengths in the vacuum of the SU(3) pure gauge theory is given downstream an analysis of the transverse profile of the chromoelectric flux tubes. Within ordinary superconductivity, a simple variational model for the magnitude of the normalized order parameter of an isolated vortex produces an analytic expression for magnetic field and supercurrent density. In the picture of SU(3) vacuum as dual superconductor, this expression provides us with the function that fits the chromoelectric field data. The smearing procedure is used in order to reduce noise.


title: On the concavity of the arithmetic volumes

abstract: In this note, we study the differentiability of the arithmetic volumes along arithmetic R-divisors, and give some equality conditions for the Brunn-Minkowski inequality for arithmetic volumes over the cone of nef and big arithmetic R-divisors.


title: Polynomial spaces reproduced by elliptic scaling functions

abstract: The Strang-Fix conditions are necessary and sufficient to reproduce spaces of algebraic polynomials up to some degree by integer shifts of compactly supported functions. W. Dahmen and Ch. Micchelli (Linear Algebra Appl. 52/3:217-234, 1983) introduced a generalization of the Strang-Fix conditions to affinely invariant subspaces of higher degree polynomials. C. de Boor (Constr. Approx. 3:199-208, 1987) raised a question on the necessity of scale-invariance of polynomial space for an arbitrary function; and he omitted the scale-invariance restriction on the space. In the paper, we present a matrix approach to determine (not necessarily scale-invariant) polynomial space contained in the span of integer shifts of a compactly supported function. Also, in the paper, we consider scaling functions that we call the elliptic scaling functions (Int. J. Wavelets Multiresolut. Inf. Process. To appear); and, using the matrix approach, we prove that the elliptic scaling functions satisfy the generalized Strang-Fix conditions and reproduce only affinely invariant polynomial spaces. Namely we prove that any algebraic polynomial contained in the span of integer shifts of a compactly supported elliptic scaling function belongs to the null-space of a homogeneous elliptic differential operator. However, in the paper, we present nonstationary elliptic scaling functions such that the scaling functions reproduce not scale-invariant (only shift-invariant) polynomial spaces.


title: Dynamic connectivity algorithms for Monte Carlo simulations of the random-cluster model

abstract: We review Sweeny's algorithm for Monte Carlo simulations of the random cluster model. Straightforward implementations suffer from the problem of computational critical slowing down, where the computational effort per edge operation scales with a power of the system size. By using a tailored dynamic connectivity algorithm we are able to perform all operations with a poly-logarithmic computational effort. This approach is shown to be efficient in keeping online connectivity information and is of use for a number of applications also beyond cluster-update simulations, for instance in monitoring droplet shape transitions. As the handling of the relevant data structures is non-trivial, we provide a Python module with a full implementation for future reference.


title: Origin of two-dimensional electron gases at oxide interfaces: insights from theory

abstract: The response of oxide thin films to polar discontinuities at interfaces and surfaces has generated an enormous activity due to the variety of interesting effects it gives rise to. A case in point is the discovery of the electron gas at the interface between LaAlO3 and SrTiO3, which has since been shown to be quasi-two-dimensional, switchable, magnetic and/or superconducting. Despite these findings, the origin of the two-dimensional electron gas is highly debated and several possible mechanisms remain. Here we review the main proposed mechanisms and attempt to model expected effects in a quantitative way with the ambition of better constraining what effects can/cannot explain the observed phenomenology. We do it in the framework of a phenomenological model for understanding electronic and/or redox screening of the chemical charge in oxide heterostructures. We also discuss the effect of intermixing, both conserving and non-conserving the total stoichiometry.