
We study an infinitehorizon discretetime optimal stopping problem under
nonexponential discounting. A new method, which we call the iterative
approach, is developed to find subgame perfect Nash equilibria. When the
discount function induces decreasing impatience, we establish the existence of
an equilibrium through fixedpoint iterations. Moreover, we show that there
exists a unique optimal equilibrium, which generates larger value than any
other equilibrium does at all times. To the best of our knowledge, this is the
first time a dominating subgame perfect Nash equilibrium is shown to exist in
the literature of timeinconsistency.

In this paper, we consider the estimation and inference of precision matrices
of a rich class of locally stationary and nonlinear time series assuming that
only one realization of the time series is observed. Using a Cholesky
decomposition technique, we show that the precision matrices can be directly
estimated via a series of least squares linear regressions with smoothly
timevarying coefficients. The method of sieves is utilized for the estimation
and is shown to be efficient and optimally adaptive in terms of estimation
accuracy and computational complexity. We establish an asymptotic theory for a
class of ${\cal L}^2$ tests based on the nonparametric sieve estimators. The
latter are used for testing whether the precision matrices are diagonal or
banded. A high dimensional Gaussian approximation result is established for a
wide class of quadratic form of nonstationary and nonlinear processes, which
is of interest by itself.

Let $\Omega$ to be one of $\X^{N+1},C[0,1],D[0,1]$: a product of Polish
spaces, space of continuous functions from $[0,1]$ to a subset of
$\mathbb{R}^d$, and space of RCLL (rightcontinuous with left limits) functions
from $[0,1]$ to $\mathbb{R}^d$ respectively. We consider the existence of a
probability measure $P$ on $\Omega$ such that $P$ has the given marginals
$\alpha$ and $\beta$ and satisfies some other convex transport constraints,
which is given by $\Gamma$. The main application we have in mind is the
martingale optimal transport problem with when the martingales are assumed to
have bounded volatility/quadratic variation. We show that such probability
measure exists if and only if the $\alpha$ average of socalled $G$expectation
of bounded uniformly continuous and bounded functions with respect to the
measures in $\Gamma$ is less than their $\beta$ average.
As a byproduct, we get a necessary and sufficient condition for the Skorokhod
embedding for bounded stopping times.

Inspired by Strotz's consistent planning strategy, we formulate the infinite
horizon meanvariance stopping problem as a subgame perfect Nash equilibrium in
order to determine time consistent strategies with no regret. Equilibria among
stopping times or randomized stopping times may not exist. This motivates us to
consider the notion of liquidation strategies, which lets the stopping right to
be divisible. We then argue that the meanstandard deviation variant of this
problem makes more sense for this type of strategies in terms of time
consistency. It turns out that an equilibrium liquidation strategy always
exists. We then analyze whether optimal equilibrium liquidation strategies
exist and whether they are unique and observe that neither may hold.

A restrictive assumption in change point analysis is "stationarity under the
null hypothesis of no changepoint", which is crucial for asymptotic theory but
not very realistic from a practical point of view. For example, if change point
analysis for correlations is performed, it is not necessarily clear that the
mean, marginal variance or higher order moments are constant, even if there is
no change in the correlation. This paper develops change point analysis for the
correlation structures under less restrictive assumptions. In contrast to
previous work, our approach does not require that the mean, variance and fourth
order joint cumulants are constant under the null hypothesis. Moreover, we also
address the problem of detecting relevant change points.

Since most of the traded options on individual stocks is of American type it
is of interest to generalize the results obtained in semistatic trading to the
case when one is allowed to statically trade American options. However, this
problem has proved to be elusive so far because of the asymmetric nature of the
positions of holding versus shorting such options. Here we provide a unified
framework and generalize the fundamental theorem of asset pricing (FTAP) and
hedging dualities in arXiv:1502.06681 (to appear in Annals of Applied
Probability) to the case where the investor can also short American options.
Following arXiv:1502.06681, we assume that the longed American options are
divisible. As for the shorted American options, we show that the divisibility
plays no role regarding arbitrage property and hedging prices. Then using the
method of enlarging probability spaces proposed in arXiv:1604.05517, we convert
the shorted American options to European options, and establish the FTAP and
sub and superhedging dualities in the enlarged space both with and without
model uncertainty.

A neutron imaging detector based on neutronsensitive microchannel plates
(MCPs) was constructed and tested at beamlines of thermal and cold neutrons.
The MCPs are made of a glass mixture containing B10 and natural Gd, which
makes the bulk of the MCP an efficient neutron converter. Contrary to the
neutron sensitive scintillator screens normally used in neutron imaging,
spatial resolution is not traded off with detection efficiency. While the best
neutron imaging scintillators have a detection efficiency around a percent, a
detection efficiency of around 50% for thermal neutrons and 70% for cold
neutrons has been demonstrated with these MCPs earlier.
Our tests show a performance similar to conventional neutron imaging
detectors, apart from the orders of magnitude better sensitivity. We
demonstrate a spatial resolution better than 150 um. The sensitivity of this
detector allows fast tomography and neutron video recording, and will make
smaller reactor sites and even portable sources suitable for neutron imaging.

In this paper, we introduce a new sparsitypromoting prior, namely, the
"normal product" prior, and develop an efficient algorithm for sparse signal
recovery under the Bayesian framework. The normal product distribution is the
distribution of a product of two normally distributed variables with zero means
and possibly different variances. Like other sparsityencouraging distributions
such as the Student's $t$distribution, the normal product distribution has a
sharp peak at origin, which makes it a suitable prior to encourage sparse
solutions. A twostage normal productbased hierarchical model is proposed. We
resort to the variational Bayesian (VB) method to perform the inference.
Simulations are conducted to illustrate the effectiveness of our proposed
algorithm as compared with other stateoftheart compressed sensing
algorithms.

We consider the superhedging price of an American option in a discretetime
market in which stocks are available for dynamic trading and European options
are available for static trading. We show that the superhedging price $\pi$ is
given by the supremum over the prices of the American option under randomized
models. That is, $\pi=\sup_{(c_i,Q_i)_i}\sum_ic_i\phi^{Q_i}$, where $c_i \in
\mathbb{R}_+$ and the martingale measure $Q^i$ are chosen such that $\sum_i
c_i=1$ and $\sum_i c_iQ_i$ prices the European options correctly, and
$\phi^{Q_i}$ is the price of the American option under the model $Q_i$. Our
result generalizes the example given in ArXiv:1604.02274 that the highest model
based price can be considered as a randomization over models.

On a filtered probability space
$(\Omega,\mathcal{F},P,\mathbb{F}=(\mathcal{F}_t)_{t=0,\dotso,T})$, we consider
stopperstopper games $\overline
V:=\inf_{\Rho\in\bT^{ii}}\sup_{\tau\in\T}\E[U(\Rho(\tau),\tau)]$ and
$\underline V:=\sup_{\Tau\in\bT^i}\inf_{\rho\in\T}\E[U(\Rho(\tau),\tau)]$ in
discrete time, where $U(s,t)$ is $\mathcal{F}_{s\vee t}$measurable instead of
$\mathcal{F}_{s\wedge t}$measurable as is often assumed in the literature,
$\T$ is the set of stopping times, and $\bT^i$ and $\bT^{ii}$ are sets of
mappings from $\T$ to $\T$ satisfying certain nonanticipativity conditions. We
convert the problems into a corresponding Dynkin game, and show that $\overline
V=\underline V=V$, where $V$ is the value of the Dynkin game. We also get the
optimal $\Rho\in\bT^{ii}$ and $\Tau\in\bT^i$ for $\overline V$ and $\underline
V$ respectively.

We consider the problem of downlink channel estimation for millimeter wave
(mmWave) MIMOOFDM systems, where both the base station (BS) and the mobile
station (MS) employ large antenna arrays for directional precoding/beamforming.
Hybrid analog and digital beamforming structures are employed in order to offer
a compromise between hardware complexity and system performance. Different from
most existing studies that are concerned with narrowband channels, we consider
estimation of wideband mmWave channels with frequency selectivity, which is
more appropriate for mmWave MIMOOFDM systems. By exploiting the sparse
scattering nature of mmWave channels, we propose a CANDECOMP/PARAFAC (CP)
decompositionbased method for channel parameter estimation (including angles
of arrival/departure, time delays, and fading coefficients). In our proposed
method, the received signal at the BS is expressed as a thirdorder tensor. We
show that the tensor has the form of a lowrank CP decomposition, and the
channel parameters can be estimated from the associated factor matrices. Our
analysis reveals that the uniqueness of the CP decomposition can be guaranteed
even when the size of the tensor is small. Hence the proposed method has the
potential to achieve substantial training overhead reduction. We also develop
CramerRao bound (CRB) results for channel parameters, and compare our proposed
method with a compressed sensingbased method. Simulation results show that the
proposed method attains mean square errors that are very close to their
associated CRBs, and presents a clear advantage over the compressed
sensingbased method in terms of both estimation accuracy and computational
complexity.

We consider the problem of uplink channel estimation for millimeter wave
(mmWave) systems, where the base station (BS) and mobile stations (MSs) are
equipped with large antenna arrays to provide sufficient beamforming gain for
outdoor wireless communications. Hybrid analog and digital beamforming
structures are employed by both the BS and the MS due to hardware constraints.
We propose a layered pilot transmission scheme and a CANDECOMP/PARAFAC (CP)
decompositionbased method for joint estimation of the channels from multiple
users (i.e. MSs) to the BS. The proposed method exploits the sparse scattering
nature of the mmWave channel and the intrinsic multidimensional structure of
the multiway data collected from multiple modes. The uniqueness of the CP
decomposition is studied and sufficient conditions for essential uniqueness are
obtained. The conditions shed light on the design of the beamforming matrix,
the combining matrix and the pilot sequences, and meanwhile provide general
guidelines for choosing system parameters. Our analysis reveals that our
proposed method can achieve a substantial training overhead reduction by
employing the layered pilot transmission scheme. Simulation results show that
the proposed method presents a clear advantage over a compressed sensingbased
method in terms of both estimation accuracy and computational complexity.

We consider a financial market where stocks are available for dynamic
trading, and European and American options are available for static trading
(semistatic trading strategies). We assume that the American options are
infinitely divisible, and can only be bought but not sold. In the first part of
the paper, we work within the framework without model ambiguity. We first get
the fundamental theorem of asset pricing (FTAP). Using the FTAP, we get the
dualities for the hedging prices of European and American options. Based on the
hedging dualities, we also get the duality for the utility maximization. In the
second part of the paper, we consider the market which admits nondominated
model uncertainty. We first establish the hedging result, and then using the
hedging duality we further get the FTAP. Due to the technical difficulty
stemming from the nondominancy of the probability measure set, we use a
discretization technique and apply the minimax theorem.

We present a comprehensive structural characterization of two different
highly pure nuclear graphites that compasses all relevant length scales from
nanometers to submm. This has been achieved by combining several experiments
and neutron techniques: Small Angle Neutron Scattering (SANS), highresolution
Spin Echo SANS (SESANS) and neutron imaging. In this way it is possible to
probe an extraordinary broad range of 6 orders of magnitude in length from
microscopic to macroscopic length scales. The results reveal a fractal
structure that extends from $ \sim $ 0.6 nm to 0.6 mm and has surface and mass
fractal dimensions both very close to 2.5, a value found for percolating
clusters and fractured ranked surfaces in 3D.

We consider multiplayer stopping games in continuous time. Unlike Dynkin
games, in our games the payoff of each player is revealed after all the players
stop. Moreover, each player can adjust her own stopping strategy by observing
other players' behaviors. Assuming the continuity of the payoff functions in
time, we show that there always exists an $\eps$Nash equilibrium in pure
stopping strategies for any $\eps>0$.

We consider twoplayer nonzerosum stopping games in discrete time. Unlike
Dynkin games, in our games the payoff of each player is revealed after both
players stop. Moreover, each player can adjust her own stopping strategy
according to the other player's action. In the first part of the paper, we
consider the game where players act simultaneously at each stage. We show that
there exists a Nash equilibrium in mixed stopping strategies. In the second
part, we assume that one player has to act first at each stage. In this case,
we show the existence of a Nash equilibrium in pure stopping strategies.

On a filtered probability space $(\Omega ,\mathcal{F},
(\mathcal{F}_t)_{t\in[0,\infty]}, \mathbb{P})$, we consider the twoplayer
nonzerosum stopping game $u^i := \mathbb{E}[U^i(\rho,\tau)],\ i=1,2$, where
the first player choose a stopping strategy $\rho$ to maximize $u^1$ and the
second player chose a stopping strategy $\tau$ to maximize $u^2$. Unlike the
Dynkin game, here we assume that $U(s,t)$ is $\mathcal{F}_{s\vee
t}$measurable. Assuming the continuity of $U^i$ in $(s,t)$, we show that there
exists an $\epsilon$Nash equilibrium for any $\epsilon>0$.

We consider a zerosum continuous time stopping game in which the payoff is
revealed in the maximum of the two stopping times instead of the minimum, which
is the case in Dynkin games.

We consider as given a discrete time financial market with a risky asset and
options written on that asset and determine both the sub and superhedging
prices of an American option in the model independent framework of
ArXiv:1305.6008. We obtain the duality of results for the sub and
superhedging prices. For the subhedging prices we discuss whether the sup and
inf in the dual representation can be exchanged (a counter example shows that
this is not true in general). For the superhedging prices we discuss several
alternative definitions and argue why our choice is more reasonable. Then
assuming that the path space is compact, we construct a discretization of the
path space and demonstrate the convergence of the hedging prices at the optimal
rate. The latter result would be useful for numerical computation of the
hedging prices. Our results generalize those of ArXiv:1304.3574 to the case
when static positions in (finitely many) European options can be used in the
hedging portfolio.

We consider the optimal stopping problem
$v^{(\eps)}:=\sup_{\tau\in\mathcal{T}_{0,T}}\mathbb{E}B_{(\tau\eps)^+}$ posed
by Shiryaev at the International Conference on Advanced Stochastic Optimization
Problems organized by the Steklov Institute of Mathematics in September 2012.
Here $T>0$ is a fixed time horizon, $(B_t)_{0\leq t\leq T}$ is the Brownian
motion, $\eps\in[0,T]$ is a constant, and $\mathcal{T}_{\eps,T}$ is the set of
stopping times taking values in $[\eps,T]$. The solution of this problem is
characterized by a path dependent reflected backward stochastic differential
equations, from which the continuity of $\eps \to v^{(\eps)}$ follows. For
large enough $\eps$, we obtain an explicit expression for $v^{(\eps)}$ and for
small $\eps$ we have lower and upper bounds. The main result of the paper is
the asymptotics of $v^{(\eps)}$ as $\eps\searrow 0$. As a byproduct, we also
obtain L\'{e}vy's modulus of continuity result in the $L^1$ sense.

An important assumption in the work on testing for structural breaks in time
series consists in the fact that the model is formulated such that the
stochastic process under the null hypothesis of "no changepoint" is
stationary. This assumption is crucial to derive (asymptotic) critical values
for the corresponding testing procedures using an elegant and powerful
mathematical theory, but it might be not very realistic from a practical point
of view.
This paper develops change point analysis under less restrictive assumptions
and deals with the problem of detecting change points in the marginal variance
and correlation structures of a nonstationary time series. A CUSUM approach is
proposed, which is used to test the "classical" hypothesis of the form $H_0:
\theta_1=\theta_2$ vs. $H_1: \theta_1 \not =\theta_2$, where $\theta_1$ and
$\theta_2$ denote second order parameters of the process before and after a
change point. The asymptotic distribution of the CUSUM test statistic is
derived under the null hypothesis. This distribution depends in a complicated
way on the dependency structure of the nonlinear nonstationary time series and
a bootstrap approach is developed to generate critical values. The results are
then extended to test the hypothesis of a {\it non relevant change point}, i.e.
$H_0:  \theta_1\theta_2  \leq \delta$, which reflects the fact that
inference should not be changed, if the difference between the parameters
before and after the changepoint is small.
In contrast to previous work, our approach does neither require the mean to
be constant nor  in the case of testing for lag $k$correlation  that the
mean, variance and fourth order joint cumulants are constant under the null
hypothesis. In particular, we allow that the variance has a change point at a
different location than the autocovariance.

We consider the fundamental theorem of asset pricing (FTAP) and hedging
prices of options under nondominated model uncertainty and portfolio
constrains in discrete time. We first show that no arbitrage holds if and only
if there exists some family of probability measures such that any admissible
portfolio value process is a local supermartingale under these measures. We
also get the nondominated optional decomposition with constraints. From this
decomposition, we get duality of the superhedging prices of European options,
as well as the sub and superhedging prices of American options. Finally, we
get the FTAP and duality of superhedging prices in a market where stocks are
traded dynamically and options are traded statically.

We show that the results of ArXiv:1305.6008 on the Fundamental Theorem of
Asset Pricing and the superhedging theorem can be extended to the case in
which the options available for static hedging (\emph{hedging options}) are
quoted with bidask spreads. In this setup, we need to work with the notion of
\emph{robust noarbitrage} which turns out to be equivalent to noarbitrage
under the additional assumption that hedging options with nonzero spread are
\emph{nonredundant}. A key result is the closedness of the set of attainable
claims, which requires a new proof in our setting.

We investigate the behavior of the Generalized Likelihood Ratio Test (GLRT)
(Fan, Zhang and Zhang [Ann. Statist. 29 (2001) 153193]) for time varying
coefficient models where the regressors and errors are nonstationary time
series and can be cross correlated. It is found that the GLRT retains the
minimax rate of local alternative detection under weak dependence and
nonstationarity. However, in general, the Wilks phenomenon as well as the
classic residual bootstrap are sensitive to either conditional
heteroscedasticity of the errors, nonstationarity or temporal dependence. An
averaged test is suggested to alleviate the sensitivity of the test to the
choice of bandwidth and is shown to be more powerful than tests based on a
single bandwidth. An alternative wild bootstrap method is proposed and shown to
be consistent when making inference of time varying coefficient models for
nonstationary time series.

We investigate the behavior of Fourier transforms for a wide class of
nonstationary nonlinear processes. Asymptotic central and noncentral limit
theorems are established for a class of nondegenerate and degenerate weighted
$V$statistics through the angle of Fourier analysis. The established theory
for $V$statistics provides a unified treatment for many important time and
spectral domain problems in the analysis of nonstationary time series, ranging
from nonparametric estimation to the inference of periodograms and spectral
densities.