
This paper studies the dimension effect of the linear discriminant analysis
(LDA) and the regularized linear discriminant analysis (RLDA) classifiers for
large dimensional data where the observation dimension $p$ is of the same order
as the sample size $n$. More specifically, built on properties of the Wishart
distribution and recent results in random matrix theory, we derive explicit
expressions for the asymptotic misclassification errors of LDA and RLDA
respectively, from which we gain insights of how dimension affects the
performance of classification and in what sense. Motivated by these results, we
propose adjusted classifiers by correcting the bias brought by the unequal
sample sizes. The biascorrected LDA and RLDA classifiers are shown to have
smaller misclassification rates than LDA and RLDA respectively. Several
interesting examples are discussed in detail and the theoretical results on
dimension effect are illustrated via extensive simulation studies.

In this paper, we focus on numerical solutions for random genetic drift
problem, which is governed by a degenerated convectiondominated parabolic
equation. Due to the fixation phenomenon of genes, Dirac delta singularities
will develop at boundary points as time evolves. Based on an energetic
variational approach (EnVarA), a balance between the maximal dissipation
principle (MDP) and least action principle (LAP), we obtain the trajectory
equation. In turn, a numerical scheme is proposed using a convex splitting
technique, with the unique solvability (on a convex set) and the energy decay
property (in time) justified at a theoretical level. Numerical examples are
presented for cases of pure drift and drift with semiselection. The remarkable
advantage of this method is its ability to catch the Dirac delta singularity
close to machine precision over any equidistant grid.

The current paper addresses the distributed guaranteedcost synchronization
problems for general highorder linear multiagent networks. Existing works on
the guaranteedcost synchronization usually require all state information of
neighboring agents and cannot give the cost budget previously. For both
leaderless and leaderfollowing interaction topologies, the current paper
firstly proposes a dynamic output feedback synchronization protocol with
guaranteedcost constraints, which can realize the tradeoff design between the
energy consumption and the synchronization regulation performance with the
given cost budget. Then, according to different structure features of
interaction topologies, leaderless and leaderfollowing guaranteedcost
synchronization analysis and design criteria are presented, respectively, and
an algorithm is proposed to deal with the impacts of nonlinear terms by using
both synchronization analysis and design criteria. Especially, an explicit
expression of the synchronization function is shown for leaderless cases, which
is independent of protocol states and the given cost budget. Finally, numerical
examples are presented to demonstrate theoretical results.

We synthesized the liquid crystal dimer and trimer members of a series of
flexible linear oligomers and characterized their microscopic and nanoscopic
properties using resonant soft xray scattering and a number of other
experimental techniques. On the microscopic scale, the twistbend phases of the
dimer and trimer appear essentially identical. However, while the liquid
crystal dimer exhibits a temperaturedependent variation of its twistbend
helical pitch varying from 100  170 {\AA} on heating, the trimer exhibits an
essentially temperatureindependent pitch of 66 {\AA}, significantly shorter
than those reported for other twistbend forming materials in the literature.
We attribute this to a specific combination of intrinsic conformational bend of
the trimer molecules and a sterically favorable intercalation of the trimers
over a commensurate fraction (twothirds) of the molecular length. We develop a
geometric model of the twistbend phase for these materials with the molecules
arranging into helical chain structures, and we fully determine their
respective geometric parameters.

We report on a method to experimentally generate ionic wind by coupling an
external high voltage electric field with an intense femtosecond laser induced
air plasma filament. The measured ionic wind velocity could be as strong as >4
m/s. It could be optimized by changing the applied electric field and the laser
induced plasma channel. The experimental observation was qualitatively
confirmed by a numerical simulation of spatial distribution of the electric
field. This technique is robust and free from sharp metallic electrodes for
coronas; it opens a way to optically generate ionic wind at a distance.

In this work, we aim at building a bridge from poor behavioral data to an
effective, quickresponse, and robust behavior model for online identity theft
detection. We concentrate on this issue in online social networks (OSNs) where
users usually have composite behavioral records, consisting of
multidimensional lowquality data, e.g., offline checkins and online user
generated content (UGC). As an insightful result, we find that there is a
complementary effect among different dimensions of records for modeling users'
behavioral patterns. To deeply exploit such a complementary effect, we propose
a joint model to capture both online and offline features of a user's composite
behavior. We evaluate the proposed joint model by comparing with some typical
models on two realworld datasets: Foursquare and Yelp. In the widelyused
setting of theft simulation (simulating thefts via behavioral replacement), the
experimental results show that our model outperforms the existing ones, with
the AUC values $0.956$ in Foursquare and $0.947$ in Yelp, respectively.
Particularly, the recall (True Positive Rate) can reach up to $65.3\%$ in
Foursquare and $72.2\%$ in Yelp with the corresponding disturbance rate (False
Positive Rate) below $1\%$. It is worth mentioning that these performances can
be achieved by examining only one composite behavior (visiting a place and
posting a tip online simultaneously) per authentication, which guarantees the
low response latency of our method. This study would give the cybersecurity
community new insights into whether and how a realtime online identity
authentication can be improved via modeling users' composite behavioral
patterns.

The progress of locationbased services has led to serious concerns on
location privacy leakage. For effective and efficient location privacy
preservation (LPP), existing methods are still not fully competent. They are
often vulnerable under the identification attack with side information, or hard
to be implemented due to the high computational complexity. In this paper, we
pursue the high protection efficacy and low computational complexity
simultaneously. We propose a \emph{scalable} LPP method based on the paradigm
of counterfeiting locations. To make fake locations extremely plausible, we
forge them through synthesizing \emph{artificial impostors} (AIs). The AIs
refer to the synthesized traces which have similar semantic features to the
actual traces, and do \emph{not} contain any target location. Two dedicated
techniques are devised: the \emph{samplingbased synthesis method} and
\emph{populationlevel semantic model}. They play significant roles in two
critical steps of synthesizing AIs. We conduct experiments on real datasets in
two cities (Shanghai, China and Asturias, Spain) to validate the high efficacy
and scalability of the proposed method. In these two datasets, the experimental
results show that our method achieves the preservation efficacy of $97.65\%$
and $96.12\%$, and its run time of building the generators is only $230.47$ and
$215.92$ seconds, respectively. This study would give the research community
new insights into improving the practicality of the stateoftheart LPP
paradigm via counterfeiting locations.

The timely provision of traffic sign information to drivers is essential for
the drivers to respond, to ensure safe driving, and to avoid traffic accidents
in a timely manner. We proposed a timely visual recognizability quantitative
evaluation method for traffic signs in largescale transportation environments.
To achieve this goal, we first address the concept of a visibility field to
reflect the visible distribution of threedimensional (3D) space and construct
a traffic sign Visibility Evaluation Model (VEM) to measure the traffic sign
visibility for a given viewpoint. Then, based on the VEM, we proposed the
concept of the Visual Recognizability Field (VRF) to reflect the visual
recognizability distribution in 3D space and established a Visual
Recognizability Evaluation Model (VREM) to measure a traffic sign visual
recognizability for a given viewpoint. Next, we proposed a Traffic Sign Timely
Visual Recognizability Evaluation Model (TSTVREM) by combining VREM, the actual
maximum continuous visual recognizable distance, and traffic big data to
measure a traffic sign visual recognizability in different lanes. Finally, we
presented an automatic algorithm to implement the TSTVREM model through traffic
sign and road marking detection and classification, traffic sign environment
point cloud segmentation, viewpoints calculation, and TSTVREM model
realization. The performance of our method for traffic sign timely visual
recognizability evaluation is tested on three road point clouds acquired by a
mobile laser scanning system (RIEGL VMX450) according to Road Traffic Signs
and Markings (GB 57681999 in China), showing that our method is feasible and
efficient.

We report here the results of searching for inelastic scattering of dark
matter (initial and final state dark matter particles differ by a small mass
splitting) with nucleon with the first 79.6day of PandaXII data (Run 9). We
set the upper limits for the spin independent WIMPnucleon scattering cross
section up to a mass splitting of 300 keV/c$^2$ at two benchmark dark matter
masses of 1 and 10 TeV/c$^2$.

We report a new search of weakly interacting massive particles (WIMPs) using
the combined low background data sets in 2016 and 2017 from the PandaXII
experiment in China. The latest data set contains a new exposure of 77.1 live
day, with the background reduced to a level of 0.8$\times10^{3}$ evt/kg/day,
improved by a factor of 2.5 in comparison to the previous run in 2016. No
excess events were found above the expected background. With a total exposure
of 5.4$\times10^4$ kg day, the most stringent upper limit on spinindependent
WIMPnucleon cross section was set for a WIMP with mass larger than 100
GeV/c$^2$, with the lowest exclusion at 8.6$\times10^{47}$ cm$^2$ at 40
GeV/c$^2$.

In this paper, we propose a recurrent neural network (RNN) with residual
attention (RRA) to learn longrange dependencies from sequential data. We
propose to add residual connections across timesteps to RNN, which explicitly
enhances the interaction between current state and hidden states that are
several timesteps apart. This also allows training errors to be directly
backpropagated through residual connections and effectively alleviates
gradient vanishing problem. We further reformulate an attention mechanism over
residual connections. An attention gate is defined to summarize the individual
contribution from multiple previous hidden states in computing the current
state. We evaluate RRA on three tasks: the adding problem, pixelbypixel MNIST
classification and sentiment analysis on the IMDB dataset. Our experiments
demonstrate that RRA yields better performance, faster convergence and more
stable training compared to a standard LSTM network. Furthermore, RRA shows
highly competitive performance to the stateoftheart methods.

The phasematching condition is a key aspect in nonlinear wavelength
conversion processes, which requires the momenta of the photons involved in the
processes to be conserved. Conventionally, nonlinear phase matching is achieved
using either birefringent or periodically poled nonlinear crystals, which
requires careful dispersion engineering and are usually narrowband. In recent
years, metasurfaces consisting of densely packed arrays of optical antennas
have been demonstrated to provide an effective optical momentum to bend light
in arbitrary ways. Here, we demonstrate that gradient metasurface structures
consisting of phased array antennas are able to circumvent the phasematching
requirement in onchip nonlinear wavelength conversion. We experimentally
demonstrate phasematchingfree second harmonic generation over many coherent
lengths in thin film lithium niobate waveguides patterned with the gradient
metasurfaces. Efficient secondharmonic generation (1660% W1cm2) in the
metasurfacebased devices was observed over a wide range of pump wavelengths
(15801650 nm).

In this paper, we study a novel secondorder energy stable Backward
Differentiation Formula (BDF) finite difference scheme for the epitaxial thin
film equation with slope selection (SS). One major challenge for the higher
oder in time temporal discretization is how to ensure an unconditional energy
stability and an efficient numerical implementation. We propose a general
framework for designing the higher order in time numerical scheme with
unconditional energy stability by using the BDF method with constant
coefficient stabilized terms. Based on the unconditional energy stability
property, we derive an $L^\infty_h (0,T; H_{h}^2)$ stability for the numerical
solution and provide an optimal the convergence analysis. To deal with the
4Laplacian solver in an $L^{2}$ gradient flow at each time step, we apply an
efficient preconditioned steepest descent algorithm and preconditioned
nonlinear conjugate gradient algorithm to solve the corresponding nonlinear
system. Various numerical simulations are present to demonstrate the stability
and efficiency of the proposed schemes and slovers.

Laser wakefield accelerators (LWFA) hold great potential to produce
highquality highenergy electron beams (e beams) and simultaneously bright
xray sources via betatron radiation, which are very promising for pumpprobe
study in ultrafast science. However, in order to obtain a highquality e beam,
electron injection and acceleration should be carefully manipulated, where a
large oscillation amplitude has to be avoided and thus the emitted xray yield
is limited. Here, we report a new scheme to experimentally enhance betatron
radiation significantly both in photon yield and photon energy by separating
electron injection and acceleration from manipulation of the ebeam transverse
oscillation in the wake via introducing a slanted thin plasma refraction slab.
Particleincell simulations indicate that the ebeam transverse oscillation
amplitude can be increased by more than 10 folds, after being steered into the
deflected laserdriven wakefield due to refraction at the slab's boundaries.
Spectral broadening of the xrays can be suppressed owing to the small
variation in the peak energy of the lowenergyspread e beam in a plasma
wiggler regime. We demonstrate that the highquality ebeam generation,
refracting and wiggling can act as a whole to realize the concurrence of
monoenergetic e beam and bright xrays in a compact LWFA.

Periodic structures of phases with orientational order of molecules, but
homogenous electron density distribution: a short pitch cholesteric, blue phase
and twistbend nematic phase, were probed by a resonant soft xray scattering
(RSoXS) at the carbon Kedge. The theoretical model shows that in case of a
simple heliconical nematic structure two resonant signals corresponding to the
full and half pitch band should be present, while only the full pitch band is
observed in experiment. This suggests that the twistbend nematic phase has
complex structure with a doublehelix, built of two interlocked, shifted
helices. We confirm that the helical pitch in the twistbend nematic phase is
in a 10 nm range, for both, the chiral and achiral materials. We also show that
the symmetry of a blue phase can unambiguously be determined through a resonant
enhancement of xray diffraction signals, by including polarization effects,
which are found to be an important indicator in phase structure determination.

In this paper, we propose a novel, thermodynamically consistent phase field
model to simulate the deformation and breakup of a ferrodroplet that is
immersed in a viscous medium and subject to an applied uniform magnetic field.
Instead of using the magnetic body force in the traditional Rosensweig model,
the key idea of this model is to propose a new magnetic energy that enables
direct effects of the magnetic field on the interface evolution. The model can
thereby be derived from the variational principle via minimizing the free
energy of the total system. This energy based modeling idea can be easily
extended to include more external fields in a convenient and consistent way for
more complicated applications involving multiple external fields. We validate
the model by performing a series of numerical simulations, including the
comparison with analytic solutions, the investigation of the effect of
different types of magnetic fields, the dynamical behaviors of the ferrodroplet
breakup under the strong magnetic field, magnetic/velocity/pressure field
distributions, the magnetic energy density, and the inertial phenomenon, etc.

Modern communication networks require high performance and scalable
electrooptic modulators that convert electrical signals to optical signals at
high speed. Existing lithium niobate modulators have excellent performance but
are bulky and prohibitively expensive to scale up. Here we demonstrate scalable
and highperformance nanophotonic electrooptic modulators made of
singlecrystalline lithium niobate microring resonators and microMachZehnder
interferometers. We show a halfwave electrooptic modulation efficiency of
1.8V$\cdot$cm and data rates up to 40 Gbps.

Integrated thinfilm lithium niobate platform has recently emerged as a
promising candidate for nextgeneration, highefficiency wavelength conversion
systems that allow dense packaging and massproduction. Here we demonstrate
efficient, phasematched second harmonic generation in lithographicallydefined
thinfilm lithium niobate waveguides with submicron dimensions. Both modal
phase matching in fixedwidth waveguides and quasiphase matching in
periodically grooved waveguides are theoretically proposed and experimentally
demonstrated. Our lowloss (~ 3.0 dB/cm) nanowaveguides possess normalized
conversion efficiencies as high as 41% W1cm2.

We describe and analyze preconditioned steepest descent (PSD) solvers for
fourth and sixthorder nonlinear elliptic equations that include pLaplacian
terms on periodic domains in 2 and 3 dimensions. The highest and lowest order
terms of the equations are constantcoefficient, positive linear operators,
which suggests a natural preconditioning strategy. Such nonlinear elliptic
equations often arise from time discretization of parabolic equations that
model various biological and physical phenomena, in particular, liquid
crystals, thin film epitaxial growth and phase transformations. The analyses of
the schemes involve the characterization of the strictly convex energies
associated with the equations. We first give a general framework for PSD in
generic Hilbert spaces. Based on certain reasonable assumptions of the linear
preconditioner, a geometric convergence rate is shown for the nonlinear PSD
iteration. We then apply the general the theory to the fourth and sixthorder
problems of interest, making use of Sobolev embedding and regularity results to
confirm the appropriateness of our preconditioners for the regularized
pLapacian problems. Our results include a sharper theoretical convergence
result for pLaplacian systems compared to what may be found in existing works.
We demonstrate rigorously how to apply the theory in the finite dimensional
setting using finite difference discretization methods. Numerical simulations
for some important physical application problems  including thin film epitaxy
with slope selection and the square phase field crystal model  are carried
out to verify the efficiency of the scheme.

In this paper we analyze and implement a secondorderintime numerical
scheme for the threedimensional phase field crystal (PFC) equation. The
numerical scheme was proposed in [46], with the unique solvability and
unconditional energy stability established. However, its convergence analysis
remains open. We present a detailed convergence analysis in this article, in
which the maximum norm estimate of the numerical solution over grid points
plays an essential role. Moreover, we outline the detailed multigrid method to
solve the highly nonlinear numerical scheme over a cubic domain, and various
threedimensional numerical results are presented, including the numerical
convergence test, complexity test of the multigrid solver and the polycrystal
growth simulation.

We present a secondorderintime finite difference scheme for the
CahnHilliardHeleShaw equations. This numerical method is uniquely solvable
and unconditionally energy stable. At each time step, this scheme leads to a
system of nonlinear equations that can be efficiently solved by a nonlinear
multigrid solver. Owing to the energy stability, we derive an $\ell^2 (0,T;
H_h^3)$ stability of the numerical scheme. To overcome the difficulty
associated with the convection term $\nabla \cdot (\phi \boldsymbol{u})$, we
perform an $\ell^\infty (0,T; H_h^1)$ error estimate instead of the classical
$\ell^\infty (0,T; \ell^2)$ one to obtain the optimal rate convergence
analysis. In addition, various numerical simulations are carried out, which
demonstrate the accuracy and efficiency of the proposed numerical scheme.

We discuss how to efficiently forward data in vehicular networks. Existing
solutions do not make full use of trajectory planning of nearby vehicles, or
social attributes. The development of onboard navigation system provides
drivers some traveling route information. The main novelty of our approach is
to envision sharing partial traveling information to the encountered vehicles
for better service. Our data forwarding algorithm utilizes this lightweight
information under the delusive paths privacy preservation together with the
social community structure in vehicular networks. We assume that data
transmission is carried by vehicles and road side units (RSUs), while cellular
network manages and coordinates relevant global information. The approximate
destination set is the set of RSUs that are often passed by the destination
vehicle. RSU importance is raised by summing encounter ratios of RSUs in the
same connected component. We first define a concept of spacetime
approachability which is derived from shared partial traveling route and
encounter information. It describes the capability of a vehicle to advance
messages toward destination. Then, we design a novel data forwarding algorithm,
called approachability based algorithm, which combines the spacetime
approachability with the social community attribute in vehicular networks. We
evaluate our approachability based algorithm on data sets from San Francisco
Cabspotting and Shanghai Taxi Movement. Results show that the partially shared
traveling information plays a positive role in data forwarding in vehicular
networks. Approachability based data forwarding algorithm achieves a better
performance than existing social based algorithms in vehicular networks.

In this paper we provide a detailed convergence analysis for fully discrete
second order (in both time and space) numerical schemes for nonlocal AllenCahn
(nAC) and nonlocal CahnHilliard (nCH) equations. The unconditional unique
solvability and energy stability ensures $\ell^4$ stability. The convergence
analysis for the nAC equation follows the standard procedure of consistency and
stability estimate for the numerical error function. For the nCH equation, due
to the complicated form of the nonlinear term, a careful expansion of its
discrete gradient is undertaken and an $H^{1}$ inner product estimate of this
nonlinear numerical error is derived to establish convergence. In addition, an
apriori $W^{1,\infty}$ bound of the numerical solution at the discrete level
is needed in the error estimate. Such a bound can be obtained by performing a
higher order consistency analysis by using asymptotic expansions for the
numerical solution. Following the technique originally proposed by Strang
(e.g., 1964), instead of the standard comparison between the exact and
numerical solutions, an error estimate between the numerical solution and the
constructed approximate solution yields an $O( s^3 + h^4)$ convergence in
$\ell^\infty (0, T; \ell^2)$ norm, which leads to the necessary bound under a
standard constraint $s \le C h$. Here, we also prove convergence of the scheme
in the maximum norm under the same constraint.

Geometric model fitting is a fundamental task in computer graphics and
computer vision. However, most geometric model fitting methods are unable to
fit an arbitrary geometric model (e.g. a surface with holes) to incomplete
data, due to that the similarity metrics used in these methods are unable to
measure the rigid partial similarity between arbitrary models. This paper hence
proposes a novel rigid geometric similarity metric, which is able to measure
both the full similarity and the partial similarity between arbitrary geometric
models. The proposed metric enables us to perform partial procedural geometric
model fitting (PPGMF). The task of PPGMF is to search a procedural geometric
model space for the model rigidly similar to a query of noncomplete point set.
Models in the procedural model space are generated according to a set of
parametric modeling rules. A typical query is a point cloud. PPGMF is very
useful as it can be used to fit arbitrary geometric models to noncomplete
(incomplete, overcomplete or hybridcomplete) point cloud data. For example,
most laser scanning data is noncomplete due to occlusion. Our PPGMF method
uses Markov chain Monte Carlo technique to optimize the proposed similarity
metric over the model space. To accelerate the optimization process, the method
also employs a novel coarsetofine model dividing strategy to reject
dissimilar models in advance. Our method has been demonstrated on a variety of
geometric models and noncomplete data. Experimental results show that the
PPGMF method based on the proposed metric is able to fit noncomplete data,
while the method based on other metrics is unable. It is also shown that our
method can be accelerated by several times via early rejection.

This paper proposes a convex optimization based distributed algorithm to
solve multiperiod optimal gaspower flow (OGPF) in coupled energy distribution
systems. At the gas distribution system side, the nonconvex Weymouth gas flow
equations is convexified as quadratic constraints. The optimal gas flow (OGF)
subproblem is solved by an iterative secondorder cone programming procedure,
whose efficiency is two orders of magnitudes higher than traditional nonlinear
methods. A convex quadratic program based initiation scheme is suggested, which
helps to find a highquality starting point. At the power distribution system
side, convex relaxation is performed on the nonconvex branch flow equations,
and the optimal power flow (OPF) subproblem gives rise to a second order cone
program. Tightness is guaranteed by the radial topology. In the proposed
distributed algorithm, OGF and OPF are solved independently, and coordinated by
the alternating direction multiplier method (ADMM). Numerical results
corroborate significant enhancements on computational robustness and efficiency
compared with existing centralized OGPF methods.