• This paper considers how to obtain MCMC quantitative convergence bounds which can be translated into tight complexity bounds in high-dimensional settings. We propose a modified drift-and-minorization 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 drift-and-minorization approach can be employed in many other specific examples to obtain complexity bounds for high-dimensional 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^2-N+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 gamma-ray burst (GRB) 111005A at 5-345 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 580-49 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 low-luminosity 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 mid-infrared observations 1-9 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 gamma-ray and radio luminosities, rapid decay, lack of a SN, and super-solar metallicity suggest that GRB 111005A represents a different rare class of GRBs than typical core-collapse events. We modelled the spectral energy distribution of the GRB 111005A host finding that it is a dwarf, moderately star-forming 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 far-infrared detection of ESO 580-49, 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 field-induced 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 magneto-topological phenomena.
  • We present the design, architecture and detailed performance of a three-dimensional (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 period-dependent 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 light-year, 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 40-m radio telescope, we found that both sources have a point-like 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 sub-pc distance accuracy with the in-beam phase-referencing very-long-baseline 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 pc-scale 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 Gallot--Kounias 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 focusing-defocusing complex short pulse equation. The bright-bright, bright-dark and dark-dark soliton solutions of the coupled focusing-defocusing 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 bright-bright soliton solution. But the interactions of bright-dark and dark-dark 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 impedance-matched 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 cross-linking 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 well-controlled 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 (SDPD-DV) for modeling and simulation of mesoscopic fluids in wall-bounded domains. The SDPD-DV method is realized with fluid particles, boundary particles and dynamically allocated virtual particles near solid boundaries. The physical domain in SDPD-DV 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 linked-cell and Verlet-list 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 SDPD-DV simulations of transient, planar, Couette liquid water flow. The top plate is moving at and separated by 10-3 m from the bottom stationary plate. The numerical results are in very good agreement with the analytical solutions. Additional SDPD-DV verification is accomplished with the simulation of a body-force driven, low-Reynolds number flow of water over a cylinder of radius . The SDPD-DV field velocity and pressure are compared with those obtained by FLUENT.
  • In this work, we review previously developed coarse-grained (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 self-assembly, vesicle fusion and membrane budding.
  • PKS 1155+251 is a radio-loud 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 sub-milli-arcsecond--resolution 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 southern-end (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 southern-end 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 high-resolution 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 many-body 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.
  • Semi-supervised learning plays an important role in large-scale 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 semi-supervised 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 state-of-the-art semi-supervised 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 hypertext-induced 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 high-mass X-ray binary (HMXB) pulsars that are known to reside in the Magellanic Clouds (MCs) have been observed by the XMM-Newton and Chandra X-ray telescopes on a regular basis for 15 years, and the XMM-Newton 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 X-ray emission states of these objects that occur when accretion to the polar caps proceeds at the smallest possible rates. These states determine the so-called 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 X-ray 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 sub-wavelength 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 right-angle bend device in air. Perforated panels with various hole-sizes 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 self-consistent 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 low-lying 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 meta-fluid 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 9-15 kHz. This finding paves the way for controlling underwater acoustics using the structured anisotropic modulus meta-fluid.
  • According to optical stellar kinematics observations, an over-massive black hole candidate has been reported by van den Bosch et al. (2012) in the normal early-type 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 early-type 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 re-correlate the data of an existing very long baseline interferometry (VLBI) experiment at these new positions. With the re-correlation 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 non-thermal. We rule out ongoing nuclear star formation and conclude that these VLBI-detected radio sources are parsec-scale 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 co-segmentation 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 pre-segmentation 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 co-segmentation. 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 state-of-the-art algorithms.
  • We consider the problem $$ \epsilon^2 \Delta u-V(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 non-degeneracy conditions with respect to the functional $\int_{\Gamma}V^{\sigma}$, where $\sigma=\frac {p+1}{p-1}-\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]).