• The multinomial model is one of the simplest statistical models. When constraints are placed on the possible values for the probabilities, however, it becomes much more difficult to deal with. Model checking and checking for prior-data conflict is considered here for such models. A theorem is proved that establishes the consistency of the check on the prior. Applications are presented to models that arise in quantum state estimation as well as the Bayesian analysis of models for ordered probabilities.
  • We improve on the Thomas-Fermi approximation for the single-particle density of fermions by introducing inhomogeneity corrections. Rather than invoking a gradient expansion, we relate the density to the unitary evolution operator for the given effective potential energy and approximate this operator by a Suzuki-Trotter factorization. This yields a hierarchy of approximations, one for each approximate factorization. For the purpose of a first benchmarking, we examine the approximate densities for a few cases with known exact densities and observe a very satisfactory, and encouraging, performance. As a bonus, we also obtain a simple fourth-order leapfrog algorithm for the symplectic integration of classical equations of motion.
  • Quantum key distribution is on the verge of real world applications, where perfectly secure information can be distributed among multiple parties. Several quantum cryptographic protocols have been theoretically proposed and independently realized in different experimental conditions. Here, we develop an experimental platform based on high-dimensional orbital angular momentum states of single photons that enables implementation of multiple quantum key distribution protocols with a single experimental apparatus. Our versatile approach allows us to experimentally survey different classes of quantum key distribution techniques, such as the 1984 Bennett & Brassard (BB84), tomographic protocols including the six-state and the Singapore protocol, and a recently introduced differential phase shift (Chau15) protocol. This enables us to experimentally compare the performance of these techniques and discuss their benefits and deficiencies in terms of noise tolerance in different dimensions. Our analysis gives an overview of the available quantum key distribution protocols for photonic orbital angular momentum and highlights the benefits of the presented schemes for different implementations and channel conditions.
  • When a photon is detected after passing through an interferometer one might wonder which path it took, and a meaningful answer can only be given if one has the means of monitoring the photon's whereabouts. We report the realization of a single-photon experiment for a two-path interferometer with path marking. In this experiment, the path of a photon ("signal") through a Mach--Zehnder interferometer becomes known by unambiguous discrimination between the two paths. We encode the signal path in the polarization state of a partner photon ("idler") whose polarization is examined by a three-outcome measurement: one outcome each for the two signal paths plus an inconclusive outcome. Our results agree fully with the theoretical predictions from a common-sense analysis of what can be said about the past of a quantum particle: The signals for which we get the inconclusive result have full interference strength, as their paths through the interferometer \emph{cannot be known}; and every photon that emerges from the dark output port of the balanced interferometer has a \emph{known path}.
  • We analyze Vaidman's three-path interferometer with weak path marking [Phys. Rev. A 87, 052104 (2013)] and find that common sense yields correct statements about the particle's path through the interferometer. This disagrees with the original claim that the particles have discontinuous trajectories at odds with common sense. In our analysis, "the particle's path" has operational meaning as acquired by a path-discriminating measurement. For a quantum-mechanical experimental demonstration of the case, one should perform a single-photon version of the experiment by Danan et al. [Phys. Rev. Lett. 111, 240402 (2013)] with unambiguous path discrimination. We present a detailed proposal for such an experiment.
  • Building on the discussion in PRA 93, 042510 (2016), we present a systematic derivation of gradient corrections to the kinetic-energy functional and the one-particle density, in particular for two-dimensional systems. We derive the leading gradient corrections from a semiclassical expansion based on Wigner's phase space formalism and demonstrate that the semiclassical kinetic-energy density functional at zero temperature cannot be evaluated unambiguously. In contrast, a density-potential functional description that effectively incorporates interactions provides unambiguous gradient corrections. Employing an averaging procedure that involves Airy functions, thereby partially resumming higher-order gradient corrections, we facilitate a smooth transition of the particle density into the classically forbidden region of arbitrary smooth potentials. We find excellent agreement of the semiclassical Airy-averaged particle densities with the exact densities for very low but finite temperatures, illustrated for a Fermi gas with harmonic potential energy. We furthermore provide criteria for the applicability of the semiclassical expansions at low temperatures. Finally, we derive a well-behaved ground-state kinetic-energy functional, which improves on the Thomas-Fermi approximation.
  • This is the documentation for generating random samples from the quantum state space in accordance with a specified distribution, associated with this webpage: http://tinyurl.com/QSampling . Ready-made samples (each with at least a million points) from various distributions are available for download, or one can generate one's own samples from a chosen distribution using the provided source codes. The sampling relies on the Hamiltonian Monte Carlo algorithm as described in New J. Phys. 17, 043018 (2015). The random samples are reposited in the hope that they would be useful for a variety of tasks in quantum information and quantum computation. Constructing credible regions for tomographic data, optimizing a function over the quantum state space with a complicated landscape, testing the typicality of entanglement among states from a multipartite quantum system, or computing the average of some quantity of interest over a subset of quantum states are but some exemplary applications among many.
  • Quantum state estimation aims at determining the quantum state from observed data. Estimating the full state can require considerable efforts, but one is often only interested in a few properties of the state, such as the fidelity with a target state, or the degree of correlation for a specified bipartite structure. Rather than first estimating the state, one can, and should, estimate those quantities of interest directly from the data. We propose the use of optimal error intervals as a meaningful way of stating the accuracy of the estimated property values. Optimal error intervals are analogs of the optimal error regions for state estimation [New J. Phys. 15, 123026 (2013)]. They are optimal in two ways: They have the largest likelihood for the observed data and the pre-chosen size, and are the smallest for the pre-chosen probability of containing the true value. As in the state situation, such optimal error intervals admit a simple description in terms of the marginal likelihood for the data for the properties of interest. Here, we present the concept and construction of optimal error intervals, report on an iterative algorithm for reliable computation of the marginal likelihood (a quantity difficult to calculate reliably), explain how plausible intervals --- a notion of evidence provided by the data --- are related to our optimal error intervals, and illustrate our methods with single-qubit and two-qubit examples.
  • When using complex Bayesian models to combine information, the checking for consistency of the information being combined is good statistical practice. Here a new method is developed for detecting prior-data conflicts in Bayesian models based on comparing the observed value of a prior to posterior divergence to its distribution under the prior predictive distribution for the data. The divergence measure used in our model check is a measure of how much beliefs have changed from prior to posterior, and can be thought of as a measure of the overall size of a relative belief function. It is shown that the proposed method is intuitive, has desirable properties, can be extended to hierarchical settings, and is related asymptotically to Jeffreys' and reference prior distributions. In the case where calculations are difficult, the use of variational approximations as a way of relieving the computational burden is suggested. The methods are compared in a number of examples with an alternative but closely related approach in the literature based on the prior predictive distribution of a minimal sufficient statistic.
  • Density functional theory (DFT) is notorious for the absence of gradient corrections to the two-dimensional (2D) Thomas-Fermi kinetic-energy functional; it is widely accepted that the 2D analog of the 3D von Weizs\"acker correction vanishes, together with all higher-order corrections. Contrary to this long-held belief, we show that the leading correction to the kinetic energy does not vanish, is unambiguous, and contributes perturbatively to the total energy. This insight emerges naturally in a simple extension of standard DFT, which has the effective potential energy as a functional variable on equal footing with the single-particle density.
  • One can implement fast two-qubit entangling gates by exploiting the Rydberg blockade. Although various theoretical schemes have been proposed, experimenters have not yet been able to demonstrate two-atom gates of high fidelity due to experimental constraints. We propose a novel scheme, which only uses a single Rydberg pulse illuminating both atoms, for the construction of neutral-atom controlled-phase gates. In contrast to the existing schemes, our approach is simpler to implement and requires neither individual addressing of atoms nor adiabatic procedures. With parameters estimated based on actual experimental scenarios, a gate fidelity higher than 0.99 is achievable.
  • The variance of difference of photocounts (VDPs) is an established measure of quantum correlations for quantum states of light. It enables us to discriminate between the classical correlation of a two-mode coherent state and the quantum correlation of a twin-beam state. We study the effect of loss and saturation of the photon-number-resolving detector on the measurement of the VDPs. An analytic function is derived for this variance, both for the coherent and the twin-beam states. It is found that the VDPs is no longer a reliable entanglement measure in the nonlinear regime of the detector response but it remains useful in some range of values of average photon numbers of the incident light. We also quantify the linear regime of the detector with saturation which will be useful for calibration of the detector quantum efficiency.
  • The entropic uncertainty relation proven by Maassen and Uffink for arbitrary pairs of two observables is known to be non-optimal. Here, we call an uncertainty relation optimal, if the lower bound can be attained for any value of either of the corresponding uncertainties. In this work we establish optimal uncertainty relations by characterising the optimal lower bound in scenarios similar to the Maassen-Uffink type. We disprove a conjecture by Englert et al. and generalise various previous results. However, we are still far from a complete understanding and, based on numerical investigation and analytical results in small dimension, we present a number of conjectures.
  • High-quality random samples of quantum states are needed for a variety of tasks in quantum information and quantum computation. Searching the high-dimensional quantum state space for a global maximum of an objective function with many local maxima or evaluating an integral over a region in the quantum state space are but two exemplary applications of many. These tasks can only be performed reliably and efficiently with Monte Carlo methods, which involve good samplings of the parameter space in accordance with the relevant target distribution. We show how the standard strategies of rejection sampling, importance sampling, and Markov-chain sampling can be adapted to this context, where the samples must obey the constraints imposed by the positivity of the statistical operator. For a comparison of these sampling methods, we generate sample points in the probability space for two-qubit states probed with a tomographically incomplete measurement, and then use the sample for the calculation of the size and credibility of the recently-introduced optimal error regions [see New J. Phys. 15 (2013) 123026]. Another illustration is the computation of the fractional volume of separable two-qubit states.
  • High-quality random samples of quantum states are needed for a variety of tasks in quantum information and quantum computation. Searching the high-dimensional quantum state space for a global maximum of an objective function with many local maxima or evaluating an integral over a region in the quantum state space are but two exemplary applications of many. These tasks can only be performed reliably and efficiently with Monte Carlo methods, which involve good samplings of the parameter space in accordance with the relevant target distribution. We show how the Markov-chain Monte Carlo method known as Hamiltonian Monte Carlo, or hybrid Monte Carlo, can be adapted to this context. It is applicable when an efficient parameterization of the state space is available. The resulting random walk is entirely inside the physical parameter space, and the Hamiltonian dynamics enable us to take big steps, thereby avoiding strong correlations between successive sample points while enjoying a high acceptance rate. We use examples of single and double qubit measurements for illustration.
  • We report an experiment in which one determines, with least tomographic effort, whether an unknown two-photon polarization state is entangled or separable. The method measures whole families of optimal entanglement witnesses. We introduce adaptive measurement schemes that greatly speed up the entanglement detection. The experiments are performed on states of different ranks, and we find good agreement with results from computer simulations.
  • Because of the constraint that the estimators be bona fide physical states, any quantum state tomography scheme - including the widely used maximum likelihood estimation - yields estimators that may have a bias, although they are consistent estimators. Schwemmer et al. (arXiv:1310.8465 [quant-ph]) illustrate this by observing a systematic underestimation of the fidelity and an overestimation of entanglement in estimators obtained from simulated data. Further, these authors argue that the simple method of linear inversion overcomes this (perceived) problem of bias, and there is the suggestion to abandon time-tested estimation procedures in favor of linear inversion. Here, we discuss the pros and cons of using biased and unbiased estimators for quantum state tomography. We conclude that the little occasional benefit from the unbiased linear-inversion estimation does not justify the high price of using unphysical estimators, which are typically the case in that scheme.
  • Quantum theory is a well-defined local theory with a clear interpretation. No "measurement problem" or any other foundational matters are waiting to be settled.
  • Disordered systems have grown in importance in the past decades, with similar phenomena manifesting themselves in many different physical systems. Because of the difficulty of the topic, theoretical progress has mostly emerged from numerical studies or analytical approximations. Here, we provide an exact, analytical solution to the problem of uniform phase disorder in a system of identical scatterers arranged with varying separations along a line. Relying on a relationship with Legendre functions, we demonstrate a simple approach to computing statistics of the transmission probability (or the conductance, in the language of electronic transport), and its reciprocal (or the resistance). Our formalism also gives the probability distribution of the conductance, which reveals features missing from previous approaches to the problem.
  • We report a controllable method for producing mixed two-photon states via Spontaneous Parametric Down-Conversion with a two-type-I crystal geometry. By using variable polarization rotators (VPRs), one obtains mixed states of various purities and degrees of entanglement depending on the parameters of the VPRs. The generated states are characterized by quantum state tomography. The experimental results are found to be in good agreement with the theory. The method can be easily implemented for various experiments which require the generation of states with controllable degrees of entanglement or mixedness.
  • In multi-level systems, the commonly used adiabatic elimination is a method for approximating the dynamics of the system by eliminating irrelevant, non-resonantly coupled levels. This procedure is, however, somewhat ambiguous and it is not clear how to improve on it systematically. We use an integro-differential equation for the probability amplitudes of the levels of interest, which is equivalent to the original Schrodinger equation for all probability amplitudes. In conjunction with a Markov approximation, the integro-differential equation is then used to generate a hierarchy of approximations, in which the zeroth order is the adiabatic-elimination approximation. It works well with a proper choice of interaction picture; the procedure suggests criteria for optimizing this choice. The first-order approximation in the hierarchy provides significant improvements over standard adiabatic elimination, without much increase in complexity, and is furthermore not so sensitive to the choice of interaction picture. We illustrate these points with several examples.
  • Rather than point estimators, states of a quantum system that represent one's best guess for the given data, we consider optimal regions of estimators. As the natural counterpart of the popular maximum-likelihood point estimator, we introduce the maximum-likelihood region---the region of largest likelihood among all regions of the same size. Here, the size of a region is its prior probability. Another concept is the smallest credible region---the smallest region with pre-chosen posterior probability. For both optimization problems, the optimal region has constant likelihood on its boundary. We discuss criteria for assigning prior probabilities to regions, and illustrate the concepts and methods with several examples.
  • Driven Raman processes --- nearly resonant two-photon transitions through an intermediate state that is non-resonantly coupled and does not acquire a sizeable population --- are commonly treated with a simplified description in which the intermediate state is removed by adiabatic elimination. While the adiabatic-elimination approximation is reliable when the detuning of the intermediate state is quite large, it cannot be trusted in other situations, and it does not allow one to estimate the population in the eliminated state. We introduce an alternative method that keeps all states in the description, without increasing the complexity by much. An integro-differential equation of Lippmann-Schwinger type generates a hierarchy of approximations, but very accurate results are already obtained in the lowest order.
  • There exists, in general, a convex set of quantum state estimators that maximize the likelihood for informationally incomplete data. We propose an estimation scheme, catered to measurement data of this kind, to search for the exact maximum-likelihood-maximum-entropy estimator using semidefinite programming and a standard multi-dimensional function optimization routine. This scheme can be used to infer the expectation values of a set of entanglement witnesses that can be used to verify the entanglement of the unknown quantum state for composite systems. Next, we establish an alternative numerical scheme that is more computationally robust for the sole purpose of maximizing the likelihood and entropy.
  • We examine the validity of the harmonic approximation, where the radio-frequency ion trap is treated as a harmonic trap, in the problem regarding the controlled collision of a trapped atom and a single trapped ion. This is equivalent to studying the effect of the micromotion since this motion must be neglected for the trapped ion to be considered as a harmonic oscillator. By applying the transformation of Cook and Shankland we find that the micromotion can be represented by two periodically oscillating operators. In order to investigate the effect of the micromotion on the dynamics of a trapped atom-ion system, we calculate (i) the coupling strengths of the micromotion operators by numerical integration and (ii) the quasienergies of the system by applying the Floquet formalism, a useful framework for studying periodic systems. It turns out that the micromotion is not negligible when the distance between the atom and the ion traps is shorter than a characteristic distance. Within this range the energy diagram of the system changes remarkably when the micromotion is taken into account, which leads to undesirable consequences for applications that are based on an adiabatic process of the trapped atom-ion system. We suggest a simple scheme for bypassing the micromotion effect in order to successfully implement a quantum controlled phase gate proposed previously and create an atom-ion macromolecule. The methods presented here are not restricted to trapped atom-ion systems and can be readily applied to studying the micromotion effect in any system involving a single trapped ion.