title: kmos: A lattice kinetic Monte Carlo framework

abstract: Kinetic Monte Carlo (kMC) simulations have emerged as a key tool for microkinetic modeling in heterogeneous catalysis and other materials applications. Systems, where site-specificity of all elementary reactions allows a mapping onto a lattice of discrete active sites, can be addressed within the particularly efficient lattice kMC approach. To this end we describe the versatile kmos software package, which offers a most user-friendly implementation, execution, and evaluation of lattice kMC models of arbitrary complexity in one- to three-dimensional lattice systems, involving multiple active sites in periodic or aperiodic arrangements, as well as site-resolved pairwise and higher-order lateral interactions. Conceptually, kmos achieves a maximum runtime performance which is essentially independent of lattice size by generating code for the efficiency-determining local update of available events that is optimized for a defined kMC model. For this model definition and the control of all runtime and evaluation aspects kmos offers a high-level application programming interface. Usage proceeds interactively, via scripts, or a graphical user interface, which visualizes the model geometry, the lattice occupations and rates of selected elementary reactions, while allowing on-the-fly changes of simulation parameters. We demonstrate the performance and scaling of kmos with the application to kMC models for surface catalytic processes, where for given operation conditions (temperature and partial pressures of all reactants) central simulation outcomes are catalytic activity and selectivities, surface composition, and mechanistic insight into the occurrence of individual elementary processes in the reaction network.


title: Classical novae and type I X-ray bursts: challenges for the 21st century

abstract: Classical nova explosions and type I X-ray bursts are the most frequent types of thermonuclear stellar explosions in the Galaxy. Both phenomena arise from thermonuclear ignition in the envelopes of accreting compact objects in close binary star systems. Detailed observations of these events have stimulated numerous studies in theoretical astrophysics and experimental nuclear physics. We discuss observational features of these phenomena and theoretical efforts to better understand the energy production and nucleosynthesis in these explosions. We also examine and summarize studies directed at identifying nuclear physics quantities with uncertainties that significantly affect model predictions.


title: Dynamical control of interference using voltage pulses in the quantum regime

abstract: As a general trend, nanoelectronics experiments are shifting toward frequencies so high that they become comparable to the device's internal characteristic time scales, resulting in new opportunities for studying the dynamical aspects of quantum mechanics. Here we theoretically study how a voltage pulse (in the quantum regime) propagates through an electronic interferometer (Fabry- Perot or Mach-Zehnder). We show that extremely fast pulses provide a conceptually new tool for manipulating quantum information: the possibility to dynamically engineer the interference pattern of a quantum system. Striking physical signatures are associated with this new regime: restoration of the interference in presence of large bias voltages; negative currents with respect to the direction of propagation of the voltage pulse; and oscillation of the total transmitted charge with the total number of injected electrons. The present findings have been made possible by the recent unlocking of our capability for simulating time-resolved quantum nanoelectronics of large systems.


title: Optimal feedback control, linear first-order PDE systems, and obstacle problems

abstract: We introduce an alternative approach for the analysis and numerical approximation of the optimal feedback control mapping. It consists in looking at a typical optimal control problem in such a way that feasible controls are mappings depending both in time and space. In this way, the feedback form of the problem is built-in from the very beginning. Optimality conditions are derived for one such optimal mapping, which by construction is the optimal feedback mapping of the problem. In formulating optimality conditions, costates in feedback form are solutions of linear, first-order transport systems, while optimal descent directions are solutions of appropriate obstacle problems. We treat situations with no constraint-sets for control and state, as well as the more general case where a constraint-set is considered for the control variable. Constraints for the state variable are deferred to a coming contribution.


title: Generic measures for geodesic flows on nonpositively curved manifolds

abstract: We study the generic invariant probability measures for the geodesic flow on connected complete nonpositively curved manifolds. Under a mild technical assumption, we prove that ergodicity is a generic property in the set of probability measures defined on the unit tangent bundle of the manifold and supported by trajectories not bounding a flat strip. This is done by showing that Dirac measures on periodic orbits are dense in that set. In the case of a compact surface, we get the following sharp result: ergod- icity is a generic property in the space of all invariant measures defined on the unit tangent bundle of the surface if and only if there are no flat strips in the universal cover of the surface. Finally, we show under suitable assumptions that generically, the invariant probability measures have zero entropy and are not strongly mixing.


title: In-medium jet shape from energy collimation in parton showers: Comparison with CMS PbPb data at 2.76 TeV

abstract: We present the medium-modified energy collimation in the leading-logarithmic approximation (LLA) and next-to-leading-logarithmic approximation (NLLA) of QCD. As a consequence of more accurate kinematic considerations in the argument of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) fragmentation functions (FFs) we find a new NLLA correction ${\cal O}(\alpha_s)$ which accounts for the scaling violation of DGLAP FFs at small $x$. The jet shape is derived from the energy collimation within the same approximations and we also compare our calculations for the energy collimation with the event generators Pythia 6 and YaJEM for the first time in this paper. The modification of jets by the medium in both cases is implemented by altering the infrared sector using the Borghini-Wiedemann model. The energy collimation and jet shapes qualitatively describe a clear broadening of showers in the medium, which is further supported by YaJEM in the final comparison of the jet shape with CMS PbPb data at center-of-mass energy 2.76 TeV. The comparison of the biased versus unbiased YaJEM jet shape with the CMS data shows a more accurate agreement for biased showers and illustrates the importance of an accurate simulation of the experimental jet-finding strategy.


title: Features of ultrasound absorption by dislocations in subgrain-free cd\boldmath$_{0.2} $hg$_{0.8}$te crystals

abstract: The temperature dependence of the ultrasound wave absorption in bulk $p$-Cd$_{0.2}$Hg$_{0.8}$Te crystals free from low-angle grain boundaries has been studied experimentally for the first time in the frequency range 10--55 MHz and the temperature interval 150--300 K, and the corresponding results of measurements are presented. The maximum value of absorption coefficient is found to increase and to shift toward higher temperatures, as the ultrasound frequency grows. The results obtained can be satisfactorily explained in the framework of the Brailsford model, which associates the ultrasound absorption with vibrations of thermally activated dislocation kinks. The characteristic parameters of this model for $p$-Cd$_{0.2}$Hg$_{0.8}$Te are determined; namely, the frequency coefficient $f_{k}\approx $~6\,$\times $\,10$^{\,9}$\,Hz and the kink diffusion activation energy $W_{k}\approx $~0.11\,eV. The dislocation concentration is also evaluated ($\alpha \approx $~2\,$\times$\,10$^{10}$\,m$^{-2}$), with the determined value being consistent with that obtained by the selective etching method (0.7\,$\times $\,10$^{10}$\,m$^{-2}$).


title: Comparaison between the two models : new approach using the $\alpha$-divergence

abstract: We propose new nonparametric accordance R\'enyi-$\alpha$ and $\alpha$-Tsallis divergence estimators for continuous distributions. We discuss this approach with a view to the selection model (on al\'etoire and autoregressive AR (1)). We lestimateur used by kernel density esttimer underlying. Nevertheless, we are able to prove that the estimators are consistent under certain conditions. We also describe how to apply these estimators and demonstrate their effectiveness through numerical experiments.


title: If the spin of gravitational field depends on spacetime dimension

abstract: In this manuscript a different perspective on gravitational field has been proposed, suggesting different features of gravity depending on spacetime dimension, which can explain preventing the formation of singularity inside blackholes and also suggest interpreting spin-1 fields emerged in theories with dimensional reduction mechanism, as different aspects of gravitational field.


title: k-Colorability is Graph Automaton Recognizable

abstract: Automata operating on general graphs have been introduced by virtue of graphoids. In this paper we construct a graph automaton that recognizes $k$-colorable graphs.


title: Neutron Production Rates by Inverse-Beta Decay in Fully Ionized Plasmas

abstract: Recently we showed that the nuclear transmutation rates are largely overestimated in the Widom-Larsen theory of the so called `Low Energy Nuclear Reactions'. Here we show that unbound plasma electrons are even less likely to initiate nuclear transmutations.


title: Graphical Interface for Visually Impaired People Based on Bi-stable Solenoids

abstract: In this paper a concept for hardware realization of graphic tactile display for visually impaired peoples is presented. For realization of tactile actuators bi-stable, solenoids and PIC based control board are used. The selected algorithm for series activation of each row of display allows using minimal number of active components to set and reset the solenoids. Finally, a program algorithm of control board is discussed. The project is funded by Bulgarian National Science Fund NSF Grant D ID 02 14, 2009 2013


title: Topology of Asymptotic Cones and S-machine

abstract: Sapir, Birget and Rips showed how to construct groups from Turing machines. To achieve such a construction they introduced the notion of S-machine. Then considering a simplified S-machine Sapir and Olshanskii showed how to construct a group such that each of its asymptotic cone is non-simply connected. Still using the notion S-machine, they constructed a group with two asymptotic cone non-homeomorphic. In this paper we show that each asymptotic cone of a group constructed following the whole method of Sapir, Birget and Rips is not simply connected.


title: On quaternionic Monge-Ampere operator, closed positive currents and Lelong-Jensen type formula on quaternionic space

abstract: In this paper, we introduce the first-order differential operators $d_0$ and $d_1$ acting on the quaternionic version of differential forms on the flat quaternionic space $\mathbb{H}^n$. The behavior of $d_0,d_1$ and $\triangle=d_0d_1$ is very similar to $\partial,\overline{\partial}$ and $\partial \overline{\partial}$ in several complex variables. The quaternionic Monge-Amp\`{e}re operator can be defined as $(\triangle u)^n$ and has a simple explicit expression. We define the notion of closed positive currents in the quaternionic case, and extend several results in complex pluripotential theory to the quaternionic case: define the Lelong number for closed positive currents, obtain the quaternionic version of Lelong-Jensen type formula, and generalize Bedford-Taylor theory, i.e., extend the definition of the quaternionic Monge-Amp\`{e}re operator to locally bounded quaternionic plurisubharmonic functions and prove the corresponding convergence theorem.


title: Non-Termination Analysis of Java Bytecode

abstract: We introduce a fully automated static analysis that takes a sequential Java bytecode program P as input and attempts to prove that there exists an infinite execution of P. The technique consists in compiling P into a constraint logic program P_CLP and in proving non-termination of P_CLP; when P consists of instructions that are exactly compiled into constraints, the non-termination of P_CLP entails that of P. Our approach can handle method calls; to the best of our knowledge, it is the first static approach for Java bytecode able to prove the existence of infinite recursions. We have implemented our technique inside the Julia analyser. We have compared the results of Julia on a set of 113 programs with those provided by AProVE and Invel, the only freely usable non-termination analysers comparable to ours that we are aware of. Only Julia could detect non-termination due to infinite recursion.


title: 3D global simulations of a cosmic-ray-driven dynamo in dwarf galaxies

abstract: Star-forming dwarf galaxies can be seen as the local proxies of the high-redshift building blocks of more massive galaxies according to the current paradigm of the hierarchical galaxy formation. They are low-mass objects, and therefore their rotation speed is very low. Several galaxies are observed to show quite strong magnetic fields. These cases of strong ordered magnetic fields seem to correlate with a high, but not extremely high, star formation rate. We investigate whether these magnetic fields could be generated by the cosmic-ray-driven dynamo. The environment of a dwarf galaxy is unfavourable for the large-scale dynamo action because of the very slow rotation that is required to create the regular component of the magnetic field. We built a 3D global model of a dwarf galaxy that consists of two gravitational components: the stars and the dark-matter halo described by the purely phenomenological profile proposed previously. We solved a system of magnetohydrodynamic (MHD) equations that include an additional cosmic-ray component described by the fluid approximation. We found that the cosmic-ray-driven dynamo can amplify the magnetic field with an exponential growth rate. The $e$-folding time is correlated with the initial rotation speed. The final mean value of the azimuthal flux for our models is of the order of few $\mu$G and the system reaches its equipartition level. The results indicate that the cosmic-ray-driven dynamo is a process that can explain the magnetic fields in dwarf galaxies.


title: Well-posedness via Monotonicity. An Overview

abstract: The idea of monotonicity (or positive-definiteness in the linear case) is shown to be the central theme of the solution theories associated with problems of mathematical physics. A "grand unified" setting is surveyed covering a comprehensive class of such problems. We elaborate the applicability of our scheme with a number examples. A brief discussion of stability and homogenization issues is also provided.


title: First order Quantum Hall Transitions in Hofstadter Butterfly in the Honeycomb Lattice

abstract: We analyze the effects of nearest neighbor repulsive interactions in the Hofstadter system in a honeycomb lattice. At low fillings, we show that, as the interaction strength is increased there are two first order transitions, a Landau transition with translational and rotational symmetries broken, followed by a topological transition with a jump in the quantized Hall conductivity. We therefore predict that in physical realizations where the interaction effects are strong, there would be translation symmetry broken states with quantized Hall conductivities that differ from those predicted by the non-interacting theory.


title: How to suppress the backscattering of conduction electrons?

abstract: It is shown theoretically that the strong coupling of electrons to a high-frequency electromagnetic field results in the nulling of electron backscattering within the Born approximation. The conditions of the effect depend only on field parameters and do not depend on concrete form of scattering potential. As a consequence, this phenomenon is of universal physical nature and can take place in various conducting systems. Since the suppression of electron backscattering results in decreasing electrical resistance, the solved quantum-mechanical problem opens a new way to control electronic transport properties of conductors by a laser-generated field. Particularly, the elaborated theory is applicable to nanostructures exposed to a strong monochromatic electromagnetic wave.


title: Navigating MazeMap: indoor human mobility, spatio-logical ties and future potential

abstract: Global navigation systems and location-based services have found their way into our daily lives. Recently, indoor positioning techniques have also been proposed, and there are several live or trial systems already operating. In this paper, we present insights from MazeMap, the first live indoor/outdoor positioning and navigation system deployed at a large university campus in Norway. Our main contribution is a measurement case study; we show the spatial and temporal distribution of MazeMap geo-location and wayfinding requests, construct the aggregated human mobility map of the campus and find strong logical ties between different locations. On one hand, our findings are specific to the venue; on the other hand, the nature of available data and insights coupled with our discussion on potential usage scenarios for indoor positioning and location-based services predict a successful future for these systems and applications.


title: Generalized Scaling Operators in White Noise Analysis and Applications to Hamiltonian Path Integrals with Quadratic Action

abstract: We give an outlook, how to realize the ideas of complex scaling from Grothaus, Vogel and Streit to phase space path integrals in the framework of White Noise Analysis. The idea of this scaling method goes back to Doss. Therefore we extend the concept complex scaling to scaling with suitable bounded operators.


title: Quantum state transfer in atom-cavity systems with uncolored Cayley interacting networks

abstract: Considering the two-photon exchange interaction between n coupled cavities each of them containing a two level atom, the atomic and photonic state transfer is investigated. In fact, n atom-cavity systems are considered to be distributed on the nodes of an uncolored Cayley network and interact with each other via the adjacency matrix of the corresponding network. Then, by employing the photon-excitation conservation and also the algebraic structure of the networks, such as irreducible characters of the groups associated with the networks, some suitable basis for the atom-cavity state space is introduced based on the corresponding generalized Fourier transform, so that the Hamiltonian of the whole system, is block diagonalized with two-dimensional blocks. Then, by solving the corresponding Schrodinger equation exactly, quantum state transfer and also entanglement generation between the atoms or the photons are discussed. For instance, the probability amplitudes associated with the photon transition between the cavities or excitation transition between the atoms are obtained in terms of the irreducibles characters of the corresponding network and the hopping parameter ? between the cavities.


title: How are identifiers named in open source software? About popularity and consistency

abstract: With the rapid increasing of software project size and maintenance cost, adherence to coding standards especially by managing identifier naming, is attracting a pressing concern from both computer science educators and software managers. Software developers mainly use identifier names to represent the knowledge recorded in source code. However, the popularity and adoption consistency of identifier naming conventions have not been revealed yet in this field. Taking forty-eight popular open source projects written in three top-ranking programming languages Java, C and C++ as examples, an identifier extraction tool based on regular expression matching is developed. In the subsequent investigation, some interesting findings are obtained. For the identifier naming popularity, it is found that Camel and Pascal naming conventions are leading the road while Hungarian notation is vanishing. For the identifier naming consistency, we have found that the projects written in Java have a much better performance than those written in C and C++. Finally, academia and software industry are urged to adopt the most popular naming conventions consistently in their practices so as to lead the identifier naming to a standard, unified and high-quality road.


title: Gravity-mediated (or Composite) Dark Matter Confronts Astrophysical Data

abstract: We consider the astrophysical bounds on a new form of dark matter, the so called Gravity-mediated Dark Matter. In this scenario, dark matter communicates with us through a mediator sector composed of gravitational resonances, namely a new scalar (radion) and a massive spin-two resonance (massive graviton). We consider specific models motivated by natural electroweak symmetry breaking or weak-scale dark matter in the context of models in warped extra-dimensions and their composite duals. The main Dark Matter annihilation mechanism is due to the interactions of KK gravitons to gauge bosons that propagate in bulk. We impose the bounds on monochromatic or continuum photons from Fermi-LAT and HESS. We also explore scenarios in which the Fermi gamma-ray line could be a manifestation of Gravity-mediated Dark Matter.


title: Quivers with loops and perverse sheaves

abstract: In this article we consider the general definition of Lustig sheaves for arbitrary quivers, possibly carrying loops. We answer a conjecture raised by Lusztig asking if the more "simple" Lusztig perverse sheaves are enough to span the whole Grothendieck group usually built with these sheaves. Our proof is based on induction, with the help of restriction and induction functors, but with non trivial first steps, consisting in the study of quivers with one vertex but possible loops. We also need to consider regularity conditions on the support of our perverse sheaves to perform efficient restrictions at imaginary vertices. From our proof emerges a new combinatorial structure on our generalized canonical basis, which is more general than the usual crystals, in that there are now more operators associated to a vertex with loops. In a second part, we construct and study a Hopf algebra which generalizes the usual quantum groups. The geometric study previously made leads to a natural definition, which includes countably infinite sets of generators at imaginary roots, with higher order Serre relations and commutativity conditions imposed by the Jordan quiver case. We finally prove that the positive part of this algebra is isomorphic to our Grothendieck group, thanks to the study of a nondegenerate Hopf pairing.


title: Slow, bursty dynamics as the consequence of quenched network topologies

abstract: Bursty dynamics of agents is shown to appear at criticality or in extended Griffiths phases, even in case of Poisson processes. I provide numerical evidence for power-law type of inter-communication time distributions by simulating the Contact Process and the Susceptible-Infected-Susceptible model. This observation suggests that in case of non-stationary bursty systems the observed non-poissonian behavior can emerge as the consequence of an underlying hidden poissonian network process, which is either critical or exhibits strong rare-region effects. On contrary, in time varying networks rare-region effects do not cause deviation from the mean-field behavior and heterogeneity induced burstyness is absent.


title: Generalized almost para-contact manifolds

abstract: In this paper, we introduce generalized almost para-contact manifolds and obtain normality conditions in terms of classical tensor fields. We show that such manifolds naturally carry certain Lie bialgebroid/quasi-Lie algebroid structures on them and we relate this new generalized manifolds with classical almost para-contact manifolds. The paper contains several examples.


title: Bounds on the ML Decoding Error Probability of RS-Coded Modulation over AWGN Channels

abstract: This paper is concerned with bounds on the maximum-likelihood (ML) decoding error probability of Reed-Solomon (RS) codes over additive white Gaussian noise (AWGN) channels. To resolve the difficulty caused by the dependence of the Euclidean distance spectrum on the way of signal mapping, we propose to use random mapping, resulting in an ensemble of RS-coded modulation (RS-CM) systems. For this ensemble of RS-CM systems, analytic bounds are derived, which can be evaluated from the known (symbol-level) Hamming distance spectrum. Also presented in this paper are simulation-based bounds, which are applicable to any specific RS-CM system and can be evaluated by the aid of a list decoding (in the Euclidean space) algorithm. The simulation-based bounds do not need distance spectrum and are numerically tight for short RS codes in the regime where the word error rate (WER) is not too low. Numerical comparison results are relevant in at least three aspects. First, in the short code length regime, RS-CM using BPSK modulation with random mapping has a better performance than binary random linear codes. Second, RS-CM with random mapping (time varying) can have a better performance than with specific mapping. Third, numerical results show that the recently proposed Chase-type decoding algorithm is essentially the ML decoding algorithm for short RS codes.


title: A Mobile Device Prototype Application for the Detection and Prediction of Node Faults in Wireless Sensor Networks

abstract: Various implementations of wireless sensor networks (i.e. personal area-, wireless body area- networks) are prone to node and network failures by such characteristics as limited node energy resources and hardware damage incurred from their surrounding environment (i.e. flooding, forest fires, a patient falling). This may jeopardize their reliability to act as early warning systems, monitoring systems for patients and athletes, and industrial and environmental observation networks. Following the current trend and widespread use of hand held, mobile communication devices, we outline an application architecture designed to detect and predict faulty nodes in wireless sensor networks. Furthermore, we implement our design as a proof of concept prototype for Android-based smartphones, which may be extended to develop other applications used for monitoring networked wireless personal area and body sensors used in other capacities. We have conducted several preliminary experiments to demonstrate the use of our design, which is capable of monitoring networks of wireless sensor devices and predicting node faults based on several localized metrics. As attributes of such networks may change over time, any models generated when the application is initialized must be updated periodically such that the applied machine learning algorithm maintains high levels of both accuracy and precision. The application is designed to discover node faults and, once identified, alert the user so that appropriate action may be taken.


title: Nuclear forces and ab initio calculations of atomic nuclei

abstract: Nuclear forces and the nuclear many-body problem have been some of Gerry Brown's main topics in his so productive life as a theoretical physicist. In this talk, I outline how Gerry's work laid the foundations of the modern theory of nuclear forces and ab initio calculations of atomic nuclei. I also present some recent developments obtained in the framework of nuclear lattice simulations.


title: Ultrafast Quenching of the Exchange Interaction in a Mott Insulator

abstract: We investigate how fast and how effective photocarrier excitation can modify the exchange interaction $J_\mathrm{ex}$ in the prototype Mott-Hubbard insulator. We demonstrate an ultrafast quenching of $J_\mathrm{ex}$ both by evaluating exchange integrals from a time-dependent response formalism and by explicitly simulating laser-induced spin precession in an antiferromagnet that is canted by an external magnetic field. In both cases, the electron dynamics is obtained from nonequilibrium dynamical mean-field theory. We find that the modified $J_\mathrm{ex}$ emerges already within a few electron hopping times after the pulse, with a reduction that is comparable to the effect of chemical doping.


title: The effect of the Basset history force on particle clustering in Homogeneous and Isotropic Turbulence

abstract: We study the effect of the Basset history force on the dynamics of small particles transported in homogeneous and isotropic turbulence and show that this term, often neglected in previous numerical studies, reduces the small-scale clustering typical of inertial particles. The contribution of this force to the total particle acceleration is, on average, responsible for about 10% of the total acceleration and particularly relevant during rare strong events. At moderate density ratios, i.e. sand or metal powder in water, its presence alters the balance of forces determining the particle acceleration.


title: Ultracompact Adiabatic Bi-sectional Tapered Coupler for the Si/III-V Heterogeneous Integration

abstract: An ultracompact bi-sectional adiabatic tapered coupler, which is suitable for an adiabatic mode transformation between a common single mode SOI wire waveguide and a Si/III-V hybrid waveguide, are proposed for Si/III-V heterogeneous integration. Since the bi-sectional tapered coupler mimics a semi-3D taper and avoids exciting the unwanted high-order modes in the thick p-InP cladding layer (which is removed in the first tapered section), the length of the adiabatic mode coupler can be dramatically shortened. Taking into account the mask design and the fabrication tolerance, we design the tapered structure only in the III/V structure, while keeping the SOI wire waveguide straight. In the proposed structure, the length of the bi-sectional tapered coupler can be 9.5 {\mu}m with a large fundamental mode-coupling ratio (over 95%) in a bandwidth of ~100 nm, and provides 100 nm tolerance to misalignment, even when the BCB layer is as thick as 50 nm.


title: Multi-Directional Multi-Level Dual-Cross Patterns for Robust Face Recognition

abstract: To perform unconstrained face recognition robust to variations in illumination, pose and expression, this paper presents a new scheme to extract "Multi-Directional Multi-Level Dual-Cross Patterns" (MDML-DCPs) from face images. Specifically, the MDMLDCPs scheme exploits the first derivative of Gaussian operator to reduce the impact of differences in illumination and then computes the DCP feature at both the holistic and component levels. DCP is a novel face image descriptor inspired by the unique textural structure of human faces. It is computationally efficient and only doubles the cost of computing local binary patterns, yet is extremely robust to pose and expression variations. MDML-DCPs comprehensively yet efficiently encodes the invariant characteristics of a face image from multiple levels into patterns that are highly discriminative of inter-personal differences but robust to intra-personal variations. Experimental results on the FERET, CAS-PERL-R1, FRGC 2.0, and LFW databases indicate that DCP outperforms the state-of-the-art local descriptors (e.g. LBP, LTP, LPQ, POEM, tLBP, and LGXP) for both face identification and face verification tasks. More impressively, the best performance is achieved on the challenging LFW and FRGC 2.0 databases by deploying MDML-DCPs in a simple recognition scheme.


title: The impact of a young radio galaxy: clues from the cosmic-ray electron population

abstract: In the framework of hierarchical structure formation AGN feedback shapes the galaxy luminosity function. Low luminosity, galaxy-scale double radio sources are ideal targets to investigate the interplay between AGN feedback and star formation. We use VLA and BIMA observations to study the radio continuum emission of NGC 3801 between 1.4 and 112.4 GHz. We find a prominent spectral break at ~10 GHz, where the spectrum steepens as expected from cosmic-ray electron (CRe) ageing. Using the equipartition magnetic field and fitting JP models locally we create a spatially resolved map of the spectral age of the CRe population. The spectral age of tau_int = 2.0 +/- 0.2 Myr agrees within a factor of two with the dynamical age of the expanding X-ray emitting shells. The spectral age varies only little across the lobes, requiring an effective mixing process of the CRe such as a convective backflow of magnetized plasma. The jet termination points have a slightly younger CRe spectral age, hinting at in-situ CRe re-acceleration. Our findings support the scenario where the supersonically expanding radio lobes heat the ISM of NGC 3801 via shock waves, and, as their energy is comparable to the energy of the ISM, are clearly able to influence the galaxy's further evolution.


title: Hot big bang or slow freeze?

abstract: We confront the big bang for the beginning of the universe with an equivalent picture of a slow freeze - a very cold and slowly evolving universe. In the freeze picture the masses of elementary particles increase and the gravitational constant decreases with cosmic time, while the Newtonian attraction remains unchanged. The freeze and big bang pictures both describe the same observations or physical reality. We present a simple "crossover model" without a big bang singularity. In the infinite past space-time is flat. Our model is compatible with present observations, describing the generation of primordial density fluctuations during inflation as well as the present transition to a dark energy dominated universe.


title: Why free markets die: An evolutionary perspective

abstract: Company mergers and acquisitions are often perceived to act as catalysts for corporate growth in free markets systems: it is conventional wisdom that those activities lead to better and more efficient markets. However, the broad adoption of this perception into corporate strategy is prone to result in a less diverse and more unstable environment, dominated by either very large or very small niche entities. We show here that ancestry, i.e. the cumulative history of mergers, is the key characteristic that encapsulates the diverse range of drivers behind mergers and acquisitions, across a range of industries and geographies. A long-term growth analysis reveals that entities which have been party to fewer mergers tend to grow faster than more highly acquisitive businesses.


title: The quaternionic commutator bracket and its implications

abstract: A quaternionic commutator bracket for position and momentum shows that the quaternionic wave function, \emph{viz.} $\widetilde{\psi}=(\frac{i}{c}\,\psi_0\,,\vec{\psi})$, represents a state of a particle with orbital angular momentum, $L=3\,\hbar$, resulting from the internal structure of the particle. This angular momentum can be attributed to spin of the particle. The vector $\vec{\psi}$, points along the direction of $\vec{L}$. When a charged particle is placed in an electromagnetic fields the interaction energy reveals that the magnetic moments interact with the electric and magnetic fields giving rise to terms similar to Aharonov-Bohm and Aharonov-Casher effects.


title: A Distributed Minimum Cut Approximation Scheme

abstract: In this paper, we study the problem of approximating the minimum cut in a distributed message-passing model, the CONGEST model. The minimum cut problem has been well-studied in the context of centralized algorithms. However, there were no known non-trivial algorithms in the distributed model until the recent work of Ghaffari and Kuhn. They gave algorithms for finding cuts of size $O(\epsilon^{-1}\lambda)$ and $(2+\epsilon)\lambda$ in $O(D)+\tilde{O}(n^{1/2+\epsilon})$ rounds and $\tilde{O}(D+\sqrt{n})$ rounds respectively, where $\lambda$ is the size of the minimum cut. This matches the lower bound they provided up to a polylogarithmic factor. Yet, no scheme that achieves $(1+\epsilon)$-approximation ratio is known. We give a distributed algorithm that finds a cut of size $(1+\epsilon)\lambda$ in $\tilde{O}(D+\sqrt{n})$ time, which is optimal up to polylogarithmic factors.


title: On the p-rank of Ext(A,B) for countable abelian groups A and B

abstract: In this note we show that the p-rank of Ext(A,B) for countable torsion-free abelian groups A and B is either countable or the size of the continuum.


title: Upper critical field as a probe for multiband superconductivity in bulk and interfacial STO

abstract: We investigate the temperature dependence of the upper critical field $H_{c 2} $ as a tool to probe the possible presence of multiband superconductivity at the interface of LAO/STO. The behaviour of $H_{c 2} $ can clearly indicate two-band superconductivity through its nontrivial temperature dependence. For the disorder scattering dominated two-dimensional LAO/STO interface we find a characteristic non-monotonic curvature of the $H_{c 2} (T) $. We also analyse the $H_{c 2} $ for multiband bulk STO and find similar behaviour.


title: Next-to-leading order improved perturbative QCD + saturation + hydrodynamics model for A+A collisions

abstract: We calculate initial conditions for the hydrodynamical evolution in ultrarelativistic heavy-ion collisions at the LHC and RHIC in an improved next-to-leading order perturbative QCD + saturation framework. Using viscous relativistic hydrodynamics, we show that we obtain a good simultaneous description of the centrality dependence of charged particle multiplicities, transverse momentum spectra and elliptic flow at the LHC and at RHIC. In particular, we discuss how the temperature dependence of the shear viscosity is constrained by these data.


title: Measurement of Anomalous Moments of Creation and Annihilation Operators

abstract: We propose a scheme for measurements of anomalous moments of creation and annihilation operators requiring only one beam splitter, only one photodetector and not many measurements with phase-shift $\varphi$ in homodyning. This has advantageous over scheme of Shchukin and Vogel [Phys. Rev. A 72, 043808 (2005)] which requires many photodetectors, many beam splitters and a large number of measurements with phase-shift $\varphi$ as a Fourier transform is calculated integrating over $\varphi$.


title: An Integer Programming Approach to UEP Coding for Multiuser Broadcast Channels

abstract: In this paper, an integer programming approach is introduced to construct Unequal Error Protection (UEP) codes for multiuser broadcast channels. We show that the optimal codes can be constructed that satisfy the integer programming bound. Based on the bound, we compute asymptotic code rate and perform throughput analysis for the degraded broadcast channel.


title: On simultaneous diophantine approximations to $\zeta(2)$ and $\zeta(3)$

abstract: We present a hypergeometric construction of rational approximations to $\zeta(2)$ and $\zeta(3)$ which allows one to demonstrate simultaneously the irrationality of each of the zeta values, as well as to estimate from below certain linear forms in 1, $\zeta(2)$ and $\zeta(3)$ with rational coefficients. A new notion of (simultaneous) diophantine exponent is introduced to formalise the arithmetic structure of these specific linear forms. Finally, the properties of this newer concept are studied and linked to the classical irrationality exponent and its generalisations given recently by S. Fischler.


title: Climate of Earth-like planets with high obliquity and eccentric orbits: implications for habitability conditions

abstract: We explore the effects of seasonal variability for the climate of Earth-like planets as determined by the two parameters polar obliquity and orbital eccentricity using a general circulation model of intermediate complexity. In the first part of the paper we examine the consequences of different values of obliquity and eccentricity for the spatio-temporal patterns of radiation and surface temperatures as well as for the main characteristics of the atmospheric circulations. In the second part we analyse the associated implications for the habitability of planets close to the outer edge of the habitable zone (HZ). This part of the paper focuses in particular on the multistability property of climate, i.e. the parallel existence of both an ice-free and an ice-covered climate state. Our results show that seasonal variability affects both the existence of and transitions between the two climate states. Moreover, our experiments reveal that planets with Earth-like atmospheres and high seasonal variability can have ice-free areas at much larger distance from the host star than planets without seasonal variability, which leads to a substantial expansion of the outer edge of the HZ. Sensitivity experiments exploring the role of azimuthal obliquity and surface heat capacity test the robustness of our results. On circular orbits, our findings obtained with a general circulation model agree well with previous studies based on one dimensional energy balance models, whereas significant differences are found on eccentric orbits.


title: Witness of mixed separable states useful for entanglement creation

abstract: Absolutely separable states form a special subset of the set of all separable states, as they remain separable under any global unitary transformation unlike other separable states. In this work we consider the set of absolutely separable bipartite states and show that it is convex and compact in any arbitrary dimensional Hilbert space. Through a generic approach of construction of suitable hermitian operators we prove the completeness of the separation axiom for identifying any separable state that is not absolutely separable. We demonstrate the action of such witness operators in different qudit systems. Examples of mixed separable systems are provided, pointing out the utility of the witness in entanglement creation using quantum gates. Decomposition of witnesses in terms of spin operators or photon polarizations facilitates their measureability for qubit states.


title: Game Semantics for Access Control

abstract: We introduce a semantic approach to the study of logics for access control and dependency analysis, based on Game Semantics. We use a variant of AJM games with explicit justification (but without pointers). Based on this, we give a simple and intuitive model of the information flow constraints underlying access control. This is used to give strikingly simple proofs of \emph{non-interference theorems} in robust, semantic versions.


title: Hydrodynamical simulations of coupled and uncoupled quintessence models II: Galaxy clusters

abstract: We study the $z=0$ properties of clusters (and large groups) of galaxies within the context of interacting and non-interacting quintessence cosmological models, using a series of adiabatic SPH simulations. Initially, we examine the average properties of groups and clusters, quantifying their differences in LCDM, uncoupled Dark Energy (\ude) and coupled Dark Energy (\cde) cosmologies. In particular, we focus upon radial profiles of the gas density, temperature and pressure, and we also investigate how the standard hydrodynamic equilibrium hypothesis holds in quintessence cosmologies. While we are able to confirm previous results about the distribution of baryons, we also find that the main discrepancy (with differences up to $20\%$) can be seen in cluster pressure profiles. We then switch attention to individual structures, mapping each halo in quintessence cosmology to its \LCDM\ counterpart. We are able to identify a series of small correlations between the coupling in the dark sector and halo spin, triaxiality and virialization ratio. When looking at spin and virialization of dark matter haloes, we find a weak ($5\%$) but systematic deviation in fifth force scenarios from \LCDM.


title: Compositional Operators in Distributional Semantics

abstract: This survey presents in some detail the main advances that have been recently taking place in Computational Linguistics towards the unification of the two prominent semantic paradigms: the compositional formal semantics view and the distributional models of meaning based on vector spaces. After an introduction to these two approaches, I review the most important models that aim to provide compositionality in distributional semantics. Then I proceed and present in more detail a particular framework by Coecke, Sadrzadeh and Clark (2010) based on the abstract mathematical setting of category theory, as a more complete example capable to demonstrate the diversity of techniques and scientific disciplines that this kind of research can draw from. This paper concludes with a discussion about important open issues that need to be addressed by the researchers in the future.


title: A microscopic quantal self-consistent cranking model for oscillations in spherical nuclei

abstract: In this article, we transform the previously-derived microscopic rotational-model Schrodinger equation into a form suitable for describing oscillations-coupled-to-intrinsic motion in spherical nuclei. The resulting equation is decomposed into two coupled cranking-type equations, one for the oscillation and another for the intrinsic motion, using a product wavefunction and a constrained variational method. The energy and cranking parameters in the coupled equations are self-consistently determined as functions of the system parameters by the solutions of the equations themselves. This self-consistency makes the two equations time-reversal invariant, unlike the conventional phenomenological cranking models. The self-consistency and time-reversal invariance accept only real solutions to the equations. For the harmonic oscillator mean-field potential, we explicitly determine these solutions and the corresponding eigenvalues, and derive the set of equations that determine self-consistently the parameters. To explore the relative importance of the various model features and approximations, we perform a preliminary scoping calculation of the excitation energy of the first excited states in the light nuclei using a sum rule to determine the oscillation frequency. The preliminary results indicate that, except in the lightest nuclei, the excitation energies are significantly overpredicted in the light nuclei due to the neglect, among other factors, of the deformation degree of freedom. The model derivation presented here serves as guide for eventually developing a corresponding model for the vibrational-rotational motion in deformed nuclei.


title: $SO(N)_2$ Braid group representations are Gaussian

abstract: We give a description of the centralizer algebras for tensor powers of spin objects in the pre-modular categories $SO(N)_2$ (for $N$ odd) and $O(N)_2$ (for $N$ even) in terms of quantum $(n-1)$-tori, via non-standard deformations of $U\mathfrak{so}_N$. As a consequence we show that the corresponding braid group representations are Gaussian representations, the images of which are finite groups. This verifies special cases of a conjecture that braid group representations coming from weakly integral braided fusion categories have finite image.


title: Study of Neural Network Algorithm for Straight-Line Drawings of Planar Graphs

abstract: Graph drawing addresses the problem of finding a layout of a graph that satisfies given aesthetic and understandability objectives. The most important objective in graph drawing is minimization of the number of crossings in the drawing, as the aesthetics and readability of graph drawings depend on the number of edge crossings. VLSI layouts with fewer crossings are more easily realizable and consequently cheaper. A straight-line drawing of a planar graph G of n vertices is a drawing of G such that each edge is drawn as a straight-line segment without edge crossings. However, a problem with current graph layout methods which are capable of producing satisfactory results for a wide range of graphs is that they often put an extremely high demand on computational resources. This paper introduces a new layout method, which nicely draws internally convex of planar graph that consumes only little computational resources and does not need any heavy duty preprocessing. Here, we use two methods: The first is self organizing map known from unsupervised neural networks which is known as (SOM) and the second method is Inverse Self Organized Map (ISOM).


title: Chandra and MMT observations of low-mass black hole active galactic nuclei accreting at low rates in dwarf galaxies

abstract: We report on Chandra X-ray observations of four candidate low-mass black hole (<10^6Msun) active galactic nuclei (AGNs) that have the estimated Eddington ratios among the lowest (~10^(-2)) found for this class. The aims are to validate the nature of their AGNs and to confirm the low Eddington ratios that are derived from the broad H_alpha line, and to explore this poorly studied regime in the AGN parameter space. Among them, two objects with the lowest significance of the broad lines are also observed with Multi-Mirror Telescope, and the high-quality optical spectra taken confirm them as Seyfert 1 AGNs and as having small black hole masses. X-ray emission is detected from the nuclei of two of the galaxies, which is variable on timescales of 10^3s, whereas no significant (or only marginal at best) detection is found for the remaining two. The X-ray luminosities are on the order of 10^(41) ergs/s or even lower, on the order of 10^(40) ergs/s for non-detections, which are among the lowest regimes ever probed for Seyfert galaxies. The low X-ray luminosities, compared to their black hole masses derived from H_alpha, confirm their low accretion rates assuming typical bolometric corrections. Our results hint at the existence of a possibly large population of under-luminous low-mass black holes in the local universe. An off-nucleus ultra-luminous X-ray source (ULX) in one of the dwarf galaxies is detected serendipitously, with a luminosity (6-9)x10^(39) ergs/s in 2-10 keV.


title: Active crystals and their stability

abstract: A recently introduced active phase field crystal model describes the formation of ordered resting and traveling crystals in systems of self-propelled particles. Increasing the active drive, a resting crystal can be forced to perform collectively ordered migration as a single traveling object. We demonstrate here that these ordered migrating structures are linearly stable. In other words, during migration, the single crystalline texture together with the globally ordered collective motion is preserved even on large length scales. Furthermore, we consider self-propelled particles on a substrate that are surrounded by a thin fluid film. We find that in this case the resulting hydrodynamic interactions can destabilize the order.


title: Simultaneous whole-animal 3D-imaging of neuronal activity using light field microscopy

abstract: 3D functional imaging of neuronal activity in entire organisms at single cell level and physiologically relevant time scales faces major obstacles due to trade-offs between the size of the imaged volumes, and spatial and temporal resolution. Here, using light-field microscopy in combination with 3D deconvolution, we demonstrate intrinsically simultaneous volumetric functional imaging of neuronal population activity at single neuron resolution for an entire organism, the nematode Caenorhabditis elegans. The simplicity of our technique and possibility of the integration into epi-fluoresence microscopes makes it an attractive tool for high-speed volumetric calcium imaging.


title: A Microkernel Architecture for Constraint Programming

abstract: This paper presents a microkernel architecture for constraint programming organized around a number of small number of core functionalities and minimal interfaces. The architecture contrasts with the monolithic nature of many implementations. Experimental results indicate that the software engineering benefits are not incompatible with runtime efficiency.


title: A contraction of the principal series by Berezin-Weyl quantization

abstract: We study a contraction of the principal series representations of a noncompact semisimple Lie group to the unitary irreducible representations of its Cartan motion group by means of the Berezin-Weyl quantization on the coadjoint orbits associated with these representations.


title: Signature, positive Hopf plumbing and the Coxeter transformation

abstract: By a theorem of A'Campo, the eigenvalues of certain Coxeter transformations are positive real or lie on the unit circle. By optimally bounding the signature of tree-like positive Hopf plumbings from below by the genus, we prove that at least two thirds of them lie on the unit circle. In contrast, we show that for divide links, the signature cannot be linearly bounded from below by the genus.


title: Understanding the gravitational and cosmological redshifts as Doppler shifts by gravitational phase factors

abstract: From the viewpoint of gauge gravitational theories, the path dependent gravitational phase factors define the Lorentz transformations between the local inertial coordinate systems of different positions. With this point we show that the spectral shifts in the curved spacetime, such as the gravitational and cosmological redshifts, can be understood as Doppler shifts. All these shifts are interpreted in a unified way as being originated from the relative motion of the free falling observers instantaneously static with the wave source and the receiver respectively. The gravitational phase factor of quantum systems in the curved spacetime is also discussed.


title: Ranking Templates for Linear Loops

abstract: We present a new method for the constraint-based synthesis of termination arguments for linear loop programs based on linear ranking templates. Linear ranking templates are parametrized, well-founded relations such that an assignment to the parameters gives rise to a ranking function. This approach generalizes existing methods and enables us to use templates for many different ranking functions with affine-linear components. We discuss templates for multiphase, piecewise, and lexicographic ranking functions. Because these ranking templates require both strict and non-strict inequalities, we use Motzkin's Transposition Theorem instead of Farkas Lemma to transform the generated $\exists\forall$-constraint into an $\exists$-constraint.


title: Complex Objects in the Polytopes of the Linear State-Space Process

abstract: A simple object (one point in $m$-dimensional space) is the resultant of the evolving matrix polynomial of walks in the irreducible aperiodic network structure of the first order DeGroot (weighted averaging) state-space process. This paper draws on a second order generalization the DeGroot model that allows complex object resultants, i.e, multiple points with distinct coordinates, in the convex hull of the initial state-space. It is shown that, holding network structure constant, a unique solution exists for the particular initial space that is a sufficient condition for the convergence of the process to a specified complex object. In addition, it is shown that, holding network structure constant, a solution exists for dampening values sufficient for the convergence of the process to a specified complex object. These dampening values, which modify the values of the walks in the network, control the system's outcomes, and any strongly connected typology is a sufficient condition of such control.


title: Energy and Information Near Black Hole Horizons

abstract: The central challenge in trying to resolve the firewall paradox is to identify excitations in the near-horizon zone of a black hole that can carry information without injuring a freely falling observer. By analyzing the problem from the point of view of a freely falling observer, I arrive at a simple proposal for the degrees of freedom that carry information out of the black hole. An infalling observer experiences the information-carrying modes as ingoing, negative energy excitations of the quantum fields. In these states, freely falling observers who fall in from infinity do not encounter a firewall, but freely falling observers who begin their free fall from a location close to the horizon are "frozen" by a flux of negative energy. When the black hole is "mined," the number of information-carrying modes increases, increasing the negative energy flux in the infalling frame without violating the equivalence principle. Finally, I point out a loophole in recent arguments that an infalling observer must detect a violation of unitarity, effective field theory, or free infall.


title: Domain Views for Constraint Programming

abstract: Views are a standard abstraction in constraint programming: They make it possible to implement a single version of each constraint, while avoiding to create new variables and constraints that would slow down propagation. Traditional constraint-programming systems provide the concept of {\em variable views} which implement a view of the type $y = f(x)$ by delegating all (domain and constraint) operations on variable $y$ to variable $x$. This paper proposes the alternative concept of {\em domain views} which only delegate domain operations. Domain views preserve the benefits of variable views but simplify the implementation of value-based propagation. Domain views also support non-injective views compositionally, expanding the scope of views significantly. Experimental results demonstrate the practical benefits of domain views.


title: Video-rate computational heterodyne holography

abstract: We present a versatile computational image rendering software of optically-acquired holograms. The reported software can process 4 Megapixel 8-bit raw frames from a sensor array acquired at a sustained rate of 80 Hz. Video-rate image rendering is achieved by streamline image processing with commodity computer graphics hardware. For time-averaged holograms acquired in off-axis optical configuration with a frequency-shifted reference beam, wide-field imaging of one tunable spectral component is permitted. This software is validated by phase-stepped hologram rendering, and non-contact monitoring of surface acoustic waves by single and dual sideband hologram rendering. It demonstrates the suitability of holography for video-rate computational laser Doppler imaging in heterodyne optical configuration.


title: Clustering South African households based on their asset status using latent variable models

abstract: The Agincourt Health and Demographic Surveillance System has since 2001 conducted a biannual household asset survey in order to quantify household socio-economic status (SES) in a rural population living in northeast South Africa. The survey contains binary, ordinal and nominal items. In the absence of income or expenditure data, the SES landscape in the study population is explored and described by clustering the households into homogeneous groups based on their asset status. A model-based approach to clustering the Agincourt households, based on latent variable models, is proposed. In the case of modeling binary or ordinal items, item response theory models are employed. For nominal survey items, a factor analysis model, similar in nature to a multinomial probit model, is used. Both model types have an underlying latent variable structure - this similarity is exploited and the models are combined to produce a hybrid model capable of handling mixed data types. Further, a mixture of the hybrid models is considered to provide clustering capabilities within the context of mixed binary, ordinal and nominal response data. The proposed model is termed a mixture of factor analyzers for mixed data (MFA-MD). The MFA-MD model is applied to the survey data to cluster the Agincourt households into homogeneous groups. The model is estimated within the Bayesian paradigm, using a Markov chain Monte Carlo algorithm. Intuitive groupings result, providing insight to the different socio-economic strata within the Agincourt region.


title: Holographic imaging of surface acoustic waves

abstract: We report on an experimental demonstration of surface acoustic waves monitoring on a thin metal plate with heterodyne optical holography. Narrowband imaging of local optical pathlength modulation is achieved with a frequency-tunable time-averaged laser Doppler holographic imaging scheme on a sensor array, at video-rate. This method enables robust and quantitative mapping of out-of-plane vibrations of nanometric amplitudes at radiofrequencies.


title: Congruences for the Fishburn Numbers

abstract: The Fishburn numbers, $\xi(n),$ are defined by a formal power series expansion $$ \sum_{n=0}^\infty \xi(n)q^n = 1 + \sum_{n=1}^\infty \prod_{j=1}^n (1-(1-q)^j). $$ For half of the primes $p$, there is a non--empty set of numbers $T(p)$ lying in $[0,p-1]$ such that if $j\in T(p),$ then for all $n\geq 0,$ $$ \xi(pn+j)\equiv 0 \pmod{p}. $$


title: Software Licensing in the Cloud Age

abstract: Cloud computing represents a major shift in information systems architecture, combining both new deployment models and new business models. Rapid provisioning, elastic scaling, and metered usage are essential characteristics of cloud services, and they require cloud resources with these same characteristics. When cloud services depend on commercial software, the licenses for that software become another resource to be managed by the cloud. This paper examines common licensing models, including open source, and how well they function in a cloud services model. It discusses creative, new, cloud-centric licensing models and how they allow providers to preserve and expand their revenue streams as their partners and customers transition to the cloud. The paper concludes by identifying the next steps to achieve standardized, cloud-friendly licensing models.


title: Linear Ranking for Linear Lasso Programs

abstract: The general setting of this work is the constraint-based synthesis of termination arguments. We consider a restricted class of programs called lasso programs. The termination argument for a lasso program is a pair of a ranking function and an invariant. We present the---to the best of our knowledge---first method to synthesize termination arguments for lasso programs that uses linear arithmetic. We prove a completeness theorem. The completeness theorem establishes that, even though we use only linear (as opposed to non-linear) constraint solving, we are able to compute termination arguments in several interesting cases. The key to our method lies in a constraint transformation that replaces a disjunction by a sum.


title: Analytical Solution of Mathieu Equation

abstract: The general solution of the homogeneous damped Mathieu equation in the analytical form, allowing its practical using in many applications, including superconductivity studies, without numerical calculations has been found.


title: Numerical study of the transverse stability of NLS soliton solutions in several classes of NLS type equations

abstract: Dispersive PDEs are important both in applications (wave phenomena e.g. in hy- drodynamics, nonlinear optics, plasma physics, Bose-Einstein condensates,...) and a mathematically very challenging class of partial differential equations, especially in the time dependent case. An important point with respect to applications is the stability of exact solutions like solitons. Whereas the linear or spectral stability can be addressed analytically in some situations, the proof of full nonlinear (in-)stability remains mostly an open question. In this paper, we numerically investi- gate the transverse (in-)stability of the solitonic solution to the one-dimensional cubic NLS equation, the well known isolated soliton, under the time evolution of several higher dimensional models, being admissible as a tranverse perturbation of the 1d cubic NLS. One of the recent work in this context [42] allowed to prove the instability of the soliton, under the flow of the classical (elliptic) 2d cubic NLS equation, for both localized or periodic perturbations. The characteristics of this instability stay however unknown. Is there a blow-up, dispersion..? We first illustrate how this instability occurs for the elliptic 2d cubic NLS equation and then show that the elliptic-elliptic Davey Stewartson system (a (2+1)-dimensional generalization of the cubic NLS equation) behaves as the former in this context. Then we investigate hyperbolic variants of the above models, for which no theory in this context is available. Namely we consider the hyperbolic 2d cubic NLS equation and the Davey-Stewartson II equations. For localized perturbations, the isolated soliton appears to be unstable for the former case, but seems to be orbitally stable for the latter. For periodic perturbations the soliton is found to be unstable for all transversally perturbed models considered.


title: Precision radial velocities of 15 M5 - M9 dwarfs

abstract: We present radial velocity measurements of a sample of M5V-M9V stars from our Red-Optical Planet Survey, ROPS, operating at 0.65-1.025 micron. Radial velocities for 15 stars, with r.m.s. precision down to 2.5 m/s over a week long time scale are achieved using Thorium-Argon reference spectra. We are sensitive to planets with m_psin(i) >= 1.5 MEarth (3 MEarth at 2-sigma) in the classical habitable zone and our observations currently rule out planets with m_psin(i) > 0.5 MJup at 0.03 AU for all our targets. A total of 9 of the 15 targets exhibit r.m.s. < 16 m/s, which enables us to rule out the presence of planets with m_psin(i) > 10 MEarth in 0.03 AU orbits. Since the mean rotation velocity is of order 8 km/s for an M6V star and 15 km/s by M9V, we avoid observing only slow rotators that would introduce a bias towards low axial inclination i << 90 deg systems, which are unfavourable for planet detection. Our targets with the highest vsini values exhibit radial velocities significantly above the photon-noise limited precision, even after accounting for vsini. We therefore monitored stellar activity via chromospheric emission from the Halpha and Ca II infrared triplet lines. A clear trend of log_10(Halpha/Lbol) with radial velocity r.m.s. is seen, implying that significant starspot activity is responsible for the observed radial velocity precision floor. The implication that most late M dwarfs are significantly spotted, and hence exhibit time varying line distortions, indicates that observations to detect orbiting planets need strategies to reliably mitigate against the effects of activity induced radial velocity variations.


title: Ranking Function Synthesis for Linear Lasso Programs

abstract: The scope of this work is the constraint-based synthesis of termination arguments for the restricted class of programs called linear lasso programs. A termination argument consists of a ranking function as well as a set of supporting invariants. We extend existing methods in several ways. First, we use Motzkin's Transposition Theorem instead of Farkas' Lemma. This allows us to consider linear lasso programs that can additionally contain strict inequalities. Existing methods are restricted to non-strict inequalities and equalities. Second, we consider several kinds of ranking functions: affine-linear, piecewise and lexicographic ranking functions. Moreover, we present a novel kind of ranking function called multiphase ranking function which proceeds through a fixed number of phases such that for each phase, there is an affine-linear ranking function. As an abstraction to the synthesis of specific ranking functions, we introduce the notion ranking function template. This enables us to handle all ranking functions in a unified way. Our method relies on non-linear algebraic constraint solving as a subroutine which is known to scale poorly to large problems. As a mitigation we formalize an assessment of the difficulty of our constraints and present an argument why they are of an easier kind than general non-linear constraints. We prove our method to be complete: if there is a termination argument of the form specified by the given ranking function template with a fixed number of affine-linear supporting invariants, then our method will find a termination argument. To our knowledge, the approach we propose is the most powerful technique of synthesis-based discovery of termination arguments for linear lasso programs and encompasses and enhances several methods having been proposed thus far.


title: Effective Casimir Conditions and Group Coherent States

abstract: Properties of group coherent states can be derived "effectively" without knowing full wave functions. The procedure is detailed in this article as an example of general methods for effective constraints. The role of constraints in the present context is played by a Casimir condition that puts states within an irreducible representation of a Lie group (or, equivalently, on a quantization of a co-adjoint orbit of the dual Lie algebra). Simplifications implied by a Casimir condition, compared with general first-class constraints, allows one to show that the correct number of degrees of freedom is obtained after imposing the condition. When combined with conditions to saturate uncertainty relations, moments of group coherent states can be derived. A detailed example in quantum cosmology (cosmic forgetfulness) illustrates the usefulness of the methods.


title: STABLAB Documentation for KdV : Numerical proof of stability of roll waves in the small-amplitude limit for inclined thin film flow

abstract: We document the MATLAB code used in the following study: Numerical proof of stability of roll waves in the small-amplitude limit for inclined thin film flow.


title: Counting Houses of Pareto Optimal Matchings in the House Allocation Problem

abstract: Let $A,B$ with $|A| = m$ and $|B| = n\ge m$ be two sets. We assume that every element $a\in A$ has a reference list over all elements from $B$. We call an injective mapping $\tau$ from $A$ to $B$ a matching. A blocking coalition of $\tau$ is a subset $A'$ of $A$ such that there exists a matching $\tau'$ that differs from $\tau$ only on elements of $A'$, and every element of $A'$ improves in $\tau'$, compared to $\tau$ according to its preference list. If there exists no blocking coalition, we call the matching $\tau$ an exchange stable matching (ESM). An element $b\in B$ is reachable if there exists an exchange stable matching using $b$. The set of all reachable elements is denoted by $E^*$. We show \[|E^*| \leq \sum_{i = 1,\ldots, m}{\left\lfloor\frac{m}{i}\right\rfloor} = \Theta(m\log m).\] This is asymptotically tight. A set $E\subseteq B$ is reachable (respectively exactly reachable) if there exists an exchange stable matching $\tau$ whose image contains $E$ as a subset (respectively equals $E$). We give bounds for the number of exactly reachable sets. We find that our results hold in the more general setting of multi-matchings, when each element $a$ of $A$ is matched with $\ell_a$ elements of $B$ instead of just one. Further, we give complexity results and algorithms for corresponding algorithmic questions. Finally, we characterize unavoidable elements, i.e., elements of $B$ that are used by all ESM's. This yields efficient algorithms to determine all unavoidable elements.


title: Recovering the chiral critical end-point via delocalization of quark interactions

abstract: We show that for the lower branch of the quark condensate and values higher than approximately $-(250 \, \mathrm{MeV})^3$ the chiral critical end-point in the Nambu--Jona-Lasinio model does not occur in the phase diagram. By using lattice motivated non-local quark interactions, we demonstrate that the critical end-point can be recovered. We study this behavior for a range of condensate values and find that the variation in the position of the critical end-point is more pronounced as the condensate is increased.


title: Random strain fluctuations as dominant disorder source for high-quality on-substrate graphene devices

abstract: We have performed systematic investigations of transport through graphene on hexagonal boron nitride (hBN) substrates, together with confocal Raman measurements and a targeted theoretical analysis, to identify the dominant source of disorder in this system. Low-temperature transport measurements on many devices reveal a clear correlation between the carrier mobility $\mu$ and the width $n^*$ of the resistance peak around charge neutrality, demonstrating that charge scattering and density inhomogeneities originate from the same microscopic mechanism. The study of weak-localization unambiguously shows that this mechanism is associated to a long-ranged disorder potential, and provides clear indications that random pseudo-magnetic fields due to strain are the dominant scattering source. Spatially resolved Raman spectroscopy measurements confirm the role of local strain fluctuations, since the line-width of the Raman 2D-peak --containing information of local strain fluctuations present in graphene-- correlates with the value of maximum observed mobility. The importance of strain is corroborated by a theoretical analysis of the relation between $\mu$ and $n^*$ that shows how local strain fluctuations reproduce the experimental data at a quantitative level, with $n^*$ being determined by the scalar deformation potential and $\mu$ by the random pseudo-magnetic field (consistently with the conclusion drawn from the analysis of weak-localization). Throughout our study, we compare the behavior of devices on hBN substrates to that of devices on SiO$_2$ and SrTiO$_3$, and find that all conclusions drawn for the case of hBN are compatible with the observations made on these other materials. These observations suggest that random strain fluctuations are the dominant source of disorder for high-quality graphene on many different substrates, and not only on hexagonal boron nitride.


title: Shape, Thermal and Surface Properties determination of a Candidate Spacecraft Target Asteroid (175706) 1996 FG3

abstract: In this paper, a 3D convex shape model of (175706) 1996 FG3, which consists of 2040 triangle facets and 1022 vertices, is derived from the known lightcurves. The best-fit orientation of the asteroid's spin axis is determined to be $\lambda =237.7^\circ$ and $\beta=-83.8^{\circ}$ considering the observation uncertainties, and its rotation period is $\sim$ 3.5935 h . Using the derived shape model, we adopt the so-called advanced thermophysical model (ATPM) to fit three published sets of mid-infrared observations of 1996 FG3 \citep{Wolters2011,Walsh2012}, so as to evaluate its surface properties. Assuming the primary and the secondary bear identical shape, albedo, thermal inertia and surface roughness, the best-fit parameters are obtained from the observations. The geometric albedo and effective diameter of the asteroid are reckoned to be $p_{\rm v}=0.045\pm0.002$, $D_{\rm eff}=1.69^{+0.05}_{-0.02}$ km. The diameters of the primary and secondary are determined to be $D_{1}=1.63^{+0.04}_{-0.03}$ km and $D_{2}=0.45^{+0.04}_{-0.03}$ km, respectively. The surface thermal inertia $\Gamma$ is derived to be a low value of $80\pm40\rm Jm^{-2}s^{-0.5}K^{-1}$ with a roughness fraction $f_{\rm R}$ of $0.8^{+0.2}_{-0.4}$. This indicates that the primary possibly has a regolith layer on its surface, which is likely to be covered by a mixture of dust, fragmentary rocky debris and sand. The minimum regolith depth is estimated to be $5\sim20\rm mm$ from the simulations of subsurface temperature distribution, indicating that 1996 FG3 could be a very suitable target for a sample return mission.


title: Experimentally-driven approach for measuring quantum phase uncertainty

abstract: We propose a new generalised formalism for estimating the quantum phase uncertainty of pure and mixed continuous-variable quantum states and compare this with the phase uncertainty given by the quantum Fisher information. In order to preserve the Hermiticity of the operators, we use the Heisenberg and Schroedinger uncertainty relations to derive expressions for the phase uncertainty from generalised Susskind-Glogower operators. This formalism not only offers the possibility of directly measuring quantum phase uncertainties in a cavity-QED experiment but also gives a significant computational saving over the quantum Fisher information approach, which requires diagonalisation of the density matrix.


title: Slowly varying control parameters, delayed bifurcations and the stability of spikes in reaction-diffusion systems

abstract: We present three examples of delayed bifurcations for spike solutions of reaction-diffusion systems. The delay effect results as the system passes slowly from a stable to an unstable regime, and was previously analysed in the context of ODE's in [P.Mandel, T.Erneux, J.Stat.Phys, 1987]. It was found that the instability would not be fully realized until the system had entered well into the unstable regime. The bifurcation is said to have been "delayed" relative to the threshold value computed directly from a linear stability analysis. In contrast, we analyze the delay effect in systems of PDE's. In particular, for spike solutions of singularly perturbed generalized Gierer-Meinhardt (GM) and Gray-Scott (GS) models, we analyze three examples of delay resulting from slow passage into regimes of oscillatory and competition instability. In the first example, for the GM model on the infinite real line, we analyze the delay resulting from slowly tuning a control parameter through a Hopf bifurcation. In the second example, we consider a Hopf bifurcation on a finite one-dimensional domain. In this scenario, as opposed to the extrinsic tuning of a system parameter through a bifurcation value, we analyze the delay of a bifurcation triggered by slow intrinsic dynamics of the PDE system. In the third example, we consider competition instabilities of the GS model triggered by the extrinsic tuning of a feed rate parameter. In all cases, we find that the system must pass well into the unstable regime before the onset of instability is fully observed, indicating delay. We also find that delay has an important effect on the eventual dynamics of the system in the unstable regime. We give analytic predictions for the magnitude of the delays as obtained through analysis of certain explicitly solvable nonlocal eigenvalue problems. The theory is confirmed by numerical solutions of the full PDE systems.


title: Positivity, Discontinuity, Finite Resources and Nonzero Error for Arbitrarily Varying Quantum Channels

abstract: This work is motivated by a quite general question: Under which circumstances are the capacities of information transmission systems continuous? The research is explicitly carried out on arbitrarily varying quantum channels (AVQCs). We give an explicit example that answers the recent question whether the transmission of messages over AVQCs can benefit from distribution of randomness between the legitimate sender and receiver in the affirmative. The specific class of channels introduced in that example is then extended to show that the deterministic capacity does have discontinuity points, while that behaviour is, at the same time, not generic: We show that it is continuous around its positivity points. This is in stark contrast to the randomness-assisted capacity, which is always continuous in the channel. Our results imply that the deterministic message transmission capacity of an AVQC can be discontinuous only in points where it is zero, while the randomness assisted capacity is nonzero. Apart from the zero-error capacities, this is the first result that shows a discontinuity of a capacity for a large class of quantum channels. The continuity of the respective capacity for memoryless quantum channels had, among others, been listed as an open problem on the problem page of the ITP Hannover for about six years before it was proven to be continuous. We also quantify the interplay between the distribution of finite amounts of randomness between the legitimate sender and receiver, the (nonzero) decoding error with respect to the average error criterion that can be achieved over a finite number of channel uses and the number of messages that can be sent. This part of our results also applies to entanglement- and strong subspace transmission. In addition, we give a new sufficient criterion for the entanglement transmission capacity with randomness assistance to vanish.


title: Search for Majorana neutrinos in $B^- \to \pi^+\mu^-\mu^-$ decays

abstract: A search for heavy Majorana neutrinos produced in the $B^- \to \pi^+\mu^-\mu^-$ decay mode is performed using 3 fb$^{-1}$ of integrated luminosity collected with the LHCb detector in $pp$ collisions at center-of-mass energies of 7 TeV and 8 TeV at the LHC. Neutrinos with masses in the range 250-5000 MeV and lifetimes from zero to 1000 ps are probed. In the absence of a signal, upper limits are set on the branching fraction ${\cal{B}}(B^- \to \pi^+\mu^-\mu^-)$ as functions of neutrino mass and lifetime. These limits are on the order of $10^{-9}$ for short neutrino lifetimes of 1 ps or less. Limits are also set on the coupling between the muon and a possible fourth-generation neutrino.


title: Cohomology of deformations

abstract: In this article we study cohomology of a group with coefficients in representations on Banach spaces and its stability under deformations. We show that small, metric deformations of the representation preserve vanishing of cohomology. As applications we obtain deformation theorems for fixed point properties on Banach spaces. In particular, our results yield fixed point theorems for affine actions in which the linear part is not uniformly bounded. Our proofs are effective and allow for quantitative estimates.


title: On Estimation of Optimal Process Parameters of the Abrasive Waterjet Machining

abstract: This paper is a commentary on the work of Zain et al. (2011a) in which they proposed the use of soft computing techniques, i.e., genetic algorithm (GA) and simulated annealing (SA), to estimate optimal process parameters of the abrasive waterjet machining. We found a fundamental problem with their work; the authors have misinterpreted the original work of Caydas and Hascalik (2008) in which the real problem is how to build predictive models for the machining process, not to optimize the model variables using the estimate model as presented by Zain et al. (2011a). In addition to this problem, the objective function to be optimized is an infinitely differentiable quadratic formula with box constraints, thus instead of using soft computing techniques, gradient methods that have better convergence property and are usually faster should be employed first. In particular, we show that the active set method outperforms GA and SA both in the machining performance and computational times (about 13 times faster than GA and 462 times faster than SA).


title: HMACA: Towards Proposing a Cellular Automata Based Tool for Protein Coding, Promoter Region Identification and Protein Structure Prediction

abstract: Human body consists of lot of cells, each cell consist of DeOxaRibo Nucleic Acid (DNA). Identifying the genes from the DNA sequences is a very difficult task. But identifying the coding regions is more complex task compared to the former. Identifying the protein which occupy little place in genes is a really challenging issue. For understating the genes coding region analysis plays an important role. Proteins are molecules with macro structure that are responsible for a wide range of vital biochemical functions, which includes acting as oxygen, cell signaling, antibody production, nutrient transport and building up muscle fibers. Promoter region identification and protein structure prediction has gained a remarkable attention in recent years. Even though there are some identification techniques addressing this problem, the approximate accuracy in identifying the promoter region is closely 68% to 72%. We have developed a Cellular Automata based tool build with hybrid multiple attractor cellular automata (HMACA) classifier for protein coding region, promoter region identification and protein structure prediction which predicts the protein and promoter regions with an accuracy of 76%. This tool also predicts the structure of protein with an accuracy of 80%.


title: Markovianity of the invariant distribution of probabilistic cellular automata on the line

abstract: We revisit the problem of finding the conditions under which synchronous probabilistic cellular automata indexed by the line $\mathbb{Z}$, or the periodic line $\cyl{n}$, depending on 2 neighbours, admit as invariant distribution the law of a space-indexed Markov chain. Our advances concerns PCA defined on a finite alphabet, where most of existing results concern size 2 alphabet. A part of the paper is also devoted to the comparison of different structures ($\mathbb{Z}$, $\cyl{n}$, and also some structures constituted with two consecutive lines of the space time diagram) with respect to the property to possess a Markovian invariant distribution.


title: Nucleon Electric Dipole Moments and the Isovector Parity- and Time-Reversal-Odd Pion-Nucleon Coupling

abstract: The isovector time-reversal- and parity-violating pion-nucleon coupling $\bar{g}_{\pi}^{(1)}$ is uniquely sensitive to dimension six interactions between right-handed light quarks and the Standard Model Higgs doublet that naturally arises in left-right symmetric models. Recent work has used the $\bar{g}_{\pi}^{(1)}$-induced one-loop contribution to the neutron electric dipole moment $d_n$, together with the present experimental $d_n$ bound, to constrain the CP-violating parameters of the left-right symmetric model. We show that this and related analyses are based on an earlier meson theory $d_n$ computation that is not consistent with the power-counting appropriate for an effective field theory. We repeat the one-loop calculation using heavy baryon chiral perturbation theory and find that the resulting $d_n$ sensitivity to $\bar{g}_{\pi}^{(1)}$ is suppressed, implying more relaxed constraints on the parameter space of the left-right symmetric model. Assuming no cancellations between this loop contribution and other contributions, such as the leading order EDM low-energy constant, the present limit on $d_n$ implies $|\bar{g}_{\pi}^{(1)}| \lesssim 1.1\times 10^{-10}$.


title: Towards a Benchmark and a Comparison Framework for Combinatorial Interaction Testing of Software Product Lines

abstract: As Software Product Lines (SPLs) are becoming a more pervasive development practice, their effective testing is becoming a more important concern. In the past few years many SPL testing approaches have been proposed, among them, are those that support Combinatorial Interaction Testing (CIT) whose premise is to select a group of products where faults, due to feature interactions, are more likely to occur. Many CIT techniques for SPL testing have been put forward; however, no systematic and comprehensive comparison among them has been performed. To achieve such goal two items are important: a common benchmark of feature models, and an adequate comparison framework. In this research-in-progress paper, we propose 19 feature models as the base of a benchmark, which we apply to three different techniques in order to analyze the comparison framework proposed by Perrouin et al. We identify the shortcomings of this framework and elaborate alternatives for further study.


title: On the equivalence of two fundamental theta identities

abstract: Two fundamental theta identities, a three-term identity due to Weierstrass and a five-term identity due to Jacobi, both with products of four theta functions as terms, are shown to be equivalent. One half of the equivalence was already proved by R.J. Chapman in 1996. The history and usage of the two identities, and some generalizations are also discussed.


title: A limit equation criterion for applying the conformal method to asymptotically cylindrical initial data sets

abstract: We prove that in a certain class of conformal data on an asymptotically cylindrical manifold, if the conformally decomposed Einstein constraint equations do not admit a solution, then one can always find a nontrivial solution to the limit equation first explored by Dahl, Gicquaud, and Humbert in [DGH11]. We also give an example of a Ricci curvature condition on the manifold which precludes the existence of a solution to this limit equation, showing that such a limit criterion can be a useful tool for studying the Einstein constraint equations on manifolds with asymptotically cylindrical ends.


title: Classification of invariant valuations on the quaternionic plane

abstract: We describe the orbit space of the action of the group $\mathrm{Sp}(2)\mathrm{Sp}(1)$ on the real Grassmann manifolds $\mathrm{Gr}_k(\mathbb{H}^2)$ in terms of certain quaternionic matrices of Moore rank not larger than $2$. We then give a complete classification of valuations on the quaternionic plane $\mathbb{H}^2$ which are invariant under the action of the group $\mathrm{Sp}(2)\mathrm{Sp}(1)$.


title: Numerical proof of stability of roll waves in the small-amplitude limit for inclined thin film flow

abstract: We present a rigorous numerical proof based on interval arithmetic computations categorizing the linearized and nonlinear stability of periodic viscous roll waves of the KdV-KS equation modeling weakly unstable flow of a thin fluid film on an incline in the small-amplitude KdV limit. The argument proceeds by verification of a stability condition derived by Bar-Nepomnyashchy and Johnson-Noble-Rodrigues-Zumbrun involving inner products of various elliptic functions arising through the KdV equation. One key point in the analysis is a bootstrap argument balancing the extremely poor sup norm bounds for these functions against the extremely good convergence properties for analytic interpolation in order to obtain a feasible computation time. Another is the way of handling analytic interpolation in several variables by a two-step process carving up the parameter space into manageable pieces for rigorous evaluation. These and other general aspects of the analysis should serve as blueprints for more general analyses of spectral stability.


title: An integrated mid-infrared, far-infrared and terahertz optical Hall effect instrument

abstract: We report on the development of the first integrated mid-infrared, far-infrared and terahertz optical Hall effect instrument, covering an ultra wide spectral range from 3 cm$^{-1}$ to 7000 cm$^{-1}$ (0.1-210 THz or 0.4-870 meV). The instrument comprises four sub-systems, where the magneto-cryostat-transfer sub-system enables the usage of the magneto-cryostat sub-system with the mid-infrared ellipsometer sub-system, and the far-infrared/terahertz ellipsometer sub-system. Both ellipsometer sub-systems can be used as variable angle-of-incidence spectroscopic ellipsometers in reflection or transmission mode, and are equipped with multiple light sources and detectors. The ellipsometer sub-systems are operated in polarizer-sample-rotating-analyzer configuration granting access to the upper left $3\times 3$ block of the normalized $4\times 4$ Mueller matrix. The closed cycle magneto-cryostat sub-system provides sample temperatures between room temperature and 1.4 K and magnetic fields up to 8 T, enabling the detection of transverse and longitudinal magnetic field-induced birefringence. We discuss theoretical background and practical realization of the integrated mid-infrared, far-infrared and terahertz optical Hall effect instrument, as well as acquisition of optical Hall effect data and the corresponding model analysis procedures. Exemplarily, epitaxial graphene grown on 6H-SiC, a tellurium doped bulk GaAs sample and an AlGaN/GaN high electron mobility transistor structure are investigated. The selected experimental datasets display the full spectral, magnetic field and temperature range of the instrument and demonstrate data analysis strategies. Effects from free charge carriers in two dimensional confinement and in a volume material, as well as quantum mechanical effects (inter-Landau-level transitions) are observed and discussed exemplarily.


title: Local a Priori Estimate on the General Scale Subdomains

abstract: The local a priori estimate for the finite element approximation is essential for underlying the local and parallel technique. It is well known that the constant coefficients in the inequality is independent of the mesh size. But it is not so clear whether the constant depends on the scale of the local subdomains. The aim of this note is to derive a new local a priori estimate on the general scale domains. We also show that the dependence of the constant appearing in the local a priori estimate on the scale of the subdomains.


title: Formation and appearance of pulsar-like white dwarfs

abstract: Accretion-driven spin-up of a magnetized white dwarf in a close binary system is discussed. We address a situation in which the magnetic field of the white dwarf is screening during the accretion phase and re-generating due to the field diffusion through the accreted material after it. We find this scenario to be effective for a formation of massive pulsar-like white dwarfs.


title: Iterative regularization for ensemble data assimilation in reservoir models

abstract: We propose the application of iterative regularization for the development of ensemble methods for solving Bayesian inverse problems. In concrete, we construct (i) a variational iterative regularizing ensemble Levenberg-Marquardt method (IR-enLM) and (ii) a derivative-free iterative ensemble Kalman smoother (IR-ES). The aim of these methods is to provide a robust ensemble approximation of the Bayesian posterior. The proposed methods are based on fundamental ideas from iterative regularization methods that have been widely used for the solution of deterministic inverse problems [21]. In this work we are interested in the application of the proposed ensemble methods for the solution of Bayesian inverse problems that arise in reservoir modeling applications. The proposed ensemble methods use key aspects of the regularizing Levenberg-Marquardt scheme developed by Hanke [16] and that we recently applied for history matching in [18]. In the case where the forward operator is linear and the prior is Gaussian, we show that the proposed IR-enLM and IR-ES coincide with standard randomized maximum likelihood (RML) and the ensemble smoother (ES) respectively. For the general nonlinear case, we develop a numerical framework to assess the performance of the proposed ensemble methods at capturing the posterior. This framework consists of using a state-of-the art MCMC method for resolving the Bayesian posterior from synthetic experiments. The resolved posterior via MCMC then provides a gold standard against to which compare the proposed IR-enLM and IR-ES. We show that for the careful selection of regularization parameters, robust approximations of the posterior can be accomplished in terms of mean and variance. Our numerical experiments showcase the advantage of using iterative regularization for obtaining more robust and stable approximation of the posterior than standard unregularized methods.


title: Remarks on geometric properties of SQG sharp fronts and $\alpha$-patches

abstract: Guided by numerical simulations, we present the proof of two results concerning the behaviour of SQG sharp fronts and $\alpha$-patches. We establish that ellipses are not rotational solutions and we prove that initially convex interfaces may lose this property in finite time.


title: Trigonometric Parallaxes of High Mass Star Forming Regions: the Structure and Kinematics of the Milky Way

abstract: Over 100 trigonometric parallaxes and proper motions for masers associated with young, high-mass stars have been measured with the BeSSeL Survey, a VLBA key science project, the EVN, and the Japanese VERA project. These measurements provide strong evidence for the existence of spiral arms in the Milky Way, accurately locating many arm segments and yielding spiral pitch angles ranging from 7 to 20 degrees. The widths of spiral arms increase with distance from the Galactic center. Fitting axially symmetric models of the Milky Way with the 3-D position and velocity information and conservative priors for the solar and average source peculiar motions, we estimate the distance to the Galactic center, Ro, to be 8.34 +/- 0.16 kpc, a circular rotation speed at the Sun, To, to be 240 +/- 8 km/s, and a rotation curve that is nearly flat (a slope of -0.2 +/- 0.4 km/s/kpc) between Galactocentric radii of 5 and 16 kpc. Assuming a "universal" spiral galaxy form for the rotation curve, we estimate the thin disk scale length to be 2.44 +/- 0.16 kpc. The parameters Ro and To are not highly correlated and are relatively insensitive to different forms of the rotation curve. Adopting a theoretically motivated prior that high-mass star forming regions are in nearly circular Galactic orbits, we estimate a global solar motion component in the direction of Galactic rotation, Vsun = 14.6 +/- 5.0 km/s. While To and Vsun are significantly correlated, the sum of these parameters is well constrained, To + Vsun = 255.2 +/- 5.1 km/s, as is the angular speed of the Sun in its orbit about the Galactic center, (To + Vsun)/Ro = 30.57 +/- 0.43 km/s/kpc. These parameters improve the accuracy of estimates of the accelerations of the Sun and the Hulse-Taylor binary pulsar in their Galactic orbits, significantly reducing the uncertainty in tests of gravitational radiation predicted by general relativity.


title: A Multigrid Method Based On Shifted-Inverse Power Technique for Eigenvalue Problems

abstract: A multigrid method is proposed in this paper to solve eigenvalue problems by the finite element method based on the shifted-inverse power iteration technique. With this scheme, solving eigenvalue problem is transformed to a series of nonsingular solutions of boundary value problems on multilevel meshes. Since replacing the difficult eigenvalue solving by the easier solution of boundary value problems, the multigrid way can improve the overall efficiency of the eigenvalue problem solving. Some numerical experiments are presented to validate the efficiency of this new method.