title: Interior Regularity Estimates in High Conductivity Homogenization and Application

abstract: In this paper, uniform pointwise regularity estimates for the solutions of conductivity equations are obtained in a unit conductivity medium reinforced by a epsilon-periodic lattice of highly conducting thin rods. The estimates are derived only at a distance epsilon^{1+tau} (for some tau>0) away from the fibres. This distance constraint is rather sharp since the gradients of the solutions are shown to be unbounded locally in L^p as soon as p>2. One key ingredient is the derivation in dimension two of regularity estimates to the solutions of the equations deduced from a Fourier series expansion with respect to the fibres direction, and weighted by the high-contrast conductivity. The dependence on powers of epsilon of these two-dimensional estimates is shown to be sharp. The initial motivation for this work comes from imaging, and enhanced resolution phenomena observed experimentally in the presence of micro-structures. We use these regularity estimates to characterize the signature of low volume fraction heterogeneities in the fibred reinforced medium assuming that the heterogeneities stay at a distance epsilon^{1+tau} away from the fibres.


title: A new class of codes for Boolean masking of cryptographic computations

abstract: We introduce a new class of rate one-half binary codes: {\bf complementary information set codes.} A binary linear code of length $2n$ and dimension $n$ is called a complementary information set code (CIS code for short) if it has two disjoint information sets. This class of codes contains self-dual codes as a subclass. It is connected to graph correlation immune Boolean functions of use in the security of hardware implementations of cryptographic primitives. Such codes permit to improve the cost of masking cryptographic algorithms against side channel attacks. In this paper we investigate this new class of codes: we give optimal or best known CIS codes of length $<132.$ We derive general constructions based on cyclic codes and on double circulant codes. We derive a Varshamov-Gilbert bound for long CIS codes, and show that they can all be classified in small lengths $\le 12$ by the building up construction. Some nonlinear permutations are constructed by using $\Z_4$-codes, based on the notion of dual distance of an unrestricted code.


title: Efficient Encoding of Watermark Numbers as Reducible Permutation Graphs

abstract: In a software watermarking environment, several graph theoretic watermark methods use numbers as watermark values, where some of these methods encode the watermark numbers as graph structures. In this paper we extended the class of error correcting graphs by proposing an efficient and easily implemented codec system for encoding watermark numbers as reducible permutation flow-graphs. More precisely, we first present an efficient algorithm which encodes a watermark number $w$ as self-inverting permutation $\pi^*$ and, then, an algorithm which encodes the self-inverting permutation $\pi^*$ as a reducible permutation flow-graph $F[\pi^*]$ by exploiting domination relations on the elements of $\pi^*$ and using an efficient DAG representation of $\pi^*$. The whole encoding process takes O(n) time and space, where $n$ is the binary size of the number $w$ or, equivalently, the number of elements of the permutation $\pi^*$. We also propose efficient decoding algorithms which extract the number $w$ from the reducible permutation flow-graph $F[\pi^*]$ within the same time and space complexity. The two main components of our proposed codec system, i.e., the self-inverting permutation $\pi^*$ and the reducible permutation graph $F[\pi^*]$, incorporate important structural properties which make our system resilient to attacks.


title: Reconstruction of Jet Properties in Pb-Pb Collisions with the ALICE-Experiment

abstract: The study of the hot and dense medium created in heavy-ion collisions via the effect of parton energy loss is one of the major goals of jet and high-\pt\ measurements in these reactions. In particular the modifications of the longitudinal and transverse jet structure. Here the subtraction of the underlying heavy-ion event, as well as a precise knowledge of the background induced fluctuations of the reconstructed jet energy is required. We extract the modification of the near- and away-side jet yield in heavy-ion collisions at the LHC and present a detailed assessment of background fluctuations for jet reconstruction.


title: Measuring the Higgs boson couplings

abstract: Once a Higgs boson has been discovered, it is also important to know what are its properties, in particular its couplings to the other particles. At the LHC, there will be many observable channels that can be used to measure the relevant parameters. Using the SFitter framework, we project these measurements onto a weak-scale effective theory with general Higgs boson couplings. Thereby, our analysis profits from the error treatment and the parameter determination tools already used for new-physics searches. We use this to study both the individual couplings and their associated errors as well as correlations between them.


title: Structural variability of 3C 111 on parsec scales

abstract: We discuss the parsec-scale structural variability of the extragalactic jet 3C 111 related to a major radio flux density outburst in 2007. The data analyzed were taken within the scope of the MOJAVE, UMRAO, and F-GAMMA programs, which monitor a large sample of the radio brightest compact extragalactic jets with the VLBA, the University of Michigan 26 m, the Effelsberg 100 m, and the IRAM 30 m radio telescopes. The analysis of the VLBA data is performed by fitting Gaussian model components in the visibility domain. We associate the ejection of bright features in the radio jet with a major flux-density outburst in 2007. The evolution of these features suggests the formation of a leading component and multiple trailing components.


title: On Joint Diagonalisation for Dynamic Network Analysis

abstract: Joint diagonalisation (JD) is a technique used to estimate an average eigenspace of a set of matrices. Whilst it has been used successfully in many areas to track the evolution of systems via their eigenvectors; its application in network analysis is novel. The key focus in this paper is the use of JD on matrices of spanning trees of a network. This is especially useful in the case of real-world contact networks in which a single underlying static graph does not exist. The average eigenspace may be used to construct a graph which represents the `average spanning tree' of the network or a representation of the most common propagation paths. We then examine the distribution of deviations from the average and find that this distribution in real-world contact networks is multi-modal; thus indicating several \emph{modes} in the underlying network. These modes are identified and are found to correspond to particular times. Thus JD may be used to decompose the behaviour, in time, of contact networks and produce average static graphs for each time. This may be viewed as a mixture between a dynamic and static graph approach to contact network analysis.


title: Factorial cluster algebras

abstract: We show that cluster algebras do not contain non-trivial units and that all cluster variables are irreducible elements. Both statements follow from Fomin and Zelevinsky's Laurent phenomenon. As an application we give a criterion for a cluster algebra to be a factorial algebra. This can be used to construct cluster algebras, which are isomorphic to polynomial rings. We also study various kinds of upper bounds for cluster algebras, and we prove that factorial cluster algebras coincide with their upper bounds.


title: The nucleon axial mass and the MiniBooNE CCQE neutrino-nucleus data

abstract: We analyze the MiniBooNE CCQE $d\sigma/dT_\mu d\cos\theta_\mu$ data using a theoretical model that has proved to be quite successful in the analysis of nuclear reactions with electron, photon and pion probes. We find that RPA and multinucleon knockout turn out to be essential for the description of the MiniBooNE data. We show these measurements are fully compatible with former determinations of nucleon axial mass $M_A$, in contrast with several previous analyses, which have suggested an anomalously large value. We find, $M_A=1.08 \pm 0.03$ GeV. We also argue that the procedure, commonly used to reconstruct the neutrino energy for QE events from the muon angle and energy, could be unreliable for a wide region of the phase space, due to the large importance of multinucleon events.


title: Implementation of PhotoZ under Astro-WISE - A photometric redshift code for large datasets

abstract: We describe the implementation of the PhotoZ code in the framework of the Astro-WISE package and as part of the Photometric Classification Server of the PanSTARRS pipeline. Both systems allow the automatic measurement of photometric redshifts for the millions of objects being observed in the PanSTARRS project or expected to be observed by future surveys like KIDS, DES or EUCLID.


title: Adaptive schemes for incomplete quantum process tomography

abstract: We propose an iterative algorithm for incomplete quantum process tomography, with the help of quantum state estimation, based on the combined principles of maximum-likelihood and maximum-entropy. The algorithm yields a unique estimator for an unknown quantum process when one has less than a complete set of linearly independent measurement data to specify the quantum process uniquely. We apply this iterative algorithm adaptively in various situations and so optimize the amount of resources required to estimate the quantum process with incomplete data.


title: Stochastic potential and quantum decoherence of heavy quarkonium in the quark-gluon plasma

abstract: We propose an open quantum systems approach to the physics of heavy quarkonia in a thermal medium, based on stochastic quantum evolution. This description emphasizes the importance of collisions with the environment and focuses on the concept of spatial decoherence of the heavy quarkonium wave function. It is shown how to determine the parameters of the dynamical evolution, i.e. the real potential and the noise strength, from a comparison with quantities to be obtained from lattice QCD. Furthermore the imaginary part of the lattice QCD heavy quark potential is found to be naturally related to the strength of the noise correlations. We discuss the time evolution of $Q\bar{Q}$ analytically in a limiting scenario for the spatial decoherence and provide a qualitative 1-dimensional numerical simulation of the real-time dynamics.


title: A supercircle description of universal three-body states in two dimensions

abstract: We consider bound states of asymmetric three-body systems confined to two dimensions. In the universal regime, two energy ratios and two mass ratios provide complete knowledge of the three-body energy measured in units of one two-body energy. The lowest number of stable bound states is produced when one mass is larger than two similar masses. We focus on selected asymmetric systems of interest in cold atom physics. The scaled three-body energy and the two scaled two-body energies are related through an equation for a supercircle whose radius increases almost linearly with three-body energy. The exponents exhibit an increasing behavior with three-body energy. The mass dependence is highly non-trivial. We give a simple relation that predicts the universal three-body energy.

categories: cond-mat.quant-gas cond-mat.mes-hall nucl-th quant-ph


title: Ising metamagnets in thin film geometry: equilibrium properties

abstract: Artificial antiferromagnets and synthetic metamagnets have attracted much attention recently due to their potential for many different applications. Under some simplifying assumptions these systems can be modeled by thin Ising metamagnetic films. In this paper we study, using both the Wang/Landau scheme and importance sampling Monte Carlo simulations, the equilibrium properties of these films. On the one hand we discuss the microcanonical density of states and its prominent features. On the other we analyze canonically various global and layer quantities. We obtain the phase diagram of thin Ising metamagnets as a function of temperature and external magnetic field. Whereas the phase diagram of the bulk system only exhibits one phase transition between the antiferromagnetic and paramagnetic phases, the phase diagram of thin Ising metamagnets includes an additional intermediate phase where one of the surface layers has aligned itself with the direction of the applied magnetic field. This additional phase transition is discontinuous and ends in a critical end point. Consequently, it is possible to gradually go from the antiferromagnetic phase to the intermediate phase without passing through a phase transition.


title: Anti-de-Sitter spacetime and its uses

abstract: This is a pedagogic account of some of the global properties of Anti-de-Sitter spacetime with a view to their application to the AdS/CFT correspondence. Particular care is taken over the distinction between Anti-de-Sitter and it's covering space. Written version of lectures given at 2nd Samos Meeting held at at Pythagoreon, Samos, Greece, 31 August - 4 September 1998 and published as Anti-de-Sitter spacetime and its uses, in Mathematical and Quantum Aspects of Relativity and Cosmology. Proceedings of the 2nd Samos Meeting on Cosmology, Geometry and Relativity, S Cotsakis and G W Gibbons eds, {\it Lecture Notes in Physics}\, {\bf 537} (2000)}


title: Comparative Wakefield Analysis of a First Prototype of a DDS Structure for CLIC Main Linac

abstract: A Damped Detuned Structure (DDS) for CLIC main linac has been proposed as an alternative to the present baseline design which is based on heavy damping. A first prototype, CLIC_DDS_A, for high power tests has been already designed and is under construction. It is also foreseen to design a further prototype, CLIC_DDS_B, to test both the wakefield suppression and high power performances. Wakefield calculations for DDS are, in the early design stage, based on single infinitely periodic cells. Though cell-to-cell interaction is taken into account to calculate the wakefields, it is important to study full structure properties using computational tools. In particular this is fundamental for defining the input parameters for the HOM coupler that is crucial for the performances of DDS. In the following a full analysis of wakefields and impedances based on simulations conducted with finite difference based electromagnetic computer code GdfidL will be presented.


title: Rotation, Scaling and Translation Analysis of Biometric Signature Templates

abstract: Biometric authentication systems that make use of signature verification methods often render optimum performance only under limited and restricted conditions. Such methods utilize several training samples so as to achieve high accuracy. Moreover, several constraints are imposed on the end-user so that the system may work optimally, and as expected. For example, the user is made to sign within a small box, in order to limit their signature to a predefined set of dimensions, thus eliminating scaling. Moreover, the angular rotation with respect to the referenced signature that will be inadvertently introduced as human error, hampers performance of biometric signature verification systems. To eliminate this, traditionally, a user is asked to sign exactly on top of a reference line. In this paper, we propose a robust system that optimizes the signature obtained from the user for a large range of variation in Rotation-Scaling-Translation (RST) and resolves these error parameters in the user signature according to the reference signature stored in the database.


title: Audio Watermarking with Error Correction

abstract: In recent times, communication through the internet has tremendously facilitated the distribution of multimedia data. Although this is indubitably a boon, one of its repercussions is that it has also given impetus to the notorious issue of online music piracy. Unethical attempts can also be made to deliberately alter such copyrighted data and thus, misuse it. Copyright violation by means of unauthorized distribution, as well as unauthorized tampering of copyrighted audio data is an important technological and research issue. Audio watermarking has been proposed as a solution to tackle this issue. The main purpose of audio watermarking is to protect against possible threats to the audio data and in case of copyright violation or unauthorized tampering, authenticity of such data can be disputed by virtue of audio watermarking.


title: Tracking Down the Highest Spindown Power Gamma-ray Pulsars

abstract: Forty six gamma-ray pulsars were reported in the First Fermi Large Area Telescope (LAT) Catalog of Gamma-ray Pulsars. Over forty more have been seen since then. A simple but effective figure-of-merit for gamma-detectability is sqrt(Edot)/d^2, where Edot is the pulsar spindown power and d the distance. We are tracking down the best gamma-ray candidates not yet seen. We present the timing and spectral analysis results of some new high spindown power, nearby gamma-ray pulsars. We also update some population distribution plots in preparation for the 2nd Fermi LAT gamma-ray Pulsar Catalog.


title: Collisional energy losses of heavy quarks in relativistic nuclear collisions at energies available at the BNL Relativistic Heavy Ion Collider

abstract: We investigate the energy losses of $b$ and $c$ quarks at hard collisions with gluons and light quarks in quark-gluon plasma produced in central $Au + Au$ collisions at RHIC energy($\sqrt{s_{NN}} = 200$ GeV). As the coupling constant at hard collision in the perturbative QCD is limited by Landau pole at momenta $Q^{2} = \Lambda^{2}$, we use the analytic models of QCD, where the Landau pole is absent. We calculate the energy losses of $b$ and $c$ quarks in analytic model of QCD at $Q^{2} \ge \Lambda^{2}$. At calculations we use the effective quasiparticle model. We use at calculations of total nuclear modification factor $R_{AA}^{b+c}$ the fragmentation of $b$ and $c$ quarks into heavy hadrons with accounting theirs energy losses. We show that this factor gives the considerably suppression of heavy $B$ and $D$ mesons in the region of middle meanings of $p_{\perp}$. The experimental data \cite{2} show, that nonphotonic electrons have the analogous suppression in the same region of $p_{\perp}$.


title: Chandra Publication Statistics

abstract: In this study we develop and propose publication metrics, based on an analysis of data from the Chandra bibliographic database, that are more meaningful and less sensitive to observatory-specific characteristics than the traditional metrics. They fall in three main categories: speed of publication; fraction of observing time published; and archival usage. Citation of results is a fourth category, but lends itself less well to definite statements. For Chandra, the median time from observation to publication is 2.36 years; after about 7 years 90% of the observing time is published; after 10 years 70% of the observing time is published more than twice; and the total annual publication output of the mission is 60-70% of the cumulative observing time available, assuming a two year lag between data retrieval and publication.


title: Hom and FP Coupler Design for the NLSF High Gradient SC Cavity

abstract: The design of both higher order mode (HOM) and fundamental power (FP) couplers for a unique New Low Surface Field (NLSF) cavity is presented. Here we present a study which uses ILC baseline couplers. The Balleyguier method of calculating external quality factor is used and the results validated using both Microwave studio and HFSS.


title: Long Horizons, High Risk Aversion, and Endogeneous Spreads

abstract: For an investor with constant absolute risk aversion and a long horizon, who trades in a market with constant investment opportunities and small proportional transaction costs, we obtain explicitly the optimal investment policy, its implied welfare, liquidity premium, and trading volume. We identify these quantities as the limits of their isoelastic counterparts for high levels of risk aversion. The results are robust with respect to finite horizons, and extend to multiple uncorrelated risky assets. In this setting, we study a Stackelberg equilibrium, led by a risk-neutral, monopolistic market maker who sets the spread as to maximize profits. The resulting endogenous spread depends on investment opportunities only, and is of the order of a few percentage points for realistic parameter values.


title: Sunspots are in many ways similar to terrestrial vortices

abstract: In this letter we identify similarities amongst sunspots and terrestrial vortices. The dark appearance of the central part of any sunspot is currently justified by an anticipated cooling effect experienced by the ionized gas. However, it cannot single-handedly reconcile the halo that surrounds the penumbra, the subsequent second dim ring that could be possibly followed by a second halo. In antithesis, light refraction due to density variations in a compressible whirl can give reason for all of these manifestations. Certain data of Wilson's depression fit better the geometric depth profile of a two-celled vortex. The last provides a hurricane equivalent manifestation for the normal and reverse Evershed effect. There is compelling evidence that alike to atmospheric vortices sunspots do also spawn meso-cyclones.


title: Higher Order Modes in Coupled Cavities of the Flash Module ACC39

abstract: We analyse the higher order modes (HOM's) in the 3.9GHz bunch shaping cavities installed in the FLASH facility at DESY. A suite of finite element computer codes (including HFSS and ACE3P) and globalised scattering matrix calculations (GSM) are used to investigate the modes in these cavities. This study is primarily focused on the dipole component of the multiband expansion of the wakefield, with the emphasis being on the development of a HOM-based BPM system for ACC39. Coupled inter-cavity modes are simulated together with a limited band of trapped modes.


title: Production and characterization of a custom-made 228Th source with reduced neutron source strength for the Borexino experiment

abstract: A custom-made 228Th source of several MBq activity was produced for the Borexino experiment for studying the external background of the detector. The aim was to reduce the unwanted neutron emission produced via (alpha,n) reactions in ceramics used typically for commercial 228Th sources. For this purpose a ThCl4 solution was converted chemically into ThO2 and embedded into a gold foil. The paper describes the production and the characterization of the custom-made source by means of gamma-activity, dose rate and neutron source strength measurements. From gamma-spectroscopic measurements it was deduced that the activity transfer from the initial solution to the final source was >91% (at 68% C.L.) and the final activity was (5.41+-0.30) MBq. The dose rate was measured by two dosimeters yielding 12.1 mSv/h and 14.3 mSv/h in 1 cm distance. The neutron source strength of the 5.41 MBq 228Th source was determined as (6.59+-0.85)/sec.


title: PT-symmetric quantum models living in an auxiliary Pontryagin space

abstract: An extension of the scope of quantum theory is proposed in a way inspired by the recent heuristic as well as phenomenological success of the use of non-Hermitian Hamiltonians which are merely required self-adjoint in a Krein space with an indefinite metric (chosen, usually, as the operator of parity). In nuce, the parity-like operators are admitted to represent the mere indefinite metric in a Pontryagin space. A constructive version of such a generalized quantization strategy is outlined and, via a non-numerical illustrative example, found feasible.


title: Describing West-3-stack-sortable permutations with permutation patterns

abstract: We describe a new method for finding patterns in permutations that produce a given pattern after the permutation has been passed once through a stack. We use this method to describe West-3-stack-sortable permutations, that is, permutations that are sorted by three passes through a stack. We also show how the method can be applied to the bubble-sort operator. The method requires the use of mesh patterns introduced by Br\"and\'en and Claesson (2011), as well as a new type of generalized pattern we call a decorated pattern.


title: Standard Quantum Teleportation of an Arbitrary N-Qubit State, Non-Existence of Magic Basis and Existence of Magic Partial Bases for 2N Entangled Qubit States with N>1

abstract: We present a simple and precise protocol for standard quantum teleportation of N-qubit state, considering the most general resource q-channel and Bell states. We find condition on these states for perfect teleportation and give explicitly the unitary transformation required to be done by Bob for achieving perfect teleportation. We discuss connection of our simple theory with the complicated related work on this subject and with character matrix, transformation, judgment and kernel operators defined in this context. We also prove that the magic basis discussed by Hill and Wootters [Phys. Rev. Lett. 78 (1997) 5022] does not exist for entangled 2N-qubit states with N > 1 but magic partial bases, similar to those discussed recently by Prakash and Maurya [Optics Commun. 284 (2011) 5024] do exist. We give explicitly all magic partial bases for N = 2.


title: Feasible Automata for Two-Variable Logic with Successor on Data Words

abstract: We introduce an automata model for data words, that is words that carry at each position a symbol from a finite alphabet and a value from an unbounded data domain. The model is (semantically) a restriction of data automata, introduced by Bojanczyk, et. al. in 2006, therefore it is called weak data automata. It is strictly less expressive than data automata and the expressive power is incomparable with register automata. The expressive power of weak data automata corresponds exactly to existential monadic second order logic with successor +1 and data value equality \sim, EMSO2(+1,\sim). It follows from previous work, David, et. al. in 2010, that the nonemptiness problem for weak data automata can be decided in 2-NEXPTIME. Furthermore, we study weak B\"uchi automata on data omega-strings. They can be characterized by the extension of EMSO2(+1,\sim) with existential quantifiers for infinite sets. Finally, the same complexity bound for its nonemptiness problem is established by a nondeterministic polynomial time reduction to the nonemptiness problem of weak data automata.


title: From travelling waves to mild chaos: a supercritical bifurcation cascade in pipe flow

abstract: We study numerically a succession of transitions in pipe Poiseuille flow that leads from simple travelling waves to waves with chaotic time-dependence. The waves at the origin of the bifurcation cascade possess a shift-reflect symmetry and are both axially and azimuthally periodic with wave numbers {\kappa} = 1.63 and n = 2, respectively. As the Reynolds number is increased, successive transitions result in a wide range of time dependent solutions that includes spiralling, modulated-travelling, modulated-spiralling, doubly-modulated-spiralling and mildly chaotic waves. We show that the latter spring from heteroclinic tangles of the stable and unstable invariant manifolds of two shift-reflect-symmetric modulated-travelling waves. The chaotic set thus produced is confined to a limited range of Reynolds numbers, bounded by the occurrence of manifold tangencies. The states studied here belong to a subspace of discrete symmetry which makes many of the bifurcation and path-following investigations presented technically feasible. However, we expect that most of the phenomenology carries over to the full state-space, thus suggesting a mechanism for the formation and break-up of invariant states that can sustain turbulent dynamics.


title: Exact deflection of a Neutral-Tachyon in the Kerr's Gravitational field

abstract: We solve in closed analytic form space-like geodesic equations in the Kerr gravitational field. Such geodesic equations describe the motion of neutral tachyons (faster than light particles) in the Kerr spacetime. More specifically we derive the closed form solution for the deflection angle of a neutral tachyon on an equatorial orbit in Kerr spacetime. The solution is expressed elegantly in terms of Lauricella's hypergeometric function F_{D}.We applied our results to three cases: first, for the calculation of the deflection angle of a neutral tachyon on an equatorial trajectory in the gravitational field of a Kerr black hole. Subsequently, we applied our exact solutions to compute the deflection angle of equatorial spacelike geodesics in the gravitational fields of Sun and Earth assuming the Kerr spacetime geometry.


title: Sections of surface bundles and Lefschetz fibrations

abstract: We investigate the possible self-intersection numbers for sections of surface bundles and Lefschetz fibrations over surfaces. When the fiber genus g and the base genus h are positive, we prove that the adjunction bound 2h-2 is the only universal bound on the self-intersection number of a section of any such genus g bundle and fibration. As a side result, in the mapping class group of a surface with boundary, we calculate the precise value of the commutator lengths of all powers of a Dehn twist about a boundary component, concluding that the stable commutator length of such a Dehn twist is 1/2. We furthermore prove that there is no upper bound on the number of critical points of genus-g Lefschetz fibrations over surfaces with positive genera admitting sections of maximal self-intersection, for g at least two.


title: Spin and pseudospin symmetry along with orbital dependency of the Dirac-Hulthen problem

abstract: The role of the Hulthen potential on the spin and pseudospin symmetry solutions is investigated systematically by solving the Dirac equation with attractive scalar S(r) and repulsive vector V(r) potentials. The spin and pseudospin symmetry along with orbital dependency (pseudospin-orbit and spin-orbit dependent couplings) of the Dirac equation are included to the solution by introducing the Hulthen-square approximation. This effective approach is based on forming the spin and pseudo-centrifugal kinetic energy term from the square of the Hulthen potential. The analytical solutions of the Dirac equation for the Hulthen potential with the spin-orbit and pseudospin-orbit-dependent couplings are obtained by using the Nikiforov-Uvarov (NU) method. The energy eigenvalue equations and wave functions for various degenerate states are presented for several spin-orbital, pseudospin-orbital and radial quantum numbers under the condition of the spin and pseudospin symmetry. Keywords: Spin and pseudospin symmetry; orbital dependency; Dirac equation; Hulthen potential; Nikiforov-Uvarov Method.


title: Absorption features caused by oscillations of electrons on the surface of a quark star

abstract: If quark stars exist, they may be enveloped in thin electron layers (electron seas), which uniformly surround the entire star. These layers will be affected by the magnetic fields of quark stars in such a way that the electron seas would transmit hydromagnetic cyclotron waves, as studied in this paper. Particular attention is devoted to vortex hydrodynamical oscillations of the electron sea. The frequency spectrum of these oscillations is derived in analytic form. If the thermal X-ray spectra of quark stars are modulated by vortex hydrodynamical vibrations, the thermal spectra of compact stars, foremost central compact objects (CCOs) and X-ray dim isolated neutron stars (XDINSs), could be used to verify the existence of these vibrational modes observationally. The central compact object 1E 1207.4-5209 appears particularly interesting in this context, since its absorption features at 0.7 keV and 1.4 keV can be comfortably explained in the framework of the hydro-cyclotron oscillation model.


title: Odd H-depth and H-separable extensions

abstract: A subring pair B < A has right depth 2n if the n+1'st relative Hochschild bar resolution group is isomorphic to a direct summand of a multiple of the n'th relative Hochschild bar resolution group as A-B-bimodules; depth 2n+1 if the same condition holds only as B-B-bimodules. It is then natural to ask what is defined if this same condition should hold as A-A-bimodules, the so-called H-depth 2n-1 condition. In particular, the H-depth 1 condition coincides with A being an H-separable extension of B. In this paper the H-depth of semisimple subalgebra pairs is derived from the transpose inclusion matrix, and for QF extensions it is derived from the odd depth of the endomorphism ring extension. For general extensions characterizations of H-depth are possible using the H-equivalence generalization of Morita theory. In certain nice categories of bimodules the minimum depth and H-depth of certain types of ring extensions are always finite.


title: Order-distance and other metric-like functions on jointly distributed random variables

abstract: We construct a class of real-valued nonnegative binary functions on a set of jointly distributed random variables, which satisfy the triangle inequality and vanish at identical arguments (pseudo-quasi-metrics). These functions are useful in dealing with the problem of selective probabilistic causality encountered in behavioral sciences and in quantum physics. The problem reduces to that of ascertaining the existence of a joint distribution for a set of variables with known distributions of certain subsets of this set. Any violation of the triangle inequality or its consequences by one of our functions when applied to such a set rules out the existence of this joint distribution. We focus on an especially versatile and widely applicable pseudo-quasi-metric called an order-distance and its special case called a classification distance.


title: BRST-Invariant Deformations of Geometric Structures in Sigma Models

abstract: We study a Lie algebra of formal vector fields $W_n$ with its application to the perturbative deformed holomorphic symplectic structure in the A-model, and a Calabi-Yau manifold with boundaries in the B-model. We show that equivalent classes of deformations are describing by a Hochschild cohomology theory of the DG-algebra ${\mathfrak A} = (A, Q)$, $Q =\bar{\partial}+\partial_{\rm deform}$, which is defined to be the cohomology of $(-1)^n Q +d_{\rm Hoch}$. Here $\bar{\partial}$ is the initial non-deformed BRST operator while $\partial_{\rm deform}$ is the deformed part whose algebra is a Lie algebra of linear vector fields ${\rm gl}_n$. We show that equivalent classes of deformations are described by a Hochschild cohomology of ${\mathfrak A}$, an important geometric invariant of the (anti)holomorphic structure on $X$. We discuss the identification of the harmonic structure $(HT^\bullet(X); H\Omega_\bullet(X))$ of affine space $X$ and the group ${\rm Ext}_{X^{2}}^n({\cO}_{\triangle}, {\cO}_{\triangle})$ (the HKR isomorphism), and bulk-boundary deformation pairing.


title: Solitary and Coupled Semiconductor Ring Lasers as Optical Spiking Neurons

abstract: We theoretically investigate the possibility of generating pulses in an excitable (asymmetric) semiconductor ring laser (SRL) using optical trigger pulses. We show that the phase difference between the injected field and the electric field inside the SRL determines the direction of the perturbation in phase space. Due to the folded shape of the excitability threshold, this has an important influence on the ability to cross it. A mechanism for exciting multiple consecutive pulses using a single trigger pulse (i.e. multi pulse excitability) is revealed. We furthermore investigate the possibility of using asymmetric SRLs in a coupled configuration, which is a first step toward an all-optical neural network using SRLs as building blocks.


title: Aoki Phases in Staggered-Wilson Fermions

abstract: We investigate the parity-broken phase (Aoki phase) for staggered-Wilson fermions by using the Gross-Neveu model and the strong-coupling lattice QCD. In the both cases the gap equations indicate the parity-broken phase exists and the pion becomes massless on the phase boundaries. We also show we can take the chiral and continuum limit in the Gross-Neveu model by tuning mass and gauge-coupling parameters. This supports the idea that the staggered-Wilson fermions can be applied to the lattice QCD simulation by taking a chiral limit, as with Wilson fermions.


title: Some examples of Baker domains

abstract: We construct entire functions with hyperbolic and simply parabolic Baker domains on which the functions are not univalent. The Riemann maps from the unit disk to these Baker domains extend continuously to certain arcs on the unit circle. The results answer questions posed by Fagella and Henriksen, Baker and Dominguez, and others.


title: Path properties of dilatively stable processes and singularity of their distributions

abstract: First, we present some results about the H\"older continuity of the sample paths of so called dilatively stable processes which are certain infinitely divisible processes having a more general scaling property than self-similarity. As a corollary, we obtain that the most important (H,delta)-dilatively stable limit processes (e.g., the LISOU and the LISCBI processes, see Igloi [4]) almost surely have a local H\"older exponent H. Next we prove that, under some slight regularity assumptions, any two dilatively stable processes with stationary increments are singular (in the sense that their distributions have disjoint supports) if their parameters H are different. We also study the more general case of not having stationary increments. Throughout the paper we specialize our results to some basic dilatively stable processes such as the above-mentioned limit processes and the fractional L\'evy process.


title: Dissipative and conservative nonlinearity in carbon nanotube and graphene mechanical resonators

abstract: Graphene and carbon nanotubes represent the ultimate size limit of one and two-dimensional nanoelectromechanical resonators. Because of their reduced dimensionality, graphene and carbon nanotubes display unusual mechanical behavior; in particular, their dynamics is highly nonlinear. Here, we review several types of nonlinear behavior in resonators made from nanotubes and graphene. We first discuss an unprecedented scenario where damping is described by a nonlinear force. This scenario is supported by several experimental facts: (i) the quality factor varies with the amplitude of the motion as a power law whose exponent coincides with the value predicted by the nonlinear damping model, (ii) hysteretic behavior (of the motional amplitude as a function of driving frequency) is absent in some of our resonators even for large driving forces, as expected when nonlinear damping forces are large, and (iii) when we quantify the linear damping force (by performing parametric excitation measurements) we find that it is significantly smaller than the nonlinear damping force. We then review parametric excitation measurements, an alternative actuation method which is based on nonlinear dynamics. Finally, we discuss experiments where the mechanical motion is coupled to electron transport through a nanotube. The coupling can be made so strong that the associated force acting on the nanotube becomes highly nonlinear with displacement and velocity. Overall, graphene and nanotube resonators hold promise for future studies on classical and quantum nonlinear dynamics.


title: Single-walled carbon nanotubes and nanocrystalline graphene reduce beam-induced movements in high-resolution electron cryo-microscopy of ice-embedded biological samples

abstract: For single particle electron cryo-microscopy (cryoEM), contrast loss due to beam-induced charging and specimen movement is a serious problem, as the thin films of vitreous ice spanning the holes of a holey carbon film are particularly susceptible to beam-induced movement. We demonstrate that the problem is at least partially solved by carbon nanotechnology. Doping ice-embedded samples with single-walled carbon nanotubes (SWNT) in aqueous suspension or adding nanocrystalline graphene supports, obtained by thermal conversion of cross-linked self-assembled biphenyl precursors, significantly reduces contrast loss in high-resolution cryoEM due to the excellent electrical and mechanical properties of SWNTs and graphene.


title: Spectroscopy of $^{18}$Na: Bridging the two-proton radioactivity of $^{19}$Mg

abstract: The unbound nucleus $^{18}$Na, the intermediate nucleus in the two-proton radioactivity of $^{19}$Mg, was studied by the measurement of the resonant elastic scattering reaction $^{17}$Ne(p,$^{17}$Ne)p performed at 4 A.MeV. Spectroscopic properties of the low-lying states were obtained in a R-matrix analysis of the excitation function. Using these new results, we show that the lifetime of the $^{19}$Mg radioactivity can be understood assuming a sequential emission of two protons via low energy tails of $^{18}$Na resonances.


title: Free Deterministic Equivalents, Rectangular Random Matrix Models, and Operator-Valued Free Probability Theory

abstract: Motivated by the asymptotic collective behavior of random and deterministic matrices, we propose an approximation (called "free deterministic equivalent") to quite general random matrix models, by replacing the matrices with operators satisfying certain freeness relations. We comment on the relation between our free deterministic equivalent and deterministic equivalents considered in the engineering literature. We do not only consider the case of square matrices, but also show how rectangular matrices can be treated. Furthermore, we emphasize how operator-valued free probability techniques can be used to solve our free deterministic equivalents. As an illustration of our methods we consider a random matrix model studied first by R. Couillet, J. Hoydis, and M. Debbah. We show how its free deterministic equivalent can be treated and we thus recover in a conceptual way their result. On a technical level, we generalize a result from scalar valued free probability, by showing that randomly rotated deterministic matrices of different sizes are asymptotically free from deterministic rectangular matrices, with amalgamation over a certain algebra of projections. In the Appendix, we show how estimates for differences between Cauchy transforms can be extended from a neighborhood of infinity to a region close to the real axis. This is of some relevance if one wants to compare the original random matrix problem with its free deterministic equivalent.


title: Transport Coefficients and nPI Methods

abstract: Transport coefficients can be obtained from 2-point correlators using the Kubo formulae. It has been shown that the full leading order result for electrical conductivity and (QCD) shear viscosity is contained in the re-summed 2-point function that is obtained from the 3-loop 3PI effective action. The theory produces all leading order contributions without the necessity for power counting, and in this sense it provides a natural framework for the calculation and suggests that one can calculate the next-to-leading contribution to transport coefficients from the 4-loop 4PI effective action. The integral equations have been derived for shear viscosity for a scalar theory with cubic and quartic interactions, with a non-vanishing field expectation value. We review these results, and explain how the calculation could be done at higher orders.


title: Entanglement entropy in all dimensions

abstract: It has long been conjectured that the entropy of quantum fields across boundaries scales as the boundary area. This conjecture has not been easy to test in spacetime dimensions greater than four because of divergences in the von Neumann entropy. Here we show that the Renyi entropy provides a convergent alternative, yielding a quantitative measure of entanglement between quantum field theoretic degrees of freedom inside and outside hypersurfaces. For the first time, we show that the entanglement entropy in higher dimensions is proportional to the higher dimensional area. We also show that the Renyi entropy diverges at specific values of the Renyi parameter (q) in each dimension, but this divergence can be tamed by introducing a mass of the quantum field.


title: On Graphs with the Smallest Eigenvalue at Least $-1-\sqrt{2}$, part II

abstract: This is a continuation of the article with the same title. In this paper, the family H is the same as in the previous paper "On Graphs with the Smallest Eigenvalue at Least $-1-\sqrt{2}$, part I". The main result is that a minimal graph which is not an H -line graph, is just isomorphic to one of the 38 graphs found by computer.


title: Critical Casimir effect in films for generic non-symmetry-breaking boundary conditions

abstract: Systems described by an O(n) symmetrical $\phi^4$ Hamiltonian are considered in a $d$-dimensional film geometry at their bulk critical points. A detailed renormalization-group (RG) study of the critical Casimir forces induced between the film's boundary planes by thermal fluctuations is presented for the case where the O(n) symmetry remains unbroken by the surfaces. The boundary planes are assumed to cause short-ranged disturbances of the interactions that can be modelled by standard surface contributions $\propto \bm{\phi}^2$ corresponding to subcritical or critical enhancement of the surface interactions. This translates into mesoscopic boundary conditions of the generic symmetry-preserving Robin type $\partial_n\bm{\phi}=\mathring{c}_j\bm{\phi}$. RG-improved perturbation theory and Abel-Plana techniques are used to compute the $L$-dependent part $f_{\mathrm{res}}$ of the reduced excess free energy per film area $A\to\infty $ to two-loop order. When $d<4$, it takes the scaling form $f_{\mathrm{res}}\approx D(c_1L^{\Phi/\nu},c_2L^{\Phi/\nu})/L^{d-1}$ as $L\to\infty$, where $c_i$ are scaling fields associated with the surface-enhancement variables $\mathring{c}_i$, while $\Phi$ is a standard surface crossover exponent. The scaling function $D(\mathsf{c}_1,\mathsf{c}_2)$ and its analogue $\mathcal{D}(\mathsf{c}_1,\mathsf{c}_2)$ for the Casimir force are determined via expansion in $\epsilon=4-d$ and extrapolated to $d=3$ dimensions. In the special case $\mathsf{c}_1=\mathsf{c}_2=0$, the expansion becomes fractional. Consistency with the known fractional expansions of D(0,0) and $\mathcal{D}(0,0)$ to order $\epsilon^{3/2}$ is achieved by appropriate reorganisation of RG-improved perturbation theory. For appropriate choices of $c_1$ and $c_2$, the Casimir forces can have either sign. Furthermore, crossovers from attraction to repulsion and vice versa may occur as $L$ increases.


title: Tunnel spin injection into graphene using Al2O3 barrier grown by atomic layer deposition on functionalized graphene surface

abstract: We demonstrate electrical tunnel spin injection from a ferromagnet to graphene through a high-quality Al2O3 grown by atomic layer deposition (ALD). The graphene surface is functionalized with a self-assembled monolayer of 3,4,9,10-perylene tetracarboxylic acid (PTCA) to promote adhesion and growth of Al2O3 with a smooth surface. Using this composite tunnel barrier of ALD-Al2O3 and PTCA, a spin injection signal of 30 ohm has been observed from non-local magnetoresistance measurements at 45 K, revealing potentially high performance of ALD-Al2O3/PTCA tunnel barrier for spin injection into graphene.


title: Quantitative Dunford-Pettis property

abstract: We investigate possible quantifications of the Dunford-Pettis property. We show, in particular, that the Dunford-Pettis property is automatically quantitative in a sense. Further, there are two incomparable mutually dual stronger versions of a quantitative Dunford-Pettis property. We prove that $L^1$ spaces and $C(K)$ spaces posses both of them. We also show that several natural measures of weak non-compactness are equal in $L^1$ spaces.


title: Energy-filtered transmission electron microscopy of biological samples on highly transparent carbon nanomembranes

abstract: Ultrathin carbon nanomembranes (CNM) comprising crosslinked biphenyl precursors have been tested as support films for energy-filtered transmission electron microscopy (EFTEM) of biological specimens. Due to their high transparency CNM are ideal substrates for electron energy loss spectroscopy (EELS) and electron spectroscopic imaging (ESI) of stained and unstained biological samples. Virtually background-free elemental maps of tobacco mosaic virus (TMV) and ferritin have been obtained from samples supported by ~ 1 nm thin CNM. Furthermore, we have tested conductive carbon nanomembranes (cCNM) comprising nanocrystalline graphene, obtained by thermal treatment of CNM, as supports for cryoEM of ice-embedded biological samples. We imaged ice-embedded TMV on cCNM and compared the results with images of ice-embedded TMV on conventional carbon film (CC), thus analyzing the gain in contrast for TMV on cCNM in a quantitative manner. In addition we have developed a method for the preparation of vitrified specimens, suspended over the holes of a conventional holey carbon film, while backed by ultrathin cCNM.


title: Getting routers out of the core: Building an optical wide area network with "multipaths"

abstract: We propose an all-optical networking solution for a wide area network (WAN) based on the notion of multipoint-to-multipoint lightpaths that, for short, we call "multipaths". A multipath concentrates the traffic of a group of source nodes on a wavelength channel using an adapted MAC protocol and multicasts this traffic to a group of destination nodes that extract their own data from the confluent stream. The proposed network can be built using existing components and appears less complex and more efficient in terms of energy consumption than alternatives like OPS and OBS. The paper presents the multipath architecture and compares its energy consumption to that of a classical router-based ISP network. A flow-aware dynamic bandwidth allocation algorithm is proposed and shown to have excellent performance in terms of throughput and delay.


title: Orientation-dependent ionization yields from strong-field ionization of fixed-in-space linear and asymmetric top molecules

abstract: The yield of strong-field ionization, by a linearly polarized probe pulse, is studied experimentally and theoretically, as a function of the relative orientation between the laser field and the molecule. Experimentally, carbonyl sulfide, benzonitrile and naphthalene molecules are aligned in one or three dimensions before being singly ionized by a 30 fs laser pulse centered at 800 nm. Theoretically, we address the behaviour of these three molecules. We consider the degree of alignment and orientation and model the angular dependence of the total ionization yield by molecular tunneling theory accounting for the Stark shift of the energy level of the ionizing orbital. For naphthalene and benzonitrile the orientational dependence of the ionization yield agrees well with the calculated results, in particular the observation that ionization is maximized when the probe laser is polarized along the most polarizable axis. For OCS the observation of maximum ionization yield when the probe is perpendicular to the internuclear axis contrasts the theoretical results.


title: Anomalous Molecular Weight Dependence of Chain Dynamics in Unentangled Polymer Blends with Strong Dynamic Asymmetry

abstract: We address the general question of how the molecular weight dependence of chain dynamics in unentangled polymers is modified by blending. By dielectric spectroscopy we measure the normal mode relaxation of polyisoprene in blends with a slow matrix of poly(ter-butylstyrene). Unentangled polyisoprene in the blend exhibits strong deviations from Rouse scaling, approaching 'entangled-like' behavior at low temperatures in concomitance with the increase of the dynamic asymmetry in the blend. The obtained results are discussed in the framework of the generalized Langevin equation formalism. On this basis, a non trivial relationship between the molecular weight dependence of the longest chain relaxation time and the nonexponentiality of the corresponding Rouse correlator is found. This result is confirmed by molecular dynamics simulations.


title: An algorithm to compute the power of Monte Carlo tests with guaranteed precision

abstract: This article presents an algorithm that generates a conservative confidence interval of a specified length and coverage probability for the power of a Monte Carlo test (such as a bootstrap or permutation test). It is the first method that achieves this aim for almost any Monte Carlo test. Previous research has focused on obtaining as accurate a result as possible for a fixed computational effort, without providing a guaranteed precision in the above sense. The algorithm we propose does not have a fixed effort and runs until a confidence interval with a user-specified length and coverage probability can be constructed. We show that the expected effort required by the algorithm is finite in most cases of practical interest, including situations where the distribution of the p-value is absolutely continuous or discrete with finite support. The algorithm is implemented in the R-package simctest, available on CRAN.


title: On $r$-colorability of random hypergraphs

abstract: The work deals with the threshold probablity for r-colorability in the binomial model H(n,k,p) of a random k-uniform hypergraph. We prove a lower bound for this threshold which improves the previously known results in the wide range of the parameters r=r(n) and k=k(n).


title: The first complete NLL BFKL study of Mueller Navelet jets at LHC

abstract: We report on the first next-to-leading BFKL study of the cross section and azimuthal decorrellation of Mueller Navelet jets. This includes next-to-leading corrections to the Green's function as well as next-to-leading corrections to the Mueller Navelet vertices. The obtained results for standard observables proposed for studies of Mueller Navelet jets show that both sources of corrections are of equal and big importance for final magnitude and final behavior of observables, in particular for the LHC kinematics investigated here in detail. The astonishing conclusion of our analysis is that the observables obtained within the complete next-lo-leading order BFKL framework of the present contribution are quite similar to the same observables obtained within next-to-leading logarithm DGLAP type treatment. The only noticeable difference is the ratio the azimuthal angular moments < cos 2 phi> / < cos phi > which still differs in both treatments.


title: NLL QCD contribution of the Electromagnetic Dipole operator to B-> X_s gamma gamma

abstract: We calculate the set of O(alpha_s) corrections to the double differential decay width dGamma_77/(ds_1*ds_2) for the process B -> X_s gamma gamma originating from diagrams involving the electromagnetic dipole operator O_7. The kinematical variables s_1 and s_2 are defined as s_i=(p_b - q_i)^2/m_b^2, where p_b, q_1, q_2 are the momenta of b-quark and two photons. While the (renormalized) virtual corrections are worked exactly for a certain range of s_1 and s_2, we retain in the gluon bremsstrahlung process only the leading power w.r.t. the (normalized) hadronic mass s_3=(p_b-q_1-q_2)^2/m_b^2 in the underlying triple differential decay width dGamma_77/(ds_1*ds_2*ds_3). The double differential decay width, based on this approximation is free of infrared- and collinear singularities when combining virtual- and bremsstrahlung corrections. The corresponding results are obtained analytically. When retaining all powers in s_3, the sum of virtual- and bremstrahlung corrections contains uncanceled 1/epsilon singularities (which are due to collinear photon emission from the s-quark) and other concepts, which go beyond perturbation theory, like parton fragmentation functions of a quark or a gluon into a photon, are needed which is beyond the scope of our paper.


title: A Multi-Wavelength Study of the Unidentified TeV Gamma-ray Source HESS J1626-490

abstract: HESS J1626-490, so far only detected with the H.E.S.S. array of imaging atmospheric Cherenkov telescopes, could not be unambiguously identified with any source seen at lower energies. Therefore, we analyzed data from an archival XMM-Newton observation, pointed towards HESS J1626-490, to classify detected point-like and extended X-ray sources according to their spectral properties. None of the detected X-ray point sources fulfills the energetic requirements to be considered as the synchrotron radiation (SR) counterpart to the VHE source assuming an Inverse Compton (IC) emission scenario. Furthermore, we did not detect any diffuse X-ray excess emission originating from the region around HESS J1626-490 above the Galactic Background. The derived upper limit for the total X-ray flux disfavors a purely leptonic emission scenario for HESS J1626-490. To characterize the Interstellar Medium surrounding HESS J1626-490 we analyzed ^{12}CO(J=1-0) molecular line data from the Nanten Galactic plane survey and HI data from the Southern Galactic Plane Survey (SGPS). We found a good morphological match between molecular and atomic gas in the -27 km/s to -18 km/s line-of-sight velocity range and HESS J1626-490. The cloud has a mass of 1.8x10^4 M_sun and is located at a mean kinematic distance of d = 1.8 kpc. Furthermore, we found a density depression in the HI gas at a similar distance which is spatially consistent with the SNR G335.2+00.1. Therefore, the most likely origin of the VHE gamma-ray emission observed with H.E.S.S. is the hadronic interaction of cosmic rays with a moderately dense molecular cloud, which we detected with Nanten. The application of a detailed hadronic model for cosmic ray transport and interaction shows that the cosmic rays could originate from the nearby SNR G335.2+00.1.


title: Superluminal group velocity of neutrinos

abstract: We suggest a possible interpretation of the recent observation by the OPERA collaboration of superluminal propagation of neutrinos. We show that it is in principle possible that the group velocity of neutrinos exceeds the speed of light without violating special relativity.


title: Random walks in cones

abstract: We study the asymptotic behavior of a multidimensional random walk in a general cone. We find the tail asymptotics for the exit time and prove integral and local limit theorems for a random walk conditioned to stay in a cone. The main step in the proof consists in constructing a positive harmonic function for our random walk under minimal moment restrictions on the increments. For the proof of tail asymptotics and integral limit theorems, we use a strong approximation of random walks by Brownian motion. For the proof of local limit theorems, we suggest a rather simple approach, which combines integral theorems for random walks in cones with classical local theorems for unrestricted random walks. We also discuss some possible applications of our results to ordered random walks and lattice path enumeration.


title: Nonlinear magneto-optical resonances for systems with J~100 observed in K2 molecules

abstract: We present the results of an experimental as well as theoretical study of nonlinear magneto-optical resonances in diatomic potassium molecules in the electronic ground state with large values of the angular momentum quantum number J~100. At zero magnetic field, the absorption transitions are suppressed because of population trapping in the ground state due to Zeeman coherences between magnetic sublevels of this state along with depopulation pumping. The destruction of such coherences in an external magnetic field was used to study the resonances in this work. K2 molecules were formed in a glass cell filled with potassium metal at a temperature above 150C. The cell was placed in an oven and was located in a homogeneous magnetic field B, which was scanned from zero to 0.7 T. Q-type and R-type transitions were excited with a tunable, single-mode diode laser with central wavelength of 660 nm. Well pronounced nonlinear Hanle effect signals were observed in the intensities of the linearly polarized components of the laser-induced fluorescence (LIF) detected in the direction parallel to the (B)-field with polarization vectors parallel (I_par) and perpendicular (I_per) to the polarization vector of the exciting laser radiation, which was orthogonal to (B). The intensities of the LIF components were detected for different experimental parameters, such as laser power density and vapor temperature, in order to compare them with numerical simulations that were based on the optical Bloch equations for the density matrix. We report good agreement of our measurements with numerical simulations. Narrow, subnatural line width dark resonances in I_per(B) were detected and explained.


title: Driving Rydberg-Rydberg transitions from a co-planar microwave waveguide

abstract: The coherent interaction between ensembles of helium Rydberg atoms and microwave fields in the vicinity of a solid-state co-planar waveguide is reported. Rydberg-Rydberg transitions, at frequencies between 25 GHz and 38 GHz, have been studied for states with principal quantum numbers in the range 30 - 35 by selective electric-field ionization. An experimental apparatus cooled to 100 K was used to reduce effects of blackbody radiation. Inhomogeneous, stray electric fields emanating from the surface of the waveguide have been characterized in frequency- and time-resolved measurements and coherence times of the Rydberg atoms on the order of 250 ns have been determined.


title: Neutrinos in Non-linear Structure Formation - a Simple SPH Approach

abstract: We present a novel method for implementing massive neutrinos in N-body simulations. Instead of sampling the neutrino velocity distribution by individual point particles we take neutrino free-streaming into account by treating it as an effective redshift dependent sound speed in a perfect isothermal fluid, and assume a relation between the sound speed and velocity dispersion of the neutrinos. Although the method fails to accurately model the true neutrino power spectrum, it is able to calculate the total matter power spectrum to the same accuracy as more complex hybrid neutrino methods, except on very small scales. We also present an easy way to update the publicly available Gadget-2 version with this neutrino approximation.


title: Estimation of drift and diffusion functions from time series data: A maximum likelihood framework

abstract: Complex systems are characterized by a huge number of degrees of freedom often interacting in a non-linear manner. In many cases macroscopic states, however, can be characterized by a small number of order parameters that obey stochastic dynamics in time. Recently techniques for the estimation of the corresponding stochastic differential equations from measured data have been introduced. This contribution develops a framework for the estimation of the functions and their respective (Bayesian posterior) confidence regions based on likelihood estimators. In succession approximations are introduced that significantly improve the efficiency of the estimation procedure. While being consistent with standard approaches to the problem this contribution solves important problems concerning the applicability and the accuracy of estimated parameters.


title: Characterizing and Improving Generalized Belief Propagation Algorithms on the 2D Edwards-Anderson Model

abstract: We study the performance of different message passing algorithms in the two dimensional Edwards Anderson model. We show that the standard Belief Propagation (BP) algorithm converges only at high temperature to a paramagnetic solution. Then, we test a Generalized Belief Propagation (GBP) algorithm, derived from a Cluster Variational Method (CVM) at the plaquette level. We compare its performance with BP and with other algorithms derived under the same approximation: Double Loop (DL) and a two-ways message passing algorithm (HAK). The plaquette-CVM approximation improves BP in at least three ways: the quality of the paramagnetic solution at high temperatures, a better estimate (lower) for the critical temperature, and the fact that the GBP message passing algorithm converges also to non paramagnetic solutions. The lack of convergence of the standard GBP message passing algorithm at low temperatures seems to be related to the implementation details and not to the appearance of long range order. In fact, we prove that a gauge invariance of the constrained CVM free energy can be exploited to derive a new message passing algorithm which converges at even lower temperatures. In all its region of convergence this new algorithm is faster than HAK and DL by some orders of magnitude.

categories: cond-mat.dis-nn cond-mat.stat-mech cs.AI cs.IT math.IT


title: Non-Standard Grain Properties, Dark Gas Reservoir, and Extended Submillimeter Excess, Probed by Herschel in the Large Magellanic Cloud

abstract: Aims: In this paper, we perform detailed modelling of the Spitzer and Herschel observations of the LMC, in order to: (i) systematically study the uncertainties and biases affecting dust mass estimates; and to (ii) explore the peculiar ISM properties of the LMC. Methods: To achieve these goals, we have modelled the spatially resolved SEDs with two alternate grain compositions, to study the impact of different submillimetre opacities on the dust mass. We have rigorously propagated the observational errors (noise and calibration) through the entire fitting process, in order to derive consistent parameter uncertainties. Results: First, we show that using the integrated SED leads to underestimating the dust mass by ~50 % compared to the value obtained with sufficient spatial resolution, for the region we studied. This might be the case, in general, for unresolved galaxies. Second, we show that Milky Way type grains produce higher gas-to-dust mass ratios than what seems possible according to the element abundances in the LMC. A spatial analysis shows that this dilemma is the result of an exceptional property: the grains of the LMC have on average a larger intrinsic submm opacity (emissivity index beta~1.7 and opacity kappa_abs(160 microns)=1.6 m2/kg) than those of the Galaxy. By studying the spatial distribution of the gas-to-dust mass ratio, we are able to constrain the fraction of unseen gas mass between ~10, and ~100 % and show that it is not sufficient to explain the gas-to-dust mass ratio obtained with Milky Way type grains. Finally, we confirm the detection of a 500 microns extended emission excess with an average relative amplitude of ~15 %, varying up to 40 %. This excess anticorrelates well with the dust mass surface density. Although we do not know the origin of this excess, we show that it is unlikely the result of very cold dust, or CMB fluctuations.


title: Mechanics: non-classical, non-quantum

abstract: A non-classical, non-quantum theory, or NCQ, is any fully consistent theory that differs fundamentally from both the corresponding classical and quantum theories, while exhibiting certain features common to both. Such theories are of interest for two primary reasons. Firstly, NCQs arise prominently in semi-classical approximation schemes. Their formal study may yield improved approximation techniques in the near-classical regime. More importantly for the purposes of this note, it may be possible for NCQs to reproduce quantum results over experimentally tested regimes while having a well defined classical limit, and hence are viable alternative theories. We illustrate an NCQ by considering an explicit class of NCQ mechanics. Here this class will be arrived at via a natural generalization of classical mechanics formulated in terms of a probability density functional.


title: Microwave sidebands for atomic physics experiments by period one oscillation in optically injected diode lasers

abstract: We show that nonlinear dynamics in diode lasers with optical injection leads to frequency tunable microwave sidebands which are suitable for atomic physics experiments. We demonstrate the applicability of the sidebands in an experiment where rubidium atoms are magneto-optically trapped with both the trap and the re-pump optical frequencies derived from one optically injected laser. We find linewidth narrowing in the optical spectrum of the injected laser at both the injection frequency and the sideband frequency. With strong optical injection which leads to frequency locking we find a complete linewidth transfer from the master to the slave. Further applications are discussed.


title: Two-Way Automata Making Choices Only at the Endmarkers

abstract: The question of the state-size cost for simulation of two-way nondeterministic automata (2NFAs) by two-way deterministic automata (2DFAs) was raised in 1978 and, despite many attempts, it is still open. Subsequently, the problem was attacked by restricting the power of 2DFAs (e.g., using a restricted input head movement) to the degree for which it was already possible to derive some exponential gaps between the weaker model and the standard 2NFAs. Here we use an opposite approach, increasing the power of 2DFAs to the degree for which it is still possible to obtain a subexponential conversion from the stronger model to the standard 2DFAs. In particular, it turns out that subexponential conversion is possible for two-way automata that make nondeterministic choices only when the input head scans one of the input tape endmarkers. However, there is no restriction on the input head movement. This implies that an exponential gap between 2NFAs and 2DFAs can be obtained only for unrestricted 2NFAs using capabilities beyond the proposed new model. As an additional bonus, conversion into a machine for the complement of the original language is polynomial in this model. The same holds for making such machines self-verifying, halting, or unambiguous. Finally, any superpolynomial lower bound for the simulation of such machines by standard 2DFAs would imply L<>NL. In the same way, the alternating version of these machines is related to L =? NL =? P, the classical computational complexity problems.


title: Lyndon-Shirshov basis and anti-commutative algebras

abstract: Chen, Fox, Lyndon 1958 \cite{CFL58} and Shirshov 1958 \cite{Sh58} introduced non-associative Lyndon-Shirshov words and proved that they form a linear basis of a free Lie algebra, independently. In this paper we give another approach to definition of Lyndon-Shirshov basis, i.e., we find an anti-commutative Gr\"{o}bner-Shirshov basis $S$ of a free Lie algebra such that $Irr(S)$ is the set of all non-associative Lyndon-Shirshov words, where $Irr(S)$ is the set of all monomials of $N(X)$, a basis of the free anti-commutative algebra on $X$, not containing maximal monomials of polynomials from $S$. Following from Shirshov's anti-commutative Gr\"{o}bner-Shirshov bases theory \cite{S62a2}, the set $Irr(S)$ is a linear basis of a free Lie algebra.


title: Bayesian multivariate mixed-scale density estimation

abstract: Although continuous density estimation has received abundant attention in the Bayesian nonparametrics literature, there is limited theory on multivariate mixed scale density estimation. In this note, we consider a general framework to jointly model continuous, count and categorical variables under a nonparametric prior, which is induced through rounding latent variables having an unknown density with respect to Lebesgue measure. For the proposed class of priors, we provide sufficient conditions for large support, strong consistency and rates of posterior contraction. These conditions allow one to convert sufficient conditions obtained in the setting of multivariate continuous density estimation to the mixed scale case. To illustrate the procedure a rounded multivariate nonparametric mixture of Gaussians is introduced and applied to a crime and communities dataset.


title: Dipping -versus Flaring in Z-track sources: resolving the controversy

abstract: We review the longterm confusion which has existed over the nature of flaring in the brightest class of low mass X-ray binary: the Z-track sources, specifically in the Cygnus X-2 sub-group. Intensity reductions in the lightcurve produce a branch in colour -colour diagrams similar to that of real flares in the Sco X-1 like group, and the nature of this branch was not clear. However, based on observations of Cygnus X-2 in which this dipping/flaring occurred it was proposed that the mass accretion rate in Z-track sources in general increases monotonically along the Z-track towards the Flaring Branch, a standard assumption widely held. It was also suggested that the Cygnus X-2 group have high inclination. Based on recent multi-wavelength observations of Cygnus X-2 we resolve these issues, showing by spectral analysis that the Dipping Branch consists of absorption events in the outer disk, unrelated to the occasional real flaring in the source. Thus motivation for Mdot increasing along the Z from Horizontal - Normal to Flaring Branch is removed, as is the idea that high inclination distinguishes the Cygnus X-2 group. Finally, the observations provide further evidence for the extended nature of the Accretion Disk Corona (ADC), and the correct modelling of the ADC Comptonized emission is crucial to the interpretation of low mass X-ray binary data.


title: Millihertz Quasi-periodic Optical Oscillations in 4U 0614+091

abstract: We report the discovery of a 1mHz optical quasi-periodic oscillation (QPO) in the candidate ultracompact low-mass X-ray binary 4U 0614+091. The ultra-low frequency QPO has no X-ray counterpart in contemporaneous RXTE/PCA data and is likely a signature of structure in the accretion disk. The QPO can be reasonably fitted with a single sine wave but with a phase jump part way through the observation, indicating that it is not coherent.We also identify a 48 min modulation, approximately consistent with the suggested orbital period of O'Brien (2005) and Shahbaz et al. (2008). If this is indeed orbital, it supports an identification of 4U 0614+091 as an ultra-compact source.


title: Note on rainbow connection number of dense graphs

abstract: An edge-colored graph $G$ is rainbow connected if any two vertices are connected by a path whose edges have distinct colors. The rainbow connection number of a connected graph $G$, denoted by $rc(G)$, is the smallest number of colors that are needed in order to make $G$ rainbow connected. Following an idea of Caro et al., in this paper we also investigate the rainbow connection number of dense graphs. We show that for $k\geq 2$, if $G$ is a non-complete graph of order $n$ with minimum degree $\delta (G)\geq \frac{n}{2}-1+log_{k}{n}$, or minimum degree-sum $\sigma_{2}(G)\geq n-2+2log_{k}{n}$, then $rc(G)\leq k$; if $G$ is a graph of order $n$ with diameter 2 and $\delta (G)\geq 2(1+log_{\frac{k^{2}}{3k-2}}{k})log_{k}{n}$, then $rc(G)\leq k$. We also show that if $G$ is a non-complete bipartite graph of order $n$ and any two vertices in the same vertex class have at least $2log_{\frac{k^{2}}{3k-2}}{k}log_{k}{n}$ common neighbors in the other class, then $rc(G)\leq k$.


title: About equality of the cosine transform to the sine transform of Fourier and the transform of Laplace

abstract: Regularity of the transform of Laplace in the opened area of 0 is proved with help of some methods of the transform of Fourier. The class of the transform of Laplace from the transform of Fourier is considered from some functions without a regularity in null. The functions are regular in the opened area of 0.


title: Quantum chaos and effective thermalization

abstract: We demonstrate effective equilibration for unitary quantum dynamics under conditions of classical chaos. Focusing on the paradigmatic example of the Dicke model, we show how a constructive description of the thermalization process is facilitated by the Glauber $Q$ or Husimi function, for which the evolution equation turns out to be of Fokker-Planck type. The equation describes a competition of classical drift and quantum diffusion in contractive and expansive directions. By this mechanism the system follows a 'quantum smoothened' approach to equilibrium, which avoids the notorious singularities inherent to classical chaotic flows.


title: Pulsed Gamma Rays from the Original Millisecond and Black Widow Pulsars: a case for Caustic Radio Emission?

abstract: We report the detection of pulsed gamma-ray emission from the fast millisecond pulsars (MSPs) B1937+21 (also known as J1939+2134) and B1957+20 (J1959+2048) using 18 months of survey data recorded by the \emph{Fermi} Large Area Telescope (LAT) and timing solutions based on radio observations conducted at the Westerbork and Nan\c{c}ay radio telescopes. In addition, we analyzed archival \emph{RXTE} and \emph{XMM-Newton} X-ray data for the two MSPs, confirming the X-ray emission properties of PSR B1937+21 and finding evidence ($\sim 4\sigma$) for pulsed emission from PSR B1957+20 for the first time. In both cases the gamma-ray emission profile is characterized by two peaks separated by half a rotation and are in close alignment with components observed in radio and X-rays. These two pulsars join PSRs J0034-0534 and J2214+3000 to form an emerging class of gamma-ray MSPs with phase-aligned peaks in different energy bands. The modeling of the radio and gamma-ray emission profiles suggests co-located emission regions in the outer magnetosphere.


title: A multiwavelength study of the supernova remnant G296.8-0.3

abstract: We report XMM-Newton observations of the Galactic supernova remnant G296.8-0.3, together with complementary radio and infrared data. The spatial and spectral properties of the X-ray emission, detected towards G296.8-0.3, was investigated in order to explore the possible evolutionary scenarios and the physical connexion with its unusual morphology detected at radio frequencies. G296.8-0.3 displays diffuse X-ray emission correlated with the peculiar radio morphology detected in the interior of the remnant and with the shell-like radio structure observed to the northwest side of the object. The X-ray emission peaks in the soft/medium energy range (0.5-3.0 keV). The X-ray spectral analysis confirms that the column density is high (NH \sim 0.64 x 10^{22} cm^{-2}) which supports a distant location (d>9 kpc) for the SNR. Its X-ray spectrum can be well represented by a thermal (PSHOCK) model, with kT \sim 0.86 keV, an ionization timescale of 6.1 x 10^{10} cm^{-3} s, and low abundance (0.12 Z_sun). The 24 microns observations show shell-like emission correlated with part of the northwest and southeast boundaries of the SNR. In addition a point-like X-ray source is also detected close to the geometrical center of the radio SNR. The object presents some characteristics of the so-called compact central objects (CCO). Its X-ray spectrum is consistent with those found at other CCOs and the value of NH is consistent with that of G296.8-0.3, which suggests a physical connexion with the SNR.


title: A Relative Value Iteration Algorithm for Non-degenerate Controlled Diffusions

abstract: The ergodic control problem for a non-degenerate controlled diffusion controlled through its drift is considered under a uniform stability condition that ensures the well-posedness of the associated Hamilton-Jacobi-Bellman (HJB) equation. A nonlinear parabolic evolution equation is then proposed as a continuous time continuous state space analog of White's `relative value iteration' algorithm for solving the ergodic dynamic programming equation for the finite state finite action case. Its convergence to the solution of the HJB equation is established using the theory of monotone dynamical systems and also, alternatively, by using the theory of reverse martingales.


title: Simulations of the Hyades

abstract: Context: Using the recent observational data of R\"oser et al. we present $N$-body simulations of the Hyades open cluster. Aims: We make an attempt to determine initial conditions of the Hyades cluster at the time of its formation in order to reproduce the present-day cumulative mass profile, stellar mass and luminosity function (LF). Methods: We performed direct $N$-body simulations of the Hyades in an analytic Milky Way potential that account for stellar evolution and include primordial binaries in a few models. Furthermore, we applied a Kroupa (2001) IMF and used extensive ensemble-averaging. Results: We find that evolved single-star King initial models with King parameters $W_0 = 6-9$ and initial particle numbers $N_0 = 3000$ provide good fits to the observational present-day cumulative mass profile within the Jacobi radius. The best-fit King model has an initial mass of $1721\ M_\odot$ and an average mass loss rate of $-2.2 \ M_\odot/\mathrm{Myr}$. The K-band LFs of models and observations show a reasonable agreement. Mass segregation is detected in both observations and models. If 33% primordial binaries are included the initial particle number is reduced by 5% as compared to the model without primordial binaries. Conclusions: The present-day properties of the Hyades can be well reproduced by a standard King or Plummer initial model when choosing appropriate initial conditions. The degeneracy of good-fitting models can be quite high due to the large dimension of the parameter space. More simulations with different Roche-lobe filling factors and primordial binary fractions are required to explore this degeneracy in more detail.


title: The near-IR $M_{bh}$ - L and $M_{bh}$ - n relations

abstract: We present near-IR surface photometry (2D-profiling) for a sample of 29 nearby galaxies for which super-massive black hole (SMBH) masses are constrained. The data is derived from the UKIDSS-LASS survey representing a significant improvement in image quality and depth over previous studies based on 2MASS data. We derive the spheroid luminosity and spheroid S\'ersic index for each galaxy with GALFIT3 and use these data to construct SMBH mass -bulge luminosity ($M_{\rm bh}$--$L$) and SMBH - S\'ersic index ($M_{\rm bh}$--$n$) relations. The best fit K-band relation for elliptical and disk galaxies is $\log(M_{\rm bh}/M_{\odot})= -0.36(\pm 0.03) (M_{\rm K} + 18) + 6.17(\pm 0.16)$ with an intrinsic scatter of 0.4$^{+0.09}_{-0.06}$dex whilst for elliptical galaxies we find $\log(M_{\rm bh}/M_{\odot})= -0.42(\pm 0.06) (M_{\rm K} + 22) + 7.5(\pm 0.15)$ with an intrinsic scatter of 0.31$^{+0.087}_{-0.047}$dex. Our revised $M_{\rm bh}$--$L$ relation agrees closely with the previous near-IR constraint by \citet{tex:G07}. The lack of improvement in the intrinsic scatter in moving to higher quality near-IR data suggests that the SMBH relations are not currently limited by the quality of the imaging data but is either intrinsic or a result of uncertainty in the precise number of required components required in the profiling process. Contrary to expectation (see \citealt{tex:GD07a}) a relation between SMBH mass and the S\'ersic index was not found at near-IR wavelengths. This latter outcome is believed to be explained by the generic inconsistencies between 1D and 2D galaxy profiling which are currently under further investigation.


title: Density of states of continuous and discrete spin models: a case study

abstract: A relation between O(n) lattice spin models and Ising models defined on the same lattice was recently put forward [L. Casetti, C. Nardini, and R. Nerattini, Phys. Rev. Lett. 106, 057208 (2011)]. Such a relation, inspired by an energy landscape analysis, implies that the density of states of an O(n) spin model on a lattice can be effectively approximated, at least close to the phase transition, in terms of the density of states of an Ising model defined on the same lattice and with the same interactions. In the present paper we show that such a relation exactly holds, albeit in a slightly modified form, in the special cases of the mean-field XY model and of the one-dimensional XY model. We also discuss the possible consequences of this result for the general case.


title: Performance improvement of the software development project using the Value Management approach

abstract: Improving performance and delivering value for customers have become a central theme in business. The software industry has become an increasingly important sector for the economy growth in Tunisia. This study aims to show how using Value Management in the Tunisian software industry for project analysis gives new insight about true project value and performance. This new approach is considered as an appropriate tool for guiding the process of making decisions. It offers tools in order to analyze the service value from the customer and organization perspectives. The results showed that the VM allows to have better performance in the software development project by linking customer satisfaction and cost analysis. The present case shows to service managers how they can benchmark project function to reduce their costs and improve resource allocation taking into consideration what customers consider important during their overall service experience. It can identify best professional practices, orient decisions to improve service value


title: Ferromagnetism in YbCu2Si2 at high pressure

abstract: We demonstrate from detailed ac-susceptibility and calorimetry studies under hydrostatic pressure that YbCu2Si2 probably orders ferromagnetically at high pressure. The (p,H,T) phase diagram, shows that the transition temperature increases with pressure, but also with an applied magnetic field. We suggest that many ytterbium systems may show a trend towards ferromagnetism and we discuss the possible reasons for this. We also examine the implications, including the potential of YbCu2Si2 and other Yb compounds for further studies of the rich physical properties that may occur near a ferromagnetic critical point.


title: Time-Reversal Invariance Violation in Heavy and in Few-body Nuclei

abstract: Time reversal invariance violating (TRIV) effects in neutron scattering are very important in a search for new physics, being complementary to neutron and atomic electric dipole moment measurements. In this relation, a sensitivity of TRIV observables to different models of CP-violation and their dependencies on nuclear structure, which can lead to new enhancement factors, are discussed.


title: Higher-order adaptive finite-element methods for orbital-free density functional theory

abstract: In the present work, we investigate the computational efficiency afforded by higher-order finite-element discretization of the saddle-point formulation of orbital-free density functional theory. We first investigate the robustness of viable solution schemes by analyzing the solvability conditions of the discrete problem. We find that a staggered solution procedure where the potential fields are computed consistently for every trial electron-density is a robust solution procedure for higher-order finite-element discretizations. We next study the numerical convergence rates for various orders of finite-element approximations on benchmark problems. We obtain close to optimal convergence rates in our studies, although orbital-free density-functional theory is nonlinear in nature and some benchmark problems have Coulomb singular potential fields. We finally investigate the computational efficiency of various higher-order finite-element discretizations by measuring the CPU time for the solution of discrete equations on benchmark problems that include large Aluminum clusters. In these studies, we use mesh coarse-graining rates that are derived from error estimates and an a priori knowledge of the asymptotic solution of the far-field electronic fields. Our studies reveal a significant 100-1000 fold computational savings afforded by the use of higher-order finite-element discretization, alongside providing the desired chemical accuracy. We consider this study as a step towards developing a robust and computationally efficient discretization of electronic structure calculations using the finite-element basis.


title: A review of the Z-track sources

abstract: The brightest class of low mass X-ray binary source: the Z-track sources are reviewed specifically with regard to the nature of the three distinct states of the sources. A physical model is presented for the Cygnus X-2 sub-group in which increasing mass accretion rate takes place on the Normal Branch resulting in high neutron star temperature and radiation pressure responsible for inner disk disruption and launching of jets. The Flaring Branch consists of unstable nuclear burning on the neutron star. It is shown that the Sco X-1 like sub-group is dominated by almost non-stop flaring consisting of both unstable burning and increase of Mdot, causing higher neutron star temperatures. Finally, results of Atoll source surveys are presented and a model for the nature of the Banana and Island states in these sources is proposed. Motion along the Banana state is caused by variation of Mdot. Measurements of the high energy cut-off of the Comptonized emission E_CO provide the electron temperature T_e of the Comptonizing ADC; above a luminosity of 2x10^37 erg/s E_CO is a few keV and T_e equals the neutron star temperature. At lower luminosities, the cut-off energy rises towards 100 keV showing heating of the corona by an unknown process. This spectral hardening is the cause of the Island state of Atoll sources. The models for Z-track and Atoll sources thus constitute a unified model for low mass X-ray binary sources.


title: Dependence of a planet's chaotic zone on particle eccentricity: the shape of debris disc inner edges

abstract: The orbit of a planet is surrounded by a chaotic zone wherein nearby particles' orbits are chaotic and unstable. Wisdom (1980) showed that the chaos is driven by the overlap of mean motion resonances which occurs within a distance (da/a)_chaos = 1.3 mu^2/7 of the planet's orbit. However, the width of mean motion resonances grows with the particles' eccentricity, which will increase the width of the chaotic zone at higher eccentricities. Here we investigate the width of the chaotic zone using the iterated encounter map and N-body integrations. We find that the classical prescription works well for particles on low-eccentricity orbits. However, above a critical eccentricity, dependent upon the mass of the planet, the width of the chaotic zone increases with eccentricity. An extension of Wisdom's analytical arguments then shows that, above the critical eccentricity, the chaotic zone width is given by (da/a)_chaos = 1.8 e^1/5 mu^1/5, which agrees well with the encounter map results. The critical eccentricity is given by e_cr = 0.21 mu^3/7. This extended chaotic zone results in a larger cleared region when a planet sculpts the inner edge of a debris disc comprised of eccentric planetesimals. Hence, the planet mass estimated from the classical chaotic zone may be erroneous. We apply this result to the HR8799 system, showing that the masses of HR8799b inferred from the truncation of the disc may vary by up to 50% depending on the disc particles' eccentricities. With a disc edge at 90AU, the necessary mass of planet b to cause the truncation would be 8-10 Jovian masses if the disc particles have low eccentricities (<0.02), but only 4-8 Jovian masses if the disc particles have higher eccentricities. Our result also has implications for the ability of a planet to feed material into an inner system, a process which may explain metal pollution in White Dwarf atmospheres.


title: New mathematical model for fluid-glucose-albumin transport in peritoneal dialysis

abstract: A mathematical model for fluid transport in peritoneal dialysis is constructed. The model is based on a three-component nonlinear system of two-dimensional partial differential equations for fluid, glucose and albumin transport with the relevant boundary and initial conditions. Non-constant steady-state solutions of the model are studied. The restrictions on the parameters arising in the model are established with the aim to obtain exact formulae for the non-constant steady-state solutions. As the result, the exact formulae for the fluid fluxes from blood to tissue and across the tissue were constructed together with two linear autonomous ODEs for glucose and albumin concentrations. The analytical results were checked for their applicability for the description of fluid-glucose-albumin transport during peritoneal dialysis.


title: On the Rate of Convergence to the Marchenko--Pastur Distribution

abstract: Let $\mathbf X=(X_{jk})$ denote $n\times p$ random matrix with entries $X_{jk}$, which are independent for $1\le j\le n,1\le k\le p$. We consider the rate of convergence of empirical spectral distribution function of the matrix $\mathbf W=\frac1p\mathbf X\mathbf X^*$ to the Marchenko--Pastur law. We assume that $\mathbf E X_{jk}=0$, $\mathbf E X_{jk}^2=1$ and that the distributions of the matrix elements $X_{jk}$ have a uniformly sub exponential decay in the sense that there exists a constant $\varkappa>0$ such that for any $1\le j \le n,\,1\le k\le p $ and any $t\ge 1$ we have $$ \mathbf{ Pr}\{|X_{jk}|>t\}\le \varkappa^{-1}\exp\{-t^{\varkappa}\}. $$ By means of a recursion argument it is shown that the Kolmogorov distance between the empirical spectral distribution of the sample covariance matrix $\mathbf W$ and the Marchenko--Pastur distribution is of order $O(n^{-1}\log^{4+\frac4{\varkappa}} n)$ with high probability.


title: Dynamics of fibers in a wide microchannel

abstract: Dynamics of single flexible non-Brownian fibers, tumbling in a Poiseuille flow between two parallel solid plane walls, is studied with the use of the hydromultipole numerical code, based on the multipole expansion of the Stokes equations, corrected for lubrication. It is shown that for a wide range of the system parameters, the migration rate towards the middle plane of the channel increases for fibers, which are closer to a wall, or are more flexible (less stiff), or are longer. The faster motion towards the channel center is accompanied by a slower translation along the flow and a larger fiber deformation.


title: Structural Optimisation: Biomechanics of the Femur

abstract: A preliminary iterative 3D meso-scale structural model of the femur was developed, in which bar and shell elements were used to represent trabecular and cortical bone respectively. The cross-sectional areas of the bar elements and the thickness values of the shell elements were adjusted over successive iterations of the model based on a target strain stimulus, resulting in an optimised construct. The predicted trabecular architecture, and cortical thickness distribution showed good agreement with clinical observations, based on the application of a single leg stance load case during gait. The benefit of using a meso-scale structural approach in comparison to micro or macro-scale continuum approaches to predictive bone modelling was achievement of the symbiotic goals of computational efficiency and structural description of the femur.


title: Using rates to measure mixed modulus-anomaly mediated supersymmetry breaking at the LHC

abstract: If SUSY is discovered at the LHC, the task will immediately turn to determining the model of SUSY breaking. Here, we employ a Mixed Modulus-Anomaly Mediated SUSY Breaking (MMAMSB) model with very similar LHC phenomenology to the more conventionally studied Constrained Minimal SUSY Model (CMSSM) and minimal Anomaly Mediated SUSY Breaking (mAMSB) models. We then study whether the models can be distinguished and measured. If we only fit to the various mass edges and mass end-points from cascade decay chains that are normally studied, a unique determination and measurement of the model is problematic without substantial amounts of LHC data. However, if event rate information is included, we can quickly distinguish and measure the correct SUSY model and exclude alternatives.


title: Probing Substellar Companions of AGB Stars through Spirals and Arcs

abstract: Recent observations of strikingly well-defined spirals in the circumstellar envelopes of asymptotic giant branch (AGB) stars point to the existence of binary companions in these objects. In the case of planet or brown dwarf mass companions, we investigate the observational properties of the spiral-onion shell wakes due to the gravitational interaction of these companions with the outflowing circumstellar matter. Three dimensional hydrodynamical simulations at high resolution show that the substellar mass objects produce detectable signatures at 100 AU distance, for the wake induced by a Jupiter to brown dwarf mass object orbiting a solar mass AGB star. In particular, the arm pattern propagates with a speed depending on the local wind and sound speeds, implying possible variations in the arm separation in the wind acceleration region and/or in a slow wind with significant temperature variation. The pattern propagation speeds of the inner and outer boundaries differ by twice the sound speed, leading to the overlap of high density boundaries in slow winds and producing a subpattern of the spiral arm feature. Vertically, the wake forms concentric arcs with angular sizes anticorrelated to the wind Mach number. We provide an empirical formula for the peak density enhancement as a function of the mass, orbital distance, and velocity of the object as well as the wind and local sound speeds. In typical condition of AGB envelopes, the arm-interarm density contrast can be greater than 30 % of the background density within a distance of ~10(M_p/M_J) AU for the object mass M_p in units of Jupiter mass M_J. These results suggest that such features may probe unseen substellar mass objects embedded in the winds of AGB stars and may be useful in planning future high sensitivity/resolution observations with ALMA.


title: Highly efficient broadband conversion of light polarization by composite retarders

abstract: Driving on an analogy with the technique of composite pulses in quantum physics, we propose highly efficient broadband polarization converters composed of sequences of ordinary retarders rotated at specific angles with respect to their fast-polarization axes.


title: Filtrations on instanton homology

abstract: In earlier work of the authors, the Khovanov complex of a knot or link appeared as the first page in a spectral sequence abutting to the instanton homology. The quantum and (co)homological gradings on Khovanov homology do not survive as gradings, but we show that they survive as filtrations.