title: Loop corrections in spin models through density consistency

abstract: Computing marginal distributions of discrete or semidiscrete Markov random fields (MRFs) is a fundamental, generally intractable problem with a vast number of applications in virtually all fields of science. We present a new family of computational schemes to approximately calculate the marginals of discrete MRFs. This method shares some desirable properties with belief propagation, in particular, providing exact marginals on acyclic graphs, but it differs with the latter in that it includes some loop corrections; i.e., it takes into account correlations coming from all cycles in the factor graph. It is also similar to the adaptive Thouless-Anderson-Palmer method, but it differs with the latter in that the consistency is not on the first two moments of the distribution but rather on the value of its density on a subset of values. The results on finite-dimensional Isinglike models show a significant improvement with respect to the Bethe-Peierls (tree) approximation in all cases and with respect to the plaquette cluster variational method approximation in many cases. In particular, for the critical inverse temperature $\beta_{c}$ of the homogeneous hypercubic lattice, the expansion of $\left(d\beta_{c}\right)^{-1}$ around $d=\infty$ of the proposed scheme is exact up to the $d^{-4}$ order, whereas the two latter are exact only up to the $d^{-2}$ order.

categories: cond-mat.stat-mech cond-mat.dis-nn cs.LG physics.comp-ph


title: The bulk-edge correspondence for continuous dislocated systems

abstract: We study topological aspects of defect modes for a family of operators $\{\mathscr{P}(t)\}_{t \in [0,2\pi]}$ on $L^2(\mathbb{R})$. $\mathscr{P}(t)$ is a periodic Schr\"odinger operator $P_0$ perturbed by a dislocated potential. This potential is periodic on the left and on the right, but acquires a phase defect $t$ from $-\infty$ relative to $+\infty$. When $t=\pi$ and the dislocation is small and adiabatic, Fefferman, Lee-Thorp and Weinstein showed in previous work that Dirac points of $P_0$ (degeneracies in the band spectrum of $P_0$) bifurcate to defect modes of $\mathscr{P}(\pi)$. We show that these modes are topologically protected at the level of the family $\{\mathscr{P}(t)\}_{t \in [0,2\pi]}$. This means that local perturbations cannot remove these states for all values of $t$, even outside the small adiabatic regime studied by Fefferman$-$Lee-Thorp$-$Weinstein. We define two topological quantities: an edge index (the signed number of eigenvalues crossing an energy gap as $t$ runs from $0$ to $2\pi$) and a bulk index (the Chern number of a Bloch eigenbundle for the periodic operator near $+\infty$). We prove that these indexes are equal to the same odd winding number. This shows the bulk-edge correspondence for the family $\{\mathscr{P}(t)\}_{t \in [0,2\pi]}$. We express this winding number in terms of the asymptotic shape of the dislocation and of the Dirac point Bloch modes. We illustrate the topological depth of our model via the computation of the bulk/edge index on a few examples.


title: Revealed Stochastic Preference: A One-Paragraph Proof and Generalization

abstract: McFadden and Richter (1991) and later McFadden (2005) show that the Axiom of Revealed Stochastic Preference characterizes rationalizability of choice probabilities through random utility models on finite universal choice spaces. This note proves the result in one short, elementary paragraph and extends it to set valued choice. The latter requires a different axiom than is reported in McFadden (2005).


title: Generation for Lagrangian cobordisms in Weinstein manifolds

abstract: We prove that Lagrangian cocores and Lagrangian linking disks of a stopped Weinstein manifold generate the Lagrangian cobordism infinity-category. As a geometric consequence, we see that any brane (after stabilization) admits a Lagrangian cobordism to a disjoint union of some standard collection of branes (cocores, linking disks, and a zero object). For example, when our stopped Weinstein manifold is a point stopped by itself, we find that any exact brane in Euclidean space admits a Lagrangian cobordism to a disjoint union of cotangent fibers and a zero object. (This is a stronger statement than one could obtain from purely Fukaya-categorical generation results.) Our methods are constructive. For example, when our Weinstein manifold is a point, after stabilization we can resolve the conormal to a compact manifold A of R^n by a sequence of cotangent fibers; the resulting filtration realizes, after passage to the wrapped Fukaya category, the Morse cochain complex of A associated to (and hence filtered by) a generic ``distance to a point'' function; the associated gradeds are the reduced homologies of the Morse attaching spheres. There is also an algebraic consequence. Lagrangian cobordism theory is conjectured (in analogue to classical cobordism theory) to be linear over a ring spectrum L controlling Lagrangian cobordisms between cotangent fibers in Euclidean spaces. Our main theorem gives strong evidence for this conjecture: The infinity-category of Lagrangians and their cobordisms in R^infinity is equivalent to a full subcategory of modules over L. We conclude by proving a \pi_0-level theorem that gives further evidence of the above conjecture: We exhibit a \pi_0-level symmetric monoidal structure compatible with the linear structure of L-modules.


title: Weak Hadamard Star Configurations and Apolarity

abstract: In [7] the authors have introduced a new construction using the Hadamard product to present star configurations of codimension $ c $ of $\mathbb{P}^n$ and which they called Hadamard star configurations. In this paper, we introduce a more general type of Hadamard star configuration. Any star configuration constructed by our approach is called a weak Hadamard star configuration. We classify weak Hadamard star configurations, and in the case $ c=n $, we investigate the existence of a (weak) Hadamard star configuration apolar to the generic homogeneous polynomials of degree $ d$.


title: Double emulsion drop evaporation and resurfacing of daughter droplet

abstract: In this study, we present experimental and theoretical analyses of double emulsion drop evaporation. After the apparent completion of evaporation of the inner phase of a double emulsion drop, surprisingly, a resurfacing of a daughter droplet is observed. We further investigated to hypothesize this phenomenon which allowed us to obtain a prolonged fixed contact line evaporation for a single phase drop along with similar occurrence of resurfacing as of the double emulsion drops.


title: Many-body localization in Landau level subbands

abstract: We explore the problem of localization in topological and non-topological nearly-flat subbands derived from the lowest Landau level, in the presence of quenched disorder and short-range interactions. We consider two models: a suitably engineered periodic potential, and randomly distributed point-like impurities. We perform numerical exact diagonalization on a torus geometry and use the mean level spacing ratio $\langle r \rangle$ as a diagnostic of ergodicity. For topological subbands, we find there is no ergodicity breaking in both the one and two dimensional thermodynamic limits. For non-topological subbands, in constrast, we find evidence of an ergodicity breaking transition at finite disorder strength in the one-dimensional thermodynamic limit. Intriguingly, indications of similar behavior in the two-dimensional thermodynamic limit are found, as well. This constitutes a novel, $\textit{continuum}$ setting for the study of the many-body localization transition in one and two dimensions.


title: Single Channel Josephson Effect in a High Transmission Atomic Contact

abstract: The Josephson effect in scanning tunneling microscopy (STM) is an excellent tool to probe the properties of the superconducting order parameter on a local scale through the Ambegaokar-Baratoff (AB) relation. Using single atomic contacts created by means of atom manipulation, we demonstrate that in the extreme case of a single transport channel through the atomic junction modifications of the current-phase relation lead to significant deviations from the linear AB formula relating the critical current to the involved gap parameters. Using the full current-phase relation for arbitrary channel transmission, we model the Josephson effect in the dynamical Coulomb blockade regime because the charging energy of the junction capacitance cannot be neglected. We find excellent agreement with the experimental data. Projecting the current-phase relation onto the charge transfer operator shows that at high transmission multiple Cooper pair tunneling may occur. These deviations become non-negligible in Josephson-STM, for example, when scanning across single adatoms.


title: Muon Colliders, Neutrino Factories, and Results from the MICE Experiment

abstract: Muon colliders and neutrino factories are attractive options for future facilities aimed at achieving the highest lepton-antilepton collision energies and precision measurements of parameters of the Higgs boson and the neutrino mixing matrix. The performance and cost of these depend on how well a beam of muons can be cooled. Recent progress in muon cooling design studies and prototype tests nourishes the hope that such facilities can be built starting in the coming decade. The status of the key technologies and their various demonstration experiments is summarized, with emphasis on recent results from the Muon Ionization Cooling Experiment (MICE).


title: Phenomenological consequences of introducing new fermions with exotic charges to $R(K^{(*)})$, muon (g-2), the primordial Lithium problem, and dark matter

abstract: In this work we show that, by introducing two $SU(3)_{c}\times SU(2)_{L}$ singlet right handed fermions carrying opposite $U(1)_{Y}$ charges and while their left handed counterparts are singlet under $SU(3)_{c}\times SU(2)_{L}$ and neutral under $U(1)_{Y}$, in the regime where the new charged lepton masses are in the electroweak scale it will be able to explain the small neutrino masses via minimum-inverse seesaw scenario (MISS) as well as the reported $R(K^{(*)})$ and muon $(g-2)$ discrepancies. Also when the charged fermion masses are well above the electro weak scale the model can not explain the reported $R(K^{(*)})$ and muon $(g-2)$ discrepancies, but in this regime the model could explain the primordial Lithium problem. The model have another interesting extension where it can produce a stable and singlet under strong interaction scalar baryon that could constitute much of the dark matter mass of the universe which lead to a simple model where ordinary matter and dark-matter having same origin. The model can also provide new annihilation channels for the scalar singlet DM as well as allowing a doubly charged scalar whose signatures could show up in HL-LHC, ILC, CEPC etc.


title: Continual Classification Learning Using Generative Models

abstract: Continual learning is the ability to sequentially learn over time by accommodating knowledge while retaining previously learned experiences. Neural networks can learn multiple tasks when trained on them jointly, but cannot maintain performance on previously learned tasks when tasks are presented one at a time. This problem is called catastrophic forgetting. In this work, we propose a classification model that learns continuously from sequentially observed tasks, while preventing catastrophic forgetting. We build on the lifelong generative capabilities of [10] and extend it to the classification setting by deriving a new variational bound on the joint log likelihood, $\log p(x; y)$.


title: Convergent and divergent beam electron holography and reconstruction of adsorbates on free-standing two-dimensional crystals

abstract: Van der Waals heterostructures have been lately intensively studied because they offer a large variety of properties that can be controlled by selecting 2D materials and their sequence in the stack. The exact arrangement of the layers as well as the exact arrangement of the atoms within the layers, both are important for the properties of the resulting device. Recently it has been demonstrated that convergent beam electron diffraction (CBED) allows quantitative three-dimensional mapping of atomic positions in three-dimensional materials from a single CBED pattern. In this study we investigate CBED in more detail by simulating and performing various CBED regimes, with convergent and divergent wavefronts, on a somewhat simplified system: a 2D monolayer crystal. In CBED, each CBED spot is in fact an in-line hologram of the sample, where in-line holography is known to exhibit high intensity contrast in detection of weak phase objects that are not detectable in conventional in-focus imaging mode. Adsorbates exhibit strong intensity contrast in zero and higher order CBED spots, whereas lattice deformation such as strain or rippling cause noticeable intensity contrast only in the first and higher order CBED spots. The individual CBED spots can be reconstructed as typical in-line holograms, and the resolution of 2.13 A can be in principle achieved in the reconstructions. We provide simulated and experimental examples of CBED of a 2D monolayer crystal. The simulations show that individual CBED spots can be treated as in-line holograms and sample distributions such as adsorbates, can be reconstructed. Individual atoms can be reconstructed from a single CBED pattern provided the later exhibits high-order CBED spots. Examples of reconstructions obtained from experimental CBED patterns, at a resolution of 2.7 A, are shown.


title: Alleviating Patch Overfitting with Automatic Test Generation: A Study of Feasibility and Effectiveness for the Nopol Repair System

abstract: Among the many different kinds of program repair techniques, one widely studied family of techniques is called test suite based repair. However, test suites are in essence input-output specifications and are thus typically inadequate for completely specifying the expected behavior of the program under repair. Consequently, the patches generated by test suite based repair techniques can just overfit to the used test suite, and fail to generalize to other tests. We deeply analyze the overfitting problem in program repair and give a classification of this problem. This classification will help the community to better understand and design techniques to defeat the overfitting problem. We further propose and evaluate an approach called UnsatGuided, which aims to alleviate the overfitting problem for synthesis-based repair techniques with automatic test case generation. The approach uses additional automatically generated tests to strengthen the repair constraint used by synthesis-based repair techniques. We analyze the effectiveness of UnsatGuided: 1) analytically with respect to alleviating two different kinds of overfitting issues; 2) empirically based on an experiment over the 224 bugs of the Defects4J repository. The main result is that automatic test generation is effective in alleviating one kind of overfitting issue--regression introduction, but due to oracle problem, has minimal positive impact on alleviating the other kind of overfitting issue--incomplete fixing.


title: Scheduling computations with provably low synchronization overheads

abstract: Work Stealing has been a very successful algorithm for scheduling parallel computations, and is known to achieve high performances even for computations exhibiting fine-grained parallelism. We present a variant of \ws\ that provably avoids most synchronization overheads by keeping processors' deques entirely private by default, and only exposing work when requested by thieves. This is the first paper that obtains bounds on the synchronization overheads that are (essentially) independent of the total amount of work, thus corresponding to a great improvement, in both algorithm design and theory, over state-of-the-art \ws\ algorithms. Consider any computation with work $T_{1}$ and critical-path length $T_{\infty}$ executed by $P$ processors using our scheduler. Our analysis shows that the expected execution time is $O\left(\frac{T_{1}}{P} + T_{\infty}\right)$, and the expected synchronization overheads incurred during the execution are at most $O\left(\left(C_{CAS} + C_{MFence}\right)PT_{\infty}\right)$, where $C_{CAS}$ and $C_{MFence}$ respectively denote the maximum cost of executing a Compare-And-Swap instruction and a Memory Fence instruction.


title: Approximate continuous data assimilation of the 2D Navier-Stokes equations via the Voigt-regularization with observable data

abstract: We propose a data assimilation algorithm for the 2D Navier-Stokes equations, based on the Azouani, Olson, and Titi (AOT) algorithm, but applied to the 2D Navier-Stokes-Voigt equations. Adapting the AOT algorithm to regularized versions of Navier-Stokes has been done before, but the innovation of this work is to drive the assimilation equation with observational data, rather than data from a regularized system. We first prove that this new system is globally well-posed. Moreover, we prove that for any admissible initial data, the $L^2$ and $H^1$ norms of error are bounded by a constant times a power of the Voigt-regularization parameter $\alpha>0$, plus a term which decays exponentially fast in time. In particular, the large-time error goes to zero algebraically as $\alpha$ goes to zero. Assuming more smoothness on the initial data and forcing, we also prove similar results for the $H^2$ norm.


title: Two body relativistic wave equations

abstract: The relativistic quantum mechanics of two interacting particles is considered. We first present a covariant formulation of kinematics and of reduced phase space, giving a short outline of the classical results. We then quantize the systems for the scalar-scalar, fermion-scalar and fermion-fermion cases. We study the spectrum and the spherical waves solutions of the free case. The interaction with central scalar and vector potentials is introduced and the explicit equations are deduced. The one particle and the non relativistic limits are recovered and the general lines for the solution of the boundary value problems are given. We make a numerical analysis of the first two cases with Coulomb interaction. For the two fermions we largely revisit the model we had previously derived in order to uniformize the description for all the three cases. In order to give a complete review we report in Appendix some of the most interesting results obtained for atomic and mesonic systems with Coulomb and Cornell potential interactions respectively.


title: The CR geometry of weighted extremal Kahler and Sasaki metrics

abstract: We establish an equivalence between conformally Einstein--Maxwell Kahler 4-manifolds (recently studied in many works) and extremal Kahler 4-manifolds (in the sense of Calabi) with nowhere vanishing scalar curvature. The corresponding pairs of Kahler metrics arise as transversal Kahler structures of Sasaki metrics compatible with the same CR structure and having commuting Sasaki--Reeb vector fields. This correspondence extends to higher dimensions using the notion of a weighted extremal Kahler metric, illuminating and uniting several explicit constructions in Kahler and Sasaki geometry. It also leads to new existence and non-existence results for extremal Sasaki metrics, suggesting a link between notions of relative weighted K-stability for a polarized variety, and relative K-stability of the Kahler cone corresponding to a Sasaki polarization.


title: Using Personal Environmental Comfort Systems to Mitigate the Impact of Occupancy Prediction Errors on HVAC Performance

abstract: Heating, Ventilation and Air Conditioning (HVAC) consumes a significant fraction of energy in commercial buildings. Hence, the use of optimization techniques to reduce HVAC energy consumption has been widely studied. Model predictive control (MPC) is one state of the art optimization technique for HVAC control which converts the control problem to a sequence of optimization problems, each over a finite time horizon. In a typical MPC, future system state is estimated from a model using predictions of model inputs, such as building occupancy and outside air temperature. Consequently, as prediction accuracy deteriorates, MPC performance--in terms of occupant comfort and building energy use--degrades. In this work, we use a custom-built building thermal simulator to systematically investigate the impact of occupancy prediction errors on occupant comfort and energy consumption. Our analysis shows that in our test building, as occupancy prediction error increases from 5\% to 20\% the performance of an MPC-based HVAC controller becomes worse than that of even a simple static schedule. However, when combined with a personal environmental control (PEC) system, HVAC controllers are considerably more robust to prediction errors. Thus, we quantify the effectiveness of PECs in mitigating the impact of forecast errors on MPC control for HVAC systems.


title: Exact instantaneous optimals in the non-geostrophic Eady problem and the detrimental effects of discretization

abstract: We derive exact analytical expressions for flow configurations that optimize the instantaneous growth rate of energy in the linear Eady problem, along with the associated growth rates. These optimal perturbations are relevant linear stability analysis, but, more importantly, they are relevant for understanding the energetics of fully nonlinear baroclinic turbulence. The optimal perturbations and their growth rates are independent of the Richardson number. The growth rates of the optimal perturbations grow linearly as the horizontal wavelength of the perturbation decreases. Perturbation energy growth at large scales is driven by extraction of potential energy from the mean flow, while at small scales it is driven by extraction of kinetic energy from the mean shear. We also analyze the effect of spatial discretization on the optimal perturbations and their growth rates. A second order energy-conserving discretization on the Arakawa B grid generally has too-weak growth rates at small scales and is less accurate than two second order discretizations on the Arakawa C grid. The two C grid discretizations, one that conserves energy and another that conserves both energy and enstrophy, yield very similar optimal perturbation growth rates that are significantly more accurate than the B grid discretization at small scales.


title: A Reliability Model for Dependent and Distributed MDS Disk Array Units

abstract: Archiving and systematic backup of large digital data generates a quick demand for multi-peta byte scale storage systems. As drive capacities continue to grow beyond the few terabytes range to address the demands of today's cloud, the likelihood of having multiple/simultaneous disk failures become a reality. Among the main factors causing catastrophic system failures, correlated disk failures and the network bandwidth are reported to be the two common source of performance degradation. The emerging trend is to use efficient/sophisticated erasure codes (EC) equipped with multiple parities and efficient repairs in order to meet the reliability/bandwidth requirements. It is known that mean time to failure and repair rates reported by the disk manufacturers cannot capture life cycle patterns of distributed storage systems. In this study, we develop failure models based on generalized Markov chains that can accurately capture correlated performance degradations with multi-parity protection schemes based on modern Maximum Distance Separable (MDS) EC. Furthermore, we use the proposed model in a distributed storage scenario to quantify two example use cases: Primarily, the common sense that adding more parity disks are only meaningful if we have a decent decorrelation between the failure domains of storage systems and the reliability of generic multiple single-dimensional EC protected storage systems.


title: The work of Mike Hochman on multidimensional symbolic dynamics and Borel dynamics

abstract: This paper is based on my April 2018 talk at the Maryland-Penn State dynamics conference, on the occasion of Mike Hochman receiving the Brin Prize.


title: Microscopy of hydrogen and hydrogen-vacancy defect structures on graphene devices

abstract: We have used scanning tunneling microscopy (STM) to investigate two types of hydrogen defect structures on monolayer graphene supported by hexagonal boron nitride (h-BN) in a gated field-effect transistor configuration. The first H-defect type is created by bombarding graphene with 1-keV ionized hydrogen and is identified as two hydrogen atoms bonded to a graphene vacancy via comparison of experimental data to first-principles calculations. The second type of H defect is identified as dimerized hydrogen and is created by depositing atomic hydrogen having only thermal energy onto a graphene surface. Scanning tunneling spectroscopy (STS) measurements reveal that hydrogen dimers formed in this way open a new elastic channel in the tunneling conductance between an STM tip and graphene.


title: Quark and glue spectroscopy of scalars and pseudoscalars in SU(3) flavor limit

abstract: Within the framework of the generalized linear sigma model with glueballs recently proposed [1,2], we study the schematic spectroscopy of scalar and pseudoscalar mesons in the SU(3) flavor limit and explore their quark and glue contents. In this framework, for both scalars and pseudoscalars, the two octet physical states are admixtures of quark-antiquark and four-quark components, and the three singlet states contain quark-antiquark, four-quark and glue components. We identify the two scalar octets with $a_0(980)$ and $a_0(1450)$ and the two pseudoscalar octets with $\pi(137)$ and $\pi(1300)$. We show that, as expected, the light pseudoscalar octet is made dominantly of quark-antiquarks whereas the light scalar octet has a reversed substructure with a dominant four-quark component. The case of singlets is more complex due to surplus of states up to around 2 GeV. We consider all 35 permutations for identifying the three pseudoscalar singlets of our model with three of the seven experimental candidates. Our numerical simulation unambiguously identifies the lightest and the heaviest pseudoscalar singlets with $\eta(547)$ and $\eta(2225)$, respectively, and favors the identification of the middle singlet with either $\eta(1295)$ or $\eta(1405)$ [or, to a lesser extent, with $\eta(1475)$] and thereby allows a probe of their substructures. We then estimate the quark and glue components and find that the three pseudoscalar singlets (from lightest to heaviest) are mainly of quark-antiquark, four-quark and glue substructure, while the corresponding three scalar singlets (from lightest to heaviest) are of four-quark, quark antiquark and glue contents. The masses of pure pseudoscalar and scalar glueballs are estimated around 2.0 and 1.6 GeV, respectively.


title: Robust Adversarial Learning via Sparsifying Front Ends

abstract: It is by now well-known that small adversarial perturbations can induce classification errors in deep neural networks. In this paper, we take a bottom-up signal processing perspective to this problem and show that a systematic exploitation of sparsity in natural data is a promising tool for defense. For linear classifiers, we show that a sparsifying front end is provably effective against $\ell_{\infty}$-bounded attacks, reducing output distortion due to the attack by a factor of roughly $K/N$ where $N$ is the data dimension and $K$ is the sparsity level. We then extend this concept to deep networks, showing that a "locally linear" model can be used to develop a theoretical foundation for crafting attacks and defenses. We also devise attacks based on the locally linear model that outperform the well-known FGSM attack. We supplement our theoretical results with experiments on the MNIST and CIFAR-10 datasets, showing the efficacy of the proposed sparsity-based defense schemes.


abstract: We study K\"ahler-Dirac fermions on Euclidean dynamical triangulations. This fermion formulation furnishes a natural extension of staggered fermions to random geometries without requring vielbeins and spin connections. We work in the quenched approximation where the geometry is allowed to fluctuate but there is no back-reaction from the matter on the geometry. By examining the eigenvalue spectrum and the masses of scalar mesons we find evidence for a four fold degeneracy in the fermion spectrum in the large volume, continuum limit. It is natural to associate this degeneracy with the well known equivalence in continuum flat space between the K\"ahler-Dirac fermion and four copies of a Dirac fermion. Lattice effects then lift this degeneracy in a manner similar to staggered fermions on regular lattices. The evidence that these discretization effects vanish in the continuum limit suggests both that lattice continuum K\"ahler-Dirac fermions are recovered at that point, and that this limit truly corresponds to smooth continuum geometries. One additional advantage of the K\"ahler-Dirac action is that it respects an exact $U(1)$ symmetry on any random triangulation. This $U(1)$ symmetry is related to continuum chiral symmetry. By examining fermion bilinear condensates we find strong evidence that this $U(1)$ symmetry is not spontaneously broken in the model at order the Planck scale. This is a necessary requirement if models based on dynamical triangulations are to provide a valid ultraviolet complete formulation of quantum gravity.


title: Streaming Graph Neural Networks

abstract: Graphs are essential representations of many real-world data such as social networks. Recent years have witnessed the increasing efforts made to extend the neural network models to graph-structured data. These methods, which are usually known as the graph neural networks, have been applied to advance many graphs related tasks such as reasoning dynamics of the physical system, graph classification, and node classification. Most of the existing graph neural network models have been designed for static graphs, while many real-world graphs are inherently dynamic. For example, social networks are naturally evolving as new users joining and new relations being created. Current graph neural network models cannot utilize the dynamic information in dynamic graphs. However, the dynamic information has been proven to enhance the performance of many graph analytic tasks such as community detection and link prediction. Hence, it is necessary to design dedicated graph neural networks for dynamic graphs. In this paper, we propose DGNN, a new {\bf D}ynamic {\bf G}raph {\bf N}eural {\bf N}etwork model, which can model the dynamic information as the graph evolving. In particular, the proposed framework can keep updating node information by capturing the sequential information of edges (interactions), the time intervals between edges and information propagation coherently. Experimental results on various dynamic graphs demonstrate the effectiveness of the proposed framework.


title: Energy-conserving Galerkin approximations for quasigeostrophic dynamics

abstract: A method is presented for constructing energy-conserving Galerkin approximations in the vertical coordinate of the full quasigeostrophic model with active surface buoyancy. The derivation generalizes the approach of Rocha \emph{et al.} (2016) to allow for general bases. Details are then presented for a specific set of bases: Legendre polynomials for potential vorticity and a recombined Legendre basis from Shen (1994) for the streamfunction. The method is tested in the context of linear baroclinic instability calculations, where it is compared to the standard second-order finite-difference method and to a Chebyshev collocation method. The Galerkin scheme is quite accurate even for a small number of degrees of freedom $N$, and growth rates converge much more quickly with increasing $N$ for the Galerkin scheme than for the finite-difference scheme. The Galerkin scheme is at least as accurate as finite differences and can in some cases achieve the same accuracy as the finite difference scheme with ten times fewer degrees of freedom. The energy-conserving Galerkin scheme is of comparable accuracy to the Chebyshev collocation scheme in most linear stability calculations, but not in the Eady problem where the Chebyshev scheme is significantly more accurate. Finally the three methods are compared in the context of a simplified version of the nonlinear equations: the two-surface model with zero potential vorticity. The Chebyshev scheme is the most accurate, followed by the Galerkin scheme and then the finite difference scheme. All three methods conserve energy with similar accuracy, despite not having any a priori guarantee of energy conservation for the Chebyshev scheme. Further nonlinear tests with non-zero potential vorticity to assess the merits of the methods will be performed in a future work.


title: Tailorable Elasticity of Cantilever Using Spatio-Angular Functionally Graded Biomimetic Scales

abstract: Cantilevered beams are of immense importance as structural and sensorial members for a number of applications. Endowing tailorable elasticity can have wide ranging engineering ramification. Such tailorability could be possible using some type of spatial gradation in the beam's material or cross section. However, these often require extensive additive and subtractive material processing or specialized casts. In this letter, we demonstrate an alternative bio inspired mechanical pathway, which is based on exploiting the nonlinearity that would arise from a functionally graded distribution of biomimetic scales on the surface using an analytical approach. This functional gradation is geometrically sourced and could arise from either spatial or angular gradation of scales. We analyze such a functionally graded cantilever beam under uniform loading. In comparison with uniformly distributed scales, we find significant differences in bending stiffness for both spatial and angular gradations. Spatial and angular functional gradation share some universality but also sharp contrasts in their effect on the underlying beam. A combination of both types of gradation in the structure can be used to alternatively increase or decrease stiffness and therefore a pathway to tailor the elasticity of a cantilever beam relatively easily. These results give rise to an architected framework for designing and optimizing the topography of leveraged solids.


title: Justification of the Lugiato-Lefever model from a damped driven $\phi^4$ equation

abstract: The Lugiato-Lefever equation is a damped and driven version of the well-known nonlinear Schr\"odinger equation. It is a mathematical model describing complex phenomena in dissipative and nonlinear optical cavities. Within the last two decades, the equation has gained a wide attention as it becomes the basic model describing optical frequency combs. Recent works derive the Lugiato-Lefever equation from a class of damped driven $\phi^4$ equations closed to resonance. In this paper, we provide a justification of the envelope approximation. From the analysis point of view, the result is novel and non-trivial as the drive yields a perturbation term that is not square integrable. The main approach proposed in this work is to decompose the solutions into a combination of the background and the integrable component. This paper is the first part of a two-manuscript series.


title: Minsum $k$-Sink Problem on Path Networks

abstract: We consider the problem of locating a set of $k$ sinks on a path network with general edge capacities that minimizes the sum of the evacuation times of all evacuees. We first present an $O(kn\log^4n)$ time algorithm when the edge capacities are non-uniform, where $n$ is the number of vertices. We then present an $O(kn\log^3 n)$ time algorithm when the edge capacities are uniform. We also present an $O(n\log n)$ time algorithm for the special case where $k=1$ and the edge capacities are non-uniform.


title: On the analysis of scheduling algorithms for structured parallel computations

abstract: Algorithms for scheduling structured parallel computations have been widely studied in the literature. For some time now, Work Stealing is one of the most popular for scheduling such computations, and its performance has been studied in both theory and practice. Although it delivers provably good performances, the effectiveness of its underlying load balancing strategy is known to be limited for certain classes of computations, particularly the ones exhibiting irregular parallelism (e.g. depth first searches). Many studies have addressed this limitation from a purely load balancing perspective, viewing computations as sets of independent tasks, and then analyzing the expected amount of work attached to each processor as the execution progresses. However, these studies make strong assumptions regarding work generation which, despite being standard from a queuing theory perspective --- where work generation can be assumed to follow some random distribution --- do not match the reality of structured parallel computations --- where the work generation is not random, only depending on the structure of a computation. In this paper, we introduce a formal framework for studying the performance of structured computation schedulers, define a criterion that is appropriate for measuring their performance, and present a methodology for analyzing the performance of randomized schedulers. We demonstrate the convenience of this methodology by using it to prove that the performance of Work Stealing is limited, and to analyze the performance of a Work Stealing and Spreading algorithm, which overcomes Work Stealing's limitation.


title: Strong laws of large numbers for arrays of random variables and stable random fields

abstract: Strong laws of large numbers are established for random fields with weak or strong dependence. These limit theorems are applicable to random fields with heavy-tailed distributions including fractional stable random fields. The conditions for SLLN are described in terms of the $p$-th moments of the partial sums of the random fields, which are convenient to verify. The main technical tool in this paper is a maximal inequality for the moments of partial sums of random fields that extends the technique of Levental, Chobanyan and Salehi \cite{chobanyan-l-s} for a sequence of random variables indexed by a one-parameter.


title: A frequency-preserving and time-invariant metamaterial-based nonlinear acoustic diode

abstract: We present the realization of an acoustic diode or rectifier, exploiting symmetry-breaking nonlinear effects like harmonic generation and wave mixing and the filtering capabilities of metamaterials. The essential difference and advantage compared with previous acoustic diode realizations is that the present is simultaneously a time invariant, frequency preserving and switchable device. This allows its application also as an on-off or amplitude-tuning switch. We evaluate its properties by means of a numerical study and demonstrate its feasibility in a preliminary experimental realization. This work may provide new opportunities for the practical realization of structural components with one-way wave propagation properties.


title: Lower Bounds for Oblivious Data Structures

abstract: An oblivious data structure is a data structure where the memory access patterns reveals no information about the operations performed on it. Such data structures were introduced by Wang et al. [ACM SIGSAC'14] and are intended for situations where one wishes to store the data structure at an untrusted server. One way to obtain an oblivious data structure is simply to run a classic data structure on an oblivious RAM (ORAM). Until very recently, this resulted in an overhead of $\omega(\lg n)$ for the most natural setting of parameters. Moreover, a recent lower bound for ORAMs by Larsen and Nielsen [CRYPTO'18] show that they always incur an overhead of at least $\Omega(\lg n)$ if used in a black box manner. To circumvent the $\omega(\lg n)$ overhead, researchers have instead studied classic data structure problems more directly and have obtained efficient solutions for many such problems such as stacks, queues, deques, priority queues and search trees. However, none of these data structures process operations faster than $\Theta(\lg n)$, leaving open the question of whether even faster solutions exist. In this paper, we rule out this possibility by proving $\Omega(\lg n)$ lower bounds for oblivious stacks, queues, deques, priority queues and search trees.


title: Compositional Set Invariance in Network Systems with Assume-Guarantee Contracts

abstract: This paper presents an assume-guarantee reasoning approach to the computation of robust invariant sets for network systems. Parameterized signal temporal logic (pSTL) is used to formally describe the behaviors of the subsystems, which we use as the template for the contract. We show that set invariance can be proved with a valid assume-guarantee contract by reasoning about individual subsystems. If a valid assume-guarantee contract with monotonic pSTL template is known, it can be further refined by value iteration. When such a contract is not known, an epigraph method is proposed to solve for a contract that is valid, ---an approach that has linear complexity for a sparse network. A microgrid example is used to demonstrate the proposed method. The simulation result shows that together with control barrier functions, the states of all the subsystems can be bounded inside the individual robust invariant sets.


title: Learning to Route Efficiently with End-to-End Feedback: The Value of Networked Structure

abstract: We introduce efficient algorithms which achieve nearly optimal regrets for the problem of stochastic online shortest path routing with end-to-end feedback. The setting is a natural application of the combinatorial stochastic bandits problem, a special case of the linear stochastic bandits problem. We show how the difficulties posed by the large scale action set can be overcome by the networked structure of the action set. Our approach presents a novel connection between bandit learning and shortest path algorithms. Our main contribution is an adaptive exploration algorithm with nearly optimal instance-dependent regret for any directed acyclic network. We then modify it so that nearly optimal worst case regret is achieved simultaneously. Driven by the carefully designed Top-Two Comparison (TTC) technique, the algorithms are efficiently implementable. We further conduct extensive numerical experiments to show that our proposed algorithms not only achieve superior regret performances, but also reduce the runtime drastically.


title: Electronic properties of the parabolic Dirac system

abstract: We investigate the electronic properties of the parabolic Dirac system. In the pres- ence of magnetic field, we discuss the thermodynamical potential and the anomalous diamagnetism due to the presence of Dirac cone. The integer Hall conductivity and the intrinsic Hall conductivity induced by the Berry curvature are also investigated. Considering the disorder-induced self-energy in self-consistent Born approximation, the RKKY interaction, long-range scattering rate and the spectral function are also calculated, which are distinct to that of the clean system (without impurity).


title: A Multilingual Study of Compressive Cross-Language Text Summarization

abstract: Cross-Language Text Summarization (CLTS) generates summaries in a language different from the language of the source documents. Recent methods use information from both languages to generate summaries with the most informative sentences. However, these methods have performance that can vary according to languages, which can reduce the quality of summaries. In this paper, we propose a compressive framework to generate cross-language summaries. In order to analyze performance and especially stability, we tested our system and extractive baselines on a dataset available in four languages (English, French, Portuguese, and Spanish) to generate English and French summaries. An automatic evaluation showed that our method outperformed extractive state-of-art CLTS methods with better and more stable ROUGE scores for all languages.


title: High-throughput Discovery of Topologically Non-trivial Materials using Spin-orbit Spillage

abstract: We present a novel methodology to identify topologically non-trivial materials based on band inversion induced by spin-orbit coupling (SOC) effect. Specifically, we compare the density functional theory (DFT) based wavefunctions with and without spin-orbit coupling and compute the spin-orbit-spillage as a measure of band-inversion. Due to its ease of calculation, without any need for symmetry analysis or dense k-point interpolation, the spillage is an excellent tool for identifying topologically non-trivial materials. Out of 30000 materials available in the JARVIS-DFT database, we applied this methodology to more than 4835 non-magnetic materials consisting of heavy atoms and low bandgaps. We found 1868 candidate materials with high-spillage (using 0.5 as a threshold). We validated our methodology by carrying out conventional Wannier-interpolation calculations for 289 candidate materials. We demonstrate that in addition to Z2 topological insulators, this screening method successfully identified many semimetals and topological crystalline insulators. Importantly, our approach is applicable to the investigation of disordered or distorted as well as magnetic materials, because it is not based on symmetry considerations. We discuss some individual example materials, as well as trends throughout our dataset, which is available at the websites: https://www.ctcms.nist.gov/~knc6/JVASP.html and https://jarvis.nist.gov/.


title: Predicting the Semantic Textual Similarity with Siamese CNN and LSTM

abstract: Semantic Textual Similarity (STS) is the basis of many applications in Natural Language Processing (NLP). Our system combines convolution and recurrent neural networks to measure the semantic similarity of sentences. It uses a convolution network to take account of the local context of words and an LSTM to consider the global context of sentences. This combination of networks helps to preserve the relevant information of sentences and improves the calculation of the similarity between sentences. Our model has achieved good results and is competitive with the best state-of-the-art systems.


title: Some Results On Relative Entropy in Quantum Field Theory

abstract: We prove that the mutual information for vacuum state as defined by Araki is finite for quantum field theory of free fermions on a Minkowski spacetime of any dimension. In the case of two dimensional chiral conformal field theory (CFT) we use our previous results for the free fermions to show that for a large class of chiral CFT the mutual information is finite. We also provide a direct relation between relative entropy and the index of a representation of conformal net.


title: Measurements of the Cross Spectra of the Cosmic Infrared and Microwave Backgrounds from 95 to 1200 GHz

abstract: We present measurements of the power spectra of cosmic infrared background (CIB) and cosmic microwave background (CMB) fluctuations in six frequency bands. Maps at the lower three frequency bands, 95, 150, and 220 GHz (3330, 2000, 1360 $\mu$m) are from the South Pole Telescope, while the upper three frequency bands, 600, 857, and 1200 GHz (500, 350, 250 $\mu$m) are observed with Herschel/SPIRE. From these data, we produce 21 angular power spectra (six auto- and fifteen cross-frequency) spanning the multipole range $600 \le \ell \le 11,000$. Our measurements are the first to cross-correlate measurements near the peak of the CIB spectrum with maps at 95 GHz, complementing and extending the measurements from Planck Collaboration et al. (2014) at 218, 550, and 857 GHz. The observed fluctuations originate largely from clustered, infrared-emitting, dusty star-forming galaxies, the CMB, and to a lesser extent radio galaxies, active galactic nuclei, and the Sunyaev-Zel'dovich effect.


title: Graph isomorphism and Gaussian boson sampling

abstract: We introduce a connection between a near-term quantum computing device, specifically a Gaussian boson sampler, and the graph isomorphism problem. We propose a scheme where graphs are encoded into quantum states of light, whose properties are then probed with photon-number-resolving detectors. We prove that the probabilities of different photon-detection events in this setup can be combined to give a complete set of graph invariants. Two graphs are isomorphic if and only if their detection probabilities are equivalent. We present additional ways that the measurement probabilities can be combined or coarse-grained to make experimental tests more amenable. We benchmark these methods with numerical simulations on the Titan supercomputer for several graph families: pairs of isospectral nonisomorphic graphs, isospectral regular graphs, and strongly regular graphs.

categories: quant-ph math-ph math.CO math.MP math.NT


title: The velocity anisotropy of the Milky Way satellite system

abstract: We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, $\beta$, as a function of Galactocentric distance and compare the MW satellite system with those of simulated MW-mass haloes from the APOSTLE and Auriga simulation suites. The anisotropy profile for the MW satellite system increases from $\beta\sim -2$ at $r\sim20$ kpc to $\beta\sim 0.5$ at $r\sim200$ kpc, indicating that satellites closer to the Galactic centre have tangentially-biased motions while those farther out have radially-biased motions. The motions of satellites around APOSTLE host galaxies are nearly isotropic at all radii, while the $\beta(r)$ profiles for satellite systems in the Auriga suite, whose host galaxies are substantially more massive in baryons than those in APOSTLE, are more consistent with that of the MW satellite system. This shape of the $\beta(r)$ profile may be attributed to the central stellar disc preferentially destroying satellites on radial orbits, or intrinsic processes from the formation of the Milky Way system.


title: Tuning from failed superconductor to failed insulator with magnetic field

abstract: Do charge modulations compete with electron pairing in high-temperature copper-oxide superconductors? We investigated this question by suppressing superconductivity in a stripe-ordered cuprate compound at low temperature with high magnetic fields. With increasing field, loss of three-dimensional superconducting order is followed by reentrant two-dimensional superconductivity and then an ultra-quantum metal phase. Circumstantial evidence suggests that the latter state is bosonic and associated with the charge stripes. These results provide experimental support to the theoretical perspective that local segregation of doped holes and antiferromagnetic spin correlations underlies the electron-pairing mechanism in cuprates.


title: Multi-level Memory for Task Oriented Dialogs

abstract: Recent end-to-end task oriented dialog systems use memory architectures to incorporate external knowledge in their dialogs. Current work makes simplifying assumptions about the structure of the knowledge base, such as the use of triples to represent knowledge, and combines dialog utterances (context) as well as knowledge base (KB) results as part of the same memory. This causes an explosion in the memory size, and makes the reasoning over memory harder. In addition, such a memory design forces hierarchical properties of the data to be fit into a triple structure of memory. This requires the memory reader to infer relationships across otherwise connected attributes. In this paper we relax the strong assumptions made by existing architectures and separate memories used for modeling dialog context and KB results. Instead of using triples to store KB results, we introduce a novel multi-level memory architecture consisting of cells for each query and their corresponding results. The multi-level memory first addresses queries, followed by results and finally each key-value pair within a result. We conduct detailed experiments on three publicly available task oriented dialog data sets and we find that our method conclusively outperforms current state-of-the-art models. We report a 15-25% increase in both entity F1 and BLEU scores.


title: Disordering, Clustering, and Laning Transitions in Particle Systems with Dispersion in the Magnus Term

abstract: We numerically examine a two-dimensional system of repulsively interacting particles with dynamics that are governed by both a damping term and a Magnus term. The magnitude of the Magnus term has one value for half of the particles and a different value for the other half of the particles. In the absence of a driving force, the particles form a triangular lattice, while when a driving force is applied, we find that there is a critical drive above which a Magnus-induced disordering transition can occur even if the difference in the Magnus term between the two particle species is as small as one percent. The transition arises due to the different Hall angles of the two species, which causes their motion to decouple at the critical drive. At higher drives, the disordered state can undergo both species and density phase separation into a density modulated stripe that is oriented perpendicular to the driving direction. We observe several additional phases that occur as a function of drive and Magnus force disparity, including a variety of density modulated diagonal laned phases. In general we find a much richer variety of states compared to systems of oppositely driven overdamped Yukawa particles. We discuss the implications of our work for skyrmion systems, where we predict that even for small skyrmion dispersities, a drive-induced disordering transition can occur along with clustering phases and pattern forming states.


title: On the Effectiveness of Type-based Control Flow Integrity

abstract: Control flow integrity (CFI) has received significant attention in the community to combat control hijacking attacks in the presence of memory corruption vulnerabilities. The challenges in creating a practical CFI has resulted in the development of a new type of CFI based on runtime type checking (RTC). RTC-based CFI has been implemented in a number of recent practical efforts such as GRSecurity Reuse Attack Protector (RAP) and LLVM-CFI. While there has been a number of previous efforts that studied the strengths and limitations of other types of CFI techniques, little has been done to evaluate the RTC-based CFI. In this work, we study the effectiveness of RTC from the security and practicality aspects. From the security perspective, we observe that type collisions are abundant in sufficiently large code bases but exploiting them to build a functional attack is not straightforward. Then we show how an attacker can successfully bypass RTC techniques using a variant of ROP attacks that respect type checking (called TROP) and also built two proof-of-concept exploits, one against Nginx web server and the other against Exim mail server. We also discuss practical challenges of implementing RTC. Our findings suggest that while RTC is more practical for applying CFI to large code bases, its policy is not strong enough when facing a motivated attacker.


title: Stationary distributions of the multi-type ASEP

abstract: We give a recursive construction of the stationary distribution of multi-type asymmetric simple exclusion processes on a finite ring or on the infinite line $Z$. The construction can be interpreted in terms of "multi-line diagrams" or systems of queues in tandem. Let $q$ be the asymmetry parameter of the system. The queueing construction generalises the one previously known for the totally asymmetric ($q=0$) case, by introducing queues in which each potential service is unused with probability $q^k$ when the queue-length is $k$. The analysis is based on the matrix product representation of Prolhac, Evans and Mallick. Consequences of the construction include: a simple method for sampling exactly from the stationary distribution for the system on a ring; results on common denominators of the stationary probabilities, expressed as rational functions of $q$ with non-negative integer coefficients; and probabilistic descriptions of "convoy formation" phenomena in large systems.


title: A gamma-ray periodic modulation in Globular Cluster 47 Tucanae

abstract: The Globular Cluster 47 Tucanae was firstly detected in gamma-rays by the Large Area Telescope (LAT) onboard the \emph{Fermi} Gamma-ray Space Telescope, and the gamma-ray emission has been widely attributed to the millisecond pulsars. In this work, we analyze the Fermi-LAT pass 8 data ranging from 2008 August to 2017 May and report the detection of a modulation with a period of $18.416\pm0.008$ hours at a significance level of $\sim 4.8\sigma$. This is the first time to detect a significant modulation with a period much longer than that of millisecond pulsars in gamma-rays from Globular Clusters. The periodic modulation signal appears in the {\it Swift}-BAT data as well. The phase-folded Chandra X-ray light curve of a point source may have provided an additional clue.


title: Masses, Sheets and Rigid SCFTs

abstract: We study mass deformations of certain three dimensional $\mathcal{N}=4$ Superconformal Field Theories (SCFTs) that have come to be called $T^\rho[G]$ theories. These are associated to tame defects of the six dimensional $(0,2)$ SCFT $X[\mathfrak{j}]$ for $\mathfrak{j}=A,D,E$. We describe these deformations using a refined version of the theory of sheets, a subject of interest in Geometric Representation Theory. In mathematical terms, we parameterize local mass-like deformations of the tamely ramified Hitchin integrable system and identify the subset of the deformations that do admit an interpretation as a mass deformation for the theories under consideration. We point out the existence of non-trivial Rigid SCFTs among these theories. We classify the Rigid theories within this set of SCFTs and give a description of their Higgs and Coulomb branches. We then study the implications for the endpoints of RG flows triggered by mass deformations in these 3d $\mathcal{N}=4$ theories. Finally, we discuss connections with the recently proposed idea of Symplectic Duality and describe some conjectures about its action.


title: Emergence of excitonic superfluid at topological-insulator surfaces

abstract: Excitons are spin integer particles that are predicted to condense into a coherent quantum state at sufficiently low temperature, and exciton condensates can be realized at much higher temperature than condensates of atoms because of strong Coulomb binding and small mass. Signatures of exciton condensation have been reported in double quantum wells1-4, microcavities5, graphene6, and transition metal dichalcogenides7. Nonetheless, transport of exciton condensates is not yet understood and it is unclear whether an exciton condensate is a superfluid8,9 or an insulating electronic crystal10,11. Topological insulators (TIs) with massless particles and unique spin textures12 have been theoretically predicted13 as a promising platform for achieving exciton condensation. Here we report experimental evidence of excitonic superfluid phase on the surface of three-dimensional (3D) TIs. We unambiguously confirmed that electrons and holes are paired into charge neutral bound states by the electric field independent photocurrent distributions. And we observed a millimetre-long transport distance of these excitons up to 40 K, which strongly suggests dissipationless propagation. The robust macroscopic quantum states achieved with simple device architecture and broadband photoexcitation at relatively high temperature are expected to find novel applications in quantum computations and spintronics.


title: Sample-Efficient Learning of Nonprehensile Manipulation Policies via Physics-Based Informed State Distributions

abstract: This paper proposes a sample-efficient yet simple approach to learning closed-loop policies for nonprehensile manipulation. Although reinforcement learning (RL) can learn closed-loop policies without requiring access to underlying physics models, it suffers from poor sample complexity on challenging tasks. To overcome this problem, we leverage rearrangement planning to provide an informative physics-based prior on the environment's optimal state-visitation distribution. Specifically, we present a new technique, Learning with Planned Episodic Resets (LeaPER), that resets the environment's state to one informed by the prior during the learning phase. We experimentally show that LeaPER significantly outperforms traditional RL approaches by a factor of up to 5X on simulated rearrangement. Further, we relax dynamics from quasi-static to welded contacts to illustrate that LeaPER is robust to the use of simpler physics models. Finally, LeaPER's closed-loop policies significantly improve task success rates relative to both open-loop controls with a planned path or simple feedback controllers that track open-loop trajectories. We demonstrate the performance and behavior of LeaPER on a physical 7-DOF manipulator in https://youtu.be/feS-zFq6J1c.


title: The SCUBA-2 Web Survey: I. Observations of CO(3-2) in hyper-luminous QSO fields

abstract: A primary goal of the SCUBA-2 Web survey is to perform tomography of the early inter-galactic medium by studying systems containing some of the brightest quasi-stellar objects (QSOs; 2.5<z<3.0) and nearby submillimetre galaxies. As a first step, this paper aims to characterize the galaxies that host the QSOs. To achieve this, a sample of 13 hyper-luminous (L_AGN>10^14 L_odot) QSOs with previous submillimetre continuum detections were followed up with CO(3-2) observations using the NOEMA interferometer. All but two of the QSOs are detected in CO(3-2); for one non-detection, our observations show a tentative 2sigma line at the expected position and redshift, and for the other non-detection we find only continuum flux density an order of magnitude brighter than the other sources. In three of the fields, a companion potentially suitable for tomography is detected in CO line emission within 25 arcsec of the QSO. We derive gas masses, dynamical masses and far-infrared luminosities, and show that the QSOs in our sample have similar properties as compared to less luminous QSOs and SMGs in the literature, despite the fact that their black-hole masses (which are proportional to L_AGN) are 1-2 orders of magnitude larger. We discuss two interpretations of these observations: this is due to selection effects, such as preferential face-on viewing angles and picking out objects in the tail ends of the scatter in host-galaxy mass and black-hole mass relationships; or the black hole masses have been overestimated because the accretion rates are super-Eddington.


title: Understand, Compose and Respond - Answering Visual Questions by a Composition of Abstract Procedures

abstract: An image related question defines a specific visual task that is required in order to produce an appropriate answer. The answer may depend on a minor detail in the image and require complex reasoning and use of prior knowledge. When humans perform this task, they are able to do it in a flexible and robust manner, integrating modularly any novel visual capability with diverse options for various elaborations of the task. In contrast, current approaches to solve this problem by a machine are based on casting the problem as an end-to-end learning problem, which lacks such abilities. We present a different approach, inspired by the aforementioned human capabilities. The approach is based on the compositional structure of the question. The underlying idea is that a question has an abstract representation based on its structure, which is compositional in nature. The question can consequently be answered by a composition of procedures corresponding to its substructures. The basic elements of the representation are logical patterns, which are put together to represent the question. These patterns include a parametric representation for object classes, properties and relations. Each basic pattern is mapped into a basic procedure that includes meaningful visual tasks, and the patterns are composed to produce the overall answering procedure. The UnCoRd (Understand Compose and Respond) system, based on this approach, integrates existing detection and classification schemes for a set of object classes, properties and relations. These schemes are incorporated in a modular manner, providing elaborated answers and corrections for negative answers. In addition, an external knowledge base is queried for required common-knowledge. We performed a qualitative analysis of the system, which demonstrates its representation capabilities and provide suggestions for future developments.


title: Gravitational Lensing by Black Holes in Einsteinian Cubic Gravity

abstract: We investigate the predictions of Einsteinian Cubic Gravity (ECG) for the lensing effects due to supermassive black holes at the center of Milky Way and other galaxies. Working in the context of spherical symmetry, we obtain the metric function from a continued fraction method and find that both time delays and the angular positions of images considerably deviate from general relativity, as large as milliarcseconds. This suggests that observational tests of ECG are indeed feasible.


title: Sports Camera Calibration via Synthetic Data

abstract: Calibrating sports cameras is important for autonomous broadcasting and sports analysis. Here we propose a highly automatic method for calibrating sports cameras from a single image using synthetic data. First, we develop a novel camera pose engine. The camera pose engine has only three significant free parameters so that it can effectively generate a lot of camera poses and corresponding edge (i.e, field marking) images. Then, we learn compact deep features via a siamese network from paired edge image and camera pose and build a feature-pose database. After that, we use a novel two-GAN (generative adversarial network) model to detect field markings in real images. Finally, we query an initial camera pose from the feature-pose database and refine camera poses using truncated distance images. We evaluate our method on both synthetic and real data. Our method not only demonstrates the robustness on the synthetic data but also achieves the state-of-the-art accuracy on a standard soccer dataset and very high performance on a volleyball dataset.


title: Combinatorial Optimization with Graph Convolutional Networks and Guided Tree Search

abstract: We present a learning-based approach to computing solutions for certain NP-hard problems. Our approach combines deep learning techniques with useful algorithmic elements from classic heuristics. The central component is a graph convolutional network that is trained to estimate the likelihood, for each vertex in a graph, of whether this vertex is part of the optimal solution. The network is designed and trained to synthesize a diverse set of solutions, which enables rapid exploration of the solution space via tree search. The presented approach is evaluated on four canonical NP-hard problems and five datasets, which include benchmark satisfiability problems and real social network graphs with up to a hundred thousand nodes. Experimental results demonstrate that the presented approach substantially outperforms recent deep learning work, and performs on par with highly optimized state-of-the-art heuristic solvers for some NP-hard problems. Experiments indicate that our approach generalizes across datasets, and scales to graphs that are orders of magnitude larger than those used during training.


title: The Case for Formation of ISP-Content Providers Consortiums by Nash Bargaining for Internet Content Delivery

abstract: The formation of consortiums of a broadband access Internet Service Provider (ISP) and multiple Content Providers (CP) is considered for large-scale content caching. The consortium members share costs from operations and investments in the supporting infrastructure. Correspondingly, the model's cost function includes marginal and fixed costs; the latter has been important in determining industry structure. Also, if Net Neutrality regulations permit, additional network capacity on the ISP's last mile may be contracted by the CPs. The number of subscribers is determined by a combination of users' price elasticity of demand and Quality of Experience. The profit generated by a coalition after pricing and design optimization determines the game's characteristic function. Coalition formation is by a bargaining procedure due to Okada (1996) based on random proposers in a non-cooperative, multi-player game-theoretic framework. A necessary and sufficient condition is obtained for the Grand Coalition to form, which bounds subsidies from large to small contributors. Caching is generally supported even under Net Neutrality regulations. The Grand Coalition's profit matches upper bounds. Numerical results illustrate the analytic results.


title: The Second-Price Knapsack Problem: Near-Optimal Real Time Bidding in Internet Advertisement

abstract: In many online advertisement (ad) exchanges, ad slots are each sold via a separate second-price auction. This paper considers the bidder's problem of maximizing the value of ads they purchase in these auctions, subject to budget constraints. This 'second-price knapsack' problem presents challenges when devising a bidding strategy because of the uncertain resource consumption: bidders win if they bid the highest amount, but pay the second-highest bid, unknown a priori. This is in contrast to the traditional online knapsack problem, where posted prices are revealed when ads arrive, and for which there exists a rich literature of primal and dual algorithms. The main results of this paper establish general methods for adapting these primal and dual online knapsack selection algorithms to the second-price knapsack problem, where the prices are revealed only after bidding. In particular, a methodology is provided for converting deterministic and randomized knapsack selection algorithms into second-price knapsack bidding strategies, that purchase ads through an equivalent set of criteria and thereby achieve the same competitive guarantees. This shows a connection between the traditional knapsack selection algorithm and second-price auction bidding algorithms, that has not previously been leveraged. Empirical analysis on real ad exchange data verifies the usefulness of this method, and gives examples where it can outperform state-of-the-art techniques.


title: Multi-Channel Auto-Encoder for Speech Emotion Recognition

abstract: Inferring emotion status from users' queries plays an important role to enhance the capacity in voice dialogues applications. Even though several related works obtained satisfactory results, the performance can still be further improved. In this paper, we proposed a novel framework named multi-channel auto-encoder (MTC-AE) on emotion recognition from acoustic information. MTC-AE contains multiple local DNNs based on different low-level descriptors with different statistics functions that are partly concatenated together, by which the structure is enabled to consider both local and global features simultaneously. Experiment based on a benchmark dataset IEMOCAP shows that our method significantly outperforms the existing state-of-the-art results, achieving $64.8\%$ leave-one-speaker-out unweighted accuracy, which is $2.4\%$ higher than the best result on this dataset.


title: TYC 5353-1137-1: an enigmatic Double Periodic Variable of semiregular amplitude

abstract: We have performed a new search for DPVs of short period in the ASAS catalog (Pojmanski, G., 1997), focusing on those stars with orbital periods between 2 to 3 days which also show variations in their brightness. From a total of 244 objects, we have found another candidate to DPV, one whose mean brightness is gradually decreasing. By fitting a 3rd order polynomial to the mean magnitude and then moving it to zero for a second analysis, a gradual decrease over 2500 days was revealed. During the last 1000 days of this decrease, a 42% increase in the variation between the minimum and maximum values of the magnitude was observed. We determined the orbital period by using the PDM IRAF software (Stellingwerf 1978) and estimated the errors for the orbital period and long cycle by visual inspection of the light curves phased with trial periods near the minimum of the periodogram given by PDM.


title: Automated Process Incorporating Machine Learning Segmentation and Correlation of Oral Diseases with Systemic Health

abstract: Imaging fluorescent disease biomarkers in tissues and skin is a non-invasive method to screen for health conditions. We report an automated process that combines intraoral fluorescent porphyrin biomarker imaging, clinical examinations and machine learning for correlation of systemic health conditions with periodontal disease. 1215 intraoral fluorescent images, from 284 consenting adults aged 18-90, were analyzed using a machine learning classifier that can segment periodontal inflammation. The classifier achieved an AUC of 0.677 with precision and recall of 0.271 and 0.429, respectively, indicating a learned association between disease signatures in collected images. Periodontal diseases were more prevalent among males (p=0.0012) and older subjects (p=0.0224) in the screened population. Physicians independently examined the collected images, assigning localized modified gingival indices (MGIs). MGIs and periodontal disease were then cross-correlated with responses to a medical history questionnaire, blood pressure and body mass index measurements, and optic nerve, tympanic membrane, neurological, and cardiac rhythm imaging examinations. Gingivitis and early periodontal disease were associated with subjects diagnosed with optic nerve abnormalities (p <0.0001) in their retinal scans. We also report significant co-occurrences of periodontal disease in subjects reporting swollen joints (p=0.0422) and a family history of eye disease (p=0.0337). These results indicate cross-correlation of poor periodontal health with systemic health outcomes and stress the importance of oral health screenings at the primary care level. Our screening process and analysis method, using images and machine learning, can be generalized for automated diagnoses and systemic health screenings for other diseases.


title: Engaging Image Captioning Via Personality

abstract: Standard image captioning tasks such as COCO and Flickr30k are factual, neutral in tone and (to a human) state the obvious (e.g., "a man playing a guitar"). While such tasks are useful to verify that a machine understands the content of an image, they are not engaging to humans as captions. With this in mind we define a new task, Personality-Captions, where the goal is to be as engaging to humans as possible by incorporating controllable style and personality traits. We collect and release a large dataset of 201,858 of such captions conditioned over 215 possible traits. We build models that combine existing work from (i) sentence representations (Mazare et al., 2018) with Transformers trained on 1.7 billion dialogue examples; and (ii) image representations (Mahajan et al., 2018) with ResNets trained on 3.5 billion social media images. We obtain state-of-the-art performance on Flickr30k and COCO, and strong performance on our new task. Finally, online evaluations validate that our task and models are engaging to humans, with our best model close to human performance.


title: Coding for Computing Arbitrary Functions Unknown to the Encoder

abstract: In this paper we consider point-to-point and distributed source coding problems where the receiver is only interested in a function of the data sent by the source encoder(s), while knowledge of the function remains unknown to the encoder(s). We find the rate region for these problems, and in particular, show that if the destination is interested in computing a non-bijective function then the rate region for the point-to-point source coding problem expands over the entropy, and the rate region over the distributed source coding problem expands over the Slepian-Wolf rate region. A novel proof technique, similar to random binning, is developed to prove these results.


title: Truncated Back-propagation for Bilevel Optimization

abstract: Bilevel optimization has been recently revisited for designing and analyzing algorithms in hyperparameter tuning and meta learning tasks. However, due to its nested structure, evaluating exact gradients for high-dimensional problems is computationally challenging. One heuristic to circumvent this difficulty is to use the approximate gradient given by performing truncated back-propagation through the iterative optimization procedure that solves the lower-level problem. Although promising empirical performance has been reported, its theoretical properties are still unclear. In this paper, we analyze the properties of this family of approximate gradients and establish sufficient conditions for convergence. We validate this on several hyperparameter tuning and meta learning tasks. We find that optimization with the approximate gradient computed using few-step back-propagation often performs comparably to optimization with the exact gradient, while requiring far less memory and half the computation time.


title: The Boca-Cobeli-Zaharescu Map Analogue for the Hecke Triangle Groups $G_q$

abstract: The Farey sequence $\mathcal{F}(Q)$ at level $Q$ is the sequence of irreducible fractions in $[0, 1]$ with denominators not exceeding $Q$, arranged in increasing order of magnitude. A simple ``next-term'' algorithm exists for generating the elements of $\mathcal{F}(Q)$ in increasing or decreasing order. That algorithm, along with a number of other properties of the Farey sequence, was encoded by F. Boca, C. Cobeli, and A. Zaharescu into what is now known as the Boca-Cobeli-Zaharescu (BCZ) map, and used to attack several problems that can be described using the statistics of subsets of the Farey sequence. In this paper, we derive the Boca-Cobeli-Zaharescu map analogue for the discrete orbits $\Lambda_q = G_q(1, 0)^T$ of the linear action of the Hecke triangle groups $G_q$ on the plane $\mathbb{R}^2$ starting with a Stern-Brocot tree analogue for the said orbits. We derive the next-term algorithm for generating the elements of $\Lambda_q$ in vertical strips in increasing order of slope, and present a number of applications to the statistics of $\Lambda_q$.


title: Model Selection using Multi-Objective Optimization

abstract: Choices in scientific research and management require balancing multiple, often competing objectives.Multiple-objective optimization (MOO) provides a unifying framework for solving multiple objective problems. Model selection is a critical component to scientific inference and prediction and concerns balancing the competing objectives of model fit and model complexity. The tradeoff between model fit and model complexity provides a basis for describing the model-selection problem within the MOO framework. We discuss MOO and two strategies for solving the MOO problem; modeling preferences pre-optimization and post-optimization. Most model selection methods are consistent with solving MOO problems via specification of preferences pre-optimization. We reconcile these methods within the MOO framework. We also consider model selection using post-optimization specification of preferences. That is, by first identifying Pareto optimal solutions, and then selecting among them. We demonstrate concepts with an ecological application of model selection using avian species richness data in the continental United States.


title: Electron dynamics in laser and quasi-static transverse electric and longitudinal magnetic fields

abstract: By deriving the 3/2 dimensional Hamiltonian equations for electrons in the intense laser radiation and quasi-static transverse electric and longitudinal magnetic fields, the electron heating mechanisms are examined both for low harmonic resonance of electron frequency in the static fields with the laser frequency and for high harmonic resonances where the overlapping of broadened resonances causes the stochastic heating. For both cases, the maximum electron kinetic energies, well beyond the ponderomotive scaling, depend only on a small parameter combining the laser amplitude, electrostatic field strength and the conserved dephasing rate.


title: Renormalization of six-dimensional Yang-Mills theory in a background gauge field

abstract: Using the background field method, we study in a general covariant gauge the renormalization of the 6-dimensional Yang-Mills theory. This requires background gauge invariant counterterms, some of which do not vanish on shell. Such counterterms occur, even off-shell, with gauge-independent coefficients. The analysis is done at one loop order and the extension to higher orders is discussed by means of the BRST identities. We examine the behaviour of the beta function, which implies that this theory is not asymptotically free.


title: First order formulation of the Yang-Mills theory in a background field

abstract: The background gauge renormalization of the first order formulation of the Yang-Mills theory is studied by using the BRST identities. Together with the background symmetry, these identities allow for an iterative proof of renormalizability to all orders in perturbation theory. However, due to the fact that certain improper diagrams which violate the BRST symmetry should be removed, the renormalizability must be deduced indirectly. The recursive method involves rescalings and mixings of the fields, which lead to a renormalized effective action for the background field theory.


title: On Restricting to One Loop Order the Radiative Effects in Quantum Gravity

abstract: The dimensionful nature of the coupling in the Einstein-Hilbert action in four dimensions implies that the theory is non-renormalizable; explicit calculation shows that beginning at two loop order, divergences arise that cannot be removed by renormalization without introducing new terms in the classical action. It has been shown that, by use of a Lagrange multiplier field to ensure that the classical equation of motion is satisfied in the path integral, radiative effects can be restricted to one loop order. We show that by use of such Lagrange multiplier fields, the Einstein-Hilbert action can be quantized without the occurrence of non-renormalizable divergences. We then apply this mechanism to a model in which there is in addition to the Einstein-Hilbert action, a fully covariant action for a self-interacting scalar field coupled to the metric. It proves possible to restrict loop diagrams involving internal lines involving the metric to one-loop order; diagrams in which the scalar field propagates occur at arbitrary high order in the loop expansion. This model also can be shown to be renormalizable. Incorporating spinor and vector fields in the same way as scalar fields is feasible, and so a fully covariant Standard Model with a dynamical metric field can also be shown to be renormalizable


title: Long Journey toward the Detection of Gravitational Waves and New Era of Gravitational Wave Astrophysics

abstract: The gravitational waves were detected directly for the first time on September 14, 2015 by two LIGO detectors at Livingston, Louisiana and Hanford, Washington, USA. Careful analysis revealed that this signal was produced by the last moment of inspiral and merger of two black holes that have been orbiting each other. Since the first detection, several other gravitational wave sources, including one neutron star merger event in August 201,7 have been observed by the LIGO together with Virgo. This article provides a very brief overview of the history toward the understanding the nature of gravitational waves, detectors, types of sources and achievements of the gravitational wave detectors. We then discuss astrophysical implications of the recent detections of black hole and neutron star binaries and provide short term and long term prospects of gravitational wave astronomy.


title: All-dielectric resonant meta-optics goes active

abstract: All-dielectric resonant nanophotonics is a rapidly developing research field driven by its exceptional application potential for low-loss nanoscale metadevices. The tight confinement of the local electromagnetic fields and interferences in resonant photonic nanostructures can boost many optical effects, thus offering novel opportunities for the subwavelength control of light-matter-interactions. Active, light-emitting nanoscale structures are of particular interest, as they offer unique opportunities for novel types of light sources and nanolasers. Here, we review the latest advances in this recently emerged and rapidly growing field of active dielectric resonant nanophotonics enabled by dipolar and multipolar Mie-type resonances. More specifically, we discuss how to employ dielectric nanostructures for the resonant control of emission from quantum dots, two-dimensional transition metal dichalcogenides, and halide perovskites. We also foresee various future research directions and applications of active all-dielectric resonant nanostructures including but not limited to lasing in topologically robust systems, light-matter interactions beyond the electric dipole limit, light-emitting surfaces, and tunable active metadevices.


title: Near/Far Side Asymmetry in the Tidally Heated Moon

abstract: Using viscoelastic mass spring model simulations to track heat distribution inside a tidally perturbed body, we measure the asymmetry of heating in the crust of a spin synchronous Moon in eccentric orbit about the Earth. With the Moon within a few Earth radii of the Earth, we find that tidal heating per unit area in a lunar crustal shell is asymmetric due to the octupole order moment in the Earth's tidal field and is 10 to 20\% higher on its near side than on its far side. Tidal heating reduces the crustal basal heat flux and the rate of magma ocean crystallization. Assuming that the local crustal growth rate depends on the local basal heat flux and the distribution of tidal heating in latitude and longitude, a heat conductivity model illustrates that a moderately asymmetric and growing lunar crust could maintain its asymmetry. Near/far side asymmetric tidal heating could serve as a trigger for uneven lunar crustal growth.


title: A Jump Stochastic Differential Equation Approach for Influence Prediction on Information Propagation Networks

abstract: We propose a novel problem formulation of continuous-time information propagation on heterogenous networks based on jump stochastic differential equations (SDE). The structure of the network and activation rates between nodes are naturally taken into account in the SDE system. This new formulation allows for efficient and stable algorithm for many challenging information propagation problems, including estimations of individual activation probability and influence level, by solving the SDE numerically. To this end, we develop an efficient numerical algorithm incorporating variance reduction; furthermore, we provide theoretical bounds for its sample complexity. Moreover, we show that the proposed jump SDE approach can be applied to a much larger class of critical information propagation problems with more complicated settings. Numerical experiments on a variety of synthetic and real-world propagation networks show that the proposed method is more accurate and efficient compared with the state-of-the-art methods.


title: Calculating Time Eigenvalues of the Neutron Transport Equation with Dynamic Mode Decomposition

abstract: A novel method to compute time eigenvalues of neutron transport problems is presented based on solutions to the time dependent transport equation. Using these solutions we use the dynamic mode decomposition (DMD) to form an approximate transport operator. This approximate operator has eigenvalues that can be directly related to the time eigenvalues of the neutron transport equation. This approach works for systems of any level of criticality and does not require the user to have estimates for the eigenvalues. Numerical results are presented for homogeneous and heterogeneous media. The numerical results indicate that the method finds the eigenvalues that are most important to the solution evolution over a given time range, and the eigenvalue with the largest real part is not necessarily important to the system evolution.


title: Atomistic origin of glass-like Zn4Sb3 thermal conductivity

abstract: Zinc antimony stands out among thermoelectrics because of its very low lattice thermal conductivity, close to the amorphous limit. Understanding the physical reason behind such an unusual crystal property is of fundamental interest for the design of new thermoelectric materials. In this work we report the results of atomistic computer simulations on experimentally determined beta-Zn4Sb3 structures. We find a remarkably anharmonic behavior of Zn atoms that could be responsible for the low thermal conductivity of Zn4Sb3: their movement, better explained as diffusive, does not contribute to thermal conduction. Moreover, phonon transport is impeded by a lack of coupling between Zn and Sb atoms in crystalline positions.


title: Thermoelectric transport properties of silicon: Towards an ab initio approach

abstract: We have combined the Boltzmann transport equation with an {\it ab initio} approach to compute the thermoelectric coefficients of semiconductors. Electron-phonon, ionized impurity, and electron-plasmon scattering rates have been taken into account. The electronic band structure and average intervalley deformation potentials for the electron-phonon coupling are obtained from the density functional theory. The linearized Boltzmann equation has then been solved numerically beyond the relaxation time approximation. Our approach has been applied to crystalline silicon. We present results for the mobility, Seebeck coefficient, and electronic contribution to the thermal conductivity, as a function of the carrier concentration and temperature. The calculated coefficients are in good quantitative agreement with experimental results.


title: Absence of Casimir regime in two-dimensional nanoribbon phonon conduction

abstract: In stark contrast with three-dimensional (3D) nanostructures, we show that boundary scattering in two-dimensional (2D) nanoribbons alone does not lead to a finite phonon mean free path. If combined with an intrinsic scattering mechanism, 2D boundary scattering does reduce the overall mean free path, however the latter does not scale proportionally to the ribbon width, unlike the well known Casimir regime occurring in 3D nanowires. We show that boundary scattering can be accounted for by a simple Mathiessen type approach for many different 3D nanowire cross sectional shapes, however this is not possible in the 2D nanoribbon case, where a complete solution of the Boltzmann transport equation is required. These facts have strong implications for the thermal conductivity of suspended nanostructures.


title: Electrostatic deflections of cantilevered metallic carbon nanotubes via charge-dipole model

abstract: We compute electrostatic fields induced deformations of cantilevered finite-length metallic carbon nanotubes, using an energy minimization method based on a charge-dipole moment interaction potential combined with an empirical many-body potential. The influence of field strength, field direction and tube geometry on the electrostatic deflection is investigated for both single and double walled tubes. These results could apply to nanoelectromechanical devices based on cantilevered carbon nanotubes.


title: Lithium ion intercalation in thin crystals of hexagonal TaSe2 gated by a polymer electrolyte

abstract: Ionic liquid gating has been used to modify properties of layered transition metal dichalcogenides (TMDCs), including two-dimensional (2D) crystals of TMDCs used extensively recently in the device work, which has led to observations of properties not seen in the bulk. The main effect comes from the electrostatic gating due to strong electric field at the interface. In addition, ionic liquid gating also leads to ion intercalation when the ion size of gate electrolyte is small compared to the interlayer spacing of TMDCs. However, the microscopic processes of ion intercalation have rarely been explored in layered TMDCs. Here, we employed a technique combining photolithography device fabrication and electrical transport measurements on the thin crystals of hexagonal TaSe2 using multiple channel devices gated by a polymer electrolyte LiClO4/PEO. The gate voltage and time dependent source-drain resistances of these thin crystals were used to obtain information on the intercalation process, the effect of ion intercalation, and the correlation between the ion occupation of allowed interstitial sites and the device characteristics. We found a gate voltage controlled modulation of the charge density waves and scattering rate of charge carriers. Our work suggests that ion intercalation can be a useful tool for layered materials engineering and 2D crystal device design.


title: Electrostatic deflections of cantilevered semiconducting single-walled carbon nanotubes

abstract: How carbon nanotubes behave in an external electric field? What will be the relation between the intensity of the electric field and the tube's deformation? What are the geometry effects on the response of carbon nanotubes to electric fields? To answer these questions, we have developed a new combined computational technique to study electrostatic field induced deformations of carbon nanotubes. In this work, we find that the deflection angle of cantilevered semiconducting single-walled carbon nanotubes is proportional to the square of the electric field strength, and the tubes can be most bent when the field angle ranges from 45 to 60 degrees. Furthermore, the deflection angle is also found to be proportional to the aspect ratio L/R. Our results provide a good qualitative agreement with those of one previous experimental study.


title: Photoproduction of J/$\psi$ in non-single-diffractive p+p collisions

abstract: Recently, significant enhancements of J/$\psi$ production at very low transverse momenta were observed by the ALICE and STAR collaboration in peripheral hadronic A+A collisions. The anomaly excesses point to evidence of coherent photon-nucleus interactions in violent hadronic heavy-ion collisions, which were conventionally studied only in ultra-peripheral collisions. Assuming that the coherent photoproduction is the underlying mechanism which is responsible for the excess observed in peripheral A+A collisions, its contribution in p+p collisions with nuclear overlap, i.e. non-single-diffractive collisions, is of particular interest. In this paper, we perform a calculation of exclusive J/$\psi$ photoproduction in non-single-diffractive p+p collisions at RHIC and LHC energies base on the pQCD motivated parametrization from world-wide experimental data, which could be further employed to improve the precision of phenomenal calculations for photoproduction in A+A collisions. The differential rapidity and transverse momentum distributions of J/$\psi$ from photoproduction are presented. In comparison with the J/$\psi$ production from hadronic interactions, we find that the contribution of photoproduction is negligible.


title: Magnetic Properties of Dust Grains, Effect of Precession and Radiative Torque Alignment

abstract: Alignment of dust grains in astrophysical environments results in the polarization of starlight as well as the polarization of radiation emitted by dust. We demonstrate the advances in grain alignment theory allow the use of linear and circular polarization to probe not only the magnetic field, but also dust composition, the dust environment, etc. We revisit the process of grain alignment by Radiative Torques (RATs) and focus on constraining magnetic susceptibility of grains via observations. We discuss the possibility of observational testing of the magnetic properties of grains as the alignment changes from being in respect to the magnetic field to being in respect to the radiation direction. This opens both a possibility of constraining the uncertain parameters of the RATs theory and provides a new way of measuring magnetic fields in interstellar medium and circumstellar regions. We provide a detailed discussion of the precession induced both by the magnetic field and the anisotropic radiation and revisit a number of key processes related to magnetic response of the grains. We consider various effects that increase the rate of magnetic relaxation both in silicate and carbonaceous grains. In particular, we find a new relaxation process related to the change of the amplitude of internal magnetization within a wobbling triaxial grain and identify a range of grain sizes in which this effect can dominate the internal alignment of angular momentum within grain axes. We show that these relaxation processes significantly change the dynamics of grains in the presence of RATs. We apply our analysis for observed grain alignment in special environments to put constraints on the enhanced magnetic properties of dust grains in the cloud near supernovae, in cometary coma, and protoplanetary disks.