title: Fully differential cross sections for $ Li^{2+}$-impact ionization of $Li(2s)$ and $Li(2p)$

abstract: A semiclassical impact parameter version of the continuum distorted wave-Eikonal initial state theory is developed to study the differential ionization of $Li$ atoms in collisions with $Li^{2+}$ ions. Both post and prior forms of the transition amplitude are considered. The fully differential cross sections are calculated for the lithium targets in their ground and their first excited states and for the projectile ions at 16~MeV impact energy. The role of the internuclear interaction as well as the significance of the post-prior discrepancy in the ejected electron spectra are investigated. The obtained results for ejection of the electron into the azimuthal plane are compared with the recent measurements and with their corresponding values obtained using a fully quantum mechanical version of the theory. In most of the cases, the consistency of the present approach with the experimental and the quantum theoretical data is reasonable. However, for $2p$-state ionization, in the cases where no experimental data exist, there is a considerable difference between the two theoretical approaches. This difference is questionable and further experiments are needed to judge which theory makes a more accurate description of the collision dynamics.


title: Twist-3 GPDs in Deeply-Virtual Compton Scattering

abstract: The sub-leading power of the scattering amplitude for deeply-virtual Compton scattering (DVCS) off the nucleon contains leading-twist and twist-3 generalized parton distributions (GPDs). We point out that in DVCS, at twist-3 accuracy, one cannot address any individual twist-3 GPD. This complication appears on top of the deconvolution issues familiar from the twist-2 DVCS amplitude. Accessible are exclusively linear combinations involving both vector and axial-vector twist-3 GPDs. This implies, in particular, that the (kinetic) orbital angular momentum of quarks can hardly be constrained by twist-3 DVCS observables. Moreover, using the quark-target model, we find that twist-3 GPDs can be discontinuous. The discontinuities however cancel in the DVCS amplitude, which further supports the hypothesis of factorization at twist-3 accuracy.


title: Cone construction II

abstract: This is the second part of two parts, titled " cone construction". In this part we prove the Lefschetz cohomologicity of the cone operator $Con$.


title: Lower-hybrid oscillations in a cold magnetized electron-positron-ion plasma

abstract: In this paper, we obtain the dispersion relation for the lower-hybrid oscillations/waves in a cold magnetized electron-positron-ion (e-p-i) plasma. It is found that the frequency of the lower-hybrid oscillations in an e-p-i plasma is approximate $1/\sqrt{3}$ times the frequency of lower-hybrid oscillations in an electron-ion (e-i) plasma, provided the plasma density is high, and the equilibrium densities of the positron and the ion are balanced. The present work may have some relevance for laboratory/astrophysical e-p-i plasmas.


title: Backlash Identification in Two-Mass Systems by Delayed Relay Feedback

abstract: Backlash, also known as mechanical play, is a piecewise differentiable nonlinearity which exists in several actuated systems, comprising, e.g., rack-and-pinion drives, shaft couplings, toothed gears, and other machine elements. Generally, the backlash is nested between the moving parts of a complex dynamic system, which handicaps its proper detection and identification. A classical example is the two-mass system which can approximate numerous mechanisms connected by a shaft (or link) with relatively high stiffness and backlash in series. Information about the presence and extent of the backlash is seldom exactly known and is rather conditional upon factors such as wear, fatigue and incipient failures in the components. This paper proposes a novel backlash identification method using one-side sensing of a two-mass system. The method is based on the delayed relay operator in feedback that allows stable and controllable limit cycles to be induced and operated within the (unknown) backlash gap. The system model, with structural transformations required for the one-side backlash measurements, is given, along with the analysis of the delayed relay in velocity feedback. Experimental evaluations are shown for a two-inertia motor bench that has coupling with backlash gap of about one degree.


title: A note on the integrability of Hamiltonian 1 : 2 : 2 resonance

abstract: We study the integrability of the Hamiltonian normal form of 1 : 2 : 2 resonance. It is known that this normal form truncated to order three is integrable. The truncated to order four normal form contains too many parameters. For a generic choice of parameters in the normal form up to order four we prove a non-integrability result using Morales-Ramis theory. We also isolate a non-trivial case of integrability.


title: Geometric ergodicity of Polya-Gamma Gibbs sampler for Bayesian logistic regression with a flat prior

abstract: The logistic regression model is the most popular model for analyzing binary data. In the absence of any prior information, an improper flat prior is often used for the regression coefficients in Bayesian logistic regression models. The resulting intractable posterior density can be explored by running Polson et al.'s (2013) data augmentation (DA) algorithm. In this paper, we establish that the Markov chain underlying Polson et al.'s (2013) DA algorithm is geometrically ergodic. Proving this theoretical result is practically important as it ensures the existence of central limit theorems (CLTs) for sample averages under a finite second moment condition. The CLT in turn allows users of the DA algorithm to calculate standard errors for posterior estimates.


title: Solar Type IV bursts at frequencies 10-30 MHz

abstract: The results of the first observations of Type IV bursts at frequencies 10-30 MHz are presented. These observations were carried out at radio telescopes UTR-2 (Kharkov, Ukraine) and URAN-2 (Poltava, Ukraine) during the period 2003-2006. Detection of Type IV bursts in wide band from 10 to 30MHz with high sensitivity and time resolution allowed to study their properties in details. These bursts have fluxes 10-2000s.f.u. at maximum phase. Their durations are about 1-2 hours and even more. Some of Type IV bursts drift from high to low frequencies with drift rates about 10kHz/s. All observed Type IV bursts have fine structures in the form of sub-bursts with durations from 2s to 20s and frequency drift rates in a majority of 1-2MHz/s. In most cases, sub-bursts with negative drift rates were registered. Sometimes sub-bursts in absorption with durations 10-200s against Type IV burst background have been observed. The Type IV burst observed on July 22, 2004 had zebra structure, in which single zebra stripes had positive, negative and infinite drift rates.


title: First-order bifurcation detection for dynamic complex networks

abstract: In this paper, we explore how network centrality and network entropy can be used to identify a bifurcation network event. A bifurcation often occurs when a network undergoes a qualitative change in its structure as a response to internal changes or external signals. In this paper, we show that network centrality allows us to capture important topological properties of dynamic networks. By extracting multiple centrality features from a network for dimensionality reduction, we are able to track the network dynamics underlying an intrinsic low-dimensional manifold. Moreover, we employ von Neumann graph entropy (VNGE) to measure the information divergence between networks over time. In particular, we propose an asymptotically consistent estimator of VNGE so that the cubic complexity of VNGE is reduced to quadratic complexity that scales more gracefully with network size. Finally, the effectiveness of our approaches is demonstrated through a real-life application of cyber intrusion detection.


title: Orthogonal gauge fixing of first order gravity

abstract: We consider the first order connection formulation of 4D general relativity in the "orthogonal" gauge. We show how the partial gauge fixing of the phase space canonical coordinates leads to the appearance of second class constraints in the theory. We employ the gauge unfixing procedure in order to successfully complete the Dirac treatment of the system. While equivalent to the inversion of the Dirac matrix, the gauge unfixing allows us to work directly with the reduced phase space and the ordinary Poisson bracket. At the same time, we explicitly derive the new set of residual first class constraints preserving the partial gauge fixing, which are linear combinations of the original constraints, and these turn out to contain nonlinear terms. While providing an explicit example of how to consistently recast general relativity in a given partial gauge, the main motivation of this classical analysis is the application of the Quantum Reduced Loop Gravity program to a Schwarzschild black hole geometry.


title: Kernel-nulling for a robust direct interferometric detection of extrasolar planets

abstract: Combining the resolving power of long-baseline interferometry with the high-dynamic range capability of nulling still remains the only technique that can directly sense the presence of structures in the innermost regions of extrasolar planetary systems. Ultimately, the performance of any nuller architecture is constrained by the partial resolution of the on-axis star whose light it attempts to cancel out. However from the ground, the effective performance of nulling is dominated by residual time-varying instrumental phase and background errors that keep the instrument off the null. Our work investigates robustness against instrumental phase. We introduce a modified nuller architecture that enables the extraction of information that is robust against piston excursions. Our method generalizes the concept of kernel, now applied to the outputs of the modified nuller so as to make them robust to second order pupil phase error. We present the general method to determine these kernel-outputs and highlight the benefits of this novel approach. We present the properties of VIKiNG: the VLTI Infrared Kernel NullinG, an instrument concept within the Hi-5 framework for the 4-UT VLTI infrastructure that takes advantage of the proposed architecture, to produce three self-calibrating nulled outputs. Stabilized by a fringe-tracker that would bring piston-excursions down to 50 nm, this instrument would be able to directly detect more than a dozen extrasolar planets so-far detected by radial velocity only, as well as many hot transiting planets and a significant number of very young exoplanets.


title: Complete Intersections of Quadrics and the Weak Lefschetz Property

abstract: We consider artinian algebras $A=\mathbb{C}[x_0,\ldots,x_m]/I$, with $I$ generated by a regular sequence of homogeneous forms of the same degree $d\geq 2$. We show that the multiplication by a general linear form from $A_{d-1}$ to $A_d$ is injective. We prove that the Weak Lefschetz Property holds for artinian complete intersection algebras as above, with $d=2$ and $m\leq 4$. Apparently, this was previously known only for $m\leq 3$. Although we are proposing only very limited progress towards the WLP conjecture for complete intersections, we hope that the methods of the present article can illustrate some geometrical aspects of the general problem.


title: FeTaSb and FeMnTiSb as promising thermoelectric materials: An ab initio approach

abstract: Thermoelectricity in principle provides a pathway to put waste heat to good use. Motivated by this we investigate thermal and electrical transport properties of two new Fe-based Heusler alloys, FeTaSb and FeMnTiSb, by a first principles approach and semiclassical Boltzmann transport theory within the constant relaxation-time approximation. We find a high power factor of \textit{p}-doped FeTaSb, competitive with best performing Heusler alloy FeNbSb at 1100 K. The obtained power factor of \textit{n}-doped FeMnTiSb at room temperature is higher than that of both FeNbSb and FeTaSb. Remarkably, FeMnTiSb can be used for both \textit{n}-type and \textit{p}-type legs in a thermoelectric module. The Seebeck coefficients of the two proposed systems are in line with those of earlier reported Heusler alloys. We also provide conservative estimates of the figure of merit for the two systems. Overall, our findings suggest a high temperature thermoelectric potential of FeTaSb while the low cost FeMnTiSb is a viable room temperature thermoelectric candidate material.


title: Electric-field control of magnetism in few-layered van der Waals magnet

abstract: Manipulating quantum state via electrostatic gating has been intriguing for many model systems in nanoelectronics. When it comes to the question of controlling the electron spins, more specifically, the magnetism of a system, tuning with electric field has been proven to be elusive. Recently, magnetic layered semiconductors have attracted much attention due to their emerging new physical phenomena. However, challenges still remain in the demonstration of a gate controllable magnetism based on them. Here, we show that, via ionic gating, strong field effect can be observed in few-layered semiconducting Cr$_{2}$Ge$_{2}$Te$_{6}$ devices. At different gate doping, micro-area Kerr measurements in the studied devices demonstrate tunable magnetization loops below the Curie temperature, which is tentatively attributed to the moment re-balance in the spin-polarized band structure. Our findings of electric-field controlled magnetism in van der Waals magnets pave the way for potential applications in new generation magnetic memory storage, sensors, and spintronics.


title: Spatiotemporal stabilization proof of concept of Broad Area Semiconductor laser sources

abstract: We provide a numerical proof of concept of a stabilization mechanism for BAS laser sources. The scheme is based on the simultaneous introduction of in-phase two-dimensional modulations on the refractive index and pump (gain). We numerically proof total stabilization in BAS laser sources, both in space and time. We also examine the interplay between the index and gain modulations and the effect of the slow relaxation of carriers on the stabilization performance.


title: REGAL: Representation Learning-based Graph Alignment

abstract: Problems involving multiple networks are prevalent in many scientific and other domains. In particular, network alignment, or the task of identifying corresponding nodes in different networks, has applications across the social and natural sciences. Motivated by recent advancements in node representation learning for single-graph tasks, we propose REGAL (REpresentation learning-based Graph ALignment), a framework that leverages the power of automatically-learned node representations to match nodes across different graphs. Within REGAL we devise xNetMF, an elegant and principled node embedding formulation that uniquely generalizes to multi-network problems. Our results demonstrate the utility and promise of unsupervised representation learning-based network alignment in terms of both speed and accuracy. REGAL runs up to 30x faster in the representation learning stage than comparable methods, outperforms existing network alignment methods by 20 to 30% accuracy on average, and scales to networks with millions of nodes each.


title: A cosmic Zevatron based on cyclotron auto-resonance

abstract: A Tevatron is an accelerator capable of imparting TeV energies to particles like a proton (1 TeV = $10^{12}$ eV). By analogy, a Zevatron is an accelerator scheme envisaged to accelerate particles to ZeV energies (1 ZeV = $10^{21}$ eV). Zevatron schemes have been proposed to explain the acceleration of ultra-high-energy-cosmic-ray (UHECR) particles detected on Earth since 1962. Here we show that nuclei of hydrogen, helium, and iron-56, may be accelerated to ZeV energies, if injected along the common direction of an existing ultra-strong magnetic field and a super-intense radiation field. This can happen if the injection speed, the magnetic field strength, the charge-to-mass ratio of the nucleus, and the radiation field frequency, are such that an auto-resonance condition is fulfilled, whereby the cyclotron frequency of the particle, around the magnetic field lines, matches the Doppler-shifted frequency of the radiation field.


title: Exact and Consistent Interpretation for Piecewise Linear Neural Networks: A Closed Form Solution

abstract: Strong intelligent machines powered by deep neural networks are increasingly deployed as black boxes to make decisions in risk-sensitive domains, such as finance and medical. To reduce potential risk and build trust with users, it is critical to interpret how such machines make their decisions. Existing works interpret a pre-trained neural network by analyzing hidden neurons, mimicking pre-trained models or approximating local predictions. However, these methods do not provide a guarantee on the exactness and consistency of their interpretation. In this paper, we propose an elegant closed form solution named $OpenBox$ to compute exact and consistent interpretations for the family of Piecewise Linear Neural Networks (PLNN). The major idea is to first transform a PLNN into a mathematically equivalent set of linear classifiers, then interpret each linear classifier by the features that dominate its prediction. We further apply $OpenBox$ to demonstrate the effectiveness of non-negative and sparse constraints on improving the interpretability of PLNNs. The extensive experiments on both synthetic and real world data sets clearly demonstrate the exactness and consistency of our interpretation.


title: A Collaborative Computer Aided Diagnosis (C-CAD) System with Eye-Tracking, Sparse Attentional Model, and Deep Learning

abstract: There are at least two categories of errors in radiology screening that can lead to suboptimal diagnostic decisions and interventions:(i)human fallibility and (ii)complexity of visual search. Computer aided diagnostic (CAD) tools are developed to help radiologists to compensate for some of these errors. However, despite their significant improvements over conventional screening strategies, most CAD systems do not go beyond their use as second opinion tools due to producing a high number of false positives, which human interpreters need to correct. In parallel with efforts in computerized analysis of radiology scans, several researchers have examined behaviors of radiologists while screening medical images to better understand how and why they miss tumors, how they interact with the information in an image, and how they search for unknown pathology in the images. Eye-tracking tools have been instrumental in exploring answers to these fundamental questions. In this paper, we aim to develop a paradigm shift CAD system, called collaborative CAD (C-CAD), that unifies both of the above mentioned research lines: CAD and eye-tracking. We design an eye-tracking interface providing radiologists with a real radiology reading room experience. Then, we propose a novel algorithm that unifies eye-tracking data and a CAD system. Specifically, we present a new graph based clustering and sparsification algorithm to transform eye-tracking data (gaze) into a signal model to interpret gaze patterns quantitatively and qualitatively. The proposed C-CAD collaborates with radiologists via eye-tracking technology and helps them to improve diagnostic decisions. The C-CAD learns radiologists' search efficiency by processing their gaze patterns. To do this, the C-CAD uses a deep learning algorithm in a newly designed multi-task learning platform to segment and diagnose cancers simultaneously.


title: Non-degeneracy of cohomological traces for general Landau-Ginzburg models

abstract: We prove non-degeneracy of the cohomological bulk and boundary traces for general open-closed Landau-Ginzburg models associated to a pair $(X,W)$, where $X$ is a non-compact complex manifold with trivial canonical line bundle and $W$ is a complex-valued holomorphic function defined on $X$, assuming only that the critical locus of $W$ is compact (but may not consist of isolated points). These results can be viewed as certain "deformed" versions of Serre duality. The first amounts to a duality property for the hypercohomology of the sheaf Koszul complex of $W$, while the second is equivalent with the statement that a certain power of the shift functor is a Serre functor on the even subcategory of the $\mathbb{Z}_2$-graded category of topological D-branes of such models.


title: Probing Patchy Saturation of Fluids in Nanoporous Media by Ultrasound

abstract: Nanoporous materials provide high surface area per unit mass and are capable of fluids adsorption. While the measurements of overall amount of fluid adsorbed by a nanopororus sample are straightforward, probing the fluid spacial distribution is non-trivial. We consider published data on adsorption and desorption of fluids in nanoporous glasses reported along with the measurements of ultrasonic waves propagation. We analyse these using Biot's theory of dynamic poroelasticity, approximating the patches as spherical shells. Our calculations show that on adsorption the patch diameter is on the order of 10-20 pore diameters, while on desorption the patch size is comparable to the sample size. Our analysis suggests that one can employ ultrasound to probe the uniformity of fluid spatial distribution in nanoporous materials.


title: Stochastic multiscale flux basis for Stokes-Darcy flows

abstract: Three algorithms are developed for uncertainty quantification in modeling coupled Stokes and Darcy flows. The porous media may consist of multiple regions with different properties. The permeability is modeled as a non-stationary stochastic variable, with its log represented as a sum of local Karhunen-Lo\`eve (KL) expansions. The problem is approximated by stochastic collocation on either tensor-product or sparse grids, coupled with a multiscale mortar mixed finite element method for the spatial discretization. A non-overlapping domain decomposition algorithm reduces the global problem to a coarse scale mortar interface problem, which is solved by an iterative solver, for each stochastic realization. In the traditional implementation, each subdomain solves a local Dirichlet or Neumann problem in every interface iteration. To reduce this cost, two additional algorithms based on deterministic or stochastic multiscale flux basis are introduced. The basis consists of the local flux (or velocity trace) responses from each mortar degree of freedom. It is computed by each subdomain independently before the interface iteration begins. The use of the multiscale flux basis avoids the need for subdomain solves on each iteration. The deterministic basis is computed at each stochastic collocation and used only at this realization. The stochastic basis is formed by further looping over all local realizations of a subdomain's KL region before the stochastic collocation begins. It is reused over multiple realizations. Numerical tests are presented to illustrate the performance of the three algorithms, with the stochastic multiscale flux basis showing significant savings in computational cost.


title: Monotonicity in inverse medium scattering on unbounded domains

abstract: We discuss a time-harmonic inverse scattering problem for the Helmholtz equation with compactly supported penetrable and possibly inhomogeneous scattering objects in an unbounded homogeneous background medium, and we develop a monotonicity relation for the far field operator that maps superpositions of incident plane waves to the far field patterns of the corresponding scattered waves. We utilize this monotonicity relation to establish novel characterizations of the support of the scattering objects in terms of the far field operator. These are related to and extend corresponding results known from factorization and linear sampling methods to determine the support of unknown scattering objects from far field observations of scattered fields. An attraction of the new characterizations is that they only require the refractive index of the scattering objects to be above or below the refractive index of the background medium locally and near the boundary of the scatterers. An important tool to prove these results are so-called localized wave functions that have arbitrarily large norm in some prescribed region while at the same time having arbitrarily small norm in some other prescribed region. We present numerical examples to illustrate our theoretical findings.


title: Statistical Link Label Modeling for Sign Prediction: Smoothing Sparsity by Joining Local and Global Information

abstract: One of the major issues in signed networks is to use network structure to predict the missing sign of an edge. In this paper, we introduce a novel probabilistic approach for the sign prediction problem. The main characteristic of the proposed models is their ability to adapt to the sparsity level of an input network. The sparsity of networks is one of the major reasons for the poor performance of many link prediction algorithms, in general, and sign prediction algorithms, in particular. Building a model that has an ability to adapt to the sparsity of the data has not yet been considered in the previous related works. We suggest that there exists a dilemma between local and global structures and attempt to build sparsity adaptive models by resolving this dilemma. To this end, we propose probabilistic prediction models based on local and global structures and integrate them based on the concept of smoothing. The model relies more on the global structures when the sparsity increases, whereas it gives more weights to the information obtained from local structures for low levels of the sparsity. The proposed model is assessed on three real-world signed networks, and the experiments reveal its consistent superiority over the state of the art methods. As compared to the previous methods, the proposed model not only better handles the sparsity problem, but also has lower computational complexity and can be updated using real-time data streams.


title: An analysis of training and generalization errors in shallow and deep networks

abstract: This paper is motivated by an open problem around deep networks, namely, the apparent absence of over-fitting despite large over-parametrization which allows perfect fitting of the training data. In this paper, we analyze this phenomenon in the case of regression problems when each unit evaluates a periodic activation function. We argue that the minimal expected value of the square loss is inappropriate to measure the generalization error in approximation of compositional functions in order to take full advantage of the compositional structure. Instead, we measure the generalization error in the sense of maximum loss, and sometimes, as a pointwise error. We give estimates on exactly how many parameters ensure both zero training error as well as a good generalization error. We prove that a solution of a regularization problem is guaranteed to yield a good training error as well as a good generalization error and estimate how much error to expect at which test data.


title: Gelation of patchy gold nanoparticles decorated by liquid-crystalline ligands: computer simulation study

abstract: We consider patchy gold nanoparticles decorated by liquid crystalline ligands. The cases of two, three, four and six symmetrically arranged patches of ligands are discussed, as well as the cases of their equatorial and uniform arrangement. A solution of decorated nanoparticles is considered within a flat pore with the solid walls and the interior filled by a polar solvent. The ligands form physical crosslinks between the nanoparticles due to strong liquid crystalline interaction, turning the solution into a gel-like structure. Gelation is done repeatedly starting each time from freshly equilibrated dispersed state of nanoparticles. The gelation dynamics and the range of network characteristics of gel are examined, depending on the type of patchy decoration and the solution density. The emphasis is given to the suitability of a gel for catalytic applications


title: McKean-Vlasov diffusion and the well-posedness of the Hookean bead-spring-chain model for dilute polymeric fluids: small-mass limit and equilibration in momentum space

abstract: We reformulate a general class of classical bead-spring-chain models for dilute polymeric fluids, with Hookean spring potentials, as McKean-Vlasov diffusion. This results in a coupled system of partial differential equations involving the unsteady incompressible linearized Navier-Stokes equations, referred to as the Oseen system, for the velocity and the pressure of the fluid, with a source term which is a nonlinear function of the probability density function, and a second-order degenerate parabolic Fokker-Planck equation, whose transport terms depend on the velocity field, for the probability density function. We show that this coupled Oseen-Fokker-Planck system has a large-data global weak solution. We then perform a rigorous passage to the limit as the masses of the beads in the bead-spring-chain converge to zero, which is shown in particular to result in equilibration in momentum space. The limiting problem is then used to perform a rigorous derivation of the Hookean bead-spring-chain model for dilute polymeric fluids, which has the interesting feature that, if the flow domain is bounded, then so is the associated configuration space domain and the associated Kramers stress tensor is defined by integration over this bounded configuration domain.


title: Initial-boundary value problems for multi-term time-fractional diffusion equations with x-dependent coefficients

abstract: In this paper, we discuss an initial-boundary value problem (IBVP) for the multi-term time-fractional diffusion equation with x-dependent coefficients. By means of the Mittag-Leffler functions and the eigenfunction expansion, we reduce the IBVP to an equivalent integral equation to show the unique existence and the analyticity of the solution for the equation. Especially, in the case where all the coefficients of the time-fractional derivatives are non-negative, by the Laplace and inversion Laplace transforms, it turns out that the decay rate of the solution for long time is dominated by the lowest order of the time-fractional derivatives. Finally, as an application of the analyticity of the solution, the uniqueness of an inverse problem in determining the fractional orders in the multi-term time-fractional diffusion equations from one interior point observation is established.


title: MAVIS: Managing Datacenters using Smartphones

abstract: Distributed monitoring plays a crucial role in managing the activities of cloud-based datacenters. System administrators have long relied on monitoring systems such as Nagios and Ganglia to obtain status alerts on their desktop-class machines. However, the popularity of mobile devices is pushing the community to develop datacenter monitoring solutions for smartphone-class devices. Here we lay out desirable characteristics of such smartphone-based monitoring and identify quantitatively the shortcomings from directly applying existing solutions to this domain. Then we introduce a possible design that addresses some of these shortcomings and provide results from an early prototype, called MAVIS, using one month of monitoring data from approximately 3,000 machines hosted by Purdue's central IT organization.


title: Lower Bounds on Sparse Spanners, Emulators, and Diameter-reducing shortcuts

abstract: We prove better lower bounds on additive spanners and emulators, which are lossy compression schemes for undirected graphs, as well as lower bounds on shortcut sets, which reduce the diameter of directed graphs. We show that any $O(n)$-size shortcut set cannot bring the diameter below $\Omega(n^{1/6})$, and that any $O(m)$-size shortcut set cannot bring it below $\Omega(n^{1/11})$. These improve Hesse's [Hesse03] lower bound of $\Omega(n^{1/17})$. By combining these constructions with Abboud and Bodwin's [AbboudB17] edge-splitting technique, we get additive stretch lower bounds of $+\Omega(n^{1/11})$ for $O(n)$-size spanners and $+\Omega(n^{1/18})$ for $O(n)$-size emulators. These improve Abboud and Bodwin's $+\Omega(n^{1/22})$ lower bounds.


title: Superconductivity induced by flexural modes in non $\sigma_{\rm h}$-symmetric Dirac-like two-dimensional materials: A theoretical study for silicene and germanene

abstract: In two-dimensional crystals that lack symmetry under reflections on the horizontal plane of the lattice (non-$\sigma_{\rm h}$-symmetric), electrons can couple to flexural modes (ZA phonons) at first order. We show that in materials of this type that also exhibit a Dirac-like electron dispersion, the strong coupling can result in electron pairing mediated by these phonons, as long as the flexural modes are not damped or suppressed by additional interactions with a supporting substrate or gate insulator. We consider several models: The weak-coupling limit, which is applicable only in the case of gapped and parabolic materials, like stanene and HfSe$_{2}$, thanks to the weak coupling; the full gap-equation, solved using the constant-gap approximation and considering statically screened interactions; its extensions to energy-dependent gap and to dynamic screening. We argue that in the case of silicene and germanene superconductivity mediated by this process can exhibit a critical temperature of a few degrees K, or even a few tens of degrees K when accounting for the effect of a high-dielectric-constant environment. We conclude that the electron/flexural-modes coupling should be included in studies of possible superconductivity in non-$\sigma_{\rm h}$-symmetric two-dimensional crystals, even if alternative forms of coupling are considered.


title: Derivatives of Eisenstein series of weight 2 and intersections of modular correspondences

abstract: We give a formula for certain values and derivatives of Siegel series and use them to compute Fourier coefficients of derivatives of the Siegel Eisenstein series of weight g/2 and genus g. When g=4, the Fourier coefficient is approximated by a certain Fourier coefficient of the central derivative of the Siegel Eisenstein series of weight 2 and genus 3, which is related to the intersection of 3 arithmetic modular correspondences. Applications include a relation between weighted averages of representation numbers of symmetric matrices.


title: Automatic Classification of Roof Shapes for Multicopter Emergency Landing Site Selection

abstract: Geographic information systems (GIS) now provide accurate maps of terrain, roads, waterways, and building footprints and heights. Aircraft, particularly small unmanned aircraft systems, can exploit additional information such as building roof structure to improve navigation accuracy and safety particularly in urban regions. This paper proposes a method to automatically label building roof shape types. Satellite imagery and LIDAR data from Witten, Germany are fed to convolutional neural networks (CNN) to extract salient feature vectors. Supervised training sets are automatically generated from pre-labeled buildings contained in the OpenStreetMap database. Multiple CNN architectures are trained and tested, with the best performing networks providing a condensed feature set for support vector machine and decision tree classifiers. Satellite and LIDAR data fusion is shown to provide greater classification accuracy than through use of either data type individually.


title: Protein motion in the nucleus: from anomalous diffusion to weak interactions

abstract: Understanding how transcription factors (TFs) regulate mammalian gene expression in space and time is a central topic in biology. To activate a gene, a TF has first to diffuse in the available space of the nucleus until it reaches a target DNA sequence or protein (target site). This eventually results in the recruitment of the whole transcriptional machinery. All these processes take place in the mammalian nucleoplasm, a highly organized and dynamic environment, in which some complexes transiently assemble and break apart, whereas others appear more stable. This diversity of dynamic behaviors arises from the number of biomolecules that make up the nucleoplasm and their pairwise interactions. Indeed, interactions energies that span several orders of magnitude, from covalent bounds to transient and dynamic interactions can shape nuclear landscapes. Thus, the nuclear environment determines how frequently and how fast a TF contacts its target site, and indirectly gene expression. How exactly transient interactions are involved in the regulation of TF diffusion is unclear, but are reflected by live cell imaging techniques such as fluorescence correlation spectroscopy, fluorescence recovery after photobleaching or single-particle tracking. Overall, the macroscopic result of these microscopic interactions is almost always anomalous diffusion, a phenomenon widely studied and modeled. Here, we review the connections between the anomalous diffusion of a TF and the microscopic organization of the nucleus, including recently described topologically associated domains and dynamic phase-separated compartments. We propose that anomalous diffusion found in single particle tracking (SPT) data result from weak and transient interactions with dynamic nuclear substructures, and that SPT data analysis would benefit form a better description of such structures.


title: Interactive Estimation of the Fractal Properties of Carbonate Rocks

abstract: Scale invariance of intrinsic patterns is an important concept in geology that can be observed in numerous geological objects and phenomena. These geological objects and phenomena are described as containing statistically selfsimilar patterns often modeled with fractal geometry. Fractal geometry has been used extensively to characterize pore space and fracture distribution of both carbonate and clastic rocks as well as the transport properties of porous media and fluid flow in reservoirs. The fractal properties are usually estimated from thin-section photomicrograph images or scanning electron microscope images. For complex rock such as carbonate rocks, automatic feature detection methods are often inaccurate. In addition, the rocks may be have been subjected to facies selective diagenesis which preferentially affect some of the rock fabric, thus increasing the difficulty in automatic detection of certain features. We present an interactive program, GeoBoxCount, for analyzing thin-section images and calculating the fractal dimension interactively. The program relies on the geologists insight in interpreting the features of interest; this significantly improves the accuracy of feature selection. The program provides two options for calculating the fractal dimension: the Hausdorff and the Minkowsi-Bouligand box-counting methods.


title: Neutron lifetime, dark matter and search for sterile neutrino

abstract: A review is focused on experimental measurements on neutron lifetime. The latest measurements with a gravitational trap (PNPI NRC KI) and a magnetic trap (LANL, USA) confirmed PNPI result of 2005. The results of measurements with storage of ultra cold neutrons are in agreement, yet, there is discrepancy with a beam experiment by $3.5{\sigma}$ (1% of decay probability), which is discussed in literature as "neutron anomaly" along with the ideas of explaining it by decay into dark matter partially. The second part of the paper is devoted to so called "reactor antineutrino anomaly", which refers to deficiency of the measured flux of antineutrino from reactor in respect to the calculated flux by $3{\sigma}$(deviation by 6.6%). Specific feature of the proposal in this paper lies in the fact that both anomalies can be accounted for by one and the same phenomenon of oscillation in baryon sector between a neutron and a neutron of dark matter $n{\rightarrow}n'$ with mass $m_{n'}$, somewhat less than mass $m_n$ of an ordinary neutron. Calculations of the proposed model require one free parameter: mass difference $m_n-m_{n'}$ if one normalizes probability of oscillations for a free neutron on "neutron anomaly" 1%, then, having succeeded to interpret 6.6% of neutron anomaly in calculations, one can determine mass difference. According to preliminary estimations, the mass difference is $m_n-m_{n'}{\approx}$ 3 MeV. However, the analysis of cumulative yields of isotopes occurs in fission fragments was performed and it does not confirm possibility of existence of additional decay channel with emission of dark matter neutron with mass difference $m_n-m_{n'}{\approx}$ 3 MeV. The result of the analysis is the conclusion that for mirror neutrons the region of the mass difference $m_n-m_{n'} {\geq}$ 3 MeV is closed. The region of the mass difference $m_n-m_{n'}{\leq}$ 2 MeV turned out to be not closed.


title: A discrete extrinsic and intrinsic Dirac operator

abstract: In differential geometry of surfaces the Dirac operator appears intrinsically as a tool to address the immersion problem as well as in an extrinsic flavour (that comes with spin transformations to comformally transfrom immersions) and the two are naturally related. In this paper we consider a corresponding pair of discrete Dirac operators, the latter on discrete surfaces with polygonal faces and normals defined on each face, and show that many key properties of the smooth theory are preserved. In particular, the corresponding spin transformations, conformal invariants for them, and the relation between this operator and its intrinsic counterpart are discussed.


title: Statistical Reasoning: Choosing and Checking the Ingredients, Inferences Based on a Measure of Statistical Evidence with Some Applications

abstract: The features of a logically sound approach to a theory of statistical reasoning are discussed. A particular approach that satisfies these criteria is reviewed. This is seen to involve selection of a model, model checking, elicitation of a prior, checking the prior for bias, checking for prior-data conflict and estimation and hypothesis assessment inferences based on a measure of evidence. A long-standing anomalous example is resolved by this approach to inference and an application is made to a practical problem of considerable importance which, among other novel aspects of the analysis, involves the development of a relevant elicitation algorithm.


title: Spatial field reconstruction with INLA: Application to IFU galaxy data

abstract: Astronomical observations of extended sources, such as cubes of integral field spectroscopy (IFS), encode auto-correlated spatial structures that cannot be optimally exploited by standard methodologies. This work introduces a novel technique to model IFS datasets, which treats the observed galaxy properties as realizations of an unobserved Gaussian Markov random field. The method is computationally efficient, resilient to the presence of low-signal-to-noise regions, and uses an alternative to Markov Chain Monte Carlo for fast Bayesian inference, the Integrated Nested Laplace Approximation (INLA). As a case study, we analyse 721 IFS data cubes of nearby galaxies from the CALIFA and PISCO surveys, for which we retrieve the maps of the following physical properties: age, metallicity, mass and extinction. The proposed Bayesian approach, built on a generative representation of the galaxy properties, enables the creation of synthetic images, recovery of areas with bad pixels, and an increased power to detect structures in datasets subject to substantial noise and/or sparsity of sampling. A snippet code to reproduce the analysis of this paper is available in the COIN toolbox, together with the field reconstructions of the CALIFA and PISCO samples.


title: Representations of the inverse hull of a 0-left cancellative semigroup

abstract: A semigroup S containing a zero element is said to be 0-left cancellative if st = sr \neq 0 implies that t = r. Given such an S we build an inverse semigroup H(S), called the inverse hull of S. Motivated by the study of certain C*-algebras associated to H(S) (a task that we will address in a subsequent article) we carry out a detailed analysis of the spectrum of the idempotent semilattice E(S) of H(S) with a special interest in identifying the ultra-characters. In order to produce examples of characters on E(S), we introduce the notion of "strings" in a semigroup, attempting to make sense of the "infinite paths" which are fundamental in the study of graph C*-algebras. Our strongest results are obtained under the assumption that S admits "least common multiples", but we also touch upon the notion of "finite alignment", motivated by the corresponding notion from the theory of higher rank graphs, and which has also appeared in recent papers by Spielberg and collaborators.


title: Large Rank-Based Models with Common Noise

abstract: For large systems of Brownian particles interacting through their ranks introduced in (Banner, Fernholz, Karatzas, 2005), the empirical cumulative distribution function satisfies a porous medium PDE. However, when we introduce a common noise, the limit is no longer deterministic. Instead, we show that this limit is a solution of a stochastic PDE related to this porous medium PDE. This stochastic PDE is somewhat similar to the equations developed for conservation laws with rough stochastic fluxes (Lions, Perthame, Souganidis, 2013).


title: Almost Sure Productivity

abstract: We define Almost Sure Productivity (ASP), a probabilistic generalization of the productivity condition for coinductively defined structures. Intuitively, a probabilistic coinductive stream or tree is ASP if it produces infinitely many outputs with probability 1. Formally, we define almost sure productivity using a final coalgebra semantics of programs inspired from Kerstan and K\"onig. Then, we introduce a core language for probabilistic streams and trees, and provide two approaches to verify ASP: a sufficient syntactic criterion, and a reduction to model-checking pCTL* formulas on probabilistic pushdown automata. The reduction shows that ASP is decidable for our core language.


title: Similarities on Graphs: Kernels versus Proximity Measures

abstract: We analytically study proximity and distance properties of various kernels and similarity measures on graphs. This helps to understand the mathematical nature of such measures and can potentially be useful for recommending the adoption of specific similarity measures in data analysis.


title: Greening Geographical Power Allocation for Cellular Networks

abstract: Harvesting energy from nature (solar, wind etc.) is envisioned as a key enabler for realizing green wireless networks. However, green energy sources are geographically distributed and the power amount is random which may not enough to power a base station by a single energy site. Burning brown energy sources such as coal and crude oil, though companied with carbon dioxide emission, provides stable power. In this paper, without sacrificing communication quality, we investigate how to perform green energy allocation to abate the dependence on brown energy with hybrid brown and green energy injected in power networks. We present a comprehensive framework to characterize the performance of hybrid green and brown energy empowered cellular network. Novel performance metric "bits/ton\ce{CO2}/Hz" is proposed to evaluate the greenness of the communication network. As green energy is usually generated from distributed geographical locations and is time varying, online geographical power allocation algorithm is proposed to maximize the greenness of communication network considering electricity transmission's physical laws i.e., Ohm's law and Kirchhoff's circuit laws. Simulations show that geographically distributed green energy sources complement each other by improving the communication capacity while saving brown energy consumption. Besides, the penetration of green energy can also help reduce power loss on the transmission breaches.


title: Nonconvex Matrix Factorization from Rank-One Measurements

abstract: We consider the problem of recovering low-rank matrices from random rank-one measurements, which spans numerous applications including covariance sketching, phase retrieval, quantum state tomography, and learning shallow polynomial neural networks, among others. Our approach is to directly estimate the low-rank factor by minimizing a nonconvex quadratic loss function via vanilla gradient descent, following a tailored spectral initialization. When the true rank is small, this algorithm is guaranteed to converge to the ground truth (up to global ambiguity) with near-optimal sample complexity and computational complexity. To the best of our knowledge, this is the first guarantee that achieves near-optimality in both metrics. In particular, the key enabler of near-optimal computational guarantees is an implicit regularization phenomenon: without explicit regularization, both spectral initialization and the gradient descent iterates automatically stay within a region incoherent with the measurement vectors. This feature allows one to employ much more aggressive step sizes compared with the ones suggested in prior literature, without the need of sample splitting.


title: Unsupervised vehicle recognition using incremental reseeding of acoustic signatures

abstract: Vehicle recognition and classification have broad applications, ranging from traffic flow management to military target identification. We demonstrate an unsupervised method for automated identification of moving vehicles from roadside audio sensors. Using a short-time Fourier transform to decompose audio signals, we treat the frequency signature in each time window as an individual data point. We then use a spectral embedding for dimensionality reduction. Based on the leading eigenvectors, we relate the performance of an incremental reseeding algorithm to that of spectral clustering. We find that incremental reseeding accurately identifies individual vehicles using their acoustic signatures.


title: Implementation of Neural Network and feature extraction to classify ECG signals

abstract: This paper presents a suitable and efficient implementation of a feature extraction algorithm (Pan Tompkins algorithm) on electrocardiography (ECG) signals, for detection and classification of four cardiac diseases: Sleep Apnea, Arrhythmia, Supraventricular Arrhythmia and Long Term Atrial Fibrillation (AF) and differentiating them from the normal heart beat by using pan Tompkins RR detection followed by feature extraction for classification purpose .The paper also presents a new approach towards signal classification using the existing neural networks classifiers.


title: Faster Algorithms for Integer Programs with Block Structure

abstract: We consider integer programming problems $\max \{ c^T x : \mathcal{A} x = b, l \leq x \leq u, x \in \mathbb{Z}^{nt}\}$ where $\mathcal{A}$ has a (recursive) block-structure generalizing "$n$-fold integer programs" which recently received considerable attention in the literature. An $n$-fold IP is an integer program where $\mathcal{A}$ consists of $n$ repetitions of submatrices $A \in \mathbb{Z}^{r \times t}$ on the top horizontal part and $n$ repetitions of a matrix $B \in \mathbb{Z}^{s \times t}$ on the diagonal below the top part. Instead of allowing only two types of block matrices, one for the horizontal line and one for the diagonal, we generalize the $n$-fold setting to allow for arbitrary matrices in every block. We show that such an integer program can be solved in time $n^2 t^2 {\phi} \cdot (rs{\Delta})^{\mathcal{O}(rs^2+ sr^2)}$ (ignoring logarithmic factors). Here ${\Delta}$ is an upper bound on the largest absolute value of an entry of $\mathcal{A}$ and ${\phi}$ is the largest binary encoding length of a coefficient of $c$. This improves upon the previously best algorithm of Hemmecke, Onn and Romanchuk that runs in time $n^3t^3 {\phi} \cdot {\Delta}^{\mathcal{O}(t^2s)}$. In particular, our algorithm is not exponential in the number $t$ of columns of $A$ and $B$. Our algorithm is based on a new upper bound on the $l_1$-norm of an element of the "Graver basis" of an integer matrix and on a proximity bound between the LP and IP optimal solutions tailored for IPs with block structure. These new bounds rely on the "Steinitz Lemma". Furthermore, we extend our techniques to the recently introduced "tree-fold IPs", where we again present a more efficient algorithm in a generalized setting.


title: TabVec: Table Vectors for Classification of Web Tables

abstract: There are hundreds of millions of tables in Web pages that contain useful information for many applications. Leveraging data within these tables is difficult because of the wide variety of structures, formats and data encoded in these tables. TabVec is an unsupervised method to embed tables into a vector space to support classification of tables into categories (entity, relational, matrix, list, and non-data) with minimal user intervention. TabVec deploys syntax and semantics of table cells, and embeds the structure of tables in a table vector space. This enables superior classification of tables even in the absence of domain annotations. Our evaluations in four real world domains show that TabVec improves classification accuracy by more than 20% compared to three state of the art systems, and that those systems require significant in domain training to achieve good results.


title: Domination of Sample Maxima and Related Extremal Dependence Measures

abstract: For a given $d$-dimensional distribution function (df) $H$ we introduce the class of dependence measures $ \mu(H,Q) = - \mathbb{E}\{ \ln H(Z_1, \ldots, Z_d)\},$ where the random vector $(Z_1, \ldots, Z_d)$ has df $Q$ which has the same marginal df's as $H$. If both $H$ and $Q$ are max-stable df's, we show that for a df $F$ in the max-domain of attraction of $H$, this dependence measure explains the extremal dependence exhibited by $F$. Moreover we prove that $\mu(H,Q)$ is the limit of the probability that the maxima of a random sample from $F$ is marginally dominated by some random vector with df in the max-domain of attraction of $Q$. We show a similar result for the complete domination of the sample maxima which leads to another measure of dependence denoted by $\lambda(Q,H)$. In the literature $\lambda(H,H)$ with $H$ a max-stable df has been studied in the context of records, multiple maxima, concomitants of order statistics and concurrence probabilities. It turns out that both $\mu(H,Q)$ and $\lambda(Q,H)$ are closely related. If $H$ is max-stable we derive useful representations for both $\mu(H,Q)$ and $\lambda(Q,H)$. Our applications include equivalent conditions for $H$ to be a product df and $F$ to have asymptotically independent components.


title: Nonparametric Estimation of Low Rank Matrix Valued Function

abstract: Let $A:[0,1]\rightarrow\mathbb{H}_m$ (the space of Hermitian matrices) be a matrix valued function which is low rank with entries in H\"{o}lder class $\Sigma(\beta,L)$. The goal of this paper is to study statistical estimation of $A$ based on the regression model $\mathbb{E}(Y_j|\tau_j,X_j) = \langle A(\tau_j), X_j \rangle,$ where $\tau_j$ are i.i.d. uniformly distributed in $[0,1]$, $X_j$ are i.i.d. matrix completion sampling matrices, $Y_j$ are independent bounded responses. We propose an innovative nuclear norm penalized local polynomial estimator and establish an upper bound on its point-wise risk measured by Frobenius norm. Then we extend this estimator globally and prove an upper bound on its integrated risk measured by $L_2$-norm. We also propose another new estimator based on bias-reducing kernels to study the case when $A$ is not necessarily low rank and establish an upper bound on its risk measured by $L_{\infty}$-norm. We show that the obtained rates are all optimal up to some logarithmic factor in minimax sense. Finally, we propose an adaptive estimation procedure based on Lepskii's method and model selection with data splitting which is computationally efficient and can be easily implemented and parallelized.


title: Black-Box Reductions for Parameter-free Online Learning in Banach Spaces

abstract: We introduce several new black-box reductions that significantly improve the design of adaptive and parameter-free online learning algorithms by simplifying analysis, improving regret guarantees, and sometimes even improving runtime. We reduce parameter-free online learning to online exp-concave optimization, we reduce optimization in a Banach space to one-dimensional optimization, and we reduce optimization over a constrained domain to unconstrained optimization. All of our reductions run as fast as online gradient descent. We use our new techniques to improve upon the previously best regret bounds for parameter-free learning, and do so for arbitrary norms.


title: Dynamics of strongly interacting unstable two-solitons for generalized Korteweg-de Vries equations

abstract: We consider the generalized Korteweg-de Vries equation $\partial_t u = -\partial_x(\partial_x^2 u + f(u))$, where $f(u)$ is an odd function of class $C^3$. Under some assumptions on $f$, this equation admits \emph{solitary waves}, that is solutions of the form $u(t, x) = Q_v(x - vt - x_0)$, for $v$ in some range $(0, v_*)$. We study pure two-solitons in the case of the same limit speed, in other words global solutions $u(t)$ such that \begin{equation} \label{eq:abstract} \tag{$\ast$} \lim_{t\to\infty}\|u(t) - (Q_v(\cdot - x_1(t)) \pm Q_v(\cdot - x_2(t)))\|_{H^1} = 0, \qquad \text{with}\quad\lim_{t \to \infty}x_2(t) - x_1(t) = \infty. \end{equation} Existence of such solutions is known for $f(u) = |u|^{p-1}u$ with $p \in \mathbb{Z} \setminus \{5\}$ and $p > 2$. We describe the~dynamical behavior of any solution satisfying \eqref{eq:abstract} under the assumption that $Q_v$ is linearly unstable (which corresponds to $p > 5$ for power nonlinearities). We prove that in this case the sign in \eqref{eq:abstract} is necessarily "$+$", which corresponds to an attractive interaction. We also prove that the~distance $x_2(t) - x_1(t)$ between the solitons equals $\frac{2}{\sqrt v}\log(\kappa t) + o(1)$ for some $\kappa = \kappa(v) > 0$.


title: Resonance line broadened quasilinear (RBQ) model for fast ion distribution relaxation due to Alfv\'enic eigenmodes

abstract: The burning plasma performance is limited by the confinement of the superalfvenic fusion products such as alpha particles and the auxiliary heating ions capable of exciting the Alfv\'enic eigenmodes (AEs). In this work the effect of AEs on fast ions is formulated within the quasi-linear (QL) theory generalized for this problem recently. The generalization involves the resonance line broadened interaction of energetic particles (EP) with AEs supplemented by the diffusion coefficients depending on EP position in the velocity space. A new resonance broadened QL code (or RBQ1D) based on this formulation allowing for EP diffusion in radial direction is built and presented in details. We reduce the wave particle interaction (WPI) dynamics to 1D case when the particle kinetic energy is nearly constant. The diffusion equation for EP distribution evolution is then one dimensional and is solved simultaneously for all particles with the equation for the evolution of the wave angular momentum. The evolution of fast ion constants of motion is governed by the QL diffusion equations which are adapted to find the fast ion distribution function. We make initial applications of the RBQ1D to DIII-D plasma with elevated q-profile where the beam ions show stiff transport properties. AE driven fast ion profile relaxation is studied for validations of the QL approach in realistic conditions of beam ion driven instabilities in DIII-D.


title: Collaborative model based design of automated and robotic agricultural vehicles in the Crescendo Tool

abstract: This paper describes a collaborative modelling approach to automated and robotic agricultural vehicle design. The Cresendo technology allows engineers from different disciplines to collaborate and produce system models. The combined models are called co-models and their execution co-simulation. To support future development efforts a template library of different vehicle and controllers types are provided. This paper describes a methodology to developing co-models from initial problem definition to deployment of the actual system. We illustrate the development methodology with an example development case from the agricultural domain. The case relates to an encountered speed controller problem on a differential driven vehicle, where we iterate through different candidate solutions and end up with an adaptive controller solution based on a combination of classical control and learning feedforward. The second case is an example of combining human control interface and co-simulation of agricultural robotic operation to illustrate collaborative development


title: Vortex interaction with a rough wall formed by a hexagonal lattice of posts

abstract: An experimental study is reported which investigates the head-on collision of a laminar vortex ring of diameter D (Re{\Gamma}= 3000) on a fakir-like surface composed of circular posts of height h*=0.068 located on a planar bed. Lattices of the posts in hexagonal and random distribution (average porosity of e = 0.94 in the layer) are compared to each other with respect to the plain wall. Prior to impact, the vortex ring develops the early state of natural azimuthal instabilities of different mode numbers N=5-7 competing with each other. While impacting with the rough wall, a secondary ring is observed which is pushed outwards and is not wrapped around the primary ring as in flat wall impact. Between both rings of opposite sign vorticity, a strong fluid rebound is induced. The hexagonal lattice causes the rapid growth of further secondary vortex structures in a regular mode number N=6 arrangement at the outer edge of the primary ring in form of six lobes which are aligned with the orientations of preferential pathways in the layer. At the outer tip of the lobes radial wall-jets are generated. Rotating the fakir geometry around the centre of impact also rotates the jets location and direction accordingly. A surface with random lattice of the posts at the same average number density is not able to repeat this observation and no regular secondary flow pattern is visible until full breakdown of the ring. The results show that a tailored arrangement of such posts can be used for near-wall flow control when patterns of preferred pathways in the posts layer lock-on with existing instability modes such as in impacting jet flows or in turbulent boundary layer flows.


title: Unimodality of the independence polynomials of non-regular caterpillars

abstract: The independence polynomial $I(G, x)$ of a graph $G$ is the polynomial in variable $x$ in which the coefficient $a_n$ on $x^n$ gives the number of independent subsets $S \subseteq V(G)$ of vertices of $G$ such that $|S| = n$. $I(G, x)$ is unimodal if there is an index $\mu$ such that that $a_0 \leq a_1 \leq$...$\leq a_{\mu-1} \leq a_{\mu} \geq a_{\mu +1} \geq$...$\geq a_{d-1} \geq a_d$ While the independence polynomials of many families of graphs with highly regular structure are known to be unimodal, little is known about less regularly structured graphs. We analyze the independence polynomials of a large infinite family of trees without regular structure and show that these polynomials are unimodal through a combinatorial analysis of the polynomials coefficients.


title: Robotic design choice overview using co-simulation

abstract: Rapid robotic system development sets a demand for multi-disciplinary methods and tools to explore and compare design alternatives. In this paper, we present collaborative modeling that combines discrete-event models of controller software with continuous-time models of physical robot components. The presented co-modeling method utilized VDM for discrete-event and 20-sim for continuous-time modeling. The collaborative modeling method is illustrated with a concrete example of collaborative model development of a mobile robot animal feeding system. Simulations are used to evaluate the robot model output response in relation to operational demands. The result of the simulations provides the developers with an overview of the impacts of each solution instance in the chosen design space. Based on the solution overview the developers can select candidates that are deemed viable to be deployed and tested on an actual physical robot.


title: Exact and Robust Conformal Inference Methods for Predictive Machine Learning With Dependent Data

abstract: We extend conformal inference to general settings that allow for time series data. Our proposal is developed as a randomization method and accounts for potential serial dependence by including block structures in the permutation scheme. As a result, the proposed method retains the exact, model-free validity when the data are i.i.d. or more generally exchangeable, similar to usual conformal inference methods. When exchangeability fails, as is the case for common time series data, the proposed approach is approximately valid under weak assumptions on the conformity score.


title: Structural Properties of Bichromatic Non-crossing Matchings

abstract: Given a set of $n$ red and $n$ blue points in the plane, we are interested in matching red points with blue points by straight line segments so that the segments do not cross. We develop a range of tools for dealing with the non-crossing matchings of points in convex position. It turns out that the points naturally partition into groups that we refer to as orbits, with a number of properties that prove useful for studying and efficiently processing the non-crossing matchings. Bottleneck matching is such a matching that minimizes the length of the longest segment. Illustrating the use of the developed tools, we solve the problem of finding bottleneck matchings of points in convex position in $O(n^2)$ time. Subsequently, combining our tools with a geometric analysis we design an $O(n)$-time algorithm for the case where the given points lie on a circle. Previously best known results were $O(n^3)$ for points in convex position, and $O(n \log n$) for points on a circle.


title: Improving the accuracy of the fast inverse square root algorithm

abstract: We present improved algorithms for fast calculation of the inverse square root for single-precision floating-point numbers. The algorithms are much more accurate than the famous fast inverse square root algorithm and have the same or similar computational cost. The main idea of our work consists in modifying the Newton-Raphson method and demanding that the maximal error is as small as possible. Such modification is possible when the distribution of Newton-Raphson corrections is not symmetric (e.g., if they are non-positive functions).


title: Links between functions and subdifferentials

abstract: A function in a class $\mathcal{F}(X)$ is said to be subdifferentially determined in $\mathcal{F}(X)$ if it is equal up to an additive constant to any function in $\mathcal{F}(X)$ with the same subdifferential. A function is said to be subdifferentially representable if it can be recovered from a subdifferential. We identify large classes of lower semicontinuous functions that possess these properties.


title: Lagrangian mean curvature flow of Whitney spheres

abstract: It is shown that an equivariant Lagrangian sphere with a positivity condition on its Ricci curvature develops a type-II singularity under the Lagrangian mean curvature flow that rescales to the product of a grim reaper with a flat Lagrangian subspace. In particular this result applies to the Whitney spheres.


title: Machine learning for Internet of Things data analysis: A survey

abstract: Rapid developments in hardware, software, and communication technologies have allowed the emergence of Internet-connected sensory devices that provide observation and data measurement from the physical world. By 2020, it is estimated that the total number of Internet-connected devices being used will be between 25 and 50 billion. As the numbers grow and technologies become more mature, the volume of data published will increase. Internet-connected devices technology, referred to as Internet of Things (IoT), continues to extend the current Internet by providing connectivity and interaction between the physical and cyber worlds. In addition to increased volume, the IoT generates Big Data characterized by velocity in terms of time and location dependency, with a variety of multiple modalities and varying data quality. Intelligent processing and analysis of this Big Data is the key to developing smart IoT applications. This article assesses the different machine learning methods that deal with the challenges in IoT data by considering smart cities as the main use case. The key contribution of this study is presentation of a taxonomy of machine learning algorithms explaining how different techniques are applied to the data in order to extract higher level information. The potential and challenges of machine learning for IoT data analytics will also be discussed. A use case of applying Support Vector Machine (SVM) on Aarhus Smart City traffic data is presented for a more detailed exploration.


title: Learning Data-Driven Objectives to Optimize Interactive Systems

abstract: Effective optimization is essential for interactive systems to provide a satisfactory user experience. However, it is often challenging to find an objective to optimize for. Generally, such objectives are manually crafted and rarely capture complex user needs in an accurate manner. We propose an approach that infers the objective directly from observed user interactions. These inferences can be made regardless of prior knowledge and across different types of user behavior. We introduce interactive system optimization, a novel algorithm that uses these inferred objectives for optimization. Our main contribution is a new general principled approach to optimizing interactive systems using data-driven objectives. We demonstrate the high effectiveness of interactive system optimization over several simulations.


title: Out-of-sample extension of graph adjacency spectral embedding

abstract: Many popular dimensionality reduction procedures have out-of-sample extensions, which allow a practitioner to apply a learned embedding to observations not seen in the initial training sample. In this work, we consider the problem of obtaining an out-of-sample extension for the adjacency spectral embedding, a procedure for embedding the vertices of a graph into Euclidean space. We present two different approaches to this problem, one based on a least-squares objective and the other based on a maximum-likelihood formulation. We show that if the graph of interest is drawn according to a certain latent position model called a random dot product graph, then both of these out-of-sample extensions estimate the true latent position of the out-of-sample vertex with the same error rate. Further, we prove a central limit theorem for the least-squares-based extension, showing that the estimate is asymptotically normal about the truth in the large-graph limit.


title: Nonparametric Testing under Random Projection

abstract: A common challenge in nonparametric inference is its high computational complexity when data volume is large. In this paper, we develop computationally efficient nonparametric testing by employing a random projection strategy. In the specific kernel ridge regression setup, a simple distance-based test statistic is proposed. Notably, we derive the minimum number of random projections that is sufficient for achieving testing optimality in terms of the minimax rate. An adaptive testing procedure is further established without prior knowledge of regularity. One technical contribution is to establish upper bounds for a range of tail sums of empirical kernel eigenvalues. Simulations and real data analysis are conducted to support our theory.


title: Learning Adversarially Fair and Transferable Representations

abstract: In this paper, we advocate for representation learning as the key to mitigating unfair prediction outcomes downstream. Motivated by a scenario where learned representations are used by third parties with unknown objectives, we propose and explore adversarial representation learning as a natural method of ensuring those parties act fairly. We connect group fairness (demographic parity, equalized odds, and equal opportunity) to different adversarial objectives. Through worst-case theoretical guarantees and experimental validation, we show that the choice of this objective is crucial to fair prediction. Furthermore, we present the first in-depth experimental demonstration of fair transfer learning and demonstrate empirically that our learned representations admit fair predictions on new tasks while maintaining utility, an essential goal of fair representation learning.


title: Characterizing and Learning Equivalence Classes of Causal DAGs under Interventions

abstract: We consider the problem of learning causal DAGs in the setting where both observational and interventional data is available. This setting is common in biology, where gene regulatory networks can be intervened on using chemical reagents or gene deletions. Hauser and B\"uhlmann (2012) previously characterized the identifiability of causal DAGs under perfect interventions, which eliminate dependencies between targeted variables and their direct causes. In this paper, we extend these identifiability results to general interventions, which may modify the dependencies between targeted variables and their causes without eliminating them. We define and characterize the interventional Markov equivalence class that can be identified from general (not necessarily perfect) intervention experiments. We also propose the first provably consistent algorithm for learning DAGs in this setting and evaluate our algorithm on simulated and biological datasets.


title: Virtual Element Method: an equilibrium-based stress recovery procedure

abstract: Within the framework of the displacement-based Virtual Element Method (VEM) for plane elasticity a significant problem is represented by an accurate evaluation of the stress field. In particular, in the classical VEM formulation, a suitable operator which maps to the strain field is introduced in order to allow the calculation of the stiffness matrix. The stress field is then computed using that strian field, by using the constitutive law. Considering for example a first-order formulation for a homogeneous material, strains are locally mapped onto constant functions, and stresses are accordingly piecewise constant. However, the virtual displacements might engender more complex strain fields for polygons which are not triangles. In this paper, Recovery by Compatibility in Patches is used in order to mitigate such an effect and, thus, enhance the accuracy of the recovered stress field. The procedure is simple, efficient and can be readily implemented in existing codes. Numerical tests confirm the soundness of the proposed approach.


title: Counting linear extensions of restricted posets

abstract: The classical 1991 result by Brightwell and Winkler states that the number of linear extensions of a poset is #P-complete. We extend this result to posets with certain restrictions. First, we prove that the number of linear extension for posets of height two is #P-complete. Furthermore, we prove that this holds for incidence posets of graphs. Finally, we prove that the number of linear extensions for posets of dimension two is #P-complete.


title: Deformation of a Half-Space from Anelastic Strain Confined in a Tetrahedral Volume

abstract: Deformation in the lithosphere-asthenosphere system can be accommodated by faulting and plastic flow. However, incorporating structural data in models of distributed deformation still represents a challenge. Here, I present solutions for the displacements and stress in a half-space caused by distributed anelastic strain confined in a tetrahedral volume. These solutions form the basis of curvilinear meshes that can adapt to realistic structural settings, such as a mantle wedge corner, a spherical shell around a magma chamber, or an aquifer. I provide computer programs to evaluate them in the cases of anti-plane strain, in-plane strain, and three-dimensional deformation. These tools may prove useful in the modeling of deformation data in tectonics, volcanology, and hydrology.


title: Autonomous Vehicle Speed Control for Safe Navigation of Occluded Pedestrian Crosswalk

abstract: Both humans and the sensors on an autonomous vehicle have limited sensing capabilities. When these limitations coincide with scenarios involving vulnerable road users, it becomes important to account for these limitations in the motion planner. For the scenario of an occluded pedestrian crosswalk, the speed of the approaching vehicle should be a function of the amount of uncertainty on the roadway. In this work, the longitudinal controller is formulated as a partially observable Markov decision process and dynamic programming is used to compute the control policy. The control policy scales the speed profile to be used by a model predictive steering controller.


title: Center of mass and K\"ahler structures

abstract: There is a sequence of positive numbers $\delta_{2n}$, such that for any connected $2n$-dimensional Riemannian manifold $M$, there are two mutually exclusive possibilities: $1)$ There is a complex structure on $M$ making it into a K\"ahler manifold, or $2)$ For any almost complex structure $J$ compatible with the metric, at every point $p\in M$, there is a smooth loop $\gamma$ at $p$ such that $dist(J_p, hol_\gamma^{-1}J_phol_\gamma)> \delta_{2n}$.


title: Projective dimension and regularity of edge ideal of some weighted oriented graphs

abstract: In this paper we provide some exact formulas for the projective dimension and the regularity of edge ideals associated to vertex weighted rooted forests and oriented cycles. As some consequences, we give some exact formulas for the depth of these ideals.


title: Spontaneous emission from radiative chiral nematic liquid crystals at the photonic band gap edge: an investigation into the role of the density of photon states near resonance

abstract: In this article, we investigate the spontaneous emission properties of radiating molecules embedded in a chiral nematic liquid crystal, under the assumption that the electronic transition frequency is close to the photonic edge mode of the structure, i.e. at resonance. We take into account the transition broadening and the decay of electromagnetic field modes supported by the so-called `mirror-less' cavity. We employ the Jaynes-Cummings Hamiltonian to describe the electron interaction with the electromagnetic field, focusing on the mode with the diffracting polarization in the chiral nematic layer. As known in these structures, the density of photon states, calculated via the Wigner method, has distinct peaks on either side of the photonic band gap, which manifests itself as a considerable modification of the emission spectrum. We demonstrate that, near resonance, there are notable differences between the behavior of the density of states and the spontaneous emission profile of these structures. In addition, we examine in some detail the case of the logarithmic peak exhibited in the density of states in 2D photonic structures and obtain analytic relations for the Lamb shift and the broadening of the atomic transition in the emission spectrum. The dynamical behavior of the atom-field system is described by a system of two first order differential equations, solved using the Green's function method and the Fourier transform. The emission spectra are then calculated and compared with experimental data.


title: Large Neighborhood-Based Metaheuristic and Branch-and-Price for the Pickup and Delivery Problem with Split Loads

abstract: We consider the multi-vehicle one-to-one pickup and delivery problem with split loads, a NP-hard problem linked with a variety of applications for bulk product transportation, bike-sharing systems and inventory re-balancing. This problem is notoriously difficult due to the interaction of two challenging vehicle routing attributes, "pickups and deliveries" and "split deliveries". This possibly leads to optimal solutions of a size that grows exponentially with the instance size, containing multiple visits per customer pair, even in the same route. To solve this problem, we propose an iterated local search metaheuristic as well as a branch-and-price algorithm. The core of the metaheuristic consists of a new large neighborhood search, which reduces the problem of finding the best insertion combination of a pickup and delivery pair into a route (with possible splits) to a resource-constrained shortest path and knapsack problem. Similarly, the branch-and-price algorithm uses sophisticated labeling techniques, route relaxations, pre-processing and branching rules for an efficient resolution. Our computational experiments on classical single-vehicle instances demonstrate the excellent performance of the metaheuristic, which produces new best known solutions for 92 out of 93 test instances, and outperforms all previous algorithms. Experimental results on new multi-vehicle instances with distance constraints are also reported. The branch-and-price algorithm produces optimal solutions for instances with up to 20 pickup-and-delivery pairs, and very accurate solutions are found by the metaheuristic.


title: Consensus in Software Engineering: A Cognitive Mapping Study

abstract: Background: Philosophers of science including Collins, Feyerabend, Kuhn and Latour have all emphasized the importance of consensus within scientific communities of practice. Consensus is important for maintaining legitimacy with outsiders, orchestrating future research, developing educational curricula and agreeing industry standards. Low consensus contrastingly undermines a field's reputation and hinders peer review. Aim: This paper aims to investigate the degree of consensus within the software engineering academic community concerning members' implicit theories of software engineering. Method: A convenience sample of 60 software engineering researchers produced diagrams describing their personal understanding of causal relationships between core software engineering constructs. The diagrams were then analyzed for patterns and clusters. Results: At least three schools of thought may be forming; however, their interpretation is unclear since they do not correspond to known divisions within the community (e.g. Agile vs. Plan-Driven methods). Furthermore, over one third of participants do not belong to any cluster. Conclusion: Although low consensus is common in social sciences, the rapid pace of innovation observed in software engineering suggests that high consensus is achievable given renewed commitment to empiricism and evidence-based practice.


title: Mass problem in the Standard Model

abstract: We propose a new $\mathrm{SU}(3)_C \otimes \mathrm{SU}(2)_{L}\otimes \mathrm{U}(1)_{Y}\otimes \mathrm{U}(1)_{X}$ gauge model which is non universal respect to the three fermion families of the Standard Model. We introduce additional one top-like quark, two bottom-like quarks and three right handed neutrinos in order to have an anomaly free theory. We also consider additional three right handed neutrinos which are singlets respect to the gauge symmetry of the model to implement see saw mechanism and give masses to the light neutrinos according to the neutrino oscillation phenomenology. In the context of this horizontal gauge symmetry we find mass ansatz for leptons and quarks. In particular, from the analysis of solar, atmospheric, reactor and accelerator neutrino oscillation experiments, we get the allow region for the Yukawa couplings for the charge and neutral lepton sectors according with the mass squared differences and mixing angles for the two neutrino hierarchy schemes, normal and inverted.


title: The Dangerous Dogmas of Software Engineering

abstract: To legitimize itself as a scientific discipline, the software engineering academic community must let go of its non-empirical dogmas. A dogma is belief held regardless of evidence. This paper analyzes the nature and detrimental effects of four software engineering dogmas - 1) the belief that software has "requirements"; 2) the division of software engineering tasks into analysis, design, coding and testing; 3) the belief that software engineering is predominantly concerned with designing "software" systems; 4) the belief that software engineering follows methods effectively. Deconstructing these dogmas reveals that they each oversimplify and over-rationalize aspects of software engineering practice, which obscures underlying phenomena and misleads researchers and practitioners. Evidenced-based practice is analyzed as a means to expose and repudiate non-empirical dogmas. This analysis results in several novel recommendations for overcoming the practical challenges of evidence-based practice.


title: Spin Solid versus Magnetic Charge Ordered State in Artificial Honeycomb Lattice of Connected Elements

abstract: The nature of magnetic correlation at low temperature in two-dimensional artificial magnetic honeycomb lattice is a strongly debated issue. While theoretical researches suggest that the system will develop a novel zero entropy spin solid state as T --> 0 K, a confirmation to this effect in artificial honeycomb lattice of connected elements is lacking. We report on the investigation of magnetic correlation in newly designed artificial permalloy honeycomb lattice of ultra-small elements, with a typical length of ~ 12 nm, using neutron scattering measurements and temperature dependent micromagnetic simulations. Numerical modeling of the polarized neutron reflectometry data elucidates the temperature dependent evolution of spin correlation in this system. As temperature reduces to ~ 7 K, the system tends to develop novel spin solid state, manifested by the alternating distribution of magnetic vortex loops of opposite chiralities. Experimental results are complemented by temperature dependent micromagnetic simulations that confirm the dominance of spin solid state over local magnetic charge ordered state in the artificial honeycomb lattice with connected elements. Our results enable a direct investigation of novel spin solid correlation in the connected honeycomb geometry of two-dimensional artificial structure.


title: New Description of Evolution of Magnetic Phases in Artificial Honeycomb Lattice

abstract: Artificial magnetic honeycomb lattice provides a two-dimensional archetypal system to explore novel phenomena of geometrically frustrated magnets. According to theoretical reports, an artificial magnetic honeycomb lattice is expected to exhibit several phase transitions to unique magnetic states as a function of reducing temperature. Experimental investigations of permalloy artificial honeycomb lattice of connected ultra- small elements, ~ 12 nm, reveal a more complicated behavior. First, upon cooling the sample to intermediate temperature, T ~ 175 K, the system manifests a non-unique state where the long range order co-exists with short-range magnetic charge order and weak spin ice state. Second, at much lower temperature, T ~ 6 K, the long-range spin solid state exhibits a re-entrant behavior. Both observations are in direct contrast to the present understanding of this system. New theoretical approaches are needed to develop a comprehensive formulation of this two dimensional magnet.


title: Fractional electrical dimensionality in the spin solid phase of artificial honeycomb lattice

abstract: Two-dimensional artificial magnetic honeycomb lattice is at the forefront of research on unconventional magnetic materials. Among the many emergent magnetic phases that are predicted to arise as a function of temperature, the low temperature spin solid phase with zero magnetization and entropy is of special importance. Here, we report an interesting perspective to the consequence of spin solid order in an artificial honeycomb lattice of ultra-small connected elements using electrical dimensionality analysis. At low temperature, $T \leq$ 30 K, the system exhibits a very strong insulating characteristic. The electrical dimensionality analysis of the experimental data reveals a fractional dimensionality of $d$ = 0.6(0.04) in the spin solid phase of honeycomb lattice at low temperature. The much smaller electrical dimension in the spin solid phase, perhaps, underscores the strong insulating behavior in this system. Also, the fractional dimensionality in an otherwise two-dimensional system suggests a non-surface-like electrical transport at low temperature in an artificial honeycomb lattice.


title: Scaling of nonlinear susceptibilities in artificial permalloy honeycomb lattice

abstract: Two-dimensional artificial magnetic honeycomb lattice is predicted to manifest thermodynamic phase transition to the spin solid order ground state at low temperature. Nonlinear susceptibilities are very sensitive to thermodynamic phase transition. We have performed the analysis of nonlinear susceptibility to explore the thermodynamic nature of spin solid phase transition in artificial honeycomb lattice of ultra-small connected permalloy (Ni$_{0.81}$Fe$_{0.19}$) elements, typical length of $\simeq$ 12 nm. The nonlinear susceptibility, $\chi_{n1}$, is found to exhibit an unusual cross-over character in both temperature and magnetic field. The higher order susceptibility $\chi_3$ changes from positive to negative as the system traverses through the spin solid phase transition at $T_s$ = 29 K. Additionally, the static critical exponents, used to test the scaling of $\chi_{n1}$, do not follow the conventional scaling relation. We conclude that the transition to the ground state is not truly thermodynamic, thus raises doubt about the validity of predicted zero entropy state in the spin solid phase.


title: Sensitivity and Bifurcation Analysis of a DAE Model for a Microbial Electrolysis Cell

abstract: Microbial electrolysis cells (MECs) are a promising new technology for producing hydrogen cheaply, efficiently, and sustainably. However, to scale up this technology, we need a better understanding of the processes in the devices. In this effort, we present a differential-algebraic equation (DAE) model of a microbial electrolysis cell with an algebraic constraint on current. We then perform sensitivity and bifurcation analysis for the DAE system. The model can be applied either to batch-cycle MECs or to continuous-flow MECs. We conduct differential-algebraic sensitivity analysis after fitting simulations to current density data for a batch-cycle MEC. The sensitivity analysis suggests which parameters have the greatest influence on the current density at particular times during the experiment. In particular, growth and consumption parameters for exoelectrogenic bacteria have a strong effect prior to the peak current density. An alternative strategy to maximizing peak current density is maintaining a long term stable equilibrium with non-zero current density in a continuous-flow MEC. We characterize the minimum dilution rate required for a stable nonzero current equilibrium and demonstrate transcritical bifurcations in the dilution rate parameter that exchange stability between several curves of equilibria. Specifically, increasing the dilution rate transitions the system through three regimes where the stable equilibrium exhibits (i) competitive exclusion by methanogens, (ii) coexistence, and (iii) competitive exclusion by exolectrogens. Positive long term current production is only feasible in the final two regimes. These results suggest how to modify system parameters to increase peak current density in a batch-cycle MEC or to increase the long term current density equilibrium value in a continuous-flow MEC.


title: Directional and Causal Information Flow in EEG for Assessing Perceived Audio Quality

abstract: In this paper, electroencephalography (EEG) measurements are used to infer change in cortical functional connectivity in response to change in audio stimulus. Experiments are conducted wherein the EEG activity of human subjects is recorded as they listen to audio sequences whose quality varies with time. A causal information theoretic framework is then proposed to measure the information flow between EEG sensors appropriately grouped into different regions of interest (ROI) over the cortex. A new causal bidirectional information (CBI) measure is defined as an improvement over standard directed information measures for the purposes of identifying connectivity between ROIs in a generalized cortical network setting. CBI can be intuitively interpreted as a causal bidirectional modification of directed information, and inherently calculates the divergence of the observed data from a multiple access channel with feedback. Further, we determine the analytical relationship between the different causal measures and compare how well they are able to distinguish between the perceived audio quality. The connectivity results inferred indicate a significant change in the rate of information flow between ROIs as the subjects listen to different audio qualities, with CBI being the best in discriminating between the perceived audio quality, compared to using standard directed information measures.


title: Minimum length RNA folding trajectories

abstract: The Kinfold and KFOLD programs for RNA folding kinetics implement the Gillespie algorithm to generate stochastic folding trajectories from an initial structure s to a target structure t, in which each intermediate secondary structure is obtained from its predecessor by the addition, removal or shift of a single base pair. Define MS2 distance between secondary structures s and t to be the minimum path length to refold s to t, where a move from MS2 is applied in each step. We describe algorithms to compute the shortest MS2 folding trajectory between any two given RNA secondary structures. These algorithms include an optimal integer programming (IP) algorithm, an accurate and efficient near-optimal algorithm, a greedy algorithm, a branch-and-bound algorithm, and an optimal algorithm if one allows intermediate structures to contain pseudoknots. Our optimal IP [resp. near-optimal IP] algorithm maximizes [resp. approximately maximizes] the number of shifts and minimizes [resp. approximately minimizes] the number of base pair additions and removals by applying integer programming to (essentially) solve the minimum feedback vertex set (FVS) problem for the RNA conflict digraph, then applies topological sort to tether subtrajectories into the final optimal folding trajectory. We prove NP-hardness of the problem to determine the minimum barrier energy over all possible MS2 folding pathways, and conjecture that computing the MS2 distance between arbitrary secondary structures is NP-hard. Since our optimal IP algorithm relies on the FVS, known to be NP-complete for arbitrary digraphs, we compare the family of RNA conflict digraphs with the following classes of digraphs (planar, reducible flow graph, Eulerian, and tournament) for which FVS is known to be either polynomial time computable or NP-hard. Source code available at http://bioinformatics.bc.edu/clotelab/MS2distance/.


title: Metrics in projective differential geometry: the geometry of solutions to the metrizability equation

abstract: Pseudo-Riemannian metrics with Levi-Civita connection in the projective class of a given torsion free affine connection can be obtained from (and are equivalent to) the maximal rank solutions of a certain overdetermined projectively invariant differential equation often called the metrizability equation. Dropping this rank assumption we study the solutions to this equation given less restrictive generic conditions on its prolonged system. In this setting we find that the solution stratifies the manifold according to the strict signature (pointwise) of the solution and does this in way that locally generalizes the stratification of a model, where the model is, in each case, a corresponding Lie group orbit decomposition of the sphere. Thus the solutions give curved generalizations of such embedded orbit structures. We describe the smooth nature of the strata and determine the geometries of each of the different strata types; this includes a metric on the open strata that becomes singular at the strata boundary, with the latter a type of projective infinity for the given metric. The approach reveals and exploits interesting highly non-linear relationships between different linear geometric partial differential equations. Apart from their direct significance, the results show that, for the metrizability equation, strong results arising for so-called normal BGG solutions, and the corresponding projective holonomy reduction, extend to a far wider class of solutions. The work also provides new results for the projective compactification of scalar-flat metrics.


title: The Category of Factorization

abstract: We introduce and investigate the category of factorization of a multiplicative, commutative, cancellative, pre-ordered monoid $A$, which we denote $\mathcal{F}(A)$. The objects of $\mathcal{F}(A)$ are factorizations of elements of $A$, and the morphisms in $\mathcal{F}(A)$ encode combinatorial similarities and differences between the factorizations. We pay particular attention to the divisibility pre-order and to the monoid $A=D\setminus\{0\}$ where $D$ is an integral domain. Among other results, we show that $\mathcal{F}(A)$ is a symmetric and strict monoidal category with weak equivalences and compute the associated category of fractions obtained by inverting the weak equivalences. Also, we use this construction to characterize various factorization properties of integral domains: atomicity, unique factorization, and so on.