title: Simplified broken Lefschetz fibrations and trisections of 4-manifolds

abstract: Shapes of four dimensional spaces can be studied effectively via maps to standard surfaces. We explain, and illustrate by quintessential examples, how to simplify such generic maps on 4-manifolds topologically, in order to derive simple decompositions into much better understood manifold pieces. Our methods not only allow us to produce various interesting families of examples, but also to establish a correspondence between simplified broken Lefschetz fibrations and simplified trisections of closed, oriented 4-manifolds.


title: An exciton-coupled electron transfer process controlled by non-Markovian environments

abstract: We theoretically investigate an exciton-coupled electron transfer (XCET) process that is conversion of an exciton into a charge transfer state. This conversion happens in an exciton transfer (XT) process, and the electron moves away in an electron transfer(ET) process in multiple environments (baths). This XCET process plays an essential role in the harvesting of solar energy in biological and photovoltaic materials. We develop a practical theoretical model to study the efficiency of XCET process that occurs either in consecutive or concerted processes under the influence of non-Markovian baths. The role of quantum coherence in the XT-ET system and the baths is investigated using reduced hierarchal equations of motion (HEOM). This model includes independent baths for each XT and ET state, in addition to a XCET bath for the conversion process. We found that, while quantum system-bath coherence is important in the XT and ET processes, coherence between the XT and ET processes must be suppressed in order to realize efficient irreversible XCET process through the weak off-diagonal interaction between the XT and ET bridge sites arises from a XCET bath.


title: Social Learning Against Data Falsification in Sensor Networks

abstract: Although surveillance and sensor networks play a key role in Internet of Things, sensor nodes are usually vulnerable to tampering due to their widespread locations. In this letter we consider data falsification attacks where an smart attacker takes control of critical nodes within the network, including nodes serving as fusion centers. In order to face this critical security thread, we propose a data aggregation scheme based on social learning, resembling the way in which agents make decisions in social networks. Our results suggest that social learning enables network resilience, even when a significant portion of the nodes have been compromised by the attacker. Finally, we show the suitability of our scheme to sensor networks by developing a low-complexity algorithm to facilitate the social learning data fusion rule in devices with restricted computational power.


title: Data-driven spectral analysis of the Koopman operator

abstract: Starting from measured data, we develop a method to compute the fine structure of the spectrum of the Koopman operator with rigorous convergence guarantees. The method is based on the observation that, in the measure-preserving ergodic setting, the moments of the spectral measure associated to a given observable are computable from a single trajectory of this observable. Having finitely many moments available, we use the classical Christoffel-Darboux kernel to separate the atomic and absolutely continuous parts of the spectrum, supported by convergence guarantees as the number of moments tends to infinity. In addition, we propose a technique to detect the singular continuous part of the spectrum as well as two methods to approximate the spectral measure with guaranteed convergence in the weak topology, irrespective of whether the singular continuous part is present or not. The proposed method is simple to implement and readily applicable to large-scale systems since the computational complexity is dominated by inverting an $N\times N$ Hermitian positive-definite Toeplitz matrix, where $N$ is the number of moments, for which efficient and numerically stable algorithms exist; in particular, the complexity of the approach is independent of the dimension of the underlying state-space. We also show how to compute, from measured data, the spectral projection on a given segment of the unit circle, allowing us to obtain a finite-dimensional approximation of the operator that explicitly takes into account the point and continuous parts of the spectrum. Finally, we describe a relationship between the proposed method and the so-called Hankel Dynamic Mode Decomposition, providing new insights into the behavior of the eigenvalues of the Hankel DMD operator. A number of numerical examples illustrate the approach, including a study of the spectrum of the lid-driven two-dimensional cavity flow.


title: Modelling chemical abundance distributions for dwarf galaxies in the Local Group: the impact of turbulent metal diffusion

abstract: We investigate stellar metallicity distribution functions (MDFs), including Fe and ${\alpha}$-element abundances, in dwarf galaxies from the Feedback in Realistic Environments (FIRE) project. We examine both isolated dwarf galaxies and those that are satellites of a Milky Way-mass galaxy. In particular, we study the effects of including a sub-grid turbulent model for the diffusion of metals in gas. Simulations that include diffusion have narrower MDFs and abundance ratio distributions, because diffusion drives individual gas and star particles toward the average metallicity. This effect provides significantly better agreement with observed abundance distributions of dwarf galaxies in the Local Group, including the small intrinsic scatter in [${\alpha}$/Fe] vs. [Fe/H] (less than 0.1 dex). This small intrinsic scatter arises in our simulations because the interstellar medium (ISM) in dwarf galaxies is well-mixed at nearly all cosmic times, such that stars that form at a given time have similar abundances to within 0.1 dex. Thus, most of the scatter in abundances at z = 0 arises from redshift evolution and not from instantaneous scatter in the ISM. We find similar MDF widths and intrinsic scatter for satellite and isolated dwarf galaxies, which suggests that environmental effects play a minor role compared with internal chemical evolution in our simulations. Overall, with the inclusion of metal diffusion, our simulations reproduce abundance distribution widths of observed low-mass galaxies, enabling detailed studies of chemical evolution in galaxy formation.


title: Lower Bounds for Numbers of Real Self-Dual Spaces in Problems of Schubert Calculus

abstract: The self-dual spaces of polynomials are related to Bethe vectors in the Gaudin model associated to the Lie algebras of types B and C. In this paper, we give lower bounds for the numbers of real self-dual spaces in intersections of Schubert varieties related to osculating flags in the Grassmannian. The higher Gaudin Hamiltonians are self-adjoint with respect to a nondegenerate indefinite Hermitian form. Our bound comes from the computation of the signature of this form.


title: Maximum number of colourings. II. 5-chromatic graphs

abstract: In 1971, Tomescu conjectured [Le nombre des graphes connexes $k$-chromatiques minimaux aux sommets \'etiquet\'es, C. R. Acad. Sci. Paris 273 (1971), 1124--1126] that every connected graph $G$ on $n$ vertices with $\chi(G) = k \geq 4$ has at most $k!(k-1)^{n-k}$ $k$-colourings, where equality holds if and only if the graph is formed from $K_k$ by repeatedly adding leaves. In this note we prove (a strengthening of) the conjecture of Tomescu when $k=5$.


title: Basic tasks of sentiment analysis

abstract: Subjectivity detection is the task of identifying objective and subjective sentences. Objective sentences are those which do not exhibit any sentiment. So, it is desired for a sentiment analysis engine to find and separate the objective sentences for further analysis, e.g., polarity detection. In subjective sentences, opinions can often be expressed on one or multiple topics. Aspect extraction is a subtask of sentiment analysis that consists in identifying opinion targets in opinionated text, i.e., in detecting the specific aspects of a product or service the opinion holder is either praising or complaining about.


title: Sim-to-Real Transfer of Robotic Control with Dynamics Randomization

abstract: Simulations are attractive environments for training agents as they provide an abundant source of data and alleviate certain safety concerns during the training process. But the behaviours developed by agents in simulation are often specific to the characteristics of the simulator. Due to modeling error, strategies that are successful in simulation may not transfer to their real world counterparts. In this paper, we demonstrate a simple method to bridge this "reality gap". By randomizing the dynamics of the simulator during training, we are able to develop policies that are capable of adapting to very different dynamics, including ones that differ significantly from the dynamics on which the policies were trained. This adaptivity enables the policies to generalize to the dynamics of the real world without any training on the physical system. Our approach is demonstrated on an object pushing task using a robotic arm. Despite being trained exclusively in simulation, our policies are able to maintain a similar level of performance when deployed on a real robot, reliably moving an object to a desired location from random initial configurations. We explore the impact of various design decisions and show that the resulting policies are robust to significant calibration error.


title: $L^p$-$L^q$ estimates for the damped wave equation and the critical exponent for the nonlinear problem with slowly decaying data

abstract: We study the Cauchy problem of the damped wave equation \begin{align*} \partial_{t}^2 u - \Delta u + \partial_t u = 0 \end{align*} and give sharp $L^p$-$L^q$ estimates of the solution for $1\le q \le p < \infty\ (p\neq 1)$ with derivative loss. This is an improvement of the so-called Matsumura estimates. Moreover, as its application, we consider the nonlinear problem with initial data in $(H^s\cap H_r^{\beta}) \times (H^{s-1} \cap L^r)$ with $r \in (1,2]$, $s\ge 0$, and $\beta = (n-1)|\frac{1}{2}-\frac{1}{r}|$, and prove the local and global existence of solutions. In particular, we prove the existence of the global solution with small initial data for the critical nonlinearity with the power $1+\frac{2r}{n}$, while it is known that the critical power $1+\frac{2}{n}$ belongs to the blow-up region when $r=1$. We also discuss the asymptotic behavior of the global solution in supercritical cases. Moreover, we present blow-up results in subcritical cases. We give estimates of lifespan and blow-up rates by an ODE argument.


title: Computing parametrized solutions for plasmonic nanogap structures

abstract: The interaction of electromagnetic waves with metallic nanostructures generates resonant oscillations of the conduction-band electrons at the metal surface. These resonances can lead to large enhancements of the incident field and to the confinement of light to small regions, typically several orders of magnitude smaller than the incident wavelength. The accurate prediction of these resonances entails several challenges. Small geometric variations in the plasmonic structure may lead to large variations in the electromagnetic field responses. Furthermore, the material parameters that characterize the optical behavior of metals at the nanoscale need to be determined experimentally and are consequently subject to measurement errors. It then becomes essential that any predictive tool for the simulation and design of plasmonic structures accounts for fabrication tolerances and measurement uncertainties. In this paper, we develop a reduced order modeling framework that is capable of real-time accurate electromagnetic responses of plasmonic nanogap structures for a wide range of geometry and material parameters. The main ingredients of the proposed method are: (i) the hybridizable discontinuous Galerkin method to numerically solve the equations governing electromagnetic wave propagation in dielectric and metallic media, (ii) a reference domain formulation of the time-harmonic Maxwell's equations to account for geometry variations; and (iii) proper orthogonal decomposition and empirical interpolation techniques to construct an efficient reduced model. To demonstrate effectiveness of the models developed, we analyze geometry sensitivities and explore optimal designs of a 3D periodic annular nanogap structure.


title: A Particle Filtering Approach for Enabling Distributed and Scalable Sharing of DSA Network Resources

abstract: Handling the massive number of devices needed in numerous applications such as smart cities is a major challenge given the scarcity of spectrum resources. Dynamic spectrum access (DSA) is seen as a potential candidate to support the connectivity and spectrum access of these devices. We propose an efficient technique that relies on particle filtering to enable distributed resource allocation and sharing for large-scale dynamic spectrum access networks. More specifically, we take advantage of the high tracking capability of particle filtering to efficiently assign the available spectrum and power resources among cognitive users. Our proposed technique maximizes the per-user throughput while ensuring fairness among users, and it does so while accounting for the different users' quality of service requirements and the channel gains' variability. Through intensive simulations, we show that our proposed approach performs well by achieving high overall throughput while improving user's fairness under different objective functions. Furthermore, it achieves higher performance when compared to state-of-the-art techniques.


title: Design Considerations of a Sub-50 {\mu}W Receiver Front-end for Implantable Devices in MedRadio Band

abstract: Emerging health-monitor applications, such as information transmission through multi-channel neural implants, image and video communication from inside the body etc., calls for ultra-low active power (<50${\mu}$W) high data-rate, energy-scalable, highly energy-efficient (pJ/bit) radios. Previous literature has strongly focused on low average power duty-cycled radios or low power but low-date radios. In this paper, we investigate power performance trade-off of each front-end component in a conventional radio including active matching, down-conversion and RF/IF amplification and prioritize them based on highest performance/energy metric. The analysis reveals 50${\Omega}$ active matching and RF gain is prohibitive for 50${\mu}$W power-budget. A mixer-first architecture with an N-path mixer and a self-biased inverter based baseband LNA, designed in TSMC 65nm technology show that sub 50${\mu}$W performance can be achieved up to 10Mbps (< 5pJ/b) with OOK modulation.


title: Asymmetric Actor Critic for Image-Based Robot Learning

abstract: Deep reinforcement learning (RL) has proven a powerful technique in many sequential decision making domains. However, Robotics poses many challenges for RL, most notably training on a physical system can be expensive and dangerous, which has sparked significant interest in learning control policies using a physics simulator. While several recent works have shown promising results in transferring policies trained in simulation to the real world, they often do not fully utilize the advantage of working with a simulator. In this work, we exploit the full state observability in the simulator to train better policies which take as input only partial observations (RGBD images). We do this by employing an actor-critic training algorithm in which the critic is trained on full states while the actor (or policy) gets rendered images as input. We show experimentally on a range of simulated tasks that using these asymmetric inputs significantly improves performance. Finally, we combine this method with domain randomization and show real robot experiments for several tasks like picking, pushing, and moving a block. We achieve this simulation to real world transfer without training on any real world data.


title: Photometric long-term variations and superhump occurrence in the Classical Nova RR Pictoris

abstract: We present an analysis of all available time-resolved photometry from the literature and new light curves obtained in 2013-2014 for the old nova RR Pictoris. The well-known hump light curve phased with the orbital period reveals significant variations over the last 42 years in shape, amplitude and other details which apparently are caused by long-term variations in the disc structure. In addition we found evidence for the presence of superhumps in 2007, with the same period (~9% longer than the orbital period), as reported earlier by other authors from observations in 2005. Possibly, superhumps arise quickly in RR Pic, but are sporadic events, because in all the other observing runs analysed no significant superhump signal was detected. We also determined an actual version of the Stolz--Schoembs relation between superhump period and orbital period, analysing separately dwarf novae, classical novae and nova-like stars, and conclude that this relation is of general validity for all superhumpers among the cataclysmic variables (CVs), in spite of small but significant differences among the sub-types mentioned above. We emphasize the importance of such a study in context with the still open question of the interrelation between the different sub-classes of CVs, crucial for our understanding of the long-term CV evolution.


title: New inertial regularized algorithm for solving strongly pseudomonotone equilibrium problems

abstract: The article introduces a new algorithm for solving a class ofequilibrium problems involving strongly pseudomonotone bifunctions with Lipschitz-type condition. We describe how to incorporate the proximal-like regularized technique with inertial effects. The main novelty of the algorithm is that it can be done without previously knowing the information on the strongly pseudomonotone and Lipschitz-type constants of cost bifunction. A reasonable explain for this is that the algorithm uses a sequence of stepsizes which is diminishing and non-summable. Theorem of strong convergence is proved. In the case, when the information on the modulus of strong pseudomonotonicity and Lispchitz-type constant is known,the rate of linear convergence of the algorithm has been established. Several of numerical experiments are performed to illustrate the convergence of the algorithm and also compare it with other algorithms.


title: Block-induced complex structures building the flare-productive solar active region 12673

abstract: Solar active region (AR) 12673 produced 4 X-class, 27 M-class, and numerous lower class flares during its passage across the visible solar disk in September 2017. Our study is to answer the questions why this AR was so flare-productive and how the X9.3 flare, the largest one of the last decade, took place. We find that there was a sunspot in the initial several days, and then two bipolar regions emerged nearby it successively. Due to the standing of the pre-existing sunspot, the movement of the bipoles was blocked, while the pre-existing sunspot maintained its quasi-circular shaped umbra only with the disappearance of a part of penumbra. Thus, the bipolar patches were significantly distorted, and the opposite polarities formed two semi-circular shaped structures. After that, two sequences of new bipolar regions emerged within the narrow semi-circular zone, and the bipolar patches separated along the curved channel. The new bipoles sheared and interacted with the previous ones, forming a complex topological system, during which numerous flares occurred. At the highly sheared region, a great deal of free energy was accumulated. On September 6, one negative patch near the polarity inversion line began to rapidly rotate and shear with the surrounding positive fields, and consequently the X9.3 flare erupted. Our results reveal that the block-induced complex structures built the flare-productive AR and the X9.3 flare was triggered by an erupting filament due to the kink instability. To better illustrate this process, a block-induced eruption model is proposed for the first time.


title: Absolute Properties of the Pulsating Post-mass Transfer Eclipsing Binary OO Draconis

abstract: OO Dra is a short-period Algol system with a $\delta$ Sct-like pulsator. We obtained time-series spectra between 2016 February and May to derive the fundamental parameters of the binary star and to study its evolutionary scenario. The radial velocity (RV) curves for both components were presented, and the effective temperature of the hotter and more massive primary was determined to be $T_{\rm eff,1}$ = 8260 $\pm$ 210 K by comparing the disentangling spectrum and the Kurucz models. Our RV measurements were solved with the $BV$ light curves of Zhang et al. (2014) using the Wilson-Devinney binary code. The absolute dimensions of each component are determined as follows: $M_1$ = 2.03 $\pm$ 0.06 M$_\odot$, $M_2$ = 0.19 $\pm$ 0.01 M$_\odot$, $R_1$ = 2.08 $\pm$ 0.03 R$_\odot$, $R_2$ = 1.20 $\pm$ 0.02 R$_\odot$, $L_1$ = 18 $\pm$ 2 L$_\odot$, and $L_2$ = 2.0 $\pm$ 0.2 L$_\odot$. Comparison with stellar evolution models indicated that the primary star resides inside the $\delta$ Sct instability strip on the main sequence, while the cool secondary component is noticeably overluminous and oversized. We demonstrated that OO Dra is an oscillating post-mass transfer R CMa-type binary; the originally more massive star became the low-mass secondary component through mass loss caused by stellar wind and mass transfer, and the gainer became the pulsating primary as the result of mass accretion. The R CMa stars, such as OO Dra, are thought to have formed by non-conservative binary evolution and ultimately to evolve into EL CVn stars.


title: The Behavior of Selected Diffuse Interstellar Bands with Molecular Fraction in Diffuse Atomic and Molecular Clouds

abstract: We study the behavior of eight diffuse interstellar bands (DIBs) in different interstellar environments, as characterized by the fraction of hydrogen in molecular form [$f$(H$_2$)], with comparisons to the corresponding behavior of various known atomic and molecular species. The equivalent widths of the five "normal" DIBs ($\lambda\lambda$5780.5, 5797.1, 6196.0, 6283.8, and 6613.6), normalized to $E(B-V)$, show a "Lambda-shaped" behavior: they increase at low $f$(H$_2$), peak at $f$(H$_2$) ~ 0.3, and then decrease. The similarly normalized column densities of Ca, Ca$^+$, Ti$^+$, and CH$^+$ also decline for $f$(H$_2$) > 0.3. In contrast, the normalized column densities of Na, K, CH, CN, and CO increase monotonically with $f$(H$_2$), and the trends exhibited by the three C$_2$ DIBs ($\lambda\lambda$4726.8, 4963.9, and 4984.8) lie between those two general behaviors. These trends with $f$(H$_2$) are accompanied by cosmic scatter, the dispersion at any given $f$(H$_2$) being significantly larger than the individual errors of measurement. The Lambda-shaped trends suggest the balance between creation and destruction of the DIB carriers differs dramatically between diffuse atomic and diffuse molecular clouds; additional processes besides ionization and shielding are needed to explain those observed trends. Except for several special cases, the highest $W$(5780)/$W$(5797) ratios, characterizing the so-called "sigma-zeta effect", occur only at $f$(H$_2$) < 0.2. We propose a sequence of DIBs based on trends in their pair-wise strength ratios with increasing $f$(H$_2$). In order of increasing environmental density, we find the $\lambda$6283.8 and $\lambda$5780.5 DIBs, the $\lambda$6196.0 DIB, the $\lambda$6613.6 DIB, the $\lambda$5797.1 DIB, and the C$_2$ DIBs.


title: Data Driven Computational Model for Bipedal Walking and Push Recovery

abstract: In this research, we have developed the data driven computational walking model to overcome the problem with traditional kinematics based model. Our model is adaptable and can adjust the parameter morphological similar to human. The human walk is a combination of different discrete sub-phases with their continuous dynamics. Any system which exhibits the discrete switching logic and continuous dynamics can be represented using a hybrid system. In this research, the bipedal locomotion is analyzed which is important for understanding the stability and to negotiate with the external perturbations. We have also studied the other important behavior push recovery. The Push recovery is also a very important behavior acquired by human with continuous interaction with environment. The researchers are trying to develop robots that must have the capability of push recovery to safely maneuver in a dynamic environment. The push is a very commonly experienced phenomenon in cluttered environment. The human beings can recover from external push up to a certain extent using different strategies of hip, knee and ankle. The different human beings have different push recovery capabilities. For example a wrestler has a better push negotiation capability compared to normal human beings. The push negotiation capability acquired by human, therefore, is based on learning but the learning mechanism is still unknown to researchers. The research community across the world is trying to develop various humanoid models to solve this mystery. Seeing all the conventional mechanics and control based models have some inherent limitations, a learning based computational model has been developed to address effectively this issue. In this research we will discuss how we have framed this problem as hybrid system.


title: Flavor-singlet spectrum in multi-flavor QCD

abstract: Studying SU(3) gauge theories with increasing number of light fermions is relevant both for understanding the strong dynamics of QCD and for constructing strongly interacting extensions of the Standard Model (e.g. UV completions of composite Higgs models). In order to contrast these many-flavors strongly interacting theories with QCD, we study the flavor-singlet spectrum as an interesting probe. In fact, some composite Higgs models require the Higgs boson to be the lightest flavor-singlet scalar in the spectrum of a strongly interacting new sector with a well defined hierarchy with the rest of the states. Moreover, introducing many light flavors at fixed number of colors can influence the dynamics of the lightest flavor-singlet pseudoscalar. We present the on-going study of these flavor-singlet channels using multiple interpolating operators on high-statistics ensembles generated by the LatKMI collaboration and we compare results with available data obtained by the Lattice Strong Dynamics collaboration. For the theory with 8 flavors, the two collaborations have generated configurations that complement each others with the aim to tackle the massless limit using the largest possible volumes.


title: Designing Planar, Ultra-Thin, Broad-Band and Material-Versatile Solar Absorbers via Bound-Electron and Exciton Absorption

abstract: Ultrathin planar absorbing layers, including semiconductor and metal films, and 2D materials, are promising building blocks for solar energy harvesting devices but poor light absorption has been a critical issue. Although interference in ultrathin absorbing layers has been studied to realize near perfect absorption at a specific wavelength, achieving high broadband absorption still remains challenging. Here, we both theoretically and experimentally demonstrated a method to tune not only reflection phase shift but also electromagnetic energy dissipation to design broadband solar absorber with simple planar structure consisting of an ultrathin absorbing layer separated from the metallic substrate by a transparent layer. We explicitly identified by deriving a new formulism that the absorbing material with refractive index proportional to the wavelength as well as extinction coefficient independent of the wavelength, is the ideal building block to create ultrathin planar broadband absorbers. To demonstrate the general strategy for naturally available absorbing materials in both high-loss (refractory metals) and weak-absorption (2D materials) regimes, we leveraged the bound-electron interband transition with a broad Lorentz oscillator peak to design a solar thermal absorber based on a ultrathin Cr layer; and leveraged the strong exciton attributed to the spin-orbit coupling for the spectrum near the band edge, and the bound-electron interband transition for shorter wavelengths, to design a solar photovoltaic absorber based on a atomically thin MoS2 layer. Furthermore, our designed ultrathin broadband solar absorbers with planar structures have comparable absorption properties compared to the absorbers with nanopatterns. Our proposed design strategies pave the way to novel nanometer thick energy harvesting and optoelectronic devices with simple planar structures.


title: Exploiting oddsmaker bias to improve the prediction of NFL outcomes

abstract: Accurately predicting the outcome of sporting events has been a goal for many groups who seek to maximize profit. What makes this challenging is that the outcome of an event can be influenced by many factors that dynamically change across time. Oddsmakers attempt to estimate these factors by using both algorithmic and subjective methods to set the spread. However, it is well-known that both human and algorithmic decision-making can be biased, so this paper explores if oddsmaker biases can be used in an exploitative manner, in order to improve the prediction of NFL game outcomes. Real-world gambling data was used to train and test different predictive models under varying assumptions. The results show that methods that leverage oddsmaker biases in an exploitative manner perform best under the conditions tested in this paper. These findings suggest that leveraging human and algorithmic decision biases in an exploitative manner may be useful for predicting the outcomes of competitive events, and could lead to increased profit for those who have financial interest in the outcomes.


title: Relaxation-Based Coarsening for Multilevel Hypergraph Partitioning

abstract: Multilevel partitioning methods that are inspired by principles of multiscaling are the most powerful practical hypergraph partitioning solvers. Hypergraph partitioning has many applications in disciplines ranging from scientific computing to data science. In this paper we introduce the concept of algebraic distance on hypergraphs and demonstrate its use as an algorithmic component in the coarsening stage of multilevel hypergraph partitioning solvers. The algebraic distance is a vertex distance measure that extends hyperedge weights for capturing the local connectivity of vertices which is critical for hypergraph coarsening schemes. The practical effectiveness of the proposed measure and corresponding coarsening scheme is demonstrated through extensive computational experiments on a diverse set of problems. Finally, we propose a benchmark of hypergraph partitioning problems to compare the quality of other solvers.


title: On a modification method of Lefschetz thimbles

abstract: The QCD at finite density is not well understood yet, where standard Monte Carlo simulation suffers from the sign problem. In order to overcome the sign problem, the method of Lefschetz thimble has been explored. Basically, the original sign problem can be less severe in a complexified theory due to the constancy of the imaginary part of an action on each thimble. However, global phase factors assigned on each thimble still remain. Their interference is not negligible in a situation where a large number of thimbles contribute to the partition function, and this could also lead to a sign problem.In this study, we propose a method to resolve this problem by modifying the structure of Lefschetz thimbles such that only a single thimble is relevant to the partition function. It can be shown that observables measured in the original and modified theories are connected by a simple identity. We exemplify that our method works well in a toy model.


title: Honk: A PyTorch Reimplementation of Convolutional Neural Networks for Keyword Spotting

abstract: We describe Honk, an open-source PyTorch reimplementation of convolutional neural networks for keyword spotting that are included as examples in TensorFlow. These models are useful for recognizing "command triggers" in speech-based interfaces (e.g., "Hey Siri"), which serve as explicit cues for audio recordings of utterances that are sent to the cloud for full speech recognition. Evaluation on Google's recently released Speech Commands Dataset shows that our reimplementation is comparable in accuracy and provides a starting point for future work on the keyword spotting task.


title: Learning Deep Context-aware Features over Body and Latent Parts for Person Re-identification

abstract: Person Re-identification (ReID) is to identify the same person across different cameras. It is a challenging task due to the large variations in person pose, occlusion, background clutter, etc How to extract powerful features is a fundamental problem in ReID and is still an open problem today. In this paper, we design a Multi-Scale Context-Aware Network (MSCAN) to learn powerful features over full body and body parts, which can well capture the local context knowledge by stacking multi-scale convolutions in each layer. Moreover, instead of using predefined rigid parts, we propose to learn and localize deformable pedestrian parts using Spatial Transformer Networks (STN) with novel spatial constraints. The learned body parts can release some difficulties, eg pose variations and background clutters, in part-based representation. Finally, we integrate the representation learning processes of full body and body parts into a unified framework for person ReID through multi-class person identification tasks. Extensive evaluations on current challenging large-scale person ReID datasets, including the image-based Market1501, CUHK03 and sequence-based MARS datasets, show that the proposed method achieves the state-of-the-art results.


title: On the action of the Steenrod algebra on the modular invariants of special linear group

abstract: We compute the action of the Steenrod algebra on generators of algebras of invariants of special linear group ${SL_n=SL(n,\mathbb Z/p)}$ in the polynomial algebra with $ p$ an odd prime number.


title: Atomic Gas in Blue Ultra Diffuse Galaxies around Hickson Compact Groups

abstract: We have found the atomic gas (HI) reservoirs of the blue ultra diffuse galaxy (UDG) candidates identified by R\'oman and Trujillo in images near Hickson Compact Groups (HCGs). We confirm that all of the objects are indeed UDGs with effective radii R_e > 1.5 kpc. Three of them are likely to be gravitationally bound to the HCG near which they project, one is plausibly gravitationally bound to the nearest HCG, and one is in the background. We measure HI masses and velocity widths for each object directly from the spectra, and use the widths together with the UDG effective radii to estimate dynamical masses and halo spin parameters. The location of the blue UDGs in the HI mass - stellar mass plane is consistent with that of the broader gas-rich galaxy population, and both their HI masses and gas richnesses are correlated with their effective radii. The blue UDGs appear to be low-mass objects with high-spin halos, although their properties are not as extreme as those of the faintest diffuse objects found in HI searches. The data presented here highlight the potential of single-dish radio observations for measuring the physical properties of blue diffuse objects detected in the optical.


title: Diverse Trends of Electron Correlation Effects for Properties with Different Radial and Angular Factors in an Atomic System: A case study in Ca$^{+}$

abstract: Atomic properties such as field shift constants, magnetic dipole and electric quadrupole hyperfine structure constants, Land\'e $g_J$ factors, and electric quadrupole moments that are described by electronic operators with different ranks and radial behaviors are studied and the role of electron correlation effects in their determination are investigated. We have adopted the Dirac-Hartree-Fock method, the second- and third-order relativistic many-body perturbation theories, and an all-order relativistic many-body method in the coupled-cluster theory framework considering only the linearized terms and also all the non-linearized terms in the singles and doubles with partial triples excitations approximation to carry out these analyses. Variations in the propagation of electron correlation effects with operators having same angular factors but different radial behaviors and with different ranks are highlighted. Corrections from the higher-order relativistic corrections due to the Breit and quantum electrodynamics interactions to all these properties are also estimated. Understanding of trends of electron correlation effects in these properties can be useful to establish accuracies in the theoretical results of different atomic properties and to substantiate validity of an approximated many-body method.

categories: physics.atom-ph physics.chem-ph physics.comp-ph quant-ph


title: A recognition algorithm for simple-triangle graphs

abstract: A simple-triangle graph is the intersection graph of triangles that are defined by a point on a horizontal line and an interval on another horizontal line. The time complexity of the recognition problem for simple-triangle graphs was a longstanding open problem, which was recently settled. This paper provides a new recognition algorithm for simple-triangle graphs to improve the time bound from $O(n^2 \overline{m})$ to $O(nm)$, where $n$, $m$, and $\overline{m}$ are the number of vertices, edges, and non-edges of the graph, respectively. The algorithm uses the vertex ordering characterization that a graph is a simple-triangle graph if and only if there is a linear ordering of the vertices containing both an alternating orientation of the graph and a transitive orientation of the complement of the graph. We also show, as a byproduct, that an alternating orientation can be obtained in $O(nm)$ time for cocomparability graphs, and it is NP-complete to decide whether a graph has an orientation that is alternating and acyclic.


title: Increasing the smoothness of vector and Hermite subdivision schemes

abstract: In this paper we suggest a method for transforming a vector subdivision scheme generating $C^{\ell}$ limits to another such scheme of the same dimension, generating $C^{\ell+1}$ limits. In scalar subdivision, it is well known that a scheme generating $C^{\ell}$ limit curves can be transformed to a new scheme producing $C^{\ell+1}$ limit curves by multiplying the scheme's symbol with the smoothing factor $\tfrac{z+1}{2}$. We extend this approach to vector and Hermite subdivision schemes, by manipulating symbols. The algorithms presented in this paper allow to construct vector (Hermite) subdivision schemes of arbitrarily high regularity from a convergent vector scheme (from a Hermite scheme whose Taylor scheme is convergent with limit functions of vanishing first component).


title: Revenue-based Attribution Modeling for Online Advertising

abstract: This paper examines and proposes several attribution modeling methods that quantify how revenue should be attributed to online advertising inputs. We adopt and further develop relative importance method, which is based on regression models that have been extensively studied and utilized to investigate the relationship between advertising efforts and market reaction (revenue). Relative importance method aims at decomposing and allocating marginal contributions to the coefficient of determination (R^2) of regression models as attribution values. In particular, we adopt two alternative submethods to perform this decomposition: dominance analysis and relative weight analysis. Moreover, we demonstrate an extension of the decomposition methods from standard linear model to additive model. We claim that our new approaches are more flexible and accurate in modeling the underlying relationship and calculating the attribution values. We use simulation examples to demonstrate the superior performance of our new approaches over traditional methods. We further illustrate the value of our proposed approaches using a real advertising campaign dataset.


title: Hydrodynamic limit and Propagation of Chaos for Brownian Particles reflecting from a Newtonian barrier

abstract: In 2001, Knight constructed a stochastic process modeling the one dimensional interaction of two particles, one being Newtonian in the sense that it obeys Newton's laws of motion, and the other particle being Brownian. We construct a multi-particle analog, using Skorohod map estimates in proving a propagation of chaos and characterizing the hydrodynamic limit as the solution to a PDE with free boundary condition. Stochastic methods are used to show existence and uniqueness for the free boundary problem, and also present an algorithm of approximating the solution.


title: Long-time asymptotics for initial-boundary value problems of integrable Fokas-Lenells equation on the half-line

abstract: We study the Schwartz class of initial-boundary value (IBV) problems for the integrable Fokas-Lenells equation on the half-line via the Deift-Zhou's nonlinear descent method analysis of the corresponding Riemann-Hilbert problem such that the asymptotics of the Schwartz class of IBV problems as t\to\infty is presented.

categories: nlin.SI math-ph math.AP math.MP quant-ph


title: Replacement AutoEncoder: A Privacy-Preserving Algorithm for Sensory Data Analysis

abstract: An increasing number of sensors on mobile, Internet of things (IoT), and wearable devices generate time-series measurements of physical activities. Though access to the sensory data is critical to the success of many beneficial applications such as health monitoring or activity recognition, a wide range of potentially sensitive information about the individuals can also be discovered through access to sensory data and this cannot easily be protected using traditional privacy approaches. In this paper, we propose a privacy-preserving sensing framework for managing access to time-series data in order to provide utility while protecting individuals' privacy. We introduce Replacement AutoEncoder, a novel algorithm which learns how to transform discriminative features of data that correspond to sensitive inferences, into some features that have been more observed in non-sensitive inferences, to protect users' privacy. This efficiency is achieved by defining a user-customized objective function for deep autoencoders. Our replacement method will not only eliminate the possibility of recognizing sensitive inferences, it also eliminates the possibility of detecting the occurrence of them. That is the main weakness of other approaches such as filtering or randomization. We evaluate the efficacy of the algorithm with an activity recognition task in a multi-sensing environment using extensive experiments on three benchmark datasets. We show that it can retain the recognition accuracy of state-of-the-art techniques while simultaneously preserving the privacy of sensitive information. Finally, we utilize the GANs for detecting the occurrence of replacement, after releasing data, and show that this can be done only if the adversarial network is trained on the users' original data.


title: Efficiency at maximum power of a quantum Carnot engine with temperature tunable baths

abstract: We investigate the efficiency at maximum power (EMP) of irreversible quantum Carnot engines that perform finite-time cycles between two temperature tunable baths. The temperature form we adopt can be experimentally realized in squeezed baths in the high temperature limit, which makes our proposal of practical relevance. Focusing on low dissipation engines, we first generalize the pervious upper as well as lower bounds for the EMP to temperature tunable cases in which they are solely determined by a generalized Carnot limit. As an illustrative example, we then consider a minimal heat engine model with a two-level spin as the working medium. It mimics a low dissipation engine as confirmed by finite time thermodynamic optimization results. The so-obtained EMP, being constrained by the generalized bounds, is well described by a generalized Curzon- Ahlborn efficiency as consequences of a left/right symmetry for a rate constant and low dissipations. Intriguing features of this minimal heat engine under optimal power output are also demonstrated.


title: On Least Squares Linear Regression Without Second Moment

abstract: If X and Y are real valued random variables such that the first moments of X, Y, and XY exist and the conditional expectation of Y given X is an affine function of X, then the intercept and slope of the conditional expectation equal the intercept and slope of the least squares linear regression function, even though Y may not have a finite second moment. As a consequence, the affine in X form of the conditional expectation and zero covariance imply mean independence.


title: Resummed tree heptagon

abstract: The form factor program for the regularized space-time S-matrix in planar maximally supersymmetric gauge theory, known as the pentagon operator product expansion, is formulated in terms of flux-tube excitations propagating on a dual two-dimensional world-sheet, whose dynamics is known exactly as a function of 't Hooft coupling. Both MHV and non-MHV amplitudes are described in a uniform, systematic fashion within this framework, with the difference between the two encoded in coupling-dependent helicity form factors expressed via Zhukowski variables. The nontrivial SU(4) tensor structure of flux-tube transitions is coupling independent and is known for any number of charged excitations from solutions of a system of Watson and Mirror equations. This description allows one to resum the infinite series of form factors and recover the space-time S-matrix exactly in kinematical variables at a given order of perturbation series. Recently, this was done for the hexagon. Presently, we successfully perform resummation for the seven-leg tree NMHV amplitude. To this end, we construct the flux-tube integrands of the fifteen independent Grassmann component of the heptagon with an infinite number of small fermion-antifermion pairs accounted for in NMHV two-channel conformal blocks.


title: A robust inverse scattering transform for the focusing nonlinear Schr\"odinger equation

abstract: We propose a modification of the standard inverse scattering transform for the focusing nonlinear Schr\"odinger equation (also other equations by natural generalization) formulated with nonzero boundary conditions at infinity. The purpose is to deal with arbitrary-order poles and potentially severe spectral singularities in a simple and unified way. As an application, we use the modified transform to place the Peregrine solution and related higher-order "rogue wave" solutions in an inverse-scattering context for the first time. This allows one to directly study properties of these solutions such as their dynamical or structural stability, or their asymptotic behavior in the limit of high order. The modified transform method also allows rogue waves to be generated on top of other structures by elementary Darboux transformations, rather than the generalized Darboux transformations in the literature or other related limit processes.


title: Liouville type theorems on the steady Navier-Stokes equations in R3

abstract: In this paper we study the Liouville type properties for solutions to the steady incompressible Navier-Stoks equations in $\mathbf{R}^{3}$. It is shown that any solution to the steady Navier-Stokes equations in $\mathbf{R}^{3}$ with finite Dirichlet integral and vanishing velocity field at far fields must be trivial. This solves an open problem. The key ingredients of the proof include a Hodge decomposition of the energy-flux and the observation that the square of the deformation matrix lies in the local Hardy space. As a by-product, we also obtain a Liouville type theorem for the steady density-dependent Navier-Stokes equations.


title: A Correspondence Between Random Neural Networks and Statistical Field Theory

abstract: A number of recent papers have provided evidence that practical design questions about neural networks may be tackled theoretically by studying the behavior of random networks. However, until now the tools available for analyzing random neural networks have been relatively ad-hoc. In this work, we show that the distribution of pre-activations in random neural networks can be exactly mapped onto lattice models in statistical physics. We argue that several previous investigations of stochastic networks actually studied a particular factorial approximation to the full lattice model. For random linear networks and random rectified linear networks we show that the corresponding lattice models in the wide network limit may be systematically approximated by a Gaussian distribution with covariance between the layers of the network. In each case, the approximate distribution can be diagonalized by Fourier transformation. We show that this approximation accurately describes the results of numerical simulations of wide random neural networks. Finally, we demonstrate that in each case the large scale behavior of the random networks can be approximated by an effective field theory.


title: Entropy-bounded solutions to the compressible Navier-Stokes equations: with far field vacuum

abstract: The entropy is one of the fundamental states of a fluid and, in the viscous case, the equation that it satisfies is highly singular in the region close to the vacuum. In spite of its importance in the gas dynamics, the mathematical analyses on the behavior of the entropy near the vacuum region, were rarely carried out; in particular, in the presence of vacuum, either at the far field or at some isolated interior points, it was unknown if the entropy remains its boundedness. The results obtained in this paper indicate that the ideal gases retain their uniform boundedness of the entropy, locally or globally in time, if the vacuum occurs at the far field only and the density decays slowly enough at the far field. Precisely, we consider the Cauchy problem to the one-dimensional full compressible Navier-Stokes equations without heat conduction, and establish the local and global existence and uniqueness of entropy-bounded solutions, in the presence of vacuum at the far field only. It is also shown that, different from the case that with compactly supported initial density, the compressible Navier-Stokes equations, with slowly decaying initial density, can propagate the regularities in the inhomogeneous Sobolev spaces.


title: FIMP Dark Matter Freeze-in Gauge Mediation and Hidden Sector

abstract: We explore the dark matter freeze-in mechanism within the gauge mediation framework, which involves a hidden Feebly Interacting Massive Particle (FIMP) coupling feebly with the messenger fields while the messengers are still in the thermal bath. The FIMP is the fermionic component of the pseudo-moduli in a generic metastable supersymmetry (SUSY) breaking model and resides in the hidden sector. The relic abundance and the mass of the FIMP are determined by the SUSY breaking scale and the feeble coupling. The gravitino, which is the canonical dark matter candidate in the gauge mediation framework, contributes to the dark matter relic abundance along with the freeze-in of the FIMP. The hidden sector thus becomes two- component with both the FIMP and gravitino lodging in the SUSY breaking hidden sector. We point out that the ratio between the FIMP and the gravitino is determined by how SUSY breaking is communicated to the messengers. In particular when the FIMP dominates the hidden sector, the gravitino becomes the minor contributor in the hidden sector. Meanwhile, the neutralino is assumed to be both the Weakly Interacting Massive Particle (WIMP) dark matter candidate in the freeze-out mechanism and the Lightest Observable SUSY Particle (LOSP). We further find out the neutralino has the sub-leading contribution to the current dark matter relic density in the parameter space of our freeze-in gauge mediation (FIGM) model. Our result links the SUSY breaking scale in the gauge mediation framework with the FIMP freeze-in production rate leading to a natural and predicting scenario for the studies of the dark matter in the hidden sector.


title: Towards an unified theory for testing statistical hypothesis: Multinormal mean with nuisance covariance matrix

abstract: Under a multinormal distribution with arbitrary unknown covariance matrix, the main purpose of this paper is to propose a framework to achieve the goal of reconciliation of Bayesian, frequentist and Fisherian paradigms for the problems of testing mean against restricted alternatives (closed convex cones). Combining Fisher's fiducial inference and Wald's decision theory via d-admissibility into an unified approach, the goal can then be achieved. To proceed, the tests constructed via the union-intersection principle are studied.


title: The Effects of Memory Replay in Reinforcement Learning

abstract: Experience replay is a key technique behind many recent advances in deep reinforcement learning. Allowing the agent to learn from earlier memories can speed up learning and break undesirable temporal correlations. Despite its wide-spread application, very little is understood about the properties of experience replay. How does the amount of memory kept affect learning dynamics? Does it help to prioritize certain experiences? In this paper, we address these questions by formulating a dynamical systems ODE model of Q-learning with experience replay. We derive analytic solutions of the ODE for a simple setting. We show that even in this very simple setting, the amount of memory kept can substantially affect the agent's performance. Too much or too little memory both slow down learning. Moreover, we characterize regimes where prioritized replay harms the agent's learning. We show that our analytic solutions have excellent agreement with experiments. Finally, we propose a simple algorithm for adaptively changing the memory buffer size which achieves consistently good empirical performance.


title: Polarized bound state in the continuum and resonances with tunable Q-factor in an anisotropic photonic crystal

abstract: We consider a one-dimensional photonic crystal composed of alternating layers of isotropic and anisotropic dielectric materials. Such a system has different band structures for different polarizations of light. We demonstrate that if an anisotropic defect layer is inserted into the structure, the crystal can support an optical bound state in the continuum. By tilting the principle dielectric axes of the defect layer relative to those of the photonic crystal we observe a long-lived resonance in the transmission spectrum. We derive an analytical expression for the decay rate of the resonance that agrees well with the numerical data by the Berreman anisotropic transfer matrix approach. An experimental set-up with a liquid crystal defect layer is proposed to tune the Q-factor of the resonance through applying an external electric field. We speculate that the set-up provides a simple and robust platform for observing optical bound states in the continuum in the form of resonances with tunable Q-factor.


title: Dirac-mode analysis for quark number density and its application for deconfinement transition

abstract: The quark number density at finite imaginary chemical potential is investigated in the lattice QCD using the Dirac-mode expansion. We find the analytical formula of the quark number density in terms of the Polyakov loop in the large quark mass regime. On the other hand, in the small quark mass region, the quark number density is investigated by using the quenched lattice QCD simulation. The quark number density is found to strongly depend on the low-lying Dirac modes while its sign does not change. This result leads to that the quark number holonomy is not sensitive to the low-lying Dirac modes. We discuss the confinement-deconfinement transition from the property of the quark number density and the quark number holonomy.


title: Monogamy relations of concurrence for any dimensional quantum systems

abstract: We study monogamy relations for arbitrary dimensional multipartite systems. Monogamy relations based on concurrence and concurrence of assistance for any dimensional $m_1\otimes m_2\otimes...\otimes m_{N}$ quantum states are derived, which give rise to the restrictions on the entanglement distributions among the subsystems. Besides, we give the lower bound of concurrence for four-partite mixed states. The approach can be readily generalized to arbitrary multipartite systems.


title: Acceleration of Gradient-based Path Integral Method for Efficient Optimal and Inverse Optimal Control

abstract: This paper deals with a new accelerated path integral method, which iteratively searches optimal controls with a small number of iterations. This study is based on the recent observations that a path integral method for reinforcement learning can be interpreted as gradient descent. This observation also applies to an iterative path integral method for optimal control, which sets a convincing argument for utilizing various optimization methods for gradient descent, such as momentum-based acceleration, step-size adaptation and their combination. We introduce these types of methods to the path integral and demonstrate that momentum-based methods, like Nesterov Accelerated Gradient and Adam, can significantly improve the convergence rate to search for optimal controls in simulated control systems. We also demonstrate that the accelerated path integral could improve the performance on model predictive control for various vehicle navigation tasks. Finally, we represent this accelerated path integral method as a recurrent network, which is the accelerated version of the previously proposed path integral networks (PI-Net). We can train the accelerated PI-Net more efficiently for inverse optimal control with less RAM than the original PI-Net.


title: Imaginary potential of moving quarkonia in a D-instanton background

abstract: The imaginary part of the inter-quarks potential of moving heavy quarkonia is investigated in a dual supergravity of the AdS background deformed by dilaton, which induces the gauge field condensate in the dual gauge theory. We analyze the quark anti-quark pair moving transverse and parallel to the plasma wind, respectively. It is shown that for both cases the presence of D-instanton density tends to increase the inter-distance and decrease the imaginary potential, reverse to the effect of the pair's velocity. Moreover, it is found that the D-instanton density has stronger effects for the parallel case rather than transverse.


title: Symmetry-unprotected nodes or gap minima in $s_{++}$ state of FeSe single crystal

abstract: We report the study on superconducting pairing mechanism of FeSe via the pair-breaking effects induced by H$^+$-irradiation combined with low-temperature specific heat measurements. A multi-gap structure with nodes or gap minima is suggested in a clean FeSe by the specific heat results. The suppression of critical temperature $T_c$ with increasing the defect density manifests a two-step behavior. When the increase in the residual resistivity is small, $\Delta\rho_0$ $<$ $\sim$4.5 $\mu\Omega$cm, $T_c$ is gradually suppressed with increasing the density of scattering centers, suggesting the presence of symmetry-unprotected nodes or gap minima. However, for $\Delta\rho_0$ $>$ $\sim$4.5 $\mu\Omega$cm, $T_c$ is almost independent of the scattering, which indicates that the nodes or gap minima are lifted and the order parameter becomes almost isotropic without sign change. Thus, the superconductivity in FeSe is found to be realized in symmetry-unprotected nodal or highly anisotropic $s_{++}$ state.


title: Studies of Single Component Fermi Gas near a $P$-wave Resonance with Lowest Order Constrained Variational Method

abstract: We study single component Fermi gas near a $p$-wave resonance with the lowest order constrained variational (LOCV) method. We obtain the energy per particle for the ground state of single component Fermi gas near a $p$-wave resonance with LOCV method. We also calculate compressibility of the single component Fermi gas near a $p$-wave resonance and it shows that near the $p$-wave resonance, the system would lose its stability and collapse. The two $p$-wave contacts are also obtained and their variation tendencies with interaction strength are consistent with recent experimental results.


title: Learning Social Image Embedding with Deep Multimodal Attention Networks

abstract: Learning social media data embedding by deep models has attracted extensive research interest as well as boomed a lot of applications, such as link prediction, classification, and cross-modal search. However, for social images which contain both link information and multimodal contents (e.g., text description, and visual content), simply employing the embedding learnt from network structure or data content results in sub-optimal social image representation. In this paper, we propose a novel social image embedding approach called Deep Multimodal Attention Networks (DMAN), which employs a deep model to jointly embed multimodal contents and link information. Specifically, to effectively capture the correlations between multimodal contents, we propose a multimodal attention network to encode the fine-granularity relation between image regions and textual words. To leverage the network structure for embedding learning, a novel Siamese-Triplet neural network is proposed to model the links among images. With the joint deep model, the learnt embedding can capture both the multimodal contents and the nonlinear network information. Extensive experiments are conducted to investigate the effectiveness of our approach in the applications of multi-label classification and cross-modal search. Compared to state-of-the-art image embeddings, our proposed DMAN achieves significant improvement in the tasks of multi-label classification and cross-modal search.


title: Distributed economic control of dynamically coupled networks

abstract: This paper investigates the synthesis of distributed economic control algorithms under which dynamically coupled physical systems are regulated to a variational equilibrium of a constrained convex game. We study two complementary cases: (i) each subsystem is linear and controllable; and (ii) each subsystem is nonlinear and in the strict-feedback form. The convergence of the proposed algorithms is guaranteed using Lyapunov analysis. Their performance is verified by two case studies on a multi-zone building temperature regulation problem and an optimal power flow problem, respectively.


title: Collective resonances near zero energy induced by a point defect in bilayer graphene

abstract: Intrinsic defects give rise to scattering processes governing the transport properties of mesoscopic systems. We investigate analytically and numerically the local density of states in Bernal stacking bilayer graphene with a point defect. With Bernal stacking structure, there are two types of lattice sites. One corresponds to connected sites, where carbon atoms from each layer stack on top of each other, and the other corresponds to disconnected sites. From our theoretical study, a picture emerges in which the pronounced zero-energy peak in the local density of states does not attribute to zero-energy impurity states associated to two different types of defects but to a collective phenomenon of the low-energy resonant states induced by the defect. To corroborate this description, we numerically show that at small system size $N$, where $N$ is the number of unit cells, the zero-energy peak near the defect scales as $1/\ln N$ for the quasi-localized zero-energy state and as $1/N$ for the delocalized zero-energy state. As the system size approaches to the thermodynamic limit, the former zero-energy peak becomes a power-law singularity $1/|E|$ in low energies, while the latter is broadened into a Lorentzian shape. A striking point is that both types of zero-energy peaks decay as $1/r^2$ away from the defect, manifesting the quasi-localized character. Based on our results, we propose a general formula for the local density of states in low-energy and in real space. Our study sheds light on this fundamental problem of defects.


title: Suppressing chemotactic blow-up through a fast splitting scenario on the plane

abstract: We revisit the question of global regularity for the Patlak-Keller-Segel (PKS) chemotaxis model. The classical 2D hyperbolic-elliptic model blows up for initial mass M>8\pi. We consider more realistic scenario which takes into account the flow of the ambient environment induced by harmonic potentials, and thus retain the identical elliptic structure as in the original PKS. Surprisingly, we find that already the simplest case of linear stationary vector field, $Ax^\top$, with large enough amplitude $A$, prevents chemotactic blow-up. Specifically, the presence of such an ambient fluid transport creates what we call a 'fast splitting scenario', which competes with the focusing effect of aggregation so that 'enough mass' is pushed away from concentration along the $x_1$-axis, thus avoiding a finite time blow-up, at least for M<16\pi. Thus, the enhanced ambient flow doubles the amount of allowable mass which evolve to global smooth solutions.


title: $F_\gamma$: a new observable for photon-hadron discrimination in hybrid air shower events

abstract: To search for ultra-high-energy photons in primary cosmic rays, air shower observables are needed that allow a good separation between primary photons and primary hadrons. We present a new observable, $F_\gamma$, which can be extracted from ground-array data in hybrid events, where simultaneous measurements of the longitudinal and the lateral shower profile are performed. The observable is based on a template fit to the lateral distribution measured by the ground array with the template taking into account the complementary information from the measurement of the longitudinal profile, i.e. the primary energy and the geometry of the shower. $F_\gamma$ shows a very good photon-hadron separation, which is even superior to the separation given by the well-known $X_\text{max}$ observable (the atmospheric depth of the shower maximum). At energies around $1\,\text{EeV}$ ($10\,\text{EeV}$), $F_\gamma$ provides a background rejection better than $97.8\,\%$ ($99.9\,\%$) at a signal efficiency of $50\,\%$. Advantages of the observable $F_\gamma$ are its technical stability with respect to irregularities in the ground array (i.e. missing or temporarily non-operating stations) and that it can be applied over the full energy range accessible to the air shower detector, down to its threshold energy. Furthermore, $F_\gamma$ complements nicely to $X_\text{max}$ such that both observables can well be combined to achieve an even better discrimination power, exploiting the rich information available in hybrid events.


title: Association and Load Optimization with User Priorities in Load-Coupled Heterogeneous Networks

abstract: In this paper, we consider the network utility maximization problem with various user priorities via jointly optimizing user association, load distribution and power control in a load-coupled heterogeneous network. In order to tackle the nonconvexity of the problem, we first analyze the problem by obtaining the optimal resource allocation strategy in closed form and characterizing the optimal base station load distribution pattern. Both observations are shown essential in simplifying the original problem and making it possible to transform the nonconvex load distribution and power control problem into convex reformulation via exponential variable transformation. An iterative algorithm with low complexity is accordingly presented to obtain a suboptimal solution to the joint optimization problem. Simulation results show that the proposed algorithm achieves better performance than conventional approaches.


title: Empirical correlation between the interfacial Dzyaloshinskii Moriya interaction and work function in metallic magnetic trilayers

abstract: The Dzyaloshinskii Moriya interaction (DMI) generates intriguing chiral magnetic objects such as magnetic skyrmions and chiral domain walls that can be used as building blocks in emerging magnetic nanodevices. To achieve better stability and functionality of these chiral objects, it is essential to achieve a larger DMI. In this paper, we report an experimental observation that in magnetic trilayer films, the DMI strength is mainly determined by the work functions of the nonmagnetic layers interfaced with the magnetic layer. The clear correlation with the intrinsic material parameters provides a guideline for material selection to engineer the DMI strength.


title: Anomalous Behaviour of Mutual Information in Finite Flocks

abstract: The existing consensus is that flocks are poised at criticality, entailing long correlation lengths and a maximal value of Shannon mutual information in the large-system limit. We show, by contrast, that for finite flocks in the long observation time limit, mutual information may not only fail to peak at criticality---as observed for other critical systems---but also diverge as noise tends to zero. This result carries implications for other finite-size, out-of-equilibrium systems, where observation times may vary widely compared to time scales of internal system dynamics; thus it may not be assumed that mutual information locates the phase transition.


title: MEDOC: a Python wrapper to load MEDLINE into a local MySQL database

abstract: Since the MEDLINE database was released, the number of documents indexed by this entity has risen every year. Several tools have been developed by the National Institutes of Health (NIH) to query this corpus of scientific publications. However, in terms of advances in big data, text-mining and data science, an option to build a local relational database containing all metadata available on MEDLINE would be truly useful to optimally exploit these resources. MEDOC (MEdline DOwnloading Contrivance) is a Python program designed to download data on an FTP and to load all extracted information into a local MySQL database. It took MEDOC 4 days and 17 hours to load the 26 million documents available on this server onto a standard computer. This indexed relational database allows the user to build complex and rapid queries. All fields can thus be searched for desired information, a task that is difficult to accomplish through the PubMed graphical interface. MEDOC is free and publicly available at https://github.com/MrMimic/MEDOC.


title: Review of Data Structures for Computationally Efficient Nearest-Neighbour Entropy Estimators for Large Systems with Periodic Boundary Conditions

abstract: Information theoretic quantities are extremely useful in discovering relationships between two or more data sets. One popular method---particularly for continuous systems---for estimating these quantities is the nearest neighbour estimators. When system sizes are very large or the systems have periodic boundary conditions issues with performance and correctness surface, however solutions are known for each problem. Here we show that these solutions are inappropriate in systems that simultaneously contain both features and discuss a lesser known alternative solution involving Vantage Point trees that is capable of addressing both issues.


title: Eigenvalue fluctuations for lattice Anderson Hamiltonians: Unbounded potentials

abstract: We consider random Schr\"odinger operators with Dirichlet boundary conditions outside lattice approximations of a smooth Euclidean domain and study the behavior of its lowest-lying eigenvalues in the limit when the lattice spacing tends to zero. Under a suitable moment assumption on the random potential and regularity of the spatial dependence of its mean, we prove that the eigenvalues of the random operator converge to those of a deterministic Schr\"odinger operator. Assuming also regularity of the variance, the fluctuation of the random eigenvalues around their mean are shown to obey a multivariate central limit theorem. This extends the authors' recent work where similar conclusions have been obtained for bounded random potentials.


title: Detecting Partial Rosettes in Tumor Histopathology Using The Cross Product

abstract: Tumors of the eye and nervous system often exhibit an arch-like arrangement of nuclei, called rosettes. Pathologists are able to identify rosettes [full circles] and the presence of partial rosettes [semi-circles] and interpret this as a sign of differentiation in a tumor. However, there is no objective method to quantitate the many partial rosettes that are obvious or not obvious to the naked eye. This paper proposes a mathematical algorithm to computationally detect the presence of obvious or non-obvious partial rosettes, henceforth referred to as Nguyen-Wu Partial Rosettes. Quantifying the degree of partial rosettes present in a tumor may allow pathologists to stratify tumors into more refined groups that may respond better to therapy or have different clinical outcomes. The Midline Cross Product [MCP] algorithm calculates the magnitude of two cross products and adds them together to obtain one value. Each of the two cross products results from [1] the line that connects the midpoints of longest lengths of two neighboring ovals, and [2] the line that extends from the midpoint from one longest length and is perpendicular to that longest length. The MCP algorithm makes nuclei that are arranged in consecutive rows and arches quantitatively distinct from nuclei that are arranged next to each other in a disorderly manner.


title: On the motivic Peterson conjecture

abstract: We show that the analogue of the Peterson conjecture on the action of Steenrod squares does not hold in motivic cohomology.


title: Variational Inference based on Robust Divergences

abstract: Robustness to outliers is a central issue in real-world machine learning applications. While replacing a model to a heavy-tailed one (e.g., from Gaussian to Student-t) is a standard approach for robustification, it can only be applied to simple models. In this paper, based on Zellner's optimization and variational formulation of Bayesian inference, we propose an outlier-robust pseudo-Bayesian variational method by replacing the Kullback-Leibler divergence used for data fitting to a robust divergence such as the beta- and gamma-divergences. An advantage of our approach is that superior but complex models such as deep networks can also be handled. We theoretically prove that, for deep networks with ReLU activation functions, the \emph{influence function} in our proposed method is bounded, while it is unbounded in the ordinary variational inference. This implies that our proposed method is robust to both of input and output outliers, while the ordinary variational method is not. We experimentally demonstrate that our robust variational method outperforms ordinary variational inference in regression and classification with deep networks.


title: The LifeV library: engineering mathematics beyond the proof of concept

abstract: LifeV is a library for the finite element (FE) solution of partial differential equations in one, two, and three dimensions. It is written in C++ and designed to run on diverse parallel architectures, including cloud and high performance computing facilities. In spite of its academic research nature, meaning a library for the development and testing of new methods, one distinguishing feature of LifeV is its use on real world problems and it is intended to provide a tool for many engineering applications. It has been actually used in computational hemodynamics, including cardiac mechanics and fluid-structure interaction problems, in porous media, ice sheets dynamics for both forward and inverse problems. In this paper we give a short overview of the features of LifeV and its coding paradigms on simple problems. The main focus is on the parallel environment which is mainly driven by domain decomposition methods and based on external libraries such as MPI, the Trilinos project, HDF5 and ParMetis. Dedicated to the memory of Fausto Saleri.


title: Fingerprinting of Compound-Specific Chlorine Isotopologue Distribution of Organochlorines by GC-HRMS

abstract: This study developed a method to measure compound-specific chlorine isotopologue distribution of organochlorines for source identification and apportionment. Complete chlorine isotopologues of individual model organochlorines were detected by gas chromatography-high resolution mass spectrometry (GC-HRMS). The measured relative abundances (RAmea), simulated relative abundances (RAsim), and relative variations between RAmea and RAsim (${\Delta}$RA) were obtained on basis of the detected MS signal intensities of individual isotopologues. The method has been partially validated in precision, injection-amount dependency and temporal drifts. The standard deviations (SDs) of RAmea of all istotopologues of perchlorethylene (PCE) and trichloroethylene (TCE) from different manufacturers were 0.002%-0.069%. The SDs of ${\Delta}$RA of the first three isotopologues of PCE and TCE from different manufacturers were 0.026%-0.155%. The ${\Delta}$RA and ${\Delta}$RA patterns of the standards of PCE and TCE from different manufacturers were able to be differentiated with statistical significance. The method has been successfully applied to an analogous case of source identification and apportionment for two trichlorodiphenyls that exhibited significant chlorine isotope fractionation on the GC-HRMS system. The results demonstrate that the ${\Delta}$RA in conjunction with isotope ratios can be applied to source identification and present high-validity identification outcomes. Moreover, the RAmea can be used to implement source apportionment for organochlorines from more sources with more reliable outcomes compared with the methods using isotope ratios only. This method opens a new way to perform fingerprinting analysis of compound-specific chlorine isotopologue distribution of organochlorines, and will be a promising high-performance approach in source delineation and apportionment for chlorinated organic compounds.


title: A Bounded Degree Lasserre Hierarchy with SOCP Relaxations for Global Polynomial Optimization and Applications

abstract: In this paper, we propose a new convergent conic programming hierarchy of relaxations involving both semi-definite cone and second-order cone constraints for solving nonconvex polynomial optimization problems to global optimality. The significance of this hierarchy is that the size and number of the semi-definite and second-order cone constraints of the relaxations are fixed and independent of the step or level of the approximation in the hierarchy. Using the Krivine-Stengle's certificate of positivity in real algebraic geometry, we establish the convergence of the hierarchy of relaxations, extending the very recent so-called bounded degree Lasserre hierarchy. In particular, we also provide a convergent bounded degree second-order cone programming (SOCP) hierarchy for solving polynomial optimization problems. We then present finite convergence at step one of the SOCP hierarchy for two classes of polynomial optimization problems: a subclass of convex polynomial optimization problems where the objective and constraint functions are SOCP-convex polynomials, defined in terms of specially structured sum of squares polynomials, and a class of polynomial optimization problems, involving polynomials with essentially non-positive coefficients. In the case of one-step convergence for problems with SOCP-convex polynomials, we show how a global solution is recovered from the relaxation via Jensen's inequality of SOCP-convex polynomials. As an application, we derive a corresponding convergent conic linear programming hierarchy for conic-convex semi-algebraic programs. Whenever the semi-algebraic set of the conic-convex program is described by convex polynomial inequalities, we show further that the values of the relaxation problems converge to the common value of the convex program and its Lagrangian dual under a constraint qualification.


title: A history of violence: insights into post-accretionary heating in carbonaceous chondrites from volatile element abundances, Zn isotopes, and water contents

abstract: Carbonaceous chondrites (CCs) may have been the carriers of water, volatile and moderately volatile elements to Earth. Investigating the abundances of these elements, their relative volatility, and isotopes of state-change tracer elements such as Zn, and linking these observations to water contents, provide vital information on the processes that govern the abundances and isotopic signatures of these species in CCs and other planetary bodies. Here we report Zn isotopic data for 28 CCs (20 CM, 6 CR, 1 C2-ung, and 1 CV3), as well as trace element data for Zn, In, Sn, Tl, Pb, and Bi in 16 samples (8 CM, 6 CR, 1 C2-ung, and 1 CV3), that display a range of elemental abundances from case-normative to intensely depleted. We use these data, water content data from literature and Zn isotopes to investigate volatile depletions and to discern between closed and open system heating. Trace element data have been used to construct relative volatility scales among the elements for the CM and CR chondrites. From least volatile to most, the scale in CM chondrites is Pb-Sn-Bi-In-Zn-Tl, and for CR chondrites it is Tl-Zn-Sn-Pb-Bi-In. These observations suggest that heated CM and CR chondrites underwent volatile loss under different conditions to one another and to that of the solar nebula, e.g. differing oxygen fugacities. Furthermore, the most water and volatile depleted samples are highly enriched in the heavy isotopes of Zn. Taken together, these lines of evidence strongly indicate that heated CM and CR chondrites incurred open system heating, stripping them of water and volatiles concomitantly, during post-accretionary shock impact(s).


title: Quantifying the AGN-driven outflows in ULIRGs (QUADROS) I: VLT/Xshooter observations of 9 nearby objects

abstract: Although now routinely incorporated into hydrodynamic simulations of galaxy evolution, the true importance of the feedback effect of the outflows driven by active galactic nuclei (AGN) remains uncertain from an observational perspective. This is due to a lack of accurate information on the densities, radial scales and level of dust extinction of the outflow regions. Here we use the unique capabilities of VLT/Xshooter to investigate the warm outflows in a representative sample of 9 local (0.06 $<$ z $<$ 0.15) ULIRGs with AGN nuclei and, for the first time, accurately quantify the key outflow properties. We find that the outflows are compact (0.05 $<$ R$_{[OIII]}$ $<$ 1.2 kpc), significantly reddened (median E(B-V)$\sim$0.5 magnitudes), and have relatively high electron densities (3.4 $<$ log$_{10}$ n$_e$ (cm$^{-3}$) $<$ 4.8). It is notable that the latter densities -- obtained using trans-auroral [SII] and [OII] emission-line ratios -- exceed those typically assumed for the warm, emission-line outflows in active galaxies, but are similar to those estimated for broad and narrow absorption line outflow systems detected in some type 1 AGN. Even if we make the most optimistic assumptions about the true (deprojected) outflow velocities, we find relatively modest mass outflow rates ($0.07 < \dot{M} < 11$ M$_{sol}$ yr$^{-1}$) and kinetic powers measured as a fraction of the AGN bolometric luminosities ($4\times10^{-4} < \dot{E}/L_{BOL} <1$%). Therefore, although warm, AGN-driven outflows have the potential to strongly affect the star formation histories in the inner bulge regions ($r \sim$ 1kpc) of nearby ULIRGs, we lack evidence that they have a significant impact on the evolution of these rapidly evolving systems on larger scales.


title: Strong Coupling of 3D Cavity Photons to Travelling Magnons At Low Temperatures

abstract: We demonstrate strong coupling between travelling magnons in an Yttrium Iron Garnet film and 3D microwave cavity photons at milli-Kelvin temperatures. The coupling strength of $350$MHz or $7.3$\% of resonance frequency is observed. The magnonic subsystem is represented by the Damon-Eshbach magnetostatic surface wave with a distribution of wave numbers giving the linewidth of 15MHz. The ways to improve this parameter are discussed. The energy gap in the spectrum given by the Zeeman energy and the shape-anisotropy energy in the film geometry give rise to a significant asymmetry of the double peak structure of the photon-magnon avoided level crossing. A structure of two parallel YIG films is investigated using the same re-entrant magnetostatic surface wave transducer revealing a higher order magnon modes existing in both films. Combination of a multi-post re-entrant cavity and multiple films is a potential base for engineering both magnon and photon spectra.


title: Evolution of the Solar Lyman-Alpha line profile during the solar cycle

abstract: Recent studies of interstellar neutral (ISN) hydrogen observed by the Interstellar Boundary Explorer (IBEX) suggested that the present understanding of the radiation pressure acting on hydrogen atoms in the heliosphere should be revised. There is a significant discrepancy between theoretical predictions of the ISN H signal using the currently used model of the solar Lyman-alpha profile by Tarnopolski et al. 2009 (TB09) and the signal due to ISN H observed by IBEX-Lo. We developed a new model of evolution of the solar Lyman-alpha profile that takes into account all available observations of the full-disk solar Lyman-alpha profiles from SUMER/SOHO, provided by Lemaire et al. 2015 (L15), covering practically the entire 23rd solar cycle. The model has three components that reproduce different features of the profile. The main shape of the emission line that is produced in the chromosphere is modeled by the kappa function; the central reversal due to absorption in the transition region is modeled by the Gauss function; the spectral background is represented by the linear function. The coefficients of all those components are linear functions of the line-integrated full-disk Lyman-alpha irradiance, which is the only free parameter of the model. The new model features potentially important differences in comparison with the model by TB09, which was based on a limited set of observations. This change in the understanding of radiation pressure, especially during low solar activity, may significantly affect the interstellar H and D distributions in the inner heliosphere and their derivative populations.


title: Self-consistent one-electron equation including exchange and correlation in terms of equivalent function and phase norm: calculation for ground and excited states in a unified way

abstract: We present an alternative one-electron equation for resolving many-electron problem to one-electron approximation and including the exchange and correlation effects in an analytical way, thereby fulfilling the requirements for ab initio calculation. To derive one-electron equation, we accept a new notion of equivalent function suggestive of the pseudo wavefunction. As a result, we reduce many-electron equation to one-electron including the exchange effect in an analytical method. Moreover we accept the notion of phase norm for two electrons to resolve the electronic correlation problem. The phase norm is used to specify the electron-approachable limit between particles. We take into consideration the electronic correlation with the help of a correlation-hole function in terms of the phase norm, by which multiplying the integrand of the operator term representing the interaction between electrons. Using the phase norm leads to analytical consideration of the electronic correlation without employing in a factitious way the additional term pertaining to correlation, so it embodies the physical essence of electronic correlation. The derived equation becomes an one-electron equation which does not include an additional term pertaining to the exchange and correlation, but takes into consideration the exchange and correlation effects in a rigorous ab initio way.


title: A novel type of rogue waves with predictability in nonlinear physics

abstract: Rogue waves named by oceanographers are ubiquitous in nature and appear in a variety of different contexts such as water waves, liquid Helium, nonlinear optics, microwave cavities, etc. In this letter, we propose a novel type of exact (2+1)-dimensional rogue waves which may be found in many physical fields described by integrable and nonintegrable models. This type of rogue waves are closely related to invisible lumps. Usually, a lump is an algebraically localized wave in space but visible at any time. Whence a lump induces a soliton, the lump will become invisible before or after a fixed time. If a bounded twin soliton is induced by a lump, the lump will become a rogue wave (or instanton) and can only be visible at an instant time. Because of the existence of the induced visible solitons, the rogue wave may be predictable in some senses. For instance, the height, the position and the arrival time of the rogue wave can be predictable. The idea is illustrated by the celebrate (2+1)-dimensional Kadomtsev-Petviashvili equation in its extended form.


title: Linear Gravity in de Sitter universe

abstract: We give in this paper an explicit construction of the covariant quantization of the rank-two massless tensor field on de Sitter space (linear covariant quantum gravity on a de Sitter background). The main ingredient of the construction is an indecomposable representation of de Sitter group. We here make the choice of a specific simple gauge fixing. We show that our gauge fixing eliminates any infrared divergence in the two-point function for the traceless part of this field. But it is not possible to do the same for the pure trace part (conformal sector). We describe the related Krein space structure and covariant field operators. This work is in the continuation of our previous ones concerning the massless minimally coupled scalar fields and the massive tensor field on de Sitter.


title: Cooperative Couplings between Octahedral Rotations and Ferroelectricity in Perovskites

abstract: The structure of ABO3 perovskites is dominated by two types of unstable modes, namely, the oxygen octahedral rotation (AFD) and ferroelectric (FE) mode. It is generally believed that such AFD and FE modes tend to compete and suppress each other. Here we use first-principles methods to show that a dual nature of the AFD-FE coupling, which turns from competitive to cooperative as the AFD mode strengthens, occurs in numerous perovskite oxides. We provide a unified model of such a dual interaction by introducing novel high-order coupling terms, and explain the atomistic origin of the resulting new form of ferroelectricity in terms of universal steric mechanisms. We also predict that such a novel form of ferroelectricity leads to atypical behaviors, such as an enhancement of the electric polarization under hydrostatic pressure.


title: Status of the continuous gravitational wave searches in the Advanced Detector Era

abstract: Periodic (almost monochromatic) gravitational waves emitted by rotating, asymmetric neutron stars are intriguing potential signals in the sensitivity band of Advanced LIGO and Advanced Virgo detectors. These signals are related to elastic and magnetic stresses in the neutron-star interior, as well as to various possible instabilities, and thus are interesting from the point of view of the largely-unknown neutron star structure. I will describe the main challenges related to these searches, the current state of the data-analysis methods and plans for the future.


title: Cell Segmentation in 3D Confocal Images using Supervoxel Merge-Forests with CNN-based Hypothesis Selection

abstract: Automated segmentation approaches are crucial to quantitatively analyze large-scale 3D microscopy images. Particularly in deep tissue regions, automatic methods still fail to provide error-free segmentations. To improve the segmentation quality throughout imaged samples, we present a new supervoxel-based 3D segmentation approach that outperforms current methods and reduces the manual correction effort. The algorithm consists of gentle preprocessing and a conservative super-voxel generation method followed by supervoxel agglomeration based on local signal properties and a postprocessing step to fix under-segmentation errors using a Convolutional Neural Network. We validate the functionality of the algorithm on manually labeled 3D confocal images of the plant Arabidopis thaliana and compare the results to a state-of-the-art meristem segmentation algorithm.


title: Supervised and Extended Restart in Random Walks for Ranking and Link Prediction in Networks

abstract: Given a real-world graph, how can we measure relevance scores for ranking and link prediction? Random walk with restart (RWR) provides an excellent measure for this and has been applied to various applications such as friend recommendation, community detection, anomaly detection, etc. However, RWR suffers from two problems: 1) using the same restart probability for all the nodes limits the expressiveness of random walk, and 2) the restart probability needs to be manually chosen for each application without theoretical justification. We have two main contributions in this paper. First, we propose Random Walk with Extended Restart (RWER), a random walk based measure which improves the expressiveness of random walks by using a distinct restart probability for each node. The improved expressiveness leads to superior accuracy for ranking and link prediction. Second, we propose SuRe (Supervised Restart for RWER), an algorithm for learning the restart probabilities of RWER from a given graph. SuRe eliminates the need to heuristically and manually select the restart parameter for RWER. Extensive experiments show that our proposed method provides the best performance for ranking and link prediction tasks, improving the MAP (Mean Average Precision) by up to 15.8% on the best competitor.


title: Can test fields destroy the event horizon in the Kerr-Taub-NUT spacetime?

abstract: In this work we investigate if the interaction of the Kerr-Taub-NUT spacetime with test scalar and neutrino fields can lead to the destruction of the event horizon. It turns out that both extremal and nearly extremal black holes can be destroyed by scalar and neutrino fields if the initial angular momentum of the spacetime is sufficiently large relative to its mass and NUT charge. This is the first example in which a classical field satisfying the null energy condition can actually destroy an extremal black hole. For scalar fields, the modes that can lead to the destruction of the horizon are restricted to a narrow range due to superradiance. Since superradiance does not occur for neutrino fields, the destruction of the horizon by neutrino fields is generic, and it cannot be fixed by backreaction effects. We also show that the extremal black holes that can be destroyed by scalar fields correspond to naked singularities in the Kerr limit, in accord with the previous results which imply that extremal Kerr black holes cannot be destroyed by scalar test fields.


title: On community structure validation in real networks

abstract: Community structure is a commonly observed feature of real networks. The term refers to the presence in a network of groups of nodes (communities) that feature high internal connectivity, but are poorly connected between each other. Whereas the issue of community detection has been addressed in several works, the problem of validating a partition of nodes as a good community structure for a real network has received considerably less attention and remains an open issue. We propose a set of indices for community structure validation of network partitions that are based on an hypothesis testing procedure that assesses the distribution of links between and within communities. Using both simulations and real data, we illustrate how the proposed indices can be employed to compare the adequacy of different partitions of nodes as community structures in a given network, to assess whether two networks share the same or similar community structures, and to evaluate the performance of different network clustering algorithms.


title: Mirror Descent and Convex Optimization Problems With Non-Smooth Inequality Constraints

abstract: We consider the problem of minimization of a convex function on a simple set with convex non-smooth inequality constraint and describe first-order methods to solve such problems in different situations: smooth or non-smooth objective function; convex or strongly convex objective and constraint; deterministic or randomized information about the objective and constraint. We hope that it is convenient for a reader to have all the methods for different settings in one place. Described methods are based on Mirror Descent algorithm and switching subgradient scheme. One of our focus is to propose, for the listed different settings, a Mirror Descent with adaptive stepsizes and adaptive stopping rule. This means that neither stepsize nor stopping rule require to know the Lipschitz constant of the objective or constraint. We also construct Mirror Descent for problems with objective function, which is not Lipschitz continuous, e.g. is a quadratic function. Besides that, we address the problem of recovering the solution of the dual problem.