title: Isoperimetry with upper mean curvature bounds and sharp stability estimates

abstract: It was proved by Almgren that among boundaries whose mean curvature is bounded from above, perimeter is uniquely minimized by balls. We obtain sharp stability estimates for Almgren's isoperimetric principle and, as an application, we deduce a sharp description of boundaries with almost constant mean curvature under a total perimeter bound which prevents bubbling.


title: Computing the generators of the truncated real radical ideal by moment matrices and SDP facial reduction

abstract: Recent breakthroughs have been made in the use of semidefinite programming and its application to real polynomial solving. For example, the real radical of a zero dimensional ideal, can be determined by such approaches as shown by Lasserre and collaborators. Some progress has been made on the determination of the real radical in positive dimension by Ma, Wang and Zhi. Such work involves the determination of maximal rank semidefinite moment matrices. Existing methods are computationally expensive and have poorer accuracy on larger examples. This paper is motivated by problems in the numerical computation of the real radical ideal in the general positive case. In this paper we give a method to compute the generators of the real radical for any given degree $d$. We combine the use of moment matrices and techniques from SDP optimization: facial reduction first developed by Borwein and Wolkowicz. In use of the semidefinite moment matrices to compute the real radical, the maximum rank property is very key, and with facial reduction, it can be guaranteed with very high accuracy. Our algorithm can be used to test the real radical membership of a given polynomial. In a special situation, we can determine the real radical ideal in the positive dimensional case.


title: On the orthogonality of the Chebyshev-Frolov lattice and applications

abstract: We deal with lattices that are generated by the Vandermonde matrices associated to the roots of Chebyshev-polynomials. If the dimension $d$ of the lattice is a power of two, i.e. $d=2^m, m \in \mathbb{N}$, the resulting lattice is an admissible lattice in the sense of Skriganov. Those are related to the Frolov cubature formulas, which recently drew attention due to their optimal convergence rates in a broad range of function spaces with mixed smoothness. We prove that the resulting lattices are orthogonal and possess a lattice representation matrix with entries not larger than $2$ (in modulus). This allows for an efficient enumeration of the Frolov cubature nodes in the $d$-cube $[-1/2,1/2]^d$ up to dimension $d=16$.


title: Minimal positive design for self-assembly of the Archimedean tilings

abstract: A challenge of molecular self-assembly is to understand how to design particles that self-assemble into a desired structure and not any of a potentially large number of undesired structures. Here we use simulation to show that a strategy of minimal positive design allows the self-assembly of networks equivalent to the 8 semiregular Archimedean tilings of the plane, structures not previously realized in simulation. This strategy consists of identifying the fewest distinct types of interparticle interaction that appear in the desired structure, and does not require enumeration of the many possible undesired structures. The resulting particles, which self-assemble into the desired networks, possess DNA-like selectivity of their interactions. Assembly of certain molecular networks may therefore require such selectivity.


title: The average singular value of a complex random matrix decreases with dimension

abstract: We obtain a recurrence relation in $d$ for the average singular value $% \alpha (d)$ of a complex valued $d\times d$\ matrix $\frac{1}{\sqrt{d}}X$ with random i.i.d., N( 0,1) entries, and use it to show that $\alpha (d)$ decreases monotonically with $d$ to the limit given by the Marchenko-Pastur distribution.\ The monotonicity of $\alpha (d)$ has been recently conjectured by Bandeira, Kennedy and Singer in their study of the Little Grothendieck problem over the unitary group $\mathcal{U}_{d}$ \cite{BKS}, a combinatorial optimization problem. The result implies sharp global estimates for $\alpha (d)$, new bounds for the expected minimum and maximum singular values, and a lower bound for the ratio of the expected maximum and the expected minimum singular value. The proof is based on a connection with the theory of Tur\'{a}n determinants of orthogonal polynomials. We also discuss some applications to the problem that originally motivated the conjecture.


title: Modelling cross-reactivity and memory in the cellular adaptive immune response to influenza infection in the host

abstract: The cellular adaptive immune response plays a key role in resolving influenza infection. Experiments where individuals are successively infected with different strains within a short timeframe provide insight into the underlying viral dynamics and the role of a cross-reactive immune response in resolving an acute infection. We construct a mathematical model of within-host influenza viral dynamics including three possible factors which determine the strength of the cross-reactive cellular adaptive immune response: the initial naive T cell number, the avidity of the interaction between T cells and the epitopes presented by infected cells, and the epitope abundance per infected cell. Our model explains the experimentally observed shortening of a second infection when cross-reactivity is present, and shows that memory in the cellular adaptive immune response is necessary to protect against a second infection.


title: On the equivalence between Kolmogorov-Smirnov and ROC curve metrics for binary classification

abstract: Binary decisions are very common in artificial intelligence. Applying a threshold on the continuous score gives the human decider the power to control the operating point to separate the two classes. The classifier,s discriminating power is measured along the continuous range of the score by the Area Under the ROC curve (AUC_ROC) in most application fields. Only finances uses the poor single point metric maximum Kolmogorov-Smirnov (KS) distance. This paper proposes the Area Under the KS curve (AUC_KS) for performance assessment and proves AUC_ROC = 0.5 + AUC_KS, as a simpler way to calculate the AUC_ROC. That is even more important for ROC averaging in ensembles of classifiers or n fold cross-validation. The proof is geometrically inspired on rotating all KS curve to make it lie on the top of the ROC chance diagonal. On the practical side, the independent variable on the abscissa on the KS curve simplifies the calculation of the AUC_ROC. On the theoretical side, this research gives insights on probabilistic interpretations of classifiers assessment and integrates the existing body of knowledge of the information theoretical ROC approach with the proposed statistical approach based on the thoroughly known KS distribution.


title: The Jackknife Estimation Method

abstract: Statistical resampling methods have become feasible for parametric estimation, hypothesis testing, and model validation now that the computer is a ubiquitous tool for statisticians. This essay focuses on the resampling technique for parametric estimation known as the Jackknife procedure. To outline the usefulness of the method and its place in the general class of statistical resampling techniques, I will quickly delineate two similar resampling methods: the bootstrap and the permutation test. I then outline the Jackknife method and show an example of its use.


title: A Comprehensive Linear Speedup Analysis for Asynchronous Stochastic Parallel Optimization from Zeroth-Order to First-Order

abstract: Asynchronous parallel optimization received substantial successes and extensive attention recently. One of core theoretical questions is how much speedup (or benefit) the asynchronous parallelization can bring us. This paper provides a comprehensive and generic analysis to study the speedup property for a broad range of asynchronous parallel stochastic algorithms from the zeroth order to the first order methods. Our result recovers or improves existing analysis on special cases, provides more insights for understanding the asynchronous parallel behaviors, and suggests a novel asynchronous parallel zeroth order method for the first time. Our experiments provide novel applications including model blending problems using the proposed asynchronous parallel zeroth order method.


title: Generalizing and Hybridizing Count-based and Neural Language Models

abstract: Language models (LMs) are statistical models that calculate probabilities over sequences of words or other discrete symbols. Currently two major paradigms for language modeling exist: count-based n-gram models, which have advantages of scalability and test-time speed, and neural LMs, which often achieve superior modeling performance. We demonstrate how both varieties of models can be unified in a single modeling framework that defines a set of probability distributions over the vocabulary of words, and then dynamically calculates mixture weights over these distributions. This formulation allows us to create novel hybrid models that combine the desirable features of count-based and neural LMs, and experiments demonstrate the advantages of these approaches.


title: Variation of transition temperatures and residual resistivity ratio in vapor-grown FeSe

abstract: The study of the iron-based superconductor FeSe has blossomed with the availability of high quality single crystals, obtained through flux/vapor-transport growth techniques below the structural transformation temperature of its tetragonal phase, T~450 C. Here, we report on the variation of sample morphology and properties due to small modifications in the growth conditions. A considerable variation of the superconducting transition temperature Tc, from 8.8 K to 3 K, which cannot be correlated with the sample composition, is observed. Instead, we point out a clear correlation between Tc and disorder, as measured by the residual resistivity ratio. Notably, the tetragonal-to-orthorhombic structural transition is also found to be quite strongly disorder dependent (Ts =72 - 90 K), and linearly correlated with Tc.


title: Effect of nonmagnetic impurities on $s_\pm$ superconductivity in the presence of incipient bands

abstract: Several Fe chalcogenide superconductors without hole pockets at the Fermi level display high temperature superconductivity, in apparent contradiction to naive spin fluctuation pairing arguments. Recently, scanning tunneling microscopy measurements have measured the influence of impurities on some of these materials, and claimed that non-magnetic impurities do not create in-gap states, leading to the conclusion that the gap must be $s_{++}$, i.e. conventional $s$ wave with no gap sign change. Here we critique this argument, and give various ways sign-changing gaps can be consistent with the absence of such bound states. In particular, we calculate the bound states for an $s_\pm$ system with a hole pocket below the Fermi level, and show that the nonmagnetic impurity bound state energy generically tracks the gap edge $E_{min}$ in the system, thereby rendering it unobservable. A failure to observe a bound state in the case of a nonmagnetic impurity can therefore not be used as an argument to exclude sign-changing pairing states.


title: Program Repair by Stepwise Correctness Enhancement

abstract: Relative correctness is the property of a program to be more-correct than another with respect to a given specification. Whereas the traditional definition of (absolute) correctness divides candidate program into two classes (correct, and incorrect), relative correctness arranges candidate programs on the richer structure of a partial ordering. In other venues we discuss the impact of relative correctness on program derivation, and on program verification. In this paper, we discuss the impact of relative correctness on program testing; specifically, we argue that when we remove a fault from a program, we ought to test the new program for relative correctness over the old program, rather than for absolute correctness. We present analytical arguments to support our position, as well as an empirical argument in the form of a small program whose faults are removed in a stepwise manner as its relative correctness rises with each fault removal until we obtain a correct program.


title: Model-based Testing of Mobile Systems -- An Empirical Study on QuizUp Android App

abstract: We present an empirical study in which model-based testing (MBT) was applied to a mobile system: the Android client of QuizUp, the largest mobile trivia game in the world. The study shows that traditional MBT approaches based on extended finite-state machines can be used to test a mobile app in an effective and efficient way. Non-trivial defects were detected on a deployed system that has millions of users and was already well tested. The duration of the overall testing effort was of three months, including the construction of the models. Maintaining a single behavioral model for the app was key in order to test it in an efficient way.


title: Using Multi-Viewpoint Contracts for Negotiation of Embedded Software Updates

abstract: In this paper we address the issue of change after deployment in safety-critical embedded system applications. Our goal is to substitute lab-based verification with in-field formal analysis to determine whether an update may be safely applied. This is challenging because it requires an automated process able to handle multiple viewpoints such as functional correctness, timing, etc. For this purpose, we propose an original methodology for contract-based negotiation of software updates. The use of contracts allows us to cleanly split the verification effort between the lab and the field. In addition, we show how to rely on existing viewpoint-specific methods for update negotiation. We illustrate our approach on a concrete example inspired by the automotive domain.


title: Monitoring Assumptions in Assume-Guarantee Contracts

abstract: Pre-deployment verification of software components with respect to behavioral specifications in the assume-guarantee form does not, in general, guarantee absence of errors at run time. This is because assumptions about the environment cannot be discharged until the environment is fixed. An intuitive approach is to complement pre-deployment verification of guarantees, up to the assumptions, with post-deployment monitoring of environment behavior to check that the assumptions are satisfied at run time. Such a monitor is typically implemented by instrumenting the application code of the component. An additional challenge for the monitoring step is that environment behaviors are typically obtained through an I/O library, which may alter the component's view of the input format. This transformation requires us to introduce a second pre-deployment verification step to ensure that alarms raised by the monitor would indeed correspond to violations of the environment assumptions. In this paper, we describe an approach for constructing monitors and verifying them against the component assumption. We also discuss limitations of instrumentation-based monitoring and potential ways to overcome it.


title: Preliminary Results Towards Contract Monitorability

abstract: This paper discusses preliminary investigations on the monitorability of contracts for web service descriptions. There are settings where servers do not guarantee statically whether they satisfy some specified contract, which forces the client (i.e., the entity interacting with the server) to perform dynamic checks. This scenario may be viewed as an instance of Runtime Verification, where a pertinent question is whether contracts can be monitored for adequately at runtime, otherwise stated as the monitorability of contracts. We consider a simple language of finitary contracts describing both clients and servers, and develop a formal framework that describes server contract monitoring. We define monitor properties that potentially contribute towards a comprehensive notion of contract monitorability and show that our simple contract language satisfies these properties.


title: The Gorenstein property for projective coordinate rings of the moduli of parabolic $\mathrm{SL}_2$-principal bundles on a smooth curve

abstract: Using combinatorial methods, we determine that a projective coordinate ring of the moduli of parabolic principal $\mathrm{SL}_2$-bundles on a marked projective curve is not Gorenstein when the genus and number of marked points are greater than $1$.


title: Open string multi-branched and Kahler potentials

abstract: We consider type II string compactifications on Calabi-Yau orientifolds with fluxes and D-branes, and analyse the F-term scalar potential that simultaneously involves closed and open string modes. In type IIA models with D6-branes such potential can be directly computed by integrating out Minkowski three-forms. The result shows a multi-branched structure along the space of lifted open string moduli, in which discrete shifts in special Lagrangian and Wilson line deformations are compensated by changes in the RR flux quanta. The same sort of discrete shift symmetries are present in the superpotential and constrain the Kahler potential. As for the latter, inclusion of open string moduli breaks the factorisation between complex structure and Kahler moduli spaces. Nevertheless, the 4d Kahler metrics display a set of interesting relations that allow to rederive the scalar potential analytically. Similar results hold for type IIB flux compactifications with D7-brane Wilson lines.


title: High Current Density 2D/3D Esaki Tunnel Diodes

abstract: The integration of two-dimensional materials such as transition metal dichalcogenides with bulk semiconductors offer interesting opportunities for 2D/3D heterojunction-based novel device structures without any constraints of lattice matching. By exploiting the favorable band alignment at the GaN/MoS2 heterojunction, an Esaki interband tunnel diode is demonstrated by transferring large area, Nb-doped, p-type MoS2 onto heavily n-doped GaN. A peak current density of 446 A/cm2 with repeatable room temperature negative differential resistance, peak to valley current ratio of 1.2, and minimal hysteresis was measured in the MoS2/GaN non-epitaxial tunnel diode. A high current density of 1 kA/cm2 was measured in the Zener mode (reverse bias) at -1 V bias. The GaN/MoS2 tunnel junction was also modeled by treating MoS2 as a bulk semiconductor, and the electrostatics at the 2D/3D interface was found to be crucial in explaining the experimentally observed device characteristics.


title: Full extremal process, cluster law and freezing for the two-dimensional discrete Gaussian Free Field

abstract: We study the local structure of the extremal process associated with the Discrete Gaussian Free Field (DGFF) in scaled-up (square-)lattice versions of bounded open planar domains subject to mild regularity conditions on the boundary. We prove that, in the scaling limit, this process tends to a Cox process decorated by independent, correlated clusters whose distribution is completely characterized. As an application, we control the scaling limit of the discrete supercritical Liouville measure, extract a Poisson-Dirichlet statistics for the limit of the Gibbs measure associated with the DGFF and establish the "freezing phenomenon" conjectured to occur in the "glassy" phase. In addition, we prove a local limit theorem for the position and value of the absolute maximum. The proofs are based on a concentric, finite-range decomposition of the DGFF and entropic-repulsion arguments for an associated random walk. Although we naturally build on our earlier work on this problem, the methods developed here are largely independent.


title: Distributed Hessian-Free Optimization for Deep Neural Network

abstract: Training deep neural network is a high dimensional and a highly non-convex optimization problem. Stochastic gradient descent (SGD) algorithm and it's variations are the current state-of-the-art solvers for this task. However, due to non-covexity nature of the problem, it was observed that SGD slows down near saddle point. Recent empirical work claim that by detecting and escaping saddle point efficiently, it's more likely to improve training performance. With this objective, we revisit Hessian-free optimization method for deep networks. We also develop its distributed variant and demonstrate superior scaling potential to SGD, which allows more efficiently utilizing larger computing resources thus enabling large models and faster time to obtain desired solution. Furthermore, unlike truncated Newton method (Marten's HF) that ignores negative curvature information by using na\"ive conjugate gradient method and Gauss-Newton Hessian approximation information - we propose a novel algorithm to explore negative curvature direction by solving the sub-problem with stabilized bi-conjugate method involving possible indefinite stochastic Hessian information. We show that these techniques accelerate the training process for both the standard MNIST dataset and also the TIMIT speech recognition problem, demonstrating robust performance with upto an order of magnitude larger batch sizes. This increased scaling potential is illustrated with near linear speed-up on upto 16 CPU nodes for a simple 4-layer network.


title: Generalized gradient expansion for inhomogeneous dynamical mean-field theory: Application to ultracold atoms in a harmonic trap

abstract: We develop a generalized gradient expansion of the inhomogeneous dynamical mean-field theory method for determining properties of ultracold atoms in a trap. This approach goes beyond the well-known local density approximation and at higher temperatures, in the normal phase, it shows why the local density approximation works so well, since the local density and generalized gradient approximations are essentially indistinguishable from each other (and from the exact solution within full inhomogeneous dynamical mean-field theory). But because the generalized gradient expansion only involves nearest-neighbor corrections, it does not work as well at low temperatures, when the systems enter into ordered phases. This is primarily due to the problem that ordered phases often satisfy some global constraints which determine the spatial ordering pattern, and the local density and generalized gradient approximations are not able to impose those kinds of constraints; they also overestimate the tendency to order. The theory is applied to phase separation of different mass fermionic mixtures represented by the Falicov-Kimball model and to determining the entropy per particle of a fermionic system represented by the Hubbard model. The generalized gradient approximation is a useful diagnostic for the accuracy of the local density approximation---when both methods agree, they are likely accurate, when they disagree, neither is likely to be correct.


title: Productive Interactions: heavy particles and non-Gaussianity

abstract: We analyze the shape and amplitude of oscillatory features in the primordial power spectrum and non-Gaussianity induced by periodic production of heavy degrees of freedom coupled to the inflaton $\phi$. We find that non-adiabatic production of particles can contribute effects which are detectable or constrainable using cosmological data even if their time-dependent masses are always heavier than the scale $\dot \phi^{1/2}$, much larger than the Hubble scale. This provides a new role for UV completion, consistent with the criteria from effective field theory for when heavy fields cannot be integrated out. This analysis is motivated in part by the structure of axion monodromy, and leads to an additional oscillatory signature in a subset of its parameter space. At the level of a quantum field theory model that we analyze in detail, the effect arises consistently with radiative stability for an interesting window of couplings up to of order $\lesssim 1$. The amplitude of the bispectrum and higher-point functions can be larger than that for Resonant Non-Gaussianity, and its signal/noise may be comparable to that of the corresponding oscillations in the power spectrum (and even somewhat larger within a controlled regime of parameters). Its shape is distinct from previously analyzed templates, but was partly motivated by the oscillatory equilateral searches performed recently by the {\it Planck} collaboration. We also make some general comments about the challenges involved in making a systematic study of primordial non-Gaussianity.


title: High-harmonic and terahertz wave spectroscopy (HATS) for aligned molecules

abstract: We present the experimental and theoretical details of our recent published letter [Phys. Rev. Lett. 115. 123002] on synchronized high-harmonic and terahertz-wave spectroscopy (HATS) from nonadiabatically aligned nitrogen molecules in dual-color laser fields. Associating the alignment-angle dependent terahertz wave generation with the synchronizing high-harmonic signal, the angular differential photoionization cross section (PICS) for molecules can be reconstructed, and the minima of the angle on PICS show great convergence between the theoretical predictions and the experimental deduced results. We also show the optimal relative phase between the dual-color laser fields for terahertz wave generation dose not change with the alignment angle at a precision of $50$ attoseconds. This all-optical method provides an alternative for investigating molecular structures and dynamics.


title: PEXSI-$\Sigma$: A Green's function embedding method for Kohn-Sham density functional theory

abstract: In this paper, we propose a new Green's function embedding method called PEXSI-$\Sigma$ for describing complex systems within the Kohn-Sham density functional theory (KSDFT) framework, after revisiting the physics literature of Green's function embedding methods from a numerical linear algebra perspective. The PEXSI-$\Sigma$ method approximates the density matrix using a set of nearly optimally chosen Green's functions evaluated at complex frequencies. For each Green's function, the complex boundary conditions are described by a self energy matrix $\Sigma$ constructed from a physical reference Green's function, which can be computed relatively easily. In the linear regime, such treatment of the boundary condition can be numerically exact. The support of the $\Sigma$ matrix is restricted to degrees of freedom near the boundary of computational domain, and can be interpreted as a frequency dependent surface potential. This makes it possible to perform KSDFT calculations with $\mathcal{O}(N^2)$ computational complexity, where $N$ is the number of atoms within the computational domain. Green's function embedding methods are also naturally compatible with atomistic Green's function methods for relaxing the atomic configuration outside the computational domain. As a proof of concept, we demonstrate the accuracy of the PEXSI-$\Sigma$ method for graphene with divacancy and dislocation dipole type of defects using the DFTB+ software package.


title: Effect of rotational-state-dependent molecular alignment on the optical dipole force

abstract: The properties of molecule-optical elements such as lenses or prisms based on the interaction of molecules with optical fields depend in a crucial way on the molecular quantum state and its alignment created by the optical field. However, in previous experimental studies, the effects of state-dependent alignment have never been included in estimates of the optical dipole force acting on the molecules while previous theoretical investigations took the state-dependent molecular alignment into account only implicitly. Herein, we consider the effects of molecular alignment explicitly and, to this end, introduce an effective polarizability which takes proper account of molecular alignment and is directly related to the alignment-dependent optical dipole force. We illustrate the significance of including molecular alignment in the optical dipole force by a trajectory study that compares previously used approximations with the present approach. The trajectory simulations were carried out for an ensemble of linear molecules subject to either propagating or standing-wave optical fields for a range of temperatures and laser intensities. The results demonstrate that the alignment-dependent effective polarizability can serve to provide correct estimates of the optical dipole force, on which a state-selection method applicable to nonpolar molecules could be based. We note that an analogous analysis of the forces acting on polar molecules subject to an inhomogeneous static electric field reveals a similarly strong dependence on molecular orientation.


title: Wormholes versus black holes: quasinormal ringing at early and late times

abstract: Recently it has been argued that the phantom thin-shell wormholes matched with the Schwarzschild space-time near the Schwarzschild radius ring like Schwarzschild black holes at early times, but differently at late times (arXiv:1602.07309). Here we consider perturbations of the wormhole which was constructed without thin-shells: the Bronnikov-Ellis wormhole supported by the phantom matter and electromagnetic field. This wormhole solution is known to be stable under specific equation of state of the phantom matter. We show that if one does not use the above thin-shell matching, the wormhole, depending on the values of its parameters, either rings as the black hole at all times or rings differently also at all times. The wormhole's spectrum, investigated here, posses a number of distinctive features. In the final part we have considered general properties of scattering around arbitrary rotating traversable wormholes. We have found that symmetric and non-symmetric (with respect to the throat) wormholes are qualitatively different in this respect: First, superradiance is allowed only if for those non-symmetric wormholes for which the asymptotic values of the rotation parameters are different on both sides from the throat. Second, the symmetric wormholes cannot mimic effectively the ringing of a black hole at a few various dominant multipoles at the same time, so that the future observations of various events should easily tell the symmetric wormhole from a black hole.


title: The challenge of engaging all students via self-paced interactive e-learning tutorials for introductory physics

abstract: As research-based self-paced e-learning tools become increasingly available, a critical issue educators encounter is implementing strategies to ensure that all students engage with them as intended. Here, we discuss the effectiveness of research-based e-learning tutorials as self-paced learning tools in large enrollment brick and mortar introductory physics courses. These interactive tutorials were developed via research in physics education and were found to be effective for a diverse group of introductory physics students in one-on-one implementation. Instructors encouraged the use of these self-paced tools in a self-paced learning environment by telling students that they would be helpful for solving the assigned homework problems and that the underlying physics principles in the tutorial problems would be similar to those in the in-class quizzes (which we call paired problems). We find that many students, who struggled in the courses in which these adaptive e-learning tutorials were assigned as a self-study tool, performed poorly on the paired problems. In contrast, a majority of student volunteers in one-on-one implementation greatly benefited from the tutorials and performed well on the paired problems. This suggests that many students enrolled in introductory physics courses did not effectively engage with the self-paced tutorials outside of class and may have only used them superficially. The findings suggest that many students in need of out-of-class remediation via self-paced learning tools may have difficulty motivating themselves and may lack the self-regulation and time-management skills to engage effectively with tools specially designed to help them learn at their own pace. We conclude by proposing a theoretical framework to help students with diverse prior preparations engage effectively with self-study tools.


title: Massively-Parallel Lossless Data Decompression

abstract: Today's exponentially increasing data volumes and the high cost of storage make compression essential for the Big Data industry. Although research has concentrated on efficient compression, fast decompression is critical for analytics queries that repeatedly read compressed data. While decompression can be parallelized somewhat by assigning each data block to a different process, break-through speed-ups require exploiting the massive parallelism of modern multi-core processors and GPUs for data decompression within a block. We propose two new techniques to increase the degree of parallelism during decompression. The first technique exploits the massive parallelism of GPU and SIMD architectures. The second sacrifices some compression efficiency to eliminate data dependencies that limit parallelism during decompression. We evaluate these techniques on the decompressor of the DEFLATE scheme, called Inflate, which is based on LZ77 compression and Huffman encoding. We achieve a 2X speed-up in a head-to-head comparison with several multi-core CPU-based libraries, while achieving a 17% energy saving with comparable compression ratios.


title: K-theory and perturbations of absolutely continuous spectra

abstract: We study the K_0 group of the commutant modulo a normed ideal of an n-tuple of commuting Hermitian operators in some of the simplest cases. In case n=1, the results, under some technical conditions are rather complete and show the key role of the absolutely continuous part when the ideal is the trace-class. For a commuting n-tuple, n>2 and the Lorentz (n, 1) ideal, we show under an absolute continuity assumption that the commutant determines a canonical direct summand in K_0. Also, certain properties involving the compact ideal, established assuming quasicentral approximate units mod the normed ideal, have weaker versions which hold assuming only finiteness of the obstruction to quasicentral approximate units.


title: Initial and Eventual Software Quality Relating to Continuous Integration in GitHub

abstract: The constant demand for new features and bug fixes are forcing software projects to shorten cycles and deliver updates ever faster, while sustaining software quality. The availability of inexpensive, virtualized, cloud-computing has helped shorten schedules, by enabling continuous integration (CI) on demand. Platforms like GitHub support CI in-the-cloud. In projects using CI, a user submitting a pull request triggers a CI step. Besides speeding up build and test, this fortuitously creates voluminous archives of build and test successes and failures. CI is a relatively new phenomenon, and these archives allow a detailed study of CI. How many problems are exposed? Where do they occur? What factors affect CI failures? Does the "initial quality" as ascertained by CI predict how many bugs will later appear ("eventual quality") in the code? In this paper, we undertake a large-scale, fine resolution study of these records, to better understand CI processes, the nature, and predictors of CI failures, and the relationship of CI failures to the eventual quality of the code. We find that: a) CI failures appear to be concentrated in a few files, just like normal bugs; b) CI failures are not very highly correlated with eventual failures; c) The use of CI in a pull request doesn't necessarily mean the code in that request is of good quality.


title: Emergent Electric Polarization by Kondo Effect in a Triangular Triple Quantum Dot

abstract: A triangular triple quantum dot gives various Kondo effects, such as the emergence of an electric polarization accompanied by a complete compensation of spin degrees of freedom. The interplay of spin and charge of electrons in quantum dots is investigated using an Anderson model with an equilateral triangular spin cluster, in which a single electron state at one site is hybridized with a conduction band in a lead. The numerical renormalization group analysis shows how a nonzero electric polarization develops with local electrons traveling in the loop and how it depends on the hybridization strength as an experimentally controllable parameter in the Kondo effect.


title: Photon-induced Background for Dilepton Searches and Measurements in pp Collisions at 13 TeV

abstract: The production of high invariant mass opposite sign lepton pairs in proton-proton collisions at the LHC is dominated by the Drell-Yan process. In addition to this photon or Z exchange mediated mechanism, gamma-gamma collisions, where photons radiated by the incoming protons collide, can produce lepton pairs. This is an important additional source of background for high mass resonant (like Z') or non-resonant (like contact interactions) searches. In this paper detailed calculations of the Drell-Yan and photon-induced cross sections in the typical acceptance of a multi-purpose LHC detector at center of mass energy 13 TeV are presented. The hint for a diphoton excess at a mass around 750 GeV, reported by the ATLAS and CMS experiments from the analysis of the 2015 data at 13 TeV, raises the possibility that such an excess could be produced through gamma-gamma collisions. A good theoretical understanding and measurements in the dilepton channels at these energies can help to elucidate this production hypothesis.


title: Revisiting Lepton Flavour Universality in B Decays

abstract: Lepton flavour universality (LFU) in B-decays is revisited by considering a class of semileptonic operators defined at a scale Lambda above the electroweak scale v. The importance of quantum effects, so far neglected in the literature, is emphasised. We construct the low-energy effective Lagrangian taking into account the running effects from Lambda down to v through the one-loop renormalization group equations (RGE) in the limit of exact electroweak symmetry and QED RGEs from v down to the 1 GeV scale. The most important quantum effects turn out to be the modification of the leptonic couplings of the vector boson Z and the generation of a purely leptonic effective Lagrangian. Large LFU breaking effects in Z and tau decays and visible lepton flavour violating (LFV) effects in the processes tau -> mu ll, tau -> mu rho, tau -> mu pi and tau -> mu eta^(') are induced.


title: Pfaffian Schur processes and last passage percolation in a half-quadrant

abstract: We study last passage percolation in a half-quadrant, which we analyze within the framework of Pfaffian Schur processes. For the model with exponential weights, we prove that the fluctuations of the last passage time to a point on the diagonal are either GSE Tracy-Widom distributed, GOE Tracy-Widom distributed, or Gaussian, depending on the size of weights along the diagonal. Away from the diagonal, the fluctuations of passage times follow the GUE Tracy-Widom distribution. We also obtain a two-dimensional crossover between the GUE, GOE and GSE distribution by studying the multipoint distribution of last passage times close to the diagonal when the size of the diagonal weights is simultaneously scaled close to the critical point. We expect that this crossover arises universally in KPZ growth models in half-space. Along the way, we introduce a method to deal with diverging correlation kernels of point processes where points collide in the scaling limit.

categories: math.PR cond-mat.stat-mech math-ph math.MP


title: Dynamical structure of magnetized dissipative accretion flow around black holes

abstract: We study the global structure of optically thin, advection dominated, magnetized accretion flow around black holes. We consider the magnetic field to be turbulent in nature and dominated by the toroidal component. With this, we obtain the complete set of accretion solutions for dissipative flows where bremsstrahlung process is regarded as the dominant cooling mechanism. We show that rotating magnetized accretion flow experiences virtual barrier around black hole due to centrifugal repulsion that can trigger the discontinuous transition of the flow variables in the form of shock waves. We examine the properties of the shock waves and find that the dynamics of the post-shock corona (PSC) is controlled by the flow parameters, namely viscosity, cooling rate and strength of the magnetic field, respectively. We separate the effective region of the parameter space for standing shock and observe that shock can form for wide range of flow parameters. We obtain the critical viscosity parameter that allows global accretion solutions including shocks. We estimate the energy dissipation at the PSC from where a part of the accreting matter can deflect as outflows and jets. We compare the maximum energy that could be extracted from the PSC and the observed radio luminosity values for several super-massive black hole sources and the observational implications of our present analysis are discussed.


title: On the work of Rodriguez Hertz on rigidity in dynamics

abstract: This paper is a survey about recent progress in measure rigidity and global rigidity of Anosov actions, and celebrates the profound contributions by Federico Rodriguez Hertz to rigidity in dynamical systems.


title: ZFIRE: The Kinematics of Star-Forming Galaxies as a Function of Environment at z~2

abstract: We perform a kinematic analysis of galaxies at $z\sim2$ in the COSMOS legacy field using near-infrared (NIR) spectroscopy from Keck/MOSFIRE as part of the ZFIRE survey. Our sample consists of 75 Ks-band selected star-forming galaxies from the ZFOURGE survey with stellar masses ranging from log(M$_{\star}$/M$_{\odot}$)$=9.0-11.0$, 28 of which are members of a known overdensity at $z=2.095$. We measure H$\alpha$ emission-line integrated velocity dispersions ($\sigma_{\rm int}$) from 50$-$230 km s$^{-1}$, consistent with other emission-line studies of $z\sim2$ field galaxies. From these data we estimate virial, stellar, and gas masses and derive correlations between these properties for cluster and field galaxies at $z\sim2$. We find evidence that baryons dominate within the central effective radius. However, we find no statistically significant differences between the cluster and the field, and conclude that the kinematics of star-forming galaxies at $z\sim2$ are not significantly different between the cluster and field environments.


title: The selection rules of high harmonic generation: the symmetries of molecules and laser fields

abstract: The selection rules of high harmonic generation (HHG) are investigated using three-dimensional time-dependent density functional theory (TDDFT). From the harmonic spectra obtained with various real molecules and different forms of laser fields, several factors that contribute to selection rules are revealed. Extending the targets to stereoscopic molecules, it is shown that the allowed harmonics are dependent on the symmetries of the projections of molecules. For laser fields, the symmetries contributing to the selection rules are discussed according to Lissajous figures and their dynamical directivities. All the phenomena are explained by the symmetry of the full time-dependent Hamiltonian under a combined transformation. We present a systematic study on the selection rules and propose an intuitive method for the judgment of allowed harmonic orders, which can be extended to more complex molecules and various forms of laser pulses.


title: On American VIX options under the generalized 3/2 and 1/2 models

abstract: In this paper, we extend the 3/2-model for VIX studied by Goard and Mazur (2013) and introduce the generalized 3/2 and 1/2 classes of volatility processes. Under these models, we study the pricing of European and American VIX options and, for the latter, we obtain an early exercise premium representation using a free-boundary approach and local time-space calculus. The optimal exercise boundary for the volatility is obtained as the unique solution to an integral equation of Volterra type. We also consider a model mixing these two classes and formulate the corresponding optimal stopping problem in terms of the observed factor process. The price of an American VIX call is then represented by an early exercise premium formula. We show the existence of a pair of optimal exercise boundaries for the factor process and characterize them as the unique solution to a system of integral equations.


title: Fast and Robust Compressive Phase Retrieval with Sparse-Graph Codes

abstract: In this paper, we tackle the compressive phase retrieval problem in the presence of noise. The noisy compressive phase retrieval problem is to recover a $K$-sparse complex signal $s \in \mathbb{C}^n$, from a set of $m$ noisy quadratic measurements: $ y_i=| a_i^H s |^2+w_i$, where $a_i^H\in\mathbb{C}^n$ is the $i$th row of the measurement matrix $A\in\mathbb{C}^{m\times n}$, and $w_i$ is the additive noise to the $i$th measurement. We consider the regime where $K=\beta n^\delta$, with constants $\beta>0$ and $\delta\in(0,1)$. We use the architecture of PhaseCode algorithm, and robustify it using two schemes: the almost-linear scheme and the sublinear scheme. We prove that with high probability, the almost-linear scheme recovers $s$ with sample complexity $\Theta(K \log(n))$ and computational complexity $\Theta(n \log(n))$, and the sublinear scheme recovers $s$ with sample complexity $\Theta(K\log^3(n))$ and computational complexity $\Theta(K\log^3(n))$. To the best of our knowledge, this is the first scheme that achieves sublinear computational complexity for compressive phase retrieval problem. Finally, we provide simulation results that support our theoretical contributions.


title: Effect of Energy Polydispersity on the Nature of Lennard-Jones Liquids

abstract: In the companion paper [T. S. Ingebrigtsen and H. Tanaka, J. Phys. Chem. B 119, 11052 (2015)] the effect of size polydispersity on the nature of Lennard-Jones (LJ) liquids, which represent most molecular liquids without hydrogen bonds, was studied. More specifically, it was shown that even highly size polydisperse LJ liquids are Roskilde-simple (RS) liquids. RS liquids are liquids with strong correlation between constant volume equilibrium fluctuations of virial and potential energy and are simpler than other types of liquids. Moreover, it was shown that size polydisperse LJ liquids have isomorphs to a good approximation. Isomorphs are curves in the phase diagram of RS liquids along which structure, dynamics, and some thermodynamic quantities are invariant in dimensionless (reduced) units. In this paper, we study the effect of energy polydispersity on the nature of LJ liquids. We show that energy polydisperse LJ liquids are RS liquids. However, a tendency of particle segregation which increases with the degree of polydispersity leads to a loss of strong virial-potential energy correlation, but is mitigated with increasing temperature and/or density. Isomorphs are a good approximation also for energy polydisperse LJ liquids, although particle-resolved quantities display a somewhat poorer scaling compared to the mean quantities along the isomorph.


title: Hierarchy, dimension, attractor and self-organization -- dynamics of mode-locked fiber lasers

abstract: Mode-locked fiber lasers are one of the most important sources of ultra-short pulses. However, A unified description for the rich variety of states and the driving forces behind the complex and diverse nonlinear behavior of mode-locked fiber lasers have yet to be developed. Here we present a comprehensive theoretical framework based upon complexity science, thereby offering a fundamentally new way of thinking about the behavior of mode-locked fiber lasers. This hierarchically structured frame work provide a model with and changeable variable dimensionality resulting in a simple and elegant view, with which numerous complex states can be described systematically. The existence of a set of new mode-locked fiber laser states is proposed for the first time. Moreover, research into the attractors' basins reveals the origin of stochasticity, hysteresis and multistability in these systems. These findings pave the way for dynamics analysis and new system designs of mode-locked fiber lasers. The paradigm will have a wide range of potential applications in diverse research fields.


title: Dynamic Control of Interference Limited Underlay D2D Network

abstract: Device-to-Device (D2D) communication appears as a key communication paradigm to realizing the vision of Internet of Things (IoT) into reality by supporting heterogeneous objects interconnection in a large scale network. These devices may be many types of objects with embedded intelligence and communication capabilities, e.g., smart phones, cars, or home appliances. The issue in in this type of communication is the interference to cellular communication caused by D2D communication. Thus, proper power control and resource allocation should be coordinated in D2D network to prevent excessive interference and drastic decrease in the throughput of the cellular system. In this paper, we consider the problem of cross-layer resource allocation in time-varying cellular wireless networks with D2D communication and incorporate average interference to cellular system as a quality-of-service constraint. Specifically, each D2D pair in the network injects packets to its queue, at rates chosen in order to maximize a global utility function, subject to network stability and interference constraints. The interference constraint enforces an arbitrarily low interference to the cellular system caused by D2D communication. We first obtain the stability region for the multiuser systems assuming that the nodes have full channel state information (CSI) of their neighbors. Then, we provide a joint flow control and scheduling scheme, which is proven to achieve a utility arbitrarily close to the maximum achievable utility. Finally, we address the consequences of practical implementation issue such as distributed scheduling by a designing algorithm, which is capable of taking advantage of diversity gain introduced by fading channels. We demonstrate the efficacy of our policies by numerical studies under various network conditions.


title: Conditions for cell size homeostasis: A stochastic hybrid systems approach

abstract: A ubiquitous feature of living cells is their growth over time followed by division into daughter cells. How isogenic cell populations maintain size homeostasis, i.e., a narrow distribution of cell size, is an intriguing fundamental problem. We model cell size using a stochastic hybrid system, where a cell grows exponentially in size (volume) over time and probabilistic division events are triggered at discrete time intervals. Moreover, whenever division events occur, size is randomly partitioned among daughter cells. We first consider a scenario, where a timer (i.e., cell-cycle clock) that measures the time since the last division event regulates both the cellular growth and division rates. Analysis reveals that such a timer-controlled system cannot achieve size homeostasis, in the sense that, the cell-to-cell size variation grows unboundedly with time. To explore biologically meaningful mechanisms for controlling size we consider two classes of regulation: a size-dependent growth rate and a size-dependent division rate. Our results show that these strategies can provide bounded intercellular variation in cell size, and exact mathematical conditions on the form of regulation needed for size homeostasis are derived. Different known forms of size control strategies, such as, the adder and the sizer are shown to be consistent with these results. Interestingly, for timer-based division mechanisms, the mean cell size depends on the noise in the cell-cycle duration but independent of errors incurred in partitioning of volume among daughter cells. In contrast, the mean cell size decreases with increasing partitioning errors for size-based division mechanisms. Finally, we discuss how organisms ranging from bacteria to mammalian cells have adopted different control approaches for maintaining size homeostasis.


title: On injective dimension of F-finite F-modules and holonomic D-modules

abstract: We investigate injective dimension of $F$-finite $F$-modules in characteristic $p$ and holonomic $D$-modules in characteristic 0. One of our main results is the following. If, either $R$ is a regular ring of finite type over an infinite field of characteristic $p>0$ and $M$ is an $F_R$-finite $F_R$-module, or $R=k[x_1,\dots,x_n]$ where $k$ is a field of characteristic 0 and $M$ is a holonomic $D(R,k)$-module, then the injective dimension of $M$ is the same as the dimension of its support.


title: Study for (anti)hypertriton and light (anti)nuclei productions in high energy collisions at 200 GeV

abstract: We use the parton and hadron cascade (PACIAE) model and the dynamically constrained phase-space coalescence (DCPC) model to investigate the productions of (anti)hypertriton and light (anti)nuclei generated by 0-10% centrality $^{12}$C+$^{12}$C, $^{24}$Mg+$^{24}$Mg, $^{40}$Ca+$^{40}$Ca, and $^{64}$Cu+$^{64}$Cu collisions at $\sqrt{S_{NN}}$ = 200 GeV with $|y|$ $<$ 1.5 and $p_{T}$ $<$ 5. We study the yield ratios of the antiparticle to particle and the rapidity distributions of the different (anti)nuclei. We find that the amounts of antimatter produced are significantly lower than that of the corresponding particles, the results of theoretical models are well consistent with the PHOBOS data. The yield ratios of the particle to antiparticle in different transverse momentum regions are also given, and we find that the ratios are raise with the increase of the transverse momentum.


title: Dictionary Learning for Robotic Grasp Recognition and Detection

abstract: The ability to grasp ordinary and potentially never-seen objects is an important feature in both domestic and industrial robotics. For a system to accomplish this, it must autonomously identify grasping locations by using information from various sensors, such as Microsoft Kinect 3D camera. Despite numerous progress, significant work still remains to be done in this field. To this effect, we propose a dictionary learning and sparse representation (DLSR) framework for representing RGBD images from 3D sensors in the context of determining such good grasping locations. In contrast to previously proposed approaches that relied on sophisticated regularization or very large datasets, the derived perception system has a fast training phase and can work with small datasets. It is also theoretically founded for dealing with masked-out entries, which are common with 3D sensors. We contribute by presenting a comparative study of several DLSR approach combinations for recognizing and detecting grasp candidates on the standard Cornell dataset. Importantly, experimental results show a performance improvement of 1.69% in detection and 3.16% in recognition over current state-of-the-art convolutional neural network (CNN). Even though nowadays most popular vision-based approach is CNN, this suggests that DLSR is also a viable alternative with interesting advantages that CNN has not.


title: Geometric embedding properties of Bestvina-Brady subgroups

abstract: We compute the relative divergence and the subgroup distortion of Bestvina-Brady subgroups. We also show that for each integer $n\geq 3$, there is a free subgroup of rank $n$ of some right-angled Artin group whose inclusion is not a quasi-isometric embedding. This result answers the question of Carr about the minimum rank $n$ such that some right-angled Artin group has a free subgroup of rank $n$ whose inclusion is not a quasi-isometric embedding. It is well-known that a right-angled Artin group $A_\Gamma$ is the fundamental group of a graph manifold whenever the defining graph $\Gamma$ is a tree. We show that the Bestvina-Brady subgroup $H_\Gamma$ in this case is a horizontal surface subgroup.


title: Multi-pretrained Deep Neural Network

abstract: Pretraining is widely used in deep neutral network and one of the most famous pretraining models is Deep Belief Network (DBN). The optimization formulas are different during the pretraining process for different pretraining models. In this paper, we pretrained deep neutral network by different pretraining models and hence investigated the difference between DBN and Stacked Denoising Autoencoder (SDA) when used as pretraining model. The experimental results show that DBN get a better initial model. However the model converges to a relatively worse model after the finetuning process. Yet after pretrained by SDA for the second time the model converges to a better model if finetuned.


title: Parallel Triangular Solvers on GPU

abstract: In this paper, we investigate GPU based parallel triangular solvers systematically. The parallel triangular solvers are fundamental to incomplete LU factorization family preconditioners and algebraic multigrid solvers. We develop a new matrix format suitable for GPU devices. Parallel lower triangular solvers and upper triangular solvers are developed for this new data structure. With these solvers, ILU preconditioners and domain decomposition preconditioners are developed. Numerical results show that we can speed triangular solvers around seven times faster.


title: Homomorphisms from Specht Modules to Signed Young Permutation Modules

abstract: We construct a class $\Theta_{\mathscr{R}}$ of homomorphisms from a Specht module $S_{\mathbb{Z}}^{\lambda}$ to a signed permutation module $M_{\mathbb{Z}}(\alpha|\beta)$ which generalises James's construction of homomorphisms whose codomain is a Young permutation module. We show that any $\phi \in \operatorname{Hom}_{{\mathbb{Z}}\mathfrak{S}_{n}}\big(S_{\mathbb{Z}}^\lambda, M_{\mathbb{Z}}(\alpha|\beta)\big)$ lies in the $\mathbb{Q}$-span of $\Theta_{\text{sstd}}$, a subset of $\Theta_{\mathscr{R}}$ corresponding to semistandard $\lambda$-tableaux of type $(\alpha|\beta)$. We also study the conditions for which $\Theta^{\mathbb{F}}_{\mathrm{sstd}}$ - a subset of $\operatorname{Hom}_{\mathbb{F}\mathfrak{S}_{n}}\big(S_{\mathbb{F}}^\lambda,M_{\mathbb{F}}(\alpha|\beta)\big)$ induced by $\Theta_{\mathrm{sstd}}$ - is linearly independent, and show that it is a basis for $\operatorname{Hom}_{\mathbb{F}\mathfrak{S}_{n}}\big(S_{\mathbb{F}}^\lambda,M_{\mathbb{F}}(\alpha|\beta)\big)$ when $\mathbb{F}\mathfrak{S}_{n}$ is semisimple.


title: On stationary solutions to the vacuum Einstein field equations

abstract: We prove that any 4-dimensional geodesically complete spacetime with a timelike Killing field satisfying the vacuum Einstein field equation $Ric(g_{M})=\lambda g_{M}$ with nonnegative cosmological constant $\lambda\geq 0$ is flat. When dim $\geq 5$, if the spacetime is assumed to be static additionally, we prove that its universal cover splits isometrically as a product of a Ricci flat Riemannian manifold and a real line.


title: Distribution of residence times as a marker to distinguish different pathways for quantum transport

abstract: Electron transport through a nanoscale system is an inherently stochastic quantum mechanical process. Electric current is a time series of electron tunnelling events separated by random intervals. Thermal and quantum noise are two sources of this randomness. In this paper, we used the quantum master equation to consider the following questions: (i) Given that an electron has tunnelled into the electronically unoccupied system from the source electrode at some particular time, how long is it until an electron tunnels out to the drain electrode to leave the system electronically unoccupied, where there were no intermediate tunnelling events ("the" tunnelling path)? (ii) Given that an electron has tunnelled into the unoccupied system from the source electrode at some particular time, how long is it until an electron tunnels out to the drain electrode to leave the system electronically unoccupied, where there were no intermediate tunnelling events ("an" tunnelling path)? (iii) What are the distributions of these times? We show that electron correlations suppress the difference between "the" and "an" electron tunnelling paths.


title: Development of Krylov and AMG linear solvers for large-scale sparse matrices on GPUs

abstract: This research introduce our work on developing Krylov subspace and AMG solvers on NVIDIA GPUs. As SpMV is a crucial part for these iterative methods, SpMV algorithms for single GPU and multiple GPUs are implemented. A HEC matrix format and a communication mechanism are established. And also, a set of specific algorithms for solving preconditioned systems in parallel environments are designed, including ILU(k), RAS and parallel triangular solvers. Based on these work, several Krylov solvers and AMG solvers are developed. According to numerical experiments, favorable acceleration performance is acquired from our Krylov solver and AMG solver under various parameter conditions.


title: Forecasting wind power - Modeling periodic and non-linear effects under conditional heteroscedasticity

abstract: In this article we present an approach that enables joint wind speed and wind power forecasts for a wind park. We combine a multivariate seasonal time varying threshold autoregressive moving average (TVARMA) model with a power threshold generalized autoregressive conditional heteroscedastic (power-TGARCH) model. The modeling framework incorporates diurnal and annual periodicity modeling by periodic B-splines, conditional heteroscedasticity and a complex autoregressive structure with non-linear impacts. In contrast to usually time-consuming estimation approaches as likelihood estimation, we apply a high-dimensional shrinkage technique. We utilize an iteratively re-weighted least absolute shrinkage and selection operator (lasso) technique. It allows for conditional heteroscedasticity, provides fast computing times and guarantees a parsimonious and regularized specification, even though the parameter space may be vast. We are able to show that our approach provides accurate forecasts of wind power at a turbine-specific level for forecasting horizons of up to 48 h (short- to medium-term forecasts).


title: Modelling discrete valued cross sectional time series with observation driven models

abstract: This paper develops computationally feasible methods for estimating random effects models in the context of regression modelling of multiple independent time series of discrete valued counts in which there is serial dependence. Given covariates, random effects and process history, the observed responses at each time in each series are independent and have an exponential family distribution. We develop maximum likelihood estimation of the mixed effects model using an observation driven generalized linear autoregressive moving average specification for the serial dependence in each series. The paper presents an easily implementable approach which uses existing single time series methods to handle the serial dependence structure in combination with adaptive Gaussian quadrature to approximate the integrals over the regression random effects required for the likelihood and its derivatives. The models and methods presented allow extension of existing mixed model procedures for count data by incorporating serial dependence which can differ in form and strength across the individual series. The structure of the model has some similarities to longitudinal data transition models with random effects. However, in contrast to that setting, where there are many cases and few to moderate observations per case, the time series setting has many observations per series and a few to moderate number of cross sectional time series. The method is illustrated on time series of binary responses to musical features obtained from a panel of listeners.


title: Large-scale Reservoir Simulations on IBM Blue Gene/Q

abstract: This paper presents our work on simulation of large-scale reservoir models on IBM Blue Gene/Q and studying the scalability of our parallel reservoir simulators. An in-house black oil simulator has been implemented. It uses MPI for communication and is capable of simulating reservoir models with hundreds of millions of grid cells. Benchmarks show that our parallel simulator are thousands of times faster than sequential simulators that designed for workstations and personal computers, and the simulator has excellent scalability.


title: Smoothing of semistable Fano varieties

abstract: Given any field $k$ (not necessarily perfect), we study the smoothing of a semistable Fano variety over $k$. In characteristic 0, the reduced semistable Fano degenerate fibers of Mori fibrations are classified. In positive characteristic, under a suitable $W_2$ lifting assumption, we prove that a semistable Fano variety always appears as a degenerate fiber in a semistable family if it has a global log structure (in the sense of Fontaine-Illusie-Kato) of semistable type. A similar smoothing result over a mixed characteristic base is also obtained.


title: Similarity transformation for equilibrium flows, including effects of blowing and suction

abstract: A similarity transformation for the mean velocity profiles is obtained in sink flow turbulent boundary layers (TBL), including effects of blowing and suction. It is based on symmetry analysis which transforms the governing partial differential equations (for mean mass and momentum) into an ordinary differential equation and yields a new result including an exact, linear relation between the mean normal ($V$) and streamwise ($U$) velocities. A characteristic length is further introduced which, under a first order expansion in wall blowing/suction velocity, leads to the similarity transformation for $U$. This transformation is shown to be a group invariant under a generalized symmetry analysis and maps different $U$ profiles under different blowing/suction conditions into a (universal) profile under no blowing/suction. Its inverse transformation enables predictions of all mean quantities in the mean mass and momentum equations - $U$, $V$ and the Reynolds shear stress - in good agreement with direct numerical simulation (DNS) data.


title: JIMWLK evolution of the odderon

abstract: We study the effects of a parity-odd "odderon" correlation in JIMWLK renormalization group evolution at high energy. Firstly we show that in the eikonal picture where the scattering is described by Wilson lines, one obtains a strict mathematical upper limit for the magnitude of the odderon amplitude compared to the parity even pomeron one. This limit increases with N_c, approaching infinity in the infinite N_c limit. We use a systematic extension of the Gaussian approximation including both 2- and 3-point correlations which enables us to close the system of equations even at finite N_c. In the large-N_c limit we recover an evolution equation derived earlier. By solving this equation numerically we confirm that the odderon amplitude decreases faster in the nonlinear case than in the linear BFKL limit. We also point out that, in the 3-point truncation at finite N_c, the presence of an odderon component introduces azimuthal angular correlations ~ cos(n phi) at all n.


title: Harbourne, Schenck and Seceleanu's Conjecture

abstract: In [HSS], Conjecture 5.5.2, Harbourne, Schenck and Seceleanu conjectured that, for $r=6$ and all $r\ge 8$, the artinian ideal $I=(\ell _1^2,\dots ,l_{r+1}^2)\subset K[x_1, \dots ,x_r]$ generated by the square of $r+1$ general linear forms $\ell _{i}$ fails the Weak Lefschetz property. This paper is entirely devoted to prove this Conjecture. It is worthwhile to point out that half of the Conjecture - namely, the case when the number of variables $r$ is even - was already proved in [mmn], Theorem 6.1.


abstract: By exploiting a link between Bergman orthogonal polynomials and the Grunsky matrix, probably first observed by K{\"u}hnau in 1985, we improve some recent results on strong asymptotics of Bergman polynomials outside the domain G of orthogonality, and on entries of the Bergman shift operator. In our proofs we suggest a new matrix approach involving the Grunsky matrix, and use well-established results in the literature relating properties of the Grunsky matrix to the regularity of the boundary of G, and the associated conformal maps. For quasiconformal boundaries, this approach allows for new insights for Bergman polynomials.


title: Octant Degeneracy, Quadrant of leptonic CPV phase at Long Baseline Neutrino Experiments and Baryogenesis

abstract: In a recent work by us, we have studied, how CP violation discovery potential can be improved at long baseline neutrino experiments (LBNE/DUNE), by combining with its ND (near detector) and reactor experiments. In this work, we discuss how this study can be further analysed to resolve entanglement of the quadrant of leptonic CPV phase and Octant of atmospheric mixing angle $ \theta_{23} $, at LBNEs. The study is done for both NH (Normal hierarchy) and IH (Inverted hierarchy), HO (Higher Octant) and LO (Lower Octant). We show how baryogenesis can enhance the effect of resolving this entanglement, and how possible values of the leptonic CP-violating phase $ \delta_{CP} $ can be predicted in this context. With respect to the latest global fit data of neutrino mixing angles, we predict the values of $ \delta_{CP} $ for different cases. In this context we present favoured values of $ \delta_{CP} $ ($ \delta_{CP} $ range at $ \geq $ 2$ \sigma $ ) constrained by the latest updated BAU range and also confront our predictions of $ \delta_{CP} $ with an up-to-date global analysis of neutrino oscillation data. We find that some region of the favoured $ \delta_{CP} $ parameter space lies within the best fit values around $ \delta_{CP} \simeq 1.3\pi-1.4 \pi $. A detailed analytic and numerical study of baryogenesis through leptogenesis is performed in this framework in a model independent way.


title: Global Entropy Solutions to Weakly Nonlinear Gas Dynamics

abstract: Entropy weak solutions with bounded periodic initial data are considered for the system of weakly nonlinear gas dynamics. Through a modified Glimm scheme, an approximate solution sequence is constructed, and then a priori estimates are provided with the methods of approximate characteristics and approximate conservation laws, which gives not only the existence and uniqueness but also the uniform total variation bounds for the entropy solutions.


title: Numerical Simulation of Multi-phase Flow in Porous Media on Parallel Computers

abstract: A parallel reservoir simulator has been developed, which is designed for large-scale black oil simulations. It handles three phases, including water, oil and gas, and three components, including water, oil and gas. This simulator can calculate traditional reservoir models and naturally fractured models. Various well operations are supported, such as water flooding, gas flooding and polymer flooding. The operation constraints can be fixed bottom-hole pressure, a fixed fluid rate, and combinations of them. The simulator is based on our in-house platform, which provides grids, cell-centred data, linear solvers, preconditioners and well modeling. The simulator and the platform use MPI for communications among computation nodes. Our simulator is capable of simulating giant reservoir models with hundreds of millions of grid cells. Numerical simulations show that our simulator matches with commercial simulators and it has excellent scalability.


title: $\gamma$-ray emission signals in the massive graviton mediated dark matter model

abstract: Dark matter may interact with Standard Model (SM) particle through the exchange of a massive spin-2 graviton producing signals that can be detected. In this work we examine the $\gamma$-ray emission signals, including the line emission and the continuous spectrum component in such a massive graviton-mediated dark matter model. The constraints of LHC data, dark matter relic density as well as the dark matter indirect detection data have been applied to narrow down the parameter space. We focus on the vector dark matter model which could produce detectable $\gamma$-ray line signal. It is found that the $\gamma$-ray line data is effective on constraining the model parameters and the ongoing and upcoming space or ground-based $\gamma$-ray experiments can constrain the model further. As for the continuous $\gamma$-ray emission, the total effective annihilation cross section is $\sim 10^{-26} {\rm cm^{3}~s^{-1}}$ except at the region where dark matter mass is around the graviton mass or half of it, which is consistent with current observational data and will be reliably probed by the upcoming CTA.


title: A characterization of almost alternating knot

abstract: Generalizing Howie and Greene's characterization of alternating knots, we give a topological characterization of almost alternating knots.


title: On Landau-Zener transitions for dephasing Lindbladians

abstract: We consider a driven open system whose evolution is described by a Lindbladian. The Lindbladian is assumed to be dephasing and its Hamiltonian part to be given by the Landau-Zener Hamiltonian. We derive a formula for the transition probability which, unlike previous results, extends the Landau-Zener formula to open systems.


title: Pomeron inspired Neutrino-nucleon Cross Section in Ultra High Energy

abstract: Studies on neutrino-nucleon ($\nu N$) cross sections at different energy scales have regained interest due to increasing importance of precision measurements, as they are needed as an ingredient in all neutrino experiments. In this paper we have calculated both charged current (CC) and neutral current (NC) $\nu$N scattering cross sections at Ultra High Energy (UHE) regime in the neutrino energy ($E_{\nu}$) region i.e. $10^{9} GeV \le E_{\nu} \le 10^{12}$ GeV using QCD inspired double asymptotic limit fit of electron-proton structure function $F_{2}^{ep}$ to low $\mathit{x}$ HERA data. The form $F_{2}^{ep} \sim x^{-\lambda(Q^{2})}$ used in our analysis, can be conjectured like a dynamic pomeron (DP)-type behaviour. We also find an analytic form of the total cross sections, $\sigma_{CC}^{\nu N}$ and $\sigma_{NC}^{\nu N}$ which appear to be of hard-pomeron exchange types. A comparative analysis of our results with those available in literature is also done. An improved understanding of $\nu N$ interactions at UHE are essentially important for future oscillation experiments. Future measurements will support/confront our predictions. \textbf{Keywords}: Neutrino cross section, Ultra High Energy, QCD, Double Asymptotic limit, dynamic pomeron, hard-pomeron.


title: Mining Software Components from Object-Oriented APIs

abstract: Object-oriented Application Programing Interfaces (APIs) support software reuse by providing pre-implemented functionalities. Due to the huge number of included classes, reusing and understanding large APIs is a complex task. Otherwise, software components are admitted to be more reusable and understandable entities than object-oriented ones. Thus, in this paper, we propose an approach for reengineering object-oriented APIs into component-based ones. We mine components as a group of classes based on the frequency they are used together and their ability to form a quality-centric component. To validate our approach, we experimented on 100 Java applications that used Android APIs.


title: Creation of two-photon states via interactions between Rydberg atoms during light storage

abstract: We propose a new method to create two-photon states in a controllable way using interaction between the Rydberg atoms during the storage and retrieval of slow light. A distinctive feature of the suggested procedure is that the slow light is stored into a superposition of two atomic coherences under conditions of electromagnetically induced transparency (EIT). Interaction between the atoms during the storage period creates entangled pairs of atoms in a superposition state that is orthogonal to the initially stored state. Restoring the slow light from this new atomic state one can produce a two photon state with a second-order correlation function determined by the atom-atom interaction and the storage time. Therefore the measurement of the restored light allows one to probe the atom-atom coupling by optical means with a sensitivity that can be increased by extending the storage time. As a realization of this idea we consider a many-body Ramsey-type technique which involves pi/2 pulses creating a superposition of Rydberg states at the beginning and the end of the storage period. In that case the regenerated light is due to the resonance dipole-dipole interaction between the atoms in the Rydberg states.


title: Presentations for rook partition monoids and algebras and their singular ideals

abstract: We obtain several presentations by generators and relations for the rook partition monoids and algebras, as well as their singular ideals. Among other results, we also calculate the minimal sizes of generating sets (some of our presentations use such minimal-size generating sets), and show that the singular part of the rook partition monoid is generated by its idempotents.


title: Percolation model with an additional source of disorder

abstract: The ranges of transmission of the mobiles in a Mobile Ad-hoc Network are not uniform in reality. They are affected by the temperature fluctuation in air, obstruction due to the solid objects, even the humidity difference in the environment, etc. How the varying range of transmission of the individual active elements affects the global connectivity in the network may be an important practical question to ask. Here a new model of percolation phenomena, with an additional source of disorder, has been introduced for a theoretical understanding of this problem. As in ordinary percolation, sites of a square lattice are occupied randomly with the probability $p$. Each occupied site is then assigned a circular disc of random value $R$ for its radius. A bond is defined to be occupied if and only if the radii $R_1$ and $R_2$ of the discs centered at the ends satisfy certain pre-defined condition. In a very general formulation, one divides the $R_1 - R_2$ plane into two regions by an arbitrary closed curve. One defines that a point within one region represents an occupied bond, otherwise it is a vacant bond. Study of three different rules under this general formulation, indicates that the percolation threshold is always larger and varies continuously. This threshold has two limiting values, one is $p_c$(sq), the percolation threshold for the ordinary site percolation on the square lattice and the other being unity. The variation of the thresholds are characterized by exponents, which are not known in the literature. In a special case, all lattice sites are occupied by discs of random radii $R \in \{0,R_0\}$ and a percolation transition is observed with $R_0$ as the control variable, similar to the site occupation probability.


title: Properties of microlensing events by wide separation planets with a moon

abstract: We investigate the properties of microlensing events caused by planetary systems where planets with a moon are widely separated from their host stars. From this investigation, we find that the moon feature generally appears as an very short-duration perturbation on the smooth asymmetric light curve of the lensing event induced by the wide separation planet; thus it can be easily discriminated from the planet feature responsible for the overall asymmetric light curve. For typical Galactic lensing events with the Einstein radius of $\sim 2$ AU, the asymmetry of the light curves due to bound planets can be noticed up to $\sim 20$AU. We also find that the perturbations of the wide planetary systems become dominated by the moon as the projected star-planet separation increases, and eventually the light curves of events produced by the systems appear as the single lensing light curve of the planet itself with a very short-duration perturbation induced by the moon, which is a representative light curve of the event induced by a star and a planet, except on the Einstein timescale of the planet. We also study the effect of a finite source star on the moon feature in the wide planetary lensing events. From this study, we find that when the lunar caustic is sufficiently separated from the planetary caustic, the lower limit on the ratio of the lunar caustic size to the source radius causing a $\geq 5\%$ lunar deviation depends mostly on the projected planet-moon separation regardless of the moon/star mass ratio, and it decreases as the planet-moon separation becomes smaller or larger than the planetary Einstein radius.


title: Zonal-mean circulation response to reduced air-sea momentum roughness

abstract: The impact of uncertainties in surface layer physics on the atmospheric general circulation is comparatively unexplored. Here the sensitivity of the zonal-mean circulation to reduced air-sea momentum roughness ($Z_{0m}$) at low flow speed is investigated with the Community Atmosphere Model (CAM3). In an aquaplanet framework with prescribed sea surface temperatures, the response to reduced $Z_{0m}$ resembles the La Ni$\tilde{\text{n}}$a minus El Ni$\tilde{\text{n}}$o response to El Ni$\tilde{\text{n}}$o Southern Oscillation variability with: i) a poleward shift of the mid-latitude westerlies extending all the way to the surface; ii) a weak poleward shift of the subtropical descent region; and iii) a weakening of the Hadley circulation, which is generally also accompanied by a poleward shift of the inter-tropical convergence zone (ITCZ) and the tropical surface easterlies. Mechanism-denial experiments show this response to be initiated by the reduction of tropical latent and sensible heat fluxes, effected by reducing $Z_{0m}$. The circulation response is elucidated by considering the effect of the tropical energy fluxes on the Hadley circulation strength, the upper tropospheric critical layer latitudes, and the lower-tropospheric baroclinic eddy forcing. The ITCZ shift is understood via moist static energy budget analysis in the tropics. The circulation response to reduced $Z_{0m}$ carries over to more complex setups with seasonal cycle, full complexity of atmosphere-ice-land-ocean interaction, and a slab ocean lower boundary condition. Hence, relatively small changes in the surface parameterization parameters can lead to a significant circulation response.


title: Circulator based on spoof surface plasmon polaritons

abstract: Circulators based on spoof surface plasmon polaritons are designed and analyzed. In the letter, we use blade structure to realize the propagation of SSPPs wave and a matching transition is used to feed energy from coplanar waveguide to the SSPPs. And the circulator shows good nonreciprocal transmission characteristics. The simulation results indicate that in the frequency band from 5 to 6.6 GHz, the isolation degree and return loss basically reaches 15dB and the insertion loss is less than 0.5dB. Moreover, the use of confinement electromagnetic waves can decrease the size of the ferrite and show a broadband characteristic.


title: Octave-spanning supercontinuum generation in a silicon-rich nitride waveguide

abstract: We experimentally show octave-spanning supercontinuum generation in a non-stoichiometric silicon-rich nitride waveguide when pumped by femtosecond pulses from an erbium fiber laser. The pulse energy and bandwidth are comparable to results achieved in stoichiometric silicon nitride waveguides, but our material platform is simpler to manufacture. We also observe wave-breaking supercontinuum generation by using orthogonal pumping in the same waveguide. Additional analysis reveals that the waveguide height is a powerful tuning parameter for generating mid-infrared dispersive waves while keeping the pump in the telecom band.


title: Fixed Point Algebras for Easy Quantum Groups

abstract: Compact matrix quantum groups act naturally on Cuntz algebras. The first author isolated certain conditions under which the fixed point algebras under this action are Kirchberg algebras. Hence they are completely determined by their $K$-groups. Building on prior work by the second author, we prove that free easy quantum groups satisfy these conditions and we compute the $K$-groups of their fixed point algebras in a general form. We then turn to examples such as the quantum permutation group $S_n^+$, the free orthogonal quantum group $O_n^+$ and the quantum reflection groups $H_n^{s+}$. Our fixed point-algebra construction provides concrete examples of free actions of free orthogonal easy quantum groups, which are related to Hopf-Galois extensions.


title: Particle-number conserving analysis of the high-spin structure of $^{159}$Ho

abstract: The high-spin rotational bands in odd-$Z$ nuclei $^{159}$Ho ($Z=67$) are investigated using the cranked shell model with the pairing correlations treated by a particle-number conserving method, in which the blocking effects are taken into account exactly. The experimental moments of inertia and alignments and their variations with the rotational frequency $\hbar\omega$ are reproduced very well by the calculations. The splitting between the signature partners of the yrast band $7/2^-[523]$ is discussed and the splitting of the excited band $7/2^+[404]$ above $\hbar\omega \sim 0.30$~MeV is predicted due to the level crossing with $1/2^+[411]$. The calculated $B(E2)$ transition probabilities are also suggested for future experiments.


title: Geometric calibration of the SND detector electromagnetic calorimeter

abstract: This paper presents the design, implementation and validation of the software alignment procedure used to perform geometric calibration of the electromagnetic calorimeter with respect to the tracking system of the SND detector which is taking data at the VEPP-2000 e^{+}e^{-}collider (BINP, Novosibirsk). This procedure is based on the mathematical model describing the relative calorimeter position. The parameter values are determined by minimizing a \chi^{2} function using the difference between particle directions reconstructed in these two subdetectors for e^{+}e^{-}\rightarrow e^{+}e^{-} scattering events. The results of the calibration applied to data and MC simulation fit the model well and give an improvement in particle reconstruction. They are used in data reconstruction and MC simulation.


title: Parametrically tunable soliton-induced resonant radiation by three-wave mixing

abstract: We show that a temporal soliton can induce resonant radiation by three-wave mixing nonlinearities. This constitutes a new class of resonant radiation whose spectral positions are parametrically tunable. The experimental verification is done in a periodically poled lithium niobate crystal, where a femtosecond near-IR soliton is excited and resonant radiation waves are observed exactly at the calculated soliton phase-matching wavelengths via the sum- and difference-frequency generation nonlinearities. This extends the supercontinuum bandwidth well into the mid-IR to span 550-5000 nm and the mid-IR edge is parametrically tunable over 1000 nm by changing the three-wave mixing phase-matching condition. The results are important for bright and broadband supercontinuum generation and for frequency comb generation in quadratic nonlinear microresonators.


title: Investigation of the two-quasiparticle bands in the doubly-odd nucleus $^{166}$Ta using a particle-number conserving cranked shell model

abstract: The high-spin rotational properties of two-quasiparticle bands in the doubly-odd ${}^{166}$Ta are analyzed using the cranked shell model with pairing correlations treated by a particle-number conserving method, in which the blocking effects are taken into account exactly. The experimental moments of inertia and alignments and their variations with the rotational frequency $\hbar\omega$ are reproduced very well by the particle-number conserving calculations, which provides a reliable support to the configuration assignments in previous works for these bands. The backbendings in these two-quasiparticle bands are analyzed by the calculated occupation probabilities and the contributions of each orbital to the total angular momentum alignments. The moments of inertia and alignments for the Gallagher-Moszkowski partners of these observed two-quasiparticle rotational bands are also predicted.


title: Inclusive Four-jet Production at 7 and 13 TeV: Azimuthal Profile in Multi-Regge Kinematics

abstract: Recently, new observables in LHC inclusive events with three tagged jets were proposed. Here, we extend that proposal to events with four tagged jets. The events are characterised by one jet in the forward direction, one in the backward direction with a large rapidity distance $Y$ from the first one and two more jets tagged in more central regions of the detector. In our setup, non-tagged associated mini-jet multiplicity is present and needs to be accounted for by the inclusion of BFKL gluon Green functions. The projection of the cross section on azimuthal-angle components opens up the opportunity for defining new ratios of correlation functions of the azimuthal angle differences among the tagged jets that can be used as probes of the BFKL dynamics.


title: Ensemble-Compression: A New Method for Parallel Training of Deep Neural Networks

abstract: Parallelization framework has become a necessity to speed up the training of deep neural networks (DNN) recently. Such framework typically employs the Model Average approach, denoted as MA-DNN, in which parallel workers conduct respective training based on their own local data while the parameters of local models are periodically communicated and averaged to obtain a global model which serves as the new start of local models. However, since DNN is a highly non-convex model, averaging parameters cannot ensure that such global model can perform better than those local models. To tackle this problem, we introduce a new parallel training framework called Ensemble-Compression, denoted as EC-DNN. In this framework, we propose to aggregate the local models by ensemble, i.e., averaging the outputs of local models instead of the parameters. As most of prevalent loss functions are convex to the output of DNN, the performance of ensemble-based global model is guaranteed to be at least as good as the average performance of local models. However, a big challenge lies in the explosion of model size since each round of ensemble can give rise to multiple times size increment. Thus, we carry out model compression after each ensemble, specialized by a distillation based method in this paper, to reduce the size of the global model to be the same as the local ones. Our experimental results demonstrate the prominent advantage of EC-DNN over MA-DNN in terms of both accuracy and speedup.


title: Coideals, quantum subgroups and idempotent states

abstract: We establish a one to one correspondence between idempotent states on a locally compact quantum group G and integrable coideals in the von Neumann algebra of bounded measurable functions on G that are preserved by the scaling group. In particular we show that there is a 1-1 correspondence between idempotent states on G and psi-expected left invariant von Neumann subalgebras of bounded measurable functions on G. We characterize idempotent states of Haar type as those corresponding to integrable normal coideals preserved by the scaling group. We also establish a one to one correspondence between open subgroups of G and convolutionally central idempotent states on the dual of G. Finally we characterize coideals corresponding to open quantum subgroups of G as those that are normal and admit an atom. As a byproduct of this study we get a number of universal lifting results for Podles condition, normality and regularity and we generalize number of results known before to hold under the coamenability assumption.


title: Source-LDA: Enhancing probabilistic topic models using prior knowledge sources

abstract: A popular approach to topic modeling involves extracting co-occurring n-grams of a corpus into semantic themes. The set of n-grams in a theme represents an underlying topic, but most topic modeling approaches are not able to label these sets of words with a single n-gram. Such labels are useful for topic identification in summarization systems. This paper introduces a novel approach to labeling a group of n-grams comprising an individual topic. The approach taken is to complement the existing topic distributions over words with a known distribution based on a predefined set of topics. This is done by integrating existing labeled knowledge sources representing known potential topics into the probabilistic topic model. These knowledge sources are translated into a distribution and used to set the hyperparameters of the Dirichlet generated distribution over words. In the inference these modified distributions guide the convergence of the latent topics to conform with the complementary distributions. This approach ensures that the topic inference process is consistent with existing knowledge. The label assignment from the complementary knowledge sources are then transferred to the latent topics of the corpus. The results show both accurate label assignment to topics as well as improved topic generation than those obtained using various labeling approaches based off Latent Dirichlet allocation (LDA).


title: On the eigenfunctions for the multi-species q-Boson system

abstract: In a previous paper a multi-species version of the q-Boson stochastic particle system is introduced and the eigenfunctions of its backward generator are constructed by using a representation of the Hecke algebra. In this article we prove a formula which expresses the eigenfunctions by means of the q-deformed bosonic operators, which are constructed from the L-operator of higher rank found in the recent work by Garbali, de Gier and Wheeler. The L-operator is obtained from the universal R-matrix of the quantum affine algebra of type A_{r}^{(1)} by the use of the q-oscillator representation. Thus our formula may be regarded as a bridge between two approaches to studying integrable stochastic systems by means of the quantum affine algebra and the affine Hecke algebra.


title: Shape dependence and anisotropic finite-size scaling of the phase coherence of three-dimensional Bose-Einstein condensed gases

abstract: We investigate the equilibrium phase-coherence properties of Bose-condensed particle systems, focusing on their shape dependence and finite-size scaling (FSS). We consider three-dimensional (3D) homogeneous systems confined to anisotropic L x L x L_a boxes, below the BEC transition temperature $T_c$. We show that the phase correlations develop peculiar anisotropic FSS for any $T<T_c$, in the large-$L$ limit keeping the ratio \lambda = L_a/L^2 fixed. This phenomenon is effectively described by the 3D spin-wave (SW) theory. Its universality is confirmed by quantum Monte Carlo simulations of the 3D Bose-Hubbard model in the BEC phase. The phase-coherence properties of very elongated BEC systems, \lambda>>1, are characterized by a coherence length \xi_a \sim A_t \rho_s/T where A_t is the transverse area and \rho_s is the superfluid density.


title: Steric effects of CO2 binding to transition metal-benzene complexes: a first-principles study

abstract: Using density functional theory (DFT) calculations, we investigated the adsorption of CO2 molecules on 3d transition metal (TM)-benzene complexes. Our calculations show that the maximum number of CO2 molecules adsorbable on Sc or Ti atoms is three, but the 18-electron rule predicts it should be four. The 18-electron rule is generally successful in predicting the maximum H2 adsorption number for TM atoms including Sc or Ti atoms. We found that the 18-electron rule fails to correctly predict CO2 binding on Sc- or Ti-benzene complexes because CO2 binding, in contrast to H2 binding, requires additional consideration for steric hindrance due to the large bond length of CO2. We calculated the occupation function for CO2 using the Tolman cone angle, which shows that three CO2 molecules fully occupy the available space around Sc- and Ti-benzene complexes. This estimation is the same maximum CO2 adsorption number predicted by DFT calculations. Therefore, we propose that the occupation function for CO2 using the Tolman cone angle is an efficient model for evaluating steric hindrance of CO2 adsorption on a surface.


title: Online Budgeted Truthful Matching

abstract: An online truthful budgeted matching problem is considered for a bipartite graph, where the right vertices are available ahead of time, and individual left vertices arrive sequentially. On arrival of a left vertex, its edge utilities (or weights) to all the right vertices and a corresponding cost (or bid) are revealed. If a left vertex is matched to any of the right vertices, then it has to be paid at least as much as its cost. The problem is to match each left vertex instantaneously and irrevocably to any one of the right vertices, if at all, to find the maximum weight matching that is truthful, under a payment budget constraint. Truthfulness condition requires that no left vertex has any incentive of misreporting its cost. Assuming that the vertices arrive in an uniformly random order (secretary model) with arbitrary utilities, a truthful algorithm is proposed that is $24\beta$-competitive (where $\beta$ is the ratio of the maximum and the minimum utility) and satisfies the payment budget constraint. Direct applications of this problem include crowdsourcing auctions, and matching wireless users to cooperative relays in device-to-device enabled cellular network.


title: Dynamics of 3D isolated thermal filaments

abstract: Simulations have been carried out to establish how electron thermal physics, introduced in the form of a dynamic electron temperature, affects isolated filament motion and dynamics in 3D. It is found that thermal effects impact filament motion in two major ways when the filament has a significant temperature perturbation compared to its density perturbation: They lead to a strong increase in filament propagation in the bi-normal direction and a significant decrease in net radial propagation. Both effects arise from the temperature dependence of the sheath current which leads to a non-uniform floating potential, with the latter effect supplemented by faster pressure loss. The reduction in radial velocity can only occur when the filament cross-section loses angular symmetry. The behaviour is observed across different filament sizes and suggests that filaments with much larger temperature perturbations than density perturbations are more strongly confined to the near SOL region.