title: Engineering non-binary Rydberg interactions via phonons in an optical lattice

abstract: Coupling electronic and vibrational degrees of freedom of Rydberg atoms held in optical tweezer arrays offers a flexible mechanism for creating and controlling atom-atom interactions. We find that the state-dependent coupling between Rydberg atoms and local oscillator modes gives rise to two- and three-body interactions which are controllable through the strength of the local confinement. This approach even permits the cancellation of two-body terms such that three-body interactions become dominant. We analyze the structure of these interactions on two-dimensional bipartite lattice geometries and explore the impact of three-body interactions on system ground state on a square lattice. Focusing specifically on a system of $ ^{87} $Rb atoms, we show that the effects of the multi-body interactions can be maximized via a tailored dressed potential within a trapping frequency range of the order of a few hundred kHz and for temperatures corresponding to a $ >90\% $ occupation of the atomic vibrational ground state. These parameters, as well as the multi-body induced time scales, are compatible with state-of-the-art arrays of optical tweezers. Our work shows a highly versatile handle for engineering multi-body interactions of quantum many-body systems in most recent manifestations on Rydberg lattice quantum simulators.


title: Single-particle shell strengths near the doubly magic nucleus $^{56}$Ni and the $^{56}$Ni(p,$\gamma$)$^{57}$Cu reaction rate in explosive astrophysical burning

abstract: Angle-integrated cross-section measurements of the $^{56}$Ni(d,n) and (d,p) stripping reactions have been performed to determine the single-particle strengths of low-lying excited states in the mirror nuclei pair $^{57}$Cu-$^{57}$Ni situated adjacent to the doubly magic nucleus $^{56}$Ni. The reactions were studied in inverse kinematics utilizing a beam of radioactive $^{56}$Ni ions in conjunction with the GRETINA $\gamma$-array. Spectroscopic factors are compared with new shell-model calculations using a full $pf$ model space with the GPFX1A Hamiltonian for the isospin-conserving strong interaction plus Coulomb and charge-dependent Hamiltonians. These results were used to set new constraints on the $^{56}$Ni(p,$\gamma$)$^{57}$Cu reaction rate for explosive burning conditions in x-ray bursts, where $^{56}$Ni represents a key waiting point in the astrophysical rp-process.


title: New Beta Cephei stars with KELT

abstract: We present the results of a search for Galactic Beta Cephei stars, which are massive pulsating stars with both pressure modes and mixed modes. Thus, these stars can serve as benchmarks for seismological studies of the interiors of massive stars. We conducted the search by performing a frequency analysis on the optical light curves of known O- and B-type stars with data from the KELT exoplanet survey. We identify 113 Beta Cephei stars, of which 86 are new discoveries, which altogether represents a ~70% increase in the number presently known. An additional 96 candidates are identified. Among our targets, we find five new eclipsing binaries and 22 stars with equal frequency spacings suggestive of rotational splitting of nonradial pulsation modes. Candidates for runaway stars among our targets and a number of interesting individual objects are discussed. Most of the known and newly discovered Beta Cephei stars will be observed by the TESS mission, providing by far the most comprehensive observational data set of massive main sequence pulsating stars of sufficient quality for detailed asteroseismic studies. Future analysis of these light curves has the potential to dramatically increase our understanding of the structure of stellar interiors and the physical processes taking place therein.


title: WASP-South hot Jupiters: WASP-178b, WASP-184b, WASP-185b & WASP-192b

abstract: We report on four new transiting hot Jupiters discovered by the WASP-South survey. WASP-178b transits a V = 9.9, A1V star with Teff = 9350 +/- 150 K, the second-hottest transit host known. It has a highly bloated radius of 1.81 +/- 0.09 Rjup, in line with the known correlation between high irradiation and large size. With an estimated temperature of 2470 +/- 60 K, the planet is one of the best targets for studying ultra-hot Jupiters that is visible from the Southern hemisphere. The three host stars WASP-184, WASP-185 and WASP-192 are all post-main-sequence G0 stars of ages 4-8 Gyr. The larger stellar radii (1.3-1.7 Msun) mean that the transits are relatively shallow (0.7-0.9%) even though the planets have moderately inflated radii of 1.2-1.3 Rjup. WASP-185b has an eccentric orbit (e = 0.24) and a relatively long orbital period of 9.4 d. A star that is 4.6 arcsec from WASP-185 and 4.4 mag fainter might be physically associated.


title: Partition of a Subset into Two Directed Cycles with Partial Degrees

abstract: Let $D=(V,A)$ be a directed graph of order $n\geq 6$. Let $W$ be a subset of $V$ with $|W|\geq 6$. Suppose that every vertex of $W$ has degree at least $(3n-3)/2$ in $D$. Then for any integer partition $|W|=n_1+n_2$ with $n_1\geq 3$ and $n_2\geq 3$, $D$ contains two disjoint directed cycles $C_1$ and $C_2$ such that $|V(C_1)\cap W|=n_1$ and $|V(C_2)\cap W|=n_2$. We conjecture that for any integer partition $|W|=n_1+n_2+\cdots +n_k$ with $k\geq 3$ and $n_i\geq 3(1\leq i\leq k)$, $D$ contains $k$ disjoint directed cycles $C_1,C_2,\ldots , C_k$ such that $|V(C_i)\cap W|=n_i$ for all $1\leq i\leq k$. The degree condition is sharp in general.


title: Subtour Elimination Constraints Imply a Matrix-Tree Theorem SDP Constraint for the TSP

abstract: De Klerk, Pasechnik, and Sotirov give a semidefinite programming constraint for the Traveling Salesman Problem (TSP) based on the matrix-tree Theorem. This constraint says that the aggregate weight of all spanning trees in a solution to a TSP relaxation is at least that of a cycle graph. In this note, we show that the semidefinite constraint holds for any weighted 2-edge-connected graph and, in particular, is implied by the subtour elimination constraints of the subtour elimination linear program. Hence, this semidefinite constraint is implied by a finite set of linear inequality constraints.


title: Global hypoellipticity for a class of periodic Cauchy operators

abstract: This note presents an investigation on the global hypoellipticity problem for Cauchy operators on $\mathbb{T}^{n+1}$ belonging to the class \linebreak $L = \prod_{j=1}^{m}\left(D_t + c_j(t) P_j(D_x)\right)$, where $P_j(D_x)$ is a pseudo-differential operator on $\mathbb{T}^n$ and $c_j = c_j(t)$, a smooth, complex valued function on $\mathbb{T}$. The main goal of this investigation consists in establishing connections between the global hypoellipticity of the operators $L$ and its normal form $L_0 = \prod_{j=1}^m \left( D_t + c_{0,j}P_j(D_x)\right)$. In order to do so, the problem is approached by combining H\"{o}rmander's and Siegel's conditions on the symbols of the operators $L_j = D_t + c_j(t) P_j(D_x)$.


title: Probing heavy dark matter decays with multi-messenger astrophysical data

abstract: We set conservative constraints on decaying dark matter particles with masses spanning a very wide range ($10^4-10^{16}$ GeV). For this we use multimessenger observations of cosmic-ray (CR) protons/antiprotons, electrons/positrons, neutrinos/antineutrinos and gamma rays. Focusing on decays into the $\bar{b}b$ channel, we simulate the spectra of dark matter yields by using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equations and the $\texttt{Pythia}$ package. We then propagate the CRs of dark matter origin till Earth by using the state-of-the-art numerical frameworks $\texttt{CRPropa}$, $\texttt{GALPROP}$ and $\texttt{HelMod}$ for the solution of the CR transport equation in the extragalactic, Galactic region and the heliosphere, respectively. Conservative limits are obtained by requiring that the predicted dark matter spectra at Earth be less than the observed CR spectra. Overall, we exclude dark matter lifetimes of $10^{28}$ s or shorter for all the masses investigated in this work. The most stringent constraints reach $10^{30}$ s for very heavy dark matter particles with masses in the range $10^{11}-10^{14}$ GeV.


title: Fairness-Efficiency Tradeoffs in Dynamic Fair Division

abstract: We investigate the tradeoffs between fairness and efficiency when allocating indivisible items over time. Suppose T items arrive over time and must be allocated upon arrival, immediately and irrevocably, to one of n agents. Agent i assigns a value v_{it} in [0,1] to the t-th item to arrive and has an additive valuation function. If the values are chosen by an adaptive adversary it is known that the algorithm that minimizes maximum pairwise envy simply allocates each item uniformly at random; the maximum pairwise envy is then sublinear in T. If the values are independently and identically drawn from an adversarially chosen distribution D, it is also known that, under some mild conditions on D, allocating to the agent with the highest value -- a Pareto efficient allocation -- is envy-free with high probability. In this paper we study fairness-efficiency tradeoffs in this setting and provide matching upper and lower bounds under a spectrum of progressively stronger adversaries. On one hand we show that, even against a non-adaptive adversary, there is no algorithm with sublinear maximum pairwise envy that Pareto dominates the simple algorithm that allocates each item uniformly at random. On the other hand, under a slightly weaker adversary regime where item values are drawn from a known distribution and are independent with respect to time, i.e. v_{it} is independent of v_{it'} but possibly correlated with v_{i't}, optimal (in isolation) efficiency is compatible with optimal (in isolation) fairness. That is, we give an algorithm that is Pareto efficient ex-post and is simultaneously optimal with respect to fairness: for each pair of agents i and j, either i envies j by at most one item (a prominent fairness notion), or $i$ does not envy j with high probability. En route, we prove a structural result about allocations of divisible items that might be of independent interest.


title: A characterization of strong iISS for time-varying impulsive systems

abstract: For general time-varying or switched (nonlinear) systems, converse Lyapunov theorems for stability are not available. In these cases, the integral input-to-state stability (iISS) property is not equivalent to the existence of an iISS-Lyapunov function but can still be characterized as the combination of global uniform asymptotic stability under zero input (0-GUAS) and uniformly bounded energy input-bounded state (UBEBS). For impulsive systems, asymptotic stability can be weak (when the asymptotic decay depends only on elapsed time) or strong (when such a decay depends also on the number of impulses that occurred). This paper shows that the mentioned characterization of iISS remains valid for time-varying impulsive systems, provided that stability is understood in the strong sense.


title: Reducing the dependence of the neural network function to systematic uncertainties in the input space

abstract: Applications of neural networks to data analyses in natural sciences are complicated by the fact that many inputs are subject to systematic uncertainties. To control the dependence of the neural network function to variations of the input space within these systematic uncertainties, several methods have been proposed. In this work, we propose a new approach of training the neural network by introducing penalties on the variation of the neural network output directly in the loss function. This is achieved at the cost of only a small number of additional hyperparameters. It can also be pursued by treating all systematic variations in the form of statistical weights. The proposed method is demonstrated with a simple example, based on pseudo-experiments, and by a more complex example from high-energy particle physics.


title: A weak criterion of bigness for toric vector bundles

abstract: We prove a weak version of a bigness criterion for equivariant vector bundles on toric varieties.


title: Helical magnetic ordering in Sr(Co1-xNix)2As2

abstract: SrCo2As2 is a peculiar itinerant magnetic system that does not order magnetically, but inelastic neutron scattering experiments observe the same stripe-type antiferromagnetic (AF) fluctuations found in many of the Fe-based superconductors along with evidence of magnetic frustration. Here we present results from neutron diffraction measurements on single crystals of Sr(Co1-xNix)2As2 that show the development of long-range AF order with Ni-doping. However, the AF order is not stripe-type. Rather, the magnetic structure consists of ferromagnetically-aligned (FM) layers (with moments laying in the layer) that are AF arranged along c with an incommensurate propagation vector of (0 0 tau), i.e. a helix. Using high-energy x-ray diffraction, we find no evidence for a temperature-induced structural phase transition that would indicate a collinear AF order. This finding supports a picture of competing FM and AF interactions within the square transition-metal layers due to flat-band magnetic instabilities. However, the composition dependence of the propagation vector suggests that far more subtle Fermi surface and orbital effects control the interlayer magnetic correlations.


title: Doublon bottleneck in the ultrafast relaxation dynamics of hot electrons in 1T-TaS_2

abstract: Employing time-resolved photoelectron spectroscopy we analyze the relaxation dynamics of hot electrons in the charge density wave / Mott material 1T-TaS_2. At 1.2 eV above the Fermi level we observe a hot electron lifetime of 12 +- 5 fs in the metallic state and of 60 +- 10 fs in the broken symmetry ground state - a direct consequence of the reduced phase space for electron-electron scattering determined by the Mott gap. Boltzmann equation calculations which account for the interaction of hot electrons in a Bloch band with a doublon-holon excitation in the Mott state provide insight into the unoccupied electronic structure in the correlated state.


title: Massively Scaling Seismic Processing on Sunway TaihuLight Supercomputer

abstract: Common Midpoint (CMP) and Common Reflection Surface (CRS) are widely used methods for improving the signal-to-noise ratio in the field of seismic processing. These methods are computationally intensive and require high performance computing. This paper optimizes these methods on the Sunway many-core architecture and implements large-scale seismic processing on the Sunway Taihulight supercomputer. We propose the following three optimization techniques: 1) we propose a software cache method to reduce the overhead of memory accesses, and share data among CPEs via the register communication; 2) we re-design the semblance calculation procedure to further reduce the overhead of memory accesses; 3) we propose a vectorization method to improve the performance when processing the small volume of data within short loops. The experimental results show that our implementations of CMP and CRS methods on Sunway achieve 3.50x and 3.01x speedup on average compared to the-state-of-the-art implementations on CPU. In addition, our implementation is capable to run on more than one million cores of Sunway TaihuLight with good scalability.


title: Well-conditioned multiproduct Hamiltonian simulation

abstract: Product formula approximations of the time-evolution operator on quantum computers are of great interest due to their simplicity, and good scaling with system size by exploiting commutativity between Hamiltonian terms. However, product formulas exhibit poor scaling with the time $t$ and error $\epsilon$ of simulation as the gate cost of a single step scales exponentially with the order $m$ of accuracy. We introduce well-conditioned multiproduct formulas, which are a linear combination of product formulas, where a single step has polynomial cost $\mathcal{O}(m^2\log{(m)})$ and succeeds with probability $\Omega(1/\operatorname{log}^2{(m)})$. Our multiproduct formulas imply a simple and generic simulation algorithm that simultaneously exploits commutativity in arbitrary systems and has a worst-case cost $\mathcal{O}(t\log^{2}{(t/\epsilon)})$ which is optimal up to poly-logarithmic factors. In contrast, prior Trotter and post-Trotter Hamiltonian simulation algorithms realize only one of these two desirable features. A key technical result of independent interest is our solution to a conditioning problem in previous multiproduct formulas that amplified numerical errors by $e^{\Omega(m)}$ in the classical setting, and led to a vanishing success probability $e^{-\Omega(m)}$ in the quantum setting.


title: Analyzing MCMC Output

abstract: Markov chain Monte Carlo (MCMC) is a sampling-based method for estimating features of probability distributions. MCMC methods produce a serially correlated, yet representative, sample from the desired distribution. As such it can be difficult to know when the MCMC method is producing reliable results. We introduce some fundamental methods for ensuring a trustworthy simulation experiment. In particular, we present a workflow for output analysis in MCMC providing estimators, approximate sampling distributions, stopping rules, and visualization tools.


title: Halving ISO(7) supergravity

abstract: Half-maximal, $\,\mathcal{N}=4\,$, sectors of $\,D=4\,$ $\,\mathcal{N}=8\,$ supergravity with a dyonic ISO(7) gauging are investigated. We focus on a half-maximal sector including three vector multiplets, that arises as a certain $\,\textrm{SO}(3)_{\textrm{R}}$-invariant sector of the full theory. We discuss the embedding of this sector into the largest half-maximal sector of the $\,\mathcal{N}=8\,$ supergravity retaining six vector multiplets. We also provide its canonical $\,\mathcal{N}=4\,$ formulation and show that, from this perspective, our model leads in its own right to a new explicit gauging of $\,\mathcal{N}=4\,$ supergravity. Finally, expressions for the restricted duality hierarchy are given and the vacuum structure is investigated. Five new non-supersymmetric AdS vacua are found numerically. The previously known $\,\mathcal{N}=2\,$ and $\,\mathcal{N}=3\,$ AdS vacua are also contained in our $\,\mathcal{N}=4\,$ model. Unlike when embedded in previously considered sectors with fewer fields, these vacua exhibit their full $\,\mathcal{N}=2\,$ and $\,\mathcal{N}=3\,$ supersymmetry within our $\,\mathcal{N}=4\,$ model.


title: On the Surface Diffusion Flow with Triple Junctions in Higher Space Dimensions

abstract: We show short time existence for the evolution of triple junction clusters driven by the surface diffusion flow. On the triple line we use the boundary conditions derived by Garcke and Novick-Cohen as the singular limit of a Cahn-Hilliard equation with degenerated mobility. These conditions are concurrency of the triple junction, angle conditions between the hypersurfaces, continuity of the chemical potentials and a flux-balance. For the existence analysis we first write the geometric problem over a fixed reference surface and then use for the resulting analytic problem an approach in a parabolic H\"older setting.


title: Broadband Parametric Impedance Matching for Small Antennas Using the Bode-Fano Limit

abstract: In this work, a parametric up-converter amplifier is introduced as a wideband impedance matching network for receiving electrically small antennas. Chu's limit restricts the minimum Q-factor of unloaded small antennas; however, the practical bandwidth of small antennas is defined by their loaded Q-factor. By connecting a small antenna to an amplifier with a real input impedance several times greater than the radiation resistance of the antenna, we propose increasing the return loss to reduce the loaded Q-factor and increase the bandwidth as a result of the Bode-Fano theorem. In addition, a parametric amplifier is used because, in comparison with transistor amplifiers, it offers low noise characteristics. The gain of the low noise parametric amplifier compensates for the loss due to the imposed mismatch. Our simulation result shows bandwidth improvements up to 32 times can be accomplished by trading 2 dB of noise figure compare to 15 dB suggested by Chu's limit for a lossy antenna.


title: On the Design of Black-box Adversarial Examples by Leveraging Gradient-free Optimization and Operator Splitting Method

abstract: Robust machine learning is currently one of the most prominent topics which could potentially help shaping a future of advanced AI platforms that not only perform well in average cases but also in worst cases or adverse situations. Despite the long-term vision, however, existing studies on black-box adversarial attacks are still restricted to very specific settings of threat models (e.g., single distortion metric and restrictive assumption on target model's feedback to queries) and/or suffer from prohibitively high query complexity. To push for further advances in this field, we introduce a general framework based on an operator splitting method, the alternating direction method of multipliers (ADMM) to devise efficient, robust black-box attacks that work with various distortion metrics and feedback settings without incurring high query complexity. Due to the black-box nature of the threat model, the proposed ADMM solution framework is integrated with zeroth-order (ZO) optimization and Bayesian optimization (BO), and thus is applicable to the gradient-free regime. This results in two new black-box adversarial attack generation methods, ZO-ADMM and BO-ADMM. Our empirical evaluations on image classification datasets show that our proposed approaches have much lower function query complexities compared to state-of-the-art attack methods, but achieve very competitive attack success rates.


title: Synchrotron x-ray scattering of UN and U2N3 epitaxial films

abstract: We examine the magnetic ordering of UN and of a closely related nitride, U2N3, by preparing thin epitaxial films and using synchrotron x-ray techniques. The magnetic configuration and subsequent coupling to the lattice are key features of the electronic structure. The well-known antiferromagnetic (AF) ordering of UN is confirmed, but the expected accompanying distortion at Tn is not observed. Instead, we propose that the strong magneto-elastic interaction at low temperature involves changes in the strain of the material. These strains vary as a function of the sample form. As a consequence, the accepted AF configuration of UN may be incorrect. In the case of cubic a-U2N3, no single crystals have been previously prepared, and we have determined the AF ordering wave-vector. The AF Tn is close to that previously reported. In addition, resonant diffraction methods have identified an aspherical quadrupolar charge contribution in U2N3 involving the 5f electrons; the first time this has been observed in an actinide compound.


title: A Tight Degree 4 Sum-of-Squares Lower Bound for the Sherrington-Kirkpatrick Hamiltonian

abstract: We show that, if $\mathbf{W} \in \mathbb{R}^{N \times N}_{\mathsf{sym}}$ is drawn from the gaussian orthogonal ensemble, then with high probability the degree 4 sum-of-squares relaxation cannot certify an upper bound on the objective $N^{-1} \cdot \mathbf{x}^\top \mathbf{W} \mathbf{x}$ under the constraints $x_i^2 - 1 = 0$ (i.e. $\mathbf{x} \in \{ \pm 1 \}^N$) that is asymptotically smaller than $\lambda_{\max}(\mathbf{W}) \approx 2$. We also conjecture a proof technique for lower bounds against sum-of-squares relaxations of any degree held constant as $N \to \infty$, by proposing an approximate pseudomoment construction.

categories: cs.DS cond-mat.stat-mech math.OC math.PR


title: Incremental Methods for Weakly Convex Optimization

abstract: Incremental methods are widely utilized for solving finite-sum optimization problems in machine learning and signal processing. In this paper, we study a family of incremental methods -- including incremental subgradient, incremental proximal point, and incremental prox-linear methods -- for solving weakly convex optimization problems. Such a problem class covers many nonsmooth nonconvex instances that arise in engineering fields. We show that the three said incremental methods have an iteration complexity of $O(\varepsilon^{-4})$ for driving a natural stationarity measure to below $\varepsilon$. Moreover, we show that if the weakly convex function satisfies a sharpness condition, then all three incremental methods, when properly initialized and equipped with geometrically diminishing stepsizes, can achieve a local linear rate of convergence. Our work is the first to extend the convergence rate analysis of incremental methods from the nonsmooth convex regime to the weakly convex regime. Lastly, we conduct numerical experiments on the robust matrix sensing problem to illustrate the convergence performance of the three incremental methods.


title: Nonlocal magnetolectric effects in diffusive conductors with spatially inhomogeneous spin-orbit coupling

abstract: We present a theoretical study of nonlocal magnetoelectric effects in diffusive hybrid structures with an intrinsic linear-in-momentum spin-orbit coupling (SOC) which is assumed to be spatially inhomogeneous. Our analysis is based on the SU(2)-covariant drift-diffusion equations from which we derive the BC at hybrid interfaces. Within this formulation, the spin current is covariantly conserved when the spin relaxation is only due to the intrinsic SOC. This conservation leads to the absence of spin Hall (SH) currents in homogeneous systems. We also consider extrinsic sources of relaxation (ESR), as magnetic impurities, which break the covariant spin conservation, and may lead to SH currents. We apply our model to describe nonlocal transport in a system with an interface separating two regions: one normal region without intrinsic SOC and one with a Rashba SOC. We first explore the inverse spin-galvanic effect, i.e., a spin polarization induced by an electric field. We demonstrate how the spatial behavior of such spin density depends on both, the direction of the electric field and the strength of the ESR rate. We also study the spin-to-charge conversion, and compute the charge current and the distribution of electrochemical potential in the Rashba region induced by a spin current injected into the normal region. In systems with an inhomogeneous SOC varying in one spatial direction, we find an interesting nonlocal reciprocity between the spin density induced by a charge current at a given point in space, and the spatially integrated current induced by a spin density injected at the same point.


title: Magnetoelectric control of topological phases in graphene

abstract: Topological antiferromagnetic (AFM) spintronics is an emerging field of research, which involves the topological electronic states coupled to the AFM order parameter known as the N$\acute{\rm e}$el vector. The control of these states is envisioned through manipulation of the N$\acute{\rm e}$el vector by spin-orbit torques driven by electric currents. Here we propose a different approach favorable for low-power AFM spintronics, where the control of the topological states in a two-dimensional material, such as graphene, is performed via the proximity effect by the voltage induced switching of the N$\acute{\rm e}$el vector in an adjacent magnetoelectric AFM insulator, such as chromia. Mediated by the symmetry protected boundary magnetization and the induced Rashba-type spin-orbit coupling at the interface between graphene and chromia, the emergent topological phases in graphene can be controlled by the N$\acute{\rm e}$el vector. Using density functional theory and tight-binding Hamiltonian approaches, we model a graphene/Cr2O3 (0001) interface and demonstrate non-trivial band gap openings in the graphene Dirac bands asymmetric between the K and K' valleys. This gives rise to an unconventional quantum anomalous Hall effect (QAHE) with a quantized value of $2e^2/h$ and an additional step-like feature at a value close to $e^2/2h$, and the emergence of the spin-polarized valley Hall effect (VHE). Furthermore, depending on the N$\acute{\rm e}$el vector orientation, we predict the appearance and transformation of different topological phases in graphene across the $180^{\circ}$ AFM domain wall, involving the QAHE, the valley-polarized QAHE and the quantum VHE (QVHE), and the emergence of the chiral edge state along the domain wall. These topological properties are controlled by voltage through magnetoelectric switching of the AFM insulator with no need for spin-orbit torques.


title: A new class of solutions to Laplace equation: Regularized multipoles of negative orders

abstract: We introduce a new class of solutions to Laplace equation, dubbed logopoles, and use them to derive a new relation between solutions in prolate spheroidal and spherical coordinates. The main novelty is that it involves spherical harmonics of the second kind, which have rarely been considered in physical problems because they are singular on the entire z axis. Logopoles, in contrast, have a finite line singularity like solid spheroidal harmonics, but are also closely related to solid spherical harmonics and can be viewed as an extension of the standard multipole ladder toward the negative multipolar orders. These new solutions may prove a fruitful alternative to either spherical or spheroidal harmonics in physical problems.

categories: math-ph math.MP physics.class-ph physics.comp-ph


title: Flocking in complex environments -- attention trade-offs in collective information processing

abstract: The ability of biological and artificial collectives to outperform solitary individuals in a wide variety of tasks depends crucially on the efficient processing of social and environmental information at the level of the collective. Here, we model collective behavior in complex environments with many potentially distracting cues. Counter-intuitively, large-scale coordination in such environments can be maximized by strongly limiting the cognitive capacity of individuals, where due to self-organized dynamics the collective self-isolates from disrupting information. We observe a fundamental trade-off between coordination and collective responsiveness to environmental cues. Our results offer important insights into possible evolutionary trade-offs in collective behavior in biology and suggests novel principles for design of artificial swarms exploiting attentional bottlenecks.


title: RoBERTa: A Robustly Optimized BERT Pretraining Approach

abstract: Language model pretraining has led to significant performance gains but careful comparison between different approaches is challenging. Training is computationally expensive, often done on private datasets of different sizes, and, as we will show, hyperparameter choices have significant impact on the final results. We present a replication study of BERT pretraining (Devlin et al., 2019) that carefully measures the impact of many key hyperparameters and training data size. We find that BERT was significantly undertrained, and can match or exceed the performance of every model published after it. Our best model achieves state-of-the-art results on GLUE, RACE and SQuAD. These results highlight the importance of previously overlooked design choices, and raise questions about the source of recently reported improvements. We release our models and code.


title: Nonthermal ion acceleration by the kink instability in nonrelativistic jets

abstract: We investigate the self-consistent particle acceleration physics associated with the development of the kink instability (KI) in nonrelativistic, electron-ion plasma jets. Using 3D fully kinetic particle-in-cell (PIC) simulations, we show that the KI efficiently converts the initial toroidal magnetic field energy into energetic ions. The accelerated ions form a nonthermal power-law tail in the energy spectrum, containing $\simeq10\%$ of the initial magnetic field energy, and with the maximum ion energy extending to the confinement energy of the jet. We find that the ions are efficiently accelerated by the concerted action of the motional electric field and highly tangled magnetic field that develop in the nonlinear phase of the KI: fast curvature drift motions of ions across magnetic field lines enable their acceleration along the electric field. We further investigate the role of Coulomb collisions on the ion acceleration efficiency, and identify the collisional threshold above which nonthermal ion acceleration is suppressed. Our results reveal how energetic ions may result from unstable nonrelativistic plasma jets in space and astrophysics, and provide constraints on the plasma conditions required to reproduce this acceleration mechanism in laboratory experiments.


title: On the Stochastic Flows on $(m+n+1)$-Dimensional Exotic Spheres

abstract: Stochastic flows of Stratonovich stochastic differential equations on exotic spheres have been studied. The consequences of the choice of exotic differential structure on stochastic processes taking place on the topological space $S^{m+n+1}$ as state space of the processes have been investigated. More precisely, we have investigated the properties of stochastic processes where the state spaces of the stochastic processes under consideration are $({m+n+1})$-dimensional differentiable manifolds which are homeomorphic but not necessarily diffeomorphic to standard ${(m+n+1)}$-dimensional sphere. The differentiable manifolds have been constructed from disjoint union $\mathbb{R}^{m+1}\times S^{n}\sqcup S^m\times \mathbb{R}^{n+1}$ by identifying every pair of its points using a map $u :\mathbb{R}^{m+1}\times S^n\rightarrow S^m\times \mathbb{R}^{n+1}$ which is constructed from a diffeomorphism $h_1\times h_2:S^m\times S^n\rightarrow S^m\times S^n$. The diffeomorphisms $h_1$ and $h_2$, therefore, can be regarded as the carriers of the "exoticism" of the constructed manifolds. For all of the above purposes, homeomorphisms $h$ from the above-constructed manifolds onto the standard sphere explicitly in term of the diffeomorphisms $h_1$ and $h_2$ have been constructed. Using the homeomorphisms $h$ and all their associated maps derived from them and expressed in terms of $h_1$ and $h_2$ as well as their derivatives, we construct the stochastic processes or flows on the above-constructed manifolds corresponding to stochastics processes on the standard sphere $S^{m+n+1}_s$. The stochastic processes yielded from the above construction on the constructed manifolds can be regarded as the same stochastic processes on $S^{m+n+1}_s$ but described in exotic differential structures on $S^{m+n+1}$.

categories: math-ph math.AP math.AT math.DG math.MP math.PR


title: Time-dependent coupled cluster theory on the Keldysh contour for non-equilibrium systems

abstract: We leverage the Keldysh formalism to extend our implementation of finite temperature coupled cluster theory [\textit{J. Chem. Theory Comput.} 2018, \textit{14}, 5690-5700] to thermal systems that have been driven out of equilibrium. The resulting Keldysh coupled cluster theory is discussed in detail. We describe the implementation of the equations necessary to perform Keldysh coupled cluster singles and doubles calculations of finite temperature dynamics, and we apply the method to some simple systems including a Hubbard model with a Peierls phase and an {\it ab initio} model of warm-dense silicon subject to an ultrafact XUV pulse.


title: Hidden Treasures: sterile neutrinos as dark matter with miraculous abundance, structure formation for different production mechanisms, and a solution to the sigma-8 problem

abstract: We discuss numerous mechanisms for production of sterile neutrinos, which can account for all or a fraction of dark matter, and which can range from warm to effectively cold dark matter, depending on the cosmological scenario. We investigate production by Higgs boson decay, $(B-L)$ gauge boson production at high temperature, as well as production via resonant and nonresonant neutrino oscillations. We calculate the effects on structure formation in these models, some for the first time. If two populations of sterile neutrinos, one warm and one cold, were produced by different mechanisms, or if sterile neutrinos account for only a fraction of dark matter, while the remainder is some other cold dark matter particle, the resulting multi-component dark matter may alleviate some problems in galaxy formation. We examine the X-ray constraints and the candidate signal at 3.5 keV. Finally, we also show that the $\sigma_8$ problem can be a signature of fractional dark matter in the form of sterile neutrinos in several mechanisms.


title: Rotor imbalance suppression by optimal control

abstract: An imbalanced rotor is considered. A system of moving balancing masses is given. We determine the optimal movement of the balancing masses to minimize the imbalance on the rotor. The optimal movement is given by an open-loop control solving an optimal control problem posed in infinite time. By methods of the Calculus of Variations, the existence of the optimum is proved and the corresponding optimality conditions have been derived. Asymptotic behavior of the control system is studied rigorously. By \L ojasiewicz inequality, convergence of the optima as time $t\to +\infty$ towards a steady configuration is ensured. An explicit estimate of the convergence rate is given. This guarantees that the optimal control stabilizes the system. In case the imbalance is below a computed threshold, the convergence occurs exponentially fast. This is proved by the Stable Manifold Theorem applied to the Pontryagin optimality system. Moreover, a closed-loop control strategy based on Reinforcement Learning is proposed. Numerical simulations have been performed, validating the theoretical results.


title: Electrical switching in a magnetically intercalated transition metal dichalcogenide

abstract: Recent advances in tuning the correlated behavior of graphene and transition-metal dichalcogenides (TMDs) have opened a new frontier in the study of many-body physics in two dimensions and promise exciting possibilities for new quantum technologies. An emerging field where these materials have yet to make a deep impact is the study of antiferromagnetic (AFM) spintronics - a relatively new research direction that promises technologies that are insensitive to external magnetic fields, fast switching times, and reduced crosstalk. In this study we present measurements on the intercalated TMD Fe1/3NbS2 which exhibits antiferromagnetic ordering below 42K. We find that current densities on the order of 10^4 A/cm^2 can reorient the magnetic order, the response of which can be detected in the sample's resistance. This demonstrates that Fe1/3NbS2 can be used as an antiferromagnetic switch with electronic "write-in" and "read-out". This switching is found to be stable over time and remarkably robust to external magnetic fields. Fe1/3NbS2 is a rare example of an AFM system that exhibits fully electronic switching behavior in single crystal form, making it appealing for low-power, low-temperature memory storage applications. Moreover, Fe1/3NbS2 is part of a much larger family of magnetically intercalated TMDs, some of which may exhibit the switching behavior at higher temperatures and form a platform from which to build tunable AFM spintronic devices.


title: The IceCube Neutrino Observatory -- Contributions to the 36th International Cosmic Ray Conference (ICRC2019)

abstract: Contributions from the IceCube Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA.


title: Investing for Discovery in Astronomy

abstract: How should we invest our available resources to best sustain astronomy's track record of discovery, established over the past few decades? Two strong hints come from (1) our history of astronomical discoveries and (2) literature citation patterns that reveal how discovery and development activities in science are strong functions of team size. These argue that progress in astronomy hinges on support for a diversity of research efforts in terms of team size, research tools and platforms, and investment strategies that encourage risk taking. These ideas also encourage us to examine the implications of the trend toward "big team science" and "survey science" in astronomy over the past few decades, and to reconsider the common assumption that progress in astronomy always means "trading up" to bigger apertures and facilities. Instead, the considerations above argue that we need a balanced set of investments in small- to large-scale initiatives and team sizes both large and small. Large teams tend to develop existing ideas, whereas small teams are more likely to fuel the future with disruptive discoveries. While large facilities are the "value" investments that are guaranteed to produce discoveries, smaller facilities are the "growth stocks" that are likely to deliver the biggest science bang per buck, sometimes with outsize returns. One way to foster the risk taking that fuels discovery is to increase observing opportunity, i.e., create more observing nights and facilitate the exploration of science-ready data.

categories: astro-ph.IM astro-ph.CO astro-ph.GA astro-ph.SR


title: Solving the one-dimensional Ising chain via mathematical induction: An intuitive approach to the transfer matrix

abstract: The aim of this work is to present a formulation to solve the one-dimensional Ising model using the elementary technique of mathematical induction. This formulation is physically clear and leads to the same partition function form as the transfer matrix method, which is a common subject in the introductory courses of statistical mechanics. In this way our formulation is a useful tool to complement the traditional more abstract transfer matrix method. The method can be straightforwardly generalized to other short-range chains, coupled chains and is also computationally friendly. These two approaches provide a more complete understanding of the system, and therefore our work can be of broad interest for undergraduate teaching in statistical mechanics.


title: Diffusion Phenomena in a Mixed Phase Space

abstract: We show that, in strongly chaotic dynamical systems, the average particle velocity can be calculated analytically by consideration of Brownian dynamics in phase space, the method of images and use of the classical diffusion equation. The method is demonstrated on the simplified Fermi-Ulam accelerator model, which has a mixed phase space with chaotic seas, invariant tori and Kolmogorov-Arnold-Moser (KAM) islands. The calculated average velocities agree well with numerical simulations and with an earlier empirical theory. The procedure can readily be extended to other systems including time-dependent billiards.


title: Action Guidance with MCTS for Deep Reinforcement Learning

abstract: Deep reinforcement learning has achieved great successes in recent years, however, one main challenge is the sample inefficiency. In this paper, we focus on how to use action guidance by means of a non-expert demonstrator to improve sample efficiency in a domain with sparse, delayed, and possibly deceptive rewards: the recently-proposed multi-agent benchmark of Pommerman. We propose a new framework where even a non-expert simulated demonstrator, e.g., planning algorithms such as Monte Carlo tree search with a small number rollouts, can be integrated within asynchronous distributed deep reinforcement learning methods. Compared to a vanilla deep RL algorithm, our proposed methods both learn faster and converge to better policies on a two-player mini version of the Pommerman game.


title: Accurate and Robust Pulmonary Nodule Detection by 3D Feature Pyramid Network with Self-supervised Feature Learning

abstract: Accurate detection of pulmonary nodules with high sensitivity and specificity is essential for automatic lung cancer diagnosis from CT scans. Although many deep learning-based algorithms make great progress for improving the accuracy of nodule detection, the high false positive rate is still a challenging problem which limits the automatic diagnosis in routine clinical practice. Moreover, the CT scans collected from multiple manufacturers may affect the robustness of Computer-aided diagnosis (CAD) due to the differences in intensity scales and machine noises. In this paper, we propose a novel self-supervised learning assisted pulmonary nodule detection framework based on a 3D Feature Pyramid Network (3DFPN) to improve the sensitivity of nodule detection by employing multi-scale features to increase the resolution of nodules, as well as a parallel top-down path to transit the high-level semantic features to complement low-level general features. Furthermore, a High Sensitivity and Specificity (HS2) network is introduced to eliminate the false positive nodule candidates by tracking the appearance changes in continuous CT slices of each nodule candidate on Location History Images (LHI). In addition, in order to improve the performance consistency of the proposed framework across data captured by different CT scanners without using additional annotations, an effective self-supervised learning schema is applied to learn spatiotemporal features of CT scans from large-scale unlabeled data. The performance and robustness of our method are evaluated on several publicly available datasets with significant performance improvements. The proposed framework is able to accurately detect pulmonary nodules with high sensitivity and specificity and achieves 90.6% sensitivity with 1/8 false positive per scan which outperforms the state-of-the-art results 15.8% on LUNA16 dataset.


title: Low-Rank Matrix Completion: A Contemporary Survey

abstract: As a paradigm to recover unknown entries of a matrix from partial observations, low-rank matrix completion (LRMC) has generated a great deal of interest. Over the years, there have been lots of works on this topic but it might not be easy to grasp the essential knowledge from these studies. This is mainly because many of these works are highly theoretical or a proposal of new LRMC technique. In this paper, we give a contemporary survey on LRMC. In order to provide better view, insight, and understanding of potentials and limitations of LRMC, we present early scattered results in a structured and accessible way. Specifically, we classify the state-of-the-art LRMC techniques into two main categories and then explain each category in detail. We next discuss issues to be considered when one considers using LRMC techniques. These include intrinsic properties required for the matrix recovery and how to exploit a special structure in LRMC design. We also discuss the convolutional neural network (CNN) based LRMC algorithms exploiting the graph structure of a low-rank matrix. Further, we present the recovery performance and the computational complexity of the state-of-the-art LRMC techniques. Our hope is that this survey article will serve as a useful guide for practitioners and non-experts to catch the gist of LRMC.


title: Vortex phase diagram and the normal state of cuprates with charge and spin orders

abstract: The phase diagram of underdoped cuprates in a magnetic field ($H$) is the key ingredient in understanding the anomalous normal state of these high-temperature superconductors. However, the upper critical field ($H_{c2}$) or the extent of superconducting phase with vortices, a type of topological excitations, and the role of charge orders that are present at high $H$, remain under debate. We address these questions by studying stripe-ordered La-214, i.e. cuprates in which charge orders are most pronounced and zero-field transition temperatures $T_{c}^{0}$ are lowest; the latter opens a much larger energy scale window to explore the vortex phases compared to previous studies. By combining linear and nonlinear transport techniques sensitive to vortex matter, we determine the $T$-$H$ phase diagram, directly detect $H_{c2}$, and reveal novel properties of the high-field ground state. Our results demonstrate that, while the vortex phase diagram of underdoped cuprates is not very sensitive to the details of the charge orders, quantum fluctuations and disorder play a key role as $T\rightarrow 0$. The presence of stripes, on the other hand, seems to alter the nature of the anomalous normal state, such that the high-field ground state is a metal, as opposed to an insulator.


title: On the difficulty of proving P equals NP in ZFC

abstract: Harvey Friedman, in his remarkable paper Finite functions and the necessary use of large cardinals, Ann. Math. 148:803-893, 1998 and in a technical report, Applications of large cardinals to graph theory, Ohio State University, 1997, presents numerous combinatorial statements with clear geometric meaning that are proved using large cardinals and shown to require them. By slightly extending some of Friedman's results, we construct an infinite class of structurally similar theorems which can be proved using the same large cardinals but also can be proved using the statement "subset sum is solvable in polynomial time." This curious connection between the P vs. NP problem and the theory of large cardinals seems to suggest that either P=NP is false or otherwise not provable in ZFC.


title: Signatures of a pair density wave at high magnetic fields in cuprates with charge and spin orders

abstract: In underdoped cuprates, the interplay of the pseudogap, superconductivity, and charge and spin ordering can give rise to exotic quantum states, including the pair density wave (PDW), in which the superconducting (SC) order parameter is oscillatory in space. However, the evidence for a PDW state remains inconclusive and its broader relevance to cuprate physics is an open question. To test the interlayer frustration, the crucial component of the PDW picture, we performed transport measurements on La$_{1.7}$Eu$_{0.2}$Sr$_{0.1}$CuO$_{4}$ and La$_{1.48}$Nd$_{0.4}$Sr$_{0.12}$CuO$_{4}$, cuprates with "striped" spin and charge orders, in perpendicular magnetic fields ($H_\perp$), and also with an additional field applied parallel to CuO$_2$ layers ($H_\parallel$). We detected several phenomena predicted to arise from the existence of a PDW, including an enhancement of interlayer SC phase coherence with increasing $H_\parallel$. Our findings are consistent with the presence of local, PDW pairing correlations that compete with the uniform SC order at $T_{c}^{0}< T<(2-6) T_{c}^{0}$, where $T_{c}^{0}$ is the $H=0$ SC transition temperature, and become dominant at intermediate $H_\perp$ as $T\rightarrow 0$. These data also provide much-needed transport signatures of the PDW in the regime where superconductivity is destroyed by quantum phase fluctuations.


title: Bernstein-Sato roots for monomial ideals in positive characteristic

abstract: Following work of Musta\c{t}\u{a} and Bitoun we recently developed a notion of Bernstein-Sato roots for arbitrary ideals, which is a prime characteristic analogue for the roots of the Bernstein-Sato polynomial. Here we prove that for monomial ideals the roots of the Bernstein-Sato polynomial (over $\mathbb{C}$) agree with the Bernstein-Sato roots of the mod-$p$ reductions of the ideal for $p$ large enough. We regard this as evidence that the characteristic-$p$ notion of Bernstein-Sato root is reasonable.


title: Attack Synthesis for Strings using Meta-Heuristics

abstract: Information leaks are a significant problem in modern computer systems and string manipulation is prevalent in modern software. We present techniques for automated synthesis of side-channel attacks that recover secret string values based on timing observations on string manipulating code. Our attack synthesis techniques iteratively generate inputs which, when fed to code that accesses the secret, reveal partial information about the secret based on the timing observations, leading to recovery of the secret at the end of the attack sequence. We use symbolic execution to extract path constraints, automata-based model counting to estimate the probability of execution paths, and meta-heuristic methods to maximize information gain based on entropy for synthesizing adaptive attack steps.


title: Deep MRI Reconstruction: Unrolled Optimization Algorithms Meet Neural Networks

abstract: Image reconstruction from undersampled k-space data has been playing an important role for fast MRI. Recently, deep learning has demonstrated tremendous success in various fields and also shown potential to significantly speed up MR reconstruction with reduced measurements. This article gives an overview of deep learning-based image reconstruction methods for MRI. Three types of deep learning-based approaches are reviewed, the data-driven, model-driven and integrated approaches. The main structure of each network in three approaches is explained and the analysis of common parts of reviewed networks and differences in-between are highlighted. Based on the review, a number of signal processing issues are discussed for maximizing the potential of deep reconstruction for fast MRI. the discussion may facilitate further development of "optimal" network and performance analysis from a theoretical point of view.

categories: eess.IV cs.CV cs.LG eess.SP physics.med-ph stat.ML


title: L^p-asymptotic stability analysis of a 1D wave equation with a nonlinear damping

abstract: This paper is concerned with the asymptotic stability analysis of a one dimensional wave equation with Dirichlet boundary conditions subject to a nonlinear distributed damping with an L p functional framework, p $\in$ [2, $\infty$]. Some well-posedness results are provided together with exponential decay to zero of trajectories, with an estimation of the decay rate. The well-posedness results are proved by considering an appropriate functional of the energy in the desired functional spaces introduced by Haraux in [11]. Asymptotic behavior analysis is based on an attractivity result on a trajectory of an infinite-dimensional linear time-varying system with a special structure, which relies on the introduction of a suitable Lyapunov functional. Note that some of the results of this paper apply for a large class of nonmonotone dampings.


title: Learning to Synthesize: Robust Phase Retrieval at Low Photon counts

abstract: The quality of inverse problem solutions obtained through deep learning [Barbastathis et al, 2019] is limited by the nature of the priors learned from examples presented during the training phase. In the case of quantitative phase retrieval [Sinha et al, 2017, Goy et al, 2019], in particular, spatial frequencies that are underrepresented in the training database, most often at the high band, tend to be suppressed in the reconstruction. Ad hoc solutions have been proposed, such as pre-amplifying the high spatial frequencies in the examples [Li et al, 2018]; however, while that strategy improves resolution, it also leads to high-frequency artifacts as well as low-frequency distortions in the reconstructions. Here, we present a new approach that learns separately how to handle the two frequency bands, low and high; and also learns how to synthesize these two bands into the full-band reconstructions. We show that this "learning to synthesize" (LS) method yields phase reconstructions of high spatial resolution and artifact-free; and it is also resilient to high-noise conditions, e.g. in the case of very low photon flux. In addition to the problem of quantitative phase retrieval, the LS method is applicable, in principle, to any inverse problem where the forward operator treats different frequency bands unevenly, i.e. is ill-posed.


title: Modular $A_4$ Symmetry Models of Neutrinos and Charged Leptons

abstract: We present a comprehensive analysis of neutrino mass and lepton mixing in theories with $A_4$ modular symmetry, where the only flavon field is the single modulus field $\tau$, and all masses and Yukawa couplings are modular forms. Similar to previous analyses, we discuss all the simplest neutrino sectors arising from both the Weinberg operator and the type I seesaw mechanism, with lepton doublets and right-handed neutrinos assumed to be triplets of $A_4$. Unlike previous analyses, we allow right-handed charged leptons to transform as all combinations of $\mathbf{1}$, $\mathbf{1}'$ and $\mathbf{1}''$ representations of $A_4$, using the simplest different modular weights to break the degeneracy, leading to ten different charged lepton Yukawa matrices, instead of the usual one. This implies ten different Weinberg models and thirty different type I seesaw models, which we analyse in detail. We find that fourteen models for both NO and IO neutrino mass ordering can accommodate the data, as compared to one in previous analyses, providing many new possibilities.


title: Tunable magneto-optical properties of single-layer tin diselenide: From GW approximation to large-scale tight-binding calculations

abstract: A parameterized tight-binding (TB) model based on the first-principles GW calculations is developed for single layer tin diselenide (SnSe$_2$) and used to study its electronic and optical properties under external magnetic field. The truncated model is derived from six maximally localized wannier orbitals on Se site, which accurately describes the quasi-particle electronic states of single layer SnSe$_2$ in a wide energy range. The quasi-particle electronic states are dominated by the hoppings between nearest wannier orbitals ($t_1$-$t_6$). Our numerical calculation shows that, due to the electron-hole asymmetry, two sets of Landau Level spectrum are obtained when a perpendicular magnetic field is applied. The Landau Level spectrum follows linear dependence on the level index and magnetic field, exhibiting properties of two-dimensional electron gas in traditional semiconductors. The optical conductivity calculation shows that the optical gap is very close to the GW value, and can be tuned by external magnetic field. Our proposed TB model can be used for further exploring the electronic, optical, and transport properties of SnSe$_2$, especially in the presence of external magnetic fields.


title: Quasicrystal vs Glass Transition: comparison structural and dynamical properties

abstract: Quasicrystals are solid structures with symmetry forbidden by crystallographic rules. Because of this some structural characteristics of quasicrystals, for instance, radial distribution function, can look similar to the ones of amorphous phases. This is of principal importance since radial distribution function is the main property to characterize the structure in molecular simulation. In the present paper we compare the radial distribution functions and dynamical properties of three systems in the vicinity of glass transition, quasicrystal formation and crystallization. We show that in spite of similarity of radial distribution functions the dynamical properties of a system in the vicinity of quasicrystal are qualitatively equivalent to the ones of crystal. Because of this combination the radial distribution functions with investigation of dynamics of the liquid allows unambiguously distinguish glass and quasicrystal.


title: Secure Distribution of Protected Content in Information-Centric Networking

abstract: The benefits of the ubiquitous caching in ICN are profound, such features make ICN promising for content distribution, but it also introduces a challenge to content protection against the unauthorized access. The protection of a content against unauthorized access requires consumer authentication and involves the conventional end-to-end encryption. However, in information-centric networking (ICN), such end-to-end encryption makes the content caching ineffective since encrypted contents stored in a cache are useless for any consumers except those who know the encryption key. For effective caching of encrypted contents in ICN, we propose a secure distribution of protected content (SDPC) scheme, which ensures that only authenticated consumers can access the content. SDPC is lightweight and allows consumers to verify the originality of the published content by using a symmetric key encryption. Moreover, SDPC naming scheme provides protection against privacy leakage. The security of SDPC was proved with the BAN logic and Scyther tool verification, and simulation results show that SDPC can reduce the content download delay.


title: Large scale continuous-time mean-variance portfolio allocation via reinforcement learning

abstract: We propose to solve large scale Markowitz mean-variance (MV) portfolio allocation problem using reinforcement learning (RL). By adopting the recently developed continuous-time exploratory control framework, we formulate the exploratory MV problem in high dimensions. We further show the optimality of a multivariate Gaussian feedback policy, with time-decaying variance, in trading off exploration and exploitation. Based on a provable policy improvement theorem, we devise a scalable and data-efficient RL algorithm and conduct large scale empirical tests using data from the S&P 500 stocks. We found that our method consistently achieves over 10% annualized returns and it outperforms econometric methods and the deep RL method by large margins, for both long and medium terms of investment with monthly and daily trading.


title: Market and Long Term Accounting Operational Performance

abstract: Following the value relevance literature, this study verifies whether the marketplace differentiates companies of high, medium, and low long-term operational performance, measured by accounting information on profitability, sales variation and indebtedness. The data comprises the Corporate Financial Statements disclosed during the period from 1996 to 2009 and stock prices of companies listed on the Sao Paulo Stock Exchange and Commodities and Futures Exchange - BM&FBOVESPA. The final sample is composed of 142 non-financial companies. Five year mobile windows were used, which resulted in ten five-year periods. After checking each company indices, the accounting variables were unified in an Index Performance Summary to synthesize the final performance for each five-year period, which allowed segregation in operational performance levels. Multiple regressions were performed using panel data techniques, fixed effects model and dummies variables, and then hypothesis tests were made. Regarding the explanatory power of each individual variable, the results show that not all behaviors are according to the research hypothesis and that the Brazilian stock market differentiates companies of high and low long-term operational performance. This distinction is not fully perceived between companies of high and medium operational performance.


title: OH maser emission in the THOR survey of the northern Milky Way

abstract: Context: OH masers trace diverse physical processes, from the expanding envelopes around evolved stars to star-forming regions or supernovae remnants. Aims: We identify the ground-state OH masers at 18cm wavelength in the area covered by ``The HI/OH/Recombination line survey of the Milky Way (THOR)''. We present a catalogue of all OH maser features and their possible associated environments. Methods: The THOR survey covers longitude and latitude ranges of 14.3<l<66.8 and b<1.25 deg. All OH ground state lines at 1612, 1665, 1667 and 1720MHz have been observed, employing the Very Large Array. The spatial resolution of the data varies between 12.5'' and 19'', the spectral resolution is 1.5km/s, and the rms sensitivity of the data is ~10mJy/beam per channel. Results: We identify 1585 individual maser spots distributed over 807 maser sites. Based on different criteria from spectral profiles to literature comparison, we try to associate the maser sites with astrophysical source types. Approximately 51\% of the sites exhibit the double-horned 1612MHz spectra typically emitted from the expanding shells of evolved stars. The separations of the two main velocity features of the expanding shells typically vary between 22 and 38km/s. In addition to this, at least 20% of the maser sites are associated with star-forming regions. While the largest fraction of 1720MHz maser spots (21 out of 53) is associated with supernova remnants, a significant fraction of the 1720MHz maser spots (17) are also associated with star-forming regions. We present comparisons to the thermal 13CO(1-0) emission as well as to other surveys of class II CH3OH and H2O maser emission. The catalogue attempts to present associations to astrophysical sources where available, and the full catalogue is available in electronic form.


title: FAKIR: An algorithm for revealing the anatomy and pose of statues from raw point sets

abstract: 3D acquisition of archaeological artefacts has become an essential part of cultural heritage research for preservation or restoration purpose. Statues, in particular, have been at the center of many projects. In this paper, we introduce a way to improve the understanding of acquired statues representing real or imaginary creatures by registering a simple and pliable articulated model to the raw point set data. Our approach performs a Forward And bacKward Iterative Registration (FAKIR) which proceeds joint by joint, needing only a few iterations to converge. We are thus able to detect the pose and elementary anatomy of sculptures, with possibly non realistic body proportions. By adapting our simple skeleton, our method can work on animals and imaginary creatures.


title: Correlation functions in ${\cal N}=2$ Supersymmetric vector matter Chern-Simons theory

abstract: We compute the two, three point function of the opearators in the spin zero multiplet of ${\cal N}=2$ Supersymmetric vector matter Chern-Simons theory at large $N$ and at all orders of 't Hooft coupling by solving the Schwinger-Dyson equation. Schwinger-Dyson method to compute four point function becomes extremely complicated and hence we use bootstrap method to solve for four point function of scaler operator $J_0^{f}=\bar\psi \psi$ and $J_0^{b}=\bar\phi \phi$. Interestingly, due to the fact that $\langle J_0^{f}J_0^{f}J_0^{b} \rangle$ is a contact term, the four point function of $ J_0^{f}$ operator looks like that of free theory up to overall coupling constant dependent factors and up to some bulk AdS contact terms. On the other hand the $J_0^{b}$ four-point function receives an additional contribution compared to the free theory expression due to the $J_0^{f}$ exchange. Interestingly, double discontinuity of this single trace operator $J_0^{f}$ vanishes and hence it only contributes to AdS-contact term.


title: Ferromagnetism in narrow bands of moir\'e superlattices

abstract: Many graphene moir\'e superlattices host narrow bands with non-zero valley Chern numbers. We provide analytical and numerical evidence for a robust spin and/or valley polarized insulator at total integer band filling in nearly flat bands of several different moir\'e materials. In the limit of a perfectly flat band, we present analytical arguments in favor of the ferromagnetic state substantiated by numerical calculations. Further, we numerically evaluate its stability for a finite bandwidth. We provide exact diagonalization results for models appropriate for ABC trilayer graphene aligned with hBN, twisted double bilayer graphene, and twisted bilayer graphene aligned with hBN. We also provide DMRG results for a honeycomb lattice with a quasi-flat band and non-zero Chern number, which extend our results to larger system sizes. We find a maximally spin and valley polarized insulator at all integer fillings when the band is sufficiently flat. We also show that interactions may induce effective dispersive terms strong enough to destabilize this state. These results still hold in the case of zero valley Chern number (for example, trivial side of TLG/hBN). We give an intuitive picture based on extended Wannier orbitals, and emphasize the role of the quantum geometry of the band, whose microscopic details may enhance or weaken ferromagnetism in moir\'e materials.


title: Fragmentation in Population III Galaxies formed through Ionizing Radiation

abstract: Population III stars forming in minihalos tend to be relatively inefficient, with each minihalo hosting one or a small number of stars which are more massive than local stars but still challenging to observe directly at high redshift. Here we explore a possible mechanisms for the generation of larger clusters of such stars: a nearby ionizing source which ionizes a late forming halo, delaying its collapse until the halo is sufficiently large that the core can self-shield and suffer runaway collapse. We use simulations with a simple but accurate model for the radiative ionizing flux and confirm the basic predictions of previous work: higher ionizing fluxes can delay the collapse to lower redshifts and higher masses, up to an order of magnitude above the atomic cooling limit. In a limited number of runs we also examine the fragmentation of the cores at even higher resolution, using both simple estimates and sink particles to show that the number of fragments is generally small, at most a handful, and that the mass accretion rate on the fragments is of order $10^{-3}$ M$_\odot$ yr$^{-1}$. This rate is sufficiently high that the descent on the main sequence (and hence the suppression of accretion) is delayed until the stellar masses are of order $100-1000$ M$_\odot$, but not high enough to produce direct collapse black holes of mass $\sim 10^5$ M$_{\odot}$. The resulting clusters are larger than those produced in minihalos, but are still likely to fall short of being easily detectable in JWST blind fields.


title: The Great Escape: Discovery of a nearby 1700 km/s star ejected from the Milky Way by Sgr A*

abstract: We present the serendipitous discovery of the fastest Main Sequence hyper-velocity star (HVS) by the Southern Stellar Stream Spectroscopic Survey (S5). The star S5-HVS1 is a $\sim 2.35$ M$_\odot$ A-type star located at a distance of $\sim 9$ kpc from the Sun and has a heliocentric radial velocity of $1017\pm 2.7$ km/s without any signature of velocity variability. The current 3-D velocity of the star in the Galactic frame is $1755\pm50$ km/s. When integrated backwards in time, the orbit of the star points unambiguously to the Galactic Centre, implying that S5-HVS1 was kicked away from Sgr A* with a velocity of $\sim 1800$ km/s and travelled for $4.8$ Myr to its current location. This is so far the only HVS confidently associated with the Galactic Centre. S5-HVS1 is also the first hyper-velocity star to provide constraints on the geometry and kinematics of the Galaxy, such as the Solar motion $V_{y,\odot}= 246.1\pm 5.3$ km/s or position $R_0=8.12\pm 0.23$ kpc. The ejection trajectory and transit time of S5-HVS1 coincide with the orbital plane and age of the annular disk of young stars at the Galactic centre, and thus may be linked to its formation. With the S5-HVS1 ejection velocity being almost twice the velocity of other hyper-velocity stars previously associated with the Galactic Centre, we question whether they have been generated by the same mechanism or whether the ejection velocity distribution has been constant over time.


title: A model for core formation in dark matter haloes and ultra diffuse galaxies by outflow episodes

abstract: We present a simple model for the response of a dissipationless spherical system to an instantaneous mass change at its center, describing the formation of flat cores in dark matter haloes and ultra-diffuse galaxies (UDGs) from feedback-driven outflow episodes in a specific mass range. This model generalizes an earlier simplified analysis of an isolated shell into a system with continuous density, velocity and potential profiles. The response is divided into an instantaneous change of potential at constant velocities due to a given mass loss or gain, followed by energy-conserving relaxation to a new Jeans equilibrium. The halo profile is modeled by a two-parameter function with a variable inner slope and an analytic potential profile (Dekel et al. 2017), which enables determining the associated kinetic energy at equilibrium. The model is tested against NIHAO cosmological zoom-in simulations, where it successfully predicts the evolution of the inner dark-matter profile between successive snapshots in about 75% of the cases, failing mainly in merger situations. This model provides a simple understanding of the formation of dark-matter halo cores and UDGs by supernova-driven outflows, and a useful analytic tool for studying such processes.


title: On the origin of magnetic driven winds and the structure of the galactic dynamo in isolated galaxies

abstract: We investigate the build-up of the galactic dynamo and subsequently the origin of a magnetic driven outflow. We use a setup of an isolated disc galaxy with a realistic circum-galactic medium (CGM). We find good agreement of the galactic dynamo with theoretical and observational predictions from the radial and toroidal components of the magnetic field as function of radius and disc scale height. We find several field reversals indicating dipole structure at early times and quadrupole structure at late times. Together with the magnetic pitch angle and the dynamo control parameters $R_{\alpha}$, $R_{\omega}$ and $D$ we present strong evidence for an $\alpha^2$-$\Omega$ dynamo. The formation of a bar in the centre leads to further amplification of the magnetic field via adiabatic compression which subsequently drives an outflow. Due to the Parker-Instability the magnetic field lines rise to the edge of the disc, break out and expand freely in the CGM driven by the magnetic pressure. Finally, we investigate the correlation between magnetic field and star formation rate. Globally, we find that the magnetic field is increasing as function of the star formation rate surface density with a slope between $0.3$ and $0.45$ in good agreement with predictions from theory and observations. Locally, we find that the magnetic field can decrease while star formation increases. We find that this effect is correlated with the diffusion of magnetic field from the spiral arms to the inter-arm regions which we explicitly include by solving the induction equation and accounting for non-linear terms.


title: Ground state phase diagram of the doped Hubbard model on the 4-leg cylinder

abstract: We study the ground state properties of the Hubbard model on a 4-leg cylinder with doped hole concentration per site $\delta\leq 12.5\%$ using density-matrix renormalization group. By keeping a large number of states for long system sizes, we find that the nature of the ground state is remarkably sensitive to the presence of next-nearest-neighbor hopping $t'$. Without $t'$ the ground state of the system corresponds with the insulating filled stripe phase with long-range charge-density-wave (CDW) order and short-range incommensurate spin correlations appears. However, for a small negative $t'$ a phase characterized by coexisting algebraic d-wave superconducting (SC)- and algebraic CDW correlations. In addition, it shows short range spin- and fermion correlations consistent with a canonical Luther-Emery (LE) liquid, except that the charge- and spin periodicities are consistent with half-filled stripes instead of the $4 k_F$ and $2 k_F$ wavevectors generic for one dimensional chains. For a small positive $t'$ yet another phase takes over showing similar SC and CDW correlations. However, the fermions are now characterized by a (near) infinite correlation length while the gapped spin system is characterized by simple staggered antiferromagnetic correlations. We will show that this is consistent with a LE formed from a weakly coupled (BCS like) d-wave superconductor on the ladder where the interactions have only the effect to stabilize a cuprate style magnetic resonance.


title: Exponential suppression of bit-flips in a qubit encoded in an oscillator

abstract: A quantum system interacts with its environment, if ever so slightly, no matter how much care is put into isolating it. As a consequence, quantum bits (qubits) undergo errors, putting dauntingly difficult constraints on the hardware suitable for quantum computation. New strategies are emerging to circumvent this problem by encoding a qubit non-locally across the phase space of a physical system. Since most sources of decoherence are due to local fluctuations, the foundational promise is to exponentially suppress errors by increasing a measure of this non-locality. Prominent examples are topological qubits which delocalize quantum information over real space and where spatial extent measures non-locality. In this work, we encode a qubit in the field quadrature space of a superconducting resonator endowed with a special mechanism that dissipates photons in pairs. This process pins down two computational states to separate locations in phase space. As we increase this separation, we measure an exponential decrease of the bit-flip rate while only linearly increasing the phase-flip rate. Since bit-flips are continuously and autonomously corrected at the single qubit level, only phase-flips are left to be corrected via a one-dimensional quantum error correction code. This exponential scaling demonstrates that resonators with non-linear dissipation are promising building blocks for universal fault-tolerant quantum computation with drastically reduced hardware overhead.


title: Luminosity-duration relation of fast radio bursts

abstract: Nature of dark energy remains unknown. Especially, to constrain the time variability of the dark-energy, a new, standardisable candle that can reach more distant Universe has been awaited. Here we propose a new distance measure using fast radio bursts (FRBs), which are a new emerging population of $\sim$ ms time scale radio bursts that can reach high-$z$ in quantity. We show an empirical positive correlation between the time-integrated luminosity (L$_{\nu}$) and rest-frame intrinsic duration ($w_{\rm int,rest}$) of FRBs. The L$_{\nu}-w_{\rm int,rest}$ correlation is with a weak strength but statistically very significant, i.e., Pearson coefficient is $\sim$ 0.5 with p-value of $\sim$0.038, despite the smallness of the current sample. This correlation can be used to measure intrinsic luminosity of FRBs from the observed $w_{\rm int,rest}$. By comparing the luminosity with observed flux, we measure luminosity distances to FRBs, and thereby construct the Hubble diagram. This FRB cosmology with the L$_{\nu}-w_{\rm int,rest}$ relation has several advantages over SNe Ia, Gamma-Ray Burst (GRB), and well-known FRB dispersion measure (DM)-$z$ cosmology; (i) access to higher redshift Universe beyond the SNe Ia, (ii) high event rate that is $\sim$ 3 order of magnitude more frequent than GRBs, and (iii) it is free from the uncertainty from intergalactic electron density models, i.e., we can remove the largest uncertainty in the well-debated DM-$z$ cosmology of FRB. Our simulation suggests that the L$_{\nu}-w_{\rm int,rest}$ relation provides us with useful constraints on the time variability of the dark energy when the next generation radio telescopes start to find FRBs in quantity.


title: Quantifying the Thermal Sunyaev-Zel'dovich Effect and Excess Millimeter Emission in Quasar Environments

abstract: In this paper we probe the hot, post-shock gas component of quasar-driven winds through the thermal Sunyaev-Zel'dovich (tSZ) effect. Combining datasets from the Atacama Cosmology Telescope, the $\textit{Herschel}$ Space Observatory, and the Very Large Array, we measure average spectral energy distributions (SEDs) of 109,829 optically-selected, radio quiet quasars from 1.4~GHz to 3000~GHz in six redshift bins between $0.3<z<3.5$. We model the emission components in the radio and far-infrared, plus a spectral distortion from the tSZ effect. At $z>1.91$, we measure the tSZ effect at $3.8\sigma$ significance with an amplitude corresponding to a total thermal energy of $3.1\times10^{60}$ ergs. If this energy is due to virialized gas, then our measurement implies quasar host halo masses are $\sim6\times10^{12}~h^{-1}$M$_\odot$. Alternatively, if the host dark matter halo masses are $\sim2\times10^{12}~h^{-1}$M$_\odot$ as some measurements suggest, then we measure a $>$90 per cent excess in the thermal energy over that expected due to virialization. If the measured SZ effect is primarily due to hot bubbles from quasar-driven winds, we find that $(5^{+1.2}_{-1.3}$) per cent of the quasar bolometric luminosity couples to the intergalactic medium over a fiducial quasar lifetime of 100 Myr. An additional source of tSZ may be correlated structure, and further work is required to separate the contributions. At $z\leq1.91$, we detect emission at 95 and 148~GHz that is in excess of thermal dust and optically thin synchrotron emission. We investigate potential sources of this excess emission, finding that CO line emission and an additional optically thick synchrotron component are the most viable candidates.


title: Bouncing cosmology in $f(R,\mathcal{G})$ gravity by order reduction

abstract: A bouncing universe is a viable candidate to solve the initial singularity problem. Here we consider bouncing solutions in the context of $f(R,\mathcal{G})$ gravity by using an order reduction technique which allows one to find solutions that are perturbatively close to General Relativity. This procedure also acts as a model selection approach. Indeed, several covariant gravitational actions leading to a bounce are directly selected by demanding that the Friedmann equation derived within such gravity theories coincides with the one emerging from Loop Quantum Cosmology.


title: Sub-kiloparsec Imaging of Lyman-alpha Emission in a Low Mass, Highly Ionized, Gravitationally Lensed Galaxy at z = 1.84

abstract: Low mass, low metallicity galaxies at low to moderate ($z\lesssim3$) redshifts offer the best opportunity for detailed examination of the interplay between massive stars, ionizing radiation and gas in sources similar to those that likely reionized the universe. We present new narrowband Hubble Space Telescope observations of Ly$\alpha$ emission and the adjacent ultraviolet (UV) continuum in the low mass ($M_{\star} = 2 \times 10^8$ M$_{\odot}$), low metallicity ($Z\sim1/20$ Z$_{\odot}$) and highly ionized gravitationally lensed galaxy SL2S J02176$-$0513 at $z=1.844$. The galaxy has strong Ly$\alpha$ emission with photometric equivalent width $W^{\rm phot}_{\rm{Ly}\alpha} = 218 \pm 12$ \AA, at odds with the Ly$\alpha$ escape fraction of 10%. However, the spectroscopic Ly$\alpha$ profile suggests the presence of broad absorption underlying the emission, and the total equivalent width is consistent with the escape fraction once this underlying absorption is included. The Ly$\alpha$ emission is more spatially extended than the UV continuum, and the 0.14" spatial resolution of HST coupled with the magnification of gravitational lensing enables us to examine the distribution of Ly$\alpha$ and the UV continuum on sub-kiloparsec scales. We find that the peaks of the Ly$\alpha$ emission and the UV continuum are offset by 650 pc, and there is no Ly$\alpha$ emission arising from the region with the strongest UV light. Our combined spectroscopic and imaging data imply a significant range in neutral hydrogen column density across the object. These observations offer indirect support for a model in which ionizing radiation escapes from galaxies through channels with low column density of neutral gas.


title: The QCD topological charge and its thermal dependence: the role of the $\eta'$

abstract: We analyze the contribution of the $\eta'(958)$ meson in the first two non-trivial moments of the QCD topological charge distribution, namely, the topological susceptibility and the fourth-order cumulant of the vacuum energy density. We perform our study within U(3) Chiral Perturbation Theory up to next-to-next-to-leading order in the combined chiral and large-$N_c$ expansion. We also describe the temperature dependence of these two quantities and compare them with previous analyses in the literature. In particular, we discuss the validity of the thermal scaling of the topological susceptibility with the quark condensate, which is intimately connected with a Ward Identity relating both quantities. We also consider isospin breaking corrections from the vacuum misalignment at leading order in the U(3) framework.


title: Entanglement Dynamics after a Quench in Ising Field Theory: A Branch Point Twist Field Approach

abstract: We extend the branch point twist field approach for the calculation of entanglement entropies to time-dependent problems in 1+1-dimensional massive quantum field theories. We focus on the simplest example: a mass quench in the Ising field theory from initial mass $m_0$ to final mass $m$. The main analytical results are obtained from a perturbative expansion of the twist field one-point function in the post-quench quasi-particle basis. The expected linear growth of the R\'enyi entropies at large times $mt\gg 1$ emerges from a perturbative calculation at second order. We also show that the R\'enyi and von Neumann entropies, in infinite volume, contain subleading oscillatory contributions of frequency $2m$ and amplitude proportional to $(mt)^{-3/2}$. The oscillatory terms are correctly predicted by an alternative perturbation series, in the pre-quench quasi-particle basis, which we also discuss. A comparison to lattice numerical calculations carried out on an Ising chain in the scaling limit shows very good agreement with the quantum field theory predictions. We also find evidence of clustering of twist field correlators which implies that the entanglement entropies are proportional to the number of subsystem boundary points.

categories: hep-th cond-mat.stat-mech math-ph math.MP


title: Calibration of the charge and energy loss per unit length of the MicroBooNE liquid argon time projection chamber using muons and protons

abstract: We describe a method used to calibrate the position- and time-dependent response of the MicroBooNE liquid argon time projection chamber anode wires to ionization particle energy loss. The method makes use of crossing cosmic-ray muons to partially correct anode wire signals for multiple effects as a function of time and position, including cross-connected TPC wires, space charge effects, electron attachment to impurities, diffusion, and recombination. The overall energy scale is then determined using fully-contained beam-induced muons originating and stopping in the active region of the detector. Using this method, we obtain an absolute energy scale uncertainty of 2\% in data. We use stopping protons to further refine the relation between the measured charge and the energy loss for highly-ionizing particles. This data-driven detector calibration improves both the measurement of total deposited energy and particle identification based on energy loss per unit length as a function of residual range. As an example, the proton selection efficiency is increased by 2\% after detector calibration.


title: Analysis of Signals via Non-Maximally Decimated Non-Uniform Filter Banks

abstract: This paper addresses the important problem of reconstructing a signal from multiple multirate observations. The observations are modeled as the output of an analysis bank, and time-domain analysis is carried out to design an optimal FIR synthesis bank. We pose this as a minimizing the mean-square problem and prove that at least one optimal solution is always possible. A parametric form for all optimal solutions is obtained for a non-maximally decimated filter bank. The necessary and sufficient conditions for an optimal solution, that results in perfect reconstruction (PR), are derived as time-domain pseudocirculant conditions. This represents a novel theoretical contribution in multirate filter bank theory. We explore PR in a more general setting. This results in the ability to design a synthesis bank with a particular delay in the reconstruction. Using this delay, one can achieve PR in cases where it might not have been possible otherwise. Further, we extend the design and analysis to non-uniform filter banks and carry out simulations to verify the derived results.


title: Reconstruction of Power System Measurements Based on Enhanced Denoising Autoencoder

abstract: This paper presents a new solution for reconstructing missing data in power system measurements. An Enhanced Denoising Autoencoder (EDAE) is proposed to reconstruct the missing data through the input vector space reconstruction based on the neighbor values correlation and Long Short-Term Memory (LSTM) networks. The proposed LSTM-EDAE is able to remove the noise, extract principle features of the dataset, and reconstruct the missing information for new inputs. The paper shows that the utilization of neighbor correlation can perform better in missing data reconstruction. Trained with LSTM networks, the EDAE is more effective in coping with big data in power systems and obtains a better performance than the neural network in conventional Denoising Autoencoder. A random data sequence and the simulated Phasor Measurement Unit (PMU) data of power system are utilized to verify the effectiveness of the proposed LSTM-EDAE.


title: A Strategy for Adaptive Sampling of Multi-fidelity Gaussian Process to Reduce Predictive Uncertainty

abstract: Multi-fidelity Gaussian process is a common approach to address the extensive computationally demanding algorithms such as optimization, calibration and uncertainty quantification. Adaptive sampling for multi-fidelity Gaussian process is a changing task due to the fact that not only we seek to estimate the next sampling location of the design variable, but also the level of the simulator fidelity. This issue is often addressed by including the cost of the simulator as an another factor in the searching criterion in conjunction with the uncertainty reduction metric. In this work, we extent the traditional design of experiment framework for the multi-fidelity Gaussian process by partitioning the prediction uncertainty based on the fidelity level and the associated cost of execution. In addition, we utilize the concept of Believer which quantifies the effect of adding an exploratory design point on the Gaussian process uncertainty prediction. We demonstrated our framework using academic examples as well as a industrial application of steady-state thermodynamic operation point of a fluidized bed process


title: Environment Probing Interaction Policies

abstract: A key challenge in reinforcement learning (RL) is environment generalization: a policy trained to solve a task in one environment often fails to solve the same task in a slightly different test environment. A common approach to improve inter-environment transfer is to learn policies that are invariant to the distribution of testing environments. However, we argue that instead of being invariant, the policy should identify the specific nuances of an environment and exploit them to achieve better performance. In this work, we propose the 'Environment-Probing' Interaction (EPI) policy, a policy that probes a new environment to extract an implicit understanding of that environment's behavior. Once this environment-specific information is obtained, it is used as an additional input to a task-specific policy that can now perform environment-conditioned actions to solve a task. To learn these EPI-policies, we present a reward function based on transition predictability. Specifically, a higher reward is given if the trajectory generated by the EPI-policy can be used to better predict transitions. We experimentally show that EPI-conditioned task-specific policies significantly outperform commonly used policy generalization methods on novel testing environments.


title: Tweet Moodifier: Towards giving emotional awareness to Twitter users

abstract: Emotional contagion in online social networks has been of great interest over the past years. Previous studies have focused mainly on finding evidence of affect contagion in homophilic atmospheres. However, these studies have overlooked users' awareness of the sentiments they share and consume online. In this paper, we present an experiment with Twitter users that aims to help them better understand which emotions they experience on this social network. We introduce Tweet Moodifier (T-Moodifier), a Google Chrome extension that enables Twitter users to filter and make explicit (through colored visual marks) the emotional content in their News Feed. We compare behavioral changes between 55 participants and 5089 of their public "friends." The comparison period spans from two weeks before installing T-Moodifier to one week thereafter. The results suggest that the use of T-Moodifier might help Twitter users increase their emotional awareness: T-Moodifier users who had access to emotional statistics about their posts produced a significantly higher percentage of neutral content. This behavioral change suggests that people could behave differently while using real-time mechanisms that increase their affect reflection. Also, post-experience, those who completed both pre- and post-surveys could assert more confidently the main emotions they shared and perceived on Twitter. This shows T-Moodifier's potential to effectively make users reflect on their News Feed.


title: A simple Newton method for local nonsmooth optimization

abstract: Superlinear convergence has been an elusive goal for black-box nonsmooth optimization. Even in the convex case, the subgradient method is very slow, and while some cutting plane algorithms, including traditional bundle methods, are popular in practice, local convergence is still sluggish. Faster variants depend either on problem structure or on analyses that elide sequences of "null" steps. Motivated by a semi-structured approach to optimization and the sequential quadratic programming philosophy, we describe a new bundle Newton method that incorporates second-order objective information with the usual linear approximation oracle. One representative problem class consists of maxima of several smooth functions, individually inaccessible to the oracle. Given as additional input just the cardinality of the optimal active set, we prove local quadratic convergence. A simple implementation shows promise on more general functions, both convex and nonconvex, and suggests first-order analogues.


title: SCATTERSEARCH: Visual Querying of Scatterplot Visualizations

abstract: Scatterplots are one of the simplest and most commonly-used visualizations for understanding quantitative, multidimensional data. However, since scatterplots only depict two attributes at a time, analysts often need to manually generate and inspect large numbers of scatterplots to make sense of large datasets with many attributes. We present a visual query system for scatterplots, SCATTERSEARCH, that enables users to visually search and browse through large collections of scatterplots. Users can query for other visualizations based on a region of interest or find other scatterplots that "look similar'' to a selected one. We present two demo scenarios, provide a system overview of SCATTERSEARCH, and outline future directions.


title: Simultaneously vanishing higher derived limits

abstract: In 1988, Sibe Marde\v{s}i\'{c} and Andrei Prasolov isolated an inverse system $\mathbf{A}$ with the property that the additivity of strong homology on any class of spaces which includes the closed subsets of Euclidean space would entail that $\lim^n\mathbf{A}$ (the $n^{\text{th}}$ derived limit of $\mathbf{A}$) vanishes for every $n >0$. Since that time, the question of whether it is consistent with the $\mathsf{ZFC}$ axioms that $\lim^n \mathbf{A}=0$ for every $n >0$ has remained open. It remains possible as well that this condition in fact implies that strong homology is additive on the category of metric spaces. We show that, assuming the existence of a weakly compact cardinal, it is indeed consistent with the $\mathsf{ZFC}$ axioms that $\lim^n \mathbf{A}=0$ for all $n >0$. We show this via a finite support iteration of Hechler forcings which is of weakly compact length. More precisely, we show that in any forcing extension by this iteration a condition equivalent to $\lim^n\mathbf{A}=0$ will hold for each $n>0$. This condition is of interest in its own right; namely, it is the triviality of every coherent $n$-dimensional family of certain specified sorts of partial functions $\mathbb{N}^2\to\mathbb{Z}$ which are indexed in turn by $n$-tuples of functions $f:\mathbb{N}\to\mathbb{N}$. The triviality and coherence in question here generalize the well-studied case of $n=1$.


title: An average of generalized Dedekind sums

abstract: We study a generalization of the classical Dedekind sum that incorporates two Dirichlet characters and develop properties that generalize those of the classical Dedekind sum. By calculating the Fourier transform of this generalized Dedekind sum, we obtain an explicit formula for its second moment. Finally, we derive upper and lower bounds for the second moment with nearly identical orders of magnitude.


title: Bias of Homotopic Gradient Descent for the Hinge Loss

abstract: Gradient descent is a simple and widely used optimization method for machine learning. For homogeneous linear classifiers applied to separable data, gradient descent has been shown to converge to the maximal margin (or equivalently, the minimal norm) solution for various smooth loss functions. The previous theory does not, however, apply to non-smooth functions such as the hinge loss which is widely used in practice. Here, we study the convergence of a homotopic variant of gradient descent applied to the hinge loss and provide explicit convergence rates to the max-margin solution for linearly separable data.


title: Determination of the QCD coupling from the static energy and the free energy

abstract: We present two determinations of the strong coupling \(\alpha_s\). The first one is from the static energy at three-loop accuracy, and may be considered an update of earlier determinations by some of us. The new analysis includes new lattice data at smaller lattice spacings, and reaches distances as short as \(0.0237\,{\rm fm}\). We present a comprehensive and detailed estimate of the error sources that contribute to the uncertainty of the final result, $\alpha_s(M_Z)=0.11660^{+0.00110}_{-0.00056}$. The second determination is based on lattice data for the singlet free energy at finite temperature up to distances as small as \(0.0081\,{\rm fm}\), from which we obtain $\alpha_s(M_Z)=0.11638^{+0.0009 5}_{-0.00087}$.