Quench dynamics and parity blocking in Majorana wires

We theoretically explore quench dynamics in a finite-sized topological fermionic p-wave superconducting wire with the goal of demonstrating that topological order can have marked effects on such non-equilibrium dynamics. In the case studied here, topological order is reflected in the presence of two (nearly) isolated Majorana fermionic end bound modes together forming an electronic state that can be occupied or not, leading to two (nearly) degenerate ground states characterized by fermion parity. Our study begins with a characterization of the static properties of the finite-sized wire, including the behavior of the Majorana end modes and the form of the tunnel coupling between them; a transfer matrix approach to analytically determine the locations of the zero energy contours where this coupling vanishes; and a Pfaffian approach to map the ground state parity in the associated phase diagram. We next study the quench dynamics resulting from initializing the system in a topological ground state and then dynamically tuning one of the parameters of the Hamiltonian. For this, we develop a dynamic quantum many-body technique that invokes a Wick's theorem for Majorana fermions, vastly reducing the numerical effort given the exponentially large Hilbert space. We investigate the salient and detailed features of two dynamic quantities-the overlap between the time-evolved state and the instantaneous ground state (adiabatic fidelity) and the residual energy. When the parity of the instantaneous ground state flips successively with time, we find that the time-evolved state can dramatically switch back and forth between this state and an excited state even when the quenching is very slow, a phenomenon that we term `parity blocking'. This parity blocking becomes prominently manifest as non-analytic jumps as a function of time in both dynamic quantities.

Non-Abelian gauge redundancy and entropic ambiguities

The von Neumann entropy of a generic quantum state is not unique unless the state can be uniquely decomposed as a sum of extremal or pure states. Therefore one reaches the remarkable possibility that there may be many entropies for a given state. We show that this happens if the GNS representation (of the algebra of observables in some quantum state) is reducible, and some representations in the decomposition occur with non-trivial degeneracy. This ambiguity in entropy, which can occur at zero temperature, can often be traced to a gauge symmetry emergent from the non-trivial topological character of the configuration space of the underlying system. We also establish the analogue of an H-theorem for this entropy by showing that its evolution is Markovian, determined by a stochastic matrix. After demonstrating this entropy ambiguity for the simple example of the algebra of 2 x 2 matrices, we argue that the degeneracies in the GNS representation can be interpreted as an emergent broken gauge symmetry, and play an important role in the analysis of emergent entropy due to non-Abelian anomalies. We work out the simplest situation with such non-Abelian symmetry, that of an ethylene molecule.

A stochasticity threshold in holography and the instability of AdS

We give strong numerical evidence that a self-interacting probe scalar field in AdS, with only a few modes turned on initially, will undergo fast thermalization only if it is above a certain energetic threshold. Below the threshold the energy stays close to constant in a few modes for a very long time instead of cascading quickly. This indicates the existence of a Strong Stochasticity Threshold (SST) in holography. The idea of SST is familiar from certain statistical mechanical systems, and we suggest that it exists also in AdS gravity. This would naturally reconcile the generic nonlinear instability of AdS observed by Bizon and Rostworowski, with the Fermi-Pasta-Ulam-Tsingou-like quasiperiodicity noticed recently for some classes of initial conditions. We show that our simple setup captures many of the relevant features of the full gravity-scalar system.

Two player game variant of the Erdos-Szekeres problem

The classical Erdos-Szekeres theorem states that a convex k-gon exists in every sufficiently large point set. This problem has been well studied and finding tight asymptotic bounds is considered a challenging open problem. Several variants of the Erdos-Szekeres problem have been posed and studied in the last two decades. The well studied variants include the empty convex k-gon problem, convex k-gon with specified number of interior points and the chromatic variant. In this paper, we introduce the following two player game variant of the Erdos-Szekeres problem: Consider a two player game where each player playing in alternate turns, place points in the plane. The objective of the game is to avoid the formation of the convex k-gon among the placed points. The game ends when a convex k-gon is formed and the player who placed the last point loses the game. In our paper we show a winning strategy for the player who plays second in the convex 5-gon game and the empty convex 5-gon game by considering convex layer configurations at each step. We prove that the game always ends in the 9th step by showing that the game reaches a specific set of configurations.

Separation index of graphs and stacked 2-spheres

In 1987, Kalai proved that stacked spheres of dimension d >= 3 are characterised by the fact that they attain equality in Barnette's celebrated Lower Bound Theorem. This result does not extend to dimension d = 2. In this article, we give a characterisation of stacked 2-spheres using what we call the separation index. Namely, we show that the separation index of a triangulated 2-sphere is maximal if and only if it is stacked. In addition, we prove that, amongst all n-vertex triangulated 2-spheres, the separation index is minimised by some n-vertex flag sphere for n >= 6. Furthermore, we apply this characterisation of stacked 2-spheres to settle the outstanding 3-dimensional case of the Lutz-Sulanke-Swartz conjecture that ``tight-neighbourly triangulated manifolds are tight''. For dimension d >= 4, the conjecture has already been proved by Effenberger following a result of Novik and Swartz. (C) 2015 Elsevier Inc. All rights reserved.

Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis

The ultimate bearing capacity of a circular footing, placed over rock mass, is evaluated by using the lower bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization. The generalized Hoek-Brown (HB) failure criterion, but by keeping a constant value of the exponent, alpha = 0.5, was used. The failure criterion was smoothened both in the meridian and pi planes. The nonlinear optimization was carried out by employing an interior point method based on the logarithmic barrier function. The results for the obtained bearing capacity were presented in a non-dimensional form for different values of GSI, m(i), sigma(ci)/(gamma b) and q/sigma(ci). Failure patterns were also examined for a few cases. For validating the results, computations were also performed for a strip footing as well. The results obtained from the analysis compare well with the data reported in literature. Since the equilibrium conditions are precisely satisfied only at the centroids of the elements, not everywhere in the domain, the obtained lower bound solution will be approximate not true. (C) 2015 Elsevier Ltd. All rights reserved.

Necessary and sufficient conditions for optimality in constrained general sum stochastic games

In this paper we first derive a necessary and sufficient condition for a stationary strategy to be the Nash equilibrium of discounted constrained stochastic game under certain assumptions. In this process we also develop a nonlinear (non-convex) optimization problem for a discounted constrained stochastic game. We use the linear best response functions of every player and complementary slackness theorem for linear programs to derive both the optimization problem and the equivalent condition. We then extend this result to average reward constrained stochastic games. Finally, we present a heuristic algorithm motivated by our necessary and sufficient conditions for a discounted cost constrained stochastic game. We numerically observe the convergence of this algorithm to Nash equilibrium. (C) 2015 Elsevier B.V. All rights reserved.

Near-Optimal Detection Thresholds for Bayesian Spectrum Sensing Under Fading

This paper considers the problem of energy-based, Bayesian spectrum sensing in cognitive radios under various fading environments. Under the well-known central limit theorem based model for energy detection, we derive analytically tractable expressions for near-optimal detection thresholds that minimize the probability of error under lognormal, Nakagami-m, and Weibull fading. For the Suzuki fading case, a generalized gamma approximation is provided, which saves on the computation of an integral. In each case, the accuracy of the theoretical expressions as compared to the optimal thresholds are illustrated through simulations.

A note on tetrablock contractions

A commuting triple of operators (A, B, P) on a Hilbert space H is called a tetrablock contraction if the closure of the set E = {(a(11),a(22),detA) : A = GRAPHICS] with parallel to A parallel to <1} is a spectral set. In this paper, we construct a functional model and produce a set of complete unitary invariants for a pure tetrablock contraction. In this construction, the fundamental operators, which are the unique solutions of the operator equations A - B* P = DPX1DP and B - A* P = DPX2DP where X-1, X-2 is an element of B(D-P) play a pivotal role. As a result of the functional model, we show that every pure tetrablock isometry (A, B, P) on an abstract Hilbert space H is unitarily equivalent to the tetrablock contraction (MG1*+G2z, MG2*+G1z, M-z) on H-DP*(2). (D), where G(1) and G(2) are the fundamental operators of (A*, B*, P*). We prove a Beurling Lax Halmos type theorem for a triple of operators (MF1*+F2z, MF2*+F1z, M-z), where epsilon is a Hilbert space and F-1, F-2 is an element of B(epsilon). We also deal with a natural example of tetrablock contraction on a functions space to find out its fundamental operators.

Deviated linear cyclic pursuit

This paper analyses deviated linear cyclic pursuit in which an agent pursues its leader with an angle of deviation in both the continuous- and discrete-time domains, while admitting heterogeneous gains and deviations for the agents. Sufficient conditions for the stability of such systems, in both the domains, are presented in this paper along with the derivation of the reachable set, which is a set of points where the agents may converge asymptotically. The stability conditions are derived based on Gershgorin's theorem. Simulations validating the theoretical results presented in this paper are provided.

Tight triangulations of closed 3-manifolds

A triangulation of a closed 2-manifold is tight with respect to a field of characteristic two if and only if it is neighbourly; and it is tight with respect to a field of odd characteristic if and only if it is neighbourly and orientable. No such characterization of tightness was previously known for higher dimensional manifolds. In this paper, we prove that a triangulation of a closed 3-manifold is tight with respect to a field of odd characteristic if and only if it is neighbourly, orientable and stacked. In consequence, the Kuhnel-Lutz conjecture is valid in dimension three for fields of odd characteristic. Next let F be a field of characteristic two. It is known that, in this case, any neighbourly and stacked triangulation of a closed 3-manifold is F-tight. For closed, triangulated 3-manifolds with at most 71 vertices or with first Betti number at most 188, we show that the converse is true. But the possibility of the existence of an F-tight, non-stacked triangulation on a larger number of vertices remains open. We prove the following upper bound theorem on such triangulations. If an F-tight triangulation of a closed 3-manifold has n vertices and first Betti number beta(1), then (n - 4) (617n - 3861) <= 15444 beta(1). Equality holds here if and only if all the vertex links of the triangulation are connected sums of boundary complexes of icosahedra. (C) 2015 Elsevier Ltd. All rights reserved.

Evaluation of the Effectiveness of Representative Methods for Determining Young's Modulus and Hardness from Instrumented Indentation Data

The effectiveness of Oliver & Pharr's (O&P's) method, Cheng & Cheng's (C&C's) method, and a new method developed by our group for estimating Young's modulus and hardness based on instrumented indentation was evaluated for the case of yield stress to reduced Young's modulus ratio (sigma(y)/E-r) >= 4.55 x 10(-4) and hardening coefficient (n) <= 0.45. Dimensional theorem and finite element simulations were applied to produce reference results for this purpose. Both O&P's and C&C's methods overestimated the Young's modulus under some conditions, whereas the error can be controlled within +/- 16% if the formulation was modified with appropriate correction functions. Similar modification was not introduced to our method for determining Young's modulus, while the maximum error of results was around +/- 13%. The errors of hardness values obtained from all the three methods could be even larger and were irreducible with any correction scheme. It is therefore suggested that when hardness values of different materials are concerned, relative comparison of the data obtained from a single standard measurement technique would be more practically useful. It is noted that the ranges of error derived from the analysis could be different if different ranges of material parameters sigma(y)/E-r and n are considered.

Uniform stability of a particle approximation of the optimal filter derivative

Sequential Monte Carlo methods, also known as particle methods, are a widely used set of computational tools for inference in non-linear non-Gaussian state-space models. In many applications it may be necessary to compute the sensitivity, or derivative, of the optimal filter with respect to the static parameters of the state-space model; for instance, in order to obtain maximum likelihood model parameters of interest, or to compute the optimal controller in an optimal control problem. In Poyiadjis et al. [2011] an original particle algorithm to compute the filter derivative was proposed and it was shown using numerical examples that the particle estimate was numerically stable in the sense that it did not deteriorate over time. In this paper we substantiate this claim with a detailed theoretical study. Lp bounds and a central limit theorem for this particle approximation of the filter derivative are presented. It is further shown that under mixing conditions these Lp bounds and the asymptotic variance characterized by the central limit theorem are uniformly bounded with respect to the time index. We demon- strate the performance predicted by theory with several numerical examples. We also use the particle approximation of the filter derivative to perform online maximum likelihood parameter estimation for a stochastic volatility model.

Hydrodynamic interaction between two vertical cylinders in water waves

The hydrodynamic interaction between two vertical cylinders in water waves is investigated based on the linearized potential flow theory. One of the two cylinders is fixed at the bottom while the other is articulated at the bottom and oscillates with small amplitudes in the direction of the incident wave. Both the diffracted wave and the radiation wave are studied in the present paper. A simple analytical expression for the velocity potential on the surface of each cylinder is obtained by means of Graf's addition theorem. The wave-excited forces and moments on the cylinders, the added masses and the radiation damping coefficients of the oscillating cylinder are all expressed explicitly in series form. The coefficients of the series are determined by solving algebraic equations. Several numerical examples are given to illustrate the effects of various parameters, such as the separation distance, the relative size of the cylinders, and the incident angle, on the first-order and steady second-order forces, the added masses and radiation-damping coefficients as well as the response of the oscillating cylinder.

Dynamical symmetry of screened Coulomb potential and isotropic harmonic oscillator

It is shown that for the screened Coulomb potential and isotropic harmonic oscillator, there exists an infinite number of closed orbits for suitable angular momentum values. At the aphelion (perihelion) points of classical orbits, an extended Runge-Lenz vector for the screened Coulomb potential and an extended quadrupole tensor for the screened isotropic harmonic oscillator are still conserved. For the screened two-dimensional (2D) Coulomb potential and isotropic harmonic oscillator, the dynamical symmetries SO3 and SU(2) are still preserved at the aphelion (perihelion) points of classical orbits, respectively. For the screened 3D Coulomb potential, the dynamical symmetry SO4 is also preserved at the aphelion (perihelion) points of classical orbits. But for the screened 3D isotropic harmonic oscillator, the dynamical symmetry SU(2) is only preserved at the aphelion (perihelion) points of classical orbits in the eigencoordinate system. For the screened Coulomb potential and isotropic harmonic oscillator, only the energy (but not angular momentum) raising and lowering operators can be constructed from a factorization of the radial Schrodinger equation.