title: The healing strategy by prioritizing minimum degree against localized attacks on interdependent spatially embedded networks

abstract: Many real infrastructure networks, such as power grids and communication networks, are not only depend on one another to function, but also embedded in space. A lot of works have been devoted to reveal the vulnerability of interdependent spatially embedded networks considers random failures. However, recently research show that they are susceptible to geographically localized attacks or failures cased by natural disasters or terrorist attacks, which affect all nodes within a given radius. In particular, small localized attacks may lead to catastrophic consequences. As a remedy of the collapse instability, one research introduced a dynamic healing strategy and the possibility of new connectivity link with a probability, called random healing, to bridge two remaining neighbors of a failed node. The random strategy is straightforward. Here, unlike previous strategy, we propose a simple but effective strategy, called healing strategy by prioritizing minimum degree(HPMD), which establishing a new link between two functioning neighbors with minimum degree during the cascading failures. In addition, we consider the distance between two nodes to avoid change the network structures. Afterwards, a comparison is made between HPMD and three healing strategies: random, degree centrality and local centrality to identify which a pair of neighbors have high priority. Simulation experiments are presented for two square lattices with the same size as interdependent spatially embedded networks under localized attacks. Results demonstrated that HPMD strategy has an outstanding performance against localized attacks, even when ratio of healing is vary. Moreover, HPMD is more timely, more applicability and costless. Our method is meaningful in practice as it can greatly enhance the resilience and robustness of interdependent spatially embedded networks against localized attacks.


title: OS Scheduling Algorithms for Improving the Performance of Multithreaded Workloads

abstract: Major chip manufacturers have all introduced multicore microprocessors. Multi-socket systems built from these processors are used for running various server applications. However to the best of our knowledge current commercial operating systems are not optimized for multi-threaded workloads running on such servers. Cache-to-cache transfers and remote memory accesses impact the performance of such workloads. This paper presents a unified approach to optimizing OS scheduling algorithms for both cache-to-cache transfers and remote DRAM accesses that also takes cache affinity into account. By observing the patterns of local and remote cache-to-cache transfers as well as local and remote DRAM accesses for every thread in each scheduling quantum and applying different algorithms, we come up with a new schedule of threads for the next quantum taking cache affinity into account. This new schedule cuts down both remote cache-to-cache transfers and remote DRAM accesses for the next scheduling quantum and improves overall performance. We present two algorithms of varying complexity for optimizing cache-to-cache transfers. One of these is a new algorithm which is relatively simpler and performs better when combined with algorithms that optimize remote DRAM accesses. For optimizing remote DRAM accesses we present two algorithms. Though both algorithms differ in algorithmic complexity we find that for our workloads they perform equally well. We used three different synthetic workloads to evaluate these algorithms. We also performed sensitivity analysis with respect to varying remote cache-to-cache transfer latency and remote DRAM latency. We show that these algorithms can cut down overall latency by up to 16.79% depending on the algorithm used.


title: The correspondence between Voigt and Reuss bounds and the decoupling constraint in a two-grid staggered solution algorithm to coupled flow and deformation in heterogeneous poroelastic media

abstract: We perform a convergence analysis of a two-grid staggered solution algorithm for the Biot system modeling coupled flow and deformation in heterogeneous poroelastic media. The algorithm first solves the flow subproblem on a fine grid using a mixed finite element method (by freezing a certain measure of the mean stress) followed by the poromechanics subproblem on a coarse grid using a conforming Galerkin method. Restriction operators map the fine scale flow solution to the coarse scale poromechanical grid and prolongation operators map the coarse scale poromechanical solution to the fine scale flow grid. The coupling iterations are repeated until convergence and Backward Euler is employed for time marching. The analysis is based on studying the equations satisfied by the difference of iterates to show that the two-grid scheme is a contraction map under certain conditions. Those conditions are used to construct the restriction and prolongation operators as well as arrive at coarse scale elastic properties in terms of the fine scale data. We show that the adjustable parameter in the measure of the mean stress is linked to the Voigt and Reuss bounds frequently encountered in computational homogenization of multiphase composites.


title: Digital holographic particle volume reconstruction using a deep neural network

abstract: This paper proposes a particle volume reconstruction directly from an in-line hologram using a deep neural network. Digital holographic volume reconstruction conventionally uses multiple diffraction calculations to obtain sectional reconstructed images from an in-line hologram, followed by detection of the lateral and axial positions, and the sizes of particles by using focus metrics. However, the axial resolution is limited by the numerical aperture of the optical system, and the processes are time-consuming. The method proposed here can simultaneously detect the lateral and axial positions, and the particle sizes via a deep neural network (DNN). We numerically investigated the performance of the DNN in terms of the errors in the detected positions and sizes. The calculation time is faster than conventional diffracted-based approaches.


title: Sphere theorems for Lagrangian and Legendrian submanifolds

abstract: In this paper, we prove some differentiable sphere theorems and topological sphere theorems for Lagrangian submanifolds in K\"ahler manifold and Legendrian submanifolds in Sasaki space form.


title: Atomic Characterizations of Weak Martingale Musielak--Orlicz Hardy Spaces and Their Applications

abstract: Let $(\Omega,\mathcal{F},\mathbb{P})$ be a probability space and $\varphi:\ \Omega\times[0,\infty)\to [0,\infty)$ be a Musielak--Orlicz function. In this article, the authors establish the atomic characterizations of weak martingale Musielak--Orlicz Hardy spaces $WH_{\varphi}^s(\Omega)$, $WH_{\varphi}^M(\Omega)$, $WH_{\varphi}^S(\Omega)$, $WP_{\varphi}(\Omega)$ and $WQ_{\varphi}(\Omega)$. Using these atomic characterizations, the authors then obtain the boundedness of sublinear operators from weak martingale Musielak--Orlicz Hardy spaces to weak Musielak--Orlicz spaces, and some martingale inequalities which further clarify the relationships among $WH_{\varphi}^s(\Omega)$, $WH_{\varphi}^M(\Omega)$, $WH_{\varphi}^S(\Omega)$, $WP_{\varphi}(\Omega)$ and $WQ_{\varphi}(\Omega)$. All these results improve and generalize the corresponding results on weak martingale Orlicz--Hardy spaces. Moreover, the authors also improve all the known results on weak martingale Musielak--Orlicz Hardy spaces. In particular, both the boundedness of sublinear operators and the martingale inequalities, for the weak weighted martingale Hardy spaces as well as for the weak weighted martingale Orlicz--Hardy spaces, are new.


title: A Method for Robust Online Classification using Dictionary Learning: Development and Assessment for Monitoring Manual Material Handling Activities Using Wearable Sensors

abstract: Classification methods based on sparse estimation have drawn much attention recently, due to their effectiveness in processing high-dimensional data such as images. In this paper, a method to improve the performance of a sparse representation classification (SRC) approach is proposed; it is then applied to the problem of online process monitoring of human workers, specifically manual material handling (MMH) operations monitored using wearable sensors (involving 111 sensor channels). Our proposed method optimizes the design matrix (aka dictionary) in the linear model used for SRC, minimizing its ill-posedness to achieve a sparse solution. This procedure is based on the idea of dictionary learning (DL): we optimize the design matrix formed by training datasets to minimize both redundancy and coherency as well as reducing the size of these datasets. Use of such optimized training data can subsequently improve classification accuracy and help decrease the computational time needed for the SRC; it is thus more applicable for online process monitoring. Performance of the proposed methodology is demonstrated using wearable sensor data obtained from manual material handling experiments, and is found to be superior to those of benchmark methods in terms of accuracy, while also requiring computational time appropriate for MMH online monitoring.


title: Highly accurate acoustic scattering: Isogeometric Analysis coupled with local high order Farfield Expansion ABC

abstract: This work is concerned with a unique combination of high order local absorbing boundary conditions (ABC) with a general curvilinear Finite Element Method (FEM) and its implementation in Isogeometric Analysis (IGA) for time-harmonic acoustic waves. The ABC employed were recently devised by Villamizar, Acosta and Dastrup [J. Comput. Phys. 333 (2017) 331] . They are derived from exact Farfield Expansions representations of the outgoing waves in the exterior of the regions enclosed by the artificial boundary. As a consequence, the error due to the ABC on the artificial boundary can be reduced conveniently such that the dominant error comes from the volume discretization method used in the interior of the computational domain. Reciprocally, the error in the interior can be made as small as the error at the artificial boundary by appropriate implementation of {\it p-} and {\it h}- refinement. We apply this novel method to cylindrical, spherical and arbitrary shape scatterers including a prototype submarine. Our numerical results exhibits spectral-like approximation and high order convergence rate. Additionally, they show that the proposed method can reduce both the pollution and artificial boundary errors to negligible levels even in very low- and high- frequency regimes with rather coarse discretization densities in the IGA. As a result, we have developed a highly accurate computational platform to numerically solve time-harmonic acoustic wave scattering in two- and three-dimensions.


title: Quadrupole Phonons in the Cadmium Isotopes

abstract: A key question concerning the spherical-vibrator attributes of states in Cadmium isotopes is addressed by means of a boson Hamiltonian encompassing U(5) partial dynamical symmetry. The U(5) symmetry is preserved in a segment of the spectrum and is broken in particular non-yrast states, and the resulting mixing with the intruder states is small. The vibrational character is thus maintained in the majority of low-lying normal states, as observed in $^{110}$Cd.


title: Type-II Seesaw Scalar Triplet Model at a 100 TeV $pp$ Collider: Discovery and Higgs Portal Coupling Determination

abstract: We investigate the collider phenomenology of the scalar triplet particles in the Type-II seesaw model at a 100 TeV $pp$ collider. Depending on triplet vacuum expectation value $v_\Delta$, the dominant discovery channels could be $H^{++} H^{--}$ and $H^{\pm\pm}H^\mp$. We find the $H^{\pm\pm}H^\mp\to W^\pm W^\pm h W^\mp/\ell^\pm \ell^\pm h W^\mp$ channels are promising for both model discovery at relatively large $v_\Delta$ and determination of the Higgs portal couplings $\lambda_4$ and $\lambda_5$. We also find that these two channels are complementary to indirect determination of $\lambda_4$ from future precise measurements on $h\to\gamma\gamma$ decay rate. Together with pair production of the doubly-charged Higgs subsequently decaying into same-sign di-leptons, the $H^{\pm\pm}H^\mp$ channels have the potential to cover a significant portion of the parameter space of the Type-II seesaw complex scalar triplet model.


title: Revised Wavelength and Spectral Response Calibrations for AKARI Near-Infrared Grism Spectroscopy: Post-Cryogenic Phase

abstract: We present a new calibration for the second-order light contamination in the near-infrared grism spectroscopy with the Infrared Camera aboard AKARI, specifically for the post-cryogenic phase of the satellite (Phase 3). Following our previous work on the cryogenic phase (Phases 1 and 2), the wavelength and spectral response calibrations were revised. Unlike Phases 1 and 2, during Phase 3, the temperature of the instrument was not stable and gradually increased from 40 to 47 K. To assess the effect of the temperature increase, we divided Phase 3 into three sub-phases and performed the calibrations separately. As in Phases 1 and 2, we confirmed that there was contamination due to the wavelength dependence of the refractive index of the grism material in every sub-phase. The wavelength calibration curves for the three sub-phases coincided with each other and did not show any significant temperature dependence. The response decreased with temperature by ~10% from the beginning to the end of Phase 3. We approximated the temperature dependence of the response as a linear relation and derived a correction factor as a function of temperature. The relative fraction of the second-order light contamination to the first-order light was found to be by 25% smaller than in Phases 1 and 2.


title: New Axion Searches at Flavor Factories

abstract: We assess the impact of searches at flavor factories for new neutral resonances that couple to both photons and gluons. These are well motivated by "heavy axion" solutions of the strong CP problem and by frameworks addressing both Dark Matter and the Higgs hierarchy problem. We use LHCb public diphoton data around the Bs mass to derive the current best limit on these resonances for masses between 4.9 and 6.3 GeV. We estimate that a future LHCb dedicated search would test an axion decay constant of O(TeV) for axion masses in the few-to-tens of GeV, being fully complementary to the low mass ATLAS and CMS searches. We also derive the impact of BABAR searches based on Upsilon decays and the future Belle-II reach.


title: Reanalyzing an Evaporating Black Hole

abstract: A coherent picture of the quantum mechanics of a collapse-formed, evaporating black hole is presented. In a distant frame, semiclassical theory in the zone describes microscopic dynamics of only the "hard modes," the modes that are hard enough to be discriminated in the timescale of Hawking emission. The thermal nature of these modes arises from microcanonical typicality of the full black hole degrees of freedom, mostly composed of the "soft modes," the modes that cannot be discriminated at the semiclassical level. The hard modes are purified by a combined system of the soft modes and early Hawking radiation, but not by either of them separately. This intrinsically tripartite structure of entanglement is general, regardless of the age of the black hole. The interior spacetime emerges only at a coarse-grained level. To describe it, an effective theory can be erected at each time, which applies only to a limited spacetime region determined by the time at which the theory is erected. The entire interior of the black hole can be described only using multiple effective theories erected at different times, realizing the idea of complementarity. We analyze implications of the entanglement structure described here for various phenomena, including Hawking evaporation and general information retrieval. For multiple entangled black holes, it implies that semiclassical objects dropped into different black holes cannot meet in the interior, although each object smoothly enters the horizon of the black hole to which it is falling. We also discuss physics in Rindler space, elucidating how it is obtained as a smooth limit of the black hole physics.


title: Constraints on the common-envelope evolution process from wide triple systems

abstract: Common envelope (CE) is an important phase in the evolution of interacting evolved binary systems. The interaction of the binary components during the CE evolution (CEE) stage gives rise to orbital inspiral and the formation of a short-period binary or a merger, on the expense of extending and/or ejecting the envelope. CEE is not well understood, as hydrodynamical simulations show that only a fraction of the CE-mass is ejected during the dynamical inspiral, in contrast with observations of post-CE binaries. Different CE models suggest different timescales are involved in the CE-ejection, and hence a measurement of the CE-ejection timescale could provide direct constraints on the CEE-process. Here we propose a novel method for constraining the mass-loss timescale of the CE, using post-CE binaries which are part of wide-orbit triple systems. The orbit/existence of a third companion constrains the CE mass-loss timescale, since rapid CE mass-loss may disrupt the triple system, while slower CE mass-loss may change the orbit of the third companion without disrupting it. As first test-cases we examine two observed post-CE binaries in wide triples, Wolf-1130 and GD-319. We follow their evolution due to mass-loss using analytic and numerical tools, and consider different mass-loss functions. We calculate a wide grid of binary parameters and mass-loss timescales in order to determine the most probable mass-loss timescale leading to the observed properties of the systems. We find that mass-loss timescales of the order of $10^{3}-10^{5}{\rm yr}$ are the most likely to explain these systems. Such long timescales are in tension with most of the CE mass-loss models, which predict shorter, dynamical timescales, but are potentially consistent with the longer timescales expected from the dust-driven winds model for CE ejection.


title: Enhanced n-body annihilation of dark matter and its indirect signatures

abstract: We examine the possible indirect signatures of dark matter annihilation processes with a non-standard scaling with the dark matter density, and in particular the case where more than two dark matter particles participate in the annihilation process. We point out that such processes can be strongly enhanced at low velocities without violating unitarity, similar to Sommerfeld enhancement in the standard case of two-body annihilation, potentially leading to visible signals in indirect searches. We study in detail the impact of such multi-body annihilations on the ionization history of the universe and consequently the cosmic microwave background, and find that unlike in the two-body case, the dominant signal can naturally arise from the end of the cosmic dark ages, after the onset of structure formation. We examine the complementary constraints from the Galactic Center, Galactic halo, and galaxy clusters, and outline the circumstances under which each search would give rise to the strongest constraints. We also show that if there is a population of ultra-compact dense dark matter clumps present in the Milky Way with sufficiently steep density profile, then it might be possible to detect point sources illuminated by multi-body annihilation, even if there is no large low-velocity enhancement. Finally, we provide a case study of a model where 3-body annihilation dominates the freezeout process, and in particular the resonant regime where a large low-velocity enhancement is naturally generated.


title: Dark Energy Survey Year 1 Results: Methods for Cluster Cosmology and Application to the SDSS

abstract: We perform the first blind analysis of cluster abundance data. Specifically, we derive cosmological constraints from the abundance and weak-lensing signal of \redmapper\ clusters of richness $\lambda\geq 20$ in the redshift range $z\in[0.1,0.3]$ as measured in the Sloan Digital Sky Survey (SDSS). We simultaneously fit for cosmological parameters and the richness--mass relation of the clusters. For a flat $\Lambda$CDM cosmological model with massive neutrinos, we find $S_8 \equiv \sigma_{8}(\Omega_m/0.3)^{0.5}=0.79^{+0.05}_{-0.04}$. This value is both consistent and competitive with that derived from cluster catalogues selected in different wavelengths. Our result is also consistent with the combined probes analyses by the Dark Energy Survey (DES) and the Kilo-Degree Survey (KiDS), and with the Cosmic Microwave Background (CMB) anisotropies as measured by \planck. We demonstrate that the cosmological posteriors are robust against variation of the richness--mass relation model and to systematics associated with the calibration of the selection function. In combination with Baryon Acoustic Oscillation (BAO) data and Big-Bang Nucleosynthesis (BBN) data, we constrain the Hubble rate to be $h=0.66\pm 0.02$, independent of the CMB. Future work aimed at improving our understanding of the scatter of the richness--mass relation has the potential to significantly improve the precision of our cosmological posteriors. The methods described in this work were developed for use in the forthcoming analysis of cluster abundances in the DES. Our SDSS analysis constitutes the first part of a staged-unblinding analysis of the full DES data set.


title: SCExAO/CHARIS Near-Infrared Direct Imaging, Spectroscopy, and Forward-Modeling of kappa And b: A Likely Young, Low-Gravity Superjovian Companion

abstract: We present SCExAO/CHARIS high-contrast imaging/$JHK$ integral field spectroscopy of $\kappa$ And b, a directly-imaged low-mass companion orbiting a nearby B9V star. We detect $\kappa$ And b at a high signal-to-noise and extract high precision spectrophotometry using a new forward-modeling algorithm for (A-)LOCI complementary to KLIP-FM developed by Pueyo (2016). $\kappa$ And b's spectrum best resembles that of a low-gravity L0--L1 dwarf (L0--L1$\gamma$). Its spectrum and luminosity are very well matched by 2MASSJ0141-4633 and several other 12.5--15 $M_{\rm J}$ free floating members of the 40 $Myr$-old Tuc-Hor Association, consistent with a system age derived from recent interferometric results for the primary, a companion mass at/near the deuterium-burning limit (13$^{+12}_{-2}$ M$_{\rm J}$), and a companion-to-primary mass ratio characteristic of other directly-imaged planets ($q$ $\sim$ 0.005$^{+0.005}_{-0.001}$). We did not unambiguously identify additional, more closely-orbiting companions brighter and more massive than $\kappa$ And b down to $\rho$ $\sim$ 0.3" (15 au). SCExAO/CHARIS and complementary Keck/NIRC2 astrometric points reveal clockwise orbital motion. Modeling points towards a likely eccentric orbit: a subset of acceptable orbits include those that are aligned with the star's rotation axis. However, $\kappa$ And b's semimajor axis is plausibly larger than 75 au and in a region where disk instability could form massive companions. Deeper $\kappa$ And high-contrast imaging and low-resolution spectroscopy from extreme AO systems like SCExAO/CHARIS and higher resolution spectroscopy from Keck/OSIRIS or, later, IRIS on the Thirty Meter Telescope could help clarify $\kappa$ And b's chemistry and whether its spectrum provides an insight into its formation environment.


title: SCExAO/CHARIS Near-IR High-Contrast Imaging and Integral Field Spectroscopy of the HIP 79977 Debris Disk

abstract: We present new, near-infrared (1.1--2.4 $\mu m$) high-contrast imaging of the bright debris disk surrounding HIP 79977 with the Subaru Coronagraphic Extreme Adaptive Optics system (SCExAO) coupled with the CHARIS integral field spectrograph. SCExAO/CHARIS resolves the disk down to smaller angular separations of (0.11"; $r \sim 14$ au) and at a higher significance than previously achieved at the same wavelengths. The disk exhibits a marginally significant east-west brightness asymmetry in $H$ band that requires confirmation. Geometrical modeling suggests a nearly edge-on disk viewed at a position angle of $\sim$ 114.6$^{o}$ east of north. The disk is best-fit by scattered-light models assuming strongly forward-scattering grains ($g$ $\sim$ 0.5--0.65) confined to a torus with a peak density at $r_{0}$ $\sim$ 53--75 au. We find that a shallow outer density power law of $\alpha_{out}=$-1-- -3 and flare index of $\beta = 1$ are preferred. Other disk parameters (e.g.~inner density power law and vertical scale height) are more poorly constrained. The disk has a slightly blue intrinsic color and its profile is broadly consistent with predictions from birth ring models applied to other debris disks. While HIP 79977's disk appears to be more strongly forward-scattering than most resolved disks surrounding 5--30 Myr-old stars, this difference may be due to observational biases favoring forward-scattering models for inclined disks vs. lower inclination, ostensibly neutral-scattering disks like HR 4796A's. Deeper, higher signal-to-noise SCExAO/CHARIS data can better constrain the disk's dust composition.


title: Fast radio burst source properties from polarization measurements

abstract: Recent polarization measurements of fast radio bursts (FRBs) provide new insights on these enigmatic sources. We show that the nearly 100% linear polarization and small variation of the polarization position angles (PAs) of multiple bursts from the same source suggest that the radiation is produced near the surface of a strongly magnetized neutron star. As the emitted radiation travels through the magnetosphere, the electric vector of the X-mode wave adiabatically rotates and stays perpendicular to the local magnetic field direction. The PA freezes at a radius where the plasma density becomes too small to be able to turn the electric vector. At the freeze-out radius, the electric field is perpendicular to the magnetic dipole moment of the neutron star projected in the plane of the sky, independent of the radiation mechanism or the orientation of the magnetic field in the emission region. We discuss a number of predictions of the model. The variation of PAs from repeating FRBs should follow the rotational period of the underlying neutron star (but the burst occurrence may not be periodic). Measuring this period will provide crucial support for the neutron star nature of the progenitors of FRBs. For FRB 121102, the small range of PA variation means that the magnetic inclination angle is less than about 20 degrees and that the observer's line of sight is outside the magnetic inclination cone. Other repeating FRBs may have a different range of PA variation from that of FRB 121102, depending on the magnetic inclination and the observer's viewing angle.


title: A Simple Non-equilibrium Feedback Model for Galaxy-Scale Star Formation: Delayed Feedback and SFR Scatter

abstract: We explore a class of simple non-equilibrium star formation models within the framework of a feedback-regulated model of the ISM, applicable to kiloparsec-scale resolved star formation relations (e.g. Kennicutt-Schmidt). Combining a Toomre-Q-dependent local star formation efficiency per free-fall time with a model for delayed feedback, we are able to match the normalization and scatter of resolved star formation scaling relations. In particular, this simple model suggests that large ($\sim$dex) variations in star formation rates (SFRs) on kiloparsec scales may be due to the fact that supernova feedback is not instantaneous following star formation. The scatter in SFRs at constant gas surface density in a galaxy then depends on the properties of feedback and when we observe its star-forming regions at various points throughout their collapse/star formation "cycles". This has the following important observational consequences: (1) the scatter and normalization of the Kennicutt-Schmidt relation are relatively insensitive to the local (small-scale) star formation efficiency, (2) but gas depletion times and velocity dispersions are; (3) the scatter in and normalization of the Kennicutt-Schmidt relation is a sensitive probe of the feedback timescale and strength; (4) even in a model where $\tilde Q_{\rm gas}$ deterministically dictates star formation locally, time evolution, variation in local conditions (e.g., gas fractions and dynamical times), and variations between galaxies can destroy much of the observable correlation between SFR and $\tilde Q_{\rm gas}$ in resolved galaxy surveys. Additionally, this model exhibits large scatter in SFRs at low gas surface densities, in agreement with observations of flat outer HI disk velocity dispersion profiles.


title: A-type stars in the Canada-France Imaging Survey II. Tracing the height of the disc at large distances with Blue Stragglers

abstract: We present the kinematics of Blue Straggler (BS) stars identified in the Canada-France-Imaging-Survey (CFIS), covering 4000 deg$^2$ on the sky in the $u$-band. The BS sample, characterised through CFIS and Pan-STARRS photometry, has been kinematically decomposed into putative halo and disc populations after cross-matching with Gaia astrometry and SDSS/SEGUE/LAMOST spectroscopy. This decomposition clearly reveals the strong flaring of the outer Milky Way disc. In particular, we show that we can detect this flaring up to a vertical height of $|Z| \simeq 8$ kpc at a Galactocentric distance of $R\sim $27 kpc. While some small level of flaring is expected for extended discs built up by radial migration, we demonstrate that this observed very strong flaring of the Milky Way disc is more likely consistent with it being dynamically heated by the repeated passage of the Sagittarius dwarf spheroidal galaxy through the midplane.


title: Precision Predictions at N$^3$LO for the Higgs Boson Rapidity Distribution at the LHC

abstract: We present precise predictions for the Higgs boson rapidity distribution at the LHC in the gluon fusion production mode. Our approach relies on the fully analytic computation of six terms in a systematic expansion of the partonic differential cross section around the production threshold of the Higgs boson at next-to-next-to-next-to leading order (N$^3$LO) in QCD perturbation theory. We observe a mild correction compared to the previous perturbative order and a significant reduction of the dependence of the cross section on the perturbative scale throughout the entire rapidity range.


title: Primordial gravitational wave phenomenology with polarized Sunyaev Zel'dovich tomography

abstract: The detection and characterization of primordial gravitational waves through their impact on the polarization anisotropies of the cosmic microwave background (CMB) is a primary science goal of current and future observations of the CMB. An ancillary dataset that will become accessible with the great leaps in sensitivity of CMB experiments is the polarized Sunyaev Zel'dovich (pSZ) effect, small-scale CMB polarization anisotropies induced by scattering from free electrons in the post-reionization Universe. The cross correlation of the pSZ effect with galaxy surveys, a technique known as pSZ tomography, can be used to reconstruct the remote quadrupole field: the CMB quadrupole observed from different locations in the Universe. Primordial gravitational waves leave a distinct imprint on the remote quadrupole field, making pSZ tomography a potential new method to characterize their properties. Building on previous work, we explore the utility of the full set of correlations between the primary CMB and the reconstructed remote quadrupole field to both provide exclusion limits on the amplitude of primordial gravitational waves, as well as to provide constraints on several phenomenological models of the tensor sector: axion gauge field inflation, general models with chiral tensors, and models with modified late-time decay of tensors. We find that relatively futuristic experimental requirements are necessary to provide competitive exclusion limits compared with the primary CMB. However, pSZ tomography can be a powerful probe of the late-time evolution of tensors and, through cross-correlations with the primary CMB, can provide mild improvements on parameter constraints in various models with chiral primordial gravitational waves.


title: Chandra high-resolution spectra of 4U~1630-47: the disappearance of the wind

abstract: We present the analysis of six {\it Chandra} X-ray high-resolution observations of the black hole low-mass X-ray binary 4U~1630-47 taken during its 2012-2013 outburst. {\rm Fe}~{\sc XXVI} K$\alpha$, K$\beta$, {\rm Fe}~{\sc XXV} K$\alpha$, K$\beta$ and {\rm Ca}~{\sc XX} K$\alpha$ blueshifted absorption lines were identified in the first four observations, which correspond to soft accretion states. The remaining observations, associated to intermediate and possibly hard accretion states, do not show significant absorption features down to equivalent width of 1 eV for both {\rm Fe}~{\sc XXVI} and {\rm Fe}~{\sc XXV}. We inferred wind launching radii between $1.2- 2.0$ ($10^{12}$ cm$/n$)$ \times 10^{11}$~cm and column densities $N({\rm H})> 10^{23}$ cm$^{-2}$. In the first four observations we found that thermal pressure is likely to be the dominant launching mechanism for the wind, although such conclusions depend on the assumed density. We used the spectral energy distributions obtained from our continuum modeling to compute thermal stability curves for all observations using the {\sc xstar} photoionization code. We found that the absence of lines in the transitional state cannot be attributed to an evolution of the plasma caused by thermal instabilities derived from the change in the continuum spectrum. In contrast, the disappearance of the wind could indicate an acceleration of the flow or that the plasma has been exhausted during the soft state.


title: Variability-selected low-luminosity active galactic nuclei candidates in the 7 Ms Chandra Deep Field-South

abstract: In deep X-ray surveys, active galactic nuclei (AGNs) with a broad range of luminosities have been identified. However, cosmologically distant low-luminosity AGN (LLAGN, $L_{\mathrm{X}} \lesssim 10^{42}$ erg s$^{-1}$) identification still poses a challenge due to significant contamination from host galaxies. Based on the 7 Ms Chandra Deep Field-South (CDF-S) survey, the longest timescale ($\sim 17$ years) deep X-ray survey to date, we utilize an X-ray variability selection technique to search for LLAGNs that remain unidentified among the CDF-S X-ray sources. We find 13 variable sources from 110 unclassified CDF-S X-ray sources. Except for one source which could be an ultraluminous X-ray source, the variability of the remaining 12 sources is most likely due to accreting supermassive black holes. These 12 AGN candidates have low intrinsic X-ray luminosities, with a median value of $7 \times10^{40}$ erg s$^{-1}$. They are generally not heavily obscured, with an average effective power-law photon index of 1.8. The fraction of variable AGNs in the CDF-S is independent of X-ray luminosity and is only restricted by the total number of observed net counts, confirming previous findings that X-ray variability is a near-ubiquitous property of AGNs over a wide range of luminosities. There is an anti-correlation between X-ray luminosity and variability amplitude for high-luminosity AGNs, but as the luminosity drops to $\lesssim 10^{42}$ erg s$^{-1}$, the variability amplitude no longer appears dependent on the luminosity. The entire observed luminosity-variability trend can be roughly reproduced by an empirical AGN variability model based on a broken power-law power spectral density function.


title: The Circular Velocity Curve of the Milky Way from $5$ to $25$ kpc

abstract: We measure the circular velocity curve $v_{\rm c}(R)$ of the Milky Way with the highest precision to date across Galactocentric distances of $5\leq R \leq 25$ kpc. Our analysis draws on the $6$-dimensional phase-space coordinates of $\gtrsim 23,000$ luminous red-giant stars, for which we previously determined precise parallaxes using a data-driven model that combines spectral data from APOGEE with photometric information from WISE, 2MASS, and Gaia. We derive the circular velocity curve with the Jeans equation assuming an axisymmetric gravitational potential. At the location of the Sun we determine the circular velocity with its formal uncertainty to be $v_{\rm c}(R_{\odot}) = (229.0\pm0.2)\rm\,km\,s^{-1}$ with systematic uncertainties at the $\sim 2-5\%$ level. We find that the velocity curve is gently but significantly declining at $(-1.7\pm0.1)\rm\,km\,s^{-1}\,kpc^{-1}$, with a systematic uncertainty of $0.46\rm\,km\,s^{-1}\,kpc^{-1}$, beyond the inner $5$ kpc. We exclude the inner $5$ kpc from our analysis due to the presence of the Galactic bar, which strongly influences the kinematic structure and requires modeling in a non-axisymmetric potential. Combining our results with external measurements of the mass distribution for the baryonic components of the Milky Way from other studies, we estimate the Galaxy's dark halo mass within the virial radius to be $M_{\rm vir} = (7.25\pm0.26)\cdot 10^{11}M_{\odot}$ and a local dark matter density of $\rho_{\rm dm}(R_{\odot}) = 0.30\pm0.03\,\rm GeV\,cm^{-3}$.


title: Twin Turtles

abstract: We present an ultraviolet extension of the Twin Higgs in which the radial mode of twin symmetry breaking is itself a pseudo-goldstone boson. This "turtle" structure raises the scale of new colored particles in exchange for additional states in the Higgs sector, making multiple Higgs-like scalars the definitive signature of naturalness in this context. We explore the parametrics and phenomenology of a concrete Twin Turtle model and demonstrate its robustness in two different supersymmetric completions. Along the way, we also introduce a new mechanism for inducing hard twin symmetry-breaking quartics via soft supersymmetry breaking.


title: Spectrophotometric parallaxes with linear models: Accurate distances for luminous red-giant stars

abstract: With contemporary infrared spectroscopic surveys like APOGEE, red-giant stars can be observed to distances and extinctions at which Gaia parallaxes are not highly informative. Yet the combination of effective temperature, surface gravity, composition, and age - all accessible through spectroscopy - determines a giant's luminosity. Therefore spectroscopy plus photometry should enable precise spectrophotometric distance estimates. Here we use the APOGEE-Gaia-2MASS-WISE overlap to train a data-driven model to predict parallaxes for red-giant branch stars with $0<\log g\leq2.2$ (more luminous than the red clump). We employ (the exponentiation of) a linear function of APOGEE spectral pixel intensities and multi-band photometry to predict parallax spectrophotometrically. The model training involves no logarithms or inverses of the Gaia parallaxes, and needs no cut on the Gaia parallax signal-to-noise ratio. It includes an L1 regularization to zero out the contributions of uninformative pixels. The training is performed with leave-out subsamples such that no star's astrometry is used even indirectly in its spectrophotometric parallax estimate. The model implicitly performs a reddening and extinction correction in its parallax prediction, without any explicit dust model. We assign to each star in the sample a new spectrophotometric parallax estimate; these parallaxes have uncertainties of a few to 15 percent, depending on data quality, which is more precise than the Gaia parallax for the vast majority of targets, and certainly any stars more than a few kpc distance. We obtain 10-percent distance estimates out to heliocentric distances of $20\,$kpc, and make global maps of the Milky Way's disk.


title: Fusing Binary Interface Defects in Topological Phases: The $\operatorname{Vec}(\mathbb{Z}/p\mathbb{Z})$ case

abstract: A binary interface defect is any interface between two (not necessarily invertible) domain walls. We compute all possible binary interface defects in Kitaev's $\mathbb{Z}/p\mathbb{Z}$ model and all possible fusions between them. Our methods can be applied to any Levin-Wen model. We also give physical interpretations for each of the defects in the $\mathbb{Z}/p\mathbb{Z}$ model. These physical interpretations provide a new graphical calculus which can be used to compute defect fusion.

categories: math.QA cond-mat.str-el math-ph math.MP quant-ph


title: SIToolBox : A package for Bayesian estimation of the isotropy violation in the CMB sky

abstract: The standard model of cosmology predicts a statistically isotropic (SI) CMB sky. However, the SI violation signals are always present in an observed sky-map. Different cosmological artifacts, measurement effects and unavoidable effects during data analysis etc. may lead to isotropy violation in an otherwise SI sky. Therefore, a proper data analysis technique should calculate all these SI violation signals, so that they can be matched with SI violation signals from the known sources and then conclude if there is any intrinsic SI violation in the CMB sky. In one of our recent works, we presented a general Bayesian formalism for measuring the isotropy violation signals in the CMB sky in presence of an idealized isotropic noise. In this paper, we have extended the mechanism and develop a software package, SIToolBox, for measuring SI violation in presence of anisotropic noise and masking.


title: Mixed WIMP-axion dark matter

abstract: We study the experimental constraints on a model of a two-component dark matter, consisting of the QCD axion, and a scalar particle, both contributing to the dark matter relic abundance of the Universe. The global Peccei-Quinn symmetry of the theory can be spontaneously broken down to a residual $\mathbb{Z}_2$ symmetry, thereby identifying this scalar as a stable weakly interacting massive particle, i.e., a dark matter candidate, in addition to the axion. We perform a comprehensive study of the model using the latest data from dark matter direct and indirect detection experiments, as well as new physics searches at the Large Hadron Collider. We find that although the model is mostly constrained by the dark matter detection experiments, it is still viable around a small region of the parameter space where the scalar dark matter is half as heavy as the Standard Model Higgs. In this allowed region, the bounds from these experiments are evaded due to a cancellation mechanism in the dark matter--Higgs coupling. The collider search results, however, are shown to impose weak bounds on the model.


title: Towards a Classification of Two-Character Rational Conformal Field Theories

abstract: We provide a simple and general construction of infinite families of consistent, modular-covariant pairs of characters satisfying the basic requirements to describe two-character RCFT. These correspond to solutions of generic second-order modular linear differential equations. To find these solutions, we first construct "quasi-characters" from the Kaneko-Zagier equation and subsequent works by Kaneko and collaborators, together with coset dual generalisations that we provide in this paper. We relate our construction to the Hecke images recently discussed by Harvey and Wu.


title: Sensitivity of dark matter haloes to their accretion histories

abstract: We apply our recently proposed "quadratic genetic modification" approach to generating and testing the effects of alternative mass accretion histories for a single $\Lambda$CDM halo. The goal of the technique is to construct different formation histories, varying the overall contribution of mergers to the fixed final mass. This enables targeted studies of galaxy and dark matter halo formation's sensitivity to the smoothness of mass accretion. Here, we focus on two dark matter haloes, each with four different mass accretion histories. We find that the concentration of both haloes systematically decreases as their merger history becomes smoother. This causal trend tracks the known correlation between formation time and concentration parameters in the overall halo population. At fixed formation time, we further establish that halo concentrations are sensitive to the order in which mergers happen. This ability to study an individual halo's response to variations in its history is highly complementary to traditional methods based on emergent correlations from an extended halo population.


title: Unified Superfluid Dark Sector

abstract: We present a novel theory of a unified dark sector, where late-time cosmic acceleration emerges from the dark matter superfluid framework. The system is described by a superfluid mixture consisting of two distinguishable states with a small energy gap, such as the ground state and an excited state of dark matter. Given their contact in the superfluid, interaction between those states can happen, converting one state into the other. This long range interaction within the superfluid couples the two superfluid phonon species through a cosine potential motivated by Josephson/Rabi interactions. As a consequence of this potential, a new dynamics of late-time accelerated expansion emerges in this system, without the need of dark energy, coming from a universe containing only this two-state DM superfluid. Because the superfluid species are non-relativistic, their sound speeds remain suitably small throughout the evolution. We calculate the expansion history and growth of linear perturbations, and compare the results to $\Lambda$CDM cosmology. For the fiducial parameters studied here, the predicted expansion and growth function are close to those of $\Lambda$CDM, but the difference in the predicted growth rate is significant at late times. The present theory nicely complements the recent proposal of dark matter superfluidity to explain the empirical success of MOdified Newtonian Dynamics (MOND) on galactic scales, thus offering a unified framework for dark matter, dark energy, and MOND phenomenology.


title: Revisiting the X-ray - Mass scaling relations of Early-type Galaxies with the Mass of their Globular Cluster Systems as a Proxy for the Total Galaxy Mass

abstract: Using globular cluster (GC) kinematics and photometry data, we calibrate the scaling relation between the total galaxy mass (M_TOT including dark matter) and total globular cluster system mass (M_GCS) in a sample of 30 early-type galaxies (ETG), confirming a nearly linear relationship between the two physical parameters. Using samples of 83 and 57 ETGs, we investigate this scaling relation in conjunction with the previously known relations between M_TOT and the ISM X-ray luminosity and temperature, respectively. We confirm that M_GCS can be effectively used as a proxy of M_TOT. We further find that the L_X,GAS - M_TOT relation is far tighter in the subsample of core ETGs, when compared to cusp ETGs. In core ETGs (old, passively evolving stellar systems) M_TOT is significantly larger than the total stellar mass M_STAR and the correlation with the hot gas properties is driven by their dark matter mass M_DM. Cusp ETGs have typically lower LX,GAS than core ETGs. In cusp ETGs, for a given M_DM, higher L_X,GAS is associated with higher M_STAR, suggesting stellar feedback as an important secondary factor for heating the ISM. Using the M_GCS- M_TOT scaling relations we compare 272 ETGs with previous estimates of the stellar-to-halo mass relation of galaxies. Our model-independent estimate of M_TOT results in a good agreement around halo masses of 10^12 Mo, but suggest higher star formation efficiency than usually assumed both at the low and at the high halo mass ends.


title: SNEMO: Improved Empirical Models for Type Ia Supernovae

abstract: Type Ia supernova cosmology depends on the ability to fit and standardize observations of supernova magnitudes with an empirical model. We present here a series of new models of Type Ia Supernova spectral time series that capture a greater amount of supernova diversity than possible with the models that are currently customary. These are entitled SuperNova Empirical MOdels (\textsc{SNEMO}\footnote{https://snfactory.lbl.gov/snemo}). The models are constructed using spectrophotometric time series from $172$ individual supernovae from the Nearby Supernova Factory, comprising more than $2000$ spectra. Using the available observations, Gaussian Processes are used to predict a full spectral time series for each supernova. A matrix is constructed from the spectral time series of all the supernovae, and Expectation Maximization Factor Analysis is used to calculate the principal components of the data. K-fold cross-validation then determines the selection of model parameters and accounts for color variation in the data. Based on this process, the final models are trained on supernovae that have been dereddened using the Fitzpatrick and Massa extinction relation. Three final models are presented here: \textsc{SNEMO2}, a two-component model for comparison with current Type~Ia models; \textsc{SNEMO7}, a seven component model chosen for standardizing supernova magnitudes which results in a total dispersion of $0.100$~mag for a validation set of supernovae, of which $0.087$~mag is unexplained (a total dispersion of $0.113$~mag with unexplained dispersion of $0.097$~mag is found for the total set of training and validation supernovae); and \textsc{SNEMO15}, a comprehensive $15$ component model that maximizes the amount of spectral time series behavior captured.


title: Monorepos: A Multivocal Literature Review

abstract: Monorepos (Monolithic Repositories) are used by large companies, such as Google and Facebook, and by popular open-source projects, such as Babel and Ember. This study provides an overview on the definition and characteristics of monorepos as well as on their benefits and challenges. Thereupon, we conducted a multivocal literature review on mostly grey literature. Our findings are fourfold. First, monorepos are single repositories that contain multiple projects, related or unrelated, sharing the same dependencies. Second, centralization and standardization are some key characteristics. Third, the main benefits include simplified dependencies, coordination of cross-project changes, and easy refactoring. Fourth, code health, codebase complexity, and tooling investments for both development and execution are considered the main challenges.


title: Multi-trace Correlators from Permutations as Moduli Space

abstract: We study the $n$-point functions of scalar multi-trace operators in the $U(N_c)$ gauge theory with adjacent scalars, such as ${\cal N}=4$ super Yang-Mills, at tree-level by using finite group methods. We derive a set of formulae of the general $n$-point functions, valid for general $n$ and to all orders of $1/N_c$. In one formula, the sum over Feynman graphs becomes a topological partition function on $\Sigma_{0,n}$ with a discrete gauge group, which resembles closed string interactions. In another formula, a new skeleton reduction of Feynman graphs generates connected ribbon graphs, which resembles open string interaction. We define the moduli space ${\cal M}_{g,n}^{\rm gauge}$ from the space of skeleton-reduced graphs in the connected $n$-point function of gauge theory. This moduli space is a proper subset of ${\cal M}_{g,n}$ stratified by the genus, and its top component gives a simple triangulation of $\Sigma_{g,n}$.


title: Single Image Haze Removal using a Generative Adversarial Network

abstract: Traditional methods to remove haze from images rely on estimating a transmission map. When dealing with single images, this becomes an ill-posed problem due to the lack of depth information. In this paper, we propose an end-to-end learning based approach which uses a modified conditional Generative Adversarial Network to directly remove haze from an image. We employ the usage of the Tiramisu model in place of the classic U-Net model as the generator owing to its higher parameter efficiency and performance. Moreover, a patch based discriminator was used to reduce artefacts in the output. To further improve the perceptual quality of the output, a hybrid weighted loss function was designed and used to train the model. Experiments on synthetic and real world hazy images demonstrates that our model performs competitively with the state of the art methods.


title: Many-body exceptional points in colliding condensates

abstract: Exceptional points describe the coalescence of the eigenmodes of a non-Hermitian matrix. When an exceptional point occurs in the unitary evolution of a many-body system, it generically leads to a dynamical instability with a finite wavevector [N. Bernier \etal, Phys. Rev. Lett. 113, 065303 (2014)]. Here, we study exceptional points in the context of the counterflow instability of colliding Bose-Einstein condensates. We show that the instability of this system is due to an exceptional point in the Bogoliubov spectrum. We further clarify the connection of this effect to the Landau criterion of superfluidity and to the scattering of classical particles. We propose an experimental set-up to directly probe this exceptional point, and demonstrate its feasibility with the aid of numerical calculations. Our work fosters the observation of exceptional points in nonequilibrium many-body quantum systems.


title: Dynamical structure factor of the three-dimensional quantum spin liquid candidate NaCaNi$_2$F$_7$

abstract: We study the spin-1 pyrochlore material NaCaNi$_2$F$_7$ with a combination of molecular dynamics simulations, stochastic dynamical theory and linear spin wave theory. The dynamical structure factor from inelastic neutron scattering is well described with a near-ideal Heisenberg Hamiltonian incorporating small anisotropic terms {and weak second-neighbor interactions}. We find that all three approaches reproduce remarkably well the momentum dependence of the scattering intensity as well as its energy dependence with the exception of the lowest energies. These results are notable in that (i) the data show a complete lack of sharp quasiparticle excitations in momentum space over much, if not all, of the energy range; (ii) linear spin-wave theory appears to apply in a regime where it would be expected to fail for a number of reasons. We elucidate what underpins these surprises, and note that basic questions about the nature of quantum spin liquidity in such systems pose themselves as a result.


title: A Note on Indexing Point Sets for Approximate Bottleneck Distance Queries

abstract: The {\em bottleneck distance} is a natural measure of the distance between two finite point sets of equal cardinality, defined as the minimum over all bijections between the point sets of the maximum distance between any pair of points put in correspondence by the bijection. In this work, we consider the problem of building a data structure $\mathbb{D}$ that indexes a collection of $m$ planar point sets (of varying sizes) and supports nearest bottleneck distance queries: given a query point set $Q$ of size $n$, we would like to find the point set(s) $P \in \mathbb{D}$ of size $n$ that are closest in terms of bottleneck distance. Without loss of generality, we assume that all point sets belong to the unit box $[0,1]^2$ in the plane and focus on the $L_\infty$ norm, although the techniques can also be used for other norms. The main contribution is a {\em trie}-based data structure finds a $6$-approximate nearest neighbor in $O(-\lg(d_B(\mathbb{D},Q)) n)$ time, where $d_B(\mathbb{D},Q)$ is the minimum bottleneck distance from $Q$ to any point set in $\mathbb{D}$.


title: Incorporation of random alloy GaBi$_x$As$_{1-x}$ barriers in InAs quantum dot molecules (I): energy levels and confined hole states

abstract: Self-assembled InAs quantum dots (QDs), which have long hole-spin coherence times and are amenable to optical control schemes, have long been explored as building blocks for qubit architectures. One such design consists of vertically stacking two QDs to create a quantum dot molecule (QDM) and using the spin-mixing properties of "molecule-like" coupled hole states for all-optical qubit manipulation. In this article, the first of two papers, we introduce the incorporation of dilute GaBiAs alloys in the barrier region between the two dots. GaBiAs is expected to increase the spin-mixing of the molecular states needed for qubit operations by raising the barrier valence band edge and spin-orbit splitting. Using an atomistic tight-binding model, we compute the properties of GaBiAs and the modification of hole states that arise when the alloy is used in the barrier of an InAs QDM. An atomistic treatment is necessary to correctly capture non-traditional alloy effects such as the band-anticrossing valence band. It also allows for the study of configurational variances and clustering effects of the alloy. We find that in InAs QDMs with a GaBiAs interdot barrier, electron states are not strongly affected by the inclusion of Bi. However, hole states are much more sensitive to the presence and configuration of Bi in the barriers. By independently studying the alloy-induced strain and electronic scattering off Bi and As orbitals, we conclude that an initial increase in QDM hole state energy at low Bi concentration is caused by the alloy-induced strain. We further find that the decrease in QDM hole energy at higher Bi concentrations can only be explained when both alloy strain and orbital effects are considered. In our second article, we use the understanding developed here to discuss how the alloyed barriers contribute to enhancement in hole spin-mixing and the implications for QDM qubit architectures.


title: Optimal choice of state tomography quorum formed by projection operators

abstract: A minimal set of measurement operators for quantum state tomography has in the non-degenerate case ideally eigenbases which are mutually unbiased. This is different for the degenerate case. Here, we consider the situation where the measurement operators are projections on individual pure quantum states. This corresponds to maximal degeneracy. We present numerically optimized sets of projectors and find that they significantly outperform those which are taken from a set of mutually unbiased bases.


title: Scaling Up Cartesian Genetic Programming through Preferential Selection of Larger Solutions

abstract: We demonstrate how efficiency of Cartesian Genetic Programming method can be scaled up through the preferential selection of phenotypically larger solutions, i.e. through the preferential selection of larger solutions among equally good solutions. The advantage of the preferential selection of larger solutions is validated on the six, seven and eight-bit parity problems, on a dynamically varying problem involving the classification of binary patterns, and on the Paige regression problem. In all cases, the preferential selection of larger solutions provides an advantage in term of the performance of the evolved solutions and in term of speed, the number of evaluations required to evolve optimal or high-quality solutions. The advantage provided by the preferential selection of larger solutions can be further extended by self-adapting the mutation rate through the one-fifth success rule. Finally, for problems like the Paige regression in which neutrality plays a minor role, the advantage of the preferential selection of larger solutions can be extended by preferring larger solutions also among quasi-neutral alternative candidate solutions, i.e. solutions achieving slightly different performance.


title: Random walks on Convergence Groups

abstract: We extend some properties of random walks on hyperbolic groups to random walks on convergence groups. In particular we prove that if a convergence group $G$ acts on a compact metrizable space $M$ with the convergence property then we can provide $G\cup M$ with a compact topology such that random walks on $G$ converge almost surely to points in $M$. Furthermore we prove that if $G$ is finitely generated and the random walk has finite entropy and finite logarithmic moment with respect to the word metric, then $M$, with the corresponding hitting measure, can be seen as a model for the Poisson boundary of $G$.


title: Diagnostic Accuracy of Content Based Dermatoscopic Image Retrieval with Deep Classification Features

abstract: Background: Automated classification of medical images through neural networks can reach high accuracy rates but lack interpretability. Objectives: To compare the diagnostic accuracy obtained by using content based image retrieval (CBIR) to retrieve visually similar dermatoscopic images with corresponding disease labels against predictions made by a neural network. Methods: A neural network was trained to predict disease classes on dermatoscopic images from three retrospectively collected image datasets containing 888, 2750 and 16691 images respectively. Diagnosis predictions were made based on the most commonly occurring diagnosis in visually similar images, or based on the top-1 class prediction of the softmax output from the network. Outcome measures were area under the ROC curve for predicting a malignant lesion (AUC), multiclass-accuracy and mean average precision (mAP), measured on unseen test images of the corresponding dataset. Results: In all three datasets the skin cancer predictions from CBIR (evaluating the 16 most similar images) showed AUC values similar to softmax predictions (0.842, 0.806 and 0.852 versus 0.830, 0.810 and 0.847 respectively; p-value>0.99 for all). Similarly, the multiclass-accuracy of CBIR was comparable to softmax predictions. Networks trained for detecting only 3 classes performed better on a dataset with 8 classes when using CBIR as compared to softmax predictions (mAP 0.184 vs. 0.368 and 0.198 vs. 0.403 respectively). Conclusions: Presenting visually similar images based on features from a neural network shows comparable accuracy to the softmax probability-based diagnoses of convolutional neural networks. CBIR may be more helpful than a softmax classifier in improving diagnostic accuracy of clinicians in a routine clinical setting.


title: Hierarchical multi-class segmentation of glioma images using networks with multi-level activation function

abstract: For many segmentation tasks, especially for the biomedical image, the topological prior is vital information which is useful to exploit. The containment/nesting is a typical inter-class geometric relationship. In the MICCAI Brain tumor segmentation challenge, with its three hierarchically nested classes 'whole tumor', 'tumor core', 'active tumor', the nested classes relationship is introduced into the 3D-residual-Unet architecture. The network comprises a context aggregation pathway and a localization pathway, which encodes increasingly abstract representation of the input as going deeper into the network, and then recombines these representations with shallower features to precisely localize the interest domain via a localization path. The nested-class-prior is combined by proposing the multi-class activation function and its corresponding loss function. The model is trained on the training dataset of Brats2018, and 20% of the dataset is regarded as the validation dataset to determine parameters. When the parameters are fixed, we retrain the model on the whole training dataset. The performance achieved on the validation leaderboard is 86%, 77% and 72% Dice scores for the whole tumor, enhancing tumor and tumor core classes without relying on ensembles or complicated post-processing steps. Based on the same start-of-the-art network architecture, the accuracy of nested-class (enhancing tumor) is reasonably improved from 69% to 72% compared with the traditional Softmax-based method which blind to topological prior.


title: OGLE Collection of Galactic Cepheids

abstract: We present here a new major part of the OGLE Collection of Variable Stars - OGLE Collection of Galactic Cepheids. The new dataset was extracted from the Galaxy Variability Survey images - a dedicated large-scale survey of the Galactic disk and outer bulge conducted by the OGLE project since 2013. The OGLE collection contains 2721 Cepheids of all types - classical, type II and anomalous. It more than doubles the number of known Galactic classical Cepheids. Due to the long-term monitoring and large number of epochs the selected sample is very pure, generally free from contaminating stars of other types often mimicking Cepheids. Its completeness is high at 90% level for classical Cepheids - tested using recent samples of Galactic Cepheids: ASAS-SN, ATLAS, Gaia DR2 and Wise catalog of variable stars. Our comparisons indicate that the completeness of the two latter datasets, Gaia DR2 and Wise catalog, is very low, at < 10% level in the magnitude range of the OGLE GVS survey (10.8 < I < 19.5 mag). Both these samples are severely contaminated by non-Cepheids (the purity is 67% and 56%, respectively). We also present several interesting objects found in the new OGLE Collection - multi-mode pulsators, a first Galactic candidate for eclipsing system containing Cepheid, a binary Cepheid candidate.


title: On norm almost periodic measures

abstract: In this paper, we study norm almost periodic measures on locally compact Abelian groups. First, we show that the norm almost periodicity of $\mu$ is equivalent to the equi-Bohr almost periodicity of $\mu*g$ for all $g$ in a fixed family of functions. Then, we show that, for absolutely continuous measures, norm almost periodicity is equivalent to the Stepanov almost periodicity of the Radon--Nikodym density.


title: Nonlinear Evolution of the Resonant Drag Instability in Magnetized Gas

abstract: We investigate, for the first time, the nonlinear evolution of the magnetized "resonant drag instabilities" (RDIs). We explore magnetohydrodynamic (MHD) simulations of gas mixed with (uniform) dust grains subject to Lorentz and drag forces, using the GIZMO code. The magnetized RDIs exhibit fundamentally different behaviour than the purely acoustic RDIs. The dust organizes into coherent structures and the system exhibits strong dust-gas separation. In the linear and early nonlinear regime, the growth rates agree with linear theory and the dust self-organizes into two-dimensional planes or "sheets." Eventually the gas develops fully nonlinear, saturated Alfv\'enic and compressible fast-mode turbulence, which fills the under-dense regions with a small amount of dust, and drives a dynamo which saturates at equipartition of kinetic and magnetic energy. The dust density fluctuations exhibit significant non-Gaussianity, and the power spectrum is strongly weighted towards the largest (box-scale) modes. The saturation level can be understood via quasi-linear theory, as the forcing and energy input via the instabilities becomes comparable to saturated tension forces and dissipation in turbulence. The magnetized simulation presented here is just one case; it is likely that the magnetic RDIs can take many forms in different parts of parameter space.

categories: astro-ph.GA astro-ph.SR physics.flu-dyn physics.plasm-ph


title: A Central Limit Theorem for the stochastic heat equation

abstract: We consider the one-dimensional stochastic heat equation driven by a multiplicative space-time white noise. We show that the spatial integral of the solution from $-R$ to $R$ converges in total variance distance to a standard normal distribution as $R$ tends to infinity, after renormalization. We also show a functional version of this central limit theorem.


title: One- and two- neutron halo at the dripline. From 11Be to 11Li and back: 10Li and parity inversion

abstract: The nuclei 11Be and 11Li provide paradigmatic examples of one-and two- neutron halo systems. Because the reaction 1H(11Li,9Li)3H is dominated by successive transfer, one can use the quantitative picture emerging from a nu- clear field theory description of the structure and reaction mechanism of the above Cooper pair transfer process and of the 2H(10Be,11Be)1H and 1H(11Be,10Be)2H reactions, to shed light on the structure of 10Li. This analysis provides important support for a parity inverted scenario with a 1/2+ virtual state at about 0.2 MeV.


title: A Bayesian Nonparametric Model for Zero-Inflated Outcomes: Prediction, Clustering, and Causal Estimation

abstract: Researchers are often interested in predicting outcomes, conducting clustering analysis to detect distinct subgroups of their data, or computing causal treatment effects. Pathological data distributions that exhibit skewness and zero-inflation complicate these tasks - requiring highly flexible, data-adaptive modeling. In this paper, we present a fully nonparametric Bayesian generative model for continuous, zero-inflated outcomes that simultaneously predicts structural zeros, captures skewness, and clusters patients with similar joint data distributions. The flexibility of our approach yields predictions that capture the joint data distribution better than commonly used zero-inflated methods. Moreover, we demonstrate that our model can be coherently incorporated into a standardization procedure for computing causal effect estimates that are robust to such data pathologies. Uncertainty at all levels of this model flow through to the causal effect estimates of interest - allowing easy point estimation, interval estimation, and posterior predictive checks verifying positivity, a required causal identification assumption. Our simulation results show point estimates to have low bias and interval estimates to have close to nominal coverage under complicated data settings. Under simpler settings, these results hold while incurring lower efficiency loss than comparator methods. Lastly, we use our proposed method to analyze zero-inflated inpatient medical costs among endometrial cancer patients receiving either chemotherapy and radiation therapy in the SEER medicare database.


title: Molecular Mott state in the deficient spinel GaV$_4$S$_8$

abstract: In this study, we investigated theoretically the Mott-insulating phase of a deficient spinel chalcogenide GaV$_4$S$_8$, which is known to form a tetrahedral V$_4$S$_4$ cluster unit that results in molecular orbitals (MOs) with a narrow bandwidth in the noninteracting limit. We used a cluster extension of charge self-consistent embedded dynamical mean-field theory to study the impact of strong intra-cluster correlations on the spectral properties as well as the structural degrees of freedom of the system. We found that the strong tetrahedral clustering renders the atomic Mott picture ineffective, and that the resulting MO picture is essential to describe the Mott phase. It was also found that, while the spectral properties can be qualitatively described by the truncation of the Hilbert space down to the lowest-energy MO, a proper description of the structural degrees of freedom requires the inclusion of multi-MO correlations that span a larger energy window. Specifically, we found that the lowest-energy MO description overemphasizes the clustering tendency, while the inclusion of the Hund's coupling between the lower- and higher-energy MOs corrects this tendency, bringing the theoretically predicted crystal structure into good agreement with the experiment.


title: Two view constraints on the epipoles from few correspondences

abstract: In general it requires at least 7 point correspondences to compute the fundamental matrix between views. We use the cross ratio invariance between corresponding epipolar lines, stemming from epipolar line homography, to derive a simple formulation for the relationship between epipoles and corresponding points. We show how it can be used to reduce the number of required points for the epipolar geometry when some information about the epipoles is available and demonstrate this with a buddy search app.


title: Comparing Two Approaches in Heteroscedastic Regression Models

abstract: Recently, a generalized test approach is proposed by Sadooghi-alvandi et al. (2016) and a fiducial approach is proposed by Xu and Li (2018) to test the equality of coefficients in several regression models with unequal variances. In this paper, it is shown that the considered quantities in these approaches are identically distributed and therefore, these approaches are same. Also, this result satisfies for the one-way ANOVA problem.


title: Univariate Mean Change Point Detection: Penalization, CUSUM and Optimality

abstract: The problem of univariate mean change point detection and localization based on a sequence of $n$ independent observations with piecewise constant means has been intensively studied for more than half century, and serves as a blueprint for change point problems in more complex settings. We provide a complete characterization of this classical problem in a general framework in which the upper bound $\sigma^2$ on the noise variance, the minimal spacing $\Delta$ between two consecutive change points and the minimal magnitude $\kappa$ of the changes, are allowed to vary with $n$. We first show that consistent localization of the change points, when the signal-to-noise ratio $\frac{\kappa \sqrt{\Delta}}{\sigma} < \sqrt{\log(n)}$, is impossible. In contrast, when $\frac{\kappa \sqrt{\Delta}}{\sigma}$ diverges with $n$ at the rate of at least $\sqrt{\log(n)}$, we demonstrate that two computationally-efficient change point estimators, one based on the solution to an $\ell_0$-penalized least squares problem and the other on the popular wild binary segmentation algorithm, are both consistent and achieve a localization rate of the order $\frac{\sigma^2}{\kappa^2} \log(n)$. We further show that such rate is minimax optimal, up to a $\log(n)$ term.


title: A Comparative Study of Fruit Detection and Counting Methods for Yield Mapping in Apple Orchards

abstract: We present new methods for apple detection and counting based on recent deep learning approaches and compare them with state-of-the-art results based on classical methods. Our goal is to quantify performance improvements by neural network-based methods compared to methods based on classical approaches. Additionally, we introduce a complete system for counting apples in an entire row. This task is challenging as it requires tracking fruits in images from both sides of the row. We evaluate the performances of three fruit detection methods and two fruit counting methods on six datasets. Results indicate that the classical detection approach still outperforms the deep learning based methods in the majority of the datasets. For fruit counting though, the deep learning based approach performs better for all of the datasets. Combining the classical detection method together with the neural network based counting approach, we achieve remarkable yield accuracies ranging from 95.56% to 97.83%.


title: Full-$f$ gyrokinetic simulation of turbulence in a helical open-field-line plasma

abstract: Curvature-driven turbulence in a helical open-field-line plasma is investigated using electrostatic five-dimensional gyrokinetic continuum simulations in an all-bad-curvature helical-slab geometry. Parameters for a National Spherical Torus Experiment scrape-off-layer plasma are used in the model. The formation and convective radial transport of plasma blobs is observed, and it is shown that the radial particle-transport levels are several times higher than diffusive Bohm-transport estimates. By reducing the strength of the poloidal magnetic field, the profile of the heat flux to the divertor plate is observed to broaden.


title: SB 796: a high-velocity RRc star

abstract: We report here on a detailed study of a c-type RR Lyrae variable (RRc variable), SB 796, serendipitously discovered in a search of the WASP public data for stars that display large photometric periodic modulation. SB 796 displays a period of $P = 0.26585$ d and semi-amplitude of $\sim$ 0.1 mag. Comparison of the modulation shape and period with the detailed analysis of LMC variables indicates that SB 796 is an RRc variable. Gaia DR2 classification corroborated our result. Radial-velocity (RV) follow-up observations revealed a periodic variation consistent with a sine modulation, with a semi-amplitude of $5.6\pm0.2$ km/s, and a minimum at phase of maximum brightness. Similar amplitude and phase were previously seen in other RRc variables. The stellar averaged RV is $\sim 250$ km/s, turning SB 796 to be a high-velocity star, while its present position, as derived from the Gaia astrometry, is at $\sim 3.5$ kpc below the Galactic plane. Integration of the stellar Galactic motion shows that SB 796 oscillates at a range of 0.5--20 kpc Galacto-centric distance, passing near the Galactic center about three times in 1 Gyr. The Galactic radial motion takes SB 796 up and down the plane to a scale height of $\sim 10$ kpc. During its $\sim10$ Gyrs estimated life time, SB 796 therefore passed $\sim 30$ times near the Galactic center.


title: How to train your MAML

abstract: The field of few-shot learning has recently seen substantial advancements. Most of these advancements came from casting few-shot learning as a meta-learning problem. Model Agnostic Meta Learning or MAML is currently one of the best approaches for few-shot learning via meta-learning. MAML is simple, elegant and very powerful, however, it has a variety of issues, such as being very sensitive to neural network architectures, often leading to instability during training, requiring arduous hyperparameter searches to stabilize training and achieve high generalization and being very computationally expensive at both training and inference times. In this paper, we propose various modifications to MAML that not only stabilize the system, but also substantially improve the generalization performance, convergence speed and computational overhead of MAML, which we call MAML++.


title: Properties of a Generalized Divergence Related to Tsallis Relative Entropy

abstract: In this paper, we investigate the partition inequality, joint convexity, and Pinsker's inequality, for a divergence that generalizes the Tsallis Relative Entropy and Kullback-Leibler divergence. The generalized divergence is defined in terms of a deformed exponential function, which replaces the Tsallis $q$-exponential. We also constructed a family of probability distributions related to the generalized divergence. We found necessary and sufficient conditions for the partition inequality to be satisfied. A sufficient condition for the joint convexity was established. We proved that the generalized divergence satisfies the partition inequality, and is jointly convex, if, and only if, it coincides with the Tsallis relative entropy. As an application of partition inequality, a criterion for the Pinsker's inequality was found.


title: Direct writing of heterostructures in single atomically precise graphene nanoribbons

abstract: Precision control of interfacial structures and electronic properties is the key to the realization of functional heterostructures. Here, utilizing the scanning tunneling microscope (STM) both as a manipulation and characterization tool, we demonstrate the fabrication of a heterostructure in a single atomically precise graphene nanoribbon (GNR) and report its electronic properties. The heterostructure is made of a seven-carbon-wide armchair GNR and a lower band gap intermediate ribbon synthesized bottom-up from a molecular precursor on an Au substrate. The short GNR segments are directly written in the ribbon with an STM tip to form atomic precision intraribbon heterostructures. Based on STM studies combined with density functional theory calculations, we show that the heterostructure has a type-I band alignment, with manifestations of quantum confinement and orbital hybridization. Our finding demonstrates a feasible strategy to create a double barrier quantum dot structure with atomic precision for novel functionalities, such as negative differential resistance devices in GNR-based nanoelectronics.


title: What Does a Successful Postdoctoral Fellowship Publication Record Look Like?

abstract: Obtaining a prize postdoctoral fellowship in astronomy and astrophysics involves a number of factors, many of which cannot be quantified. One criterion that can be measured is the publication record of an applicant. The publication records of past fellowship recipients may, therefore, provide some quantitative guidance for future prospective applicants. We investigated the publication patterns of recipients of the NASA prize postdoctoral fellowships in the Hubble, Einstein, and Sagan programs from 2014 through 2017, using the NASA ADS reference system. We tabulated their publications at the point where fellowship applications were submitted, and we find that the 133 fellowship recipients in that time frame had a median of 6 +/- 2 first-author publications, and 14 +/- 6 co-authored publications. The full range of first author papers is 1 to 15, and for all papers ranges from 2 to 76, indicating very diverse publication patterns. Thus, while fellowship recipients generally have strong publication records, the distribution of both first-author and co-authored papers is quite broad; there is no apparent threshold of publications necessary to obtain these fellowships. We also examined the post-PhD publication rates for each of the three fellowship programs, between male and female recipients, across the four years of the analysis and find no consistent trends. We hope that these findings will prove a useful reference to future junior scientists.


title: Automatically Detecting Self-Reported Birth Defect Outcomes on Twitter for Large-scale Epidemiological Research

abstract: In recent work, we identified and studied a small cohort of Twitter users whose pregnancies with birth defect outcomes could be observed via their publicly available tweets. Exploiting social media's large-scale potential to complement the limited methods for studying birth defects, the leading cause of infant mortality, depends on the further development of automatic methods. The primary objective of this study was to take the first step towards scaling the use of social media for observing pregnancies with birth defect outcomes, namely, developing methods for automatically detecting tweets by users reporting their birth defect outcomes. We annotated and pre-processed approximately 23,000 tweets that mention birth defects in order to train and evaluate supervised machine learning algorithms, including feature-engineered and deep learning-based classifiers. We also experimented with various under-sampling and over-sampling approaches to address the class imbalance. A Support Vector Machine (SVM) classifier trained on the original, imbalanced data set, with n-grams, word clusters, and structural features, achieved the best baseline performance for the positive classes: an F1-score of 0.65 for the "defect" class and 0.51 for the "possible defect" class. Our contributions include (i) natural language processing (NLP) and supervised machine learning methods for automatically detecting tweets by users reporting their birth defect outcomes, (ii) a comparison of feature-engineered and deep learning-based classifiers trained on imbalanced, under-sampled, and over-sampled data, and (iii) an error analysis that could inform classification improvements using our publicly available corpus. Future work will focus on automating user-level analyses for cohort inclusion.


title: Identifying Tidal Disruption Events via Prior Photometric Selection of Their Preferred Hosts

abstract: A nuclear transient detected in a post-starburst galaxy or other quiescent galaxy with strong Balmer absorption is likely to be a Tidal Disruption Event (TDE). Identifying such galaxies within the planned survey footprint of the Large Synoptic Survey Telescope (LSST)---before a transient is detected---will make TDE classification immediate and follow-up more efficient. Unfortunately, spectra for identifying most such galaxies are unavailable, and simple photometric selection is ineffective; cutting on "green valley" UV/optical/IR colors produces samples that are highly contaminated and incomplete. Here we propose a new strategy using only photometric optical/UV/IR data from large surveys. Applying a machine learning Random Forest classifier to a sample of ~400k SDSS galaxies with GALEX and WISE photometry, including 13,592 quiescent Balmer-strong galaxies, we achieve 53-61% purity and 8-21% completeness, given the range in redshift. For the subset of 1299 post-starburst galaxies, we achieve 63-73% purity and 5-12% completeness. Given these results, the range of likely TDE and supernova rates, and that 36-75% of TDEs occur in quiescent Balmer-strong hosts, we estimate that 13-99% of transients observed in photometrically-selected host galaxies will be TDEs and that we will discover 119-248 TDEs per year with LSST. Using our technique, we present a new catalog of 67,484 candidate galaxies expected to have a high TDE rate, drawn from the SDSS, Pan-STARRS, DES, and WISE photometric surveys. This sample is 3.5x larger than the current SDSS sample of similar galaxies, thereby providing a new path forward for transient science and galaxy evolution studies.


title: A pragmatic approach to formal fundamental physics

abstract: A minimal practical formal structure for a fundamental theory is suggested. A mechanism that produces such a structure is reviewed. The proposed mechanism has possibilities of producing non-canonical phenomena in SU(2) and SU(3) quantum gauge theories. These might provide testable conditional predictions. One possibility is a vacuum condensate of SU(2) gauge fields associated to certain trajectories of the SU(2) Yang-Mills flow. Contents 1 Formal fundamental physics 1.1 Against Quantum Gravity 1.2 Against mathematical idealizations 2 A minimal practical formal structure 2.1 An effective QFT for distances > L and an effective S-matrix for distances < L, for observers at every scale L >> l_P 2.2 QFT renormalization group operates from smaller distance L to larger; S-matrix renormalization group operates from larger L to smaller 2.3 An S-matrix does not imply a hamiltonian 3 A mechanism that produces such a formal structure 3.1 Summary 3.2 2d-QFT of the string worldsheet 3.3 Effective string S-matrix with IR cutoff L 3.4 Effective 2d coupling constants 3.5 Implement the S-matrix renormalization group 3.6 Production of an effective QFT with UV cutoff L 3.7 Possible non-canonical degrees of freedom and couplings in SU(2) and SU(3) quantum gauge theory 3.8 2d winding modes and 2d instantons 3.9 Vacuum condensate of SU(2) Yang-Mills flow defects 4 To do Appendix. Notes on the line of thought A.1 Search for a mechanism that produces QFT A.2 Pragmatism and the S-matrix philosophy


title: Plane curves which are quantum homogeneous spaces

abstract: Let $\mathcal{C}$ be a decomposable plane curve over an algebraically closed field $k$ of characteristic 0. That is, $\mathcal{C}$ is defined in $k^2$ by an equation of the form $g(x) = f(y)$, where $g$ and $f$ are polynomials of degree at least 2. We use this data to construct 3 pointed Hopf algebras, $A(x,a,g)$, $A(y,b,f)$ and $A(g,f)$, in the first two of which $g$ [resp. $f$] are skew primitive central elements, and the third being a factor of the tensor product of the first two. We conjecture that $A(g,f)$ contains the coordinate ring $\mathcal{O}(\mathcal{C})$ of $\mathcal{C}$ as a quantum homogeneous space, and prove this when each of $g$ and $f$ has degree at most 5 or is a power of the variable. We obtain many properties of these Hopf algebras, and show that, for small degrees, they are related to previously known algebras. For example, when $g$ has degree 3 $A(x,a,g)$ is a PBW deformation of the localisation at powers of a generator of the downup algebra $A(-1,-1,0)$.


title: Monitoring & Mitigation of Delayed Voltage Recovery using {\mu}PMU Measurements with Reduced Distribution System Model

abstract: This paper proposes a new method to monitor and mitigate fault induced delayed voltage recovery (FIDVR) phenomenon in distribution systems using {\mu}PMU measurements in conjunction with a Reduced Distribution System Model (RDSM). The recovery time estimated from a dynamic analysis of the FIDVR is used to monitor its behavior and a linear optimization is formulated to control air conditioner loads and DER reactive power injection to mitigate the FIDVR severity. The RDSM is made up of several sub-models, each of which is analogous to the Composite Load Model (CLM) with selected parameters. The linear formulation in combination with the RDSM reduces the computation time, enabling online execution. Simulated {\mu}PMU measurements from the IEEE 37 node distribution system connected to the IEEE 9 bus system under various fault scenarios are used to evaluate the proposed methodology. The resulting mitigation schemes are validated using combined transmission-distribution system simulations, thereby demonstrating that {\mu}PMU measurements along with the RDSM enable FIDVR mitigation by optimal control of reactive power injection from DERs with minimal load disconnection


title: Monitoring Long Term Voltage Instability due to Distribution & Transmission Interaction using Unbalanced {\mu}PMU & PMU Measurements

abstract: This paper extends the idea of the Thevenin equivalent into unbalanced 3{\phi} circuits and proposes a 3{\phi} long-term voltage stability indicator (VSI), that can identify critical loads in a system. Furthermore, in order to identify whether the voltage stability limit is due to the transmission network or a distribution network, a transmission-distribution distinguishing index (TDDI) is proposed. The novelty in the proposed indices is that they can account for the unbalance in the lines and loads, enabling them to use unbalanced phasor measurements naturally. This is supported by mathematical derivations and numerical results. A convex optimization formulation to estimate the 3{\phi} Thevenin equivalent using PMU & {\mu}PMU measurements is proposed, making it possible to calculate VSI and TDDI in an online model-free manner. Numerical simulations performed using co-simulation between Pypower and GridlabD are presented for the IEEE 9 bus and the 30 bus transmission networks combined with several modified IEEE 13 node and 37 node distribution networks. These case studies validate the proposed 3{\phi}-VSI and TDDI over a wide range of scenarios and demonstrate the importance of {\mu}PMU measurements in identifying the regions causing long term voltage instability.


title: Constraints and degrees of freedom in Lorentz-violating field theories

abstract: Many current models which "violate Lorentz symmetry" do so via a vector or tensor field which takes on a vacuum expectation value, thereby spontaneously breaking the underlying Lorentz symmetry of the Lagrangian. To obtain a tensor field with this behavior, one can posit a smooth potential for this field, in which case it would be expected to lie near the minimum of its potential. Alternately, one can enforce a non-zero tensor value via a Lagrange multiplier. The present work explores the relationship between these two types of theories in the case of vector models. In particular, the na\"ive expectation that a Lagrange multiplier "kills off" one degree of freedom via its constraint does not necessarily hold for vector models that already contain primary constraints. It is shown that a Lagrange multiplier can only reduce the degrees of freedom of a model if the field-space function defining the vacuum manifold commutes with the primary constraints.


title: Transient stratification force on particles crossing a density interface

abstract: We perform a series of experiments to measure Lagrangian trajectories of settling and rising particles as they traverse a density interface of thickness $h$ using an index-matched water-salt-ethanol solution. The experiments confirm the substantial deceleration that particles experience as a result of the additional force exerted on the particle due to the sudden change in density. This stratification force is calculated from the measurement data for all particle trajectories. In the absence of suitable parameterizations in the literature, a simple phenomenological model is developed which relies on parameterizations of the effective wake volume and recovery time scale. The model accurately predicts the particle trajectories obtained in our experiments and those of \cite{Fernando1999}. Furthermore, the model demonstrates that the problem depends on four key parameters, namely the entrance Reynolds number $Re_1$, entrance Froude number $Fr$, particle to fluid density ratio $\rho_p/\rho_f$, and relative interface thickness $h/a$.


title: Smart Network Field Theory: The Technophysics of Blockchain and Deep Learning

abstract: The aim of this paper is to propose a theoretical construct, smart network field theory, for the characterization, monitoring, and control of smart network systems. Smart network systems are intelligent autonomously-operating networks, a new form of global computational infrastructure that includes blockchains, deep learning, and autonomous-strike UAVs. These kinds of large-scale networks are a contemporary reality with thousands, millions, and billions of constituent elements, and entail a foundational and theoretically-robust model for their design and operation. Hence this work proposes smart network field theory, drawing from statistical physics, effective field theories, and model systems, for criticality detection and fleet-many item orchestration in smart network systems. Smart network field theory falls within the broader concern of technophysics (the application of physics to the study of technology), in which a key objective is deriving standardized methods for assessing system criticality and phase transition, and defining interim system structure between the levels of microscopic noise and macroscopic labels. The farther implications of this work include the possibility of recasting the P/NP computational complexity schema as one no longer based on traditional time (concurrency) and space constraints, due to the availability of smart network computational resources.


title: Post-earthquake modelling of transportation networks using an agent-based model

abstract: Surface transportation systems are an essential part of urban transportation infrastructure and are susceptible to damage from earthquakes. This damage, along with the lack of prior warning of earthquake events, may lead to severe and unexpected disruption of normal traffic patterns, which may seriously impair post-disaster response. Accordingly, it is important to understand and model the performance of urban transportation systems immediately following an earthquake, to evaluate its capability to support emergency response, e.g., the movement of firefighters, search and rescue teams and medical personnel, and the transportation of injured people to emergency treatment facilities. For this purpose, a scenario-based methodology is developed to model the performance of a transportation network immediately following an earthquake using an agent-based model. This methodology accounts for the abrupt changes in destination, irrational behavior of drivers, unavailability of traffic information and the impairment to traffic capacity due to bridge damage and building debris. An illustration using the road network of Tangshan City, China shows that the method can capture the traffic flow characteristics immediately after an earthquake and can determine the capability of the transportation network to transfer injured people to hospitals considering the above factors. Thus, it can provide rational support for evaluating the performance of the surface transportation system under immediate post-disaster emergency conditions.


title: Universal origin of boson peak vibrational anomalies in ordered crystals and in amorphous materials

abstract: The vibrational spectra of solids, both ordered and amorphous, in the low-energy regime, control the thermal and transport properties of materials, from heat capacity to heat conduction, electron-phonon couplings, conventional superconductivity etc. The old Debye model of vibrational spectra at low energy gives the vibrational density of states (VDOS) as proportional to the frequency squared, but in many materials the spectrum departs from this law which results in a peak upon normalizing the VDOS by frequency squared, which is known as the "boson peak". A description of the VDOS of solids (both crystals and glasses) is presented starting from first principles. Without using any assumptions whatsoever about the existence and nature of "disorder" in the material, it is shown that the boson peak in the VDOS of both ordered crystals and glasses arises naturally from the competition between elastic mode propagation and viscous damping. The theory explains the recent experimental observations of boson peak in perfectly ordered crystals, which cannot be explained based on previous theoretical frameworks. The theory also explains, for the first time, how the vibrational spectrum changes with the atomic density of the solid, and explains recent experimental observations of this effect.

categories: cond-mat.soft cond-mat.dis-nn cond-mat.stat-mech hep-th


title: A High Accuracy and High Sensitivity System Architecture for Electrical Impedance Tomography System

abstract: A high accuracy and high sensitivity system architecture is proposed for the read-out circuit of electrical impedance tomography system-on-chip. The switched ratiometric technique is applied in the proposed architecture. The proposed system architecture minimizes the device noise by processing signals from both read-out electrodes and the stimulus. The quantized signals are post-processed in the digital processing unit for proper signal demodulation and impedance ratio calculation. Our proposed architecture improves the sensitivity of the read-out circuit, cancels out the gain fluctuations in the system, and overcomes the effects of motion artifacts on measurements.