title: Design Analysis for Optimal Calibration of Diffusivity in Reactive Multilayers

abstract: Calibration of the uncertain Arrhenius diffusion parameters for quantifying mixing rates in Zr-Al nanolaminate foils was performed in a Bayesian setting [Vohra et al., 2014]. The parameters were inferred in a low temperature regime characterized by homogeneous ignition and a high temperature regime characterized by self-propagating reactions in the multilayers. In this work, we extend the analysis to find optimal experimental designs that would provide the best data for inference. We employ a rigorous framework that quantifies the expected in- formation gain in an experiment, and find the optimal design conditions using numerical techniques of Monte Carlo, sparse quadrature, and polynomial chaos surrogates. For the low temperature regime, we find the optimal foil heating rate and pulse duration, and confirm through simulation that the optimal design indeed leads to sharper posterior distributions of the diffusion parameters. For the high temperature regime, we demonstrate potential for increase in the expected information gain of the posteriors by increasing sample size and reducing uncertainty in measurements. Moreover, posterior marginals are also produced to verify favorable experimental scenarios for this regime.


title: Convolution identities for Tribonacci numbers with symmetric formulae

abstract: Many kinds of convolution identities have been considered about several numbers, including Bernoulli, Euler, Genocchi, Cauchy, Stirling, and Fibonacci numbers. The well-known basic result about Bernoulli numbers is due to Euler. The convolution identities have been studied. In this paper, we give convolution identities for Tribonacci numbers without binomial coefficients and with binomial coefficients. Convolution identities of Fiboancci numbers or Lucas numbers can be expressed in the form of linear combinations of Fibonacci numbers and Lucas numbers only. Fibonacci numbers and Lucas numbers are in pairs with different values. Covolution identities of Tribonacci numbers can be expressed in the linear combination of several Tribonacci numbers with different values. Symmetric formulas are basic tools for these expressions.


title: Calculation of entanglement in graph states up to five-qubit based on generalized concurrence

abstract: We propose a new classification for the entanglement in graph states based on generalized con- currence. The numerical results indicate that the eight different three-qubit graph states in three categories, 64 four-qubit graph states in five categories and 1024 five-qubit graph states are in ten classes. We also compare this classification with equivalence classes of these graph states under local complementation (LC) operator, and the obtained result suggests that classification by generalized concurrence is not in contradiction with the LC-rule.


title: A boundary theory approach to de Finetti's theorem

abstract: We show that boundary theory for transient Markov chains, as initiated by Doob, can be used to prove de Finetti's classical representation result for exchangeable random sequences. We also include the relevant parts of the theory, with full proofs.


title: Some families of generalized Mathieu-type power series, associated probability distributions and related functional inequalities involving complete monotonicity and log-convexity

abstract: By making use of the familiar Mathieu series and its generalizations, the authors derive a number of new integral representations and present a systematic study of probability density functions and probability distributions associated with some generalizations of the Mathieu series. In particular, the mathematical expectation, variance and the characteristic functions, related to the probability density functions of the considered probability distributions are derived. As a consequence, some interesting inequalities involving complete monotonicity and log-convexity are derived.


title: Diabolical entropy

abstract: Milnor and Thurston's famous paper proved monotonicity of the topological entropy for the real quadratic family. Guckenheimer showed that it is H\"older continuous. We obtain a precise formula for the H\"older exponent at almost every quadratic parameter. Furthermore, the entropy of most parameters is proven to be in a set of Hausdorff dimension smaller than one, while most values of the entropy arise from a set of parameters of dimension smaller than one.


title: Anomalous mass dimension of multi-flavor QCD

abstract: Models of strongly interacting theories with a large mass anomalous dimension ($\gamma_m$) provide an interesting possibility for the dynamical origin of the electroweak symmetry breaking. A laboratory for these models is QCD with many flavors, which may present a non-trivial fixed point associated to a conformal region. Studies based on conformal field theories and on Schwinger-Dyson equations have suggested the existence of bounds on the mass anomalous dimension at the fixed points of these models. In this note we discuss $\gamma_m$ values of multi-flavor QCD exhibiting a non-trivial fixed point and affected by relevant four-fermion interactions.


title: High-dynamic-range water window ptychography

abstract: Ptychographic imaging with soft X-rays, especially in the water window energy range, suffers from limited detector dynamic range that directly influences the maximum spatial resolution achievable. High-dynamic-range data can be obtained by multiple exposures. By this approach we have increased the dynamic range of a ptychographic data set by a factor of 76 and obtained diffraction signal till the corners of the detector. The real space half period resolution was improved from 50 nm for the single exposure data to 18 nm for the high-dynamic-range data.


title: Performance Bounds for Remote Estimation under Energy Harvesting Constraints

abstract: Remote estimation with an energy harvesting sensor with a limited data and energy buffer is considered. The sensor node observes an unknown Gaussian field and communicates its observations to a remote fusion center using the energy it harvested. The fusion center employs minimum mean-square error (MMSE) estimation to reconstruct the unknown field. The distortion minimization problem under the online scheme, where the sensor has access to only the statistical information for the future energy packets is considered. We provide performance bounds on the achievable distortion under a slotted block transmission scheme, where at each transmission time slot, the data and the energy buffer are completely emptied. Our bounds provide insights to the trade-offs between the buffer sizes, the statistical properties of the energy harvesting process and the achievable distortion. In particular, these trade-offs illustrate the insensitivity of the performance to the buffer sizes for signals with low degree of freedom and suggest performance improvements with increasing buffer size for signals with relatively higher degree of freedom. Depending only on the mean, variance and finite support of the energy arrival process, these results provide practical insights for the battery and buffer sizes for deployment in future energy harvesting wireless sensing systems.


title: Mining the Web for Pharmacovigilance: the Case Study of Duloxetine and Venlafaxine

abstract: Adverse reactions caused by drugs following their release into the market are among the leading causes of death in many countries. The rapid growth of electronically available health related information, and the ability to process large volumes of them automatically, using natural language processing (NLP) and machine learning algorithms, have opened new opportunities for pharmacovigilance. Survey found that more than 70% of US Internet users consult the Internet when they require medical information. In recent years, research in this area has addressed for Adverse Drug Reaction (ADR) pharmacovigilance using social media, mainly Twitter and medical forums and websites. This paper will show the information which can be collected from a variety of Internet data sources and search engines, mainly Google Trends and Google Correlate. While considering the case study of two popular Major depressive Disorder (MDD) drugs, Duloxetine and Venlafaxine, we will provide a comparative analysis for their reactions using publicly-available alternative data sources.


title: Stable components of sound fields in the ocean

abstract: A method is proposed for finding the wave field components which are weakly sensitive to the sound speed perturbation in the ocean acoustic waveguides. Such a component is formed by a narrow beam of rays whose spread in vertical direction, up to the observation range, remains less than the vertical scale of perturbation. These rays pass through practically the same inhomogeneities and therefore their phases acquire close increments. If the ray amplitudes vary insignificantly, then (i) the stable component of the monochromatic field in the perturbed and unperturbed waveguide differ by only a constant phase factor, and (ii) in the case of transient wave field the perturbation causes only an additional time delay of the stable component as a whole. It is shown how the stable components can be selected from the total wave field using the field expansions in coherent states or in normal modes. The existence of stable components is demonstrated by numerical simulation of sound field in a deep water waveguide. It turns out, that even though the assumptions (i) and (ii) are not met exactly, the stable components in the perturbed and unperturbed waveguide are quite close.


title: Why do only some Radio Galaxies acquire giant Sizes?

abstract: We study the morphology of the radio emission of giant radio galaxies (GRGs), a rare type of active galaxies, in order to find new clues for the reasons of their large size. Using radio images from two sky surveys at 1.4 GHz we quantified their radio morphology by measuring the geometry (armlength, bending angle, etc.) and flux symmetry for 58 such objects. Preliminary analysis suggests that radio source symmetry does not differ between galaxies and quasars, that there is no evidence for a decrease in linear size with redshift, and only a marginal trend for increasing symmetry with larger size. A merging with data available for other such objects is expected to yield new clues on the possible reasons for their large size. We also searched radio survey images around the positions of 1059 galaxies and quasars for further giant radio sources, and present a preliminary sample of six of these.


title: Real-time Error Control for Surgical Simulation

abstract: Objective: To present the first real-time a posteriori error-driven adaptive finite element approach for real-time simulation and to demonstrate the method on a needle insertion problem. Methods: We use corotational elasticity and a frictional needle/tissue interaction model. The problem is solved using finite elements within SOFA. The refinement strategy relies upon a hexahedron-based finite element method, combined with a posteriori error estimation driven local $h$-refinement, for simulating soft tissue deformation. Results: We control the local and global error level in the mechanical fields (e.g. displacement or stresses) during the simulation. We show the convergence of the algorithm on academic examples, and demonstrate its practical usability on a percutaneous procedure involving needle insertion in a liver. For the latter case, we compare the force displacement curves obtained from the proposed adaptive algorithm with that obtained from a uniform refinement approach. Conclusions: Error control guarantees that a tolerable error level is not exceeded during the simulations. Local mesh refinement accelerates simulations. Significance: Our work provides a first step to discriminate between discretization error and modeling error by providing a robust quantification of discretization error during simulations.


title: Transitions in a magnetized quasi-laminar spherical Couette Flow

abstract: First results of a new spherical Couette experiment are presented. The liquid metal flow in a spherical shell is exposed to a homogeneous axial magnetic field. For a Reynolds number Re=1000, we study the effect of increasing Hartmann number Ha. The resulting flow structures are inspected by ultrasound Doppler velocimetry. With a weak applied magnetic field, we observe an equatorially anti-symmetric jet instability with azimuthal wave number m=3. As the magnetic field strength increases, this instability vanishes. When the field is increased further, an equatorially symmetric return flow instability arises. Our observations are shown to be in good agreement with linear stability analysis and non-linear flow simulations.


title: A Search for double-lobed radio emission from Galactic Stars and Spiral Galaxies

abstract: We present a systematic search for two types of very unusual astronomical objects: Galactic stars and spiral galaxies with double radio lobes, i.e. radio emission on opposite sides of the optical object, suggesting the ejection of jets from them. We designed an algorithm to search for pairs of radio sources straddling objects from two unprecedented samples of 878,031 Galactic stars from the Sloan Digital Sky Survey and 675,874 spiral galaxy candidates drawn from the recent literature. We found three new examples of double-lobed radio stars, while for the spiral galaxies we only rediscovered one known such double source, confirming that the latter objects are extremely rare.


Proceedings of the 1st International Workshop on Robot Learning and Planning (RLP 2016)


title: Remarks on the maximum atom-bond connectivity index of graphs with given parameters

abstract: The atom-bond connectivity (ABC) index is a degree-based molecular structure descriptor that can be used for modelling thermodynamic properties of organic chemical compounds. Motivated by its applicable potential, a series of investigations have been carried out in the past several years. In this note we first consider graphs with given edge-connectivity that attain the maximum ABC index. In particular, we give an affirmative answer to the conjecture about the structure of graphs with edge-connectivity equal to one that maximize the ABC index, which was recently raised by Zhang, Yang, Wang and Zhang~\cite{zywz mabciggp-2016}. In addition, we provide supporting evidence for another conjecture posed by the same authors which concerns graphs that maximize the ABC index among all graphs with chromatic number equal to some fixed $\chi \geq 3$. Specifically, we confirm this conjecture in the case where the order of the graph is divisible by $\chi$.


title: Properties of Healthcare Teaming Networks as a Function of Network Construction Algorithms

abstract: Network models of healthcare systems can be used to examine how providers collaborate, communicate, refer patients to each other. Most healthcare service network models have been constructed from patient claims data, using billing claims to link patients with providers. The data sets can be quite large, making standard methods for network construction computationally challenging and thus requiring the use of alternate construction algorithms. While these alternate methods have seen increasing use in generating healthcare networks, there is little to no literature comparing the differences in the structural properties of the generated networks. To address this issue, we compared the properties of healthcare networks constructed using different algorithms and the 2013 Medicare Part B outpatient claims data. Three different algorithms were compared: binning, sliding frame, and trace-route. Unipartite networks linking either providers or healthcare organizations by shared patients were built using each method. We found that each algorithm produced networks with substantially different topological properties. Provider networks adhered to a power law, and organization networks to a power law with exponential cutoff. Censoring networks to exclude edges with less than 11 shared patients, a common de-identification practice for healthcare network data, markedly reduced edge numbers and greatly altered measures of vertex prominence such as the betweenness centrality. We identified patterns in the distance patients travel between network providers, and most strikingly between providers in the Northeast United States and Florida. We conclude that the choice of network construction algorithm is critical for healthcare network analysis, and discuss the implications for selecting the algorithm best suited to the type of analysis to be performed.


title: Schroedinger electrons interacting with optical gratings: quantum mechanical study of the inverse Smith Purcell effect

abstract: Slow swift electrons with low self-inertia interact differently with matter and light in comparison with their relativistic counterparts: they are easily recoiled, reflected, and also diffracted form optical gratings and nanostructures. As a consequence, they can be also better manipulated into the desired shape. For example, they get bunched quite fast in interaction with acceleration gratings in presence of an external electromagnetic radiation, a phenomenon which can be desirable in development of superradiant coherent light sources. Here, I examine the spatiotemporal behavior of pulsed electron wave packets at low energies interacting with pulsed light and optical gratings, using a quantum-mechanical self-consistent numerical toolbox which is introduced here. It will be shown that electron pulses are accelerated very fast in interaction with the near-field of the grating, demanding that a synchronicity condition is met. To prevent the electrons to be transversely deflected from the grating a symmetric double-grating configuration is necessary. It is found that even in this configuration, diffraction due to the interaction of the electron with the standing-wave light inside the gap between the gratings, is a source of defocusing. Moreover, the longitudinal broadening of the electron pulse directly affects the final shape of the electron wave packet due to the occurrence of multiple electron-photon scatterings. These investigations pave the way towards the design of more efficient electron-driven photon sources and accelerators.


title: The Tayler instability at low magnetic Prandtl numbers: Chiral symmetry breaking and synchronizable helicity oscillations

abstract: The current-driven, kink-type Tayler instability (TI) is a key ingredient of the Tayler-Spruit dynamo model for the generation of stellar magnetic fields, but is also discussed as a mechanism that might hamper the up-scaling of liquid metal batteries. Under some circumstances, the TI involves a helical flow pattern which goes along with some alpha effect. Here we focus on the chiral symmetry breaking and the related impact on the alpha effect that would be needed to close the dynamo loop in the Tayler-Spruit model. For low magnetic Prandtl numbers, we observe intrinsic oscillations of the alpha effect. These oscillations serve then as the basis for a synchronized Tayler-Spruit dynamo model, which could possibly link the periodic tidal forces of planets with the oscillation periods of stellar dynamos.


title: Defect tolerance: fundamental limits and examples

abstract: This paper addresses the problem of adding redundancy to a collection of physical objects so that the overall system is more robust to failures. In contrast to its information counterpart, which can exploit parity to protect multiple information symbols from a single erasure, physical redundancy can only be realized through duplication and substitution of objects. We propose a bipartite graph model for designing defect-tolerant systems in which defective objects are replaced by judiciously connected redundant objects. The fundamental limits of this model are characterized under various asymptotic settings and both asymptotic and finite-size systems that approach these limits are constructed. Among other results, we show that simple modular redundancy is in general suboptimal. As we develop, this combinatorial problem of defect tolerant system design has a natural interpretation as one of graph coloring, and the analysis is significantly different from that traditionally used in information redundancy for error-control codes.


title: Crafting GBD-Net for Object Detection

abstract: The visual cues from multiple support regions of different sizes and resolutions are complementary in classifying a candidate box in object detection. Effective integration of local and contextual visual cues from these regions has become a fundamental problem in object detection. In this paper, we propose a gated bi-directional CNN (GBD-Net) to pass messages among features from different support regions during both feature learning and feature extraction. Such message passing can be implemented through convolution between neighboring support regions in two directions and can be conducted in various layers. Therefore, local and contextual visual patterns can validate the existence of each other by learning their nonlinear relationships and their close interactions are modeled in a more complex way. It is also shown that message passing is not always helpful but dependent on individual samples. Gated functions are therefore needed to control message transmission, whose on-or-offs are controlled by extra visual evidence from the input sample. The effectiveness of GBD-Net is shown through experiments on three object detection datasets, ImageNet, Pascal VOC2007 and Microsoft COCO. This paper also shows the details of our approach in wining the ImageNet object detection challenge of 2016, with source code provided on \url{https://github.com/craftGBD/craftGBD}.


title: First results from a microwave cavity axion search at 24 micro-eV

abstract: We report on the first results from a new microwave cavity search for dark matter axions with masses above $20~\mu\text{eV}$. We exclude axion models with two-photon coupling $g_{a\gamma\gamma} \gtrsim 2\times10^{-14}~\text{GeV}^{-1}$ over the range $23.55~\mu\text{eV} < m_a < 24.0~\mu\text{eV}$. These results represent two important achievements. First, we have reached cosmologically relevant sensitivity an order of magnitude higher in mass than any existing limits. Second, by incorporating a dilution refrigerator and Josephson parametric amplifier, we have demonstrated total noise approaching the standard quantum limit for the first time in an axion search.


title: Variance-based regularization with convex objectives

abstract: We develop an approach to risk minimization and stochastic optimization that provides a convex surrogate for variance, allowing near-optimal and computationally efficient trading between approximation and estimation error. Our approach builds off of techniques for distributionally robust optimization and Owen's empirical likelihood, and we provide a number of finite-sample and asymptotic results characterizing the theoretical performance of the estimator. In particular, we show that our procedure comes with certificates of optimality, achieving (in some scenarios) faster rates of convergence than empirical risk minimization by virtue of automatically balancing bias and variance. We give corroborating empirical evidence showing that in practice, the estimator indeed trades between variance and absolute performance on a training sample, improving out-of-sample (test) performance over standard empirical risk minimization for a number of classification problems.


title: Contraction Principles in $M_s$-metric Spaces

abstract: In this paper, we give an interesting extension of the partial S-metric space which was introduced [4] to the M_s-metric space. Also, we prove the existence and uniqueness of a fixed point for a self mapping on an Ms-metric space under different contraction principles.


title: A Gentle Tutorial of Recurrent Neural Network with Error Backpropagation

abstract: We describe recurrent neural networks (RNNs), which have attracted great attention on sequential tasks, such as handwriting recognition, speech recognition and image to text. However, compared to general feedforward neural networks, RNNs have feedback loops, which makes it a little hard to understand the backpropagation step. Thus, we focus on basics, especially the error backpropagation to compute gradients with respect to model parameters. Further, we go into detail on how error backpropagation algorithm is applied on long short-term memory (LSTM) by unfolding the memory unit.


title: Perspectives on the mathematics of biological patterning and morphogenesis

abstract: A central question in developmental biology is how size and position are determined. The genetic code carries instructions on how to control these properties in order to regulate the pattern and morphology of structures in the developing organism. Transcription and protein translation mechanisms implement these instructions. However, this cannot happen without some manner of sampling of epigenetic information on the current patterns and morphological forms of structures in the organism. Any rigorous description of space- and time-varying patterns and morphological forms reduces to one among various classes of spatio-temporal partial differential equations. Reaction-transport equations represent one such class. Starting from simple Fickian diffusion, the incorporation of reaction, phase segregation and advection terms can represent many of the patterns seen in the animal and plant kingdoms. Morphological form, requiring the development of three-dimensional structure, also can be represented by these equations of mass transport, albeit to a limited degree. The recognition that physical forces play controlling roles in shaping tissues leads to the conclusion that (nonlinear) elasticity governs the development of morphological form. In this setting, inhomogeneous growth drives the elasticity problem. The combination of reaction-transport equations with those of elasto-growth makes accessible a potentially unlimited spectrum of patterning and morphogenetic phenomena in developmental biology. This perspective communication is a survey of the partial differential equations of mathematical physics that have been proposed to govern patterning and morphogenesis in developmental biology. Several numerical examples are included to illustrate these equations and the corresponding physics, with the intention of providing physical insight wherever possible.


title: HD 164492C: a rapidly-rotating, H$\alpha$-bright, magnetic early B star associated with a 12.5d spectroscopic binary

abstract: We employ high resolution spectroscopy and spectropolarimetry to derive the physical properties and magnetic characteristics of the multiple system HD 164492C, located in the young open cluster M20. The spectrum reveals evidence of 3 components: a broad-lined early B star (HD 164492C1), a narrow-lined early B star (HD 164492C2), and a late B star (HD 164492C3). Components C2 and C3 exhibit significant ($>100$ km/s) bulk radial velocity variations with a period of $12.5351(7)$ d that we attribute to eccentric binary motion around a common centre-of-mass. Component C1 exhibits no detectable radial velocity variations. Using constraints derived from modeling the orbit of the C2+C3 binary and from synthesis of the combined spectrum, we determine the approximate physical characteristics of the components. We conclude that a coherent evolutionary solution consistent with the published age of M20 implies a distance to the system of $0.9\pm 0.2$ kpc, corresponding to the smallest published values. We confirm the detection of a strong magnetic field in the combined spectrum. The field is clearly associated with the broad-lined C1 component of the system. Repeated measurement of the longitudinal magnetic field allows the derivation of the rotation period of the magnetic star, $P_{\rm rot}=1.36986(6)$ d. We derive the star's magnetic geometry, finding $i=63\pm 6$ deg, $\beta=33\pm 6$ deg and a dipole polar strength $B_{\rm d}=7.9^{+1.2}_{-1.0}$ kG. Strong emission - varying according to the magnetic period - is detected in the H$\alpha$ profile. This is consistent with the presence of a centrifugal magnetosphere surrounding the rapidly rotating magnetic C1 component.


title: Baryon asymmetry from primordial black holes

abstract: We propose a new scenario of the baryogenesis from primordial black holes (PBH). Assuming presence of microscopic baryon (or lepton) number violation and a CP violating operator such as $\partial_\alpha F(\mathcal{R_{....}} ) J^\alpha$, where $F(\mathcal{R_{....}})$ is a scalar function of the Riemann tensor, time evolution of an evaporating black hole generates baryonic (leptonic) chemical potential at the horizon; consequently PBH enumerates asymmetric Hawking radiation between baryons (leptons) and anti-baryons (leptons). Though the operator is higher dimensional and largely suppressed by a high mass scale $M_*$, we show that sufficient amount of asymmetry can be generated for a wide range of parameters of the PBH mass $M_{\rm PBH}$, its abundance $\Omega_{\rm PBH}$, and the scale $M_*$.


title: Partially Observed Optimal Control for Mean-Field SDEs

abstract: In this paper, we are concerned with a stochastic optimal control problem of mean-field type under partial observation, where the state equation is governed by the controlled nonlinear mean-field stochastic differential equation, moreover the observation noise is allowed to enter into the state equation and the observation coefficients may depend not only on the control process and but also on its probability distribution. Under standard assumptions on the coefficients, by dual analysis and convex variation, we establish the maximum principle for optimal control in a strong sense as well as a weak one, respectively. As an application, a partially observed linear quadratic control problem of mean-field type is studied detailed and the corresponding dual characterization and state feedback presentation of the partially observed optimal control are obtained by the stochastic maximum principles and the classic technique of completing squares.


title: Iterative proportional scaling revisited: a modern optimization perspective

abstract: This paper revisits the classic iterative proportional scaling (IPS) from a modern optimization perspective. In contrast to the criticisms made in the literature, we show that based on a coordinate descent characterization, IPS can be slightly modified to deliver coefficient estimates, and from a majorization-minimization standpoint, IPS can be extended to handle log-affine models with features not necessarily binary-valued or nonnegative. Furthermore, some state-of-the-art optimization techniques such as block-wise computation, randomization and momentum-based acceleration can be employed to provide more scalable IPS algorithms, as well as some regularized variants of IPS for concurrent feature selection.


title: Modeling, implementation and performance analysis of mobility load balancing for LTE downlink data transmission

abstract: We propose in this paper a simulation implementation of self-organizing Networks optimization related to mobility load balancing (MLB) for lte systems using ns-3. the implementation is achieved toward two MLB algorithms dynamically adjusting handover (HO) parameters based on the reference signal received power measurements. such adjustments are done with respect to loads of both an overloaded cell and its cell neighbours having enough available resources enabling to achieve load balancing. numerical investigations through selected key performance indicators of the proposed MLB algorithms when compared with another HO algorithm (already implemented in ns3) based on A3 event highlight the significant MLB gains provided in terms global network throughput, packet loss rate and the number of successful HO without incurring significant overhead.


title: Indirect Gaussian Graph Learning beyond Gaussianity

abstract: This paper studies how to capture dependency graph structures from real data which may not be Gaussian. Starting from marginal loss functions not necessarily derived from probability distributions, we utilize an additive over-parametrization with shrinkage to incorporate variable dependencies into the criterion. An iterative Gaussian graph learning algorithm is proposed with ease in implementation. Statistical analysis shows that the estimators achieve satisfactory accuracy with the error measured in terms of a proper Bregman divergence. Real-life examples in different settings are given to demonstrate the efficacy of the proposed methodology.


title: Solving Marginal MAP Problems with NP Oracles and Parity Constraints

abstract: Arising from many applications at the intersection of decision making and machine learning, Marginal Maximum A Posteriori (Marginal MAP) Problems unify the two main classes of inference, namely maximization (optimization) and marginal inference (counting), and are believed to have higher complexity than both of them. We propose XOR_MMAP, a novel approach to solve the Marginal MAP Problem, which represents the intractable counting subproblem with queries to NP oracles, subject to additional parity constraints. XOR_MMAP provides a constant factor approximation to the Marginal MAP Problem, by encoding it as a single optimization in polynomial size of the original problem. We evaluate our approach in several machine learning and decision making applications, and show that our approach outperforms several state-of-the-art Marginal MAP solvers.


title: An Introduction to Geometric Topology

abstract: This book provides a self-contained introduction to the topology and geometry of surfaces and three-manifolds. The main goal is to describe Thurston's geometrisation of three-manifolds, proved by Perelman in 2002. The book is divided into three parts: the first is devoted to hyperbolic geometry, the second to surfaces, and the third to three-manifolds. It contains complete proofs of Mostow's rigidity, the thick-thin decomposition, Thurston's classification of the diffeomorphisms of surfaces (via Bonahon's geodesic currents), the prime and JSJ decomposition, the topological and geometric classification of Seifert manifolds, and Thurston's hyperbolic Dehn filling Theorem.


title: Superlinearity of geodesic length in 2$D$ critical first-passage percolation

abstract: First-passage percolation is the study of the metric space $(\mathbb{Z}^d,T)$, where $T$ is a random metric defined as the weighted graph metric using random edge-weights $(t_e)_{e\in \mathcal{E}^d}$ assigned to the nearest-neighbor edges $\mathcal{E}^d$ of the $d$-dimensional cubic lattice. We study the so-called critical case in two dimensions, in which $\mathbb{P}(t_e=0)=p_c$, where $p_c$ is the threshold for two-dimensional bond percolation. In contrast to the standard case $(<p_c)$, the distance $T(0,x)$ in the critical case grows sub linearly in $x$ and geodesics are expected to have Euclidean length which is superlinear. We show a strong version of this super linearity, namely that there is $s>1$ such that with probability at least $1-e^{-\|x\|_1^c}$, the minimal length geodesic from $0$ to $x$ has at least $\|x\|_1^s$ number of edges. Our proofs combine recent ideas to bound $T$ for general critical distributions, and modifications of techniques of Aizenman-Burchard to estimate the Hausdorff dimension of random curves.


title: Integrable nonlocal nonlinear equations

abstract: A nonlocal nonlinear Schr\"odinger (NLS) equation was recently found by the authors and shown to be an integrable infinite dimensional Hamiltonian equation. Unlike the classical (local) case, here the nonlinearly induced "potential" is $PT$ symmetric thus the nonlocal NLS equation is also $PT$ symmetric. In this paper, new {\it reverse space-time} and {\it reverse time} nonlocal nonlinear integrable equations are introduced. They arise from remarkably simple symmetry reductions of general AKNS scattering problems where the nonlocality appears in both space and time or time alone. They are integrable infinite dimensional Hamiltonian dynamical systems. These include the reverse space-time, and in some cases reverse time, nonlocal nonlinear Schr\"odinger, modified Korteweg-deVries (mKdV), sine-Gordon, $(1+1)$ and $(2+1)$ dimensional three-wave interaction, derivative NLS, "loop soliton", Davey-Stewartson (DS), partially $PT$ symmetric DS and partially reverse space-time DS equations. Linear Lax pairs, an infinite number of conservation laws, inverse scattering transforms are discussed and one soliton solutions are found. Integrable reverse space-time and reverse time nonlocal discrete nonlinear Schr\"odinger type equations are also introduced along with few conserved quantities. Finally, nonlocal Painlev\'e type equations are derived from the reverse space-time and reverse time nonlocal NLS equations.


title: Constraint Control of Nonholonomic Mechanical Systems

abstract: We derive an optimal control formulation for a nonholonomic mechanical system using the nonholonomic constraint itself as the control. We focus on Suslov's problem, which is defined as the motion of a rigid body with a vanishing projection of the body frame angular velocity on a given direction $\boldsymbol{\xi}$. We derive the optimal control formulation, first for an arbitrary group, and then in the classical realization of Suslov's problem for the rotation group $SO(3)$. We show that it is possible to control the system using the constraint $\boldsymbol{\xi}(t)$ and demonstrate numerical examples in which the system tracks quite complex trajectories such as a spiral.


title: The Discrete Adjoint Method for Exponential Integration

abstract: The implementation of the discrete adjoint method for exponential time differencing (ETD) schemes is considered. This is important for parameter estimation problems that are constrained by stiff time-dependent PDEs when the discretized PDE system is solved using an exponential integrator. We also discuss the closely related topic of computing the action of the sensitivity matrix on a vector, which is required when performing a sensitivity analysis. The PDE system is assumed to be semi-linear and can be the result of a linearization of a nonlinear PDE, leading to exponential Rosenbrock-type methods. We discuss the computation of the derivatives of the $\varphi$-functions that are used by ETD schemes and find that the derivatives strongly depend on the way the $\varphi$-functions are evaluated numerically. A general adjoint exponential integration method, required when computing the gradients, is developed and its implementation is illustrated by applying it to the Krogstad scheme. The applicability of the methods developed here to pattern formation problems is demonstrated using the Swift-Hohenberg model.


title: Unveiling the link between fractional Schr\"odinger equation and light propagation in honeycomb lattice

abstract: We establish a link between the fractional Schr\"odinger equation (FSE) and light propagation in the honeycomb lattice (HCL) - the Dirac-Weyl equation (DWE). The fractional Laplacian in FSE causes a modulation of the dispersion relation of the system, which in the limiting case becomes linear. In the HCL, the dispersion relation is already linear around the Dirac point, suggesting a possible connection with the FSE. Here, we demonstrate this connection by describing light propagation in both FSE and HCL, using DWE. Thus, we propagate Gaussian beams according to FSE, HCL around the Dirac point, and DWE, to discover very similar behavior - the conical diffraction. However, if an additional potential is brought into the system, the link between FSE and HCL is broken, because the added potential serves as a perturbation, which breaks the translational periodicity of HCL and destroys Dirac cones in the dispersion relation.


title: Physical model for turbulent friction on rough surfaces

abstract: We present a physical model for turbulent friction on rough surfaces with regularly distributed roughness elements. Wall shear stresses are expressed as functions of physical quantities. Surfaces with varying roughness densities and roughness elements with different aspect ratios are considered. We propose a straight forward method based on the conservation of momentum to deduce the drag on elements by expressing it as functions of the maximum drag and drag reductions ratios, as the drag on individual elements decreases as packing density increases. A drag reduction effect of momentum redistribution is proposed and the mutual sheltering effect is studied. These two drag reduction mechanisms for individual elements are significant for sparse and dense surfaces, respectively. Reduction ratios for redistribution effect and mutual sheltering effect are deduced, for the two different types of rough surfaces. The shear stress on elements and the total wall shear stress are obtained as the result of the drag analysis. The estimated wall shear stresses of the proposed model are consistent with classical experimental measurements.


title: The Jacobian ideal of a commutative ring and annihilators of cohomology

abstract: It is proved that for a ring $R$ that is either an affine algebra over a field, or an equicharacteristic complete local ring, some power of the Jacobian ideal of $R$ annihilates $\mathrm{Ext}^{d+1}_{R}(-,-)$, where $d$ is the Krull dimension of $R$. Sufficient conditions are identified under which the Jacobian ideal itself annihilates these Ext-modules, and examples are provided that show that this is not always the case. A crucial new idea is to consider a derived version of the Noether different of an algebra.


title: Relative Morita equivalence of Cuntz--Krieger algebras and flow equivalence of topological Markov shifts

abstract: In this paper, we will introduce notions of relative version of imprimitivity bimodules and relative version of strong Morita equivalence for pairs of $C^*$-algebras $(\mathcal{A}, \mathcal{D})$ such that $\mathcal{D}$ is a $C^*$-subalgebra of $\mathcal{A}$ with certain conditions. We will then prove that two pairs $(\mathcal{A}_1, \mathcal{D}_1)$ and $(\mathcal{A}_2, \mathcal{D}_2)$ are relatively Morita equivalent if and only if their relative stabilizations are isomorphic. In particularly, for two pairs $(\mathcal{O}_A, \mathcal{D}_A)$ and $(\mathcal{O}_B, \mathcal{D}_B)$ of Cuntz--Krieger algebras with their canonical masas, they are relatively Morita equivalent if and only if their underlying two-sided topological Markov shifts $(\bar{X}_A,\bar{\sigma}_A)$ and $(\bar{X}_B,\bar{\sigma}_B)$ are flow equivalent. We also introduce a relative version of the Picard group ${\operatorname{Pic}}(\mathcal{A}, \mathcal{D})$ for the pair $(\mathcal{A}, \mathcal{D})$ of $C^*$-algebras and study them for the Cuntz--Krieger pair $(\mathcal{O}_A, \mathcal{D}_A)$.


title: Stellar Velocity Dispersion in Mergers: The Effects of Dust and Star Formation

abstract: We investigate the effects of stellar evolution and dust on measurements of stellar velocity dispersion in mergers of disk galaxies. $N$-body simulations and radiative transfer analysis software are used to obtain mass-weighted and flux-weighted measurements of stellar velocity dispersion. We find that the distribution of dust with respect to the distribution of young stars in such systems is more important than the total degree of attenuation. The presence of dust typically causes flux-weighted measurements of stellar velocity dispersion to be elevated with respect to mass-weighted measurements because dust preferentially obscures young stars, which tend to be dynamically cooler than older stellar populations in such systems. In exceptional situations, in which young stars are not preferentially obscured by dust, flux-weighted velocity dispersion measurements tend to be negatively offset with respect to mass-weighted measurements because the dynamically cool young stellar populations are more luminous, per unit mass, than older stellar populations. Our findings provide a context for comparing observationally-obtained measurements of velocity dispersion with measurements of velocity dispersion obtained from galaxy merger simulations.


title: Finite Rank Isopairs

abstract: An algebraic isopair is a commuting pair of pure isometries that is annihilated by a polynomial defining a distinguished variety $\mathcal{V}$. The notion of the rank of a pure algebraic isopair with finite bimultiplicity is introduced. For $\mathcal{V} $, a union of $s$ irreducible varieties $\mathcal{V}_j$, the rank is a $s$-tuple $\alpha=(\alpha_1,...,\alpha_s)$ of natural numbers. A pure algebraic isopair of finite bimultiplicity with rank $\alpha$ is described as a restriction of a $\max\{\alpha_1,...,\alpha_s\}$-cyclic pure algebraic isopair to a finite codimensional invariant subspace. The restriction of a pure algebraic isopair of finite bimultiplicity with rank $\alpha$ to a finite codimensional invariant subspace is at least $\max\{\alpha_1,...,\alpha_s\}$-cyclic and there is a $\max\{\alpha_1,...,\alpha_s\}$-cyclic finite codimensional invariant subspace.


title: Existence of a highest wave in a fully dispersive two-way shallow water model

abstract: We consider the existence of periodic traveling waves in a bidirectional Whitham equation, combining the full two-way dispersion relation from the incompressible Euler equations with a canonical shallow water nonlinearity. Of particular interest is the existence of a highest, cusped, traveling wave solution, which we obtain as a limiting case at the end of the main bifurcation branch of $2\pi$-periodic traveling wave solutions continuing from the zero state. Unlike the unidirectional Whitham equation, containing only one branch of the full Euler dispersion relation, where such a highest wave behaves like $|x|^{1/2}$ near its crest, the cusped waves obtained here behave like $|x\log|x||$. Although the linear operator involved in this equation can be easily represented in terms of an integral operator, it maps continuous functions out of the H\"older and Lipschitz scales of function spaces by introducing logarithmic singularities. Since the nonlinearity is also of higher order than that of the unidirectional Whitham equation, several parts of our proofs and results deviate from those of the corresponding unidirectional equation, with the analysis of the logarithmic singularity being the most subtle component. This paper is part of a longer research programme for understanding the interplay between nonlinearities and dispersion in the formation of large-amplitude waves and their singularities.


title: Relating Chromophoric and Structural Disorder in Conjugated Polymers

abstract: The optoelectronic properties of amorphous conjugated polymers are sensitive to conformational disorder and spectroscopy provides the means for structural characterization of the fragments of the chain which interact with light - "chromophores". A faithful interpretation of spectroscopic conformational signatures, however, presents a key challenge. We investigate the relationship between the ground state optical gaps, the properties of the excited states, and the structural features of chromophores of a single molecule poly(3-hexyl)-thiophene (P3HT), using quantum-classical atomistic simulations. Our results demonstrate that chromophoric disorder reflects an interplay between excited state de-localization and electron-hole polarization, and is controlled by torsional disorder that is specifically associated with the presence of side chains. Within this conceptual framework, we predict and explain a counter-intuitive spectral signature of P3HT: a red-shifted absorption, despite shortening of chromophores, with increasing temperature.


title: Exact time evolution of the asymmetric Hubbard dimer

abstract: We examine the time evolution of an asymmetric Hubbard dimer, which has a different on-site interaction on the two sites. The Hamiltonian has a time-dependent hopping term, which can be employed to describe an electric field (which creates a Hamiltonian with complex matrix elements), or it can describe a modulation of the lattice (which has real matrix elements). By examining the symmetries under spin and pseudospin, we show that the former case involves at most a 3 x 3 block---it can be mapped onto the time evolution of a time-independent Hamiltonian, so the dynamics can be evaluated analytically and exactly (by solving a nontrivial cubic equation). We also show that the latter case reduces to at most 2 x 2 blocks, and hence the time evolution for a single Trotter step can be determined exactly, but the time evolution generically requires a Trotter product.


title: Doing Moore with Less -- Leapfrogging Moore's Law with Inexactness for Supercomputing

abstract: Energy and power consumption are major limitations to continued scaling of computing systems. Inexactness, where the quality of the solution can be traded for energy savings, has been proposed as an approach to overcoming those limitations. In the past, however, inexactness necessitated the need for highly customized or specialized hardware. The current evolution of commercial off-the-shelf(COTS) processors facilitates the use of lower-precision arithmetic in ways that reduce energy consumption. We study these new opportunities in this paper, using the example of an inexact Newton algorithm for solving nonlinear equations. Moreover, we have begun developing a set of techniques we call reinvestment that, paradoxically, use reduced precision to improve the quality of the computed result: They do so by reinvesting the energy saved by reduced precision.


title: The External Calibrator for Hydrogen Observatories

abstract: Multiple instruments are pursuing constraints on dark energy, observing reionization and opening a window on the dark ages through the detection and characterization of the 21cm hydrogen line across the redshift spectrum, from nearby to z=25. These instruments, including CHIME in the sub-meter and HERA in the meter bands, are wide-field arrays with multiple-degree beams, typically operating in transit mode. Accurate knowledge of their primary beams is critical for separation of bright foregrounds from the desired cosmological signals, but difficult to achieve through astronomical observations alone. Previous beam calibration work has focused on model verification and does not address the need of 21cm experiments for routine beam mapping, to the horizon, of the as-built array. We describe the design and methodology of a drone-mounted calibrator, the External Calibrator for Hydrogen Observatories (ECHO), that aims to address this need. We report on a first set of trials to calibrate low-frequency dipoles and compare ECHO measurements to an established beam-mapping system based on transmissions from the Orbcomm satellite constellation. We create beam maps of two dipoles at a 9-degree resolution and find sample noise ranging from 1 to 2%. Assuming this sample noise represents the error in the measurement, the higher end of this range is roughly twice the desired requirement. The overall performance of ECHO suggests that the desired precision and angular coverage is achievable in practice with modest improvements. We identify the main sources of systematic error and uncertainty in our measurements and describe the steps needed to overcome them.


title: Research and Education in Computational Science and Engineering

abstract: Over the past two decades the field of computational science and engineering (CSE) has penetrated both basic and applied research in academia, industry, and laboratories to advance discovery, optimize systems, support decision-makers, and educate the scientific and engineering workforce. Informed by centuries of theory and experiment, CSE performs computational experiments to answer questions that neither theory nor experiment alone is equipped to answer. CSE provides scientists and engineers of all persuasions with algorithmic inventions and software systems that transcend disciplines and scales. Carried on a wave of digital technology, CSE brings the power of parallelism to bear on troves of data. Mathematics-based advanced computing has become a prevalent means of discovery and innovation in essentially all areas of science, engineering, technology, and society; and the CSE community is at the core of this transformation. However, a combination of disruptive developments---including the architectural complexity of extreme-scale computing, the data revolution that engulfs the planet, and the specialization required to follow the applications to new frontiers---is redefining the scope and reach of the CSE endeavor. This report describes the rapid expansion of CSE and the challenges to sustaining its bold advances. The report also presents strategies and directions for CSE research and education for the next decade.


title: Learning Human-Robot Handovers Through $\pi$-STAM: Policy Improvement With Spatio-Temporal Affordance Maps

abstract: Human-robot handovers are characterized by high uncertainty and poor structure of the problem that make them difficult tasks. While machine learning methods have shown promising results, their application to problems with large state dimensionality, such as in the case of humanoid robots, is still limited. Additionally, by using these methods and during the interaction with the human operator, no guarantees can be obtained on the correct interpretation of spatial constraints (e.g., from social rules). In this paper, we present Policy Improvement with Spatio-Temporal Affordance Maps -- $\pi$-STAM, a novel iterative algorithm to learn spatial affordances and generate robot behaviors. Our goal consists in generating a policy that adapts to the unknown action semantics by using affordances. In this way, while learning to perform a human-robot handover task, we can (1) efficiently generate good policies with few training episodes, and (2) easily encode action semantics and, if available, enforce prior knowledge in it. We experimentally validate our approach both in simulation and on a real NAO robot whose task consists in taking an object from the hands of a human. The obtained results show that our algorithm obtains a good policy while reducing the computational load and time duration of the learning process.


title: Visual Closed-Loop Control for Pouring Liquids

abstract: Pouring a specific amount of liquid is a challenging task. In this paper we develop methods for robots to use visual feedback to perform closed-loop control for pouring liquids. We propose both a model-based and a model-free method utilizing deep learning for estimating the volume of liquid in a container. Our results show that the model-free method is better able to estimate the volume. We combine this with a simple PID controller to pour specific amounts of liquid, and show that the robot is able to achieve an average 38ml deviation from the target amount. To our knowledge, this is the first use of raw visual feedback to pour liquids in robotics.


title: On the zeros of random harmonic polynomials: the naive model

abstract: A complex harmonic polynomial is the sum of a complex polynomial and a conjugated complex polynomial, of degrees $n$ and $m$ respectively. Li and Wei (2009) presented a formula for the expected number of zeros of a random harmonic polynomial in $\mathbb{C}$. In this paper we prove that if $m$ is a fixed number, this expectation is asymptotically $n$ as $n\rightarrow\infty$, and if $m=n$ we find a lower and upper bound of order $n\log n$ for this expectation for sufficiently large $n$.


title: Generalized Boltzmann Distribution for Systems out of Equilibrium

abstract: The Boltzmann distribution predicts the collective behavior of systems at thermodynamic equilibrium as a function of their constituent parts. Yet most systems in nature are not at equilibrium, and a unified theory of their behavior does not currently exist. Here, I show that the Boltzmann distribution is a special case of a general distribution that governs all stochastic systems, even if far from equilibrium. The generalized Boltzmann distribution is explained as an analog of the voltage equation in electronics, where resistors, batteries, node voltages, and path currents correspond to equilibrium rate constants, driven rate constants, probabilities, and probability flows, respectively. The general distribution recapitulates known properties of weakly driven systems and enables new closed-form solutions for strongly-driven systems. These solutions provide insight into fundamental limits on system performance, and experimental data show that living systems can operate at those limits. The formal mapping between non-equilibrium systems and electronic circuits may provide a unified framework to simplify, understand, and ultimately control the behavior of complex non-equilibrium systems.


title: Theoretical description of pump/probe experiments in electron mediated charge-density-wave insulators

abstract: In this review, we develop the formalism employed to describe charge-density-wave insulators in pump/probe experiments using ultra short driving pulses of light. The theory emphasizes exact results in the simplest model for a charge-density wave insulator (given by a noninteracting systems with two bands and a gap) and by employing nonequilibrium dynamical mean-field theory to solve the Falicov-Kimball model in its ordered phase. We show both how to develop the formalism and how the solutions behave. Care is taken to describe the details behind these calculations and to show how to verify their accuracy via sum-rule constraints.


title: Tunneling Dynamics Between Atomic Bright Solitons

abstract: We investigate tunneling behavior between two bright solitons in a Bose-Einstein condensate with attractive contact interactions between atoms. The explicit tunneling properties including tunneling particles and oscillation period are described analytically, which indicates that the periodic tunneling form is a nonlinear Josephson type oscillation. The results suggest that the breathing behavior of solitons comes from the tunneling mechanism in an effective double-well potential, which is quite different from the modulational instability mechanism for Akhmediev breather and K-M breather. Furthermore, we obtain a phase diagram for two soliton interaction which admits tunneling property, particle-like property, interference property, and a resonant interaction case. The explicit conditions for them are clarified based on the defined critical distance $d_c$ and spatial interference period $D$.


title: A note on homological properties of Nakayama algebras

abstract: Using the resolution quiver for a connected Nakayama algebra, a fast algorithm is given to decide whether its global dimension is finite or not and whether it is Gorenstein or not. The latter strengthens a result of Ringel.


title: Learning Spatial-Semantic Context with Fully Convolutional Recurrent Network for Online Handwritten Chinese Text Recognition

abstract: Online handwritten Chinese text recognition (OHCTR) is a challenging problem as it involves a large-scale character set, ambiguous segmentation, and variable-length input sequences. In this paper, we exploit the outstanding capability of path signature to translate online pen-tip trajectories into informative signature feature maps using a sliding window-based method, successfully capturing the analytic and geometric properties of pen strokes with strong local invariance and robustness. A multi-spatial-context fully convolutional recurrent network (MCFCRN) is proposed to exploit the multiple spatial contexts from the signature feature maps and generate a prediction sequence while completely avoiding the difficult segmentation problem. Furthermore, an implicit language model is developed to make predictions based on semantic context within a predicting feature sequence, providing a new perspective for incorporating lexicon constraints and prior knowledge about a certain language in the recognition procedure. Experiments on two standard benchmarks, Dataset-CASIA and Dataset-ICDAR, yielded outstanding results, with correct rates of 97.10% and 97.15%, respectively, which are significantly better than the best result reported thus far in the literature.


title: Time-Average Optimization with Non-Convex Decision Set and Its Convergence

abstract: This paper considers time-average optimization, where a decision vector is chosen every time step within a (possibly non-convex) set, and the goal is to minimize a convex function of the time averages subject to convex constraints on these averages. Such problems have applications in networking, multi-agent systems, and operations research, where decisions are constrained to a discrete set and the decision average can represent average bit rates or average agent actions. This time-average optimization extends traditional convex formulations to allow a non-convex decision set. This class of problems can be solved by Lyapunov optimization. A simple drift-based algorithm, related to a classical dual subgradient algorithm, converges to an $\epsilon$-optimal solution within $O(1/\epsilon^2)$ time steps. Further, the algorithm is shown to have a transient phase and a steady state phase which can be exploited to improve convergence rates to $O(1/\epsilon)$ and $O(1/{\epsilon^{1.5}})$ when vectors of Lagrange multipliers satisfy locally-polyhedral and locally-smooth assumptions respectively. Practically, this improved convergence suggests that decisions should be implemented after the transient period.


title: Collider Phenomenology of $e^{-}e^{-}\to W^{-}W^{-}$

abstract: Whether in nature exists Majorana neutrino is one of the most fundamental questions in physics. It is directly related to the violation of accidental $U(1)$-lepton number symmetry(LNV). Enormous efforts had been put into such test and among them, one conventional experiment is the neutrinoless double-beta decay ($0\nu\beta\beta$). On the other hand, there have been proposals of future electron-positron colliders such as International Linear Collider (ILC) for precision measurements of Higgs boson or new discovered beyond SM particles and it can potentially be converted into an electron-electron collider. Such collider enables a new probe of Majorana neutrino via the inverse $0\nu\beta\beta$ process($e^{-}e^{-}\to W^- W^-$) as an alternative and complementary test to the conventional neutrinoless double-beta decay experiments. In this paper, we illustrate the searching strategies of $e^{-}e^{-}\to W^- W^-$ at future electron colliders, and find that the channels of pure hadronic, semi-leptonic with muon and pure leptonic with muons have the most discovery potential.


title: Skeletal Geometric Complexes and Their Symmetries

abstract: Skeletal polyhedra and polygonal complexes are finite or infinite periodic structures in 3-space with interesting geometric, combinatorial, and algebraic properties. These structures can be viewed as finite or infinite periodic graphs (nets) equipped with a polyhedral superstructure imposed by the faces, allowed to be skew, zigzag, or helical. The article describes skeletal structures with maximal symmetry.


title: Reconstruction of signals from their autocorrelation and cross-correlation vectors, with applications to phase retrieval and blind channel estimation

abstract: We consider the problem of reconstructing two signals from the autocorrelation and cross-correlation measurements. This inverse problem is a fundamental one in signal processing, and arises in many applications, including phase retrieval and blind channel estimation. In a typical phase retrieval setup, only the autocorrelation measurements are obtainable. We show that, when the measurements are obtained using three simple "masks", phase retrieval reduces to the aforementioned reconstruction problem. The classic solution to this problem is based on finding common factors between the $z$-transforms of the autocorrelation and cross-correlation vectors. This solution has enjoyed limited practical success, mainly due to the fact that it is not sufficiently stable in the noisy setting. In this work, inspired by the success of convex programming in provably and stably solving various quadratic constrained problems, we develop a semidefinite programming-based algorithm and provide theoretical guarantees. In particular, we show that almost all signals can be uniquely recovered by this algorithm (up to a global phase). Comparative numerical studies demonstrate that the proposed method significantly outperforms the classic method in the noisy setting.


title: Tilting modules of affine quasi-hereditary algebras

abstract: We discuss tilting modules of affine quasi-hereditary algebras. We present an existence theorem of indecomposable tilting modules when the algebra has a large center and use it to deduce a criterion for an exact functor between two affine highest weight categories to give an equivalence. As an application, we prove that the Arakawa-Suzuki functor [Arakawa-Suzuki, J. of Alg. 209 (1998)] gives a fully faithful embedding of a block of the deformed BGG category of $\mathfrak{gl}_{m}$ into the module category of a suitable completion of degenerate affine Hecke algebra of $GL_{n}$.


title: Hamilton-Souplet-Zhang type gradient estimates for porous medium type equations on Riemannian manifolds

abstract: In this paper, by employ the cutoff function and the maximum principle, some Hamilton-Souplet-Zhang type gradient estimates for porous medium type equation are deduced. As a special case, an Hamilton-Souplet-Zhang type gradient estimates of the heat equation is derived which is different from the result of Souplet-Zhang. Furthermore, our results generalize those of Zhu. As application, some Livillous theorems for ancient solution are derived.


title: Singularity-free Numerical Scheme for the Stationary Wigner Equation

abstract: For the stationary Wigner equation with inflow boundary conditions, its numerical convergence with respect to the velocity mesh size are deteriorated due to the singularity at velocity zero. In this paper, using the fact that the solution of the stationary Wigner equation is subject to an algebraic constraint, we prove that the Wigner equation can be written into a form with a bounded operator $\mathcal{B}[V]$, which is equivalent to the operator $\mathcal{A}[V]=\Theta[V]/v$ in the original Wigner equation under some conditions. Then the discrete operators discretizing $\mathcal{B}[V]$ are proved to be uniformly bounded with respect to the mesh size. Based on the therectical findings, a signularity-free numerical method is proposed. Numerical reuslts are proivded to show our improved numerical scheme performs much better in numerical convergence than the original scheme based on discretizing $\mathcal{A}[V]$.


title: Beam halo study on ATF damping ring

abstract: Halo distribution is a key topic for background study. This paper has developed an analytical method to give an estimation of ATF beam halo distribution. The equilibrium particle distribution of the beam tail in the ATF damping ring is calculated analytically with different emittance and different vacuum degree. The analytical results agree the measurements very well. This is a general method which can be applied to any electron rings.


title: Constant Matrix Element Approximation to Time-Resolved Angle-Resolved Photoemission Spectroscopy

abstract: We describe different approximations associated with employing a constant matrix element for the coupling of light to multiband electrons in the context of time-resolved angle-resolved photoemission spectroscopy (TR-ARPES). In particular, we demonstrate that the constant matrix approximation only holds for one of the coordinate representations, and changing to other bases, requires including nonconstant corrections to the matrix element. We also discuss some simplifying approximations, where a constant matrix element is employed in multiple bases, and the consequences of this further approximation (especially with respect to the TR-ARPES signal no longer being nonnegative). We also discuss issues related to gauge invariance of the final spectra.


title: Eccentricity and Spin-Orbit Misalignment in Short-Period Stellar Binaries as a Signpost of Hidden Tertiary Companions

abstract: Eclipsing binaries are observed to have a range of eccentricities and spin-orbit misalignments (stellar obliquities). Whether such properties are primordial, or arise from post-formation dynamical interactions remains uncertain. This paper considers the scenario in which the binary is the inner component of a hierarchical triple stellar system, and derives the requirements that the tertiary companion must satisfy in order to raise the eccentricity and obliquity of the inner binary. Through numerical integrations of the secular octupole-order equations of motion of stellar triples, coupled with the spin precession of the oblate primary star due to the torque from the secondary, we obtain a simple, robust condition for producing spin-orbit misalignment in the inner binary: In order to excite appreciable obliquity, the precession rate of the stellar spin axis must be smaller than the orbital precession rate due to the tertiary companion. This yields quantitative requirements on the mass and orbit of the tertiary. We also present new analytic expressions for the maximum eccentricity and range of inclinations allowing eccentricity excitation (Lidov-Kozai window), for stellar triples with arbitrary masses and including the non-Keplerian potentials introduced by general relativity, stellar tides and rotational bulges. The results of this paper can be used to place constraints on unobserved tertiary companions in binaries that exhibit high eccentricity and/or spin-orbit misalignment, and will be helpful in guiding efforts to detect external companions around stellar binaries. As an application, we consider the eclipsing binary DI Herculis, and identify the requirements that a tertiary companion must satisfy to produce the observed spin-orbit misalignment.


title: Learning Deep Generative Spatial Models for Mobile Robots

abstract: We propose a new probabilistic framework that allows mobile robots to autonomously learn deep, generative models of their environments that span multiple levels of abstraction. Unlike traditional approaches that combine engineered models for low-level features, geometry, and semantics, our approach leverages recent advances in Sum-Product Networks (SPNs) and deep learning to learn a single, universal model of the robot's spatial environment. Our model is fully probabilistic and generative, and represents a joint distribution over spatial information ranging from low-level geometry to semantic interpretations. Once learned, it is capable of solving a wide range of tasks: from semantic classification of places, uncertainty estimation, and novelty detection, to generation of place appearances based on semantic information and prediction of missing data in partial observations. Experiments on laser-range data from a mobile robot show that the proposed universal model obtains performance superior to state-of-the-art models fine-tuned to one specific task, such as Generative Adversarial Networks (GANs) or SVMs.


title: Generating Isolated Terawatt-Attosecond X-ray Pulses via a Chirped Laser Enhanced High-Gain Free-electron Laser

abstract: A feasible method is proposed to generate isolated attosecond terawatt x-ray radiation pulses in high-gain free-electron lasers. In the proposed scheme, a frequency chirped laser pulse is employed to generate a gradually-varied spacing current enhancement of the electron beam and a series of spatiotemporal shifters are applied between the undulator sections to amplify a chosen ultra-short radiation pulse from self-amplified spontaneous emission. Three-dimensional start-to-end simulations have been carried out and the calculation results demonstrated that 0.15 nm x-ray pulses with peak power over 1TW and duration of several tens of attoseconds could be achieved by using the proposed technique.


title: Three Lectures of Flavor and CP violation within and Beyond the Standard Model

abstract: These notes are based on the lectures I gave at the 2015 European School of High-Energy Physics (ESHEP2015). I discuss 1) flavor physics within the Standard Model, 2) effective field theories and Minimal Flavor Violation, 3) flavor physics in theories beyond the Standard Model and "high energy" flavor transitions of the top quark and of the Higgs boson. As a bi-product, I present the most updated constraints from the measurements of B_s -> mu^+mu^-, as well as I discuss the most recent development in the LHC searches for top flavor changing couplings.


title: Direct Meissner Effect Observation of Superconductivity in Compressed H2S

abstract: Recently, an extremely high superconducting temperature (Tc) of ~200 K has been reported in the sulfur hydride system above 100 GPa. This result is supported by theoretical predictions and verified experimentally. The crystal structure of the superconducting phase was also identified experimentally, confirming the theoretically predicted structure as well as a decomposition mechanism from H2S to H3S+S. Even though nuclear resonant scattering has been successfully used to provide magnetic evidence for a superconducting state, a direct measurement of the important Meissner effect is still lacking. Here we report in situ alternating-current magnetic susceptibility measurements on compressed H2S under high pressures. It is shown that superconductivity suddenly appears at 117 GPa and that Tc reaches 183 K at 149 GPa before decreasing monotonically with a further increase in pressure. This evolution agrees with both theoretical calculations and earlier experimental measurements. The idea of conventional high temperature superconductivity in hydrogen-dominant compounds has thus been realized in the sulfur hydride system under hydrostatic pressure, opening further exciting perspectives for possibly realizing room temperature superconductivity in hydrogen-based compounds.


title: Vortex oscillations around a hemisphere-cylinder body with a high fineness ratio

abstract: The vortex unsteadiness around a hemisphere-cylinder body at AOAs of 10 to 80 deg was studied using Large Eddy Simulation (LES) and Dynamic Mode Decomposition (DMD). The Reynolds number (Re) based on the cylinder diameter of the body is 22000. The results show that vortex oscillations exist over the forebody at the whole range of AOAs. The oscillation is characterized by alternate oscillations of a forebody leeward vortex pair up and down and in-phase swings from side to side. The vortex shedding can be found at the afterbody as AOAs more than 20o, and the shedding region moves forwards gradually with AOAs increasing, and accordingly the region of vortex oscillations contracts and eventually only exists near the nose as AOAs sufficiently high. The vortex oscillation and shedding all induce fluctuating side forces along the body, but the ones from vortex oscillations are larger. The frequencies of vortex oscillations are similar to the ones of vortex shedding at the AOAs of 10o-40o with St=0.085-0.12, in which the flow fields over the afterbody are dominated by vortex shedding and forebody-vortex wakes together; while at AOAs of 50o-80o, the frequencies along the body are apparently divided into two regimes in which the vortex oscillations over the forebody have the frequencies of 0.053-0.064, and the vortex shedding over the afterbody has the frequencies of 0.16-0.2. The DMD analysis shows that except the mean flows, the most energetic modes correspond to the vortex oscillation at the forebody and to the vortex shedding at the afterbody respectively, and the frequencies from DMD are identical to the ones of the side forces obtained by fast Fourier transform. The time-averaged flow fields for the vortex oscillations are symmetric, and no apparent asymmetry exists. Therefore, the vortex pair over the forebody oscillates around a symmetric mean flow field.


title: Planetary model of sunspot emergence: A spectral and autocorrelation analysis

abstract: This paper is concerned with intermediate range periodicity in the sunspot area spectrum. An empirical model of sunspot area emergence based on Mercury planet conjunctions was developed and the spectra of the model variation and the sunspot area variation compared. By including solar cycle amplitude modulation and the effect of solar magnetic field reversal the model was able to predict fine detail in the sunspot area spectrum. As Mercury planet conjunctions occur predictably it was possible to compare the time variation of band limited components of sunspot area with the corresponding component variations in the model. When the model component variation was stable corresponding components of sunspot area lagged the model variation by a few tens of days. When a 180 degree phase change occurred in the model variation the corresponding component of sunspot area followed the change over an interval of a few hundred days, first by decreasing to zero and then emerging in phase with the model variation. Where periodicities in sunspot area did not match periods of Fourier components of the model the autocorrelation function of the model variation provided an indirect numerical link to these periodicities. Included in these indirect periodicities were components at periods of about 138, 235, 355, 472, 530, 605, 833, 1190 and 2020 days. Direct periodicities in the model, those linked directly to Fourier components of sunspot area, included components at about 88, 116, 176, 200, 292 and 405 days. The model was developed without allowing for a variable threshold for sunspot emergence. However, it was found that with a variable threshold above which the model variation is effective in triggering sunspots, weak Fourier components at the indirect periods were generated. The model with threshold then encompasses all of the significant intermediate range periodicities in sunspot area.


title: Enabling Medical Translation for Low-Resource Languages

abstract: We present research towards bridging the language gap between migrant workers in Qatar and medical staff. In particular, we present the first steps towards the development of a real-world Hindi-English machine translation system for doctor-patient communication. As this is a low-resource language pair, especially for speech and for the medical domain, our initial focus has been on gathering suitable training data from various sources. We applied a variety of methods ranging from fully automatic extraction from the Web to manual annotation of test data. Moreover, we developed a method for automatically augmenting the training data with synthetically generated variants, which yielded a very sizable improvement of more than 3 BLEU points absolute.


title: On Deductive Systems of AC Semantics for Rough Sets

abstract: Antichain based semantics for general rough sets were introduced recently by the present author. In her paper two different semantics, one for general rough sets and another for general approximation spaces over quasi-equivalence relations, were developed. These semantics are improved and studied further from a lateral algebraic logic perspective in this research. The main results concern the structure of the algebras and deductive systems in the context.


title: $KMR$ $k_t$-factorization procedure for the description of the $LHCb$ forward hadron-hadron $Z^0$ production at $\sqrt{s}=13\;TeV$

abstract: Quit recently, two sets of new experimental data from the $LHCb$ and the $CMS$ collaborations have been published, concerning the production of the $Z^0$ vector boson in hadron-hadron collisions with the center-of-mass energy $E_{CM}= \sqrt{s}=13\;TeV$. On the other hand, in our recent work, we have conducted a set of $NLO$ calculations for the production of the electroweak gauge vector bosons, utilizing the unintegrated parton distribution functions ($UPDF$) in the frameworks of $Kimber$-$Martin$-$Ryskin$ ($KMR$) or $Martin$-$Ryskin$-$Watt$ ($MRW$) and the $k_t$-factorization formalism, concluding that the results of the $KMR$ scheme are arguably better in describing the existing experimental data, coming from $D0$, $CDF$, $CMS$ and $ATLAS$ collaborations. In the present work, we intend to follow the same $NLO$ formalism and calculate the rate of the production of the $Z^0$ vector boson, utilizing the $UPDF$ of $KMR$ within the dynamics of the recent data. It will be shown that our results are in good agreement with the new measurements of the $LHCb$ and the $CMS$ collaborations.


title: Construction of Regular Black Holes in General Relativity

abstract: We present a general procedure for constructing exact black hole solutions with electric or magnetic charges in General Relativity coupled to a nonlinear electrodynamics. We obtain a variety of two-parameter family spherically symmetric black hole solutions. In particular, the singularity at the central of the space-time can be cancelled in the parameters space and the black hole solutions become regular everywhere in the space-time. We study the global properties of the solutions and derive the first law of thermodynamics. We also generalize the procedure to include a cosmological constant and construct regular black hole solutions that are asymptotic to anti-de Sitter space-time.


title: Two- and Multi-phase Quadrature Surfaces

abstract: In this paper we shall initiate the study of the two- and multi-phase quadrature surfaces (QS), which amounts to a two/multi-phase free boundary problems of Bernoulli type. The problem is studied mostly from a potential theoretic point of view that (for two-phase case) relates to integral representation $$ \int_{\partial \Omega^+} g h (x) \ d\sigma_x - \int_{\partial \Omega^-} g h (x) \ d\sigma_x= \int h d\mu \ , $$ where $d\sigma_x$ is the surface measure, $\mu= \mu^+ - \mu^-$ is given measure with support in (a priori unknown domain) $\Omega$, $g$ is a given smooth positive function, and the integral holds for all functions $h$, which are harmonic on $\overline \Omega$. Our approach is based on minimization of the corresponding two- and multi-phase functional and the use of its one-phase version as a barrier. We prove several results concerning existence, qualitative behavior, and regularity theory for solutions. A central result in our study states that three or more junction points do not appear.


title: A higher rank Racah algebra and the $\mathbb{Z}_2^{n}$ Laplace-Dunkl operator

abstract: A higher rank generalization of the (rank one) Racah algebra is obtained as the symmetry algebra of the Laplace-Dunkl operator associated to the $\mathbb{Z}_2^n$ root system. This algebra is also the invariance algebra of the generic superintegrable model on the $n$-sphere. Bases of Dunkl harmonics are constructed explicitly using a Cauchy-Kovalevskaia theorem. These bases consist of joint eigenfunctions of maximal Abelian subalgebras of the higher rank Racah algebra. A method to obtain expressions for both the connection coefficients between these bases and the action of the symmetries on these bases is presented.


title: Recent Developments in Neutrino-Nucleus Scattering (Theory)

abstract: Selected theoretical developments in neutrino-nucleus scattering in 2015-2016 are reviewed.


title: Lusztig data of Kashiwara-Nakashima tableaux in types $B$ and $C$

abstract: We provide an explicit combinatorial description of the embedding of the crystal of Kashiwara-Nakashima tableaux in types $B$ and $C$ into that of $\bi$-Lusztig data associated to a family of reduced expressions $\bi$ of the longest element $w_0$. The description of the embedding is simple and elementary using only the Sch\"{u}tzenberger's jeu de taquin and RSK algorithm. A spinor model for classical crystals plays an important role as an intermediate object connecting Kashiwara-Nakashima tableaux and Lusztig data.


title: On the dimension of Furstenberg measure for $SL_2(R)$ random matrix products

abstract: Let $\mu$ be a measure on $SL_{2}(\mathbb{R})$ generating a non-compact and totally irreducible subgroup, let $\chi>0$ denote its Lyapunov exponent, and let $\nu$ be the associated stationary (Furstenberg) measure for the action on the projective line. We prove that if $\mu$ is supported on finitely many matrices with algebraic entries, then \[ \dim\nu=\min\{1,\frac{h_{\textrm{RW}}(\mu)}{2\chi}\} \] where $h_{\textrm{RW}}(\mu)$ is the random walk entropy of $\mu$, and $\dim$ denotes pointwise dimension. In particular, for every $\delta>0$, there is a neighborhood $U$ of the identity in $SL_{2}(\mathbb{R})$ such that if a measure $\mu\in\mathcal{P}(U)$ is supported on algebraic matrices with all atoms of size at least $\delta$, and generates a group which is non-compact and totally irreducible, then its stationary measure $\nu$ satisfies $\dim\nu=1$.


title: Left-right model for dark matter

abstract: We argue that dark matter can automatically arise from a gauge theory that possesses a non-minimal left-right gauge symmetry, SU(3)_C \otimes SU(M)_L \otimes SU(N)_R \otimes U(1)_X, for (M,N) = (2,3), (3,2), (3,3), \cdots, and (5,5).


title: A model independent safeguard for unbinned Likelihood

abstract: We present a universal method to include residual un-modeled background shape uncertainties in likelihood based statistical tests for high energy physics and astroparticle physics. This approach provides a simple and natural protection against mismodeling, thus lowering the chances of a false discovery or of an over constrained confidence interval, and allows a natural transition to unbinned space. Unbinned likelihood allows optimal usage of information for the data and the models, and enhances the sensitivity. We show that the asymptotic behavior of the test statistic can be regained in cases where the model fails to describe the true background behavior, and present 1D and 2D case studies for model-driven and data-driven background models. The resulting penalty on sensitivities follows the actual discrepancy between the data and the models, and is asymptotically reduced to zero with increasing knowledge.


title: Quenching of Satellite Galaxies at the Outskirts of Galaxy Clusters

abstract: We find, using cosmological simulations of galaxy clusters, that the hot X-ray emitting intra-cluster medium (ICM) enclosed within the outer accretion shock extends out to $R_{\rm shock}\sim(2 - 3) R_{\rm vir}$, where $R_{\rm vir}$ is the standard virial radius of the halo. Using a simple analytic model for satellite galaxies in the cluster, we evaluate the effect of ram pressure stripping on the gas in the inner discs and in the haloes at different distances from the cluster centre. We find that significant removal of star-forming disc gas occurs only at $r \lesssim 0.5 R_{\rm vir}$, while gas removal from the satellite halo is more effective and can occur when the satellite is found between $ R_{\rm vir}$ and $R_{\rm shock}$. Removal of halo gas sets the stage for quenching of the star formation by starvation over $2\textrm{--}3\,\mathrm{Gyr}$, prior to the satellite entry to the inner cluster halo. This scenario explains the presence of quenched galaxies, preferentially discs, at the outskirts of galaxy clusters, and the delayed quenching of satellites compared to central galaxies.


title: On $L$-packets and depth for $SL_2(K)$ and its inner form

abstract: We consider the group $SL_2(K)$, where $K$ is a local non-archimedean field of characteristic two. We prove that the depth of any irreducible representation of $SL_2 (K)$ is larger than the depth of the corresponding Langlands parameter, with equality if and only if the L-parameter is essentially tame. We also work out a classification of all $L$-packets for $SL_2 (K)$ and for its non-split inner form, and we provide explicit formulae for the depths of their $L$-parameters.


title: Coprimeness-preserving non-integrable extension to the two-dimensional discrete Toda lattice equation

abstract: We introduce a so-called `coprimeness-preserving non-integrable' extension (another terminology is `quasi-integrable' extension) to the two-dimensional Toda lattice equation. We believe that this equation is the first example of such discrete equation defined over a three-dimensional lattice. We prove that all the iterates of the equation are irreducible Laurent polynomials of the initial data and that every pair of two iterates is co-prime, which indicate confined singularities of the equation. By reducing the equation to two- or one-dimensional lattices, we obtain coprimeness-preserving non-integrable extensions to the one-dimensional Toda lattice equation and the Somos-4 recurrence.


title: No Percolation in low temperature spin glass

abstract: We consider the Edwards-Anderson Ising Spin Glass model for non negative temperatures T: We define the natural notion of Boltzmann- Gibbs measure for the Edwards-Anderson spin glass at a given temperature, and of unsatisfied edges. We prove that for low enough temperatures, in almost every spin configuration the graph formed by the unsatisfied edges is made of finite connected components. In other words, the unsatisfied edges do not percolate.


title: Multi-wavelength observations of IGR J17544-2619 from quiescence to outburst

abstract: In this paper we report on a long multi-wavelength observational campaign of the supergiant fast X-ray transient prototype IGR J17544-2619. A 150 ks-long observation was carried out simultaneously with XMM-Newton and NuSTAR, catching the source in an initial faint X-ray state and then undergoing a bright X-ray outburst lasting about 7 ks. We studied the spectral variability during outburst and quiescence by using a thermal and bulk Comptonization model that is typically adopted to describe the X-ray spectral energy distribution of young pulsars in high mass X-ray binaries. Although the statistics of the collected X-ray data were relatively high we could neither confirm the presence of a cyclotron line in the broad-band spectrum of the source (0.5-40 keV), nor detect any of the previously reported tentative detection of the source spin period. The monitoring carried out with Swift/XRT during the same orbit of the system observed by XMM-Newton and NuSTAR revealed that the source remained in a low emission state for most of the time, in agreement with the known property of all supergiant fast X-ray transients being significantly sub-luminous compared to other supergiant X-ray binaries. Optical and infrared observations were carried out for a total of a few thousands of seconds during the quiescence state of the source detected by XMM-Newton and NuSTAR. The measured optical and infrared magnitudes were slightly lower than previous values reported in the literature, but compatible with the known micro-variability of supergiant stars. UV observations obtained with the UVOT telescope on-board Swift did not reveal significant changes in the magnitude of the source in this energy domain compared to previously reported values.


title: A new selection strategy for selective cluster ensemble based on Diversity and Independency

abstract: This research introduces a new strategy in cluster ensemble selection by using Independency and Diversity metrics. In recent years, Diversity and Quality, which are two metrics in evaluation procedure, have been used for selecting basic clustering results in the cluster ensemble selection. Although quality can improve the final results in cluster ensemble, it cannot control the procedures of generating basic results, which causes a gap in prediction of the generated basic results' accuracy. Instead of quality, this paper introduces Independency as a supplementary method to be used in conjunction with Diversity. Therefore, this paper uses a heuristic metric, which is based on the procedure of converting code to graph in Software Testing, in order to calculate the Independency of two basic clustering algorithms. Moreover, a new modeling language, which we called as "Clustering Algorithms Independency Language" (CAIL), is introduced in order to generate graphs which depict Independency of algorithms. Also, Uniformity, which is a new similarity metric, has been introduced for evaluating the diversity of basic results. As a credential, our experimental results on varied different standard data sets show that the proposed framework improves the accuracy of final results dramatically in comparison with other cluster ensemble methods.


title: Strongly scale-dependent CMB dipolar asymmetry from super-curvature fluctuations

abstract: We reconsider the observed CMB dipolar asymmetry in the context of open inflation, where a supercurvature mode might survive the bubble nucleation. If such a supercurvature mode modulates the amplitude of the curvature power spectrum, it would easily produce an asymmetry in the power spectrum. We show that current observational data can be accommodated in a three-field model, with simple quadratic potentials and a non-trivial field-space metric. Despite the presence of three fields, we believe this model is so far the simplest that can match current observations. We are able to match the observed strong scale dependence of the dipolar asymmetry, without a fine tuning of initial conditions, breaking slow roll or adding a feature to the evolution of any field.


title: Zero Shot Hashing

abstract: This paper provides a framework to hash images containing instances of unknown object classes. In many object recognition problems, we might have access to huge amount of data. It may so happen that even this huge data doesn't cover the objects belonging to classes that we see in our day to day life. Zero shot learning exploits auxiliary information (also called as signatures) in order to predict the labels corresponding to unknown classes. In this work, we attempt to generate the hash codes for images belonging to unseen classes, information of which is available only through the textual corpus. We formulate this as an unsupervised hashing formulation as the exact labels are not available for the instances of unseen classes. We show that the proposed solution is able to generate hash codes which can predict labels corresponding to unseen classes with appreciably good precision.


title: Finding dominant transition pathways via global optimization of action

abstract: We present a new computational approach, Action-CSA, to sample multiple reaction pathways with fixed initial and final states through global optimization of the Onsager-Machlup action using the conformational space annealing method. This approach successfully samples not only the most dominant pathway but also many other possible paths without initial guesses on reaction pathways. Pathway space is efficiently sampled by crossover operations of a set of paths and preserving the diversity of sampled pathways. The sampling ability of the approach is assessed by finding pathways for the conformational changes of alanine dipeptide and hexane. The benchmarks demonstrate that the rank order and the transition time distribution of multiple pathways identified by the new approach are in good agreement with those of long molecular dynamics simulations. We also show that the lowest action folding pathway of the mini-protein FSD-1 identified by the new approach is consistent with previous molecular dynamics simulations and experiments.


title: Lasso-based forecast combinations for forecasting realized variances

abstract: Volatility forecasts are key inputs in financial analysis. While lasso based forecasts have shown to perform well in many applications, their use to obtain volatility forecasts has not yet received much attention in the literature. Lasso estimators produce parsimonious forecast models. Our forecast combination approach hedges against the risk of selecting a wrong degree of model parsimony. Apart from the standard lasso, we consider several lasso extensions that account for the dynamic nature of the forecast model. We apply forecast combined lasso estimators in a comprehensive forecasting exercise using realized variance time series of ten major international stock market indices. We find the lasso extended "ordered lasso" to give the most accurate realized variance forecasts. Multivariate forecast models, accounting for volatility spillovers between different stock markets, outperform univariate forecast models for longer forecast horizons.


title: Nonlinear optical properties of undoped and doped with Zr and Nb KTiOPO4 crystals

abstract: The structure of the ferroelectric phase of undoped KTiOPO$_4$ and its solid solutions with zirconium and niobium is studied from first principles within the density functional theory. The second-order nonlinear susceptibility tensor and the spontaneous polarization of these materials are obtained. It is shown that an improvement in nonlinear optical properties of KTiOPO$_4$ upon its doping with Zr and Nb cannot be explained by a systematic change in the composition of crystals and is apparently associated with the occurrence of defects. Possible structures of such defects are discussed.


title: Higher iterated Hilbert coefficients of the graded components of bigraded modules

abstract: Let $S=K[x_1,\ldots,x_n]$ be the polynomial ring over the field $K$, and let $I\subset S$ be a graded ideal. It is shown that the higher iterated Hilbert coefficients of the graded $S$-modules $\Tor_i^S(M,I^k)$ and $\Ext^i_S(M,I^k)$ are polynomial functions in $k$, and an upper bound for their degree is given. These results are derived by considering suitable bigraded modules.