title: Generalized Papapetrou's equations of motion for an extended test body within static and isotropic metrics

abstract: Applying Dixon's general equations of motion for extended bodies, we compute the Papapetrou's equations for an extended test body on static and isotropic metrics. We incorporate the force and the torque terms which involve multipole moments, beyond dipole moment, from the energy-momentum tensor. We obtain the vector form equations for both Corinaldesi-Papapetrou and Tulczyjew-Dixon spin supplementary conditions. An expanded effective mass, including interactions between the structure of the body and the gravitational field, is also found.


title: Design and construction of a high-energy photon polarimeter

abstract: We report on the design and construction of a high-energy photon polarimeter for measuring the degree of polarization of a linearly-polarized photon beam. The photon polarimeter uses the process of pair production on an atomic electron (triplet production). The azimuthal distribution of scattered atomic electrons following triplet production yields information regarding the degree of linear polarization of the incident photon beam. The polarimeter, operated in conjunction with a pair spectrometer, uses a silicon strip detector to measure the recoil electron distribution resulting from triplet photoproduction in a beryllium target foil. The analyzing power $\Sigma_A$ for the device using a 75 $\rm{\mu m}$ beryllium converter foil is about 0.2, with a relative systematic uncertainty in $\Sigma_A$ of 1.5%.


title: An iBeacon based Proximity and Indoor Localization System

abstract: Indoor localization and Location Based Services (LBS) can greatly benefit from the widescale proliferation of communication devices. The basic requirements of a system that can provide the aforementioned services are energy efficiency, scalability, lower costs, wide reception range, high localization accuracy and availability. Different technologies such as WiFi, UWB, RFID have been leveraged to provide LBS and Proximity Based Services (PBS), however they do not meet the aforementioned requirements. Apple's Bluetooth Low Energy (BLE) based iBeacon solution primarily intends to provide Proximity Based Services (PBS). However, it suffers from poor proximity detection accuracy due to its reliance on Received Signal Strength Indicator (RSSI) that is prone to multipath fading and drastic fluctuations in the indoor environment. Therefore, in this paper, we present our iBeacon based accurate proximity and indoor localization system. Our two algorithms Server-Side Running Average (SRA) and Server-Side Kalman Filter (SKF) improve the proximity detection accuracy of iBeacons by 29% and 32% respectively, when compared with Apple's current moving average based approach. We also present our novel cascaded Kalman Filter-Particle Filter (KFPF) algorithm for indoor localization. Our cascaded filter approach uses a Kalman Filter (KF) to reduce the RSSI fluctuation and then inputs the filtered RSSI values into a Particle Filter (PF) to improve the accuracy of indoor localization. Our experimental results, obtained through experiments in a space replicating real-world scenario, show that our cascaded filter approach outperforms the use of only PF by 28.16% and 25.59% in 2-Dimensional (2D) and 3-Dimensional (3D) environments respectively, and achieves a localization error as low as 0.70 meters in 2D environment and 0.947 meters in 3D environment.


title: Star Formation Under the Outflow: The Discovery of a Non-Thermal Jet from OMC-2 FIR 3 and its Relationship to the Deeply Embedded FIR 4 Protostar

abstract: We carried out multiwavelength (0.7-5 cm), multiepoch (1994-2015) Very Large Array (VLA) observations toward the region enclosing the bright far-IR sources FIR 3 (HOPS 370) and FIR 4 (HOPS 108) in OMC-2. We report the detection of 10 radio sources, seven of them identified as young stellar objects. We image a well-collimated radio jet with a thermal free-free core (VLA 11) associated with the Class I intermediate-mass protostar HOPS 370. The jet presents several knots (VLA 12N, 12C, 12S) of non-thermal radio emission (likely synchrotron from shock-accelerated relativistic electrons) at distances of ~7,500-12,500 au from the protostar, in a region where other shock tracers have been previously identified. These knots are moving away from the HOPS 370 protostar at ~ 100 km/s. The Class 0 protostar HOPS 108, which itself is detected as an independent, kinematically decoupled radio source, falls in the path of these non-thermal radio knots. These results favor the previously proposed scenario where the formation of HOPS 108 has been triggered by the impact of the HOPS 370 outflow with a dense clump. However, HOPS 108 presents a large proper motion velocity of ~ 30 km/s, similar to that of other runaway stars in Orion, whose origin would be puzzling within this scenario. Alternatively, an apparent proper motion could result because of changes in the position of the centroid of the source due to blending with nearby extended emission, variations in the source shape, and /or opacity effects.


title: Perturbation Theory for BAO reconstructed fields: one-loop results in real-space matter density field

abstract: We compute the power spectrum at one-loop order in standard perturbation theory for the matter density field to which a standard Lagrangian Baryonic acoustic oscillation (BAO) reconstruction technique is applied. The BAO reconstruction method corrects the bulk motion associated with the gravitational evolution using the inverse Zel'dovich approximation (ZA) for the smoothed density field. We find that the overall amplitude of one-loop contributions in the matter power spectrum substantially decrease after reconstruction. The reconstructed power spectrum thereby approaches the initial linear spectrum when the smoothed density field is close enough to linear, i.e., the smoothing scale $R_s$ larger than around 10$h^{-1}$Mpc. On smaller $R_s$,however, the deviation from the linear spectrum becomes significant on large scales ($k\lt R_s^{-1}$) due to the nonlinearity in the smoothed density field, and the reconstruction is inaccurate. Compared with N-body simulations, we show that the reconstructed power spectrum at one loop order agrees with simulations better than the unreconstructed power spectrum. We also calculate the tree-level bispectrum in standard perturbation theory to investigate non-Gaussianity in the reconstructed matter density field. We show that the amplitude of the bispectrum significantly decreases for small $k$ after reconstruction and that the tree-level bispectrum agrees well with N-body results in the weakly nonlinear regime.


title: How to avoid the curse of dimensionality: scalability of particle filters with and without importance weights

abstract: Particle filters are a popular and flexible class of numerical algorithms to solve a large class of nonlinear filtering problems. However, standard particle filters with importance weights have been shown to require a sample size that increases exponentially with the dimension D of the state space in order to achieve a certain performance, which precludes their use in very high-dimensional filtering problems. Here, we focus on the dynamic aspect of this curse of dimensionality (COD) in continuous time filtering, which is caused by the degeneracy of importance weights over time. We show that the degeneracy occurs on a time-scale that decreases with increasing D. In order to soften the effects of weight degeneracy, most particle filters use particle resampling and improved proposal functions for the particle motion. We explain why neither of the two can prevent the COD in general. In order to address this fundamental problem, we investigate an existing filtering algorithm based on optimal feedback control that sidesteps the use of importance weights. We use numerical experiments to show that this Feedback Particle Filter (FPF) by Yang et al. (2013) does not exhibit a COD.


title: Ultrafast Hot Carrier Dynamics in GaN and its Impact on the Efficiency Droop

abstract: GaN is a key material for lighting technology. Yet, the carrier transport and ultrafast dynamics that are central in GaN light emitting devices are not completely understood. We present first-principles calculations of carrier dynamics in GaN, focusing on electron-phonon (e-ph) scattering and the cooling and nanoscale dynamics of hot carriers. We find that e-ph scattering is significantly faster for holes compared to electrons, and that for hot carriers with an initial 0.5$-$1 eV excess energy, holes take a significantly shorter time ($\sim$0.1 ps) to relax to the band edge compared to electrons, which take $\sim$1 ps. The asymmetry in the hot carrier dynamics is shown to originate from the valence band degeneracy, the heavier effective mass of holes compared to electrons, and the details of the coupling to different phonon modes in the valence and conduction bands. We show that the slow cooling of hot electrons and their long ballistic mean free paths (over 3 nm) are a possible cause of efficiency droop in GaN light emitting diodes. Taken together, our work sheds light on the ultrafast dynamics of hot carriers in GaN and the nanoscale origin of efficiency droop.


title: The $K$-theory of toric schemes over regular rings of mixed characteristic

abstract: We show that if $X$ is a toric scheme over a regular commutative ring $k$ then the direct limit of the $K$-groups of $X$ taken over any infinite sequence of nontrivial dilations is homotopy invariant. This theorem was previously known for regular commutative rings containing a field. The affine case of our result was conjectured by Gubeladze. We prove analogous results when $k$ is replaced by an appropriate $K$-regular, not necessarily commutative $k$-algebra.


title: Invited review: Fluctuation-induced transport. From the very small to the very large scales

abstract: The study of fluctuation-induced transport is concerned with the directed motion of particles on a substrate when subjected to a fluctuating external field. Work over the last two decades provides now precise clues on how the average transport depends on three fundamental aspects: the shape of the substrate, the correlations of the fluctuations and the mass, geometry, interaction and density of the particles. These three aspects, reviewed here, acquire additional relevance because the same notions apply to a bewildering variety of problems at very different scales, from the small nano or micro-scale, where thermal fluctuations effects dominate, up to very large scales including ubiquitous cooperative phenomena in granular materials.


title: Analysis of STM images with pure and CO-functionalized tips: A first-principles and experimental study

abstract: We describe a first principles method to calculate scanning tunneling microscopy (STM) images, and compare the results to well-characterized experiments combining STM with atomic force microscopy (AFM). The theory is based on density functional theory (DFT) with a localized basis set, where the wave functions in the vacuum gap are computed by propagating the localized-basis wave functions into the gap using a real-space grid. Constant-height STM images are computed using Bardeen's approximation method, including averaging over the reciprocal space. We consider copper adatoms and single CO molecules adsorbed on Cu(111), scanned with a single-atom copper tip with and without CO functionalization. The calculated images agree with state-of-the-art experiments, where the atomic structure of the tip apex is determined by AFM. The comparison further allows for detailed interpretation of the STM images.


title: Asymptotic $U(1)$ charges at spatial infinity

abstract: Large gauge symmetries in Minkowski spacetime are often studied in two distinct regimes: either at asymptotic (past or future) times or at spatial infinity. By working in harmonic gauge, we provide a unified description of large gauge symmetries (and their associated charges) that applies to both regimes. At spatial infinity the charges are conserved and interpolate between those defined at the asymptotic past and future. This explains the equality of asymptotic past and future charges, as recently proposed in connection with Weinberg's soft photon theorem.


title: Can scalars have asymptotic symmetries?

abstract: Recently it has been understood that certain soft factorization theorems for scattering amplitudes can be written as Ward identities of new asymptotic symmetries. This relationship has been established for soft particles with spins $s > 0$, most notably for soft gravitons and photons. Here we study the remaining case of soft scalars. We show that a class of Yukawa-type theories, where a massless scalar couples to massive particles, have an infinite number of conserved charges. This raises the question as to whether one can associate asymptotic symmetries to scalars.


title: Robust Kronecker-Decomposable Component Analysis for Low-Rank Modeling

abstract: Dictionary learning and component analysis are part of one of the most well-studied and active research fields, at the intersection of signal and image processing, computer vision, and statistical machine learning. In dictionary learning, the current methods of choice are arguably K-SVD and its variants, which learn a dictionary (i.e., a decomposition) for sparse coding via Singular Value Decomposition. In robust component analysis, leading methods derive from Principal Component Pursuit (PCP), which recovers a low-rank matrix from sparse corruptions of unknown magnitude and support. However, K-SVD is sensitive to the presence of noise and outliers in the training set. Additionally, PCP does not provide a dictionary that respects the structure of the data (e.g., images), and requires expensive SVD computations when solved by convex relaxation. In this paper, we introduce a new robust decomposition of images by combining ideas from sparse dictionary learning and PCP. We propose a novel Kronecker-decomposable component analysis which is robust to gross corruption, can be used for low-rank modeling, and leverages separability to solve significantly smaller problems. We design an efficient learning algorithm by drawing links with a restricted form of tensor factorization. The effectiveness of the proposed approach is demonstrated on real-world applications, namely background subtraction and image denoising, by performing a thorough comparison with the current state of the art.


title: Combining Neural Networks and Tree Search for Task and Motion Planning in Challenging Environments

abstract: We consider task and motion planning in complex dynamic environments for problems expressed in terms of a set of Linear Temporal Logic (LTL) constraints, and a reward function. We propose a methodology based on reinforcement learning that employs deep neural networks to learn low-level control policies as well as task-level option policies. A major challenge in this setting, both for neural network approaches and classical planning, is the need to explore future worlds of a complex and interactive environment. To this end, we integrate Monte Carlo Tree Search with hierarchical neural net control policies trained on expressive LTL specifications. This paper investigates the ability of neural networks to learn both LTL constraints and control policies in order to generate task plans in complex environments. We demonstrate our approach in a simulated autonomous driving setting, where a vehicle must drive down a road in traffic, avoid collisions, and navigate an intersection, all while obeying given rules of the road.


title: On the structure of elliptic curves over finite extensions of $\mathbb{Q}_p$ with additive reduction

abstract: Let $p$ be a prime and let $K$ be a finite extension of $\mathbb{Q}_p$. Let $E/K$ be an elliptic curve with additive reduction. In this paper, we study the topological group structure of the set of points of good reduction of $E(K)$. In particular, if $K/\mathbb{Q}_p$ is unramified, we show how one can read off the topological group structure from the Weierstrass coefficients defining $E$.


title: Large-scale asymmetry between spin patterns of spiral galaxies

abstract: Spin patterns of spiral galaxies can be broadly separated into galaxies with clockwise (Z-wise) patterns and galaxies with counterclockwise (S-wise) spin patterns. While the differences between these patterns are visually noticeable, they are a matter of the perspective of the observer, and therefore in a sufficiently large universe no other differences are expected between galaxies with Z-wise and S-wise patterns. Here large datasets of spiral galaxies separated by their spin patterns are used to show that spiral galaxies with Z-wise spin patterns are photometrically different from spiral galaxies with S-wise patterns. That asymmetry changes based on the direction of observation, such that the observed asymmetry in one hemisphere is aligned with the inverse observed asymmetry in the opposite hemisphere. The results are consistent across different sky surveys (SDSS and PanSTARRS) and analysis methods. The proximity of the most probable asymmetry axis to the galactic pole suggests that the asymmetry might be driven by relativistic beaming. Annotated data from SDSS and PanSTARRS are publicly available.


title: User Experience of the CoSTAR System for Instruction of Collaborative Robots

abstract: How can we enable novice users to create effective task plans for collaborative robots? Must there be a tradeoff between generalizability and ease of use? To answer these questions, we conducted a user study with the CoSTAR system, which integrates perception and reasoning into a Behavior Tree-based task plan editor. In our study, we ask novice users to perform simple pick-and-place assembly tasks under varying perception and planning capabilities. Our study shows that users found Behavior Trees to be an effective way of specifying task plans. Furthermore, users were also able to more quickly, effectively, and generally author task plans with the addition of CoSTAR's planning, perception, and reasoning capabilities. Despite these improvements, concepts associated with these capabilities were rated by users as less usable, and our results suggest a direction for further refinement.


title: Width Hierarchies for Quantum and Classical Ordered Binary Decision Diagrams with Repeated Test

abstract: We consider quantum, nondterministic and probabilistic versions of known computational model Ordered Read-$k$-times Branching Programs or Ordered Binary Decision Diagrams with repeated test ($k$-QOBDD, $k$-NOBDD and $k$-POBDD). We show width hierarchy for complexity classes of Boolean function computed by these models and discuss relation between different variants of $k$-OBDD.


title: Random unitaries, amenable linear groups and Jordan's theorem

abstract: It is well known that a dense subgroup $G$ of the complex unitary group $U(d)$ cannot be amenable as a discrete group when $d>1$. When $d$ is large enough we give quantitative versions of this phenomenon in connection with certain estimates of random Fourier series on the compact group $\bar G$ that is the closure of $G$. Roughly, we show that if $\bar G$ covers a large enough part of $U(d)$ in the sense of metric entropy then $G$ cannot be amenable. The results are all based on a version of a classical theorem of Jordan that says that if $G$ is finite, or amenable as a discrete group, then $G$ contains an Abelian subgroup with index $e^{o(d^2)}$.


title: Heavy neutral pseudoscalar decays into gauge bosons in the Littlest Higgs Model

abstract: We study two-body decays of a new neutral pseudoscalar into gauge bosons within the context of the Littlest Higgs model. Concretely, the $\Phi^P \to WW, VV, gg$ processes induced at the one-loop level, with $V=\gamma, Z$, are considered. Since the branching ratios of the $\Phi^P \to VV$ decays result very suppressed, only the $\Phi^P \to WW, gg$ processes are thoroughly studied. The branching ratios for the $\Phi^P \to gg$ and $\Phi^P \to WW$ decays are of the order of $10^{-4}$ and $10^{-6}$, respectively, for $f$ around 2 TeV, which represents the global symmetry breaking scale of the theory. The production cross section of the $\Phi^P$ boson via gluon fusion at LHC is estimated.


title: The cosmic microwave background Cold Spot anomaly: the impact of sky masking and the expected contribution from the Integrated Sachs-Wolfe effect

abstract: We re-analyse the cosmic microwave background (CMB) Cold Spot (CS) anomaly with particular focus on understanding the bias a mask (contaminated by Galactic and point sources) may introduce. We measure the coldest spot, found by applying the Spherical Mexican Hat Wavelet transform on 100 000 cut-sky (masked) and full-sky CMB simulated maps. The CS itself is barely affected by the mask; we estimate a 94 per cent probability that the CS is the full-sky temperature minimum. However, approximately 48 per cent (masked fraction of the mask) of full-sky minima are obscured by the mask. Since the observed minima are slightly hotter than the full-sky ensemble of minima, a cut-sky analysis would have found the CS to be significant at approximately 2.2 sigma with a wavelet angular scale of R = 5 degrees. None the less, comparisons to full-sky minima show the CS significance to be only approximately 1.9 sigma and less than 2 sigma for all R. The CS on the last scattering surface may be hotter due to the integrated Sachs-Wolfe effect in the line of sight. However, our simulations show that this on average only approximately 10 per cent (about 10 micro K but consistent with zero) of the CS temperature profile. This is consistent with Lambda and cold dark matter reconstructions of this effect based on observed line-of-sight voids.


title: Planetesimal formation by the streaming instability in a photoevaporating disk

abstract: Recent years have seen growing interest in the streaming instability as a candidate mechanism to produce planetesimals. However, these investigations have been limited to small-scale simulations. We now present the results of a global protoplanetary disk evolution model that incorporates planetesimal formation by the streaming instability, along with viscous accretion, photoevaporation by EUV, FUV, and X-ray photons, dust evolution, the water ice line, and stratified turbulence. Our simulations produce massive (60-130 $M_\oplus$) planetesimal belts beyond 100 au and up to $\sim 20 M_\oplus$ of planetesimals in the middle regions (3-100 au). Our most comprehensive model forms 8 $M_\oplus$ of planetesimals inside 3 au, where they can give rise to terrestrial planets. The planetesimal mass formed in the inner disk depends critically on the timing of the formation of an inner cavity in the disk by high-energy photons. Our results show that the combination of photoevaporation and the streaming instability are efficient at converting the solid component of protoplanetary disks into planetesimals. Our model, however, does not form enough early planetesimals in the inner and middle regions of the disk to give rise to giant planets and super-Earths with gaseous envelopes. Additional processes such as particle pileups and mass loss driven by MHD winds may be needed to drive the formation of early planetesimal generations in the planet forming regions of protoplanetary disks.


title: Lipschitz-free spaces over compact subsets of superreflexive spaces are weakly sequentially complete

abstract: Let $M$ be a compact subset of a superreflexive Banach space. We prove that the Lipschitz-free space $\mathcal{F}(M)$, the predual of the Banach space of Lipschitz functions on $M$, has the Pe{\l}czy\'nski's property ($V^\ast$). As a consequence, the Lipschitz-free space $\mathcal{F}(M)$ is weakly sequentially complete.


title: Past, Present and Future of UHECR Observations

abstract: Great advances have been made in the study of ultra-high energy cosmic rays (UHECR) in the past two decades. These include the discovery of the spectral cut-off near 5 x 10^19 eV and complex structure at lower energies, as well as increasingly precise information about the composition of cosmic rays as a function of energy. Important improvements in techniques, including extensive surface detector arrays and high resolution air fluorescence detectors, have been instrumental in facilitating this progress. We discuss the status of the field, including the open questions about the nature of spectral structure, systematic issues related to our understanding of composition, and emerging evidence for anisotropy at the highest energies. We review prospects for upgraded and future observatories including Telescope Array, Pierre Auger and JEM-EUSO and other space-based proposals, and discuss promising new technologies based on radio emission from extensive air showers produced by UHECR.


title: Homological mirror symmetry without corrections

abstract: Let $X$ be a closed symplectic manifold equipped a Lagrangian torus fibration over a base $Q$. A construction first considered by Kontsevich and Soibelman produces from this data a rigid analytic space $Y$, which can be considered as a variant of the $T$-dual introduced by Strominger, Yau, and Zaslow. We prove that the Fukaya category of tautologically unobstructed graded Lagrangians in $X$ embeds fully faithfully in the derived category of (twisted) coherent sheaves on $Y$, under the technical assumption that $\pi_2(Q)$ vanishes (all known examples satisfy this assumption). The main new tool is the construction and computation of Floer cohomology groups of Lagrangian fibres equipped with topological infinite rank local systems that correspond, under mirror symmetry, to the affinoid rings introduced by Tate, equipped with their natural topologies as Banach algebras.


title: On the proof of the Thin Sandwich Conjecture in arbitrary dimensions

abstract: In this paper we show the validity, under certain geometric conditions, of Wheeler's thin sandwich conjecture for higher dimensional theories of gravity. We extend the results shown by R. Bartnik and G. Fodor for the 3-dimensional case in two ways. On the one hand, we show that the results obtained by the mentioned authors are valid in arbitrary dimensions, and on the other hand we show that the geometric hypotheses needed for the proofs can always be satisfied, which constitutes in itself a new result for the 3-dimensional case. In this way, we show that on any compact n-dimensional manifold, n greater or equal to 3, there is an open set in the space of all possible initial data where the thin sandwich problem is well-posed.


title: Conditional stable soliton resolution for a semi-linear Skyrme equation

abstract: We study a semi-linear version of the Skyrme system due to Adkins and Nappi. The objects in this system are maps from $(1+3)$-dimensional Minkowski space into the $3$-sphere and 1-forms on $\mathbb{R}^{1+3}$, coupled via a Lagrangian action. Under a co-rotational symmetry reduction we establish the existence, uniqueness, and unconditional asymptotic stability of a family of stationary solutions $Q_n$, indexed by the topological degree $n \in \mathbb{N} \cup \{0\}$ of the underlying map. We also prove that an arbitrarily large equivariant perturbation of $Q_n$ leads to a globally defined solution that scatters to $Q_n$ in infinite time as long as the critical norm for the solution remains bounded on the maximal interval of existence given by the local Cauchy theory. We remark that the evolution equations are super-critical with respect to the conserved energy.


title: The SIC Question: History and State of Play

abstract: Recent years have seen significant advances in the study of symmetric informationally complete (SIC) quantum measurements, also known as maximal sets of complex equiangular lines. Previously, the published record contained solutions up to dimension 67, and was with high confidence complete up through dimension 50. Computer calculations have now furnished solutions in all dimensions up to 151, and in several cases beyond that, as large as dimension 844. These new solutions exhibit an additional type of symmetry beyond the basic definition of a SIC, and so verify a conjecture of Zauner in many new cases. The solutions in dimensions 68 through 121 were obtained by Andrew Scott, and his catalogue of distinct solutions is, with high confidence, complete up to dimension 90. Additional results in dimensions 122 through 151 were calculated by the authors using Scott's code. We recap the history of the problem, outline how the numerical searches were done, and pose some conjectures on how the search technique could be improved. In order to facilitate communication across disciplinary boundaries, we also present a comprehensive bibliography of SIC research.


title: Nonlinear damped wave equations for the sub-Laplacian on the Heisenberg group and for Rockland operators on graded Lie groups

abstract: In this paper we study the Cauchy problem for the semilinear damped wave equation for the sub-Laplacian on the Heisenberg group. In the case of the positive mass, we show the global in time well-posedness for small data for power like nonlinearities. We also obtain similar well-posedness results for the wave equations for Rockland operators on general graded Lie groups. In particular, this includes higher order operators on $\mathbb R^n$ and on the Heisenberg group, such as powers of the Laplacian or of the sub-Laplacian. In addition, we establish a new family of Gagliardo-Nirenberg inequalities on graded Lie groups that play a crucial role in the proof but which are also of interest on their own: if $G$ is a graded Lie group of homogeneous dimension $Q$ and $a>0$, $1<r<\frac{Q}{a},$ and $1\leq p\leq q\leq \frac{rQ}{Q-ar},$ then we have the following Gagliardo-Nirenberg type inequality $$ \|u\|_{L^{q}(G)}\lesssim \|u\|_{\dot{L}_{a}^{r}(G)}^{s} \|u\|_{L^{p}(G)}^{1-s} $$ for $s=\left(\frac1p-\frac1q\right) \left(\frac{a}Q+\frac1p-\frac1r\right)^{-1}\in [0,1]$ provided that $\frac{a}Q+\frac1p-\frac1r\not=0$, where $\dot{L}_{a}^{r}$ is the homogeneous Sobolev space of order $a$ over $L^r$. If $\frac{a}Q+\frac1p-\frac1r=0$, we have $p=q=\frac{rQ}{Q-ar}$, and then the above inequality holds for any $0\leq s\leq 1$.


title: Central limit theorem for Fourier transform and periodogram of random fields

abstract: In this paper we show that the limiting distribution of the real and the imaginary part of the double Fourier transform of a stationary random field is almost surely an independent vector with Gaussian marginal distributions, whose variance is, up to a constant, the field's spectral density. The dependence structure of the random field is general and we do not impose any restrictions on the speed of convergence to zero of the covariances, or smoothness of the spectral density. The only condition required is that the variables are adapted to a commuting filtration and are regular in some sense. The results go beyond the Bernoulli fields and apply to both short range and long range dependence. They can be easily applied to derive the asymptotic behavior of the periodogram associated to the random field. The method of proof is based on new probabilistic methods involving martingale approximations and also on borrowed and new tools from harmonic analysis.


title: Cross-Validation with Confidence

abstract: Cross-validation is one of the most popular model selection methods in statistics and machine learning. Despite its wide applicability, traditional cross validation methods tend to select overfitting models, due to the ignorance of the uncertainty in the testing sample. We develop a new, statistically principled inference tool based on cross-validation that takes into account the uncertainty in the testing sample. This new method outputs a set of highly competitive candidate models containing the best one with guaranteed probability. As a consequence, our method can achieve consistent variable selection in a classical linear regression setting, for which existing cross-validation methods require unconventional split ratios. When used for regularizing tuning parameter selection, the method can provide a further trade-off between prediction accuracy and model interpretability. We demonstrate the performance of the proposed method in several simulated and real data examples.


title: A finite simple group is CCA if and only if it has no element of order four

abstract: A Cayley graph for a group $G$ is CCA if every automorphism of the graph that preserves the edge-orbits under the regular representation of $G$ is an element of the normaliser of $G$. A group $G$ is then said to be CCA if every connected Cayley graph on $G$ is CCA. We show that a finite simple group is CCA if and only if it has no element of order 4. We also show that "many" 2-groups are non-CCA.


title: Decomposition theory of modules: the case of Kronecker algebra

abstract: Let $A$ be a finite-dimensional algebra over an algebraically closed field $\Bbbk$. For any finite-dimensional $A$-module $M$ we give a general formula that computes the indecomposable decomposition of $M$ without decomposing it, for which we use the knowledge of AR-quivers that are already computed in many cases. The proof of the formula here is much simpler than that in a prior literature by Dowbor and Mr\'oz. As an example we apply this formula to the Kronecker algebra $A$ and give an explicit formula to compute the indecomposable decomposition of $M$, which enables us to make a computer program.


title: Robust Polynomial Reconstruction via Chinese Remainder Theorem in the Presence of Small Degree Residue Errors

abstract: Based on unique decoding of the polynomial residue code with non-pairwise coprime moduli, a polynomial with degree less than that of the least common multiple (lcm) of all the moduli can be accurately reconstructed when the number of residue errors is less than half the minimum distance of the code. However, once the number of residue errors is beyond half the minimum distance of the code, the unique decoding may fail and lead to a large reconstruction error. In this paper, assuming that all the residues are allowed to have errors with small degrees, we consider how to reconstruct the polynomial as accurately as possible in the sense that a reconstructed polynomial is obtained with only the last $\tau$ number of coefficients being possibly erroneous, when the residues are affected by errors with degrees upper bounded by $\tau$. In this regard, we first propose a multi-level robust Chinese remainder theorem (CRT) for polynomials, namely, a trade-off between the dynamic range of the degree of the polynomial to be reconstructed and the residue error bound $\tau$ is formulated. Furthermore, a simple closed-form reconstruction algorithm is also proposed.


title: Minimum Degree-Weighted Distance Decoding for Polynomial Residue Codes with Non-Pairwise Coprime Moduli

abstract: This paper presents a new decoding for polynomial residue codes, called the minimum degree-weighted distance decoding. The newly proposed decoding is based on the degree-weighted distance and different from the traditional minimum Hamming distance decoding. It is shown that for the two types of minimum distance decoders, i.e., the minimum degree-weighted distance decoding and the minimum Hamming distance decoding, one is not absolutely stronger than the other, but they can complement each other from different points of view.


title: Data Driven Exploratory Attacks on Black Box Classifiers in Adversarial Domains

abstract: While modern day web applications aim to create impact at the civilization level, they have become vulnerable to adversarial activity, where the next cyber-attack can take any shape and can originate from anywhere. The increasing scale and sophistication of attacks, has prompted the need for a data driven solution, with machine learning forming the core of many cybersecurity systems. Machine learning was not designed with security in mind, and the essential assumption of stationarity, requiring that the training and testing data follow similar distributions, is violated in an adversarial domain. In this paper, an adversary's view point of a classification based system, is presented. Based on a formal adversarial model, the Seed-Explore-Exploit framework is presented, for simulating the generation of data driven and reverse engineering attacks on classifiers. Experimental evaluation, on 10 real world datasets and using the Google Cloud Prediction Platform, demonstrates the innate vulnerability of classifiers and the ease with which evasion can be carried out, without any explicit information about the classifier type, the training data or the application domain. The proposed framework, algorithms and empirical evaluation, serve as a white hat analysis of the vulnerabilities, and aim to foster the development of secure machine learning frameworks.


title: Bidirectional-Convolutional LSTM Based Spectral-Spatial Feature Learning for Hyperspectral Image Classification

abstract: This paper proposes a novel deep learning framework named bidirectional-convolutional long short term memory (Bi-CLSTM) network to automatically learn the spectral-spatial feature from hyperspectral images (HSIs). In the network, the issue of spectral feature extraction is considered as a sequence learning problem, and a recurrent connection operator across the spectral domain is used to address it. Meanwhile, inspired from the widely used convolutional neural network (CNN), a convolution operator across the spatial domain is incorporated into the network to extract the spatial feature. Besides, to sufficiently capture the spectral information, a bidirectional recurrent connection is proposed. In the classification phase, the learned features are concatenated into a vector and fed to a softmax classifier via a fully-connected operator. To validate the effectiveness of the proposed Bi-CLSTM framework, we compare it with several state-of-the-art methods, including the CNN framework, on three widely used HSIs. The obtained results show that Bi-CLSTM can improve the classification performance as compared to other methods.


title: Variation of entanglement entropy and mutual information in fermion-fermion scattering

abstract: We study the behavior of entanglement between different degrees of freedom of scattering fermions, based on an exemplary QED scattering process $e^+e^-\longrightarrow\mu^+\mu^-$. The variation of entanglement entropy between two fermions from an initial state to the final state was computed, with respect to different entanglement between the ingoing particles. This variation of entanglement entropy is found to be proportional to an area quantity, the total cross section. We also study the spin-momentum and helicity-momentum entanglements within one particle in the aforementioned scattering process. The calculations of the relevant variations of mutual information in the same inertial frame reveals that, for a maximally entangled initial state, the scattering between the particles does not affect the degree of both of these entanglements of one particle in the final state. It is also found that the increasing degree of entanglement between two ingoing particles would restrict the generation of entanglement between spin (helicity) and momentum of one outgoing particle. And the entanglement between spin and momentum within one particle in the final state is shown to always be stronger than that for helicity-momentum for a general initial entanglement state, implying significantly distinct properties of entanglement for the helicity and spin perceived by an inertial observer.


title: Toward Traffic Patterns in High-speed Railway Communication Systems: Power Allocation and Antenna Selection

abstract: In high-speed railway (HSR) communication systems, distributed antenna is usually employed to support frequent handover and enhance the signal to noise ratio to user equipments. In this case, dynamic time-domain power allocation and antenna selection (PAWAS) could be jointly optimized to improve the system performances. This paper consider this problem in such a simple way where dynamic switching between multiple-input-multiple-output (MIMO) and single-input-multiple-output (SIMO) is allowed and exclusively utilized, while the channel states and traffic demand are taken into account. The channel states includes sparse and rich scattering terrains, and the traffic patterns includes delay-sensitive and delay-insensitive as well as hybrid. Some important results are obtained in theory. In sparse scattering terrains, for delay-sensitive traffic, the PAWAS can be viewed as the generalization of channel-inversion associated with transmit antenna selection. On the contrary, for delay-insensitive traffic, the power allocation with MIMO can be viewed as channel-inversion, but with SIMO, it is traditional water-filling. For the hybrid traffic, the PAWAS can be partitioned as delay-sensitive and delay-insensitive parts by some specific strategies. In rich scattering terrains, the corresponding PAWAS is derived by some amendments in sparse scattering terrains and similar results are then presented.


title: Molecular Hydrogen: a benchmark system for Near Threshold Resonances at higher partial waves

abstract: Benchmark reactions involving molecular hydrogen, such as H$_2$+D or H$_2$+Cl, provide the ideal platforms to investigate the effect of Near Threshold Resonances (NTR) on scattering processes. Due to the small reduced mass of those systems, shape resonances due to particular partial waves can provide features at scattering energies up to a few Kelvins, reachable in recent experiments. We explore the effect of NTRs on elastic and inelastic scattering for higher partial waves $\ell$ in the case of H$_2$+Cl for $s$-wave and H$_2$+D for $p$-wave scattering, and find that NTRs lead to a different energy scaling of the cross sections as compared to the well known Wigner threshold regime. We give a theoretical analysis based on Jost functions for short range interaction potentials. To explore higher partial waves, we adopt a three channel model that incorporates all key ingredients, and explore how the NTR scaling is affected by $\ell$. A short discussion on the effect of the long-range form of the interaction potential is also provided.


title: Why do similarity matching objectives lead to Hebbian/anti-Hebbian networks?

abstract: Modeling self-organization of neural networks for unsupervised learning using Hebbian and anti-Hebbian plasticity has a long history in neuroscience. Yet, derivations of single-layer networks with such local learning rules from principled optimization objectives became possible only recently, with the introduction of similarity matching objectives. What explains the success of similarity matching objectives in deriving neural networks with local learning rules? Here, using dimensionality reduction as an example, we introduce several variable substitutions that illuminate the success of similarity matching. We show that the full network objective may be optimized separately for each synapse using local learning rules both in the offline and online settings. We formalize the long-standing intuition of the rivalry between Hebbian and anti-Hebbian rules by formulating a min-max optimization problem. We introduce a novel dimensionality reduction objective using fractional matrix exponents. To illustrate the generality of our approach, we apply it to a novel formulation of dimensionality reduction combined with whitening. We confirm numerically that the networks with learning rules derived from principled objectives perform better than those with heuristic learning rules.


title: Perspective: Energy Landscapes for Machine Learning

abstract: Machine learning techniques are being increasingly used as flexible non-linear fitting and prediction tools in the physical sciences. Fitting functions that exhibit multiple solutions as local minima can be analysed in terms of the corresponding machine learning landscape. Methods to explore and visualise molecular potential energy landscapes can be applied to these machine learning landscapes to gain new insight into the solution space involved in training and the nature of the corresponding predictions. In particular, we can define quantities analogous to molecular structure, thermodynamics, and kinetics, and relate these emergent properties to the structure of the underlying landscape. This Perspective aims to describe these analogies with examples from recent applications, and suggest avenues for new interdisciplinary research.

categories: stat.ML cond-mat.dis-nn cs.CV cs.LG hep-th


title: MoS2-graphene in-plane contact for high interfacial thermal conduction

abstract: Recent studies showed that the in-plane and inter-plane thermal conductivities of two-dimensional (2D) MoS2 are low, posing a significant challenge in heat management in MoS2-based electronic devices. To address this challenge, we design the interfaces between MoS2 and graphene by fully utilizing graphene, a 2D material with an ultra-high thermal conduction. We first perform ab initio atomistic simulations to understand the bonding nature and structure stability of the interfaces. Our results show that the designed interfaces, which are found to be connected together by strong covalent bonds between Mo and C atoms, are energetically stable. We then perform molecular dynamics simulations to investigate the interfacial thermal conductance. It is found surprisingly that the interface thermal conductance is high, comparable to that of graphene-metal covalent-bonded interfaces. Importantly, each interfacial Mo-C bond serves as an independent thermal channel, enabling the modulation of interfacial thermal conductance by controlling Mo vacancy concentration at the interface. The present work provides a viable route for heat management in MoS2 based electronic devices.


title: Distribution of captured planetesimals in circumplanetary gas disks and implications for accretion of regular satellites

abstract: Regular satellites of giant planets are formed by accretion of solid bodies in circumplanetary disks. Planetesimals that are moving on heliocentric orbits and are sufficiently large to be decoupled from the flow of the protoplanetary gas disk can be captured by gas drag from the circumplanetary disk. In the present work, we examine the distribution of captured planetesimals in circumplanetary disks using orbital integrations. We find that the number of captured planetesimals reaches an equilibrium state as a balance between continuous capture and orbital decay into the planet. The number of planetesimals captured into retrograde orbits is much smaller than those on prograde orbits, because the former ones experience strong headwind and spiral into the planet rapidly. We find that the surface number density of planetesimals at the current radial location of regular satellites can be significantly enhanced by gas drag capture, depending on the velocity dispersions of planetesimals and the width of the gap in the protoplanetary disk. Using a simple model, we also examine the ratio of the surface densities of dust and captured planetesimals in the circumplanetary disk, and find that solid material at the current location of regular satellites can be dominated by captured planetesimals when the velocity dispersion of planetesimals is rather small and a wide gap is not formed in the protoplanetary disk. In this case, captured planetesimals in such a region can grow by mutual collision before spiraling into the planet, and would contribute to the growth of regular satellites.


title: Boson-vortex duality in compressible spin-orbit coupled BECs

abstract: Using a (1+2)-dimensional boson-vortex duality between non-linear electrodynamics and a two-component compressible Bose-Einstein condensate (BEC) with spin-orbit (SO) coupling, we obtain generalised versions of the hydrodynamic continuity and Euler equations where the phase defect and non-defect degrees of freedom enter separately. We obtain the generalised Magnus force on vortices under SO coupling, and associate the linear confinement of vortices due to SO coupling with instanton fluctuations of the dual theory.


title: Glass-based charged particle detector performance for Horizon-T EAS detector system

abstract: An implementation of a novel of glass-based detector with fast response and wide detection range is needed to increase resolution for ultra-high energy cosmic rays detection. Such detector has been designed and built for the Horizon-T detector system at Tien Shan high-altitude Science Station. The main characteristics, such as design, duration of the detector pulse and calibration of a single particle response are discussed.


title: Changing Fashion Cultures

abstract: The paper presents a novel concept that analyzes and visualizes worldwide fashion trends. Our goal is to reveal cutting-edge fashion trends without displaying an ordinary fashion style. To achieve the fashion-based analysis, we created a new fashion culture database (FCDB), which consists of 76 million geo-tagged images in 16 cosmopolitan cities. By grasping a fashion trend of mixed fashion styles,the paper also proposes an unsupervised fashion trend descriptor (FTD) using a fashion descriptor, a codeword vetor, and temporal analysis. To unveil fashion trends in the FCDB, the temporal analysis in FTD effectively emphasizes consecutive features between two different times. In experiments, we clearly show the analysis of fashion trends and fashion-based city similarity. As the result of large-scale data collection and an unsupervised analyzer, the proposed approach achieves world-level fashion visualization in a time series. The code, model, and FCDB will be publicly available after the construction of the project page.


title: Null Geodesics and Red-Blue Shifts of Photons Emitted from Geodesic Particles Around a Non-Commutative Black Hole Spacetime

abstract: We investigate the geodesic motion of massless test particles in the background of a noncommutative geometry inspired Schwarzschild black hole. The behaviour of effective potential is analysed in the equatorial plane and the possible motions of massless particles (i.e. photons) for different values of impact parameter are discussed accordingly. We have also calculated the frequency shift of photons in this spacetime. Further, the mass parameter of a non-commutative inspired Schwarzschild black hole is computed in terms of the measurable redshift of photons emitted by massive particles moving along circular geodesics in equatorial plane. It is observed that the the gravitational field of a non-commutative inspired Schwarzschild black hole is more attractive than the Schwarzschild black hole in General Relativity.


title: Local action of curl-less vector potential on vortex dynamics in superconductors

abstract: Study of the Abrikosov vortex motion in superconductors based on time-dependent Ginzburg-Landau equations reveals an opportunity to locally detect the values of the Aharonov-Bohm type curl-less vector potentials.


title: Constraint on cosmological parameters by Hubble parameter from gravitational wave standard sirens of neutron star binary system

abstract: In this paper, we present a new method of measuring Hubble parameter($H(z)$), making use of the anisotropy of luminosity distance($d_{L}$), and the analysis of gravitational wave(GW) of neutron star(NS) binary system. The method has never been put into practice before due to the lack of the ability of detecting GW. LIGO's success in detecting GW of black hole(BH) binary system merger announced the possibility of this new method. We apply this method to several GW detecting projects, including Advanced LIGO(Adv-LIGO), Einstein Telescope(ET) and DECIGO, finding that the $H(z)$ by Adv-LIGO and ET is of bad accuracy, while the $H(z)$ by DECIGO shows a good accuracy. We use the error information of $H(z)$ by DECIGO to simulate $H(z)$ data at every 0.1 redshift span, and put the mock data into the forecasting of cosmological parameters. Compared with the available 38 observed $H(z)$ data(OHD), mock data shows an obviously tighter constraint on cosmological parameters, and a concomitantly higher value of Figure of Merit(FoM). For a 3-year-observation by standard sirens of DECIGO, the FoM value is as high as 834.9. If a 10-year-observation is launched, the FoM could reach 2783.1. For comparison, the FoM of 38 actual observed $H(z)$ data is 9.3. These improvement indicates that the new method has great potential in further cosmological constraints.


title: Symmetric functions and wavefunctions of the six-vertex model by Izergin-Korepin analysis

abstract: We analyze wavefunctions of the six-vertex model by extending the Izergin-Korepin analysis on the domain wall boundary partition functions. We particularly focus on the case with triangular boundary. By using the $U_q(sl_2)$ $R$-matrix and a special class of the triangular $K$-matrix, we first introduce an analogue of the wavefunctions of the integrable six-vertex model with triangular boundary. We first give a characterization of the wavefunctions by extending our recent work of the Izergin-Korepin analysis of the domain wall boundary partition function with triangular boundary, and then determine the explicit form of the symmetric functions representing the wavefunctions by showing that it satisfies all the required properties. We also illustrate the Izergin-Korepin analysis for the case of ordinary wavefunctions as it is the basic case.


title: On integrability aspects of the supersymmetric sine-Gordon equation

abstract: In this paper we study certain integrability properties of the supersymmetric sine-Gordon equation. We construct Lax pairs with their zero-curvature representations which are equivalent to the supersymmetric sine-Gordon equation. From the fermionic linear spectral problem, we derive coupled sets of super Riccati equations and the auto-B\"acklund transformation of the supersymmetric sine-Gordon equation. In addition, a detailed description of the associated Darboux transformation is presented and non-trivial super multisoliton solutions are constructed. These integrability properties allow us to provide new explicit geometric characterizations of the bosonic supersymmetric version of the Sym--Tafel formula for the immersion of surfaces in a Lie superalgebra. These characterizations are expressed only in terms of the independent bosonic and fermionic variables.


title: Two-dimensional grating magneto-optical trap

abstract: We demonstrate a two-dimensional grating magneto-optical trap (2D GMOT) with a single input cooling laser beam and a planar diffraction grating using $^{87}$Rb. This configuration increases experimental access when compared with a traditional 2D MOT. As described in the paper, the output flux is several hundred million rubidium atoms/s at a mean velocity of $16.5(9)$ m/s and a velocity distribution of $4(3)$ m/s standard deviation. We use the atomic beam from the 2D GMOT to demonstrate loading of a three dimensional grating MOT (3D GMOT) with $2.46(7)\times 10^8$ atoms. Methods to improve output flux are discussed.


title: Collective excitations of self-bound droplets of a dipolar quantum fluid

abstract: We calculate the collective excitations of a dipolar Bose-Einstein condensate in the regime where it self-binds into droplets stabilized by quantum fluctuations. We show that the filament-shaped droplets act as a quasi-one-dimensional waveguide along which low angular momentum phonons propagate. The evaporation (unbinding) threshold occurring as the atom number $N$ is reduced to the critical value $N_c$ is associated with a monopole-like excitation going soft as $\epsilon_0\sim(N-N_c)^{1/4}$. Considering the system in the presence of a trapping potential, we quantify the crossover from a trap-bound condensate to a self-bound droplet.


title: Self corrective Perturbations for Semantic Segmentation and Classification

abstract: Convolutional Neural Networks have been a subject of great importance over the past decade and great strides have been made in their utility for producing state of the art performance in many computer vision problems. However, the behavior of deep networks is yet to be fully understood and is still an active area of research. In this work, we present an intriguing behavior: pre-trained CNNs can be made to improve their predictions by structurally perturbing the input. We observe that these perturbations - referred as Guided Perturbations - enable a trained network to improve its prediction performance without any learning or change in network weights. We perform various ablative experiments to understand how these perturbations affect the local context and feature representations. Furthermore, we demonstrate that this idea can improve performance of several existing approaches on semantic segmentation and scene labeling tasks on the PASCAL VOC dataset and supervised classification tasks on MNIST and CIFAR10 datasets.


title: Diversification-Based Learning in Computing and Optimization

abstract: Diversification-Based Learning (DBL) derives from a collection of principles and methods introduced in the field of metaheuristics that have broad applications in computing and optimization. We show that the DBL framework goes significantly beyond that of the more recent Opposition-based learning (OBL) framework introduced in Tizhoosh (2005), which has become the focus of numerous research initiatives in machine learning and metaheuristic optimization. We unify and extend earlier proposals in metaheuristic search (Glover, 1997, Glover and Laguna, 1997) to give a collection of approaches that are more flexible and comprehensive than OBL for creating intensification and diversification strategies in metaheuristic search. We also describe potential applications of DBL to various subfields of machine learning and optimization.


title: Polynomial expressions of $p$-ary auction functions

abstract: Let $\mathbb{F}_p$ be the finite field of prime order $p$. For any function $f \colon \mathbb{F}_p{}^n \to \mathbb{F}_p$, there exists a unique polynomial over $\mathbb{F}_p$ having degree at most $p-1$ with respect to each variable which coincides with $f$. We call it the minimal polynomial of $f$. It is in general a non-trivial task to find a concrete expression of the minimal polynomial of a given function, which has only been worked out for limited classes of functions in the literature. In this paper, we study minimal polynomial expressions of several functions that are closely related to some practically important procedures such as auction and voting.


title: Vertical Structure of the Transition Zone from Infalling Rotating Envelope to Disk in the Class 0 Protostar, IRAS04368+2557

abstract: We have resolved for the first time the radial and vertical structure of the almost edge-on envelope/disk system of the low-mass Class 0 protostar L1527. For that, we have used ALMA observations with a spatial resolution of 0.25$^{\prime\prime}$$\times$0.13$^{\prime\prime}$ and 0.37$^{\prime\prime}$$\times$0.23$^{\prime\prime}$ at 0.8 mm and 1.2 mm, respectively. The L1527 dust continuum emission has a deconvolved size of 78 au $\times$ 21 au, and shows a flared disk-like structure. A thin infalling-rotating envelope is seen in the CCH emission outward of about 150 au, and its thickness is increased by a factor of 2 inward of it. This radius lies between the centrifugal radius (200 au) and the centrifugal barrier of the infalling-rotating envelope (100 au). The gas stagnates in front of the centrifugal barrier and moves toward vertical directions. SO emission is concentrated around and inside the centrifugal barrier. The rotation speed of the SO emitting gas is found to be decelerated around the centrifugal barrier. A part of the angular momentum could be extracted by the gas which moves away from the mid-plane around the centrifugal barrier. If this is the case, the centrifugal barrier would be related to the launching mechanism of low velocity outflows, such as disk winds.


title: Family-Personalized Dietary Planning with Temporal Dynamics

abstract: Poor diet and nutrition in the United States has immense financial and health costs, and development of new tools for diet planning could help families better balance their financial and temporal constraints with the quality of their diet and meals. This paper formulates a novel model for dietary planning that incorporates two types of temporal constraints (i.e., dynamics on the perishability of raw ingredients over time, and constraints on the time required to prepare meals) by explicitly incorporating the relationship between raw ingredients and selected food recipes. Our formulation is a diet planning model with integer-valued decision variables, and so we study the problem of designing approximation algorithms (i.e, algorithms with polynomial-time computation and guarantees on the quality of the computed solution) for our dietary model. We develop a deterministic approximation algorithm that is based on a deterministic variant of randomized rounding, and then evaluate our deterministic approximation algorithm with numerical experiments of dietary planning using a database of about 2000 food recipes and 150 raw ingredients.


title: Speeding up adiabatic state conversion in optomechanical system

abstract: Speeding up adiabatic method has attracted much attention with the wide applications in quantum information processing. In this paper, two kinds of methods, Lewis-Riesenfeld invariant-based inverse engineering and transitionless quantum driving are applied to implement speeding up adiabatic state conversion in optomechanical system. The perfect population transfer can be achieved within a short time. At last, the energetic cost is analysed for the transitionless quantum driving.


title: Integral control of port-Hamiltonian systems: non-passive outputs without coordinate transformation

abstract: In this paper we present a method for the addition of integral action to non-passive outputs of a class of port-Hamiltonian systems. The proposed integral controller is a dynamic extension, constructed from the open loop system, such that the closed loop preserves the port-Hamiltonian form. It is shown that the controller is able to reject the effects of both matched and unmatched disturbances, preserving the regulation of the non-passive outputs. Previous solutions to this problem have relied on a change of coordinates whereas the presented solution is developed using the original state vector and, therefore, retains its physical interpretation. In addition, the resulting closed loop dynamics have a natural interpretation as a Control by Interconnection scheme.


title: New method for torque magnetometry using a commercially available membrane-type surface-stress sensor

abstract: We present a new method for torque magnetometry by using a commercially available membrane-type surface-stress sensor (MSS). This sensor has a silicon membrane supported by four beams in which piezoresistive paths are integrated. Although originally developed as a gas sensor, it can be used for torque measurement by modifying its on-chip wiring. We demonstrate the magnetic-torque measurement of submillimeter-sized crystals at a low temperature and in strong magnetic fields. This MSS can observe de-Haas-van-Alphen oscillation, which confirms that it can be an alternative tool for self-sensitive microcantilevers.


title: Synchronization of coupled metronomes on two layers

abstract: Coupled metronomes serve as a paradigmatic model for exploring the collective behaviors of complex dynamical systems, as well as a classical setup for classroom demonstrations of synchronization phenomena. Whereas previous studies of metronome synchronization have been concentrating on symmetric coupling schemes, here we consider the asymmetric case by adopting the scheme of layered metronomes. Specifically, we place two metronomes on each layer, and couple two layers by placing one on top of the other. By varying the initial conditions of the metronomes and adjusting the friction between the two layers, a variety of synchronous patterns are observed in experiment, including the splay synchronization (SS) state, the generalized splay synchronization (GSS) state , the anti-phase synchronization (APS) state, the in-phase delay synchronization (IPDS) state, and the in-phase synchronization (IPS) state. In particular, the IPDS state, in which the metronomes on each layer are synchronized in phase but are of a constant phase delay to metronomes on the other layer, is observed for the first time. In addition, a new technique based on audio signals is proposed for pattern detection, which is more convenient and easier to apply than the existing acquisition techniques. Furthermore, a theoretical model is developed to explain the experimental observations, and is employed to explore the dynamical properties of the patterns, including the basin distributions and the pattern transitions. Our study sheds new lights on the collective behaviors of coupled metronomes, and the developed setup can be used in the classroom for demonstration purposes.


title: Spectral Analysis of Piezoelectric Tensors

abstract: A third order real tensor is called a piezoelectric-type tensor if it is partially symmetric with respect to its last two indices. The piezoelectric tensor is a piezoelectric-type tensor of dimension three. We introduce C-eigenvalues and C-eigenvectors for piezoelectric-type tensors. Here, "C" names after Curie brothers, who first discovered the piezoelectric effect. We show that C-eigenvalues always exist, they are invariant under orthonormal transformations, and for a piezoelectric-type tensor, the largest C-eigenvalue and its C-eigenvectors form the best rank-one piezoelectric-type approximation of that tensor. This means that for the piezoelectric tensor, its largest C-eigenvalue determines the highest piezoelectric coupling constant. We further show that for the piezoelectric tensor, the largest C-eigenvalue corresponds to the electric displacement vector with the largest $2$-norm in the piezoelectric effect under unit uniaxial stress, and the strain tensor with the largest $2$-norm in the converse piezoelectric effect under unit electric field vector. Thus, C-eigenvalues and C-eigenvectors have concrete physical meanings in piezoelectric effect and converse piezoelectric effect. Finally, we apply C-eigenvalues and associated C-eigenvectors for various piezoelectric crystals with different symmetries.


title: Planar Object Tracking in the Wild: A Benchmark

abstract: Planar object tracking is an actively studied problem in vision-based robotic applications. While several benchmarks have been constructed for evaluating state-of-the-art algorithms, there is a lack of video sequences captured in the wild rather than in constrained laboratory environment. In this paper, we present a carefully designed planar object tracking benchmark containing 210 videos of 30 planar objects sampled in the natural environment. In particular, for each object, we shoot seven videos involving various challenging factors, namely scale change, rotation, perspective distortion, motion blur, occlusion, out-of-view, and unconstrained. The ground truth is carefully annotated semi-manually to ensure the quality. Moreover, eleven state-of-the-art algorithms are evaluated on the benchmark using two evaluation metrics, with detailed analysis provided for the evaluation results. We expect the proposed benchmark to benefit future studies on planar object tracking.


title: Recurrent Multimodal Interaction for Referring Image Segmentation

abstract: In this paper we are interested in the problem of image segmentation given natural language descriptions, i.e. referring expressions. Existing works tackle this problem by first modeling images and sentences independently and then segment images by combining these two types of representations. We argue that learning word-to-image interaction is more native in the sense of jointly modeling two modalities for the image segmentation task, and we propose convolutional multimodal LSTM to encode the sequential interactions between individual words, visual information, and spatial information. We show that our proposed model outperforms the baseline model on benchmark datasets. In addition, we analyze the intermediate output of the proposed multimodal LSTM approach and empirically explain how this approach enforces a more effective word-to-image interaction.


title: Unsupervised Basis Function Adaptation for Reinforcement Learning

abstract: When using reinforcement learning (RL) algorithms it is common, given a large state space, to introduce some form of approximation architecture for the value function (VF). The exact form of this architecture can have a significant effect on an agent's performance, however, and determining a suitable approximation architecture can often be a highly complex task. Consequently there is currently interest among researchers in the potential for allowing RL algorithms to adaptively generate (i.e. to learn) approximation architectures. One relatively unexplored method of adapting approximation architectures involves using feedback regarding the frequency with which an agent has visited certain states to guide which areas of the state space to approximate with greater detail. In this article we will: (a) informally discuss the potential advantages offered by such methods; (b) introduce a new algorithm based on such methods which adapts a state aggregation approximation architecture on-line and is designed for use in conjunction with SARSA; (c) provide theoretical results, in a policy evaluation setting, regarding this particular algorithm's complexity, convergence properties and potential to reduce VF error; and finally (d) test experimentally the extent to which this algorithm can improve performance given a number of different test problems. Taken together our results suggest that our algorithm (and potentially such methods more generally) can provide a versatile and computationally lightweight means of significantly boosting RL performance given suitable conditions which are commonly encountered in practice.


title: The Two Higgs-Doublet Extension of the Standard Model and S3 Symmetry

abstract: (This is a recompiled version of $Chin. J. of Phys., Vol. 26, No. 4, 180$ since the previous available electronic document on web was unreadable) We impose a permutation symmetry S3 upon the two Higgs doublet extension of the SU(2) x U( 1) standard model, and then break S3 symmetry spontaneously. A mass relation is obtained from which the mass of the top quark is predicted to be 39 Gev at 1 Gev renormalization scale, if the relative strength between the coupling of the first Higgs field with the up-type quarks is assumed to be the same as that with the down-type quarks. We find that in our model the flavor changing neutral currents vanish automatically.


title: Complex Balancing Reconstructed to the Asymptotic Stability of Mass-action Chemical Reaction Networks with Conservation Laws

abstract: Motivated by the fact that the pseudo-Helmholtz function is a valid Lyapunov function for characterizing asymptotic stability of complex balanced mass action systems (MASs), this paper develops the generalized pseudo-Helmholtz function for stability analysis for more general MASs assisted with conservation laws. The key technique is to transform the original network into two different MASs, defined by reconstruction and reverse reconstruction, with an important aspect that the dynamics of the original network for free species is equivalent to that of the reverse reconstruction. Stability analysis of the original network is then conducted based on an analysis of how stability properties are retained from the original network to the reverse reconstruction. We prove that the reverse reconstruction possesses only an equilibrium in each positive stoichiometric compatibility class if the corresponding reconstruction is complex balanced. Under this complex balanced reconstruction strategy, the asymptotic stability of the reverse reconstruction, which also applies to the original network, is thus reached by taking the generalized pseudo-Helmholtz function as the Lyapunov function. To facilitate applications, we further provide a systematic method for computing complex balanced reconstructions assisted with conservation laws. Some representative examples are presented to exhibit the validity of the complex balanced reconstruction strategy.


title: Role of zero synapses in unsupervised feature learning

abstract: Synapses in real neural circuits can take discrete values, including zero (silent or potential) synapses. The computational role of zero synapses in unsupervised feature learning of unlabeled noisy data is still unclear, thus it is important to understand how the sparseness of synaptic activity is shaped during learning and its relationship with receptive field formation. Here, we formulate this kind of sparse feature learning by a statistical mechanics approach. We find that learning decreases the fraction of zero synapses, and when the fraction decreases rapidly around a critical data size, an intrinsically structured receptive field starts to develop. Further increasing the data size refines the receptive field, while a very small fraction of zero synapses remain to act as contour detectors. This phenomenon is discovered not only in learning a handwritten digits dataset, but also in learning retinal neural activity measured in a natural-movie-stimuli experiment.

categories: q-bio.NC cond-mat.dis-nn cond-mat.stat-mech cs.LG cs.NE


title: Effective joint distribution of eigenvalues of Hecke operators

abstract: In 1997, Serre proved that the eigenvalues of normalised $p$-th Hecke operator $T^{'}_p$ acting on the space of cusp forms of weight $k$ and level $N$ are equidistributed in $[-2,2]$ with respect to a measure that converge to the Sato-Tate measure, whenever $N+k \to \infty$. In 2009, Murty and Sinha proved the effective version of Serre's theorem. In 2011, using Kuznetsov trace formula, Lau and Wang derived the effective joint distribution of eigenvalues of normalized Hecke operators acting on the space of primitive cusp forms of weight $k$ and level $1$. In this paper, we extend the result of Lau and Wang to space of cusp forms of higher level. Here we use Eichler-Selberg trace formula instead of Kuznetsov trace formula to deduce our result.


title: Unconventional Large Linear Magnetoresistance in Cu$_{2-x}$Te

abstract: We report a large linear magnetoresistance in Cu$_{2-x}$Te, reaching $\Delta\rho/\rho(0)$ = 250\% at 2 K in a 9 T field. This is observed for samples with $x$ in the range 0.13 to 0.22, and the results are comparable to the effects observed in Ag$_2 X$ materials, although in this case the results appear for a much wider range of bulk carrier density. Examining the magnitude vs. crossover field from low-field quadratic to high-field linear behavior, we show that models based on classical transport behavior best explain the observed results. The effects are traced to misdirected currents due to topologically inverted behavior in this system, such that stable surface states provide the high mobility transport channels. The resistivity also crosses over to a $T^2$ dependence in the temperature range where the large linear MR appears, an indicator of electron-electron interaction effects within the surface states. Thus this is an example of a system in which these interactions dominate the low-temperature behavior of the surface states.


title: Exact recursive updating of uncertainty sets for discrete-time plants with a lag

abstract: In a recent paper [arXiv:1612.04918] there are new results concerning the polytopic set of possible states of a linear discrete-time SISO system subject to bounded disturbances from measurements corrupted by bounded noise. Using these results we construct an algorithm which, for the special case of a plant with a lag, recursively updates these polytopic sets when new measurements arrive.


title: Quantitative homogenization in nonlinear elasticity for small loads

abstract: We study quantitative periodic homogenization of integral functionals in the context of non-linear elasticity. Under suitable assumptions on the energy densities (in particular frame indifference; minimality, non-degeneracy and smoothness at the identity; $p\geq d$-growth from below; and regularity of the microstructure), we show that in a neighborhood of the set of rotations, the multi-cell homogenization formula of non-convex homogenization reduces to a single-cell formula. The latter can be expressed with help of correctors. We prove that the homogenized integrand admits a quadratic Taylor expansion in an open neighborhood of the rotations -- a result that can be interpreted as the fact that homogenization and linearization commute close to the rotations. Moreover, for small applied loads, we provide an estimate on the homogenization error in terms of a quantitative two-scale expansion.


title: Fast Stochastic Variance Reduced Gradient Method with Momentum Acceleration for Machine Learning

abstract: Recently, research on accelerated stochastic gradient descent methods (e.g., SVRG) has made exciting progress (e.g., linear convergence for strongly convex problems). However, the best-known methods (e.g., Katyusha) requires at least two auxiliary variables and two momentum parameters. In this paper, we propose a fast stochastic variance reduction gradient (FSVRG) method, in which we design a novel update rule with the Nesterov's momentum and incorporate the technique of growing epoch size. FSVRG has only one auxiliary variable and one momentum weight, and thus it is much simpler and has much lower per-iteration complexity. We prove that FSVRG achieves linear convergence for strongly convex problems and the optimal $\mathcal{O}(1/T^2)$ convergence rate for non-strongly convex problems, where $T$ is the number of outer-iterations. We also extend FSVRG to directly solve the problems with non-smooth component functions, such as SVM. Finally, we empirically study the performance of FSVRG for solving various machine learning problems such as logistic regression, ridge regression, Lasso and SVM. Our results show that FSVRG outperforms the state-of-the-art stochastic methods, including Katyusha.


title: Anonymized Local Privacy

abstract: In this paper, we introduce the family of Anonymized Local Privacy mechanisms. These mechanisms have an output space of three values "Yes", "No", or "$\perp$" (not participating) and leverage the law of large numbers to generate linear noise in the number of data owners to protect privacy both before and after aggregation yet preserve accuracy. We describe the suitability in a distributed on-demand network and evaluate over a real dataset as we scale the population.


title: Failures of Gradient-Based Deep Learning

abstract: In recent years, Deep Learning has become the go-to solution for a broad range of applications, often outperforming state-of-the-art. However, it is important, for both theoreticians and practitioners, to gain a deeper understanding of the difficulties and limitations associated with common approaches and algorithms. We describe four types of simple problems, for which the gradient-based algorithms commonly used in deep learning either fail or suffer from significant difficulties. We illustrate the failures through practical experiments, and provide theoretical insights explaining their source, and how they might be remedied.


title: Sums of quadratic functions with two discriminants

abstract: Zagier in [4] discusses a construction of a function $F_{k,D}(x)$ defined for an even integer $k \geq 2$, and a positive discriminant $D$. This construction is intimately related to half-integral weight modular forms. In particular, the average value of this function is a constant multiple of the $D$-th Fourier coefficient of weight $k+1/2$ Eisenstein series constructed by H. Cohen in \cite{Cohen}. In this note we consider a construction which works both for even and odd positive integers $k$. Our function $F_{k,D,d}(x)$ depends on two discriminants $d$ and $D$ with signs sign$(d)=$ sign$(D)=(-1)^k$, degenerates to Zagier's function when $d=1$, namely, \[ F_{k,D,1}(x)=F_{k,D}(x), \] and has very similar properties. In particular, we prove that the average value of $F_{k,D,d}(x)$ is again a Fourier coefficient of H. Cohen's Eisenstein series of weight $k+1/2$, while now the integer $k \geq 2$ is allowed to be both even and odd.


title: Efficient and Robust Recovery of Sparse Signal and Image Using Generalized Nonconvex Regularization

abstract: This work addresses the robust reconstruction problem of a sparse signal from compressed measurements. We propose a robust formulation for sparse reconstruction which employs the $\ell_1$-norm as the loss function for the residual error and utilizes a generalized nonconvex penalty for sparsity inducing. The $\ell_1$-loss is less sensitive to outliers in the measurements than the popular $\ell_2$-loss, while the nonconvex penalty has the capability of ameliorating the bias problem of the popular convex LASSO penalty and thus can yield more accurate recovery. To solve this nonconvex and nonsmooth minimization formulation efficiently, we propose a first-order algorithm based on alternating direction method of multipliers (ADMM). A smoothing strategy on the $\ell_1$-loss function has been used in deriving the new algorithm to make it convergent. Further, a sufficient condition for the convergence of the new algorithm has been provided for generalized nonconvex regularization. In comparison with several state-of-the-art algorithms, the new algorithm showed better performance in numerical experiments in recovering sparse signals and compressible images. The new algorithm scales well for large-scale problems, as often encountered in image processing.


title: Fredholm conditions on non-compact manifolds: theory and examples

abstract: We give explicit Fredholm conditions for classes of pseudodifferential operators on suitable singular and non-compact spaces. In particular, we include a "user's guide" to Fredholm conditions on particular classes of manifolds including asymptotically hyperbolic manifolds, asymptotically Euclidean (or conic) manifolds, and manifolds with poly-cylindrical ends. The reader interested in applications should be able read right away the results related to those examples, beginning with Section 5. Our general, theoretical results are that an operator adapted to the geometry is Fredholm if, and only if, it is elliptic and all its limit operators, in a sense to be made precise, are invertible. Central to our theoretical results is the concept of a Fredholm groupoid, which is the class of groupoids for which this characterization of the Fredholm condition is valid. We use the notions of exhaustive and strictly spectral families of representations to obtain a general characterization of Fredholm groupoids. In particular, we introduce the class of the so-called groupoids with Exel's property as the groupoids for which the regular representations are exhaustive. We show that the class of "stratified submersion groupoids" has Exel's property, where stratified submersion groupoids are defined by glueing fibered pull-backs of bundles of Lie groups. We prove that a stratified submersion groupoid is Fredholm whenever its isotropy groups are amenable. Many groupoids, and hence many pseudodifferential operators appearing in practice, fit into this framework. This fact is explored to yield Fredholm conditions not only in the above mentioned classes, but also on manifolds that are obtained by desingularization or by blow-up of singular sets.

categories: math.OA math-ph math.AP math.FA math.MP math.SP


title: Quantum criticality at the superconductor to insulator transition revealed by specific heat measurements

abstract: The superconductor-insulator transition (SIT) is considered an excellent example of a quantum phase transition which is driven by quantum fluctuations at zero temperature. The quantum critical point is characterized by a diverging correlation length and a vanishing energy scale. Low energy fluctuations near quantum criticality may be experimentally detected by specific heat, $c_{\rm p}$, measurements. Here, we use a unique highly sensitive experiment to measure $c_{\rm p}$ of two-dimensional granular Pb films through the SIT. The specific heat shows the usual jump at the mean field superconducting transition temperature $T_{\rm c}^{\rm {mf}}$ marking the onset of Cooper pairs formation. As the film thickness is tuned toward the SIT, $T_{\rm c}^{\rm {mf}}$ is relatively unchanged, while the magnitude of the jump and low temperature specific heat increase significantly. This behaviour is taken as the thermodynamic fingerprint of quantum criticality in the vicinity of a quantum phase transition.


title: Dimensional-invariance principles in coupled dynamical systems-- A unified analysis and applications

abstract: In this paper we study coupled dynamical systems and investigate dimension properties of the subspace spanned by solutions of each individual system. Relevant problems on \textit{collinear dynamical systems} and their variations are discussed recently by Montenbruck et. al. in \cite{collinear2017SCL}, while in this paper we aim to provide a unified analysis to derive the dimensional-invariance principles for networked coupled systems, and to generalize the invariance principles for networked systems with more general forms of coupling terms. To be specific, we consider two types of coupled systems, one with scalar couplings and the other with matrix couplings. Via the \textit{rank-preserving flow theory}, we show that any scalar-coupled dynamical system (with constant, time-varying or state-dependent couplings) possesses the dimensional-invariance principles, in that the dimension of the subspace spanned by the individual systems' solutions remains invariant. For coupled dynamical systems with matrix coefficients/couplings, necessary and sufficient conditions (for constant, time-varying and state-dependent couplings) are given to characterize dimensional-invariance principles. The proofs via a rank-preserving matrix flow theory in this paper simplify the analysis in \cite{collinear2017SCL}, and we also extend the invariance principles to the cases of time-varying couplings and state-dependent couplings. Furthermore, subspace-preserving property and signature-preserving flows are also developed for coupled networked systems with particular coupling terms. These invariance principles provide insightful characterizations to analyze transient behaviors and solution evolutions for a large family of coupled systems, such as multi-agent consensus dynamics, distributed coordination systems, formation control systems, among others.


title: Generalized fermion symmetry, its currents algebra and Ward-Takahashi identities

abstract: We introduce a new local symmetry into the fermion sector of a gauge invariant Lagrangian which may or may not contain a scalar or spontaneous symmetry breaking. The Standard Model in the unitary gauge and QCD are particular cases where this symmetry may apply. We determine the associated vector and axial vector currents and their conservation laws. We show that a single current conservation law may lead to multiple Ward-Takahashi identities. Our results can potentially have important consequences for effective models of low-energy QCD and hadron structure. As an specific example, we discuss the construction of tetraquark states within a generalized linear sigma model and show that this new symmetry probes the tetraquarks in a manner that is consistent with the large $N_c$ limit of QCD.


title: Robust SfM with Little Image Overlap

abstract: Usual Structure-from-Motion (SfM) techniques require at least trifocal overlaps to calibrate cameras and reconstruct a scene. We consider here scenarios of reduced image sets with little overlap, possibly as low as two images at most seeing the same part of the scene. We propose a new method, based on line coplanarity hypotheses, for estimating the relative scale of two independent bifocal calibrations sharing a camera, without the need of any trifocal information or Manhattan-world assumption. We use it to compute SfM in a chain of up-to-scale relative motions. For accuracy, we however also make use of trifocal information for line and/or point features, when present, relaxing usual trifocal constraints. For robustness to wrong assumptions and mismatches, we embed all constraints in a parameterless RANSAC-like approach. Experiments show that we can calibrate datasets that previously could not, and that this wider applicability does not come at the cost of inaccuracy.


title: All-Optical Nanopositioning of High-Q Silica Microspheres

abstract: A tunable, all-optical, coupling method has been realized for a high-\textit{Q} silica microsphere and an optical waveguide. By means of a novel optical nanopositioning method, induced thermal expansion of an asymmetric microsphere stem for laser powers up to 171~mW has been observed and used to fine tune the microsphere-waveguide coupling. Microcavity displacements ranging from (0.612~$\pm$~0.13) -- (1.5 $\pm$ 0.13) $\mu$m and nanometer scale sensitivities varying from (2.81 $\pm$ 0.08) -- (7.39 $\pm$ 0.17) nm/mW, with an apparent linear dependency of coupling distance on stem laser heating, were obtained. Using this method, the coupling was altered such that different coupling regimes could be explored for particular samples. This tunable coupling method, in principle, could be incorporated into lab-on-a-chip microresonator systems, photonic molecule systems, and other nanopositioning frameworks.


title: The rapidity dependence of the average transverse momentum in p+p and p+Pb collisions - revisited

abstract: We revisit the rapidity dependence of the average transverse momentum $\left\langle p_{T}\right\rangle$ in a $pp$ and $pA$ collision, using the high energy factorization (HEF). We update previous predictions for the $\left\langle p_{T}\right\rangle$ both in central and forward rapidity region using parton densities following from extended BK and BFKL evolution equations to account for corrections of higher orders as well as one obtained by Sudakov resummation. Furthermore, we demonstrate that in the mid rapidity region the saturation based formalisms predict increase of the transversal momentum while in the forward rapidity region the $\left\langle p_{T}(y)\right\rangle$ decreases.


title: Distributed algorithm for controlling scale-free polygonal formations

abstract: This paper presents a distributed algorithm for controlling the deployment of a team of mobile agents in formations whose shapes can be characterized by a broad class of polygons, including regular ones, where each agent occupies a corner of the polygon. The algorithm shares the appealing properties of the popular distance-based rigid formation control, but with the additional advantage of requiring the control of fewer pairs of neighboring agents. Furthermore, the scale of the polygon can be controlled by only one pair of neighboring agents. We also exploit the exponential stability of the controlled formation in order to steer the formation as a whole with translations and rotations in a prescribed way. We provide both theoretical analysis and illustrative simulations.


title: On the Hilbert coefficients, depth of associated graded rings and reduction numbers

abstract: Let $(R,\mathfrak{m})$ be a $d$-dimensional Cohen-Macaulay local ring, $I$ an $\mathfrak{m}$-primary ideal of $R$ and $J=(x_1,...,x_d)$ a minimal reduction of $I$. We show that if $J_{d-1}=(x_1,...,x_{d-1})$ and $\sum\limits_{n=1}^\infty\lambda{({I^{n+1}\cap J_{d-1}})/({J{I^n} \cap J_{d-1}})=i}$ where i=0,1, then depth $G(I)\geq{d-i-1}$. Moreover, we prove that if $e_2(I) = \sum_{n=2}^\infty (n-1) \lambda (I^n/JI^{n-1})-2;$ or if $I$ is integrally closed and $e_2(I) = \sum_{n=2}^\infty (n-1)\lambda({{I^{n}}}/JI^{n-1})-i$ where $i=3,4$, then $e_1(I) = \sum_{n=1}^\infty \lambda(I^n / JI^{n-1})-1.$ In addition, we show that $r(I)$ is independent. Furthermore, we study the independence of $r(I)$ with some other conditions.


title: A decomposition result for Kirchhoff plate bending problems and a new discretization approach

abstract: A new approach is introduced for deriving a mixed variational formulation for Kirchhoff plate bending problems with mixed boundary conditions involving clamped, simply supported, and free boundary parts. Based on a regular decomposition of an appropriate nonstandard Sobolev space for the bending moments, the fourth-order problem can be equivalently written as a system of three (consecutively to solve) second-order problems in standard Sobolev spaces. This leads to new discretization methods, which are flexible in the sense, that any existing and well-working discretization method and solution strategy for standard second-order problems can be used as a modular building block of the new method. Similar results for the first biharmonic problem have been obtained in our previous work [W. Krendl, K. Rafetseder and W. Zulehner, A decomposition result for biharmonic problems and the Hellan-Herrmann-Johnson method, ETNA, 2016]. The extension to more general boundary conditions encounters several difficulties including the construction of an appropriate nonstandard Sobolev space, the verification of Brezzi's conditions, and the adaptation of the regular decomposition.


title: A Donsker-type Theorem for Log-likelihood Processes

abstract: Let $(\Omega, \mathcal{F}, (\mathcal{F})_{t\ge 0}, P)$ be a complete stochastic basis, $X$ a semimartingale with predictable compensator $(B, C, \nu)$. Consider a family of probability measures $\mathbf{P}=( {P}^{n, \psi}, \psi\in \Psi, n\ge 1)$, where $\Psi$ is an index set, $ {P}^{n, \psi}\stackrel {loc} \ll{P}$, and denote the likelihood ratio process by $Z_t^{n, \psi} =\frac{dP^{n, \psi}|_{\mathcal{F}_t}}{d P|_{\mathcal{F}_t}}$. Under some regularity conditions in terms of logarithm entropy and Hellinger processes, we prove that $\log Z_t^{n}$ converges weakly to a Gaussian process in $\ell^\infty(\Psi)$ as $n\rightarrow\infty$ for each fixed $t>0$.


title: Minimum Cuts and Shortest Cycles in Directed Planar Graphs via Noncrossing Shortest Paths

abstract: Let $G$ be an $n$-node simple directed planar graph with nonnegative edge weights. We study the fundamental problems of computing (1) a global cut of $G$ with minimum weight and (2) a~cycle of $G$ with minimum weight. The best previously known algorithm for the former problem, running in $O(n\log^3 n)$ time, can be obtained from the algorithm of \Lacki, Nussbaum, Sankowski, and Wulff-Nilsen for single-source all-sinks maximum flows. The best previously known result for the latter problem is the $O(n\log^3 n)$-time algorithm of Wulff-Nilsen. By exploiting duality between the two problems in planar graphs, we solve both problems in $O(n\log n\log\log n)$ time via a divide-and-conquer algorithm that finds a shortest non-degenerate cycle. The kernel of our result is an $O(n\log\log n)$-time algorithm for computing noncrossing shortest paths among nodes well ordered on a common face of a directed plane graph, which is extended from the algorithm of Italiano, Nussbaum, Sankowski, and Wulff-Nilsen for an undirected plane graph.


title: Convergence analysis of energy conserving explicit local time-stepping methods for the wave equation

abstract: Local adaptivity and mesh refinement are key to the efficient simulation of wave phenomena in heterogeneous media or complex geometry. Locally refined meshes, however, dictate a small time-step everywhere with a crippling effect on any explicit time-marching method. In [18] a leap-frog (LF) based explicit local time-stepping (LTS) method was proposed, which overcomes the severe bottleneck due to a few small elements by taking small time-steps in the locally refined region and larger steps elsewhere. Here a rigorous convergence proof is presented for the fully-discrete LTS-LF method when combined with a standard conforming finite element method (FEM) in space. Numerical results further illustrate the usefulness of the LTS-LF Galerkin FEM in the presence of corner singularities.


title: On Bernstein Type Inequalities for Stochastic Integrals of Multivariate Point Processes

abstract: We consider the stochastic integrals of multivariate point processes and study their concentration phenomena. In particular, we obtain a Bernstein type of concentration inequality through Dol\'eans-Dade exponential formula and a uniform exponential inequality using a generic chaining argument. As applications, we obtain a upper bound for a sequence of discrete time martingales indexed by a class of functionals, and so derive the rate of convergence for nonparametric maximum likelihood estimators, which is an improvement of earlier work of van de Geer.


title: Nonconvex Regularization Based Sparse Recovery and Demixing with Application to Color Image Inpainting

abstract: This work addresses the recovery and demixing problem of signals that are sparse in some general dictionary. Involved applications include source separation, image inpainting, super-resolution, and restoration of signals corrupted by clipping, saturation, impulsive noise, or narrowband interference. We employ the $\ell_q$-norm ($0 \le q < 1$) for sparsity inducing and propose a constrained $\ell_q$-minimization formulation for the recovery and demixing problem. This nonconvex formulation is approximately solved by two efficient first-order algorithms based on proximal coordinate descent and alternative direction method of multipliers (ADMM), respectively. The new algorithms are convergent in the nonconvex case under some mild conditions and scale well for high-dimensional problems. A convergence condition of the new ADMM algorithm has been derived. Furthermore, extension of the two algorithms for multi-channels joint recovery has been presented, which can further exploit the joint sparsity pattern among multi-channel signals. Various numerical experiments showed that the new algorithms can achieve considerable performance gain over the $\ell_1$-regularized algorithms.