
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
priordata 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 ThomasFermi approximation for the singleparticle 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
SuzukiTrotter 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 fourthorder 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 highdimensional 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 sixstate 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 singlephoton experiment for a twopath interferometer with path marking. In
this experiment, the path of a photon ("signal") through a MachZehnder
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 threeoutcome measurement: one
outcome each for the two signal paths plus an inconclusive outcome. Our results
agree fully with the theoretical predictions from a commonsense 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 threepath 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 pathdiscriminating measurement. For a quantummechanical
experimental demonstration of the case, one should perform a singlephoton
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 kineticenergy functional and the
oneparticle density, in particular for twodimensional systems. We derive the
leading gradient corrections from a semiclassical expansion based on Wigner's
phase space formalism and demonstrate that the semiclassical kineticenergy
density functional at zero temperature cannot be evaluated unambiguously. In
contrast, a densitypotential functional description that effectively
incorporates interactions provides unambiguous gradient corrections. Employing
an averaging procedure that involves Airy functions, thereby partially
resumming higherorder 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
Airyaveraged 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 wellbehaved
groundstate kineticenergy functional, which improves on the ThomasFermi
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 . Readymade 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 prechosen size, and are the smallest for the prechosen
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 singlequbit and twoqubit 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 priordata 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 twodimensional (2D) ThomasFermi kineticenergy functional;
it is widely accepted that the 2D analog of the 3D von Weizs\"acker correction
vanishes, together with all higherorder corrections. Contrary to this
longheld 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 singleparticle density.

One can implement fast twoqubit entangling gates by exploiting the Rydberg
blockade. Although various theoretical schemes have been proposed,
experimenters have not yet been able to demonstrate twoatom 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
neutralatom controlledphase 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 twomode coherent state and the quantum
correlation of a twinbeam state. We study the effect of loss and saturation of
the photonnumberresolving detector on the measurement of the VDPs. An
analytic function is derived for this variance, both for the coherent and the
twinbeam 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 nonoptimal. 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 MaassenUffink 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.

Highquality random samples of quantum states are needed for a variety of
tasks in quantum information and quantum computation. Searching the
highdimensional 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 Markovchain 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 twoqubit states probed with a
tomographically incomplete measurement, and then use the sample for the
calculation of the size and credibility of the recentlyintroduced optimal
error regions [see New J. Phys. 15 (2013) 123026]. Another illustration is the
computation of the fractional volume of separable twoqubit states.

Highquality random samples of quantum states are needed for a variety of
tasks in quantum information and quantum computation. Searching the
highdimensional 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 Markovchain 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 twophoton 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 [quantph])
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 timetested 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 linearinversion estimation does not justify the high price of
using unphysical estimators, which are typically the case in that scheme.

Quantum theory is a welldefined 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 twophoton states via
Spontaneous Parametric DownConversion with a twotypeI 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 multilevel systems, the commonly used adiabatic elimination is a method
for approximating the dynamics of the system by eliminating irrelevant,
nonresonantly coupled levels. This procedure is, however, somewhat ambiguous
and it is not clear how to improve on it systematically. We use an
integrodifferential 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
integrodifferential equation is then used to generate a hierarchy of
approximations, in which the zeroth order is the adiabaticelimination
approximation. It works well with a proper choice of interaction picture; the
procedure suggests criteria for optimizing this choice. The firstorder
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 maximumlikelihood point estimator, we
introduce the maximumlikelihood regionthe 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 regionthe smallest
region with prechosen 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 twophoton transitions through an
intermediate state that is nonresonantly 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
adiabaticelimination 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 integrodifferential equation of
LippmannSchwinger 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 maximumlikelihoodmaximumentropy estimator using semidefinite
programming and a standard multidimensional 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
radiofrequency 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 atomion 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 atomion 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 atomion macromolecule. The
methods presented here are not restricted to trapped atomion systems and can
be readily applied to studying the micromotion effect in any system involving a
single trapped ion.