title: Privacy of Dependent Users Against Statistical Matching

abstract: Modern applications significantly enhance user experience by adapting to each user's individual condition and/or preferences. While this adaptation can greatly improve a user's experience or be essential for the application to work, the exposure of user data to the application presents a significant privacy threat to the users\textemdash even when the traces are anonymized\textemdash since the statistical matching of an anonymized trace to prior user behavior can identify a user and their habits. Because of the current and growing algorithmic and computational capabilities of adversaries, provable privacy guarantees as a function of the degree of anonymization and obfuscation of the traces are necessary. Our previous work has established the requirements on anonymization and obfuscation in the case that data traces are independent between users. However, the data traces of different users will be dependent in many applications, and an adversary can potentially exploit such. In this paper, we consider the impact of dependency between user traces on their privacy. First, we demonstrate that the adversary can readily identify the association graph of the obfuscated and anonymized version of the data, revealing which user data traces are dependent. Next, we demonstrate that the adversary can use this association graph to break user privacy with significantly shorter traces than in the case of independent users, and that obfuscating data traces independently across users is often insufficient to remedy such leakage. Finally, we discuss how users can improve privacy by employing joint obfuscation that removes or reduces the data dependency.


title: Water/Cosolvent Attraction Induced Phase Separation: a Molecular Picture of Cononsolvency

abstract: Cononsolvency is a phenomenon for which the solubility of a macromolecule decreases or even vanishes in the mixture of two good solvents. Although it has been widely applied in physicochemical, green chemical and pharmaceutical industry, its origin is still under active debate. Here, by using combined neutron total scattering, deuterium-labelling and all-atom molecular dynamic simulations, we demonstrated that it is the strong water/cosolvent attraction that leads to the cononsolvency. The combined approach presented here has opened a new route for investigating the most probable all-atom structure in macromolecular solutions and the thermodynamic origin of solubilities.


title: Superradiant Instabilities in a Class of Scalar-Tensor Horndeski Theory

abstract: We study the superradiance effect in a class of scalar-tensor Horndeski theory. We first study the dynamics of a massive charged scalar wave scattered off the horizon of a Reissner-Nordstr\"om black hole which except its canonical coupling to gravity it is also coupled kinetically to curvature. We find that a trapping potential is formed outside the horizon of a Reissner-Nordstr\"om black hole, due to this coupling, and as the strength of the new coupling is increased, the scattered wave is superradiantly amplified, resulting to the instability of the Reissner-Nordstr\"om spacetime. We then considered the backreacting effect of the scalar field coupled to curvature interacting with the background metric and we study the superradiant effect and find the superradiance conditions of a massive charged wave scattered off the horizon of a Horndeski black hole.


title: Pulsar Timing Constraints on the Fermi Massive Black-Hole Binary Blazar Population

abstract: Blazars are a sub-population of quasars whose jets are nearly aligned with the line-of-sight, which tend to exhibit multi-wavelength variability on a variety of timescales. Quasi-periodic variability on year-like timescales has been detected in a number of bright sources, and has been connected to the orbital motion of a putative massive black hole binary. If this were indeed the case, those blazar binaries would contribute to the nanohertz gravitational-wave stochastic background. We test the binary hypothesis for the blazar population observed by the \textit{Fermi} Gamma-Ray Space Telescope, which consists of BL Lacertae objects and flat-spectrum radio quasars. Using mock populations informed by the luminosity functions for BL Lacertae objects and flat-spectrum radio quasars with redshifts $z \le 2$, we calculate the expected gravitational wave background and compare it to recent pulsar timing array upper limits. The two are consistent only if a fraction $\lesssim 10^{-3}$ of blazars hosts a binary with orbital periods $<5$ years. We therefore conclude that binarity cannot significantly explain year-like quasi-periodicity in blazars.


title: Formation of LISA Black Hole Binaries in Merging Dwarf Galaxies: the Imprint of Dark Matter

abstract: Theoretical models for the expected merger rates of intermediate-mass black holes (IMBHs) are vital for planned gravitational-wave detection experiments such as the Laser Interferometer Space Antenna (LISA). Using collisionless $N$-body simulations of dwarf galaxy (DG) mergers, we examine how the orbital decay of IMBHs and the efficiency of IMBH binary formation depend on the central dark matter (DM) density profile of the merging DGs. Specifically, we explore various asymptotic inner slopes $\gamma$ of the DG's DM density distribution, ranging from steep cusps ($\gamma=1$) to shallower density profiles ($\gamma<1$), motivated by well-known baryonic-feedback effects as well as by DM models that differ from cold DM at the scales of DGs. We find that the inner DM slope is crucial for the formation (or lack thereof) of an IMBH binary; only mergers between DGs with cuspy DM profiles ($\gamma=1$) are favourable to forming a hard IMBH binary, whereas when $\gamma<1$ the IMBHs stall at a separation of 50-100 pc. Consequently, the rate of LISA signals from IMBH coalescence will be determined by the fraction of DGs with a cuspy DM profile. Conversely, the LISA event rates at IMBH mass scales offer in principle a novel way to place constraints on the inner structure of DM halos in DGs and address the core-cusp controversy. We also show that, with spatial resolutions of $\sim$0.1 kpc, as often adopted in cosmological simulations, all IMBHs stall, independent of $\gamma$. This suggests caution in employing cosmological simulations of galaxy formation to study BH dynamics in DGs.


title: The role of multiple images and model priors in measuring $H_0$ from supernova Refsdal in galaxy cluster MACS J1149.5+2223

abstract: Multiple image gravitational lensing systems with measured time delays provide a promising one-step method for determining $H_0$. MACS J1149, which lenses SN Refsdal into a quad S1-S4, and two other widely separated images, SX and SY, is a perfect candidate. If time delays are pinned down, the remaining uncertainty arises from the mass distribution in the lens. In MACS J1149, the mass in the relevant lens plane region can be constrained by (i) many multiple images, (ii) the mass of the galaxy splitting S1-S4 (which, we show, is correlated with $H_0$), (iii) magnification of SX (also correlated with $H_0$), and (iv) prior assumptions on the mass distribution. Our goal is not to estimate $H_0$, but to understand its error budget, i.e., estimate uncertainties associated with each of these constraints. Using multiple image positions alone, yields very large uncertainty, despite the fact that the position of SX is recovered to within $\!\le\!0.036$" (rms $\!\le\!0.36$") by GRALE lens inversion. Fixing the mass of the galaxy that splits S1-S4 reduces $1\sigma$ uncertainties to $\sim 23\%$, while fixing the magnification of SX yields $1\sigma$ uncertainties of $32\%$. We conclude that smaller uncertainties, of order few percent, are a consequence of imposing prior assumptions on the shapes of the galaxy and cluster mass distributions, which may or may not apply in a highly non-equilibrium environment of a merging cluster. We propose that if a measurement of $H_0$ is to be considered reliable, it must be supported by a wide range of lens inversion methods.


title: The highly polarized dusty emission core of Cygnus A

abstract: We report the detection of linearly polarized emission at 53 and 89 $\mu$m, from the radio-loud active galactic nucleus (AGN) Cygnus A using HAWC+ onboard SOFIA. We measure a highly polarized core of $11\pm3$% and $9\pm2$% with a position angle (P.A.) of polarization of $43\pm8^{\circ}$ and $39\pm7^{\circ}$ at 53 and 89 $\mu$m, respectively. We find (1) a synchrotron dominated core with a flat spectrum ($+0.21\pm0.05$) and a turn-over at $543\pm120$ $\mu$m, which implies synchrotron emission is insignificant in the infrared (IR), and (2) a $2-500$ $\mu$m bump peaking at $\sim40$ $\mu$m described by a blackbody component with color temperature of $107\pm9$ K. The polarized SED has the same shape as the IR bump of the total flux SED. We observe a change in the P.A. of polarization of $\sim20^{\circ}$ from 2 to 89 $\mu$m, which suggests a change of polarization mechanisms. The ultraviolet, optical and near-IR polarization has been convincingly attributed to scattering by polar dust, consistent with the usual torus scenario, though this scattered component can only be directly observed from the core in the near-IR. By contrast, the gradual rotation by $\sim20^{\circ}$ towards the far-IR, and the near-perfect match between the total and polarized IR bumps, indicate that dust emission from aligned dust grains becomes dominant at $10-100$ $\mu$m, with a large polarization of 10% at a nearly constant P.A. This result suggests that a coherent dusty and magnetic field structure dominates the $10-100$ $\mu$m emission around the AGN.


title: Self-regulating jets during the Common Envelope phase

abstract: Jets launched from a compact object (CO) during a common envelope (CE) may play a key role in the evolution of the system, and may also be an efficient removal channel for its material. In this work we study, through a large set of three-dimensional hydrodynamic simulations, the effects that jets launched from either a black hole (BH) or a neutron star (NS) have during a CE phase. The jet is self-consistently powered by a fraction ($\eta$) of the mass accretion rate that reaches the CO. For low mass accretion efficiencies ($\eta$<0.1\%), the jet is not able to drill through the material accreting onto the CO and forms a bulge around it. For higher efficiencies, the jet is able to drill through the bulge and in intermediate efficiencies ($\eta \sim$1-5\%) may produce an oscillating mass accretion rate behavior. We find that the jet may deposit enough energy to unbind the outer layers of the CE. The self-regulated jets present variability in size and orientation while their cocoons expand smoothly over the CE. The mass accretion rate initially decreases due to the ram pressure of the cocoon. If the launched jet is not able to drill through the bulge, the accretion is such that it may convert a NS into a BH within a decade or double the mass of a BH in a few years. If a jet is present, it deposits enough energy to unbind the outer layers of the CE and the system evolves into a grazing envelope, configuration.


title: Higher-Order Topology, Monopole Nodal Lines, and the Origin of Large Fermi Arcs in Transition Metal Dichalcogenides XTe$_2$ (X=Mo,W)

abstract: In recent years, transition metal dichalcogenides (TMDs) have garnered great interest as topological materials -- monolayers of centrosymmetric $\beta$-phase TMDs have been identified as 2D topological insulators (TIs), and bulk crystals of noncentrosymmetric $\gamma$-phase MoTe$_2$ and WTe$_2$ have been identified as type-II Weyl semimetals. However, ARPES and STM probes of these TMDs have revealed huge, "arc-like" surface states that overwhelm, and are sometimes mistaken for, the much smaller topological surface Fermi arcs of bulk type-II Weyl points. In this letter, we use first-principles calculations and (nested) Wilson loops to analyze the bulk and surface electronic structure of both $\beta$- and $\gamma$-MoTe$_2$, finding that $\beta$-MoTe$_2$ ($\gamma$-MoTe$_2$ gapped with symmetry-preserving distortion) is an inversion-symmetry-indicated $\mathbb{Z}_{4}$-nontrivial ($noncentrosymmetric, non$-$symmetry$-$indicated$) higher-order TI (HOTI) driven by double band inversion. Both structural phases of MoTe$_2$ exhibit the same surface features as WTe$_2$, revealing that the large Fermi arcs are in fact not topologically trivial, but are rather the characteristic split and gapped fourfold surface states of a HOTI. We also show that, when the effects of SOC are neglected, $\beta$-MoTe$_2$ is a nodal-line semimetal with $\mathbb{Z}_{2}$-nontrivial monopole nodal lines (MNLSM). This finding confirms that MNLSMs driven by double band inversion are the weak-SOC limit of HOTIs, implying that MNLSMs are higher-order topological $semimetals$ with flat-band-like hinge states, which we find to originate from the corner modes of 2D "fragile" TIs.


title: Physical Bayesian modelling of the non-linear matter distribution: new insights into the Nearby Universe

abstract: Accurate analyses of present and next-generation galaxy surveys require new ways to handle effects of non-linear gravitational structure formation in data. To address these needs we present an extension of our previously developed algorithm for Bayesian Origin Reconstruction from Galaxies to analyse matter clustering at non-linear scales in observations. This is achieved by incorporating a numerical particle mesh model of structure formation into our Bayesian inference framework. The algorithm simultaneously infers the 3D primordial matter fluctuations from which present non-linear observations formed and provides reconstructions of velocity fields and structure formation histories. The physical forward modelling approach automatically accounts for non-Gaussian features in evolved matter density fields and addresses the redshift space distortion problem associated with peculiar motions of galaxies. Our algorithm employs a hierarchical Bayes approach to jointly account for observational effects, such as galaxy biases, selection effects, and noise. Corresponding parameters are marginalized out via a sophisticated Markov Chain Monte Carlo approach relying on a combination of a multiple block sampling framework and a Hamiltonian Monte Carlo sampler. We demonstrate the performance of the method by applying it to the 2M++ galaxy compilation, tracing the matter distribution of the Nearby Universe. We show accurate and detailed inferences of the 3D non-linear dark matter distribution of the Nearby Universe. As exemplified in the case of the Coma cluster, we provide mass estimates that are compatible with those obtained from weak lensing and X-ray observations. For the first time, we reconstruct the vorticity of the non-linear velocity field from observations. In summary, our method provides plausible and detailed inferences of dark matter and velocity fields of our cosmic neighbourhood.


title: Globular Clusters Formed within Dark Halos I: present-day abundance, distribution and kinematics

abstract: We explore a scenario where metal poor globular clusters (GCs) form at the centres of their own dark matter halos in the early universe before reionization. This hypothesis leads to predictions about the abundance, distribution and kinematics of GCs today that we explore using cosmological N-body simulations and analytical modelling. We find that selecting the massive tail of collapsed objects at $z\gtrsim 9$ as GC formation sites leads to four main predictions: i) a highly clustered population of GCs around galaxies today, ii) a natural scaling between number of GCs and halo virial mass that follows roughly the observed trend, iii) a very low number of free floating GCs outside massive halos and iv) GCs should be embedded within massive and extended dark matter (sub)halos. We find that the strongest constraint to the model is given by the combination of (i) and (ii): a mass cut to tagged GCs halos which accounts for the number density of metal poor GCs today predicts a radial distribution that is too extended compared to recent observations. On the other hand, a mass cut sufficient to match the observed half number radius could only explain 60% of the metal poor population. In all cases, observations favour early redshifts for GC formation ($z\geq 15$).


title: Searching for Decaying and Annihilating Dark Matter with Line Intensity Mapping

abstract: The purpose of line-intensity mapping (IM), an emerging tool for extragalactic astronomy and cosmology, is to measure the integrated emission along the line of sight from spectral lines emitted from galaxies and the intergalactic medium. The observed frequency of the line then provides a distance determination allowing the three-dimensional distribution of the emitters to be mapped. Here we discuss the possibility to use these measurements to seek radiative decays or annihilations from dark-matter particles. The photons from monoenergetic decays will be correlated with the mass distribution, which can be determined from galaxy surveys, weak-lensing surveys, or the IM mapping experiments themselves. We discuss how to seek this cross-correlation and then estimate the sensitivity of various IM experiments in the dark-matter mass-lifetime parameter space. We find prospects for improvements of ten orders of magnitude in sensitivity to decaying/annihilating dark matter in the frequency bands targeted for IM experiments.


title: Dissecting the roles of mass and environment quenching in galaxy evolution with EAGLE

abstract: We exploit the pioneering cosmological hydrodynamical simulation, EAGLE, to study how the connection between halo mass (M_halo), stellar mass (M*) and star-formation rate (SFR) evolves across redshift. Using Principal Component Analysis we identify the key axes of correlation between these physical quantities, for the full galaxy sample and split by satellite/central and low/high halo mass. The first principal component of the z=0 EAGLE galaxy population is a positive correlation between M_halo, M* and SFR. This component is particularly dominant for central galaxies in low mass haloes. The second principal component, most significant in high mass haloes, is a negative correlation between M_halo and SFR, indicative of environmental quenching. For galaxies above M*~10^10M_solar, however, the SFR is seen to decouple from the M_halo-M* correlation; this result is found to be independent of environment, suggesting that mass quenching effects are also in operation. We find extremely good agreement between the EAGLE principal components and those of SDSS galaxies; this lends confidence to our conclusions. Extending our study to EAGLE galaxies in the range z=0-4, we find that, although the relative numbers of galaxies in the different subsamples change, their principal components do not change significantly with redshift. This indicates that the physical processes that govern the evolution of galaxies within their dark matter haloes act similarly throughout cosmic time. Finally, we present halo occupation distribution model fits to EAGLE galaxies and show that one flexible 6-parameter functional form is capable of fitting a wide range of different mass- and SFR-selected subsamples.


title: The Star Formation Efficiency per Free Fall Time in Nearby Galaxies

abstract: We estimate the star formation efficiency per gravitational free fall time, $\epsilon_{\rm ff}$, from observations of nearby galaxies with resolution matched to the typical size of a Giant Molecular Cloud. This quantity, $\epsilon_{\rm ff}$, is theoretically important but so far has only been measured for Milky Way clouds or inferred indirectly in a few other galaxies. Using new, high resolution CO imaging from the PHANGS-ALMA survey, we estimate the gravitational free-fall time at 60 to 120 pc resolution, and contrast this with the local molecular gas depletion time to estimate $\epsilon_{\rm ff}$. Assuming a constant thickness of the molecular gas layer ($H = 100$ pc) across the whole sample, the median value of $\epsilon_{\rm ff}$ in our sample is $0.7\%$. We find a mild scale-dependence, with higher $\epsilon_{\rm ff}$ measured at coarser resolution. Individual galaxies show different values of $\epsilon_{\rm ff}$, with the median $\epsilon_{\rm ff}$ ranging from $0.3\%$ to $2.6\%$. We find the highest $\epsilon_{\rm ff}$ in our lowest mass targets, reflecting both long free-fall times and short depletion times, though we caution that both measurements are subject to biases in low mass galaxies. We estimate the key systematic uncertainties, and show the dominant uncertainty to be the estimated line-of-sight depth through the molecular gas layer and the choice of star formation tracers.


title: Scalar singlet dark matter in non-standard cosmologies

abstract: We study production of dark matter (DM) in models with a non-standard expansion history. We consider both freeze-out and freeze-in mechanisms for producing the observed DM abundance in a model where the DM consists of scalar singlet particles coupled to the Standard Model sector via the Higgs portal. We show that a non-standard expansion phase can lead to a significant change in the DM abundance and therefore to observational ramifications. For example, for DM freeze-in the required portal coupling can be much larger, whereas for DM freeze-out much smaller values become allowed. We evaluate the relevant constraints and discuss prospects for direct detection of such DM.


title: Quantum localization bounds Trotter errors in digital quantum simulation

abstract: A fundamental challenge in digital quantum simulation (DQS) is the control of inherent errors. These appear when discretizing the time evolution generated by the Hamiltonian of a quantum many-body system as a sequence of quantum gates, called Trotterization. Here, we show that quantum localization-by constraining the time evolution through quantum interference-strongly bounds these errors for local observables. Consequently, for generic quantum many-body Hamiltonians, Trotter errors can become independent of system size and total simulation time. For local observables, DQS is thus intrinsically much more robust than what one might expect from known error bounds on the global many-body wave function. This robustness is characterized by a sharp threshold as a function of the Trotter step size. The threshold separates a regular region with controllable Trotter errors, where the system exhibits localization in the space of eigenstates of the time-evolution operator, from a quantum chaotic regime where the trajectory is quickly scrambled throughout the entire Hilbert space. Our findings show that DQS with comparatively large Trotter steps can retain controlled Trotter errors for local observables. It is thus possible to reduce the number of quantum gate operations required to represent the desired time evolution faithfully, thereby mitigating the effects of imperfect individual gate operations


title: One-loop quantum gravitational backreaction on the local Hubble rate

abstract: We determine corrections to the Hubble rate due to graviton loops in a cosmological background spacetime of constant deceleration parameter. The corrections are gauge-invariant, based on a recent proposal for all-order gauge-invariant observables in perturbative quantum gravity. We find explicit expressions for the cases of matter- and radiation-dominated eras and slow-roll inflation with vanishing second slow-roll parameter. Interestingly, in the latter case the corrections can be described by a quantum-corrected first slow-roll parameter, which brings the spacetime closer to de Sitter space.


title: The 3CR $Chandra$ extragalactic survey at 1.0$<z<$1.5

abstract: The aim of this paper is to present an analysis of newly acquired X-ray observations of 16 extragalactic radio sources, listed in the Third Cambridge Revised (3CR) catalog, and not previously observed by Chandra. Observations were performed during Chandra Cycle 17, extending X-ray coverage for the 3CR extragalactic catalog up to $z$=1.5. Among the 16 targets, two lie at $z<$0.5 (i.e., 3CR27, at $z$=0.184 and 3CR69, at $z$=0.458), all the remaining 14 have redshifts between 1.0 and 1.5. In the current sample there are three compact steep spectrum (CSS) sources, three quasars and an FRI radio galaxy, while the other nine are FRII radio galaxies. All radio sources have an X-ray counterpart. We measured nuclear X-ray fluxes as well as X-ray emission associated with radio jet knots, hotspots or lobes in three energy bands: soft (0.5-1 keV), medium (1-2 keV) and hard (2-7 keV). We also performed standard X-ray spectral analysis for the four brightest nuclei. We discovered X-ray emission associated with: the radio lobe of 3CR124; a hotspot of the quasar 3CR220.2; another hotspot of the radio galaxy 3CR238; and the jet knot of 3CR297. We also detected extended X-ray emission around the nuclear region of 3CR124 and 3CR297 on scales of several tens of kpc. Finally, we present an update on the X-ray observations performed with Chandra and XMM-Newton on the entire 3CR extragalactic catalog.


title: Far-IR dust properties of highly dust obscured AGNs from the AKARI and WISE all-sky surveys

abstract: The combination of the AKARI and WISE infrared all-sky surveys provides an unique opportunity to identify and characterize the most highly dust obscured AGNs in the universe. Dust-obscured AGNs are not easily detectable and potentially underrepresented in extragalactic surveys due to their high optical extinction, but are readily found in the WISE catalog due to their extremely red mid-IR colors. Combining these surveys with photometry from Pan-STARRS and Herschel, we use SED modeling to characterize the extinction and dust properties of these AGNs. From mid-IR WISE colors, we are able to compute bolometric corrections to AGN luminosities. Using AKARI's far-IR wavelength photometry and broadband AGN/galaxy spectral templates, we estimate AGN dust mass and temperature using simple analytic models with 3-4 parameters. Even without spectroscopic data, we can determine a number of AGN dust properties only using SED analysis. These methods, combined with the abundance of archival photometric data publically available, will be valuable for large-scale studies of dusty, IR-luminous AGNs.


title: A Parsec-scale Bipolar H$_2$ Outflow in the Massive Star Forming Infrared Dark Cloud Core MSXDC G053.11+00.05 MM1

abstract: We present a parsec-scale molecular hydrogen (H$_2$ 1-0 S(1) at 2.12~\micron) outflow discovered from the UKIRT Widefield Infrared Survey for H$_2$. The outflow is located in the infrared dark cloud core MSXDC G053.11+00.05 MM1 at 1.7 kpc and likely associated with two young stellar objects (YSOs) at the center. The overall morphology of the outflow is bipolar along the NE-SW direction with a brighter lobe to the southwest, but the detailed structure consists of several flows and knots. With the total length of $\sim$1 pc, the outflow luminosity is fairly high with $L_{\rm H_{2}} > 6~L_{\sun}$, implying a massive outflow-driving YSO if the entire outflow is driven by a single source. The two putative driving sources, located at the outflow center, show photometric variability of $\gtrsim$1 mag in {\it H}- and {\it K}-bands. This, with their early evolutionary stage from spectral energy distribution (SED) fitting, indicates that both are capable of ejecting outflows and may be eruptive variable YSOs. The YSO masses inferred from SED fitting are $\sim$10~$M_{\sun}$ and $\sim$5~$M_{\sun}$, suggesting the association of the outflow with massive YSOs. The geometrical morphology of the outflow is well explained by the lower mass YSO by assuming a single source origin, but without kinematic information, the contribution from the higher mass YSO cannot be ruled out. Considering star formation process by fragmentation of a high-mass core into several lower mass stars, we also suggest the possible presence of another, yet-undetected driving source deeply embedded in the core.


title: A NUMA-Aware Provably-Efficient Task-Parallel Platform Based on the Work-First Principle

abstract: Task parallelism is designed to simplify the task of parallel programming. When executing a task parallel program on modern NUMA architectures, it can fail to scale due to the phenomenon called work inflation, where the overall processing time that multiple cores spend on doing useful work is higher compared to the time required to do the same amount of work on one core, due to effects experienced only during parallel executions such as additional cache misses, remote memory accesses, and memory bandwidth issues. It's possible to mitigate work inflation by co-locating the computation with the data, but this is nontrivial to do with task parallel programs. First, by design, the scheduling for task parallel programs is automated, giving the user little control over where the computation is performed. Second, the platforms tend to employ work stealing, which provides strong theoretical guarantees, but its randomized protocol for load balancing does not discern between work items that are far away versus ones that are closer. In this work, we propose NUMA-WS, a NUMA-aware task parallel platform engineering based on the work-first principle. By abiding by the work-first principle, we are able to obtain a platform that is work efficient, provides the same theoretical guarantees as the classic work stealing scheduler, and mitigates work inflation. Furthermore, we implemented a prototype platform by modifying Intel's Cilk Plus runtime system and empirically demonstrate that the resulting system is work efficient and scalable.


title: Critical Kaluza-Klein black holes and black strings in D = 10

abstract: We construct static vacuum localized black holes and non-uniform black strings in ten spacetime dimensions, where one of the dimension is compactified on circle. We study the phase diagram of black objects with these boundary conditions, especially near the critical point where localized black holes and non-uniform black strings merge. Remarkably, we find that the merger happens at a cusp in the phase diagram. We verify that the critical geometry is controlled by a Ricci-flat double cone as previously predicted. However, unlike the lower dimensional cases, we find that physical quantities approach to their critical values according to a power law plus a logarithmic correction. We extract the critical exponents and find very good agreement with the predictions from the double cone geometry. According to holography, localized black holes and black strings are dual to thermal states of (1 + 1)-dimensional SU($N$) maximal Super-Yang Mills theory compactified on a circle; we recover and extend the details of the (recently found) 1st order phase transition in this system from the gravity side.


title: Thermal fluctuations in antiferromagnetic nanostructures

abstract: A theoretical model is developed that can accurately analyze the effects of thermal fluctuations in antiferromagnetic (AFM) nano-particles. The approach is based on Fourier series representation of the random effective field with cut-off frequencies of physical origin at low and high limits while satisfying the fluctuation-dissipation theorem at the same time. When coupled with the formalism of a Langevin dynamical equation, it can describe the stochastic N\'eel vector dynamics with the AFM parameters, circumventing the arbitrariness of the commonly used treatments in the micro-magnetic simulations. Subsequent application of the model to spontaneous N\'eel vector switching provides a thermal stability analysis of the AFM states. The numerical simulation shows that the AFM states are much less prone to the thermally induced accidental flips than the ferromagnetic counterparts, suggesting a longer retention time for the former.


title: Stability in the high-dimensional cohomology of congruence subgroups

abstract: We prove a representation stability result for the codimension-one cohomology of the level three congruence subgroup of $\mathbf{SL}_n(\mathbb{Z})$. This is a special case of a question of Church-Farb-Putman which we make more precise. Our methods involve proving several finiteness properties of the Steinberg module for the group $\mathbf{SL}_n(K)$ for $K$ a field. This also lets us give a new proof of Ash-Putman-Sam's homological vanishing theorem for the Steinberg module. We also prove an integral refinement of Church-Putman's homological vanishing theorem for the Steinberg module for the group $\mathbf{SL}_n(\mathbb{Z})$.


title: FC-groups with finitely many automorphism orbits

abstract: Let $G$ be a group. The orbits of the natural action of $Aut(G)$ on $G$ are called "automorphism orbits" of $G$, and the number of automorphism orbits of $G$ is denoted by $\omega(G)$. In this paper we prove that if $G$ is an FC-group with finitely many automorphism orbits, then the derived subgroup $G'$ is finite and $G$ admits a decomposition $G = Tor(G) \times D$, where $Tor(G)$ is the torsion subgroup of $G$ and $D$ is a divisible characteristic subgroup of $Z(G)$. We also show that if $G$ is an infinite FC-group with $\omega(G) \leqslant 8$, then either $G$ is soluble or $G \cong A_5 \times H$, where $H$ is an infinite abelian group with $\omega(H)=2$. Moreover, we describe the structure of the infinite non-soluble FC-groups with at most eleven automorphism orbits.


title: 3D mapping of the neutral X-ray absorption in the local interstellar medium: The Gaia and XMM-Newton synergy

abstract: We present a three-dimensional map of the hydrogen density distribution in the Galactic interstellar medium. The hydrogen equivalent column densities were obtained from the Exploring the X-ray Transient and variable Sky project ({\sc EXTraS}) which provides equivalent $N_{\rm H}$ values from X-ray spectral fits of observations within the {\it XMM-Newton} Data Release. {\sc EXTraS} include multiple fits for each source, allowing an accurate determination of the equivalent column densities, which depends on the continuum modeling of the spectra. A cross-correlation between the {\sc EXTraS} catalogue and the first {\it Gaia} Data Release was performed in order to obtain accurate parallax and distance measurements. We use a Bayesian method explained in \citet{rez17} in order to predict the most probable distribution of the density at any arbitrary point, even for lines of sight along which there are no initial observation. The resulting map shows small-scale density structures which can not been modeled by using analytic density profiles. In this paper we present a proof of concept of the kind of science possible with the synergy of these catalogs. However, given the systematic uncertainties connected to the source identification and to the dependence of $N_{\rm H}$ on the spectral model, the present maps should be considered qualitatively at this point.


title: Killing Horizons: Negative Temperatures and Entropy Super-Additivity

abstract: Many discussions in the literature of spacetimes with more than one Killing horizon note that some horizons have positive and some have negative surface gravities, but assign to all a positive temperature. However, the first law of thermodynamics then takes a non-standard form. We show that if one regards the Christodoulou and Ruffini formula for the total energy or enthalpy as defining the Gibbs surface, then the rules of Gibbsian thermodynamics imply that negative temperatures arise inevitably on inner horizons, as does the conventional form of the first law. We provide many new examples of this phenomenon, including black holes in STU supergravity. We also give a discussion of left and right temperatures and entropies, and show that both the left and right temperatures are non-negative. The left-hand sector contributes exactly half the total energy of the system, and the right-hand sector contributes the other half. Both the sectors satisfy conventional first laws and Smarr formulae. For spacetimes with a positive cosmological constant, the cosmological horizon is naturally assigned a negative Gibbsian temperature. We also explore entropy-product formulae and a novel entropy-inversion formula, and we use them to test whether the entropy is a super-additive function of the extensive variables. We find that super-additivity is typically satisfied, but we find a counterexample for dyonic Kaluza-Klein black holes.

categories: hep-th cond-mat.stat-mech gr-qc quant-ph


title: A generalized Newton iteration for computing the solution of the inverse Henderson problem

abstract: We develop a generalized Newton scheme IHNC for the construction of effective pair potentials for systems of interacting point-like particles.The construction is made in such a way that the distribution of the particles matches a given radial distribution function. The IHNC iteration uses the hypernetted-chain integral equation for an approximate evaluation of the inverse of the Jacobian of the forward operator. In contrast to the full Newton method realized in the Inverse Monte Carlo (IMC) scheme, the IHNC algorithm requires only a single molecular dynamics computation of the radial distribution function per iteration step, and no further expensive cross-correlations. Numerical experiments are shown to demonstrate that the method is as efficient as the IMC scheme, and that it easily allows to incorporate thermodynamical constraints.


title: Splash singularities for a general Oldroyd model with finite Weissenberg number

abstract: In this paper we study a 2D Oldroyd free-boundary model which describes the evolution of a viscoelastic fluid. We prove existence of splash singularities, namely points where the boundary remains smooth but self-intersects. This paper extends the previous results obtained for infinite Weissenberg number to the case of any finite Weissenberg number. The main difficulty of this paper is due to the non linear balance law of the elastic tensor which cannot be reduced, as in the case of infinite Weissenberg, to the transport equations for the deformation gradient. Our strategy in accurate local existence result depending on the Weissenberg number and the combination of Conformal and Lagrangian transformations. The existence of splash singularities is guarantee by a suitable choice of the initial data combined with stability estimates.


title: Deep Learning Based Instance Segmentation in 3D Biomedical Images Using Weak Annotation

abstract: Instance segmentation in 3D images is a fundamental task in biomedical image analysis. While deep learning models often work well for 2D instance segmentation, 3D instance segmentation still faces critical challenges, such as insufficient training data due to various annotation difficulties in 3D biomedical images. Common 3D annotation methods (e.g., full voxel annotation) incur high workloads and costs for labeling enough instances for training deep learning 3D instance segmentation models. In this paper, we propose a new weak annotation approach for training a fast deep learning 3D instance segmentation model without using full voxel mask annotation. Our approach needs only 3D bounding boxes for all instances and full voxel annotation for a small fraction of the instances, and uses a novel two-stage 3D instance segmentation model utilizing these two kinds of annotation, respectively. We evaluate our approach on several biomedical image datasets, and the experimental results show that (1) with full annotated boxes and a small amount of masks, our approach can achieve similar performance as the best known methods using full annotation, and (2) with similar annotation time, our approach outperforms the best known methods that use full annotation.


title: Charge transport properties of the Majorana zero mode induced noncollinear spin selective Andreev reflection

abstract: We study the charge transport of the spin-selective Andreev reflection(SSAR) effect between a spin polarized scanning tunneling microscope(STM) tip and a Majorana zero mode(MZM). Considering both the MZM and the excited states, we calculate the conductance and the shot noise power of the noncollinear SSAR using scattering theory. We find the influence of first excited states cannot be avoided when the spin polarization direction of the STM tip and the MZM are not collinear. In this case, the first excited states give rise to inside peaks and change the conductance peak value at zero energy. Moreover, we numerically calculate the shot noise power and the Fano factor of the SSAR effect. Our calculation shows that the shot noise power and the Fano factor are related to the angle between the spin polarization direction of the STM tip and that of the MZM. These transport properties of the SSAR effect provide additional characteristics to detect the MZM via SSAR.


title: Comment on Demystifying the constancy of the Ermakov-Lewis invariant for a time-dependent oscillator

abstract: We show that a simple modification of the Lagrangian proposed by Padmanabhan in the paper [Mod. Phys. Lett. A 33, 1830005 (2018), arXiv:1712.07328] leads to the most general dynamical invariant in [Ray and Reid, Phys. Lett. A 71, 317 (1979)]


title: Interfacial Structure of SrZr$_{x}$Ti$_{1-x}$O$_3$ films on Ge

abstract: The interfacial structure of SrZr$_{x}$Ti$_{1-x}$O$_3$ films grown on semiconducting Ge substrates are investigated by synchrotron X-ray diffraction and first-principles density functional theory. By systematically tuning Zr content x, the effects of bonding at the interface and epitaxial strain on the physical structure of the film can be distinguished. The interfacial perovskite layers are found to be polarized as a result of cation-anion ionic displacements perpendicular to the perovskite/semiconductor interface. We find a correlation between the observed buckling and valence band offsets at the SrZr$_{x}$Ti$_{1-x}$O$_3$/Ge interface. The theoretical valence band offsets for the polar structures are in agreement with reported X-ray photoelectron spectroscopy measurements. These results have important implications for the integration of functional oxide materials with established semiconductor based technologies.


title: Two-dimensional hydrogen-like atom in a uniform magnetic field. Large-order perturbation theory

abstract: The well known hypervirial perturbation method (HPM)\ based on hypervirial relations and the Hellmann-Feynman theorem is suitable for the calculation of perturbation corrections of large order for the two-dimensional hydrogen-like atom in a uniform magnetic field. We show analytical results in terms of the quantum numbers and large order corrections for particular states. This approach appears to be simpler and more efficient than the recently proposed one based on Green functions.


title: Bloch surface waves at the telecommunication wavelength with Lithium Niobate as top layer for integrated optics

abstract: Lithium niobate (LN) based devices are widely used in integrated and nonlinear optics. This material is robust and resistive to high temperatures, which makes the LN-based devices table, but challenging to fabricate. In this work we report on the design, manufacturing and characterization of engineered dielectric media with thin film lithium niobate (TFLN) on top for the coupling and propagation of electromagnetic surface waves at the telecommunication wavelengths. The designed one-dimentional photonic crystal (1DPhC) sustains Bloch surface waves at the multilayer/air interface at 1550 nm wavelength with a propagation detected over a distance of 3 mm. The working wavelength and improved BSW propagation parameters open the way for exploration of nonlinear properties of BSW based devices. It is also expected that these novel devices will modify BSW propagation and coupling by external thermal/electrical stimuli due to the improved quality of the TFLN top layer of 1DPhC.


title: Low-energy excitations and non-BCS superconductivity in Nb$_x$-Bi$_2$Se$_3$

abstract: When certain elemental metals like Cu, Sr and Nb are intercalated between the layers of Bi$_2$Se$_3$, a topological insulator, the intercalated systems superconduct with critical temperatures around 3 K. Naturally, in all these cases, the possibility of topological superconductivity was suggested and explored. However, in cases of Cu and Sr intercalated systems, the low-temperature scanning tunneling microscopy (STM) experiments revealed fully formed gaps where no signature of low-energy states, a requisite for topological superconductivity, was found. Here, through STM spectroscopy down to 400 mK we show that in Nb$_x$-Bi$_2$Se$_3$ ($x$ = 0.25), the spectra deviate from a BCS-like behavior and the spectral weight at low-bias is large. Our observations are consistent with the idea that the order parameter of Nb$_x$-Bi$_2$Se$_3$ is nodal. Therefore, our results conclude that compared to other members of the family, Nb$_x$-Bi$_2$Se$_3$ has a stronger possibility of being a topological superconductor.


title: Noncommutative Landau problem and shifted energy levels

abstract: In this paper we discuss a Landau levels problem within the framework of noncommutative configuration space and phase space. We show that the associated energy levels are being shifted in terms of the noncommutative parameter and can be directly obtained from its associated generalized Landau Hamiltonian. A straight forward way of constructing the shifted energy levels is just via a proper defined Bopp shift applied onto the spatial and momenta operators that are used to describe the Landau levels in symmetric gauge potential.


title: Solar chemical composition in the hot gas of cool-core ellipticals, groups, and clusters of galaxies

abstract: The hot intracluster medium (ICM) pervading galaxy clusters and groups is rich in metals, which were synthesised by billions of supernovae and have accumulated in cluster gravitational wells for several Gyrs. Since the products of both Type Ia and core-collapse supernovae - expected to explode over different time scales - are found in the ICM, constraining accurately the chemical composition these hot atmospheres can provide invaluable information on the history of the enrichment of large-scale structures. Recently, Hitomi observations reported solar abundance ratios in the core of the Perseus cluster, in tension with previous XMM-Newton measurements obtained for 44 cool-core clusters, groups, and massive ellipticals (the CHEERS sample). In this work, we revisit the CHEERS results by using an updated version of the spectral code used to fit the data (SPEXACT v3), the same as was used to obtain the Hitomi measurements. Despite limitations in the spectral resolution, the average Cr/Fe and Ni/Fe ratios are now found to be remarkably consistent with unity and in excellent agreement with the Hitomi results. Our updated measurements suggest that the solar composition of the ICM of Perseus is a very common feature in nearby cool-core systems.


title: Adversarial Reprogramming of Neural Networks

abstract: Deep neural networks are susceptible to \emph{adversarial} attacks. In computer vision, well-crafted perturbations to images can cause neural networks to make mistakes such as confusing a cat with a computer. Previous adversarial attacks have been designed to degrade performance of models or cause machine learning models to produce specific outputs chosen ahead of time by the attacker. We introduce attacks that instead {\em reprogram} the target model to perform a task chosen by the attacker---without the attacker needing to specify or compute the desired output for each test-time input. This attack finds a single adversarial perturbation, that can be added to all test-time inputs to a machine learning model in order to cause the model to perform a task chosen by the adversary---even if the model was not trained to do this task. These perturbations can thus be considered a program for the new task. We demonstrate adversarial reprogramming on six ImageNet classification models, repurposing these models to perform a counting task, as well as classification tasks: classification of MNIST and CIFAR-10 examples presented as inputs to the ImageNet model.


title: Bispectrum as Baryon Acoustic Oscillation Interferometer

abstract: The galaxy bispectrum, measuring excess clustering of galaxy triplets, offers a probe of dark energy via baryon acoustic oscillations (BAOs). However up to now it has been severely underused due to the combinatorically explosive number of triangles. Here we exploit interference in the bispectrum to identify triangles that amplify BAOs. This approach reduces the computational cost of estimating covariance matrices, offers an improvement in BAO constraints equivalent to lengthening BOSS by 30%, and simplifies adding bispectrum BAO information to future large-scale redshift survey analyses.


title: Augmented Hilbert series of numerical semigroups

abstract: A numerical semigroup $S$ is a subset of the non-negative integers containing $0$ that is closed under addition. The Hilbert series of $S$ (a formal power series equal to the sum of terms $t^n$ over all $n \in S$) can be expressed as a rational function in $t$ whose numerator is characterized in terms of the topology of a simplicial complex determined by membership in $S$. In this paper, we obtain analogous rational expressions for the related power series whose coefficient of $t^n$ equals $f(n)$ for one of several semigroup-theoretic invariants $f:S \to \mathbb R$ known to be eventually quasipolynomial.


title: Statistical Description of Transport in Multimode Fibers with Mode-Dependent Loss

abstract: We analyze coherent wave transport in a new physical setting associated with multimode wave systems where reflection is completely suppressed and mode-dependent losses together with mode-mixing are dictating the wave propagation. An additional physical constraint is the fact that in realistic circumstances the access to the scattering (or transmission) matrix is incomplete. We have addressed all these challenges by providing a statistical description of wave transport which fuses together a free probability theory approach with a Filtered Random Matrix ensemble. Our theoretical predictions have been tested successfully against experimental data of light transport in multimode fibers.


title: Syntax Error Recovery in Parsing Expression Grammars

abstract: Parsing Expression Grammars (PEGs) are a formalism used to describe top-down parsers with backtracking. As PEGs do not provide a good error recovery mechanism, PEG-based parsers usually do not recover from syntax errors in the input, or recover from syntax errors using ad-hoc, implementation-specific features. The lack of proper error recovery makes PEG parsers unsuitable for using with Integrated Development Environments (IDEs), which need to build syntactic trees even for incomplete, syntactically invalid programs. We propose a conservative extension, based on PEGs with labeled failures, that adds a syntax error recovery mechanism for PEGs. This extension associates recovery expressions to labels, where a label now not only reports a syntax error but also uses this recovery expression to reach a synchronization point in the input and resume parsing. We give an operational semantics of PEGs with this recovery mechanism, and use an implementation based on such semantics to build a robust parser for the Lua language. We evaluate the effectiveness of this parser, alone and in comparison with a Lua parser with automatic error recovery generated by ANTLR, a popular parser generator.


title: Knots and solenoids that cannot be attractors of self-homeomorphisms of $\mathbb{R}^3$

abstract: As a first step to understand how complicated attractors for dynamical systems can be, one may consider the following realizability problem: given a continuum $K \subseteq \mathbb{R}^3$, decide when $K$ can be realized as an attractor for a homeomorphism of $\mathbb{R}^3$. In this paper we introduce toroidal sets as those continua $K \subseteq \mathbb{R}^3$ that have a neighbourhood basis comprised of solid tori and, generalizing the classical notion of genus of a knot, give a natural definition of the genus of toroidal sets and study some of its properties. Using these tools we exhibit knots and solenoids for which the answer to the realizability problem stated above is negative.


title: Confronting SUSY GUT with Dark Matter, Sparticle Spectroscopy and Muon $(g-2)$

abstract: We explore the implications of LHC and cold dark matter searches for supersymmetric particle mass spectra in two different grand unified models with left-right symmetry, $SO(10)$ and $ SU(4)_c \times SU(2)_L \times SU(2)_R$ (4-2-2). We identify characteristic differences between the two scenarios, which imply distinct correlations between experimental measurements and the particular structure of the GUT group. The gauge structure of 4-2-2 enhances significantly the allowed parameter space as compared to $SO(10)$, giving rise to a variety of coannihilation scenarios compatible with the LHC data, LSP dark matter and the ongoing muon g-2 experiment.


title: Evaluation of adaptive treatment strategies in an observational study where time-varying covariates are not monitored systematically

abstract: In studies based on electronic health records (EHR), the frequency of covariate monitoring can vary by covariate type, across patients, and over time. This can lead to major challenges: first, the difference in monitoring protocols may invalidate the extrapolation of study results obtained in one population to the other, and second, monitoring can act as a time-varying confounder of the causal effect of a time-varying treatment on the outcomes of interest. This paper demonstrates how to account for non-systematic covariate monitoring when evaluating dynamic treatment interventions, and how to evaluate joint dynamic treatment-censoring and static monitoring interventions, in a real world, EHR-based, comparative effectiveness research (CER) study of patients with type II diabetes mellitus. First, we show that the effects of dynamic treatment-censoring regimes can be identified by including indicators of monitoring events in the adjustment set. Second, we demonstrate the poor performance of the standard inverse probability weighting (IPW) estimator of the effects of joint treatment-censoring-monitoring interventions, due to a large decrease in data support resulting in a large increase in standard errors and concerns over finite-sample bias from near-violations of the positivity assumption for the monitoring process. Finally, we detail an alternate IPW estimator of the effects of these interventions using the No Direct Effect assumption. We demonstrate that this estimator can result in improved efficiency but at the cost of increased bias concerns over structural near-violations of the positivity assumption for the treatment process. To conclude, this paper develops and illustrates new tools that researchers can exploit to appropriately account for non-systematic covariate monitoring in CER, and to ask new causal questions about the joint effects of treatment and monitoring interventions.


title: Toward modern educational IT-ecosystems: from learning management systems to digital platforms

abstract: The development of a learning management system (LMS) as a key service seems to be very effective for creation of educational digital platforms. Such platforms for both higher education institutions and various companies can provide the opportunities for networked forms of educational communication, improve the quality of the perception of innovative technologies and support tools for progress of talented youth as well as knowledge transfer. An example of such LMS is presented. The paper focuses on the demand for further development of learning management systems, their integration with modern digital platforms and potential exploitation as key services of such platforms in the context of the current educational trends of Industry 4.0 and the global trend towards a transition to a digital economy. The implementation of artificial intelligence technologies into the educational process is mentioned as an innovative way to form IT-ecosystems of modern education.


title: The Harish-Chandra integral: An introduction with examples

abstract: This expository paper introduces the theory of Harish-Chandra integrals, a family of special functions that express the integral of an exponential function over the adjoint orbits of a compact Lie group. Originally studied in the context of harmonic analysis on Lie algebras, Harish-Chandra integrals now have diverse applications in many areas of mathematics and physics. We review a number of these applications, present several different proofs of Harish-Chandra's celebrated exact formula for the integrals, and give detailed derivations of the specific integral formulae for all compact classical groups. These notes are intended for mathematicians and physicists who are familiar with the basics of Lie groups and Lie algebras but who may not be specialists in representation theory or harmonic analysis.


title: Proton Form Factor Ratio, $\mu_p G_E^p/G_M^p$ from Double Spin Asymmetry

abstract: The ratio of the electric and magnetic form factor of the proton, $\mu_p G_E^p/G_M^p$, has been measured for elastic electron-proton scattering with polarized beam and target up to four-momentum transfer squared, $Q^2=5.66$ (GeV/c)$^2$ using the double spin asymmetry for target spin orientation aligned nearly perpendicular to the beam momentum direction. This measurement of $\mu_p G_E^p/G_M^p$ agrees with the $Q^2$ dependence of previous recoil polarization data and reconfirms the discrepancy at high $Q^2$ between the Rosenbluth and the polarization-transfer method with a different measurement technique and systematic uncertainties uncorrelated to those of the recoil-polarization measurements. The form factor ratio at $Q^2$=2.06 (GeV/c)$^2$ has been measured as $\mu_p G_E^p/G_M^p = 0.720 \pm 0.176_{stat} \pm 0.039_{sys}$, which is in agreement with an earlier measurement with the polarized target technique at similar kinematics. The form factor ratio at $Q^2$=5.66 (GeV/c)$^2$ has been determined as $\mu_p G_E^p/G_M^p=0.244\pm0.353_{stat}\pm0.013_{sys}$, which represents the highest $Q^2$ reach with the double spin asymmetry with polarized target to date.


title: UNRUH EFFECT: Introductory Notes to Quantum Effects for Accelerated Observers

abstract: These pedagogical notes are dedicated to a derivation of the Unruh effect. There is special emphasis on the transparency of the arguments and the exhibition of detailed calculations. We assume the reader has a basic knowledge of quantum mechanics and relativity. The notes can hopefully be used in a self-contained way by advanced undergraduate or beginning graduate students.


title: Further investigation of changes in cometary rotation

abstract: Samarasinha & Mueller (2013) related changes of cometary rotation to other physical parameters for four Jupiter family comets defining a parameter $X$, which is approximately constant within a factor of two irrespective of the active fraction of a comet. Two additional comets are added to this sample in this paper and the claim of a nearly constant parameter $X$ for these six comets is confirmed, albeit with a larger scatter. Taking the geometric mean of $X$ for all the comets above excluding 2P/Encke (as $X$ for each comet was determined with respect to that of 2P/Encke), the expected changes in the rotation periods for a sample of 24 periodic comets are derived. We identify comets from this sample that are most likely to show observationally detectable changes in their rotation periods. Using this sample and including the six comets used to determine $X$, we find a correlation between the parameter $\zeta$ (i.e. the total water production per unit surface area per orbit approximated by that inside of 4 au) and the perihelion distance $q$; specifically we derive $\zeta$ $\propto$ $q^{-0.8}$ and provide a theoretical basis for this in Appendix A. This relationship between $\zeta$ and $q$ enables ready comparisons of activity due to insolation between comets. Additionally, a relationship between the nuclear radius $R$ and the rotation period $P$ is found. Specifically, we find that on average smaller nuclei have smaller rotation periods compared to the rotation periods of larger nuclei. This is consistent with expectations for rotational evolution and spin-up of comet nuclei, providing strong observational evidence for sublimation-driven rotational changes in comets.


title: A1-invariants in Galois cohomology and a claim of Morel

abstract: We establish a variant of the splitting principle of Garibaldi-Merkurjev-Serre for invariants taking values in a strictly homotopy invariant sheaf. As an application, we prove the folklore result of Morel that pi_0 of the motivic localization of the group completion of the stack of finite \'etale schemes is given by the sheaf of unramified Grothendieck-Witt groups.


title: Tracking black hole kicks from gravitational wave observations

abstract: Coalescing binary black holes emit anisotropic gravitational radiation. This causes a net emission of linear momentum that produces a gradual acceleration of the source. As a result, the final remnant black hole acquires a characteristic velocity known as recoil velocity or gravitational kick. The symmetries of gravitational wave emission are reflected in the interactions of the gravitational wave modes emitted by the binary. In this work we make use of the rich information encoded in the higher-order modes of the gravitational wave emission to infer the component of the kick along the line-of-sight (or \textit{radial kick}). We do this by performing parameter inference on simulated signals given by numerical relativity waveforms for non-spinning binaries using numerical relativity templates of aligned-spin (non-precessing) binary black holes. We find that for suitable sources, namely those with mass ratio $q\geq 2$ and total mass $M \sim 100M_\odot$, and for modest radial kicks of $120km/s$, the $90\%$ credible intervals of our posterior probability distributions can exclude a zero kick at a signal-to-noise ratio of $15$; using a single Advanced LIGO detector working at its early sensitivity. The measurement of a non-zero radial kick component would provide the first observational signature of net transport of linear momentum by gravitational waves away from their source.


title: Late afterglows of GW/GRB 170817A

abstract: The gamma-ray burst that followed the first detection of gravitational waves from the merger of a Binary Neutron Stars and its low energy counterparts were in many respects unusual and challenge our understanding of mechanisms involved in their production. In a previous work we used a phenomenological formulation of relativistic shocks and synchrotron emission to analyse the prompt gamma-ray emission of GW/GRB 170817A. Here we use the same model to analyse late afterglows of this event. The main goal is to see whether synchrotron emission alone can explain late afterglows. We find that collision between a mildly relativistic outflow from the merger with a Lorentz factor of $\sim 1.2-3$ and the ISM/circumburst material can explain observations, if the synchrotron self-absorption of radio emission and extinction of optical/IR photons are taken into account. Alternatively, an additional source of X-ray is necessary to explain the data. These conclusions are independent of the model used here and can be deduced directly from data. We also show that at the time of its encounter with circumburst material the outflow could have been still mildly magnetized. The origin for optical extinction could be a dust rich old faint star cluster surrounding the BNS, which had also helped its formation and merger. Such environment evades present observational constraints and is consistent with our conclusions about properties and evolution of the progenitor neutron stars obtained from the prompt gamma-ray. If the synchrotron emission was produced by collision between density shells, the extinction might have occurred inside the outflow itself rather than externally. A plausible source of additional X-ray is decay of medium and heavy isotopes produced by the kilonova, including r-processes, and recombination of cooled electrons. Contribution of these processes should be quantified in future works.


title: Fluctuation Dynamo in a Collisionless, Weakly Magnetized Plasma

abstract: Results from a numerical study of fluctuation dynamo in a collisionless, weakly magnetized plasma are presented. The key difference between this dynamo and its magnetohydrodynamic (MHD) counterpart is the adiabatic production of magnetic-field-aligned pressure anisotropy by the amplification of a weak seed field. This in turn drives kinetic instabilities on the ion-Larmor scale---namely, firehose and mirror---which sever the adiabatic link between the thermal and magnetic pressures, thereby allowing the dynamo to proceed. After an initial phase of rapid growth driven by these instabilities, the magnetic energy grows exponentially and exhibits a $k^{3/2}$ spectrum that peaks near the resistive scale, similar to the large-magnetic-Prandtl-number ($\mathrm{Pm}\gg{1}$) MHD dynamo. The magnetic field self-organizes into a folded-sheet topology, with direction reversals at the resistive scale and field lines curved at the parallel scale of the flow. The effective $\mathrm{Pm}$ is determined by whether the ion-Larmor scale is above or below the field-reversing scale: in the former case, particles undergo Bohm-like diffusion; in the latter case, particles scatter primarily off firehose fluctuations residing at the ends of the magnetic folds, and the viscosity becomes anisotropic. The magnetic field ultimately saturates at dynamical strengths, with its spectral peak migrating towards larger scales. This feature, along with an anti-correlation of magnetic-field strength and field-line curvature and a gradual thinning of magnetic sheets into ribbons, resemble the saturated state of the large-$\mathrm{Pm}$ dynamo, the primary differences manifesting in firehose/mirror-unstable regions. These results have implications for magnetic-field growth in the weakly collisional intracluster medium of galaxy clusters.


title: Slow Cooling in Low Metallicity Clouds: An Origin of Globular Cluster Bimodality?

abstract: We explore the relative role of small-scale fragmentation and global collapse in low-metallicity clouds, pointing out that in such clouds the cooling time may be longer than the dynamical time, allowing the cloud to collapse globally before it can fragment. This, we suggest, may help to explain the formation of the low-metallicity globular cluster population, since such dense stellar systems need a large amount of gas to be collected in a small region (without significant feedback during the collapse). To explore this further, we carry out numerical simulations of low-metallicity Bonner-Ebert stable gas clouds, demonstrating that there exists a critical metallicity (between 0.001 and 0.01 $Z_\odot$) below which the cloud collapses globally without fragmentation. We also run simulations including a background radiative heating source, showing that this can also produce clouds that do not fragment, and that the critical metallicity -- which can exceed the no-radiation case -- increases with the heating rate.


title: Neutron star bulk viscosity, "spin-flip" and GW emission of newly born magnetars

abstract: The viscosity-driven "spin-flip" instability in newborn magnetars with interior toroidal magnetic fields is re-examined. We calculate the bulk viscosity coefficient ($\zeta$) of cold, $npe \mu$ matter in neutron stars (NS), for selected values of the nuclear symmetry energy and in the regime where $\beta$-equilibration is slower than characteristic oscillation periods. We show that: i) $\zeta$ is larger than previously assumed and the instability timescale correspondingly shorter; ii) for a magnetically-induced ellipticity $\epsilon_B \lesssim 4 \times 10^{-3}$, typically expected in newborn magnetars, spin-flip occurs for initial spin periods $\lesssim 2-3$ ms, with some dependence on the NS equation of state (EoS). We then calculate the detectability of GW signals emitted by newborn magnetars subject to "spin-flip", by accounting also for the reduction in range resulting from realistic signal searches. For an optimal range of $\epsilon_B \sim (1-5) \times 10^{-3}$, and birth spin period $\lesssim 2$ ms, we estimate an horizon of $\gtrsim 4$ Mpc, and $\gtrsim 30$ Mpc, for Advanced and third generation interferometers at design sensitivity, respectively. A supernova (or a kilonova) is expected as the electromagnetic counterpart of such GW events. Outside of the optimal range for GW emission, EM torques are more efficient in extracting the NS spin energy, which may power even brighter EM transients.


title: Spherical harmonics method for computing the image stress due to a spherical void

abstract: We develop an efficient numerical method for calculating the image stress field induced by spherical voids in materials. The method is applied to dislocation-void interactions in dislocation dynamics simulations. We obtain a complete set of basis functions for the solution of the image stress problem, as well as their corresponding displacement, stress, and traction fields in terms of the linear combination of spherical harmonics. Using the fast transformation between the real and spherical-harmonics spaces provided by the SHTOOLS package, the method is more efficient than other image stress solvers such as the finite-element method. This method can readily be extended for solving elasticity problems involving inclusions and inhomogeneities, contact between spheres, as well as other differential equations with spherical boundaries beyond elasticity.


title: Secrecy Beamforming for SWIPT MISO Heterogeneous Cellular Networks

abstract: In this paper, we consider the secure transmission design for a multiple-input single-output Femtocell overlaid with a Macrocell in co-channel deployment. The Femtocell base station sends confidential messages to information receiving Femtocell users (FUs) and energy signals to energy receiving (ER) FUs while limiting the interference to Macrocell users (MUs). The ER FUs have the potential to wiretap the confidential messages. By taking fairness into account, we propose a sum logarithmic secrecy rate maximization beamforming design problem under the interference constraints for MUs and energy harvesting (EH) constraints for ER FUs. The formulated design problem is nontrivial to solve due to the nonconvexity which lies in the objective and the constraints. To tackle the design problem, a semidefinite relaxation and successive convex approximation based algorithm is proposed. Simulation results demonstrate the effectiveness of the proposed beamforming design.


title: The symbiotic contact process: phase transitions, hysteresis cycles, and bistability

abstract: We performed Monte Carlo simulations of the symbiotic contact process on different spatial dimensions ($d$). On the complete and random graphs (infinite dimension), we observe hysteresis cycles and bistable regions, what is consistent with the discontinuous absorbing-state phase transition predicted by mean-field theory. By contrast, on a regular square lattice, we find no signs of bistability or hysteretic behavior. This result suggests that the transition in two dimensions is rather continuous. Based on our numerical observations, we conjecture that the nature of the transition changes at the upper critical dimension ($d_c$), from continuous ($d<d_c$) to discontinuous ($d>d_c$).


title: Domestic Corpuscular Inflaton

abstract: The aim of this paper is to provide a more precise description of the paradigm of corpuscular slow-roll inflation, which was previously introduced by Casadio et al. in [1]. Specifically, we start by expanding the Starobinsky theory on a curved background and then infer the number and nature of the propagating degrees of freedom, both in the true inflationary phase and in a quasi-de Sitter approximation. We correctly find that the particle spectrum contains a transverse trace-free mode and a scalar one. The scalar mode displays a tachyonic nature during the slow-roll phase, due to the instability of the system, whereas it acquires the appropriate oscillatory behavior as the background approaches a critical value of the curvature. These results confirm the fact that the Einstein-Hilbert term acts as a perturbation to the quadratic one, and is responsible for driving the early universe out of the inflationary phase, thus realising the inflaton field in terms of pure (corpuscular) gravity.


title: 3D Normal Coordinate Systems for Cortical Areas

abstract: A surface-based diffeomorphic algorithm to generate 3D coordinate grids in the cortical ribbon is described. In the grid, normal coordinate lines are generated by the diffeomorphic evolution from the grey/white (inner) surface to the grey/csf (outer) surface. Specifically, the cortical ribbon is described by two triangulated surfaces with open boundaries. Conceptually, the inner surface sits on top of the white matter structure and the outer on top of the gray matter. It is assumed that the cortical ribbon consists of cortical columns which are orthogonal to the white matter surface. This might be viewed as a consequence of the development of the columns in the embryo. It is also assumed that the columns are orthogonal to the outer surface so that the resultant vector field is orthogonal to the evolving surface. Then the distance of the normal lines from the vector field such that the inner surface evolves diffeomorphically towards the outer one can be construed as a measure of thickness. Applications are described for the auditory cortices in human adults and cats with normal hearing or hearing loss. The approach offers great potential for cortical morphometry.


title: GenerationMania: Learning to Semantically Choreograph

abstract: Beatmania is a rhythm action game where players must reproduce some of the sounds of a song by pressing specific controller buttons at the correct time. In this paper we investigate the use of deep neural networks to automatically create game stages - called charts - for arbitrary pieces of music. Our technique uses a multi-layer feed-forward network trained on sound sequence summary statistics to predict which sounds in the music are to be played by the player and which will play automatically. We use another neural network along with rules to determine which controls should be mapped to which sounds. We evaluated our system on the ability to reconstruct charts in a held-out test set, achieving an $F_1$-score that significantly beats LSTM baselines.


title: Real-Time Optimal Power Flow under Wind Energy Penetration-Part I: Approach

abstract: Real-time optimal power flow (RT-OPF) under wind energy penetration is highly desired but extremely difficult to realize. This is basically due to the conflict between the fast changes in wind power generation and the slow response from the optimization computation. This paper (Part I) presents a prediction-updating approach to address this challenge. We consider essential scenarios around forecasted data of wind power that would probably happen during the computation time required for solving a large-scale complex optimal power flow problem. Parallel computing is used to solve the individual OPF problems corresponding to these scenarios. This provides for the forecasted time horizon probable reference operations in the form of a lookup-table. One of these operations will be selected based on the actual wind power and realized to the grid for the current time interval, thus leading to a RT-OPF framework. The proposed approach is implemented in Part II of this paper using a 41-bus medium-voltage distribution network as a case study.


title: Real-Time Optimal Power Flow under Wind Energy Penetration-Part II: Implementation

abstract: In this paper (Part II) we implement the prediction-updating approach developed in Part I to address fast changes in wind power generation when solving a complex real-time optimal power flow (RT-OPF) problem. The approach considers essential scenarios around forecasted wind power values in a moving prediction horizon (120 seconds). The individual optimal power flow problems corresponding to these scenarios are solved in parallel using a multi-processor server. Then the operation strategy is updated in a short sampling time (every 20 seconds) considering real wind power values. The RT-OPF problem is formulated considering both technical and economic aspects simultaneously. The RT-OPF is implemented on a 41-bus medium-voltage distribution network with two wind stations. The results show the benefits of the proposed approach and highlight further challenges of RT-OPF.


title: Revisiting the science case for near-UV spectroscopy with the VLT

abstract: In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at far-blue visible wavelengths for the foreseeable future. High-efficiency (>20%) observations of the ground UV (300-400 nm) at medium resolving power (R~20,000) are required to address a number of exciting topics in stellar astrophysics, while also providing new insights in extragalactic science. Anticipating strong demand to better exploit this diagnostic-rich wavelength region, we revisit the science case and instrument requirements previously assembled for the CUBES concept for the Very Large Telescope.


title: Notes on ten-dimensional localized black holes and deconfined states in two-dimensional SYM

abstract: We numerically construct static localized black holes in ten spacetime dimensions with one compact periodic dimension. In particular, we investigate the critical regime in which the poles of the localized black hole are about to merge. When approaching the critical region, the behavior of physical quantities is described by a single real valued exponent giving rise to a logarithmic scaling of the thermodynamic quantities, in agreement with the theoretical prediction derived from the double-cone metric. As a peculiarity, the localized black hole solution in ten dimensions can be related to the spatially deconfined phase of two dimensional N=(8,8) super Yang-Mills theory (SYM) on a spatial circle. We use the localized black hole solutions to determine the SYM phase diagram. In particular, we compute the location of the first order phase confinement/deconfinement transition and the related latent heat to unprecedented accuracy.


title: Combining Radon transform and Electrical Capacitance Tomography for a $2d+1$ imaging device

abstract: This paper describes a coplanar non invasive non destructive capacitive imaging device. We first introduce a mathematical model for its output, and discuss some of its theoretical capabilities. We show that the data obtained from this device can be interpreted as a weighted Radon transform of the electrical permittivity of the measured object near its surface. Image reconstructions from experimental data provide good surface resolution as well as short depth imaging, making the apparatus a $2d+1$ imager. The quality of the images leads us to expect that excellent results can be delivered by \emph{ad-hoc} optimized inversion formulas. There are also interesting, yet unexplored, theoretical questions on imaging that this sensor will allow to test.


title: Remarks on Dark Matter Constituents with Many Solar Masses

abstract: Dark matter constituents of many solar masses will accrete normal matter which emits X-rays that can be downgraded to microwaves which may distort the precisely-measured black-body spectrum of the Cosmic Microwave Background. However, it is known from elsewhere that spherical models of accretion vastly overestimate the amount accreted and consequently the emitted X-rays. Therefore, exclusion plots based on spherical accretion for the allowed fraction of the dark matter versus the MACHO mass give upper limits on intermediate-mass MACHOs which are too severe, sometimes by orders of magnitude.


title: ExoMol line lists XXV: A hot line list for silicon sulphide, SiS

abstract: SiS has long been observed in the circumstellar medium of the carbon-rich star IRC+10216 CW Leo. Comprehensive and accurate rotation-vibrational line lists and partition functions are computed for 12 isotopologues of silicon sulphide ($^{28}$Si$^{32}$S, $^{28}$Si$^{34}$S, $^{29}$Si$^{32}$S, $^{28}$Si$^{33}$S, $^{30}$Si$^{32}$S, $^{29}$Si$^{34}$S, $^{30}$Si$^{34}$S, $^{28}$Si$^{36}$S, $^{29}$Si$^{33}$S, $^{29}$Si$^{36}$S, $^{30}$Si$^{33}$S, $^{30}$Si$^{36}$S) in its ground (X $^1\Sigma^+$) electronic state. The calculations employ an existing spectroscopically accurate potential energy curve (PEC) derived from experimental measurements and a newly-computed ab initio dipole moment curve (DMC). The $^{28}$Si$^{32}$S line list includes 10104 states and 91715 transitions. These line lists are available from the ExoMol website www.exomol.com and the CDS database.


title: Scoring Alternative Forecast Distributions: Completing the Kullback Distance Complex

abstract: We develop two surprising new results regarding the use of proper scoring rules for evaluating the predictive quality of two alternative sequential forecast distributions. Both of the proponents prefer to be awarded a score derived from the other's distribution rather than a score awarded on the basis of their own. A Pareto optimal exchange of their scoring outcomes provides the basis for a comparison of forecast quality that is preferred by both forecasters, and also evades a feature of arbitrariness inherent in using the forecasters' own achieved scores. The well-known Kullback divergence, used as a measure of information, is evaluated via the entropies in the two forecast distributions and the two cross-entropies between them. We show that Kullback's symmetric measure needs to be appended by three component measures if it is to characterise completely the information content of the two asserted probability forecasts. Two of these do not involve entropies at all. The resulting 'Kullback complex' supported by the 4-dimensional measure is isomorphic to an equivalent vector measure generated by the forecasters' expectations of their scores, each for one's own score and for the other's score. We foreshadow the results of a sophisticated application of the Pareto relative scoring procedure for actual sequentional observations, and we propose a standard format for evaluation.


title: Astrophysical Nuclear Reactions: from Hydrogen Burning to Supernovae Explosions

abstract: The work reported in this dissertation will concern the study of kinetics, cycles and quantum mechanic processes of nuclear reactions involved in massive stars together with explosive nucleosynthesis phenomena having a key role in supernovae explosion. Four main fields will be explored: physics of nuclear reaction, hydrogen burning, CNO cycles and explosive supernovae phenomena. The first chapter will provide an overview of the background of the nuclear physics and of the potential approaches to the field as quantum mechanic aspects, cross section and Maxwell-Boltzmann distribution. The next chapter will describe the hydrogen burning processes for stars with low mass (our sun), the CNO and hot CNO cycles involved for higher mass stars and other cycles as sodium and magnesium. The third chapter will review the helium burning process involved in Red Giants as result of the exhausted hydrogen fuel. A short overview of other helium processes will be briefly discussed. The fourth chapter will provide a summary of the accepted model of explosive burning process involved in supernovae explosion. Explosive nucleosynthesis together with heavy elements as silicon burning processes will be also described.


title: Open class determinacy is preserved by forcing

abstract: The principle of open class determinacy is preserved by pre-tame class forcing, and after such forcing, every new class well-order is isomorphic to a ground-model class well-order. Similarly, the principle of elementary transfinite recursion $\text{ETR}_\Gamma$ for a fixed class well-order $\Gamma$ is preserved by pre-tame class forcing. The full principle ETR itself is preserved by countably strategically closed pre-tame class forcing, and after such forcing, every new class well-order is isomorphic to a ground-model class well-order. Meanwhile, it remains open whether ETR is preserved by all forcing, including the forcing merely to add a Cohen real.


title: Space-time duality and quantum temporal imaging

abstract: Using the space-time analogy, we compare the performance of quantum temporal imaging with its classical counterpart. We consider a temporal imaging scheme, based on the sum-frequency generation (SFG) time lens, but our results can be applied to other temporal imaging schemes such as, for instance, four-wave mixing. Extending the theory presented in our previous publications, in this paper we take into account the finite time aperture of the imaging system, characterized by its pupil function. Using the quantum theory, we obtain a unitary transformation of the quantum field from the input to the output of the imaging scheme and identify the contribution of the vacuum fluctuations missing in the classical theory. This contribution plays a key role in the quantum temporal imaging of nonclassical temporal waveforms, characterized by nonclassical fluctuations of the electromagnetic field. As an example, we consider quantum temporal imaging of broadband squeezed light and formulate the criteria for conservation of its squeezing properties at the output of the system.


title: Tuning of MC generator MPI models

abstract: MC models of multiple partonic scattering inevitably introduce many free parameters, either fundamental to the models or from their integration with MC treatments of primary-scattering evolution. This non-perturbative and non-factorisable physics in particular cannot currently be constrained from theoretical principles, and hence parameter optimisation against experimental data is required. This process is commonly referred to as MC tuning. We summarise the principles, problems and history of MC tuning, and the still-evolving modern approach to both model optimisation and estimation of modelling uncertainties.


title: Cross-Discourse and Multilingual Exploration of Textual Corpora with the DualNeighbors Algorithm

abstract: Word choice is dependent on the cultural context of writers and their subjects. Different words are used to describe similar actions, objects, and features based on factors such as class, race, gender, geography and political affinity. Exploratory techniques based on locating and counting words may, therefore, lead to conclusions that reinforce culturally inflected boundaries. We offer a new method, the DualNeighbors algorithm, for linking thematically similar documents both within and across discursive and linguistic barriers to reveal cross-cultural connections. Qualitative and quantitative evaluations of this technique are shown as applied to two cultural datasets of interest to researchers across the humanities and social sciences. An open-source implementation of the DualNeighbors algorithm is provided to assist in its application.


title: Modeling the growth of organisms validates a general relation between metabolic costs and natural selection

abstract: Metabolism and evolution are closely connected: if a mutation incurs extra energetic costs for an organism, there is a baseline selective disadvantage that may or may not be compensated for by other adaptive effects. A long-standing, but to date unproven, hypothesis is that this disadvantage is equal to the fractional cost relative to the total resting metabolic expenditure. This hypothesis has found a recent resurgence as a powerful tool for quantitatively understanding the strength of selection among different classes of organisms. Our work explores the validity of the hypothesis from first principles through a generalized metabolic growth model, versions of which have been successful in describing organismal growth from single cells to higher animals. We build a mathematical framework to calculate how perturbations in maintenance and synthesis costs translate into contributions to the selection coefficient. This allows us to show that the hypothesis is an approximation to the actual baseline selection coefficient. Moreover we can derive the correct prefactor in its functional form, as well as analytical bounds on its accuracy for any given realization of the model. We illustrate our general framework using a special case of the growth model, which we show provides a quantitative description of metabolic expenditures in data collected from a wide array of unicellular organisms (both prokaryotes and eukaryotes). In all these cases we demonstrate that the hypothesis is an excellent approximation, allowing estimates of baseline selection coefficients to within 15% of their actual values. Even in a broader biological parameter range, covering growth data from multicellular organisms, the hypothesis continues to work well, always within an order of magnitude of the correct result. Our work thus justifies its use as a versatile tool, setting the stage for its wider deployment.


title: Intra-pixel response characterization of a HgCdTe near infrared detector with a pronounced crosshatch pattern

abstract: The "crosshatch" pattern is a recurring "feature" of HgCdTe arrays, specifically the Teledyne HAWAII family of near infrared detectors. It is a fixed pattern of high frequency QE variations along 3 directions generally thought to be related to the crystal structure of HgCdTe. The pattern is evident in detectors used in Hubble WFC3, WISE, JWST, and in candidate detectors for Euclid and WFIRST. Using undersampled point sources projected onto a HAWAII-2RG detector, we show that the pattern induces photometric variations that are not removed by a flat-field calibration, thus demonstrating that the QE variations occur on scales smaller than the 18 micron pixels. Our testbed is the Precision Projector Laboratory's astronomical scene generator, which can rapidly characterize the full detector by scanning thousands of undersampled spots. If not properly calibrated, detectors showing strong crosshatch may induce correlated errors in photometry, astrometry, spectroscopy, and shape measurements.


title: A New Angle on L2 Regularization

abstract: Imagine two high-dimensional clusters and a hyperplane separating them. Consider in particular the angle between: the direction joining the two clusters' centroids and the normal to the hyperplane. In linear classification, this angle depends on the level of L2 regularization used. Can you explain why?


title: Neural-net-induced Gaussian process regression for function approximation and PDE solution

abstract: Neural-net-induced Gaussian process (NNGP) regression inherits both the high expressivity of deep neural networks (deep NNs) as well as the uncertainty quantification property of Gaussian processes (GPs). We generalize the current NNGP to first include a larger number of hyperparameters and subsequently train the model by maximum likelihood estimation. Unlike previous works on NNGP that targeted classification, here we apply the generalized NNGP to function approximation and to solving partial differential equations (PDEs). Specifically, we develop an analytical iteration formula to compute the covariance function of GP induced by deep NN with an error-function nonlinearity. We compare the performance of the generalized NNGP for function approximations and PDE solutions with those of GPs and fully-connected NNs. We observe that for smooth functions the generalized NNGP can yield the same order of accuracy with GP, while both NNGP and GP outperform deep NN. For non-smooth functions, the generalized NNGP is superior to GP and comparable or superior to deep NN.


title: A Proof-of-principle for Time-Of-Flight Positron Emission Tomography Imaging

abstract: Time-Of-Flight (TOF) is a noble technique that is used in Positron Emission Tomography (PET) imaging worldwide. The scintillator based imaging system that is being used around the world for TOF-PET is very expensive. Multi-gap Resistive Plate Chambers (MRPCs) are gaseous detectors which are easy to fabricate, inexpensive and have excellent position and timing resolution. They can be used as a suitable alternative to highly expensive scintillators. For the sole purpose of TOF-PET, a pair of 18 cm $\times$ 18 cm, 5 gap, glass-based MRPC modules have been fabricated. Our main aim was to determine the shift in the position of the source (Na-22) with these fabricated MRPCs. In this document, the details of the experimental results will be presented.


title: A hybrid deep learning approach for medical relation extraction

abstract: Mining relationships between treatment(s) and medical problem(s) is vital in the biomedical domain. This helps in various applications, such as decision support system, safety surveillance, and new treatment discovery. We propose a deep learning approach that utilizes both word level and sentence-level representations to extract the relationships between treatment and problem. While deep learning techniques demand a large amount of data for training, we make use of a rule-based system particularly for relationship classes with fewer samples. Our final relations are derived by jointly combining the results from deep learning and rule-based models. Our system achieved a promising performance on the relationship classes of I2b2 2010 relation extraction task.


title: Fully Distributed Cooperative Charging for Plug-in Electric Vehicles in Constrained Power Networks

abstract: Plug-in Electric Vehicles (PEVs) play a pivotal role in transportation electrification. The flexible nature of PEVs' charging demand can be utilized for reducing charging cost as well as optimizing the operating cost of power and transportation networks. Utilizing charging flexibilities of geographically spread PEVs requires design and implementation of efficient optimization algorithms. To this end, we propose a fully distributed algorithm to solve the PEVs' Cooperative Charging with Power constraints (PEV-CCP). Our solution considers the electric power limits that originate from physical characteristics of charging station, such as on-site transformer capacity limit, and allows for containing charging burden of PEVs on the electric distribution network. Our approach is also motivated by the increasing load demand at the distribution level due to additional PEV charging demand. Our proposed approach distributes computation among agents (PEVs) to solve the PEV-CCP problem in a distributed fashion through an iterative interaction between neighboring agents. The structure of each agent's update functions ensures an agreement on a price signal while enforcing individual PEV constraints. In addition to converging towards the globally-optimum solution, our algorithm ensures the feasibility of each PEV's decision at each iteration. We have tested performance of the proposed approach using a fleet of PEVs.


title: CR-GAN: Learning Complete Representations for Multi-view Generation

abstract: Generating multi-view images from a single-view input is an essential yet challenging problem. It has broad applications in vision, graphics, and robotics. Our study indicates that the widely-used generative adversarial network (GAN) may learn "incomplete" representations due to the single-pathway framework: an encoder-decoder network followed by a discriminator network. We propose CR-GAN to address this problem. In addition to the single reconstruction path, we introduce a generation sideway to maintain the completeness of the learned embedding space. The two learning pathways collaborate and compete in a parameter-sharing manner, yielding considerably improved generalization ability to "unseen" dataset. More importantly, the two-pathway framework makes it possible to combine both labeled and unlabeled data for self-supervised learning, which further enriches the embedding space for realistic generations. The experimental results prove that CR-GAN significantly outperforms state-of-the-art methods, especially when generating from "unseen" inputs in wild conditions.


title: Constraints on Upper Cutoffs in the Mass Functions of Young Star Clusters

abstract: We test claims that the power-law mass functions of young star clusters (ages $\lesssim\mbox{few}\times10^8$~yr) have physical upper cutoffs at $M_*\sim10^5~M_{\odot}$. Specifically, we perform maximum-likelihood fits of the Schechter function, $\psi(M)=dN/dM\propto M^{\beta}~\mbox{exp}(-M/M_*)$, to the observed cluster masses in eight well-studied galaxies (LMC, SMC, NGC 4214, NGC 4449, M83, M51, Antennae, and NGC 3256). In most cases, we find that a wide range of cutoff mass is permitted ($10^5~M_\odot \lesssim M_* < \infty$). We find a weak detection at $M_* \sim 10^5~M_\odot$ in one case (M51) and strong evidence against this value in two cases. However, when we include realistic errors in cluster masses in our analysis, the constraints on $M_*$ become weaker and there are no significant detections (even for M51). Our data are generally consistent with much larger cutoffs, at $M_*\sim\mbox{few}\times10^6~M_{\odot}$. This is the predicted cutoff from dynamical models in which old globular clusters and young clusters observed today formed by similar physical processes with similar initial mass functions.


title: Stationary double black hole without naked ring singularity

abstract: Recently double black hole vacuum and electrovacuum metrics attracted attention as exact solutions suitable for visualization of ultra-compact objects beyond the Kerr paradigm. However, many of the proposed systems are plagued with ring curvature singularities. Here we present a new simple solution of this type which is asymptotically Kerr, has zero electric and magnetic charges, but is endowed with magnetic dipole moment and electric quadrupole moment. It is manifestly free of ring singularities, and contains only a mild string-like singularity on the axis corresponding to a distributional energy-momentum tensor.Its main constituents are two extreme co-rotating black holes carrying equal electric andopposite magnetic and NUT charges.

categories: gr-qc astro-ph.HE hep-th math-ph math.MP