Goal-based h-adaptivity of the 1-D diamond difference discrete ordinate method.

The quantity of interest (QoI) associated with a solution of a partial differential equation (PDE) is not, in general, the solution itself, but a functional of the solution. Dual weighted residual (DWR) error estimators are one way of providing an estimate of the error in the QoI resulting from the discretisation of the PDE. This paper aims to provide an estimate of the error in the QoI due to the spatial discretisation, where the discretisation scheme being used is the diamond difference (DD) method in space and discrete ordinate (SNSN) method in angle. The QoI are reaction rates in detectors and the value of the eigenvalue (Keff)(Keff) for 1-D fixed source and eigenvalue (KeffKeff criticality) neutron transport problems respectively. Local values of the DWR over individual cells are used as error indicators for goal-based mesh refinement, which aims to give an optimal mesh for a given QoI

General CMB and Primordial Bispectrum Estimation I: Mode Expansion, Map-Making and Measures of f_NL

We present a detailed implementation of two bispectrum estimation methods which can be applied to general non-separable primordial and CMB bispectra. The method exploits bispectrum mode decompositions on the domain of allowed wavenumber or multipole values. Concrete mode examples constructed from symmetrised tetrahedral polynomials are given, demonstrating rapid convergence for known bispectra. We use these modes to generate simulated CMB maps of high resolution (l > 2000) given an arbitrary primordial power spectrum and bispectrum or an arbitrary late-time CMB angular power spectrum and bispectrum. By extracting coefficients for the same separable basis functions from an observational map, we are able to present an efficient and general f_NL estimator for a given theoretical model. The estimator has two versions comparing theoretical and observed coefficients at either primordial or late times, thus encompassing a wider range of models, including secondary anisotropies, lensing and cosmic strings. We provide examples and validation of both f_NL estimation methods by direct comparison with simulations in a WMAP-realistic context. In addition, we show how the full bispectrum can be extracted from observational maps using these mode expansions, irrespective of the theoretical model under study. We also propose a universal definition of the bispectrum parameter F_NL for more consistent comparison between theoretical models. We obtain WMAP5 estimates of f_NL for the equilateral model from both our primordial and late-time estimators which are consistent with each other, as well as with results already published in the literature. These general bispectrum estimation methods should prove useful for the analysis of nonGaussianity in the Planck satellite data, as well as in other contexts.Comment: 41 pages, 17 figure

A time-splitting pseudospectral method for the solution of the Gross-Pitaevskii equations using spherical harmonics with generalised-Laguerre basis functions

We present a method for numerically solving a Gross-Pitaevskii system of equations with a harmonic and a toroidal external potential that governs the dynamics of one- and two-component Bose-Einstein condensates. The method we develop maintains spectral accuracy by employing Fourier or spherical harmonics in the angular coordinates combined with generalised-Laguerre basis functions in the radial direction. Using an error analysis, we show that the method presented leads to more accurate results than one based on a sine transform in the radial direction when combined with a time-splitting method for integrating the equations forward in time. In contrast to a number of previous studies, no assumptions of radial or cylindrical symmetry is assumed allowing the method to be applied to 2D and 3D time-dependent simulations. This is accomplished by developing an efficient algorithm that accurately performs the generalised-Laguerre transforms of rotating Bose-Einstein condensates for different orders of the Laguerre polynomials. Using this spatial discretisation together with a second order Strang time-splitting method, we illustrate the scheme on a number of 2D and 3D computations of the ground state of a non-rotating and rotating condensate. Comparisons between previously derived theoretical results for these ground state solutions and our numerical computations show excellent agreement for these benchmark problems. The method is further applied to simulate a number of time-dependent problems including the Kelvin-Helmholtz instability in a two-component rotating condensate and the motion of quantised vortices in a 3D condensate

The Radial Masa in a Free Group Factor is Maximal Injective

The radial (or Laplacian) masa in a free group factor is the abelian von Neumann algebra generated by the sum of the generators (of the free group) and their inverses. The main result of this paper is that the radial masa is a maximal injective von Neumann subalgebra of a free group factor. We also investigate tensor products of maximal injective algebras. Given two inclusions $B_i\subset M_i$ of type $\mathrm{I}$ von Neumann algebras in finite von Neumann algebras such that each $B_i$ is maximal injective in $M_i$, we show that the tensor product $B_1\ \bar{\otimes}\ B_2$ is maximal injective in $M_1\ \bar{\otimes}\ M_2$ provided at least one of the inclusions satisfies the asymptotic orthogonality property we establish for the radial masa. In particular it follows that finite tensor products of generator and radial masas will be maximal injective in the corresponding tensor product of free group factors.Comment: 25 Pages, Typos corrected and exposition improve

Discrete Quantum Mechanics

A comprehensive review of the discrete quantum mechanics with the pure imaginary shifts and the real shifts is presented in parallel with the corresponding results in the ordinary quantum mechanics. The main subjects to be covered are the factorised Hamiltonians, the general structure of the solution spaces of the Schroedinger equation (Crum's theorem and its modification), the shape invariance, the exact solvability in the Schroedinger picture as well as in the Heisenberg picture, the creation/annihilation operators and the dynamical symmetry algebras, the unified theory of exact and quasi-exact solvability based on the sinusoidal coordinates, the infinite families of new orthogonal (the exceptional) polynomials. Two new infinite families of orthogonal polynomials, the X_\ell Meixner-Pollaczek and the X_\ell Meixner polynomials are reported.Comment: 61 pages, 1 figure. Comments and references adde

Quantum states cannot be transmitted efficiently classically

We show that any classical two-way communication protocol with shared randomness that can approximately simulate the result of applying an arbitrary measurement (held by one party) to a quantum state of $n$ qubits (held by another), up to constant accuracy, must transmit at least $\Omega(2^n)$ bits. This lower bound is optimal and matches the complexity of a simple protocol based on discretisation using an $\epsilon$-net. The proof is based on a lower bound on the classical communication complexity of a distributed variant of the Fourier sampling problem. We obtain two optimal quantum-classical separations as easy corollaries. First, a sampling problem which can be solved with one quantum query to the input, but which requires $\Omega(N)$ classical queries for an input of size $N$. Second, a nonlocal task which can be solved using $n$ Bell pairs, but for which any approximate classical solution must communicate $\Omega(2^n)$ bits.Comment: 24 pages; v3: accepted version incorporating many minor corrections and clarification