title: Towards an informed quest for accretion disc winds in quasars: the intriguing case of Ton 28

abstract: We report on the detection of a blueshifted Fe K absorption feature in two consecutive XMM-Newton observations of the luminous blue quasar Ton 28, at the 4$\sigma$ cumulative significance. The rest energy of 9.2 keV implies the presence of an accretion disc wind with bulk outflow velocity of $\sim$0.28$c$, while the kinetic power is most likely a few per cent of the quasar luminosity. Remarkably, Ton 28 had been specifically selected as an optimal target to reveal an ultra-fast X-ray wind based on its total luminosity ($L_{\rm bol} > 10^{46}$ erg s$^{-1}$) and [OIII] $\lambda$5007\AA\ equivalent width (EW < 6 \AA), suggestive of high accretion rate and low inclination, respectively. Other peculiar optical/UV emission-line properties include narrow H$\beta$, strong FeII and blueshifted CIV. These are key parameters in the Eigenvector 1 formalism, and are frequently found in active galaxies with ongoing accretion disc winds, hinting at a common physical explanation. Provided that the effectiveness of our selection method is confirmed with similar sources, this result could represent the first step towards the characterization of black-hole winds through multiwavelength indicators in the absence of high-quality X-ray spectra.


title: Double-Degenerate Carbon-Oxygen and Oxygen-Neon White Dwarf Mergers: A New Mechanism for Faint and Rapid Type Ia Supernovae

abstract: Type Ia supernovae (SNe Ia) originate from the thermonuclear explosion of carbon-oxygen white dwarfs (CO WDs), giving rise to luminous optical transients. A relatively common variety of subluminous SNe Ia events, referred to as SNe Iax, are believed to arise from the failed detonation of a CO WD. In this paper, we explore failed detonation SNe Ia in the context of the double-degenerate channel of merging white dwarfs. In particular, we have carried out the first fully three-dimensional simulation of the merger of a ONe WD with a CO WD. While the hot, tidally-disrupted carbon-rich disk material originating from the CO WD secondary is readily susceptible to detonation in such a merger, the ONe WD primary core is not. This merger yields a failed detonation, resulting in the ejection of a small amount of mass, and leaving behind a kicked, super-Chandrasekhar ONe WD remnant enriched by the fallback of the products of nuclear burning. The resulting outburst is a rapidly-fading optical transient with a small amount of radioactive $^{56}$Ni powering the light curve. Consequently, the ONe-CO WD merger naturally produces a very faint and rapidly-fading transient, fainter even than the faintest Type Iax events observed to date, such as SN $2008$ha and SN 2010ae. More massive ONe primaries than considered here may produce brighter and longer-duration transients.


title: Entanglement of Hard-Core Bosons in Bipartite Lattices

abstract: The entanglement of hard-core bosons in square and honeycomb lattices with nearest-neighbor interactions is estimated by means of quantum Monte Carlo simulations and spin-wave analysis. The particular U(1)-invariant form of the concurrence is used to establish a connection with observables such as density and superfluid density. For specific regimes the concurrence is expressed as a combination of boson density and superfluid density.


title: On the gaps of the spectrum of volumes of trades

abstract: A pair $\{T_0,T_1\}$ of disjoint collections of $k$-subsets (blocks) of a set $V$ of cardinality $v$ is called a $t$-$(v,k)$ trade or simply a $t$-trade if every $t$-subset of $V$ is included in the same number of blocks of $T_0$ and $T_1$. The cardinality of $T_0$ is called the volume of the trade. Using the weight distribution of the Reed--Muller code, we prove the conjecture that for every $i$ from $2$ to $t$, there are no $t$-trades of volume greater than $2^{t+1}-2^i$ and less than $2^{t+1}-2^{i-1}$ and derive restrictions on the $t$-trade volumes that are less than $2^{t+1}+2^{t-1}$.


title: Flavor symmetries and the topology of special K\"ahler structures at rank 1

abstract: We propose a method for determining the flavor charge lattice of the continuous flavor symmetry of rank-1 4d N = 2 superconformal field theories (SCFTs) and IR free gauge theories from topological invariants of the special K\"ahler structure of the mass-deformed Coulomb branches (CBs) of the theories. The method is based on the middle homology of the total space of the elliptic fibration over the CB, and is a generalization of the F-theory string web description of flavor charge lattices. The resulting lattices, which we call 'string web lattices', contain not only information about the flavor symmetry of the SCFT but also additional information encoded in the lattice metric derived from the middle homology intersection form. This additional information clearly reflects the low energy electric and magnetic charges of BPS states on the CB, but there are other properties of the string web lattice metric which we have not been able to understand in terms of properties of the BPS spectrum. We compute the string web lattices of all rank-1 SCFTs and IR free gauge theories. We find agreement with results obtained by other methods, and find in a few cases that the string web lattice gives additional information on the flavor symmetry.


title: Aligned and Spontaneous Flavor Violation

abstract: We present a systematic spurion setup called Aligned Flavor Violation (AFV) that allows for new physics couplings to quarks that are aligned with the Standard Model (SM) Yukawas, but do not necessarily share their hierarchies nor are family universal. Additionally, we show that there is an important subset of AFV called Spontaneous Flavor Violation (SFV), which naturally arises from UV completions where the quark family number and CP groups are spontaneously broken. Flavor-changing neutral currents are strongly suppressed in SFV extensions of the SM. We study SFV from an effective field theory perspective and demonstrate that SFV new physics with significant and preferential couplings to first or second generation quarks may be close to the TeV scale.


title: Nonlinear Color-Metallicity Relations of Globular Clusters. VIII. Reproducing Color Distributions of Individual Globular Cluster Systems in the Virgo and Fornax Galaxy Clusters

abstract: The color distributions of globular clusters (GCs) in individual early-type galaxies show great diversity in their morphology. Based on the conventional "linear" relationship between the colors and metallicities of GCs, the GC metallicity distributions inferred from colors and in turn their formation histories, should be as diverse as they appear. In contrast, here we show that an alternative scenario rooted in the "nonlinear" nature of the metallicity-to-color transformation points to a simpler and more coherent picture. Our simulations of the color distributions for $\sim$80 GC systems in early-type galaxies from the ACS Virgo and Fornax Cluster Surveys suggest that the majority ($\sim$70$\%$) of early-type galaxies have old ($\sim$13 Gyr) and coeval GCs. Their variety in the color distribution morphology stems mainly from one parameter, the mean metallicity of a GC system. Furthermore, the color distributions of the remaining ($\sim$30$\%$) GC systems are also explained by the nonlinearity scenario, assuming additional young or intermediate-age GCs with a number fraction of $\sim$20$\%$ of underlying old GCs. Our results reinforce the nonlinearity explanation for the GC color bimodality and provide a new perspective on early-type galaxy formation in the cluster environment, such as the Virgo and Fornax galaxy clusters.


title: On the Entropy Cone for Large Regions at Late Times

abstract: In this note we show that the holographic entanglement entropy inequalities that hold for constant time slices are also valid for large regions and times in covariant spacetimes containing collapsing black branes, where the leading part of the entropy is computed by the simplified membrane theory of arXiv:1803.10244.


title: Anisotropic and controllable Gilbert-Bloch dissipation in spin valves

abstract: Spin valves form a key building block in a wide range of spintronic concepts and devices from magnetoresistive read heads to spin-transfer-torque oscillators. We elucidate the dependence of the magnetic damping in the free layer on the angle its equilibrium magnetization makes with that in the fixed layer. The spin pumping-mediated damping is anisotropic and tensorial, with Gilbert- and Bloch-like terms. Our investigation reveals a mechanism for tuning the free layer damping in-situ from negligible to a large value via the orientation of fixed layer magnetization, especially when the magnets are electrically insulating. Furthermore, we expect the Bloch contribution that emerges from the longitudinal spin accumulation in the non-magnetic spacer to play an important role in a wide range of other phenomena in spin valves.


title: Multi-scale Quantum Criticality driven by Kondo-lattice Coupling in Pyrochlore Systems

abstract: Pyrochlore systems ($A_{2}B_{2}O_{7}$) with $A$-site rare-earth local moments and $B$-site $5d$ conduction electrons offer excellent material platforms for the discovery of exotic quantum many-body ground states. Notable examples include U(1) quantum spin liquid of the local moments and semimetallic non-Fermi liquid of the conduction electrons. Here we investigate emergent quantum phases and their transitions driven by the Kondo lattice coupling between such highly entangled quantum ground states. Using the renormalization group method, it is shown that weak Kondo lattice coupling is irrelevant, leading to a fractionalized semimetal phase with decoupled local moments and conduction electrons. Upon increasing the Kondo lattice coupling, this phase is unstable to the formation of broken symmetry states. Particularly important is the opposing influence of the Kondo lattice coupling and long-range Coulomb interaction. The former prefers to break the particle-hole symmetry while the latter tends to restore it. The characteristic competition leads to possibly multiple phase transitions, first from a fractionalized semimetal phase to a fractionalized Fermi surface state with particle-hole pockets, followed by the second transition to a fractionalized ferromagnetic state. Multiscale quantum critical behaviors appear at nonzero temperatures and with external magnetic field near such quantum phase transitions. We discuss the implication of these results to the experiments on Pr$_2$Ir$_2$O$_7$.


title: Roll of the Dice: A Stochastically Sampled IMF Alters the Stellar Content of Simulated Dwarf Galaxies

abstract: Cosmological simulations are reaching the resolution necessary to study ultra-faint dwarf galaxies. Observations indicate that in small populations, the stellar initial mass function (IMF) is not fully populated; rather, stars are sampled in a way that can be approximated as coming from an underlying probability density function. To ensure the accuracy of cosmological simulations in the ultra-faint regime, we present an improved treatment of the IMF. We implement a self-consistent, stochastically populated IMF in cosmological hydrodynamic simulations. We test our method using high-resolution simulations of a Milky Way halo, run to $z=6$, yielding a sample of nearly 100 galaxies. We also use an isolated dwarf galaxy to investigate the resulting systematic differences in galaxy properties. We find that a stochastic IMF in simulations makes feedback burstier, strengthening feedback, and quenching star formation earlier in small dwarf galaxies. For galaxies in halos with mass $\lesssim10^{8.5}$ M$_\odot$, a stochastic IMF typically leads to lower stellar mass compared to a continuous IMF, sometimes by more than an order of magnitude. We show that existing methods of ensuring discrete supernovae incorrectly determine the mass of the star particle and its associated feedback. This leads to overcooling of surrounding gas, with at least ${\sim}$10 per cent higher star formation and ${\sim}$30 per cent higher cold gas content. Going forward, to accurately model dwarf galaxies and compare to observations, it will be necessary to incorporate a stochastically populated IMF that samples the full spectrum of stellar masses.


title: Origin of 1I/'Oumuamua. II. An ejected exo-Oort cloud object?

abstract: 1I/'Oumuamua is the first detected interstellar interloper. We test the hypothesis that it is representative of a background population of exo-Oort cloud objects ejected under the effect of post-main sequence mass loss and stellar encounters. We do this by comparing the cumulative number density of interstellar objects inferred from the detection of 1I/'Oumuamua to that expected from these two clearing processes. We consider the 0.08--8 M$_{\odot}$ mass range, take into account the dependencies with stellar mass, Galactocentric distance, and evolutionary state, and consider a wide range of size distributions for the ejected objects. Our conclusion is that 1I/'Oumuamua is likely not representative of this background population, even though there are large uncertainties in the masses and size distributions of the exo-Oort Clouds. We discuss whether the number density of free-floating, planetary-mass objects derived from gravitational microlensing surveys could be used as a discriminating measurement regarding 1I/'Oumuamua's origin (given their potential common origin). We conclude that this is challenged by the mass limitation of the surveys and the resulting uncertainty of the mass distribution of the free floaters. The detection of interlopers may be one of the few observational constraints of the small end of this population, with the caveat that, as we conclude here and in Moro-Mart\'{\i}n (2018), in the case of 1I/'Oumuamua, it might not be appropriate to assume this object is representative of an isotropic background population, which makes the derivation of a number density very challenging.


title: The Cosmological Bootstrap: Inflationary Correlators from Symmetries and Singularities

abstract: Scattering amplitudes at weak coupling are highly constrained by Lorentz invariance, locality and unitarity, and depend on model details only through coupling constants and particle content. In this paper, we develop an understanding of inflationary correlators which parallels that of flat-space scattering amplitudes. Specifically, we study slow-roll inflation with weak couplings to extra massive particles, for which all correlators are controlled by an approximate conformal symmetry on the boundary of the spacetime. After classifying all possible contact terms in de Sitter space, we derive an analytic expression for the four-point function of conformally coupled scalars mediated by the tree-level exchange of massive scalars. Conformal symmetry implies that the correlator satisfies a pair of differential equations with respect to spatial momenta, encoding bulk time evolution in purely boundary terms. The absence of unphysical singularities completely fixes this correlator. A spin-raising operator relates it to the correlators associated with the exchange of particles with spin, while weight-shifting operators map it to the four-point function of massless scalars. We explain how these de Sitter four-point functions can be perturbed to obtain inflationary three-point functions. We reproduce many classic results in the literature and provide a complete classification of all inflationary three- and four-point functions arising from weakly broken conformal symmetry. The inflationary bispectrum associated with the exchange of particles with arbitrary spin is completely characterized by the soft limit of the simplest scalar-exchange four-point function of conformally coupled scalars and a series of contact terms. Finally, we demonstrate that the inflationary correlators contain flat-space scattering amplitudes via a suitable analytic continuation of the external momenta.


title: ZFOURGE: Extreme 5007$\AA$ emission may be a common early-lifetime phase for star-forming galaxies at $z > 2.5$

abstract: Using the \prospector\ spectral energy distribution (SED) fitting code, we analyze the properties of 19 Extreme Emission Line Galaxies (EELGs) identified in the bluest composite SED in the \zfourge\ survey at $2.5 \leq z \leq 4$. \prospector\ includes a physical model for nebular emission and returns probability distributions for stellar mass, stellar metallicity, dust attenuation, and nonparametric star formation history (SFH). The EELGs show evidence for a starburst in the most recent 50 Myr, with the median EELG having a specific star formation rate (sSFR) of 4.6 Gyr$^{-1}$ and forming 15\% of its mass in this short time. For a sample of more typical star-forming galaxies (SFGs) at the same redshifts, the median SFG has a sSFR of 1.1 Gyr$^{-1}$ and forms only $4\%$ of its mass in the last 50 Myr. We find that virtually all of our EELGs have rising SFHs, while most of our SFGs do not. From our analysis, we hypothesize that many, if not most, star-forming galaxies at $z \geq 2.5$ undergo an extreme H$\beta$+$[\hbox{{\rm O}\kern 0.1em{\sc iii}}]$ emission line phase early in their lifetimes. In a companion paper, we obtain spectroscopic confirmation of the EELGs as part of our {\sc MOSEL} survey. In the future, explorations of uncertainties in modeling the UV slope for galaxies at $z>2$ are needed to better constrain their properties, e.g. stellar metallicities.


title: Determining the Nonperturbative Collins-Soper Kernel From Lattice QCD

abstract: At small transverse momentum $q_T$, transverse-momentum dependent parton distribution functions (TMDPDFs) arise as genuinely nonperturbative objects that describe Drell-Yan like processes in hadron collisions as well as semi-inclusive deep-inelastic scattering. TMDPDFs naturally depend on the hadron momentum, and the associated evolution is determined by the Collins-Soper equation. For $q_T \sim \Lambda_\mathrm{QCD}$ the corresponding evolution kernel (or anomalous dimension) is nonperturbative and must be determined as an independent ingredient in order to relate TMDPDFs at different scales. We propose a method to extract this kernel using lattice QCD and the Large-Momentum Effective Theory, where the physical TMD correlation involving light-like paths is approximated by a quasi TMDPDF, defined using equal-time correlation functions with a large-momentum hadron state. The kernel is determined from a ratio of quasi TMDPDFs extracted at different hadron momenta.


title: Fraunhofer response and supercurrent spin switching in black phosphorus with strain and disorder

abstract: We develop theory models for both ballistic and disordered superconducting monolayer black phosphorus devices in the presence of magnetic exchange field and strain. The ballistic case is studied through a microscopic Bogoliubov-de Gennes formalism while for the disordered case we formulate a quasiclassical model. Utilizing the two models, we theoretically study the response of supercurrent to an externally applied magnetic field in two-dimensional black phosphorus Josephson junctions. Our results demonstrate that the response of the supercurrent to a perpendicular magnetic field in ballistic samples can deviate from the standard Fraunhofer interference pattern when the Fermi level and mechanical strain are varied. This finding suggests the combination of chemical potential and strain is an efficient external knob to control the current response in highly sensitive strain-effect transistors and superconducting quantum interference devices. We also study the supercurrent in a superconductor-ferromagnet-ferromagnet-superconductor junction where the magnetizations of the two adjacent magnetized regions are uniform with misaligned orientations. We show that the magnetization misalignment can control the excitation of harmonics higher than the first harmonic sin(phi) (in which phi is the phase difference between the superconductors) in supercurrent and constitutes a full spin switching current element. Finally, we discuss possible experimental implementations of our findings. We foresee our models and discussions could provide guidelines to experimentalists in designing devices and future investigations.


title: Ballistic Dark Matter oscillates above $\Lambda$CDM

abstract: Dark matter may have been relativistic and collisional until relatively late times and become cold and collisionless after a phase transition before the matter-radiation equality of the standard $\Lambda$CDM cosmology. We show that such a dark matter has large peculiar velocities due to acoustic oscillations before the phase transition, and evolves ballistically after the phase transition in the collisionless phase until the initial acoustic velocities are redshifted away. We show that this Ballistic Dark Matter (BDM) results in new non-trivial interesting features in the cosmological observables. In particular, the linear matter power spectrum exhibits acoustic oscillations on scales smaller than the Hubble scale at the time of phase transition, and for fast transitions the power at the acoustic peaks in the matter power spectrum exceeds that in a $\Lambda$CDM cosmology. If BDM only forms a part of the total dark matter, an odd vs. even acoustic peak asymmetry becomes prominent. We give an approximate analytical treatment of the linear perturbations in BDM, explaining these features. We also discuss the possibility to constrain BDM using cosmic microwave background and large scale structure data.


title: Evolution of entanglement spectra under generic quantum dynamics

abstract: We characterize the early stages of the approach to equilibrium in isolated quantum systems through the evolution of the entanglement spectrum. We find that the entanglement spectrum of a subsystem evolves with at least three distinct timescales. First, on an o(1) timescale, independent of system or subsystem size and the details of the dynamics, the entanglement spectrum develops nearest-neighbor level repulsion. The second timescale sets in when the light-cone has traversed the subsystem. Between these two times, the density of states of the reduced density matrix takes a universal, scale-free 1/f form; thus, random-matrix theory captures the local statistics of the entanglement spectrum but not its global structure. The third time scale is that on which the entanglement saturates; this occurs well after the light-cone traverses the subsystem. Between the second and third times, the entanglement spectrum compresses to its thermal Marchenko-Pastur form. These features hold for chaotic Hamiltonian and Floquet dynamics as well as a range of quantum circuit models.


title: Detections of [OIII] 88 $\mu$m in Two Quasars in the Reionization Epoch

abstract: With the Atacama Large Millimeter/sub-millimeter Array (ALMA), we report detections of the far-infrared (FIR) [OIII] 88 $\mu$m line and the underlying dust continuum in the two quasars in the reionization epoch, J205406.48-000514.8 (hereafter J2054-0005) at $z=6.0391\pm0.0002$ and J231038.88+185519.7 (hereafter J2310+1855) at $z=6.0035\pm0.0007$. The [OIII] luminosity of J2054-0005 and J2310+1855 are $L_{\rm [OIII]} = (6.8\pm0.6) \times 10^{9}$ and $(2.4\pm0.6) \times 10^{9}$ $L_{\rm \odot}$, corresponding to $\approx$ 0.05\%\ and 0.01\%\ of the total infrared luminosity, $L_{\rm TIR}$, respectively. Combining these [OIII] luminosities with [CII] 158 $\mu$m luminosities in the literature, we find that J2054-0005 and J2310+1855 have the [OIII]-to-[CII] luminosity ratio of $2.1\pm0.4$ and $0.3\pm0.1$, respectively, the latter of which is the lowest among objects so far reported at $z>6$. Combining [OIII] observations at $z\approx6-9$ from our study and the literature, we identify the [OIII] line deficit: objects with larger $L_{\rm TIR}$ have lower $L_{\rm [OIII]}$-to-$L_{\rm TIR}$ ratios. Furthermore, we also find that the anti-correlation is shifted toward higher $L_{\rm TIR}$ value when compared to the local [OIII] line deficit.


title: Neutrinoless Double-$\beta$ Decay with Nonstandard Majoron Emission

abstract: We present a novel mode of neutrinoless double beta decay with emission of a light Majoron-like scalar particle phi. We assume it couples via an effective seven-dimensional operator with a (V+A) lepton current and (V+-A) quark currents leading to a long-range contribution that is unsuppressed by the light neutrino mass. We calculate the total double beta decay rate and determine the fully differential shape for this mode. We find that future double beta decay searches are sensitive to scales of the order Lambda_NP ~ 1 TeV for the effective operator and a light scalar m_phi < 0.2 MeV, based on ordinary double beta decay Majoron searches. The angular and energy distributions can deviate considerably from that of two-neutrino double beta decay, which is the main background. We point out possible ultraviolet completions where such an effective operator can emerge.


title: Strings in Ramond-Ramond Backgrounds from the Neveu-Schwarz-Ramond Formalism

abstract: We treat RR flux backgrounds of type II string theory in the framework of closed superstring field theory based on the NSR formalism, focusing on two examples: (1) the pp-wave background supported by 5-form flux, and (2) $AdS_3\times S^3\times M_4$ supported by mixed 3-form fluxes. In both cases, we analyze the classical string field solution perturbatively, and compute the correction to the dispersion relation of string states to quadratic order in the RR flux. In the first example, our result is in a delicate way consistent with that obtained from lightcone quantization of the Green-Schwarz string. In the second example, we will obtain numerically the mass corrections to pulsating type IIB strings in $AdS_3\times S^3\times M_4$. Our results, valid at finite $AdS$ radius, agree with previously known answers in the semiclassical limit and in the BMN limit respectively.


title: Dark Matter Bars in Spinning Halos

abstract: We study nonlinear response of spinning dark matter (DM) halos to dynamic and secular evolution of stellar bars in the embedded galactic disks, using high-resolution numerical simulations. For a sequence of halos with the cosmological spin parameter lambda=0-0.09, and a representative angular momentum distribution, we analyze evolution of induced DM bars amplitude and quantify parameters of the response as well as trapping of DM orbits and angular momentum transfer by the main and secondary resonances. We find that (1) maximal amplitude of DM bars depends strongly on `lambda', while that of the stellar bars is indifferent to `lambda'; (2) Efficiency of resonance trapping of DM orbits by the bar increases with `lambda', and so is the mass and the volume of DM bars; (3) Contribution of resonance transfer of angular momentum to the DM halo increases with `lambda', and for larger spin, the DM halo `talks' to itself, by moving the angular momentum to larger radii --- this process is maintained by resonances; (4) Prograde and retrograde DM orbits play different roles in angular momentum transfer. The `active' part of the halo extends well beyond the bar region, up to few times the bar length in equatorial plane and away from this plane. (5) We model evolution of diskless DM halos and halos with frozen disks, and found them to be perfectly stable to any Fourier modes. Finally, further studies adopting a range of mass and specific angular momentum distributions of the DM halo will generalize the dependence of DM response on the halo spin and important implications for direct detection of DM and that of the associated stellar tracers, such as streamers.


title: A Catalogue of Galactic Supernova Remnants in the far-infrared: revealing ejecta dust in pulsar wind nebulae

abstract: We search for far-infrared (FIR) counterparts of known supernova remnants (SNRs) in the Galactic plane (10 deg <| l |< 60 deg) at 70-500 micron using the Herschel Infrared Galactic Plane Survey (Hi-GAL). Of 71 sources studied, we find that 29 (41 per cent) SNRs have a clear FIR detection of dust emission associated with the SNR. Dust from 8 of these is in the central region, and 4 indicate pulsar wind nebulae (PWNe) heated ejecta dust. A further 23 have dust emission in the outer shell structures which is potentially related to swept up material. Many Galactic SNe have dust signatures but we are biased towards detecting ejecta dust in young remnants and those with a heating source (shock or PWN). We estimate the dust temperature and mass contained within three PWNe, G11.2-0.3, G21.5-0.9, and G29.7-0.3 using modified blackbody fits. To more rigorously analyse the dust properties at various temperatures and dust emissivity index beta, we use point process mapping (PPMAP). We find significant quantities of cool dust (at 20-40 K) with dust masses of Md = 0.34 +/- 0.14 solar mass, Md = 0.29 +/- 0.08 solar mass, and Md = 0.51 +/- 0.13 solar mass for G11.2-0.3, G21.5-0.9, and G29.7-0.3 respectively. We derive the dust emissivity index for the PWN ejecta dust in G21.5-0.3 to be beta = 1.4 +/- 0.5 compared to dust in the surrounding medium where beta = 1.8 +/- 0.1.


title: TOPoS V: Abundance ratios in a sample of very metal-poor turn-off stars

abstract: Extremely metal-poor stars are keys to understand the early evolution of our Galaxy. The ESO large programme TOPoS has been tailored to analyse a new set of metal-poor turn-off stars, whereas most of the previously known extremely metal-poor stars are giant stars. Sixty five turn-off stars (preselected from SDSS spectra) have been observed with the X-Shooter spectrograph at the ESO VLT Unit Telescope 2, to derive accurate and detailed abundances of magnesium, silicon, calcium, iron, strontium and barium. We analysed medium-resolution spectra (R ~ 10 000) obtained with the ESO X-Shooter spectrograph and computed the abundances of several alpha and neutron-capture elements using standard one-dimensional local thermodynamic equilibrium (1D LTE) model atmospheres. Our results confirms the super-solar [Mg/Fe] and [Ca/Fe] ratios in metal-poor turn-off stars as observed in metal-poor giant stars. We found a significant spread of the [alpha/Fe] ratios with several stars showing sub-solar [Ca/Fe] ratios. We could measure the abundance of strontium in 12 stars of the sample, leading to abundance ratios [Sr/Fe] around the Solar value. We detected barium in two stars of the sample. One of the stars (SDSS J114424-004658) shows both very high [Ba/Fe] and [Sr/Fe] abundance ratios (>1 dex).


title: Demonstration of Einstein-Podolsky-Rosen Steering Using Hybrid Continuous- and Discrete-Variable Entanglement of Light

abstract: Einstein-Podolsky-Rosen steering is known to be a key resource for one-sided device-independent quantum information protocols. Here we demonstrate steering using hybrid entanglement between continuous- and discrete-variable optical qubits. To this end, we report on suitable steering inequalities and detail the implementation and requirements for this demonstration. Steering is experimentally certified by observing a violation by more than 5 standard deviations. Our results illustrate the potential of optical hybrid entanglement for applications in heterogeneous quantum networks that would interconnect disparate physical platforms and encodings.


title: Modulation of the Work Function by the Atomic Structure of Strong Organic Electron Acceptors on H-Si(111)

abstract: Advances in hybrid organic/inorganic architectures for optoelectronics can be achieved by understanding how the atomic and electronic degrees of freedom cooperate or compete to yield the desired functional properties. Here we show how work-function changes are modulated by the structure of the organic components in model hybrid systems. We consider two cyano-quinodimethane derivatives (F4-TCNQ and F6-TCNNQ), which are strong electron-acceptor molecules, adsorbed on H-Si(111). From systematic structure searches employing range-separated hybrid HSE06 functional including many body van der Waals contributions, we predict that despite their similar composition, these molecules adsorb with significantly different densely-packed geometries in the first layer, due to strong intermolecular interaction. F6-TCNNQ shows a much stronger intralayer interaction (primarily due to van der Waals contributions) than F4-TCNQ in multilayered structures. The densely-packed geometries induce a large interface-charge rearrangement that result in a work-function increase of 1.11 and 1.76 eV for F4-TCNQ and F6-TCNNQ, respectively. Nuclear fluctuations at room temperature produce a wide distribution of work-function values, well modeled by a normal distribution with {\sigma}=0.17 eV. We corroborate our findings with experimental evidence of pronounced island formation for F6-TCNNQ on H-Si(111) and with the agreement of trends between predicted and measured work-function changes.


title: E$_{7(7)}$ Exceptional Field Theory in Superspace

abstract: We formulate the locally supersymmetric E$_{7(7)}$ exceptional field theory in a $(4+56|32)$ dimensional superspace, corresponding to a 4D $N\!=\!8$ "external" superspace augmented with an "internal" 56-dimensional space. This entails the unification of external diffeomorphisms and local supersymmetry transformations into superdiffeomorphisms. The solutions to the superspace Bianchi identities lead to on-shell duality equations for the $p$-form field strengths for $p\leq 4$. The reduction to component fields provides a complete description of the on-shell supersymmetric theory. As an application of our results, we perform a generalized Scherk-Schwarz reduction and obtain the superspace formulation of maximal gauged supergravity in four dimensions parametrized by an embedding tensor.


title: Sign-changing blowing-up solutions for the critical nonlinear heat equation

abstract: Let $\Omega$ be a smooth bounded domain in $\mathbb{R}^n$ and denote the regular part of the Green's function on $\Omega$ with Dirichlet boundary condition as $H(x,y)$. Assume that $q \in \Omega$ and $n\geq 5$. We prove that there exists an integer $k_0$ such that for any integer $k\geq k_0$ there exist initial data $u_0$ and smooth parameter functions $\xi(t)\to q$, $0<\mu(t)\to 0$ as $t\to +\infty$ such that the solution $u_q$ of the critical nonlinear heat equation \begin{equation*} \begin{cases} u_t = \Delta u + |u|^{\frac{4}{n-2}}u\text{ in } \Omega\times (0, \infty),\\ u = 0\text{ on } \partial \Omega\times (0, \infty),\\ u(\cdot, 0) = u_0 \text{ in }\Omega, \end{cases} \end{equation*} has the form \begin{equation*} u_q(x, t) \approx \mu(t)^{-\frac{n-2}{2}}\left(Q_k\left(\frac{x-\xi(t)}{\mu(t)}\right) - H(x, q)\right), \end{equation*} where the profile $Q_k$ is the non-radial sign-changing solution of the Yamabe equation \begin{equation*} \Delta Q + |Q|^{\frac{4}{n-2}}Q = 0\text{ in }\mathbb{R}^n, \end{equation*} constructed in \cite{delpinomussofrankpistoiajde2011}. In dimension 5 and 6, we also prove the stability of $u_q(x, t)$.


title: Interaction of surface acoustic waves with electronic excitations in graphene

abstract: This article reviews the main theoretical and experimental advances regarding the interaction between surface acoustic waves (SAWs) and electronic excitations in graphene. The coupling of the graphene electron gas to the SAW piezoelectric field can modify the propagation properties of the SAW, and even amplify the intensity of SAWs traveling along the graphene layer. Conversely, the periodic electric and strain fields of the SAW can be used to modify the graphene Dirac cone and to couple light into graphene plasmons. Finally, SAWs can generate acousto-electric currents in graphene. These increase linearly with the SAW frequency and power but, in contrast to conventional currents, they depend non-monotonously on the graphene electric conductivity. Most of these functionalities have been reported in graphene transferred to the surface of strong piezoelectric insulators. The recent observation of acousto-electric currents in epitaxial graphene on SiC opens the way to the large-scale fabrication of graphene-based acousto-electric devices patterned directly on a semi-insulating wafer.


title: The Stellar Contents of Intermediate Mass Disk Galaxies in the Virgo Cluster. I. GMOS Spectra

abstract: The stellar contents of six intermediate mass early-type disk galaxies in the Virgo cluster are examined using deep long slit spectra. The isophotal and photometric properties of the galaxies at near- and mid-infrared wavelengths are also examined. Characteristic ages and metallicities are found by making comparisons with the strengths of lines measured from model spectra. The light from the central regions of these galaxies at visible/red wavelengths is not dominated by old populations. Rather, the central regions of four galaxies (NGC 4305, NGC 4306, NGC 4497, and NGC 4620) are dominated by populations with ages ~ 1.5 - 3 Gyr. Centrally-concentrated line emission is found in two of the galaxies (NGC 4491 and NGC 4584), and the relative strengths of Halpha and [SII]6746 are consistent with this emission originating in star-forming regions. The disks of these galaxies are dominated by populations that are at least 1 Gyr older than those near the centers, indicating that the quenching of star formation progressed from large radii inwards, and did not occur over a short timescale.


title: Duality and normalization, variations on a theme of Serre and Reid

abstract: We discuss the naive duality theory of coherent, torsion free, $S_2$ sheaves on schemes. Version 2: A solution of the original Conjecture 6.2 is added as an appendix by Hailong Dao.


title: Cosmological Hydrodynamic Simulations with Suppressed Variance in the Lyman-$\alpha$ Forest Power Spectrum

abstract: We test a method to reduce unwanted sample variance when predicting Lyman-$\alpha$ (ly$\alpha$) forest power spectra from cosmological hydrodynamical simulations. Sample variance arises due to sparse sampling of modes on large scales and propagates to small scales through non-linear gravitational evolution. To tackle this, we generate initial conditions in which the density perturbation amplitudes are {\it fixed} to the ensemble average power spectrum -- and are generated in {\it pairs} with exactly opposite phases. We run $50$ such simulations ($25$ pairs) and compare their performance against $50$ standard simulations by measuring the ly$\alpha$ 1D and 3D power spectra at redshifts $z=2$, 3, and 4. Both ensembles use periodic boxes of $40$ Mpc/h containing $512^3$ particles each of dark matter and gas. As a typical example of improvement, for wavenumbers $k=0.25$ h/Mpc at $z=3$, we find estimates of the 1D and 3D power spectra converge $34$ and $12$ times faster in a paired-fixed ensemble compared with a standard ensemble. We conclude that, by reducing the computational time required to achieve fixed accuracy on predicted power spectra, the method frees up resources for exploration of varying thermal and cosmological parameters -- ultimately allowing the improved precision and accuracy of statistical inference.


title: Casimir-Polder forces in inhomogeneous backgrounds

abstract: Casimir-Polder interactions are considered in an inhomogeneous, dispersive background. We consider both the interaction between a polarizable atom and a perfectly conducting wall, and between such an atom and a plane interface between two different inhomogeneous dielectric media. Renormalization is achieved by subtracting the interaction with the local inhomogeneous medium by itself. The results are expressed in general form, and generalize the Dzyaloshinskii-Lifshitz-Pitaevskii interaction between a spatially homogeneous dielectric interface and a polarizable atom.


title: Quantum-like modeling of the order effect in decision making: POVM viewpoint on the Wang-Busemeyer QQ-equality

abstract: In recent years, quantum mechanics has been actively used in areas outside of physics, such as psychology, sociology, theory of decision-making, game theory, and others. In particular, quantum mechanics is used to explain the paradoxes arising in cognitive psychology and decision making. Wang and Busemeyer invented a quantum model and approach as well as non-parametric equality (so-called QQ-equality), explaining the questions order effect. The primary objective of this note is to test the possibility to expand the Wang-Busemeyer model by considering questions which are mathematically represented by positive operator valued measures. We found that, for such observables, the QQ-equality can be violated. But, we also showed that, in principle, it is possible to reduce expanded model to the original Wang-Busemeyer model by expanding the context of the questions. This version of preprint is aimed to point out to annoying miscalculation in version 1. This miscalculation might mislead a reader who is not experienced in operating with POVMs. Otherwise the main line of construction and reasoning presented in version 1 is right and it can be easily completed by the reader on the basis of version 1 and the correction remark in version 2.


title: Modified gravity revealed along geodesic tracks

abstract: The study of the dynamics of a two-body system in modified gravity constitutes a more complex problem than in Newtonian gravity. Numerical methods are typically needed to solve the equations of geodesics. Despite the complexity of the problem, the study of a two-body system in $f(R)$ gravity leads to a new exciting perspective hinting the right strategy to adopt in order to probe modified gravity. Our results point out some differences between the {\em semiclassical} (Newtonian) approach, and the {\em relativistic} (geodesic) one thus suggesting that the latter represents the best strategy for future tests of modified theories of gravity. { Finally, we have also highlighted the capability of forthcoming observations to serve as smoking gun of modified gravity revealing a departure from GR or further reducing the parameter space of $f(R)$ gravity}. \keywords{$f(R)$ gravity \and binary system \and geodesics \and precession


title: Intrinsic curvature and topology of shadows in Kerr spacetime

abstract: From the viewpoint of differential geometry and topology, we investigate the characterization of the shadows in a Kerr spacetime. Two new quantities, the length of the shadow boundary and the local curvature radius are introduced. Each shadow can be uniquely determined by these two quantities. For the black hole case, the result shows that we can constrain the black hole spin and the angular coordinate of the observer only by measuring the maximum and minimum of the curvature radius. While for the naked singularity case, we adopt the length parameter and the maximum of the curvature radius. This technique is completely independent of the coordinate system and the location of the shadow, and is expected to uniquely determine the parameters of the spacetime. Moreover, we propose a topological covariant quantity to measure and distinguish different topological structures of the shadows.


title: Density-matrix embedding theory study of the one-dimensional Hubbard-Holstein model

abstract: We present a density-matrix embedding theory (DMET) study of the one-dimensional Hubbard-Holstein model, which is paradigmatic for the interplay of electron-electron and electron-phonon interactions. Analyzing the single-particle excitation gap, we find a direct Peierls insulator to Mott insulator phase transition in the adiabatic regime of slow phonons in contrast to a rather large intervening metallic phase in the anti-adiabatic regime of fast phonons. We benchmark the DMET results for both on-site energies and excitation gaps against density-matrix renormalization group (DMRG) results and find excellent agreement of the resulting phase boundaries. We also compare the fully quantum treatment of phonons against the standard Born-Oppenheimer (BO) approximation. The BO approximation gives qualitatively similar results to DMET in the adiabatic regime, but fails entirely in the anti-adiabatic regime, where BO predicts a sharp direct transition from Mott to Peierls insulator, whereas DMET correctly shows a large intervening metallic phase. This highlights the importance of quantum fluctuations in the phononic degrees of freedom for metallicity in the one-dimensional Hubbard-Holstein model.


title: On subfields of the second generalization of the GK maximal function field

abstract: The second generalized GK function fields $K_n$ are a recently found family of maximal function fields over the finite field with $q^{2n}$ elements, where $q$ is a prime power and $n \ge 1$ an odd integer. In this paper we construct many new maximal function fields by determining various Galois subfields of $K_n$. In case $\gcd(q+1,n)=1$ and either $q$ is even or $q \equiv 1 \pmod{4},$ we find a complete list of Galois subfields of $K_n.$ Our construction adds several previously unknown genera to the genus spectrum of maximal curves.


title: Linking gravitational waves and X-ray phenomena with joint LISA and Athena observations

abstract: The evolution of cosmic structures, the formation and growth of the first black holes and the connection to their baryonic environment are key unsolved problems in astrophysics. The X-ray Athena mission and the gravitational-wave Laser Interferometer Space Antenna (LISA) offer independent and complementary angles on these problems. We show that up to 10 black hole binaries in the mass range 10^5 - 10^8 Msun discovered by LISA at redshift <~ 3.5 could be detected by Athena in an exposure time up to 100 ks, if prompt X-ray emission of ~ 1% - 10% of the Eddington luminosity is present. Likewise, if any LISA-detected extreme mass ratio inspirals occur in accretion disks, Athena can detect associated electromagnetic emission out to redshift ~ 1. Finally, warned by LISA, Athena can point in advance and stare at stellar-mass binary black hole mergers at redshift <~ 0.1. These science opportunities emphasise the vast discovery space of simultaneous observations from the two observatories, which would be missed if they were operated in different epochs.


title: Generating Texts with Integer Linear Programming

abstract: Concept-to-text generation typically employs a pipeline architecture, which often leads to suboptimal texts. Content selection, for example, may greedily select the most important facts, which may require, however, too many words to express, and this may be undesirable when space is limited or expensive. Selecting other facts, possibly only slightly less important, may allow the lexicalization stage to use much fewer words, or to report more facts in the same space. Decisions made during content selection and lexicalization may also lead to more or fewer sentence aggregation opportunities, affecting the length and readability of the resulting texts. Building upon on a publicly available state of the art natural language generator for Semantic Web ontologies, this article presents an Integer Linear Programming model that, unlike pipeline architectures, jointly considers choices available in content selection, lexicalization, and sentence aggregation to avoid greedy local decisions and produce more compact texts, i.e., texts that report more facts per word. Compact texts are desirable, for example, when generating advertisements to be included in Web search results, or when summarizing structured information in limited space. An extended version of the proposed model also considers a limited form of referring expression generation and avoids redundant sentences. An approximation of the two models can be used when longer texts need to be generated. Experiments with three ontologies confirm that the proposed models lead to more compact texts, compared to pipeline systems, with no deterioration or with improvements in the perceived quality of the generated texts.


title: Some New Layer Architectures for Graph CNN

abstract: While convolutional neural networks (CNNs) have recently made great strides in supervised classification of data structured on a grid (e.g. images composed of pixel grids), in several interesting datasets, the relations between features can be better represented as a general graph instead of a regular grid. Although recent algorithms that adapt CNNs to graphs have shown promising results, they mostly neglect learning explicit operations for edge features while focusing on vertex features alone. We propose new formulations for convolutional, pooling, and fully connected layers for neural networks that make more comprehensive use of the information available in multi-dimensional graphs. Using these layers led to an improvement in classification accuracy over the state-of-the-art methods on benchmark graph datasets.


title: DEEPGONET: Multi-label Prediction of GO Annotation for Protein from Sequence Using Cascaded Convolutional and Recurrent Network

abstract: The present gap between the amount of available protein sequence due to the development of next generation sequencing technology (NGS) and slow and expensive experimental extraction of useful information like annotation of protein sequence in different functional aspects, is ever widening, which can be reduced by employing automatic function prediction (AFP) approaches. Gene Ontology (GO), comprising of more than 40, 000 classes, defines three aspects of protein function names Biological Process (BP), Cellular Component (CC), Molecular Function (MF). Multiple functions of a single protein, has made automatic function prediction a large-scale, multi-class, multi-label task. In this paper, we present DEEPGONET, a novel cascaded convolutional and recurrent neural network, to predict the top-level hierarchy of GO ontology. The network takes the primary sequence of protein as input which makes it more useful than other prevailing state-of-the-art deep learning based methods with multi-modal input, making them less applicable for proteins where only primary sequence is available. All predictions of different protein functions of our network are performed by the same architecture, a proof of better generalization as demonstrated by promising performance on a variety of organisms while trained on Homo sapiens only, which is made possible by efficient exploration of vast output space by leveraging hierarchical relationship among GO classes. The promising performance of our model makes it a potential avenue for directing experimental protein functions exploration efficiently by vastly eliminating possible routes which is done by the exploring only the suggested routes from our model. Our proposed model is also very simple and efficient in terms of computational time and space compared to other architectures in literature.


title: AdS$_6$ T-duals and AdS$_6 \times S^2$ geometries in Type IIB

abstract: We summarise the recent results in 1810.08093, where the Type IIB AdS$_6$ solutions T-dual to the Brandhuber-Oz solution to massive Type IIA supergravity are shown to fit within an extension to allow for delocalized 5- and 7-branes of the global AdS$_6 \times S^2$ solutions constructed by D'Hoker, Gutperle and Uhlemann. We show that the Abelian T-dual solution provides an explicit example of a Riemann surface with the topology of an annulus.


title: Simulations of gaussian systems in Minkowski time

abstract: Many research programs aiming to deal with the sign problem were proposed since the advent of lattice field theory. Several of these try to achieve this by exploiting properties of analytic functions. This is also the case for our study. There auxiliary complex variables are introduced and desired weight is obtained after integrating them out. In this note we clarify conceptual difficulties with this procedure encountered in previous works. In the process we observe an exciting connection with thimbles and discover an interesting hidden symmetry present in the problem. Problem of negative eigenvalues of the action will be revisited and considered from a different perspective. As a byproduct we perform simulations of simple quantum systems directly in Minkowski time.


title: A Mixture of Expert Approach for Low-Cost Customization of Deep Neural Networks

abstract: The ability to customize a trained Deep Neural Network (DNN) locally using user-specific data may greatly enhance user experiences, reduce development costs, and protect user's privacy. In this work, we propose to incorporate a novel Mixture of Experts (MOE) approach to accomplish this goal. This architecture comprises of a Global Expert (GE), a Local Expert (LE) and a Gating Network (GN). The GE is a trained DNN developed on a large training dataset representative of many potential users. After deployment on an embedded edge device, GE will be subject to customized, user-specific data (e.g., accent in speech) and its performance may suffer. This problem may be alleviated by training a local DNN (the local expert, LE) on a small size customized training data to correct the errors made by GE. A gating network then will be trained to determine whether an incoming data should be handled by GE or LE. Since the customized dataset is in general very small, the cost of training LE and GN would be much lower than that of re-training of GE. The training of LE and GN thus can be performed at local device, properly protecting the privacy of customized training data. In this work, we developed a prototype MOE architecture for handwritten alphanumeric character recognition task. We use EMNIST as the generic dataset, LeNet5 as GE, and handwritings of 10 users as the customized dataset. We show that with the LE and GN, the classification accuracy is significantly enhanced over the customized dataset with almost no degradation of accuracy over the generic dataset. In terms of energy and network size, the overhead of LE and GN is around 2.5% compared to those of GE.


title: Multidimensional staggered grid residual distribution scheme for Lagrangian hydrodynamics

abstract: We present the second-order multidimensional Staggered Grid Hydrodynamics Residual Distribution (SGH RD) scheme for Lagrangian hydrodynamics. The SGH RD scheme is based on the staggered finite element discretizations as in [Dobrev et al., SISC, 2012]. However, the advantage of the residual formulation over classical FEM approaches consists in the natural mass matrix diagonalization which allows one to avoid the solution of the linear system with the global sparse mass matrix while retaining the desired order of accuracy. This is achieved by using Bernstein polynomials as finite element shape functions and coupling the space discretization with the deferred correction type timestepping method. Moreover, it can be shown that for the Lagrangian formulation written in non-conservative form, our residual distribution scheme ensures the exact conservation of the mass, momentum and total energy. In this paper we also discuss construction of numerical viscosity approximations for the SGH RD scheme allowing to reduce the dissipation of the numerical solution. Thanks to the generic formulation of the staggered grid residual distribution scheme, it can be directly applied to both single- and multimaterial and multiphase models. Finally, we demonstrate computational results obtained with the proposed residual distribution scheme for several challenging test problems.


title: The origin of the two populations of blue stragglers in M30

abstract: We analyze the position of the two populations of blue stragglers in the globular cluster M30 in the Hertzsprung-Russell diagram. Both populations of blue stragglers are brighter than the cluster's turn-off, but one population (the blue blue-stragglers) align along the zero-age main-sequence whereas the (red) population is elevated in brightness (or colour) by $\sim 0.75$ mag. Based on stellar evolution and merger simulations we argue that the red population, which composes about 40\% of the blue stragglers in M 30, is formed at a constant rate of $\sim 2.8$ blue stragglers per Gyr over the last $\sim 10$ Gyr. The blue population is formed in a burst that started $\sim 3.2$ Gyr ago at a peak rate of $30$ blue stragglers per Gyr$^{-1}$ with an e-folding time scale of $0.93$ Gyr. We speculate that the burst resulted from the core collapse of the cluster at an age of about 9.8 Gyr, whereas the constantly formed population is the result of mass transfer and mergers through binary evolution. In that case about half the binaries in the cluster effectively result in a blue straggler.


title: Visualizing Intramolecular Distortions as the Origin of Transverse Magnetic Anisotropy

abstract: The magnetic properties of metal-organic complexes are strongly influenced by conformational changes in the ligand. The flexibility of Fe-tetra-pyridyl-porphyrin molecules leads to different adsorption configurations on a Au(111) surface. By combining low-temperature scanning tunneling spectroscopy and atomic force microscopy, we resolve a correlation of the molecular configuration with different spin states and magnitudes of magnetic anisotropy. When the macrocycle exhibits a laterally-undistorted saddle shape, the molecules lie in a S=1 state with axial anisotropy arising from a square-planar ligand field. If the symmetry in the molecular ligand field is reduced by a lateral distortion of the molecule, we find a finite contribution of transverse anisotropy. Some of the distorted molecules lie in a S=2 state, again exhibiting substantial transverse anisotropy.


title: On the Complexity of Properties of Transformation Semigroups

abstract: We investigate the computational complexity for determining various properties of a finite transformation semigroup given by generators. We introduce a simple framework to describe transformation semigroup properties that are decidable in $\mathsf{AC^0}$. This framework is then used to show that the problems of deciding whether a transformation semigroup is a group, commutative or a semilattice are in $\mathsf{AC^0}$. Deciding whether a semigroup has a left (resp.right) zero is shown to be $\mathsf{NL}$-complete, as are the problems of testing whether a transformation semigroup is nilpotent, $\mathcal{R}$-trivial or has central idempotents. We also give $\mathsf{NL}$ algorithms for testing whether a transformation semigroup is idempotent, orthodox, completely regular, Clifford or has commuting idempotents. Some of these algorithms are direct consequences of the more general result that arbitrary fixed semigroup equations can be tested in~$\mathsf{NL}$. Moreover, we show how to compute left and right identities of a transformation semigroup in polynomial time. Finally, we show that checking whether an element is regular is $\mathsf{PSPACE}$-complete. \


title: Complex Crystals from Size-disperse Spheres

abstract: Colloids are rarely perfectly uniform but follow a distribution of sizes, shapes, and charges. This dispersity can be inherent (static) or develop and change over time (dynamic). Despite a long history of research, the conditions under which non-uniform particles crystallize and which crystal forms is still not well understood. Here, we demonstrate that hard spheres with Gaussian radius distribution and dispersity up to 19% always crystallize if compressed slow enough, and they do so in surprisingly complex ways. This result is obtained by accelerating event-driven simulations with particle swap moves for static dispersity and particle resize moves for dynamic dispersity. Above 6% dispersity, AB$_2$ Laves, AB$_{13}$, and a region of Frank-Kasper phases are found. The Frank-Kasper region includes a quasicrystal approximant with Pearson symbol oS276. Our findings are relevant for ordering phenomena in soft matter and alloys.

categories: cond-mat.soft cond-mat.mes-hall cond-mat.mtrl-sci nlin.AO physics.comp-ph


title: An Interdisciplinary Comparison of Sequence Modeling Methods for Next-Element Prediction

abstract: Data of sequential nature arise in many application domains in forms of, e.g. textual data, DNA sequences, and software execution traces. Different research disciplines have developed methods to learn sequence models from such datasets: (i) in the machine learning field methods such as (hidden) Markov models and recurrent neural networks have been developed and successfully applied to a wide-range of tasks, (ii) in process mining process discovery techniques aim to generate human-interpretable descriptive models, and (iii) in the grammar inference field the focus is on finding descriptive models in the form of formal grammars. Despite their different focuses, these fields share a common goal - learning a model that accurately describes the behavior in the underlying data. Those sequence models are generative, i.e, they can predict what elements are likely to occur after a given unfinished sequence. So far, these fields have developed mainly in isolation from each other and no comparison exists. This paper presents an interdisciplinary experimental evaluation that compares sequence modeling techniques on the task of next-element prediction on four real-life sequence datasets. The results indicate that machine learning techniques that generally have no aim at interpretability in terms of accuracy outperform techniques from the process mining and grammar inference fields that aim to yield interpretable models.


title: Compact G2 holonomy spaces from SU(3) structures

abstract: We construct novel classes of compact G2 spaces from lifting type IIA flux backgrounds with O6 planes. There exists an extension of IIA Calabi-Yau orientifolds for which some of the D6 branes (required to solve the RR tadpole) are dissolved in $F_2$ fluxes. The backreaction of these fluxes deforms the Calabi-Yau manifold into a specific class of SU(3)-structure manifolds. The lift to M-theory again defines compact G2 manifolds, which in case of toroidal orbifolds are a twisted generalisation of the Joyce construction. This observation also allows a clear identification of the moduli space of a warped compactification with fluxes. We provide a few explicit examples, of which some can be constructed from T-dualising known IIB orientifolds with fluxes. Finally we discuss supersymmetry breaking in this context and suggest that the purely geometric picture in M-theory could provide a simpler setting to address some of the consistency issues of moduli stabilisation and de Sitter uplifting.


title: A representation of joint moments of CUE characteristic polynomials in terms of Painleve functions

abstract: We establish a representation of the joint moments of the characteristic polynomial of a CUE random matrix and its derivative in terms of a solution of the sigma-Painleve V equation. The derivation involves the analysis of a formula for the joint moments in terms of a determinant of generalised Laguerre polynomials using the Riemann-Hilbert method. We use this connection with the sigma-Painleve V equation to derive explicit formulae for the joint moments and to show that in the large-matrix limit the joint moments are related to a solution of the sigma-Painleve III equation. Using the conformal block expansion of the tau-functions associated with the sigma-Painleve V and the sigma-Painleve III equations leads to general conjectures for the joint moments.


title: Globular clusters as a probe for Weyl Conformal Gravity

abstract: Eventual flattening of velocity dispersion profiles of some galactic globular clusters in the Milky Way cannot be explained in the framework of Newtonian gravity and hence in general theory of relativity in the weak field limit, without resorting to the occurrence of tidal effects. We explore the possibility of explaining such deviation from expected Keplerian fall-off in dispersion profiles within the context of Weyl conformal gravity. We choose a set of 20 globular clusters for which recent kinematic measurements are available. We model the globular clusters with approximate Hernquist mass profiles and choose a constant mass-to-light ratio throughout the cluster as the only free parameter in the model. Our analysis finds reasonable Weyl gravity fits to the observed dispersion profiles, exhibiting both Keplerian decline and eventual flattening, with acceptable mass-to-light ratios. We further recover a Tully-Fisher like scaling relation in globular clusters through Weyl gravity.


title: Aligning Very Small Parallel Corpora Using Cross-Lingual Word Embeddings and a Monogamy Objective

abstract: Count-based word alignment methods, such as the IBM models or fast-align, struggle on very small parallel corpora. We therefore present an alternative approach based on cross-lingual word embeddings (CLWEs), which are trained on purely monolingual data. Our main contribution is an unsupervised objective to adapt CLWEs to parallel corpora. In experiments on between 25 and 500 sentences, our method outperforms fast-align. We also show that our fine-tuning objective consistently improves a CLWE-only baseline.


title: Three-charge superstrata with internal excitations

abstract: We construct a new family of three-charge $\frac{1}{8}-$BPS smooth solutions that have the same charges as the supersymmetric D1D5P Black Hole and are non-invariant under rotations of the compact manifold. We work in type IIB string theory on $T^4$ and we show how the supergravity and BPS equations reduce to a linear system, arranged in two `layers' of partial differential equations. We then build two solutions of our system of equations: the first is a superdescendant three-charge solution, obtained by acting with rigid symmetries on a seed two-charge solution, on which we can perform a non-trivial check for the system of equations; the second is a new superstratum solution that has both internal and external excitations. We then describe the CFT heavy states dual to these new geometries.


title: Accurate gradient computations for shape optimization via discrete adjoints in CFD-related multiphysics problems

abstract: As more and more multiphysics effects are entering the field of CFD simulations, this raises the question how they can be accurately captured in gradient computations for shape optimization. The latter has been successfully enriched over the last years by the use of (discrete) adjoints. One can think of them as Lagrange multipliers to the flow field problem linked to an objective function that depends on quantities like pressure or momentums, and they will set also the framework for this paper. It is split into two main parts: First, we show how one can compute coupled discrete adjoints using automatic differentiation in an effective way that is still easily extendable for all kinds of other couplings. Second, we suppose that a valuable first application are so-called conjugate heat transfer problems which are gaining more and more interest from the automobile and aeronautics industry. Therefore we present an implementation for this capability within the open-source solver SU2 as well as for the generic adjoint computation algorithm.


title: On the existence and approximation of a dissipating feedback

abstract: Given a matrix $A\in \R^{n\times n}$ and a tall rectangular matrix $B \in \R^{n\times q}$, $q < n$, we consider the problem of making the pair $(A,B)$ dissipative, that is the determination of a {\it feedback} matrix $K \in \R^{q\times n}$ such that the field of values of $A-B K$ lies in the left half open complex plane. We review and expand classical results available in the literature on the existence and parameterization of the class of dissipating matrices, and we explore new matrix properties associated with the problem. In addition, we discuss various computational strategies for approximating the minimal Frobenius norm dissipating $K$.


title: Effective Feature Representation for Clinical Text Concept Extraction

abstract: Crucial information about the practice of healthcare is recorded only in free-form text, which creates an enormous opportunity for high-impact NLP. However, annotated healthcare datasets tend to be small and expensive to obtain, which raises the question of how to make maximally efficient uses of the available data. To this end, we develop an LSTM-CRF model for combining unsupervised word representations and hand-built feature representations derived from publicly available healthcare ontologies. We show that this combined model yields superior performance on five datasets of diverse kinds of healthcare text (clinical, social, scientific, commercial). Each involves the labeling of complex, multi-word spans that pick out different healthcare concepts. We also introduce a new labeled dataset for identifying the treatment relations between drugs and diseases.


title: Valley selective optical control of excitons in 2D semiconductors using Chiral metasurface

abstract: Recent advances in condensed matter physics have shown that the valley degree of freedom of electrons in 2D materials with hexagonal symmetry, such as graphene, h-BN, and TMDs, can be efficiently exploited, leading to the emergent field of valleytronics, which offers unique opportunities for efficient data transfer, computing and storage. The ability to couple the valley degree of freedom of electrons with light can further expand the ways one manipulate this degree of freedom, thus envisioning a new class of solid-state-photonic interfaces and devices. Besides this expansion of control of valley by light-waves, coupling of photons with valley-polarized electrons can dramatically expand the landscape of available optical responses, which may bring new means of controlling light in photonic devices. In this work we design such hybrid solid-state photonic metasurface integrating 2D TMD and photonic all-dielectric metasurface. While TMD is naturally endowed with the property of valley to optical-polarization coupling, the photonic metasurface is designed to produce chiral field which selectively couples to the valley degree of freedom of solid-state TMD component. We experimentally demonstrate that such coupling leads to controlled valley polarization due to the coupling of 2D materials with the chiral photonic metasurface. The measured emission from valley excitons in this hybrid system yields the preferential emission of specific helicity.

categories: physics.optics cond-mat.mtrl-sci cond-mat.other physics.app-ph


title: Light escape cones in local reference frames of Kerr-de Sitter black hole spacetimes and related black hole shadows

abstract: We construct the light escape cones of isotropic spot sources of radiation residing in special classes of reference frames in the Kerr-de Sitter (KdS) black hole spacetimes, namely, in the fundamental class of 'non-geodesic' locally non-rotating reference frames (LNRFs), and two classes of 'geodesic' frames, the radial geodesic frames (RGFs), both falling and escaping, and the frames related to the circular geodesic orbits (CGFs). We compare the cones constructed in a given position for the LNRFs, RGFs, and CGFs. We have shown that the photons locally counter-rotating relative to LNRFs with positive impact parameter and negative covariant energy are confined to the ergosphere region. Finally, we demonstrate that the light escaping cones govern the shadows of black holes located in front of a radiating screen, as seen by the observers in the considered frames. For shadows related to distant static observers the LNRFs are relevant.


title: Deep Generative Model with Beta Bernoulli Process for Modeling and Learning Confounding Factors

abstract: While deep representation learning has become increasingly capable of separating task-relevant representations from other confounding factors in the data, two significant challenges remain. First, there is often an unknown and potentially infinite number of confounding factors coinciding in the data. Second, not all of these factors are readily observable. In this paper, we present a deep conditional generative model that learns to disentangle a task-relevant representation from an unknown number of confounding factors that may grow infinitely. This is achieved by marrying the representational power of deep generative models with Bayesian non-parametric factor models, where a supervised deterministic encoder learns task-related representation and a probabilistic encoder with an Indian Buffet Process (IBP) learns the unknown number of unobservable confounding factors. We tested the presented model in two datasets: a handwritten digit dataset (MNIST) augmented with colored digits and a clinical ECG dataset with significant inter-subject variations and augmented with signal artifacts. These diverse data sets highlighted the ability of the presented model to grow with the complexity of the data and identify the absence or presence of unobserved confounding factors.


title: Efficient Collection of Connected Vehicles Data with Precision Guarantees

abstract: Connected vehicles disseminate detailed data, including their position and speed, at a very high frequency. Such data can be used for accurate real-time analysis, prediction and control of transportation systems. The outstanding challenge for such analysis is how to continuously collect and process extremely large volumes of data. To address this challenge, efficient collection of data is critical to prevent overburdening the communication systems and overreaching computational and memory capacity. In this work, we propose an efficient data collection scheme that selects and transmits only a small subset of data to alleviate data transmission burden. As a demonstration, we have used the proposed approach to select data points to be transmitted from 10,000 connected vehicles trips available in the Safety Pilot Model Deployment dataset. The presented results show that collection ratio can be as small as 0.05 depending on the required precision. Moreover, a simulation study was performed to evaluate the travel time estimation accuracy using the proposed data collection approach. Results show that the proposed data collection approach can significantly improve travel time estimation accuracy.


title: The UEA multivariate time series classification archive, 2018

abstract: In 2002, the UCR time series classification archive was first released with sixteen datasets. It gradually expanded, until 2015 when it increased in size from 45 datasets to 85 datasets. In October 2018 more datasets were added, bringing the total to 128. The new archive contains a wide range of problems, including variable length series, but it still only contains univariate time series classification problems. One of the motivations for introducing the archive was to encourage researchers to perform a more rigorous evaluation of newly proposed time series classification (TSC) algorithms. It has worked: most recent research into TSC uses all 85 datasets to evaluate algorithmic advances. Research into multivariate time series classification, where more than one series are associated with each class label, is in a position where univariate TSC research was a decade ago. Algorithms are evaluated using very few datasets and claims of improvement are not based on statistical comparisons. We aim to address this problem by forming the first iteration of the MTSC archive, to be hosted at the website www.timeseriesclassification.com. Like the univariate archive, this formulation was a collaborative effort between researchers at the University of East Anglia (UEA) and the University of California, Riverside (UCR). The 2018 vintage consists of 30 datasets with a wide range of cases, dimensions and series lengths. For this first iteration of the archive we format all data to be of equal length, include no series with missing data and provide train/test splits.


title: Large Tournament Games

abstract: We consider a stochastic tournament game in which each player is rewarded based on her rank in terms of the completion time of her own task and is subject to cost of effort. When players are homogeneous and the rewards are purely rank dependent, the equilibrium has a surprisingly explicit characterization, which allows us to conduct comparative statics and obtain explicit solution to several optimal reward design problems. In the general case when the players are heterogenous and payoffs are not purely rank dependent, we prove the existence, uniqueness and stability of the Nash equilibrium of the associated mean field game, and the existence of an approximate Nash equilibrium of the finite-player game. Our results have some potential economic implications; e.g., they lend support to government subsidies for R and D, assuming the profits to be made are substantial.


title: Void distributions reveal structural link between jammed packings and protein cores

abstract: Dense packing of hydrophobic residues in the cores of globular proteins determines their stability. Recently, we have shown that protein cores possess packing fraction $\phi \approx 0.56$, which is the same as dense, random packing of amino acid-shaped particles. In this article, we compare the structural properties of protein cores and jammed packings of amino acid-shaped particles in much greater depth by measuring their local and connected void regions. We find that the distributions of surface Voronoi cell volumes and local porosities obey similar statistics in both systems. We also measure the probability that accessible, connected void regions percolate as a function of the size of a spherical probe particle and show that both systems possess the same critical probe size. By measuring the critical exponent $\tau$ that characterizes the size distribution of connected void clusters at the onset of percolation, we show that void percolation in packings of amino acid-shaped particles and protein cores belong to the same universality class, which is different from that for void percolation in jammed sphere packings. We propose that the connected void regions of proteins are a defining feature of proteins and can be used to differentiate experimentally observed proteins from decoy structures that are generated using computational protein design software. This work emphasizes that jammed packings of amino acid-shaped particles can serve as structural and mechanical analogs of protein cores, and could therefore be useful in modeling the response of protein cores to cavity-expanding and -reducing mutations.


title: On Single-Channel Speech Enhancement and On Non-Linear Modulation-Domain Kalman Filtering

abstract: This report focuses on algorithms that perform single-channel speech enhancement. The author of this report uses modulation-domain Kalman filtering algorithms for speech enhancement, i.e. noise suppression and dereverberation, in [1], [2], [3], [4] and [5]. Modulation-domain Kalman filtering can be applied for both noise and late reverberation suppression and in [2], [1], [3] and [4], various model-based speech enhancement algorithms that perform modulation-domain Kalman filtering are designed, implemented and tested. The model-based enhancement algorithm in [2] estimates and tracks the speech phase. The short-time-Fourier-transform-based enhancement algorithm in [5] uses the active speech level estimator presented in [6]. This report describes how different algorithms perform speech enhancement and the algorithms discussed in this report are addressed to researchers interested in monaural speech enhancement. The algorithms are composed of different processing blocks and techniques [7]; understanding the implementation choices made during the system design is important because this provides insights that can assist the development of new algorithms. Index Terms - Speech enhancement, dereverberation, denoising, Kalman filter, minimum mean squared error estimation.


title: Predictive Modeling of Biomedical Signals Using Controlled Spatial Transformation

abstract: An important paradigm in smart health is developing diagnosis tools and monitoring a patient's heart activity through processing Electrocardiogram (ECG) signals is a key example, sue to high mortality rate of heart-related disease. However, current heart monitoring devices suffer from two important drawbacks: i) failure in capturing inter-patient variability, and ii) incapability of identifying heart abnormalities ahead of time to take effective preventive and therapeutic interventions. This paper proposed a novel predictive signal processing method to solve these issues. We propose a two-step classification framework for ECG signals, where a global classifier recognizes severe abnormalities by comparing the signal against a universal reference model. The seemingly normal signals are then passed through a personalized classifier, to recognize mild but informative signal morphology distortions. The key idea is to develop a novel deviation analysis based on a controlled nonlinear transformation to capture significant deviations of the signal towards any of predefined abnormality classes. Here, we embrace the proven but overlooked fact that certain features of ECG signals reflect underlying cardiac abnormalities before the occurrences of cardiac disease. The proposed method achieves a classification accuracy of 96.6% and provides a unique feature of predictive analysis by providing warnings before critical heart conditions. In particular, the chance of observing a severe problem (a red alarm) is raised by about 5% to 10% after observing a yellow alarm of the same type. Although we used this methodology to provide early precaution messages to elderly and high-risk heart-patients, the proposed method is general and applicable to similar bio-medical signal processing applications.


title: Manifold Learning of Four-dimensional Scanning Transmission Electron Microscopy

abstract: Four-dimensional scanning transmission electron microscopy (4D-STEM) of local atomic diffraction patterns is emerging as a powerful technique for probing intricate details of atomic structure and atomic electric fields. However, efficient processing and interpretation of large volumes of data remain challenging, especially for two-dimensional or light materials because the diffraction signal recorded on the pixelated arrays is weak. Here we employ data-driven manifold leaning approaches for straightforward visualization and exploration analysis of the 4D-STEM datasets, distilling real-space neighboring effects on atomically resolved deflection patterns from single-layer graphene, with single dopant atoms, as recorded on a pixelated detector. These extracted patterns relate to both individual atom sites and sublattice structures, effectively discriminating single dopant anomalies via multi-mode views. We believe manifold learning analysis will accelerate physics discoveries coupled between data-rich imaging mechanisms and materials such as ferroelectric, topological spin and van der Waals heterostructures.

categories: eess.IV cond-mat.mtrl-sci physics.data-an stat.ML


title: Spatial structure of the color field in the SU(3) flux tube

abstract: We report on the chromoelectric and chromomagnetic fields generated by a static quark-antiquark pair at zero temperature in pure gauge SU(3). From the spatial structure of chromoelectric field we extract its nonperturbative part and discuss its properties.


title: Tiling-based models of perimeter and area

abstract: We consider polygonal tilings of certain regions and use these to give intuitive definitions of tiling-based perimeter and area. We apply these definitions to rhombic tilings of Elnitsky polygons, computing sharp bounds and average values for perimeter tiles in convex centrally symmetric 2n-gons. These bounds and values have implications for the combinatorics of reduced decompositions of permutations. We also classify the permutations whose polygons gave minimal perimeter, defined in two different ways. We conclude by looking at some of these questions in the context of domino tilings, giving a recursive formula and generating function for one family, and describing a family of minimal-perimeter regions.


title: Regularizing effect of absorption terms in singular problems

abstract: We prove existence of solutions to problems whose model is $$\begin{cases} \displaystyle -\Delta_p u + u^q = \frac{f}{u^\gamma} & \text{in}\ \Omega, \newline u\ge0 &\text{in}\ \Omega,\newline u=0 &\text{on}\ \partial\Omega, \end{cases}$$ where $\Omega$ is an open bounded subset of $\mathbb{R}^N$ ($N\ge2$), $\Delta_p u$ is the $p$-laplacian operator for $1\le p <N$, $q>0$, $\gamma\ge 0$ and $f$ is a nonnegative function in $L^m(\Omega)$ for some $m\ge1$. In particular we analyze the regularizing effect produced by the absorption term in order to infer the existence of finite energy solutions in case $\gamma\le 1$. We also study uniqueness of these solutions as well as examples which show the optimality of the results. Finally, we find local $W^{1,p}$-solutions in case $\gamma>1$.


title: On the Deuring Polynomial for Drinfeld Modules in Legendre Form

abstract: We study a family $\psi^{\lambda}$ of $\mathbb F_q[T]$-Drinfeld modules, which is a natural analog of Legendre elliptic curves. We then find a surprising recurrence giving the corresponding Deuring polynomial $H_{p(T)}(\lambda)$ characterising supersingular Legendre Drinfeld modules $\psi^{\lambda}$ in characteristic $p(T)$.


title: The radio-bright accreting millisecond X-ray pulsar IGR J17591-2342

abstract: IGR J17591$-$2342 is a 527-Hz accreting millisecond X-ray pulsar that was discovered in outburst in 2018 August. In this paper, we present quasi-simultaneous radio and X-ray monitoring of this source during the early part of the outburst. IGR J17591$-$2342 is highly absorbed in X-rays, with an equivalent hydrogen absorption along the line of sight, $N_{\rm H}$, of $\approx$4.4$\times$10$^{22}$\,cm$^{-2}$, where the Galactic column density is expected to be $\approx$1 -- 2$\times$10$^{22}$\,cm$^{-2}$. The high absorption suggests that the source is either relatively distant ($>$6\,kpc), or the X-ray emission is strongly absorbed by material local to the system. Radio emission detected by the Australia Telescope Compact Array shows that, for a given X-ray luminosity and for distances greater than 3\,kpc, this source was exceptionally radio loud when compared to other accreting neutron stars in outburst ($L_{\rm X} > 10^{33}$\,erg\,s$^{-1}$). For most reasonable distances, IGR J17591$-$2342 appeared as radio luminous as actively accreting, stellar-mass black hole X-ray binaries.


title: Lattice Hydrodynamics

abstract: Using the combinatorics of two interpenetrating face centered cubic lattices together with the part of calculus naturally encoded in combinatorial topology, we construct from first principles a lattice model of 3D incompressible hydrodynamics on triply periodic three space. Actually the construction applies to every dimension, but has special duality features in dimension three.


title: Band dispersion of graphene with structural defects

abstract: We study the band dispersion of graphene with randomly distributed structural defects using two complementary methods, exact diagonalization of the tight-binding Hamiltonian and implementing a self-consistent T matrix approximation. We identify three distinct types of impurities resulting in qualitatively different spectra in the vicinity of the Dirac point. First, resonant impurities, such as vacancies or 585 defects, lead to stretching of the spectrum at the Dirac point with a finite density of localized states. This type of spectrum has been observed in epitaxial graphene by photoemission spectroscopy and discussed extensively in the literature. Second, nonresonant (weak) impurities, such as paired vacancies or Stone-Wales defects, do not stretch the spectrum but provide a line broadening that increases with energy. Finally, disorder that breaks sublattice symmetry, such as vacancies placed in only one sublattice, open a gap around the Dirac point and create an impurity band in the middle of this gap. We find good agreement between the results of the two methods and also with the experimentally measured spectra.


title: Applications of the moduli continuity method to log K-stable pairs

abstract: The 'moduli continuity method' permits an explicit algebraisation of the Gromov-Hausdorff compactification of K\"ahler-Einstein metrics on Fano manifolds in some fundamental examples. In this paper, we apply such method in the 'log setting' to describe explicitly some compact moduli spaces of K-polystable log Fano pairs. We focus on situations when the angle of singularities is perturbed in an interval sufficiently close to one, by considering constructions arising from Geometric Invariant Theory. More precisely, we discuss the cases of pairs given by cubic surfaces with anticanonical sections, and of projective space with non-Fano hypersurfaces, and we show ampleness of the CM line bundle on their good moduli space (in the sense of Alper). Finally, we introduce a conjecture relating K-stability (and degenerations) of log pairs formed by a fixed Fano variety and pluri-anticanonical sections to certain natural GIT quotients.


title: Structure refinement from 'Digital' Large Angle Convergent Beam Electron Diffraction Patterns

abstract: We use semi-automated data acquisition and processing to produce digital large angle CBED (D-LACBED) patterns. We demonstrate refinements of atomic coordinates and isotropic Debye-Waller factors for well-known materials using simulations produced with a neutral, spherical independent atom model. We find that atomic coordinate refinements in Al2O3 have sub-pm precision and accuracy. Isotropic DWFs are accurate for Cu, a simple fcc metal, but do not agree with X-ray measurements of GaAs or Al2O3. This indicates that models of bonding will be essential to fully interpret D-LACBED data.


title: SDRL: Interpretable and Data-efficient Deep Reinforcement Learning Leveraging Symbolic Planning

abstract: Deep reinforcement learning (DRL) has gained great success by learning directly from high-dimensional sensory inputs, yet is notorious for the lack of interpretability. Interpretability of the subtasks is critical in hierarchical decision-making as it increases the transparency of black-box-style DRL approach and helps the RL practitioners to understand the high-level behavior of the system better. In this paper, we introduce symbolic planning into DRL and propose a framework of Symbolic Deep Reinforcement Learning (SDRL) that can handle both high-dimensional sensory inputs and symbolic planning. The task-level interpretability is enabled by relating symbolic actions to options.This framework features a planner -- controller -- meta-controller architecture, which takes charge of subtask scheduling, data-driven subtask learning, and subtask evaluation, respectively. The three components cross-fertilize each other and eventually converge to an optimal symbolic plan along with the learned subtasks, bringing together the advantages of long-term planning capability with symbolic knowledge and end-to-end reinforcement learning directly from a high-dimensional sensory input. Experimental results validate the interpretability of subtasks, along with improved data efficiency compared with state-of-the-art approaches.


title: Proton-transfer Ferroelectricity / Multiferroicity in Rutile Oxyhydroxides

abstract: Oxyhydroxide minerals like FeOOH have been a research focus in geology for studying the Earth interior, and also in chemistry for studying oxygen electrocatalysis activity. In this paper we provide first-principles evidence of a new class of ferroelectrics or multiferroics among them:GaOOH,InOOH,CrOOH,FeOOH, which are earth-abundant minerals and have been experimentally verified to possess distorted rutile structures, are ferroelectric with considerable polarizations(up to 24 muC/cm2) and piezoelectric coefficients. Their atomic-thick layer may possess vertical polarization robust against depolarizing field due to the formation of O-H O bonds that can hardly be symmetrized. Moreover,CrOOH (guyanaite) is revealed to be a combination of high-Tc in-plane type-I multiferroics and vertical type-II multiferroics, which is strain-tunable and may render a desirable coupling between magnetism and ferroelectricity. Supported by experimental evidence on reversible conversion between metal oxyhydroxides and dioxides and their nice lattice match that renders convenient epitaxial growth, heterostructure composed of oxyhydroxides and prevalent metal dioxides (e.g. TiO2, SnO2 and CrO2) may be constructed for various applications like ferroelectric field-effect transistors and multiferroic tunneling junctions.


title: Network models for characterization of trabecular bone

abstract: Trabecular bone is a lightweight, compliant material organized as a web of struts and rods (trabeculae) that erode with age and the onset of bone diseases like osteoporosis, leading to increased fracture risk. The traditional diagnostic marker of osteoporosis, bone mineral density (BMD), has been shown in ex vivo experiments to correlate poorly with fracture resistance when considered on its own, while structural features in conjunction with BMD can explain more of the variation in trabecular bone strength. We develop a network-based model of trabecular bone by creating graphs from micro-computed tomography images of human bone, with weighted links representing trabeculae and nodes representing branch points. These graphs enable calculation of quantitative network metrics to characterize trabecular structure. We also create finite element models of the networks in which each link is represented by a beam, facilitating analysis of the mechanical response of the bone samples to simulated loading. We examine the structural and mechanical properties of trabecular bone at the scale of individual trabeculae (of order 0.1 mm) and at the scale of selected volumes of interest (approximately a few mm), referred to as VOIs. At the VOI scale, we find significant correlations between the stiffness of VOIs and 10 different structural metrics. Individually, the volume fraction of each VOI is most strongly correlated to the stiffness of the VOI. We use multiple linear regression to identify the smallest subset of variables needed to capture the variation in stiffness. In a linear fit, we find that node degree, weighted node degree, Z-orientation, weighted Z-orientation, trabecular spacing, link length, and the number of links are the structural metrics that are most significant (p < 0.05) in capturing the variation of stiffness in trabecular networks.


title: Supertranslations in Higher Dimensions Revisited

abstract: In this paper, we revisit the question of identifying Soft Graviton theorem in higher (even) dimensions with Ward identities associated with Asymptotic symmetries. Building on the prior work of \cite{strominger}, we compute, from first principles, the (asymptotic) charge associated to Supertranslation symmetry in higher even dimensions and show that (i) these charges are non-trivial, finite and (ii) the corresponding Ward identities are indeed the soft graviton theorems.


title: Unusual Ferroelectricity of Trans-Unitcell Ion-Displacement and Multiferroic Soliton in Sodium and Potassium Hydroxides

abstract: We show the first-principles evidence of a hitherto unreported type of ferroelectricity with ultra-long ion-displacement in sodium and potassium hydroxides. Even a small amount of proton vacancies can completely change the mode of proton-transfer from intra-unitcell to trans-unitcell, giving rise to multiferroic soliton with mobile magnetism and a tremendous polarization that can be orders of magnitude higher compared with most perovskite ferroelectrics. Their vertical polarizations of thin-film are robust against depolarizing field, rendering various designs of 2D ferroelectric field-transistors with non-destructive readout and ultra-high on/off ratio via sensing the switchable metallic/insulating state.


title: Latest results of the OPERA experiment on nu-tau appearance in the CNGS neutrino beam

abstract: OPERA is a long-baseline experiment designed to search for $\nu_{\mu}\to\nu_{\tau}$ oscillations in appearance mode. It was based at the INFN Gran Sasso laboratory (LNGS) and took data from 2008 to 2012 with the CNGS neutrino beam from CERN. After the discovery of $\nu_\tau$ appearance in 2015, with $5.1\sigma$ significance, the criteria to select $\nu_\tau$ candidates have been extended and a multivariate approach has been used for events identification. In this way the statistical uncertainty in the measurement of the oscillation parameters and of $\nu_\tau$ properties has been improved. Results are reported.


title: doFORC tool for calculating first-order reversal curve diagrams of noisy scattered data

abstract: First-order reversal curves (FORC) diagram method is one of the most successful characterization techniques used to characterize complex hysteretic phenomena not only in magnetism, but also in other areas of science like in ferroelectricity, geology, archeology, light-induced and pressure hysteresis in spin-transition materials, etc. Because the definition of the FORC diagram involves a second-order derivative, the main problem in their numerical calculation is that the derivative of a function for which only discrete noise-contaminated data values are available magnifies the noise that is inevitably present in measurements. In this paper we present doFORC tool for calculating FORC diagrams of noise scattered data. It can provide both a smooth approximation of the measured magnetization and all its partial derivatives. doFORC is a free, portable application working on various operating systems, with an easy to use graphical interface, with four regression methods implemented to obtain a smooth approximation of the data which may then be differentiated to obtain approximations for derivatives. In order to perform the diagnostics and goodness of fit doFORC computes residuals to characterize the difference between observed and predicted values, generalized cross-validation to measure the predictive performance, two information criteria to quantify the information that is lost by using an approximate model, and three degrees of freedom to compare different amounts of smoothing being performed by different smoothing methods. Based on these doFORC can perform automatic smoothing parameter selection.


title: An Evolutionary Algorithm with Crossover and Mutation for Model-Based Clustering

abstract: An evolutionary algorithm (EA) is developed as an alternative to the EM algorithm for parameter estimation in model-based clustering. This EA facilitates a different search of the fitness landscape, i.e., the likelihood surface, utilizing both crossover and mutation. Furthermore, this EA represents an efficient approach to "hard" model-based clustering and so it can be viewed as a sort of generalization of the k-means algorithm, which is itself equivalent to a restricted Gaussian mixture model. The EA is illustrated on several datasets, and its performance is compared to other hard clustering approaches and model-based clustering via the EM algorithm.