
This paper considers how to obtain MCMC quantitative convergence bounds which
can be translated into tight complexity bounds in highdimensional settings. We
propose a modified driftandminorization approach, which establishes a
generalized drift condition defined in a subset of the state space. The subset
is called the ``large set'' and is chosen to rule out some ``bad'' states which
have poor drift property when the dimension gets large. Using the ``large set''
together with a ``centered'' drift function, a quantitative bound can be
obtained which can be translated into a tight complexity bound. As a
demonstration, we analyze a certain realistic Gibbs sampler algorithm and
obtain a complexity upper bound for the mixing time, which shows that the
number of iterations required for the Gibbs sampler to converge is constant
under certain conditions on the observed data and the initial state. It is our
hope that this modified driftandminorization approach can be employed in many
other specific examples to obtain complexity bounds for highdimensional Markov
chains.

In this paper, lower bounds on the probability of a finite union of events
are considered, i.e. $P\left(\bigcup_{i=1}^N A_i\right)$, in terms of the
individual event probabilities $\{P(A_i), i=1,\ldots,N\}$ and the sums of the
pairwise event probabilities, i.e., $\{\sum_{j:j\neq i} P(A_i\cap A_j),
i=1,\ldots,N\}$. The contribution of this paper includes the following: (i) in
the class of all lower bounds that are established in terms of only the
$P(A_i)$'s and $\sum_{j:j\neq i} P(A_i\cap A_j)$'s, the optimal lower bound is
given numerically by solving a linear programming (LP) problem with $N^2N+1$
variables; (ii) a new analytical lower bound is proposed based on a relaxed LP
problem, which is at least as good as the bound due to Kuai, et al.; (iii)
numerical examples are provided to illustrate the performance of the bounds.

We report the detection of the radio afterglow of a long gammaray burst
(GRB) 111005A at 5345 GHz, including the very long baseline interferometry
observations with the positional error of 0.2 mas. The afterglow position is
coincident with the disk of a galaxy ESO 58049 at z= 0.01326 (~1" from its
center), which makes GRB 111005A the second closest GRB known to date, after
GRB 980425. The radio afterglow of GRB 111005A was an order of magnitude less
luminous than those of local lowluminosity GRBs, and obviously than those of
cosmological GRBs. The radio flux was approximately constant and then
experienced an unusually rapid decay a month after the GRB explosion. Similarly
to only two other GRBs, we did not find the associated supernovae (SN), despite
deep near and midinfrared observations 19 days after the GRB explosion,
reaching ~20 times fainter than other SNe associated with GRBs. Moreover, we
measured twice solar metallicity for the GRB location. The low gammaray and
radio luminosities, rapid decay, lack of a SN, and supersolar metallicity
suggest that GRB 111005A represents a different rare class of GRBs than typical
corecollapse events. We modelled the spectral energy distribution of the GRB
111005A host finding that it is a dwarf, moderately starforming galaxy,
similar to the host of GRB 980425. The existence of two local GRBs in such
galaxies is still consistent with the hypothesis that the GRB rate is
proportional to the cosmic star formation rate (SFR) density, but suggests that
the GRB rate is biased towards low SFRs. Using the farinfrared detection of
ESO 58049, we conclude that the hosts of both GRBs 111005A and 980425 exhibit
lower dust content than what would be expected from their stellar masses and
optical colours.

We report the experimental realization of Dirac semimetal state in NdSb, a
material with antiferromagnetic ground state. The occurrence of topological
semimetal state has been well supported by our band structure calculations and
the experimental observation of chiral anomaly induced negative
magnetoresistance. A fieldinduced Fermi surface reconstruction is observed, in
response to the change of spin polarization. The observation of topological
semimetal state in a magnetic material provides an opportunity to investigate
the magnetotopological phenomena.

We present the design, architecture and detailed performance of a
threedimensional (3D) underwater acoustic carpet cloak (UACC). The proposed
system of the 3D UACC is an octahedral pyramid which is composed of periodical
steel strips. This underwater acoustic device, placed over the target to hide,
is able to manipulate the scattered wavefront to mimic a reflecting plane. The
effectiveness of the prototype is experimentally demonstrated in an anechoic
tank. The measured acoustic pressure distributions show that the 3D UACC can
work in all directions in a wide frequency range. This experimental
verification of 3D device paves the way for guidelines on future practical
applications.

Newton's gravitational constant $G$ may vary with time at an extremely low
level. The time variability of $G$ will affect the orbital motion of a
millisecond pulsar in a binary system and cause a tiny difference between the
orbital perioddependent measurement of the kinematic distance and the direct
measurement of the annual parallax distance. PSR J0437$$4715 is the nearest
millisecond pulsar and the brightest at radio. To explore the feasibility of
achieving a parallax distance accuracy of one lightyear, comparable to the
recent timing result, with the technique of differential astrometry, we
searched for compact radio sources quite close to PSR J0437$$4715. Using
existing data from the Very Large Array and the Australia Telescope Compact
Array, we detected two sources with flat spectra, relatively stable flux
densities of 0.9 and 1.0 mJy at 8.4 GHz and separations of 13 and 45 arcsec.
With a network consisting of the Long Baseline Array and the Kunming 40m radio
telescope, we found that both sources have a pointlike structure and a
brightness temperature of $\geq$10$^7$ K. According to these radio inputs and
the absence of counterparts in the other bands, we argue that they are most
likely the compact radio cores of extragalactic active galactic nuclei rather
than Galactic radio stars. The finding of these two radio active galactic
nuclei will enable us to achieve a subpc distance accuracy with the inbeam
phasereferencing verylongbaseline interferometric observations and provide
one of the most stringent constraints on the time variability of $G$ in the
near future.

The $\gamma$ray bright blazar CTA 102 is studied using imaging (new 15 GHz
and archival 43 GHz Very Long Baseline Array, VLBA data) and time variable
optical flux density, polarization degree and electric vector position angle
(EVPA) spanning between 2015 June 1 and 2016 October 1, covering a prominent
$\gamma$ray flare during 2016 January. The pcscale jet indicates expansion
with oscillatory features upto 17 mas. Component proper motions are in the
range 0.04  0.33 mas/yr with acceleration upto 1.2 mas followed by a slowing
down beyond 1.5 mas. A jet bulk Lorentz factor $\geq$ 17.5, position angle of
128.3 degrees, inclination angle $\leq$ 6.6 degrees and intrinsic half opening
angle $\leq$ 1.8 degrees are derived from the VLBA data. These inferences are
employed in a helical jet model to infer long term variability in flux density,
polarization degree, EVPA and a rotation of the Stokes Q and U parameters. A
core distance of $r_{\rm core,43 \ GHz}$ = 22.9 pc, and a magnetic field
strength at 1 pc and the core location of 1.57 G and 0.07 G respectively are
inferred using the core shift method. The study is useful in the context of
estimating jet parameters and in offering clues to distinguish mechanisms
responsible for variability over different timescales.

This is a short survey on existing upper and lower bounds on the probability
of the union of a finite number of events using partial information given in
terms of the individual or pairwise event probabilities (or their sums). New
proofs for some of the existing bounds are provided and new observations
regarding the existing GallotKounias bound are given.

Nonlinear Schr\"odinger equation, short pulse equation and complex short
pulse equation have important application in nonlinear optics. They can be
derived from the Maxwell equation. In this paper, we investigate a coupled
focusingdefocusing complex short pulse equation.
The brightbright, brightdark and darkdark soliton solutions of the coupled
focusingdefocusing complex short pulse equation are given.
Then the breathers are derived from the dark soliton solution. The rogue wave
solutions are also constructed.
The dynamics and the asymptotic behavior of the soliton solutions are
analyzed, which reveals that there exist the elastic or inelastic collision in
brightbright soliton solution. But the interactions of brightdark and
darkdark soliton solutions are both elastic.

We design an impedance matching underwater acoustic bend with pentamode
microstructure. The proposed bend is assembled by pentamode lattice. The
effective density and compressive mod ulus of each unit cell can be tuned
simultaneously, which are modulated to guarantee both the bending effect and
high transmission. The standard deviations of transmitted phase are calculated
to quantitatively evaluate the degree of the distortion of the transmitted
wavefront, while the trans mission is calculated to appraise the degree of
acoustic impedance matching. The low standard deviations and high transmission
indicate that the designed bend has a nice broadband bending effect and is
impedancematched to water. This design has potential applications in
underwater communication and underwater detection.

Sound absorption or dissipation principally involves joint interactions
between sound waves, material morphology and the air medium. How these elements
work most efficiently for sound absorption remains elusive to date. In this
paper, we suggest a fundamental relation concisely crosslinking the three
elements, which reveals that optimal sound absorption efficiency occurs when
the pore size of the material is twice the thickness of the viscous boundary
layer of the acoustic air medium. The study is validated by microlattice
materials comprising of wellcontrolled regular structures that absorb sound in
a tunable manner. Optimized material morphology in terms of pore size and
porosity is determined to provide a robust guidance for optimizing sound
absorbing materials.

This work presents the mathematical and computational aspects of a smooth
dissipative particle dynamics with dynamic virtual particle allocation method
(SDPDDV) for modeling and simulation of mesoscopic fluids in wallbounded
domains. The SDPDDV method is realized with fluid particles, boundary
particles and dynamically allocated virtual particles near solid boundaries.
The physical domain in SDPDDV contains external and internal solid boundaries,
periodic inlets and outlets, and the fluid region. The solid boundaries of the
domain are represented with boundary particles which have an assigned position,
wall velocity, and temperature upon initialization. The fluid domain is
discretized with fluid particles placed in a global index. The algorithm for
nearest neighbor particle search is based on a combination of the linkedcell
and Verletlist approaches and utilizes large rectangular cells for
computational efficiency. The density model of a fluid particle in the
proximity of a solid boundary includes the contribution from the virtual
particles in its truncated support domain. The thermodynamic properties of a
virtual particle are identical to those of the corresponding fluid particle.
Additional verification involves SDPDDV simulations of transient, planar,
Couette liquid water flow. The top plate is moving at and separated by 103 m
from the bottom stationary plate. The numerical results are in very good
agreement with the analytical solutions. Additional SDPDDV verification is
accomplished with the simulation of a bodyforce driven, lowReynolds number
flow of water over a cylinder of radius . The SDPDDV field velocity and
pressure are compared with those obtained by FLUENT.

In this work, we review previously developed coarsegrained (CG) particle
models for biological membrane and red blood cells (RBCs) and discuss the
advantages of the CG particle method over the continuum and atomic simulations
on modeling biological phenomena. CG particle models can largely increase the
length scale and time scale of atomic simulations by eliminating fast degrees
of freedom while preserving the mesoscopic structures and properties of the
simulated system. One the other hand, CG particle models can be used to capture
microstructural alternations in diseased RBCs and simulate topological changes
of biological membrane and RBCs, which are major challenges to typical
continuum representations of membrane and RBCs. The power and versatility of
the CG particle methods are demonstrated through simulating the dynamical
processes involving significant topological changes, such as lipid
selfassembly, vesicle fusion and membrane budding.

PKS 1155+251 is a radioloud quasar source at z=0.203. Observations using
very long baseline interferometry (VLBI) at ~2, 5, 8 and 15 GHz show that the
structure of the radio source is quite complicated on parsec scales and that
the outer hot spots are apparently undergoing a significant contraction.
Because these results cannot be fully explained based on the compact symmetric
object (CSO) scenario with a radio core located between the northern and
southern complexes, we made observations with the Very Long Baseline Array
(VLBA) at 24 and 43 GHz to search for compact substructures and alternative
interpretations. The results show that the radio core revealed in the previous
VLBI observations remains compact with a flat spectrum in our
submilliarcsecondresolution images; the northern lobe emission becomes
faint at 24 GHz and is mostly resolving out at 43 GHz; the southern complex is
more bright but has been resolved into the brightest southernend (S1) and jet
or tail alike components westwards. Explaining the southern components aligned
westward with a standard CSO scenario alone remains a challenge. As for the
flatter spectral index of the southernend component S1 between 24 and 43 GHz
in our observations and the significant 15 GHz VLBA flux variability of S1, an
alternative scenario is that the southern complex may be powered by a secondary
black hole residing at S1. But more sensitive and highresolution VLBI
monitoring is required to discriminate the CSO and the binary black hole
scenarios.

We show that the transition between a ETH phase and a manybody localized
phase is marked by the different finite size scaling behaviour of the decay of
the Loschmidt Echo and its temporal fluctuations  after a quantum quench  in
the infinite time limit, despite the fact that the finite time behaviour of
such quantities is dramatically different approach the MBL phase, so that
temporal fluctuations cannot be inferred from the infinite time average of the
Loschmidt Echo. We also show the different scrambling powers of ETH and MBL
Hamiltonians as a probe to the different approaches to equilibrium.

Semisupervised learning plays an important role in largescale machine
learning. Properly using additional unlabeled data (largely available nowadays)
often can improve the machine learning accuracy. However, if the machine
learning model is misspecified for the underlying true data distribution, the
model performance could be seriously jeopardized. This issue is known as model
misspecification. To address this issue, we focus on generative models and
propose a criterion to detect the onset of model misspecification by measuring
the performance difference between models obtained using supervised and
semisupervised learning. Then, we propose to automatically modify the
generative models during model training to achieve an unbiased generative
model. Rigorous experiments were carried out to evaluate the proposed method
using two image classification data sets PASCAL VOC'07 and MIR Flickr. Our
proposed method has been demonstrated to outperform a number of
stateoftheart semisupervised learning approaches for the classification
task.

An effective way to suppress the cascading failure risk is the branch
capacity upgrade, whose optimal decision making, however, may incur high
computational burden. A practical way is to find out some critical branches as
the candidates in advance. This paper proposes a simulation data oriented
approach to identify the critical branches with higher importance in cascading
failure propagation. First, a concept of cascading failure chain (CFC) is
introduced and numerous samples of CFC are generated with an AC power flow
based cascading failure simulator. Then, a directed weighted graph is
constructed, whose edges denotes the severities of branch interactions. Third,
the weighted hypertextinduced topic search (HITS) algorithm is used to rate
and rank this graph's vertices,through which the critical branches can be
identified accordingly. Validations on IEEE 118bus and RTS96 systems show that
the proposed approach can identify critical branches whose capacity upgrades
suppress cascading failure risk more greatly. Moreover, it is also shown that
structural importance of a branch does not agree with its importance in
cascading failure, which indicates the effectiveness of the proposed approach
compared with structure vulnerabilities based identifying methods.

The 80 highmass Xray binary (HMXB) pulsars that are known to reside in the
Magellanic Clouds (MCs) have been observed by the XMMNewton and Chandra Xray
telescopes on a regular basis for 15 years, and the XMMNewton and Chandra
archives contain nearly complete information about the duty cycles of the
sources with spin periods P_S < 100 s. We have rerprocessed the archival data
from both observatories and we combined the output products with all the
published observations of 31 MC pulsars with P_S < 100 s in an attempt to
investigate the faintest Xray emission states of these objects that occur when
accretion to the polar caps proceeds at the smallest possible rates. These
states determine the socalled propeller lines of the accreting pulsars and
yield information about the magnitudes of their surface magnetic fields. We
have found that the faintest states of the pulsars segregate into five discrete
groups which obey to a high degree of accuracy the theoretical relation between
spin period and Xray luminosity. So the entire population of these pulsars can
be described by just five propeller lines and the five corresponding magnetic
moments (0.29, 0.53, 1.2, 2.9, and 7.3, in units of 10^30 G cm^3).

We propose the design of an impedance matching acoustic bend in this article.
The bending structure is composed of subwavelength unit cells with perforated
plates and side pipes, whose mass density and bulk modulus can be tuned
simultaneously. So the refraction index and the impedance of the acoustic bend
can be modulated simultaneously to guarantee both the bending effect and the
high transmission. The simulation results of sound pressure field distribution
show that the bending effect of the impedance matching acoustic bend is very
good. Transmission spectra of the impedance matching acoustic bend and the
acoustic bend composed of perforated plates only are both calculated for
comparison. The results indicate that the impedance matching acoustic bend is
successful in improving the impedance ratio and increasing the transmission
obviously.

In this paper, we design, fabricate and experimentally characterize a
broadband acoustic rightangle bend device in air. Perforated panels with
various holesizes are used to construct the bend structure. Both the simulated
and the experimental results verify that acoustic beam can be rotated
effectively through the acoustic bend in a wide frequency range. This model may
have potential applications in some areas such as sound absorption and acoustic
detection in pipeline.

We present a new second order complete active space selfconsistent field
implementation to converge wavefunctions for both large active spaces and large
atomic orbital (AO) bases. Our algorithm decouples the active space
wavefunction solver from the orbital optimization in the microiterations, and
thus may be easily combined with various modern active space solvers. We also
introduce efficient approximate orbital gradient and Hessian updates, and step
size determination. We demonstrate its capabilities by calculating the
lowlying states of the Fe(\Roman{2})porphine complex with modest resources
using a density matrix renormalization group solver in a CAS(22,27) active
space and a 3000 AO basis.

Application of transformation theory to underwater acoustics has been a
challenging task because highly anisotropic density is unachievable in water. A
possible strategy is to exploit anisotropic modulus rather than density, while
has not been experimentally demonstrated. We present an annular underwater
acoustic cloak designed from particular graded solid microstructures. The
geometry tailored microstructures mimics metafluid with highly anisotropic
modulus through substantially suppressed shear wave. Transient wave experiments
are conducted with the cloak in a designed 2D underwater waveguide system and
proved excellent cloaking performance for enclosed target over broadband
frequency 915 kHz. This finding paves the way for controlling underwater
acoustics using the structured anisotropic modulus metafluid.

According to optical stellar kinematics observations, an overmassive black
hole candidate has been reported by van den Bosch et al. (2012) in the normal
earlytype galaxy NGC 1277. This galaxy is located in the central region of the
Perseus cluster. Westerbork Synthesis Radio Telescope (WSRT) observations have
shown that NGC 1277 and other earlytype galaxies in the neighbourhood have
radio counterparts. These nuclear radio sources have stable flux densities on
time scale of years. In order to investigate the origin of the radio emission
from these normal galaxies, we selected five sources (NGC 1270, NGC 1272, NGC
1277, NGC 1278 and VZw 339) residing in the central 10 arcminute region of the
Perseus cluster and requested to recorrelate the data of an existing very long
baseline interferometry (VLBI) experiment at these new positions. With the
recorrelation data provided by the European VLBI Network (EVN), we imaged the
five sources with a resolution of about eight milliarcseconds and detected all
of them with a confidence level above 5{\sigma} at 1.4 GHz. They show compact
structure and brightness temperatures above $10^7$ K, which implies that the
radio emission is nonthermal. We rule out ongoing nuclear star formation and
conclude that these VLBIdetected radio sources are parsecscale jet activity
associated with the supermassive black holes in low luminosity active galactic
nuclei (LLAGNs), although there are no clear signs of nuclear activity observed
in the optical and infrared bands. Using the fundamental plane relation in
black holes, we find no significant evidence for or against an extremely
massive black hole hiding in NGC 1277.

We present an unsupervised method for cosegmentation of a set of 3D shapes
from the same class with the aim of segmenting the input shapes into consistent
semantic parts and establishing their correspondence across the set. Starting
from meaningful presegmentation of all given shapes individually, we construct
the correspondence between same candidate parts and obtain the labels via
functional maps. And then, we use these labels to mark the input shapes and
obtain results of cosegmentation. The core of our algorithm is to seek for an
optimal correspondence between semantically similar parts through functional
maps and mark such shape parts. Experimental results on the benchmark datasets
show the efficiency of this method and comparable accuracy to the
stateoftheart algorithms.

We consider the problem $$
\epsilon^2 \Delta
uV(y)u+u^p\,=\,0,~~u>0~~\quad\mbox{in}\quad\Omega,~~\quad\frac {\partial
u}{\partial \nu}\,=\,0\quad\mbox{on}~~~\partial \Omega, $$ where $\Omega$ is a
bounded domain in $\mathbb R^2$ with smooth boundary, the exponent $p>1$,
$\epsilon>0$ is a small parameter, $V$ is a uniformly positive, smooth
potential on $\bar{\Omega}$, and $\nu$ denotes the outward normal of $\partial
\Omega$. Let $\Gamma$ be a curve intersecting orthogonally with $\partial
\Omega$ at exactly two points and dividing $\Omega$ into two parts. Moreover,
$\Gamma$ satisfies stationary and nondegeneracy conditions with respect to the
functional $\int_{\Gamma}V^{\sigma}$, where $\sigma=\frac {p+1}{p1}\frac 12$.
We prove the existence of a solution $u_\epsilon$ concentrating along the whole
of $\Gamma$, exponentially small in $\epsilon$ at any positive distance from
it, provided that $\epsilon$ is small and away from certain critical numbers.
In particular, this establishes the validity of the two dimensional case of a
conjecture by A. Ambrosetti, A. Malchiodi and W.M. Ni(p.327, [4]).