title: Modelling a Hydrodynamic Instability in Freely Settling Colloidal Gels

abstract: Attractive colloidal dispersions, suspensions of fine particles which aggregate and frequently form a space spanning elastic gel are ubiquitous materials in society with a wide range of applications. The colloidal networks in these materials can exist in a mode of free settling when the network weight exceeds its compressive yield stress. An equivalent state occurs when the network is held fixed in place and used as a filter through which the suspending fluid is pumped. In either scenario, hydrodynamic instabilities leading to loss of network integrity occur. Experimental observations have shown that the loss of integrity is associated with the formation of eroded channels, so-called streamers, through which the fluid flows rapidly. However, the dynamics of growth and subsequent mechanism of collapse remain poorly understood. Here, a phenomenological model is presented that describes dynamically the radial growth of a streamer due to erosion of the network by rapid fluid back flow. The model exhibits a finite-time blowup -- the onset of catastrophic failure in the gel -- due to activated breaking of the inter-colloid bonds. Brownian dynamics simulations of hydrodynamically interacting and settling colloids in dilute gels are employed to examine the initiation and propagation of this instability, which is in good agreement with the theory. The model dynamics are also shown to accurately replicate measurements of streamer growth in two different experimental systems. The predictive capabilities and future improvements of the model are discussed and a stability-state diagram is presented providing insight into engineering strategies for avoiding settling instabilities in networks meant to have long shelf lives.


title: A Sharp Free Surface Finite Volume Method Applied to Gravity Wave Flows

abstract: This paper presents a sharp free surface method for fully nonlinear flow of two immiscible phases for wave propagation problems in the Finite Volume framework. The method resolves a sharp interface between two phases by combining the geometric reconstruction Volume-of-Fluid scheme isoAdvector for accurate advection of the free surface with the Ghost Fluid Method for the consistent treatment of density and pressure gradient discontinuities at the free surface. The method uses a compact computational stencil irrespective of cell shape and is formally second-order accurate in time and space. The primary focus of this work is to present the combined method and verify and validate it for wave-related problems in ocean sciences, marine and coastal engineering, by considering the following test cases: i) wave propagation of a two-dimensional wave with moderate steepness, ii) green water (water-on-deck) simulations for a ship model with violent free surface flow patterns. The method is implemented in OpenFOAM, an open source software for Computational Fluid Dynamics.


title: Electron mirror branch: Observational evidence from "historical" AMPTE-IRM and Equator-S measurements

abstract: Based on now historical magnetic and plasma data and available wave spectra from the AMPTE-IRM spacecraft, and on as well historical Equator-S high-cadence magnetic field observations of mirror modes in the magnetosheath near the dayside magnetopause, we present some observational evidence for a recent theoretical evaluation by Noreen et al. (2017) of the contribution of a global electron temperature anisotropy to the evolution of mirror modes in the high-temperature anisotropic collisionless plasma of the magnetosheath causing a separate electron mirror branch. These old data most probably indicate that signatures of this electron effect on mirror modes had indeed been observed already long ago in magnetic and wave data though had not been recognised as such. Unfortunately either poor time resolution or complete lack of plasma data would have inhibited the confirmation of the notoriously required pressure balance in the electron branch for unambiguous confirmation of a separate electron mirror mode. If confirmed by future high-resolution observations, in both cases the large mirror mode amplitudes suggest that mirror modes escape quasilinear saturation being in a state of weak kinetic plasma turbulence.


title: The Impact of Feedback in Massive Star Formation. II. Lower Star Formation Efficiency at Lower Metallicity

abstract: We conduct a theoretical study of the formation of massive stars over a wide range of metallicities from 1e-5 to 1Zsun and evaluate the star formation efficiencies (SFEs) from prestellar cloud cores taking into account multiple feedback processes. Unlike for simple spherical accretion, in the case of disk accretion feedback processes do not set upper limits on stellar masses. At solar metallicity, launching of magneto-centrifugally-driven outflows is the dominant feedback process to set SFEs, while radiation pressure, which has been regarded to be pivotal, has only minor contribution even in the formation of over-100Msun stars. Photoevaporation becomes significant in over-20Msun star formation at low metallicities of <1e-2Zsun, where dust absorption of ionizing photons is inefficient. We conclude that if initial prestellar core properties are similar, then massive stars are rarer in extremely metal-poor environments of 1e-5 - 1e-3Zsun. Our results give new insight into the high-mass end of the initial mass function and its potential variation with galactic and cosmological environments.


title: GRB: Greedy Routing Protocol with Backtracking for Mobile Ad Hoc Network (Extended Version)

abstract: Routing protocols for Mobile Ad Hoc Networks (MANETs) have been extensively studied for more than fifteen years. Position-based routing protocols route packets towards the destination using greedy forwarding (i.e., an intermediate node forwards packets to a neighbor that is closer to the destination than itself). Different position-based protocols use different strategies to pick the neighbor to forward the packet. If a node has no neighbor that is closer to the destination than itself, greedy forwarding fails. In this case, we say there is void (no neighboring nodes) in the direction of the destination. Different position-based routing protocols use different methods for dealing with voids. In this paper, we use a simple backtracking technique to deal with voids and design a position-based routing protocol called "Greedy Routing Protocol with Backtracking (GRB)". We compare the performance of our protocol with the well known Greedy Perimeter Stateless Routing (GPSR) routing and the Ad-Hoc On-demand Distance Vector (AODV) routing protocol as well as the Dynamic Source Routing (DSR) protocol. Our protocol needs much less routing-control packets than those needed by DSR, AODV, and GPSR. Simulation results also show that our protocol has a higher packet-delivery ratio, lower end-to-end delay, and less hop count on average than AODV.


title: Gamma-widths, lifetimes and fluctuations in the nuclear quasi-continuum

abstract: Statistical $\gamma$-decay from highly excited states is determined by the nuclear level density (NLD) and the $\gamma$-ray strength function ($\gamma$SF). These average quantities have been measured for several nuclei using the Oslo method. For the first time, we exploit the NLD and $\gamma$SF to evaluate the $\gamma$-width in the energy region below the neutron binding energy, often called the quasi-continuum region. The lifetimes of states in the quasi-continuum are important benchmarks for a theoretical description of nuclear structure and dynamics at high temperature. The lifetimes may also have impact on reaction rates for the rapid neutron-capture process, now demonstrated to take place in neutron star mergers.


title: A hybrid inverse problem in the fluorescence ultrasound modulated optical tomography in the diffusive regime

abstract: We investigate a hybrid inverse problem in fluorescence ultrasound modulated optical tomography (fUMOT) in the diffusive regime. We prove that the absorption coefficient of the fluorophores at the excitation frequency and the quantum efficiency coefficient can be uniquely and stably reconstructed from boundary measurement of the photon currents, provided that some background medium parameters are known. Reconstruction algorithms are proposed and numerically implemented as well.


title: Anomalous enhancement of the isospin-violating $\Lambda(1405)$ production by a triangle singularity in $\Lambda_c\rightarrow\pi^+\pi^0\pi^0\Sigma^0$

abstract: The decay of $\Lambda_c^+$ into $\pi^+\pi^0\Lambda(1405)$ with the $\Lambda(1405)$ decay into $\pi^0\Sigma^0$ through a triangle diagram is studied. This process is initiated by $ \Lambda_c^+ \to \pi^+ \bar{K}^*N $, then the $\bar{K}^*$ decays into $\bar{K} \pi$ and $\bar{K} N$ produce the $\Lambda(1405)$ through a triangle loop containing $\bar{K}^*N\bar{K}$ which develops a singularity around $1890$ MeV. This process is prohibited by the isospin symmetry, but the decay into this channel is enhanced by the contribution of the triangle diagram, which is sensitive to the mass of the internal particles. We find a narrow peak in the $\pi^0\Sigma^0$ invariant mass distribution, which originates from the $\Lambda(1405)$ amplitude, but is tied to the mass differences between the charged and neutral $\bar{K}$ or $N$ states. The observation of the unavoidable peak of the triangle singularity in the isospin-violating $\Lambda(1405)$ production would provide further support for hadronic molecular picture of the $\Lambda(1405)$ and further information on the $\bar{K} N $ interaction.


title: Electric Dipole Moments from CP-Violating Scalar Leptoquark Interactions

abstract: We analyze the implications of CP-violating scalar leptoquark (LQ) interactions for experimental probes of parity- and time-reversal violating properties of polar molecules. These systems are predominantly sensitive to the electric dipole moment (EDM) of the electron and nuclear-spin-independent (NSID) electron-nucleon interaction. The LQ model can generate both a tree-level NSID interaction as well as the electron EDM at one-loop order. Including both interactions, we find that the NSID interaction can dominate the molecular response. For moderate values of couplings, the current experimental results give roughly two orders of magnitude stronger limits on the electron EDM than one would otherwise infer from a sole-source analysis.


title: Designing a Micro-Benchmark Suite to Evaluate gRPC for TensorFlow: Early Experiences

abstract: Remote procedure call (RPC) is the backbone of many modern distributed systems. Google's gRPC is one of the most popular open source RPC frameworks available in the community. gRPC is the main communication engine for Google's Deep Learning framework TensorFlow. TensorFlow primarily uses gRPC for communicating tensors and administrative tasks among different processes. Tensor updates during the training phase are communication intensive and thus TensorFlow's performance is heavily dependent on the underlying network and the efficacy of the communication engine. Training deep learning models on TensorFlow can take significant time ranging from several minutes to several hours, even several days. Thus system researchers need to devote a lot of time to understand the impact of communication on the overall performance. Clearly, there is lack of benchmarks available for system researchers. Therefore, we propose TF-gRPC-Bench micro-benchmark suite that enables system researchers to quickly understand the impact of the underlying network and communication runtime on deep learning workloads. To achieve this, we first analyze the characteristics of TensorFlow workload over gRPC by training popular deep learning models. Then, we propose three micro-benchmarks that take account these workload characteristics. In addition, we comprehensively evaluate gRPC with TF-gRPC-Bench micro-benchmark suite on different clusters over Ethernet, IPoIB, and RDMA, and present the results.


title: Phase Retrieval in $\ell_2(\RR)$

abstract: We will review the major results in finite dimensional real phase retrieval for vectors and projections. We then (1)prove that many of these theorems hold in infinite dimensions, (2) give counter-examples to show that many others fail in infinite dimensions, (3)list finite dimensional results are unknown for $\ell_2$.


title: A singular Darboux type theorem and non-integrable projective distributions of degree one

abstract: We prove a singular Darboux type theorem for homogeneous polynomial closed $2$-forms of degree one on $\mathbb{C}^n$. As application, we classify non-integrable codimension one distributions, of degree one, and arbitrary classes on projective spaces.


title: Artificial Bone and Teeth through Controlled Ice Growth in Colloidal Suspensions

abstract: The formation of regular patterns is a common feature of many solidification processes involving cast materials. We describe here how regular patterns can be obtained in porous alumina and hydroxyapatite (HAP) by controlling the freezing of ceramic slurries followed by subsequent ice sublimation and sintering, leading to multilayered porous ceramic structures with homogeneous and well-defined architecture. These porous materials can be infiltrated with a second phase of choice to yield biomimetic nacre-like composites with improved mechanical properties, which could be used for artificial bone and teeth applications. Proper control of the solidification patterns provides powerful means of control over the final functional properties. We discuss the relationships between the experimental results, ice growth fundamentals, the physics of ice and the interaction between inert particles and the solidification front during directional freezing.


title: A Modified Image Comparison Algorithm Using Histogram Features

abstract: This article discuss the problem of color image content comparison. Particularly, methods of image content comparison are analyzed, restrictions of color histogram are described and a modified method of images content comparison is proposed. This method uses the color histograms and considers color locations. Testing and analyzing of based and modified algorithms are performed. The modified method shows 97% average precision for a collection containing about 700 images without loss of the advantages of based method, i.e. scale and rotation invariant.


title: Matlab code for Lyapunov exponents of fractional order systems

abstract: In this paper the Benettin-Wolf algorithm to determine all Lyapunov exponents for a class of fractional-order systems modeled by Caputo's derivative and the corresponding Matlab code are presented. First it is proved that the considered class of fractional-order systems admits the necessary variational system necessary to find the Lyapunov exponents. The underlying numerical method to solve the extended system of fractional order, composed of the initial value problem and the variational system, is the predictor-corrector Adams-Bashforth-Moulton for fractional differential equations. The Matlab program prints and plots the Lyapunov exponents as function of time. Also, the programs to obtain Lyapunov exponents as function of the bifurcation parameter and as function of the fractional order are described. The Matlab program for Lyapunov exponents is developed from an existing Matlab program for Lyapunov exponents of integer order. To decrease the computing time, a fast Matlab program which implements the Adams-Bashforth-Moulton method, is utilized. Four representative examples are considered.


title: Self-Organization and Artificial Life: A Review

abstract: Self-organization has been an important concept within a number of disciplines, which Artificial Life (ALife) also has heavily utilized since its inception. The term and its implications, however, are often confusing or misinterpreted. In this work, we provide a mini-review of self-organization and its relationship with ALife, aiming at initiating discussions on this important topic with the interested audience. We first articulate some fundamental aspects of self-organization, outline its usage, and review its applications to ALife within its soft, hard, and wet domains. We also provide perspectives for further research.


title: The role of the general relativity on icy body reservoirs under the effects of an inner eccentric Jupiter

abstract: Recent studies have analyzed the dynamical evolution of outer small body populations under the effects of an eccentric inner massive perturber. Such outer reservoirs are composed of particles on prograde and retrograde orbits, as well as particles whose orbit flips from prograde to retrograde and back again showing a coupling between the inclination i and the ascending node longitude \Omega (Type-F particles). We analyze the role of the General Relativity (GR) on the dynamics of outer particles under the influence of an inner eccentric Jupiter-mass planet produced by a planetary scattering event. In particular, we study how the GR affects the dynamical evolution of the outer Type-F particles, which experience an eccentric Lidov-Kozai mechanism. We carry out N-body simulations with and without GR effects. When the GR is included, the extreme values of \Omega are obtained for retrograde inclinations, while the minimum and maximum inclinations allowed for Type-F particles increase in comparison with that derived without GR effects. According to this, if the GR is included in the simulations, the range of prograde (retrograde) inclinations of the libration region is reduced (increased) respect to that obtained in absence of GR. We find two new class of particles when the GR effects are included in the simulations. On the one hand, particles whose orbital plane flips from prograde to retrograde and back again without experiencing a coupling between i and \Omega. On the other hand, retrograde particles that show a strong coupling between i and \Omega. We infer that the GR may significantly modify the dynamical properties of the outer reservoirs that evolve under the effects of an eccentric inner perturber.


title: Comparing the Max and Noisy-Or Pooling Functions in Multiple Instance Learning for Weakly Supervised Sequence Learning Tasks

abstract: Many sequence learning tasks require the localization of certain events in sequences. Because it can be expensive to obtain strong labeling that specifies the starting and ending times of the events, modern systems are often trained with weak labeling without explicit timing information. Multiple instance learning (MIL) is a popular framework for learning from weak labeling. In a common scenario of MIL, it is necessary to choose a pooling function to aggregate the predictions for the individual steps of the sequences. In this paper, we compare the "max" and "noisy-or" pooling functions on a speech recognition task and a sound event detection task. We find that max pooling is able to localize phonemes and sound events, while noisy-or pooling fails. We provide a theoretical explanation of the different behavior of the two pooling functions on sequence learning tasks.


title: High-order encoding schemes for floodlight quantum key distribution

abstract: Floodlight quantum key distribution (FL-QKD) has realized a 1.3\,Gbit/s secret-key rate (SKR) over a 10-dB-loss channel against a frequency-domain collective attack [Quantum~Sci.~Technol. {\bf 3}, 025007 (2018)]. It achieved this remarkable SKR by means of binary phase-shift keying (BPSK) of multiple optical modes. Moreover, it did so with available technology, and without space-division or wavelength-division multiplexing. In this paper we explore whether replacing FL-QKD's BPSK modulation with a high-order encoding can further increase that protocol's SKR. First, we show that going to $K$-ary phase-shift keying with $K = 32$ doubles---from 2.0 to 4.5\,Gbit/s---the theoretical prediction from [Phys.~Rev.~A {\bf 94}, 012322 (2016)] for FL-QKD's BPSK SKR on a 50-km-long fiber link. Second, we show that $2d\times 2d$ quadrature amplitude modulation (QAM) does not offer any SKR improvement beyond what its $d=1$ case---which is equivalent to quadrature phase-shift keying---provides.


title: Formation of Planetary Populations I: Metallicity & Envelope Opacity Effects

abstract: We present a comprehensive body of simulations of the formation of exoplanetary populations that incorporate the role of planet traps in slowing planetary migration. The traps we include in our model are the water ice line, the disk heat transition, and the dead zone outer edge. We reduce our model parameter set to two physical parameters: the opacity of the accreting planetary atmospheres ($\kappa_{\rm{env}}$) and a measure of the efficiency of planetary accretion after gap opening ($f_{\rm{max}}$). We perform planet population synthesis calculations based on the initial observed distributions of host star and disk properties - their disk masses, lifetimes, and stellar metallicities. We find the frequency of giant planet formation scales with disk metallicity, in agreement with the observed Jovian planet frequency-metallicity relation. We consider both X-ray and cosmic ray disk ionization models, whose differing ionization rates lead to different dead zone trap locations. In both cases, Jovian planets form in our model out to 2-3 AU, with a distribution at smaller radii dependent on the disk ionization source and the setting of envelope opacity. We find that low values of $\kappa_{\rm{env}}$ (0.001-0.002 cm$^2$ g$^{-1}$) and X-ray disk ionization are necessary to obtain a separation between hot Jupiters near 0.1 AU, and warm Jupiters outside 0.6 AU, a feature present in the data. Our model also produces a large number of super Earths, but the majority are outside of 2 AU. As our model assumes a constant dust to gas ratio, we suggest that radial dust evolution must be taken into account to reproduce the observed super Earth population.


title: Music Genre Classification using Machine Learning Techniques

abstract: Categorizing music files according to their genre is a challenging task in the area of music information retrieval (MIR). In this study, we compare the performance of two classes of models. The first is a deep learning approach wherein a CNN model is trained end-to-end, to predict the genre label of an audio signal, solely using its spectrogram. The second approach utilizes hand-crafted features, both from the time domain and the frequency domain. We train four traditional machine learning classifiers with these features and compare their performance. The features that contribute the most towards this multi-class classification task are identified. The experiments are conducted on the Audio set data set and we report an AUC value of 0.894 for an ensemble classifier which combines the two proposed approaches.


title: Detection of anisotropic particles in levitated optomechanics

abstract: We discuss the detection of an anisotropic particle trapped by an elliptically polarized focused Gaussian laser beam. We obtain the full rotational and translational dynamics, as well as, the measured photo-current in a general-dyne detection. As an example, we discuss a toy model of homodyne detection, which captures the main features typically found in experimental setups.


title: Size and Temperature Transferability of Direct and Local Deep Neural Networks for Atomic Forces

abstract: A direct and local deep learning (DL) model for atomic forces is presented. We demonstrate the model performance in bulk aluminum, sodium, and silicon; and show that its errors are comparable to those found in state-of-the-art machine learning and DL models. We then analyze the model's performance as a function of the number of neighbors included and show that one can ascertain physical attributes of the system from the analysis of the deep learning model's behavior. Finally, we test the size scaling performance of the model, and the transferability between different temperatures, and show that our model performs well in both scaling to larger systems and high-to-low temperature predictability.


title: Linear forms of a given Diophantine type and lattice exponents

abstract: In this paper we prove an existence theorem concerning linear forms of a given Diophantine type and apply it to study the structure of the spectrum of lattice exponents.


title: Inferring telescope polarization properties through spectral lines without linear polarization

abstract: We present a technique to determine the polarization properties of a telescope through observations of spectral lines that have no intrinsic linear polarization signals. For such spectral lines, any observed linear polarization must be induced by the telescope optics. We apply the technique to observations taken with the SPINOR at the DST and demonstrate that we can retrieve the characteristic polarization properties of the DST at three wavelengths of 459, 526, and 615 nm. We determine the amount of crosstalk between the intensity Stokes I and the linear and circular polarization states Stokes Q, U, and V, and between Stokes V and Stokes Q and U. We fit a set of parameters that describe the polarization properties of the DST to the observed crosstalk values. The values for the ratio of reflectivities X and the retardance tau match those derived with the telescope calibration unit within the error bars. Residual crosstalk after applying a correction for the telescope polarization stays at a level of 3-10%. We find that it is possible to derive the parameters that describe the polarization properties of a telescope from observations of spectral lines without intrinsic linear polarization signal. Such spectral lines have a dense coverage (about 50 nm separation) in the visible part of the spectrum (400-615 nm), but none were found at longer wavelengths. Using spectral lines without intrinsic linear polarization is a promising tool for the polarimetric calibration of current or future solar telescopes such as DKIST.


title: 3D network modelling of fracture processes in fibre-reinforced geomaterials

abstract: The width of fracture process zones in geomaterials is commonly assumed to depend on the type of heterogeneity of the material. Still, very few techniques exist, which link the type of heterogeneity to the width of the fracture process zone. Here, fracture processes in geomaterials are numerically investigated with structural network approaches, whereby the heterogeneity in the form of large aggregates and low volume fibres is modelled geometrically as poly-dispersed ellipsoids and mono-dispersed line segments, respectively. The influence of aggregates, fibres and combinations of both on fracture processes in direct tensile tests of periodic cells is investigated. For all studied heterogeneities, the fracture process zone localises at the start of the softening regime into a rough fracture. For aggregates, the width of the fracture process zone is greater than for analyses without aggregates. Fibres also increase the initial width of the fracture process zone and, in addition, result in a widening of this zone due to fibre pull out.


title: Socioeconomic Dependencies of Linguistic Patterns in Twitter: A Multivariate Analysis

abstract: Our usage of language is not solely reliant on cognition but is arguably determined by myriad external factors leading to a global variability of linguistic patterns. This issue, which lies at the core of sociolinguistics and is backed by many small-scale studies on face-to-face communication, is addressed here by constructing a dataset combining the largest French Twitter corpus to date with detailed socioeconomic maps obtained from national census in France. We show how key linguistic variables measured in individual Twitter streams depend on factors like socioeconomic status, location, time, and the social network of individuals. We found that (i) people of higher socioeconomic status, active to a greater degree during the daytime, use a more standard language; (ii) the southern part of the country is more prone to use more standard language than the northern one, while locally the used variety or dialect is determined by the spatial distribution of socioeconomic status; and (iii) individuals connected in the social network are closer linguistically than disconnected ones, even after the effects of status homophily have been removed. Our results inform sociolinguistic theory and may inspire novel learning methods for the inference of socioeconomic status of people from the way they tweet.

categories: cs.CL cs.CY cs.SI physics.soc-ph stat.ML


title: Coulomb Excitation of Deuteron in Peripheral Collisions with a Heavy Ion

abstract: We develop an ab initio, non-perturbative, time-dependent Basis Function (tBF) method to solve the nuclear structure and scattering problems in a unified manner. We apply this method to a test problem: the Coulomb excitation of a trapped deuteron by an impinging heavy ion. The states of the deuteron system are obtained by the ab initio nuclear structure calculation implementing a realistic inter-nucleon interaction with a weak external trap to localize the center of mass and to discretize the continuum. The evolution of the internal state of the deuteron system is directly solved using the equation of motion for the scattering. We analyze the excitation mechanism of the deuteron system by investigating its internal transition probabilities and observables as functions of the exposure time and the incident speed. In this investigation, the dynamics of the Coulomb excitation are revealed by the time evolution of the system's internal charge distribution.


title: Approximation Capabilities of Immersed Finite Element Spaces for Elasticity Interface Problems

abstract: We construct and analyze a group of immersed finite element (IFE) spaces formed by linear, bilinear and rotated Q1 polynomials for solving planar elasticity equation involving interface. The shape functions in these IFE spaces are constructed through a group of approximate jump conditions such that the unisolvence of the bilinear and rotated Q1 IFE shape functions are always guaranteed regardless of the Lam\`e parameters and the interface location. The boundedness property and a group of identities of the proposed IFE shape functions are established. A multi-point Taylor expansion is utilized to show the optimal approximation capabilities for the proposed IFE spaces through the Lagrange type interpolation operators.


title: Performance of interFoam on the simulation of progressive waves

abstract: The performance of interFoam (a widely-used solver within the popular open source CFD package OpenFOAM) in simulating the propagation of a nonlinear (stream function solution) regular wave is investigated in this work, with the aim of systematically documenting its accuracy. It is demonstrated that over time there is a tendency for surface elevations to increase, wiggles to appear in the free surface, and crest velocities to become (severely) overestimated. It is shown that increasing the temporal and spatial resolution can mitigate these undesirable effects, but that a relatively small Courant number is required. It is further demonstrated that the choice of discretization schemes and solver settings (often treated as a "black box" by users) can have a major impact on the results. This impact is documented, and it is shown that obtaining a "diffusive balance" is crucial to accurately propagate a surface wave over long distances without requiring exceedingly high temporal and spatial resolutions. Finally, the new code isoAdvector is compared to interFoam, which is demonstrated to produce comparably accurate results, while maintaining a sharper surface. It is hoped that the systematic documentation of the performance of the interFoam solver will enable its more accurate and optimal use, as well as increase awareness of potential shortcomings, by CFD researchers interested in the general CFD simulation of free surface waves.


title: Crystal Loss and Quality Pooling for Unconstrained Face Verification and Recognition

abstract: In recent years, the performance of face verification and recognition systems based on deep convolutional neural networks (DCNNs) has significantly improved. A typical pipeline for face verification includes training a deep network for subject classification with softmax loss, using the penultimate layer output as the feature descriptor, and generating a cosine similarity score given a pair of face images or videos. The softmax loss function does not optimize the features to have higher similarity score for positive pairs and lower similarity score for negative pairs, which leads to a performance gap. In this paper, we propose a new loss function, called Crystal Loss, that restricts the features to lie on a hypersphere of a fixed radius. The loss can be easily implemented using existing deep learning frameworks. We show that integrating this simple step in the training pipeline significantly improves the performance of face verification and recognition systems. We achieve state-of-the-art performance for face verification and recognition on challenging LFW, IJB-A, IJB-B and IJB-C datasets over a large range of false alarm rates (10-1 to 10-7).


title: Analytical Infrared Limit of Nonlinear Thomson Scattering Including Radiation Reaction

abstract: If an electric charge is accelerated by a sufficiently intense electromagnetic field, the effects of the radiation emitted by the charge on the charge dynamics (radiation reaction) cannot be ignored. Here we show that classical radiation-reaction effects alter qualitatively and quantitatively the infrared behavior of the spectrum of the radiation emitted by an electron in the presence of an intense electromagnetic plane wave (nonlinear Thomson scattering). An analytical expression of the infrared limit of nonlinear Thomson scattering is provided, which includes radiation-reaction effects and is valid for an arbitrary plane wave. Apart from their own conceptual importance and as a signature of classical radiation reaction, these results provide the limiting expression of the corresponding and yet unknown exact infrared behavior of strong-field QED in an intense plane wave.


title: 'Divide-and-conquer' semiclassical molecular dynamics: An application to water clusters

abstract: We present an investigation of vibrational features in water clusters performed by means of our recently established divide-and-conquer semiclassical approach [M. Ceotto, G. Di Liberto, and R. Conte, Phys. Rev. Lett. 119, 010401 (2017)]. This technique allows us to simulate quantum vibrational spectra of high-dimensional systems starting from full-dimensional classical trajectories and projection of the semiclassical propagator onto a set of lower dimensional subspaces. The potential energy surface employed is a many-body representation up to three-body terms, in which monomers and two-body interactions are described by the high level Wang-Huang-Braams-Bowman (WHBB) water potential, while, for three-body interactions, calculations adopt a fast permutationally invariant ab initio surface at the same level of theory of the WHBB 3-body potential. Applications range from the water dimer up to the water decamer, a system made of 84 vibrational degrees of freedom. Results are generally in agreement with previous variational estimates in the literature. This is particularly true for the bending and the high-frequency stretching motions, while estimates of modes strongly influenced by hydrogen bonding are red shifted, in a few instances even substantially, as a consequence of the dynamical and global picture provided by the semiclassical approach.


title: Optical properties of coupled silicon nanowires and unusual mechanical inductions

abstract: A recent study of the photonic coupling between metallic nanowires has revealed new degrees of freedom in the system. Unexpected spin torques were induced on dimers when illuminated with linearly polarized plane-waves. As near-field observables, the spectra of torques showed more resolved resonances than typical far-field spectra. Here the study is extended to silicon dimers. Strong forces and torques are exerted by light under both polarizations s and p, contrary to plasmonic systems where the resonant strong forces are found only for p-polarization. The systems made of high-dielectric possess volume resonances that induce the forces differently than in plasmonic systems, which have surface resonances. The asymmetry in strong near-fields is responsible for the unusual mechanics of the system. Some consequences of that may include the breaking of the action-reaction principle or the appearance of pulling forces. The numerical study is based on an exact method. The work is thought for the design of nanorotators and nanodetectors. It suggests a new viewpoint about optical forces: the resultant dynamics of topological variations of electromagnetic fields.


title: Searching for Supersolidity in Ultracold Atomic Bose Condensates with Rashba Spin-Orbit Coupling

abstract: We developed functional integral formulation for the stripe phase of a spinor Bose-Einstein condensates with Rashba spin-orbit coupling. The excitation spectrum is found to exhibit double gapless band structures, identified to be two Goldstone modes resulting from spontaneously broken internal gauge symmetry and translational invariance symmetry. The sound velocities display anisotropic behaviors with the lower branch vanishes in the direction perpendicular to the stripe in the x-y plane. At the transition point between the plane wave phase and the stripe phase, physical quantities such as fluctuation correction to the ground state energy and quantum depletion of the condensates exhibit discontinuity, characteristic of the first order phase transition. Despite strong quantum fluctuations induced by Rashba spin-orbit coupling, we show that the supersolid phase is stable against quantum depletion. Finally we extend our formulation to finite temperatures to account for interactions between excitations.


title: Bohr sets in triple products of large sets in amenable groups

abstract: We answer a question of Hegyv\'ari and Ruzsa concerning effective estimates of the Bohr-regularity of certain triple sums of sets with positive upper Banach densities in the integers. Our proof also works for any discrete amenable group, and it does not require all addends in the triple products we consider to have positive (left) upper Banach densities; one of the addends is allowed to only have positive upper asymptotic density with respect to a (possibly very sparse) ergodic sequence.


title: The UV spectrum of the Ultracool Dwarf LSR J1835+3259 observed with the Hubble Space Telescope

abstract: An interesting question about ultracool dwarfs recently raised in the literature is whether their emission is purely internally driven or partially powered by external processes similar to planetary aurora known from the solar system. In this work we present Hubble Space Telescope observations of the energy fluxes of the M8.5 ultracool dwarf LSR J1835+3259 throughout the UV. The obtained spectra reveal that the object is generally UV-fainter compared to other earlier-type dwarfs. We detect the \ion{Mg}{2} doublet at 2800 \A and constrain an average flux throughout the Near-UV. In the Far-UV without Lyman alpha, the ultracool dwarf is extremely faint with an energy output at least a factor of 1000 smaller as expected from auroral emission physically similar to that on Jupiter. We also detect the red wing of the Lyman alpha emission. Our overall finding is that the observed UV spectrum of LSR J1835+3259 resembles the spectrum of mid/late-type M-dwarf stars relatively well, but it is distinct from a spectrum expected from Jupiter-like auroral processes.


title: Driving Solar Giant Cells through the Self-Organization of Near-Surface Plumes

abstract: Global 3D simulations of solar giant-cell convection have provided significant insight into the processes which yield the Sun's observed differential rotation and cyclic dynamo action. However, as we move to higher resolution simulations a variety of codes have encountered what has been termed the convection conundrum. As these simulations increase in resolution and hence the level of turbulence achieved, they tend to produce weak or even anti-solar differential rotation patterns associated with a weak rotational influence (high Rossby number) due to large convective velocities. One potential culprit for this convection conundrum is the upper boundary condition applied in most simulations which is generally impenetrable. Here we present an alternative stochastic plume boundary condition which imposes small-scale convective plumes designed to mimic near-surface convective downflows, thus allowing convection to carry the majority of the outward solar energy flux up to and through our simulated upper boundary. The use of a plume boundary condition leads to significant changes in the convective driving realized in the simulated domain and thus to the convective energy transport, the dominant scale of the convective enthalpy flux, and the relative strength of the strongest downflows, the downflow network, and the convective upflows. These changes are present even far from the upper boundary layer. Additionally, we demonstrate that in spite of significant changes, giant cell morphology in the convective patterns is still achieved with self-organization of the imposed boundary plumes into downflow lanes, cellular patterns, and even rotationally-aligned banana cells in equatorial regions. This plume boundary presents an alternative pathway for 3D global convection simulations where driving is non-local and may provide a new approach towards addressing the convection conundrum.


title: Excitation properties of the divacancy in 4H-SiC

abstract: We investigate the quenching of the photoluminescence (PL) from the divacancy defect in 4H-SiC consisting of a nearest-neighbour silicon and carbon vacancies. The quenching occurs only when the PL is excited below certain photon energies (thresholds), which differ for the four different inequivalent divacancy configurations in 4H-SiC. Refined theoretical ab initio calculation for the charge-transfer levels of the divacancy show very good agreement between the position of the (0/-) level with respect to the conduction band for each divacancy configurations and the corresponding experimentally observed threshold, allowing us to associate the PL decay with conversion of the divacancy from neutral to negative charge state due to capture of electrons photoionized from other defects (traps) by the excitation. Electron paramagnetic resonance measurements are conducted in dark and under excitation similar to that used in the PL experiments and shed light on the possible origin of traps in the different samples. A simple model built on this concept agrees well with the experimentally-observed decay curves.


title: Cokernels of the Cartan Matrix and Stratifying Systems

abstract: We study the cokernel of the application given by the Cartan Matrix $C_\Lambda$ of a finite dimensional $k$-algebra $\Lambda.$ This produces a finitely generated abelian group, the Cartan group $G_\Lambda,$ which is invariant under derived equivalences. We are interested in the case when $G_\Lambda$ is finite. For a standardly stratified algebra, it is shown that this group is always finite and some interesting connections with the standard modules are found. As a consequence, it is got that $G_\Lambda$ can be seen as a measure of how far is a standardly stratified algebra $\Lambda$ to be quasi-hereditary. Finally, it is also shown that any finite abelian group can be realized as the Cartan group of some standardly stratified algebra.


title: Mobility and Clustering of Barium Ions and Dications in High Pressure Xenon Gas

abstract: The clustering and drift properties of barium ions in xenon gas are explored theoretically, using density functional theory and computational ion mobility theory, with the goal of better understanding barium ion transport for neutrinoless double beta decay. We derive the equilibrium conformations, energies and entropies of molecular ions in the Ba$^{+}$-Xe and Ba$^{++}$-Xe systems, which yield a predictive model of cluster formation in high pressure gas. We calculate ion-neutral interaction potential curves for these species and use them to predict effective molecular ion mobilities. Our calculation consistently reproduces experimental data on effective mobility and molecular ion formation for the Ba$^+$ system, and predicts strong cluster formation in the Ba$^{++}$ system, dominated by stable [BaXe$_6$]$^{++}$,[BaXe$_7$]$^{++}$, [BaXe$_8$]$^{++}$ and [BaXe$_9$]$^{++}$ complexes in the range of interest. Some implications for barium tagging in gas-phase neutrinoless double beta decay experiments are discussed, and the first predictions of pressure-dependent mobility of the doubly charged Ba$^{++}$ species are presented.

categories: physics.ins-det hep-ex nucl-ex physics.chem-ph


title: Femtometer Displacement Resolution with Phase-Insensitive Doppler Sensing

abstract: The measurement of extremely small displacements is of utmost importance, both for fundamental studies [1-4], and practical applications [5-7]. One way to estimate a small displacement is to measure the Doppler shift generated in light reflected off an object moving with a known periodic frequency. This remote sensing technique converts a displacement measurement into a frequency measurement, and has been considerably successful [8-14]. The displacement sensitivity of this technique is limited by the Doppler frequency noise floor and by the velocity of the moving object. Other primary limitations are hours of integration time [12,13] and optimal operation only in a narrow Doppler frequency range. Here we show a sensitive device capable of measuring $\mu$Hz/$\sqrt{\text{Hz}}$ Doppler frequency shifts corresponding to tens of fm displacements for a mirror oscillating at 2 Hz. While the Doppler shift measured is comparable to other techniques [12.13], the position sensitivity is orders of magnitude better, and operates over several orders of magnitude of Doppler frequency range. In addition, unlike other techniques which often rely on interferometric methods, our device is phase insensitive, making it unusually robust to noise.


title: Structural Dynamics and Evolution of Capsule Endoscopy (Pill Camera) Technology in Gastroenterologist Assertion

abstract: This research paper examined and re-evaluates the technological innovation, theory, structural dynamics and evolution of Pill Camera (Capsule Endoscopy) technology in redirecting the response manner of small bowel (intestine) examination in human.


title: Applying Computer Algebra Systems with SAT Solvers to the Williamson Conjecture

abstract: We employ tools from the fields of symbolic computation and satisfiability checking---namely, computer algebra systems and SAT solvers---to study the Williamson conjecture from combinatorial design theory and increase the bounds to which Williamson matrices have been enumerated. In particular, we completely enumerate all Williamson matrices of even order up to and including 70 which gives us deeper insight into the behaviour and distribution of Williamson matrices. We find that, in contrast to the case when the order is odd, Williamson matrices of even order are quite plentiful and exist in every even order up to and including 70. As a consequence of this and a new construction for 8-Williamson matrices we construct 8-Williamson matrices in all odd orders up to and including 35. We additionally enumerate all Williamson matrices whose orders are divisible by 3 and less than 70, finding one previously unknown set of Williamson matrices of order 63.


title: Comment on "Minimum Action Path Theory Reveals the Details of Stochastic Transitions Out of Oscillatory States"

abstract: De la Cruz et al. [Phys. Rev. Lett. 120, 128102 (2018); arXiv:1705.08683] studied a noise-induced transition in an oscillating stochastic population undergoing birth- and death-type reactions. They applied the Freidlin-Wentzell WKB formalism to determine the most probable path to the noise-induced escape from a limit cycle predicted by deterministic theory, and to find the probability distribution of escape time. Here we raise a number of objections to their calculations.


title: Towards Deep Learning based Hand Keypoints Detection for Rapid Sequential Movements from RGB Images

abstract: Hand keypoints detection and pose estimation has numerous applications in computer vision, but it is still an unsolved problem in many aspects. An application of hand keypoints detection is in performing cognitive assessments of a subject by observing the performance of that subject in physical tasks involving rapid finger motion. As a part of this work, we introduce a novel hand key-points benchmark dataset that consists of hand gestures recorded specifically for cognitive behavior monitoring. We explore the state of the art methods in hand keypoint detection and we provide quantitative evaluations for the performance of these methods on our dataset. In future, these results and our dataset can serve as a useful benchmark for hand keypoint recognition for rapid finger movements.


title: Algebraic aspects of when and how a Feynman diagram reduces to simpler ones

abstract: The method of Symmetries of Feynman Integrals defines for any Feynman diagram a set of partial differential equations. On some locus in parameter space the equations imply that the diagram can be reduced to a linear combination of simpler diagrams. This paper provides a systematic method to determine this locus and the associated reduction through an algebraic method involving factorization of maximal minors.


title: Looking at Hands in Autonomous Vehicles: A ConvNet Approach using Part Affinity Fields

abstract: In the context of autonomous driving, where humans may need to take over in the event where the computer may issue a takeover request, a key step towards driving safety is the monitoring of the hands to ensure the driver is ready for such a request. This work, focuses on the first step of this process, which is to locate the hands. Such a system must work in real-time and under varying harsh lighting conditions. This paper introduces a fast ConvNet approach, based on the work of original work of OpenPose for full body joint estimation. The network is modified with fewer parameters and retrained using our own day-time naturalistic autonomous driving dataset to estimate joint and affinity heatmaps for driver & passenger's wrist and elbows, for a total of 8 joint classes and part affinity fields between each wrist-elbow pair. The approach runs real-time on real-world data at 40 fps on multiple drivers and passengers. The system is extensively evaluated both quantitatively and qualitatively, showing at least 95% detection performance on joint localization and arm-angle estimation.


title: Statistical Properties of Scale-Invariant Helical Magnetic Fields and Applications to Cosmology

abstract: We investigate the statistical properties of isotropic, stochastic, Gaussian distributed, helical magnetic fields characterized by different shapes of the energy spectra at large length scales and study the associated realizability condition. We discuss smoothed magnetic fields that are commonly used when the primordial magnetic field is constrained by observational data. We are particularly interested in scale-invariant magnetic fields that can be generated during the inflationary stage by quantum fluctuations. We determine the correlation length of such magnetic fields and relate it to the infrared cutoff of perturbations produced during inflation. We show that this scale determines the observational signatures of the inflationary magnetic fields on the cosmic microwave background. At smaller scales, the scale-invariant spectrum changes with time. It becomes a steeper weak-turbulence spectrum at progressively larger scales. We show numerically that the critical length scale where this happens is the turbulent-diffusive scale, which increases with the square root of time.


title: Enhancing terahertz generation from a two-color plasma using optical parametric amplifier waste light

abstract: We show experimentally that the terahertz (THz) emission of a plasma, generated in air by a two-color laser pulse (containing a near IR frequency and its second harmonic), can be enhanced by the addition of an 800-nm pulse. We observed enhancements of the THz electric field by a factor of up to 30. This provides a widely accessible means for researchers using optical parametric amplifiers (OPA) to increase their THz yields by simply adding the residual pump beam of the OPA to the plasma generating beam. We investigate the dependence of the THz electric field enhancement factor on the powers of the two-color beam as well as the 800-nm enhancement beam. Numerical calculations using the well-known photocurrent model are in excellent agreement with the experimental observations.


title: Hypersurfaces of space forms carrying a totally geodesic foliation

abstract: In this paper we give a complete local parametric classification of the hypersurfaces with dimension at least three of a space form that carry a totally geodesic foliation of codimension one. A classification under the assumption that the leaves of the foliation are complete was given in \cite{drt} for Euclidean hypersurfaces. We prove that there exists exactly one further class of local examples in Euclidean space, all of which have rank two. We also extend the classification under the global assumption of completeness of the leaves for hypersurfaces of the sphere and show that there exist plenty of examples in hyperbolic space.


title: Steering random spin systems to speed up the quantum adiabatic algorithm

abstract: A general time-dependent quantum system can be driven fast from its initial ground state to its final ground state without generating transitions by adding a steering term to the Hamiltonian. We show how this technique can be modified to improve on the standard quantum adiabatic algorithm by making a single-particle and cluster approximation to the steering term. The method is applied to a one-dimensional Ising model in a random field. For the limit of strong disorder, the correction terms significantly enhance the probability for the whole system to remain in the ground state for the proposed non-stoquastic annealing protocol. We demonstrate that even when transitions occur for stronger interaction between qubits, the most probable quantum state is one of the lower energy states of the final Hamiltonian. Since the method can be applied to any model, and more sophisticated approximations to the steering term are possible, the new technique opens up an avenue for the improvement of the quantum adiabatic algorithm.


title: Query Shortest Paths Amidst Growing Discs

abstract: The determination of collision-free shortest paths among growing discs has previously been studied for discs with fixed growing rates. Here, we study a more general case of this problem, where: (1) the speeds at which the discs are growing are polynomial functions of degree $\dd$, and (2) the source and destination points are given as query points. We show how to preprocess the $n$ growing discs so that, for two given query points $s$ and $d$, a shortest path from $s$ to $d$ can be found in $O(n^2 \log (\dd n))$ time. The preprocessing time of our algorithm is $O(n^2 \log n + k \log k)$ where $k$ is the number of intersections between the growing discs and the tangent paths (straight line paths which touch the boundaries of two growing discs). We also prove that $k \in O(n^3\dd)$.


title: Predictive and Prescriptive Analytics for Location Selection of Add-on Retail Products

abstract: In this paper, we study an analytical approach to selecting expansion locations for retailers selling add-on products whose demand is derived from the demand of another base product. Demand for the add-on product is realized only as a supplement to the demand of the base product. In our context, either of the two products could be subject to spatial autocorrelation where demand at a given location is impacted by demand at other locations. Using data from an industrial partner selling add-on products, we build predictive models for understanding the derived demand of the add-on product and establish an optimization framework for automating expansion decisions to maximize expected sales. Interestingly, spatial autocorrelation and the complexity of the predictive model impact the complexity and the structure of the prescriptive optimization model. Our results indicate that the models formulated are highly effective in predicting add-on product sales, and that using the optimization framework built on the predictive model can result in substantial increases in expected sales over baseline policies.


title: Quantum teleportation in vacuum using only Unruh-DeWitt detectors

abstract: We consider entanglement extraction into two two-level Unruh-DeWitt detectors from a vacuum of a neutral massless quantum scalar field in a four-dimensional spacetime, where the general monopole coupling to the scalar field is assumed. Based on the reduced density matrix of the two detectors derived within the perturbation theory, we show that the single copy of the entangled pair of the detectors can be utilized in quantum teleportation even when the detectors are separated acausally, while we observe no violation of the Bell-CHSH inequality. In the case of the Minkowski vacuum, in particular, we find that entanglement usable in quantum teleportation is extracted due to the special relativistic effect when the detectors are in a relative inertial motion, while it is not when they are comoving inertially and the switching of the detectors is executed adiabatically at infinite past and future.


title: A parameterization analysis for acoustic full-waveform inversion of sub-wavelength anomalies

abstract: In the case of multi-parameter full-waveform inversion, the computation of the additional Hessian terms that contain derivatives with respect to more than one type of parameter is necessary. If a simple gradient-based minimization is used, different choices of parameterization can be interpreted as different preconditioners that change the condition number of the Hessian. If the non-linear inverse problem is well-posed, then the inversion should converge to a band-limited version of the true solution irrespective of the parameterization choice, provided we start sufficiently close to the global minimum. However, the choice of parameterization will affect the rate of convergence to the exact solution and the best choice of parameterization is the one with the fastest rate. In this paper, we search for the best choice for acoustic multi-parameter full-waveform inversion, where 1. anomalies with a size less than a quarter of the dominant wavelength have to be estimated without the risk of converging to a local minimum; 2. the scattered wavefield is recorded at all the scattering angles; 3. a steepest-descent minimization scheme is used. Our examples suggest that the best choice of parameterization depends on the contrast of the subsurface scatterer that the inversion tries to estimate. Based on the results, we observe that there is no best parameterization choice for full-waveform inversion. We also observe that a parameterization using the acoustic impedance and mass density has the worst convergence rate. Finally, we also show that the parameterization analysis during a hierarchical inversion, where the data have limited scattering angles, only helps to select a subspace for mono-parameter inversion. For multi-parameter hierarchical inversion, the search for the best parameterization in terms of the convergence speed might be obfuscated by non-uniqueness problems.


title: A Mixture Model to Detect Edges in Sparse Co-expression Graphs

abstract: In the early days of gene expression data, researchers have focused on gene-level analysis, and particularly on finding differentially expressed genes. This usually involved making a simplifying assumption that genes are independent, which made likelihood derivations feasible and allowed for relatively simple implementations. In recent years, the scope has expanded to include pathway and `gene set' analysis in an attempt to understand the relationships between genes. We develop a method to recover a gene network's structure from co-expression data, measured in terms of normalized Pearson's correlation coefficients between gene pairs. We treat these co-expression measurements as weights in the complete graph in which nodes correspond to genes. To decide which edges exist in the gene network, we fit a three-component mixture model such that the observed weights of `null edges' follow a normal distribution with mean 0, and the non-null edges follow a mixture of two lognormal distributions, one for positively- and one for negatively-correlated pairs. We show that this so-called L2N mixture model outperforms other methods in terms of power to detect edges, and it allows us to control the false discovery rate. Importantly, the method makes no assumptions about the true network structure.


title: Neural-Guided Deductive Search for Real-Time Program Synthesis from Examples

abstract: Synthesizing user-intended programs from a small number of input-output examples is a challenging problem with several important applications like spreadsheet manipulation, data wrangling and code refactoring. Existing synthesis systems either completely rely on deductive logic techniques that are extensively hand-engineered or on purely statistical models that need massive amounts of data, and in general fail to provide real-time synthesis on challenging benchmarks. In this work, we propose Neural Guided Deductive Search (NGDS), a hybrid synthesis technique that combines the best of both symbolic logic techniques and statistical models. Thus, it produces programs that satisfy the provided specifications by construction and generalize well on unseen examples, similar to data-driven systems. Our technique effectively utilizes the deductive search framework to reduce the learning problem of the neural component to a simple supervised learning setup. Further, this allows us to both train on sparingly available real-world data and still leverage powerful recurrent neural network encoders. We demonstrate the effectiveness of our method by evaluating on real-world customer scenarios by synthesizing accurate programs with up to 12x speed-up compared to state-of-the-art systems.


title: On the Existence of Uniformly Most Powerful Bayesian Tests With Application to Non-Central Chi-Squared Tests

abstract: Uniformly most powerful Bayesian tests (UMPBT's) are an objective class of Bayesian hypothesis tests that can be considered the Bayesian counterpart of classical uniformly most powerful tests. Because the rejection regions of UMPBT's can be matched to the rejection regions of classical uniformly most powerful tests (UMPTs), UMPBT's provide a mechanism for calibrating Bayesian evidence thresholds, Bayes factors, classical significance levels and p-values. The purpose of this article is to expand the application of UMPBT's outside the class of exponential family models. Specifically, we introduce sufficient conditions for the existence of UMPBT's and propose a unified approach for their derivation. An important application of our methodology is the extension of UMPBT's to testing whether the non-centrality parameter of a chi-squared distribution is zero. The resulting tests have broad applicability, providing default alternative hypotheses to compute Bayes factors in, for example, Pearson's chi-squared test for goodness-of-fit, tests of independence in contingency tables, and likelihood ratio, score and Wald tests.


title: Hospital Readmission Prediction - Applying Hierarchical Sparsity Norms for Interpretable Models

abstract: Hospital readmissions have become one of the key measures of healthcare quality. Preventable readmissions have been identified as one of the primary targets for reducing costs and improving healthcare delivery. However, most data driven studies for understanding readmissions have produced black box classification and predictive models with moderate performance, which precludes them from being used effectively within the decision support systems in the hospitals. In this paper we present an application of structured sparsity-inducing norms for predicting readmission risk for patients based on their disease history and demographics. Most existing studies have focused on hospital utilization, test results, etc., to assign a readmission label to each episode of hospitalization. However, we focus on assigning a readmission risk label to a patient based on their disease history. Our emphasis is on interpreting the models to improve the understanding of the readmission problem. To achieve this, we exploit the domain induced hierarchical structure available for the disease codes which are the features for the classification algorithm. We use a tree based sparsity-inducing regularization strategy that explicitly uses the domain hierarchy. The resulting model not only outperforms standard regularization procedures but is also highly sparse and interpretable. We analyze the model and identify several significant factors that have an effect on readmission risk. Some of these factors conform to existing beliefs, e.g., impact of surgical complications and infections during hospital stay. Other factors, such as the impact of mental disorder and substance abuse on readmission, provide empirical evidence for several pre-existing but unverified hypotheses. The analysis also reveals previously undiscovered connections such as the influence of socioeconomic factors like lack of housing and malnutrition.


title: Real-Time Prediction of the Duration of Distribution System Outages

abstract: This paper addresses the problem of predicting duration of unplanned power outages, using historical outage records to train a series of neural network predictors. The initial duration prediction is made based on environmental factors, and it is updated based on incoming field reports using natural language processing to automatically analyze the text. Experiments using 15 years of outage records show good initial results and improved performance leveraging text. Case studies show that the language processing identifies phrases that point to outage causes and repair steps.


title: Different approach to the modeling of nonfree particle diffusion

abstract: A new approach to the modeling of nonfree particle diffusion is presented. The approach uses a general setup based on geometric graphs (networks of curves), which means that particle diffusion in anything from arrays of barriers and pore networks to general geometric domains can be considered and that the (free random walk) central limit theorem can be generalized to cover also the nonfree case. The latter gives rise to a continuum-limit description of the diffusive motion where the effect of partially absorbing barriers is accounted for in a natural and non-Markovian way that, in contrast to the traditional approach, quantifies the absorptivity of a barrier in terms of a dimensionless parameter in the range 0 to 1. The generalized theorem gives two general analytic expressions for the continuum-limit propagator: an infinite sum of Gaussians and an infinite sum of plane waves. These expressions entail the known method-of-images and Laplace eigenfunction expansions as special cases and show how the presence of partially absorbing barriers can lead to phenomena such as line splitting and band gap formation in the plane wave wave-number spectrum.


title: Supersymmetric Tools in Yang-Mills Theories at Strong Coupling: the Beginning of a Long Journey

abstract: Development of holomorphy-based methods in super-Yang-Mills theories started in the early 1980s and lead to a number of breakthrough results. I review some results in which I participated. The discovery of Seiberg's duality and the Seiberg-Witten solution of ${\cal N}=2$ Yang-Mills were the milestones in the long journey of which, I assume, much will be said in other talks. I will focus on the discovery (2003) of non-Abelian vortex strings with various degree of supersymmetry, supported in some four-dimensional Yang-Mills theories and some intriguing implications of this discovery. One of the recent results is the observation of a soliton string in the bulk ${\cal N}=2$ theory with the $U(2)$ gauge group and four flavors, which can become critical in a certain limit. This is the case of a "reverse holography," with a very transparent physical meaning.


title: Thermal theory of aluminum particle ignition in continuum, free-molecular, and transition heat transfer regimes

abstract: Most studies on nano- and micro- sized aluminum particle ignition have been focused on the processes occuring inside particles. In the current paper, thermal ignition of an aluminum particle in the air is simulated with different heat transfer models: continuum, free-molecular and Fuchs model. A single parabolic oxidation law is assumed in the particle size range from nano- to millimeter diameters. A particle is considered ignited when it reaches the oxide melting point. The criterion defining the limits of validity for each model is the ratio of continuum and free-molecular heat transfer rates. The dependence of ignition temperature $T_i$ on particle size is in a qualitative agreement with the experimental trends: $T_i$ can have values in the range 700--1500 K for nanoparticles due to the dominating contribution of a free-molecular heat transfer, and sharp growth of $T_i$ with the particle size in the range of 1-100 $\mu$m diameter is due to the transitional character of heat transfer. For small values of accommodation coefficient, ignition may occur in the critical ignition mode with the thermal runaway. The results suggest the importance of heat transfer and, in particular, energy accommodation in ignition of nano- and micro- sized particles.


title: The Logical Essentials of Bayesian Reasoning

abstract: This chapter offers an accessible introduction to the channel-based approach to Bayesian probability theory. This framework rests on algebraic and logical foundations, inspired by the methodologies of programming language semantics. It offers a uniform, structured and expressive language for describing Bayesian phenomena in terms of familiar programming concepts, like channel, predicate transformation and state transformation. The introduction also covers inference in Bayesian networks, which will be modelled by a suitable calculus of string diagrams.


title: Depth Pooling Based Large-scale 3D Action Recognition with Convolutional Neural Networks

abstract: This paper proposes three simple, compact yet effective representations of depth sequences, referred to respectively as Dynamic Depth Images (DDI), Dynamic Depth Normal Images (DDNI) and Dynamic Depth Motion Normal Images (DDMNI), for both isolated and continuous action recognition. These dynamic images are constructed from a segmented sequence of depth maps using hierarchical bidirectional rank pooling to effectively capture the spatial-temporal information. Specifically, DDI exploits the dynamics of postures over time and DDNI and DDMNI exploit the 3D structural information captured by depth maps. Upon the proposed representations, a ConvNet based method is developed for action recognition. The image-based representations enable us to fine-tune the existing Convolutional Neural Network (ConvNet) models trained on image data without training a large number of parameters from scratch. The proposed method achieved the state-of-art results on three large datasets, namely, the Large-scale Continuous Gesture Recognition Dataset (means Jaccard index 0.4109), the Large-scale Isolated Gesture Recognition Dataset (59.21%), and the NTU RGB+D Dataset (87.08% cross-subject and 84.22% cross-view) even though only the depth modality was used.


title: Probabilistic Contraction Analysis of Iterated Random Operators

abstract: In many branches of engineering, Banach contraction mapping theorem is employed to establish the convergence of certain deterministic algorithms. Randomized versions of these algorithms have been developed that have proved useful in data-driven problems. In a class of randomized algorithms, in each iteration, the contraction map is approximated with an operator that uses independent and identically distributed samples of certain random variables. This leads to iterated random operators acting on an initial point in a complete metric space, and it generates a Markov chain. In this paper, we develop a new stochastic dominance based proof technique, called probabilistic contraction analysis, for establishing the convergence in probability of Markov chains generated by such iterated random operators in certain limiting regime. The methods developed in this paper provides a general framework for understanding convergence of a wide variety of Monte Carlo methods in which contractive property is present. We apply the convergence result to conclude the convergence of fitted value iteration and fitted relative value iteration in continuous state and continuous action Markov decision problems as representative applications of the general framework developed here.


title: Determinants of Binomial-Related Circulant Matrices

abstract: Due to their rich algebraic structures and various applications, circulant matrices have been of interest and continuously studied. In this paper, the notions of Binomial-related matrices have been introduced. Such matrices are circulant matrices whose the first row is the coefficients of $(x+zy)^n$, where $z$ is a complex number of norm $1$ and $n$ is a positive integer. In the case where $z\in \{1,-1,i,-i\}$, the explicit formula for the determinant of such matrices are completely determined. Known results on the determinants of binomial circulant matrices can be viewed as special cases. Finally, some open problems are discussed.


title: Development of composite control-variate stratified sampling approach for efficient stochastic calculation of molecular integrals

abstract: In this work, the composite control-variate stratified sampling (CCSS) method is presented for calculation of MO integrals without transformation of AO integrals. The central idea of this approach is to obtain the 2-electron MO integrals by direct integration of 2-electron coordinates. This method does not require or use pre-computed AO integrals and the value of the MOs at any point in space is obtained directly from the linear combination of AOs. The integration over the electronic coordinates was performed using stratified sampling Monte Carlo method. This approach was implemented by dividing the integration region into a set of non-overlapping segments and performing Monte Carlo calculations on each segment. The Monte Carlo sampling points for each segment were optimized to minimize the total variance of the sample mean. Additional variance reduction of the overall calculations was achieved by introducing control-variate in the stratified sampling scheme. The composite aspect of the CCSS allows for simultaneous computation of multiple MO integrals during the stratified sampling evaluation. The main advantage of the CCSS method is that unlike rejection sampling Monte Carlo methods such as Metropolis algorithm, the stratified sampling uses all instances of the calculated functions for the evaluation of the sample mean. The CCSS method is designed to be used for large systems where AO-to-MO transformation is computationally prohibitive. Because it is based on numerical integration, the CCSS method can be applied to a wide variety of integration kernels and does not require \textit{a priori} knowledge of analytical integrals. In this work, the developed CCSS method was applied for calculation of excitonic properties in CdSe quantum dots using electron-hole explicitly correlated Hartree-Fock (eh-XCHF) and geminal-screened electron-hole interaction kernel (GSIK) methods.


title: New recursive approximations for variable-order fractional operators with applications

abstract: To broaden the range of applicability of variable-order fractional differential models, reliable numerical approaches are needed to solve the model equation. In this paper, we develop Laguerre spectral collocation methods for solving variable-order fractional initial value problems on the half line. Specifically, we derive three-term recurrence relations to efficiently calculate the variable-order fractional integrals and derivatives of the modified generalized Laguerre polynomials, which lead to the corresponding fractional differentiation matrices that will be used to construct the collocation methods. Comparison with other existing methods shows the superior accuracy of the proposed spectral collocation methods.


title: Entropic uncertainty relations under localizations on discrete quantum groups

abstract: The uncertainty principle has been established within the framework of locally compact quantum groups in recent years. This paper demonstrates that entropic uncertainty relations can be strengthened under localizations on discrete quantum groups, which is the case if the dual compact quantum group $\mathbb{G}$ is the free orthogonal quantum group $O_N^+$ with $N\geq 3$ or if $\mathbb{G}$ admits an infinite $\Lambda(p)$ set with $p>2$. On the other hand, this paper explains the reason why such phenomena do not appear when $\mathbb{G}$ is one of the connected semisimple compact Lie groups, $O_2^+$ and the quantum $SU(2)$ groups. Also, we discuss the divergence of entropic uncertainty relations together with some explicit explanations.


title: Admissible level $\mathfrak{osp}(1|2)$ minimal models and their relaxed highest weight modules

abstract: The minimal model $\mathfrak{osp}(1|2)$ vertex operator superalgebras are the simple quotients of affine vertex operator superalgebras constructed from the affine Lie super algebra $\widehat{\mathfrak{osp}}(1|2)$ at certain rational values of the level $k$. We classify all isomorphism classes of $\mathbb{Z}_2$-graded simple relaxed highest weight modules over the minimal model $\mathfrak{osp}(1|2)$ vertex operator superalgebras in both the Neveu-Schwarz and Ramond sectors. To this end, we combine free field realisations, screening operators and the theory of symmetric functions in the Jack basis to compute explicit presentations for the Zhu algebras in both the Neveu-Schwarz and Ramond sectors. Two different free field realisations are used depending on the level. For $k<-1$, the free field realisation resembles the Wakimoto free field realisation of affine $\mathfrak{sl}(2)$ and is originally due to Bershadsky and Ooguri. It involves 1 free boson (or rank 1 Heisenberg vertex algebra), one $\beta\gamma$ bosonic ghost system and one $bc$ fermionic ghost system. For $k>-1$, the argument presented here requires the bosonisation of the $\beta\gamma$ system by embedding it into an indefinite rank 2 lattice vertex algebra.


title: Variable selection using pseudo-variables

abstract: Penalized regression has become a standard tool for model building across a wide range of application domains. Common practice is to tune the amount of penalization to tradeoff bias and variance or to optimize some other measure of performance of the estimated model. An advantage of such automated model-building procedures is that their operating characteristics are well-defined, i.e., completely data-driven, and thereby they can be systematically studied. However, in many applications it is desirable to incorporate domain knowledge into the model building process; one way to do this is to characterize each model along the solution path of a penalized regression estimator in terms of an operating characteristic that is meaningful within a domain context and then to allow domain experts to choose from among these models using these operating characteristics as well as other factors not available to the estimation algorithm. We derive an estimator of the false selection rate for each model along the solution path using a novel variable addition method. The proposed estimator applies to both fixed and random designs and allows for $p \gg n$. The proposed estimator can be used to estimate a model with a pre-specified false selection rate or can be overlaid on the solution path to facilitate interactive model exploration. We characterize the asymptotic behavior of the proposed estimator in the case of a linear model under a fixed design; however, simulation experiments show that the proposed estimator provides consistently more accurate estimates of the false selection rate than competing methods across a wide range of models.


title: Lateral heterostructures formed by thermally converting n-type SnSe2 to p-type SnSe

abstract: Different two-dimensional materials, when combined together to form heterostructures, can exhibit exciting properties that do not exist in individual components. Therefore, intensive research efforts have been devoted to their fabrication and characterization. Previously, vertical and in-plane two-dimensional heterostructures have been formed by mechanical stacking and chemical vapor deposition. Here we report a new material system that can form in-plane p-n junctions by thermal conversion of n-type SnSe2 to p-type SnSe. Through scanning tunneling microscopy and density functional theory studies, we find that these two distinctively different lattices can form atomically sharp interfaces and have a type II to nearly type III band alignment. We also demonstrate that this method can be used to create micron sized in-plane p-n junctions at predefined locations. These findings pave the way for further exploration of the intriguing properties of the SnSe2-SnSe heterostructure.


title: An integration of fast alignment and maximum-likelihood methods for electron subtomogram averaging and classification

abstract: Motivation: Cellular Electron CryoTomography (CECT) is an emerging 3D imaging technique that visualizes subcellular organization of single cells at submolecular resolution and in near-native state. CECT captures large numbers of macromolecular complexes of highly diverse structures and abundances. However, the structural complexity and imaging limits complicate the systematic de novo structural recovery and recognition of these macromolecular complexes. Efficient and accurate reference-free subtomogram averaging and classification represent the most critical tasks for such analysis. Existing subtomogram alignment based methods are prone to the missing wedge effects and low signal-to-noise ratio (SNR). Moreover, existing maximum-likelihood based methods rely on integration operations, which are in principle computationally infeasible for accurate calculation. Results: Built on existing works, we propose an integrated method, Fast Alignment Maximum Likelihood method (FAML), which uses fast subtomogram alignment to sample sub-optimal rigid transformations. The transformations are then used to approximate integrals for maximum-likelihood update of subtomogram averages through expectation-maximization algorithm. Our tests on simulated and experimental subtomograms showed that, compared to our previously developed fast alignment method (FA), FAML is significantly more robust to noise and missing wedge effects with moderate increases of computation cost.Besides, FAML performs well with significantly fewer input subtomograms when the FA method fails. Therefore, FAML can serve as a key component for improved construction of initial structural models from macromolecules captured by CECT.


title: Remarks on singular Cayley graphs and vanishing elements of simple groups

abstract: Let $\Gamma$ be a finite graph and let $A(\Gamma)$ be its adjacency matrix. Then $\Gamma$ is {\it singular} if $A(\Gamma)$ is singular. The singularity of graphs is of certain interest in graph theory and algebraic combinatorics. Here we investigate this problem for Cayley graphs ${\rm Cay}(G,H)$ when $G$ is a finite group and when the connecting set $H$ is a union of conjugacy classes of $G.$ In this situation the singularity problem reduces to finding an irreducible character $\chi$ of $G$ for which $\sum_{h\in H}\,\chi(h)=0.$ At this stage we focus on the case when $H$ is a single conjugacy class $h^G$ of $G.$ Here the above equality is equivalent to $\chi(h)=0$. Much is known in this situation, with essential information coming from the block theory of representations of finite groups. An element $h\in G$ is called vanishing if $\chi(h)=0$ for some irreducible character $\chi$ of $G.$ We study vanishing elements mainly in finite simple groups and in alternating groups in particular. We suggest some approaches for constructing singular Cayley graphs.


title: Interval partition evolutions with emigration related to the Aldous diffusion

abstract: We construct a stationary Markov process corresponding to the evolution of masses and distances of subtrees along the spine from the root to a branch point in a conjectured stationary, continuum random tree-valued diffusion that was proposed by David Aldous. As a corollary this Markov process induces a recurrent extension, with Dirichlet stationary distribution, of a Wright-Fisher diffusion for which zero is an exit boundary of the coordinate processes. This extends previous work of Pal who argued a Wright-Fisher limit for the three-mass process under the conjectured Aldous diffusion until the disappearance of the branch point. In particular, the construction here yields the first stationary, Markovian projection of the conjectured diffusion. Our construction follows from that of a pair of interval partition-valued diffusions that were previously introduced by the current authors as continuum analogues of down-up chains on ordered Chinese restaurants with parameters (1/2,1/2) and (1/2,0). These two diffusions are given by an underlying Crump-Mode-Jagers branching process, respectively with or without immigration. In particular, we adapt the previous construction to build a continuum analogue of a down-up ordered Chinese restaurant process with the unusual parameters (1/2,-1/2), for which the underlying branching process has emigration.


title: Erratum: Investigating Neutron Polarizabilities through Compton Scattering on ${}^3$He

abstract: We provide updated predictions for elastic \gamma ${}^3$He cross sections and asymmetries that correct erroneous results we published in Phys. Rev. Lett. 98, 232303 (2007) and Nucl. Phys. A 819, 98 (2009).


title: A Euclidean Algorithm for Binary Cycles with Minimal Variance

abstract: The problem is considered of arranging symbols around a cycle, in such a way that distances between different instances of a same symbol be as uniformly distributed as possible. A sequence of moments is defined for cycles, similarly to the well-known praxis in statistics and including mean and variance. Mean is seen to be invariant under permutations of the cycle. In the case of a binary alphabet of symbols, a fast, constructive, sequencing algorithm is introduced, strongly resembling the celebrated Euclidean method for greatest common divisor computation, and the cycle returned is characterized in terms of symbol distances. A minimal variance condition is proved, and the proposed Euclidean algorithm is proved to satisfy it, thus being optimal. Applications to productive systems and information processing are briefly discussed.


title: Should We Adjust for the Test for Pre-trends in Difference-in-Difference Designs?

abstract: The common practice in difference-in-difference (DiD) designs is to check for parallel trends prior to treatment assignment, yet typical estimation and inference does not account for the fact that this test has occurred. I analyze the properties of the traditional DiD estimator conditional on having passed (i.e. not rejected) the test for parallel pre-trends. When the DiD design is valid and the test for pre-trends confirms it, the typical DiD estimator is unbiased, but traditional standard errors are overly conservative. Additionally, there exists an alternative unbiased estimator that is more efficient than the traditional DiD estimator under parallel trends. However, when in population there is a non-zero pre-trend but we fail to reject the hypothesis of parallel pre-trends, the DiD estimator is generally biased relative to the population DiD coefficient. Moreover, if the trend is monotone, then under reasonable assumptions the bias from conditioning exacerbates the bias relative to the true treatment effect. I propose new estimation and inference procedures that account for the test for parallel trends, and compare their performance to that of the traditional estimator in a Monte Carlo simulation.


title: Neptune's 5:2 Resonance in the Kuiper Belt

abstract: Observations of Kuiper belt objects (KBOs) in Neptune's 5:2 resonance present two puzzles: this third order resonance hosts a surprisingly large population, comparable to the prominent populations of Plutinos and Twotinos in the first order 3:2 and 2:1 resonances, respectively; secondly, their eccentricities are concentrated near $0.4$. To shed light on these puzzles, we investigate the phase space near this resonance with use of Poincar\'e sections of the circular planar restricted three body model. We find several transitions in the phase space structure with increasing eccentricity, which we explain with the properties of the resonant orbit relative to Neptune's. The resonance width is narrow for very small eccentricities, but widens dramatically for $e\gtrsim0.2$, reaching a maximum near $e\approx0.4$, where it is similar to the maximum widths of the 2:1 and 3:2 resonances. We confirm these results with N-body numerical simulations, including the effects of all four giant planets and a wide range of orbital inclinations of the KBOs. We find that the boundaries of the stable resonance zone are not strongly sensitive to inclination and remain very similar to those found with the simplified three body model, with the caveat that orbits of eccentricity above $\sim0.53$ are unstable; higher eccentricity orbits are phase-protected from destabilizing encounters with Neptune but not with Uranus. These results show that the 5:2 resonant KBOs are not more puzzling than the Plutinos and Twotinos; however, detailed understanding of the origins of eccentric, inclined resonant KBOs remains a challenge.


title: A Segmentation-aware Deep Fusion Network for Compressed Sensing MRI

abstract: Compressed sensing MRI is a classic inverse problem in the field of computational imaging, accelerating the MR imaging by measuring less k-space data. The deep neural network models provide the stronger representation ability and faster reconstruction compared with "shallow" optimization-based methods. However, in the existing deep-based CS-MRI models, the high-level semantic supervision information from massive segmentation-labels in MRI dataset is overlooked. In this paper, we proposed a segmentation-aware deep fusion network called SADFN for compressed sensing MRI. The multilayer feature aggregation (MLFA) method is introduced here to fuse all the features from different layers in the segmentation network. Then, the aggregated feature maps containing semantic information are provided to each layer in the reconstruction network with a feature fusion strategy. This guarantees the reconstruction network is aware of the different regions in the image it reconstructs, simplifying the function mapping. We prove the utility of the cross-layer and cross-task information fusion strategy by comparative study. Extensive experiments on brain segmentation benchmark MRBrainS validated that the proposed SADFN model achieves state-of-the-art accuracy in compressed sensing MRI. This paper provides a novel approach to guide the low-level visual task using the information from mid- or high-level task.


title: Constructions for optimal Ferrers diagram rank-metric codes

abstract: Optimal rank-metric codes in Ferrers diagrams can be used to construct good subspace codes. Such codes consist of matrices having zeros at certain fixed positions. This paper generalizes the known constructions for Ferrers diagram rank-metric (FDRM) codes. Via a criteria for linear maximum rank distance (MRD) codes, an explicit construction for a class of systematic MRD codes is presented, which is used to produce new optimal FDRM codes. By exploring subcodes of Gabidulin codes, if each of the rightmost $\delta-1$ columns in Ferrers diagram $\cal F$ has at least $n-r$ dots, where $r$ is taken in a range, then the conditions that an FDRM code in $\cal F$ is optimal are established. The known combining constructions for FDRM code are generalized by introducing the concept of proper combinations of Ferrers diagrams.


title: Classification of Vehicles Based on Audio Signals using Quadratic Discriminant Analysis and High Energy Feature Vectors

abstract: The focus of this paper is on classification of different vehicles using sound emanated from the vehicles. In this paper,quadratic discriminant analysis classifies audio signals of passing vehicles to bus, car, motor, and truck categories based on features such as short time energy, average zero cross rate, and pitch frequency of periodic segments of signals. Simulation results show that just by considering high energy feature vectors, better classification accuracy can be achieved due to the correspondence of low energy regions with noises of the background. To separate these elements, short time energy and average zero cross rate are used simultaneously.In our method,we have used a few features which are easy to be calculated in time domain and enable practical implementation of efficient classifier. Although, the computation complexity is low, the classification accuracy is comparable with other classification methods based on long feature vectors reported in literature for this problem.


title: Transcendental Semigroup that has Simply Connected Fatou Components

abstract: We prove that there exists a non-trivial transcendental semigroup S such that the wandering (or pre-periodic or periodic) components of Fatou set F(S) has at least a simply connected domain D.


title: A duality for labeled graphs and factorizations with applications to graph embeddings and Hurwitz enumeration

abstract: The set of factorizations of permutations in to $m$ transpositions of some symmetric group $\mathcal{S}_n$ is naturally in bijection with the set of graphs of order $n$ and size $m$ with both edges and vertices labeled. We define a notion of duality (the \emph{mind-body duality}) for factorizations and such labeled graphs and interpret it in terms of Properly Embedded Graphs, a class of graphs embedded in a bounded compact oriented surface with all the vertices lying in the boundary, and show a close connection of this duality with the Hurwitz action of the Braid Group. Connections with the theory of Cellularly Embedded Graphs are highlighted and hints of possible applications are given. In this paper we focus on developing the necessary theory, leaving specific applications and further developments for future projects.


title: Quantum Tunneling Of Electron Snake States In An Inhomogeneous Magnetic Field

abstract: In a two dimensional free electron gas (2DEG) subjected to a perpendicular spatially varying magnetic field, the classical paths of electrons are snake-like trajectories that weave along the line where the field crosses zero. But quantum mechanically this system is described by a symmetric double well potential which, for low excitations, leads to very different electron behavior. We compute the spectrum, as well as the wavefunctions, for states of definite parity in the limit of nearly degenerate states, i.e. for electrons sufficiently far from the $B_z=0$ line. Transitions between the states are shown to give rise to a tunneling current. If the well is made asymmetrical by a time-dependent parity breaking perturbation then Rabi-like oscillations between parity states occur. Resonances can be excited and used to stimulate the transfer of electrons from one side of the potential barrier to the other through quantum tunneling.


title: Rigidity of the $\operatorname{Sine}_{\beta}$ process

abstract: We show that the $\operatorname{Sine}_{\beta}$ point process, defined as the scaling limit of the Circular Beta Ensemble when the dimension goes to infinity, and generalizing the determinantal sine-kernel process, is rigid in the sense of Ghosh and Peres: the number of points in a given bounded Borel set $B$ is almost surely equal to a measurable function of the position of the points outside $B$.


title: Gene Co-expression Network analysis of Lung Squamous Cell Carcinoma data

abstract: We performed a gene co-expression analysis on Lung Squamous Cell Carcinoma data to find modules (groups) of genes that may highly impact the growth of these type of tumors. Additionally, we used cancer survival data to relate modules to prognostic significance in terms of survival time. Analysis on RNA-seq data revealed modules which are significant in gene enrichment analysis. Specifically, two genes - RFC4 and ECT2 - have been found to be significant in terms of survival time. We also performed a second gene co-expression analysis on a second dataset of microarray data, and many significant genes found in this analysis could also be found in the RNA-seq data implying that these genes might indeed play a crucial role in Lung Squamous Cancer. All the R code for the analysis can be found at: \url{https://github.com/nafizh/Gene_coexpression_analysis_lung_cancer}