title: Same-different problems strain convolutional neural networks

abstract: The robust and efficient recognition of visual relations in images is a hallmark of biological vision. We argue that, despite recent progress in visual recognition, modern machine vision algorithms are severely limited in their ability to learn visual relations. Through controlled experiments, we demonstrate that visual-relation problems strain convolutional neural networks (CNNs). The networks eventually break altogether when rote memorization becomes impossible, as when intra-class variability exceeds network capacity. Motivated by the comparable success of biological vision, we argue that feedback mechanisms including attention and perceptual grouping may be the key computational components underlying abstract visual reasoning.\


title: Uneven flows: On cosmic bulk flows, local observers, and gravity

abstract: Using N-body simulations we study the impact of various systematic effects on the bulk flow (BF) and the Cosmic Mach Number (CMN). We consider two types of systematics: those related to survey properties and those induced by observer's location in the Universe. In the former category we model sparse sampling, velocity errors, and survey incompleteness. In the latter, we consider Local Group (LG) analogue observers, placed in a specific location within the Cosmic Web, satisfying various observational criteria. We differentiate such LG observers from Copernican ones, who are at random locations. We report strong systematic effects on the measured BF and CMN induced by sparse sampling, velocity errors and radial incompleteness. For BF most of these effects exceed 10\% for scales $R\leq100h^{-1}$Mpc. For CMN some of these systematics can be catastrophically large ($>50\%$) also on bigger scales. Moreover, we find that the position of the observer in the Cosmic Web significantly affects the locally measured BF (CMN), with effects as large as $\sim20\%$ ($30\%)$ at $R\leq50h^{-1}$Mpc for a LG-like observer as compared to a random one. This effect is comparable to the sample variance. To highlight the importance of these systematics, we additionally study a model of modified gravity (MG) with $\sim15\%$ enhanced growth rate. We found that the systematic effects can mimic the modified gravity signal. The worst-case scenario is realized for a case of a LG-like observer, when the effects induced by local structures are degenerate with the enhanced growth rate fostered by MG. Our results indicate that dedicated constrained simulations and realistic mock galaxy catalogs will be absolutely necessary to fully benefit from the statistical power of the forthcoming peculiar velocity data from surveys such as TAIPAN, WALLABY, Cosmic Flows-4 and SKA.


title: On domination perfect graphs

abstract: Let $\gamma(G)$ and $\beta(G)$ denote the domination number and the covering number of a graph $G$, respectively. A connected non-trivial graph $G$ is said to be $\gamma\beta$-{perfect} if $\gamma(H)=\beta(H)$ for every non-trivial induced connected subgraph $H$ of $G$. In this note we present an elementary proof of a characterization of the $\gamma\beta$-perfect graphs.


title: A Study of WhatsApp Usage Patterns and Prediction Models without Message Content

abstract: Internet social networks have become a ubiquitous application allowing people to easily share text, pictures, and audio and video files. Popular networks include WhatsApp, Facebook, Reddit and LinkedIn. We present an extensive study of the usage of the WhatsApp social network, an Internet messaging application that is quickly replacing SMS messaging. In order to better understand people's use of the network, we provide an analysis of over 6 million messages from over 100 users, with the objective of building demographic prediction models using activity data. We performed extensive statistical and numerical analysis of the data and found significant differences in WhatsApp usage across people of different genders and ages. We also inputted the data into the Weka data mining package and studied models created from decision tree and Bayesian network algorithms. We found that different genders and age demographics had significantly different usage habits in almost all message and group attributes. We also noted differences in users' group behavior and created prediction models, including the likelihood a given group would have relatively more file attachments, if a group would contain a larger number of participants, a higher frequency of activity, quicker response times and shorter messages. We were successful in quantifying and predicting a user's gender and age demographic. Similarly, we were able to predict different types of group usage. All models were built without analyzing message content. We present a detailed discussion about the specific attributes that were contained in all predictive models and suggest possible applications based on these results.


title: Mimetic gravity as DHOST theories

abstract: We show that theories of mimetic gravity can be viewed as degenerate higher-order scalar-tensor (DHOST) theories that admit an extra local (gauge) symmetry in addition to the usual diffeomorphism invariance. We reformulate and classify mimetic theories in this perspective. Using the effective theory of dark energy, recently extended to include DHOST theories, we then investigate the linear perturbations about a homogeneous and isotropic background for all mimetic theories. We also include matter, in the form of a $k$-essence scalar field, and we derive the quadratic action for linear perturbations in this case.


title: Bootstrap validation of links of a minimum spanning tree

abstract: We describe two different bootstrap methods applied to the detection of a minimum spanning tree obtained from a set of multivariate variables. We show that two different bootstrap procedures provide partly distinct information that can be highly informative about the investigated complex system. Our case study, based on the investigation of daily returns of a portfolio of stocks traded in the US equity markets, shows the degree of robustness and completeness of the information extracted with popular information filtering methods such as the minimum spanning tree and the planar maximally filtered graph. The first method performs a "row bootstrap" whereas the second method performs a "pair bootstrap". We show that the parallel use of the two methods is suggested especially for complex systems presenting both a nested hierarchical organization together with the presence of global feedback channels.


title: Predicting Customer Churn: Extreme Gradient Boosting with Temporal Data

abstract: Accurately predicting customer churn using large scale time-series data is a common problem facing many business domains. The creation of model features across various time windows for training and testing can be particularly challenging due to temporal issues common to time-series data. In this paper, we will explore the application of extreme gradient boosting (XGBoost) on a customer dataset with a wide-variety of temporal features in order to create a highly-accurate customer churn model. In particular, we describe an effective method for handling temporally sensitive feature engineering. The proposed model was submitted in the WSDM Cup 2018 Churn Challenge and achieved first-place out of 575 teams.


title: Network organization of coopetitive genetic influences on cortical morphologies

abstract: Brain can be represented as a network, where regions are the nodes and relations between the regions are edges. Within a network, co-existence of cooperative and competitive relationships between different nodes is called coopetition. Inter-regional genetic influences on morphological phenotypes (cortical thickness, surface area) of cortex display such coopetitive relationships. Here, we have represented these genetic influences as a network and shown that cooperative and competitive genetic influences on cortical morphological phenotypes follow distinct organization principles. Utilizing the theory of structural balance, we have shown that the pattern of collective regulation of cortical morphological phenotypes by cooperative and competitive genetic influences are overall bilaterally symmetric and such patterns of collective genetic regulation are similar to the principal modes of population variation of cortical morphological phenotypes. Finally, we have observed that the maximally and minimally imbalanced regions corresponding to the collective genetic regulation partially overlap with the cortical structural network hubs.


title: Precision determination of the strong coupling constant within a global PDF analysis

abstract: We present a determination of the strong coupling constant $\alpha_s(m_Z)$ based on the NNPDF3.1 determination of parton distributions, which for the first time includes constraints from jet production, top-quark pair differential distributions, and the $Z$ $p_T$ distributions using exact NNLO theory. Our result is based on a novel extension of the NNPDF methodology - the correlated replica method - which allows for a simultaneous determination of $\alpha_s$ and the PDFs with all correlations between them fully taken into account. We study in detail all relevant sources of experimental, methodological and theoretical uncertainty. At NNLO we find $\alpha_s(m_Z) = 0.1185 \pm 0.0005^\text{(exp)}\pm 0.0001^\text{(meth)}$, showing that methodological uncertainties are negligible. We conservatively estimate the theoretical uncertainty due to missing higher order QCD corrections (N$^3$LO and beyond) from half the shift between the NLO and NNLO $\alpha_s$ values, finding $\Delta\alpha^{\rm th}_s =0.0011$.


title: DarkSUSY 6 : An Advanced Tool to Compute Dark Matter Properties Numerically

abstract: The nature of dark matter remains one of the key science questions. Weakly Interacting Massive Particles (WIMPs) are among the best motivated particle physics candidates, allowing to explain the measured dark matter density by employing standard big-bang thermodynamics. We introduce here a radically new version of the widely used DarkSUSY package, which allows to compute the properties of such dark matter particles numerically. With DarkSUSY 6 one can accurately predict a large variety of astrophysical signals from dark matter, such as direct detection rates in low-background counting experiments and indirect detection signals through antiprotons, antideuterons, gamma rays and positrons from the Galactic halo, or high-energy neutrinos from the center of the Earth or of the Sun. For thermally produced dark matter like WIMPs, high-precision tools are provided for the computation of the relic density in the Universe today, as well as for the size of the smallest dark matter protohalos. Furthermore, the code allows to calculate dark matter self-interaction rates, which may affect the distribution of dark matter at small cosmological scales. Compared to earlier versions, DarkSUSY 6 introduces many significant physics improvements and extensions. The most fundamental new feature of this release, however, is that the code has been completely re-organized and brought into a highly modular and flexible shape. Switching between different pre-implemented dark matter candidates has thus become straight-forward, just as adding new - WIMP or non-WIMP - particle models or replacing any given functionality in a fully user-specified way. In this article, we describe the physics behind the computer package, along with the main structure and philosophy of this major revision of DarkSUSY. A detailed manual is provided together with the public release at www.darksusy.org.


title: Extended and broad Lyman alpha emission around a BAL quasar at z~5

abstract: In this work we report deep MUSE observations of a Broad Absorption Line (BAL) quasar at z ~ 5, revealing a Lyman alpha nebula with a maximum projected linear size of ~ 60 kpc around the quasar (down to our 2-sigma SB limit per layer of ~ 9e-19 erg/s/cm^2/arcsec^2 for a 1 arcsec^2 aperture). After correcting for the cosmological surface brightness dimming, we find that our nebula, at z ~ 5, has an intrinsically less extended Lyman alpha emission than nebulae at lower redshift. However, such a discrepancy is greatly reduced when referring to comoving distances, which take into account the cosmological growth of dark matter (DM) haloes, suggesting a positive correlation between the size of Lyman alpha nebulae and the sizes of DM haloes/structures around quasars. Differently from the typical nebulae around radio-quiet non-BAL quasars, in the inner regions (~ 10 kpc) of the circumgalactic medium (CGM) of our source, the velocity dispersion of the Lyman alpha emission is very high (FWHM > 1000 km/s), suggesting that in our case we may be probing outflowing material associated with the quasar.


title: Anomalous thermalization and transport in disordered interacting Floquet systems

abstract: Local observables in generic periodically driven closed quantum systems are known to relax to values described by periodic infinite temperature ensembles. At the same time, ergodic static systems exhibit anomalous thermalization of local observables and satisfy a modified version of the eigenstate thermalization hypothesis (ETH), when disorder is present. This raises the question, how does the introduction of disorder affect relaxation in periodically driven systems? In this work, we analyze this problem by numerically studying transport and thermalization in an archetypal example. We find that thermalization is anomalous and is accompanied by subdiffusive transport with a disorder dependent dynamical exponent. Distributions of matrix elements of local operators in the eigenbases of a family of effective time-independent Hamiltonians, which describe the stroboscopic dynamics of such systems, show anomalous departures from predictions of ETH signaling that only a modified version of ETH is satisfied. The dynamical exponent is shown to be related to the scaling of the variance of these distributions with system size.


title: Bow shock in merging cluster A520: the edge of the radio halo and the electron-ion equilibration timescale

abstract: We studied the prominent bow shock in the merging galaxy cluster A520 using a deep Chandra X-ray observation and archival VLA radio data. This shock is a useful diagnostic tool, owing to its clear geometry and relatively high Mach number. At the "nose" of the shock, we measure a Mach number of $M=2.4_{-0.2}^{+0.4}$. The shock becomes oblique away from the merger axis, with the Mach number falling to $\simeq$1.6 around 30$^{\circ}$ from the nose. The electron temperature immediately behind the shock nose is consistent with that from the Rankine-Hugoniot adiabat, and is higher (at a 95% confidence) than expected for adiabatic compression of electrons followed by Coulomb electron-proton equilibration, indicating the presence of equilibration mechanisms faster than Coulomb collisions. This is similar to an earlier finding for the Bullet cluster. We also combined four archival VLA datasets to obtain a better image of the cluster's giant radio halo at 1.4 GHz. An abrupt edge of the radio halo traces the shock front, and no emission is detected in the pre-shock region. If the radio edge were due only to adiabatic compression of relativistic electrons in pre-shock plasma, we would expect a pre-shock radio emission detectable in this radio dataset; however, an interferometric artifact dominates the uncertainty, so we cannot rule this model out. Other interesting features of the radio halo include a peak at the remnant of the cool core, suggesting that the core used to have a radio minihalo, and a peak marking a possible region of high turbulence.


title: Gauge-Invariant Formalism of Cosmological Weak Lensing

abstract: We present the gauge-invariant formalism of cosmological weak lensing, accounting for all the relativistic effects due to the scalar, vector, and tensor perturbations at the linear order. While the light propagation is fully described by the geodesic equation, the relation of the photon wavevector to the physical quantities requires the specification of the frames, where they are defined. By constructing the local tetrad bases at the observer and the source positions, we clarify the relation of the weak lensing observables such as the convergence, the shear, and the rotation to the physical size and shape defined in the source rest-frame and the observed angle and redshift measured in the observer rest-frame. Compared to the standard lensing formalism, additional relativistic effects contribute to all the lensing observables. We explicitly verify the gauge-invariance of the lensing observables and compare our results to previous work. In particular, we demonstrate that even in the presence of the vector and tensor perturbations, the physical rotation of the lensing observables vanishes at the linear order, while the tetrad basis rotates along the light propagation compared to a FRW coordinate. Though the latter is often used as a probe of primordial gravitational waves, the rotation of the tetrad basis is indeed not a physical observable. We further clarify its relation to the E-B decomposition in weak lensing. Our formalism provides a transparent and comprehensive perspective of cosmological weak lensing.


title: Prospects for resolving the Hubble constant tension with standard sirens

abstract: The Hubble constant ($H_0$) estimated from the local Cepheid-supernova (SN) distance ladder is in 3-$\sigma$ tension with the value extrapolated from cosmic microwave background (CMB) data assuming the standard cosmological model. Whether this tension represents new physics or systematic effects is the subject of intense debate. Here, we investigate how new, independent $H_0$ estimates can arbitrate this tension, assessing whether the measurements are consistent with being derived from the same model using the posterior predictive distribution (PPD). We show that, with existing data, the inverse distance ladder formed from BOSS baryon acoustic oscillation measurements and the Pantheon SN sample yields an $H_0$ posterior near-identical to the Planck CMB measurement. The observed local distance ladder value is a very unlikely draw from the resulting PPD. Turning to the future, we find that a sample of $\sim50$ binary neutron star "standard sirens" (detectable within the next decade) will be able to adjudicate between the local and CMB estimates.


title: Particle-without-Particle: a practical pseudospectral collocation method for linear partial differential equations with distributional sources

abstract: Partial differential equations with distributional sources---in particular, involving (derivatives of) delta distributions---have become increasingly ubiquitous in numerous areas of physics and applied mathematics. It is often of considerable interest to obtain numerical solutions for such equations, but any singular ("particle"-like) source modeling invariably introduces nontrivial computational obstacles. A common method to circumvent these is through some form of delta function approximation procedure on the computational grid; however, this often carries significant limitations on the efficiency of the numerical convergence rates, or sometimes even the resolvability of the problem at all. In this paper, we present an alternative technique for tackling such equations which avoids the singular behavior entirely: the "Particle-without-Particle" method. Previously introduced in the context of the self-force problem in gravitational physics, the idea is to discretize the computational domain into two (or more) disjoint pseudospectral (Chebyshev-Lobatto) grids such that the "particle" is always at the interface between them; thus, one only needs to solve homogeneous equations in each domain, with the source effectively replaced by jump (boundary) conditions thereon. We prove here that this method yields solutions to any linear PDE the source of which is any linear combination of delta distributions and derivatives thereof supported on a one-dimensional subspace of the problem domain. We then implement it to numerically solve a variety of relevant PDEs: hyperbolic (with applications to neuroscience and acoustics), parabolic (with applications to finance), and elliptic. We generically obtain improved convergence rates relative to typical past implementations relying on delta function approximations.

categories: physics.comp-ph gr-qc math.NA q-bio.NC q-fin.CP


title: LOFAR 150-MHz observations of SS 433 and W 50

abstract: We present LOFAR high-band data over the frequency range 115-189 MHz for the X-ray binary SS 433, obtained in an observing campaign from 2013 February - 2014 May. Our results include a deep, wide-field map, allowing a detailed view of the surrounding supernova remnant W 50 at low radio frequencies, as well as a light curve for SS 433 determined from shorter monitoring runs. The complex morphology of W 50 is in excellent agreement with previously published higher-frequency maps; we find additional evidence for a spectral turnover in the eastern wing, potentially due to foreground free-free absorption. Furthermore, SS 433 is tentatively variable at 150 MHz, with both a debiased modulation index of 11 per cent and a $\chi^2$ probability of a flat light curve of $8.2 \times 10^{-3}$. By comparing the LOFAR flux densities with contemporaneous observations carried out at 4800 MHz with the RATAN-600 telescope, we suggest that an observed $\sim$0.5-1 Jy rise in the 150-MHz flux density may correspond to sustained flaring activity over a period of approximately six months at 4800 MHz. However, the increase is too large to be explained with a standard synchrotron bubble model. We also detect a wealth of structure along the nearby Galactic plane, including the most complete detection to date of the radio shell of the candidate supernova remnant G 38.7-1.4. This further demonstrates the potential of supernova remnant studies with the current generation of low-frequency radio telescopes.


title: Supergravity in the group-geometric framework: a primer

abstract: We review the group-geometric approach to supergravity theories, in the perspective of recent developments and applications. Usual diffeomorphisms, gauge symmetries and supersymmetries are unified as superdiffeomorphisms in a supergroup manifold. Integration on supermanifolds is briefly revisited, and used as a tool to provide a bridge between component and superspace actions. As an illustration of the constructive techniques, the cases of $d=3,4$ off-shell supergravities and $d=5$ Chern-Simons supergravity are discussed in detail. A cursory account of $d=10+2$ supergravity is also included. We recall a covariant canonical formalism, well adapted to theories described by Lagrangians $d$-forms, that allows to define a form hamiltonian and to recast constrained hamiltonian systems in a covariant form language. Finally, group geometry and properties of spinors and gamma matrices in $d=s+t$ dimensions are summarized in Appendices.


title: On the Computational Complexity of Curing the Sign Problem

abstract: Quantum many-body systems whose Hamiltonians are non-stoquastic, i.e., have positive off-diagonal matrix elements in a given basis, are known to pose severe limitations on the efficiency of Quantum Monte Carlo algorithms designed to simulate them, due to the infamous sign problem. We study the computational complexity associated with `curing' non-stoquastic Hamiltonians, i.e., transforming them into sign-problem-free ones. We prove that if such transformations are limited to single-qubit Clifford group elements or general single-qubit orthogonal matrices, finding the curing transformation is NP-complete. We discuss the implications of this result.


title: Deflections of ultra-high energy cosmic rays by the Milky Way magnetic field: how well can they be corrected?

abstract: Locating the sources of ultra-high energy cosmic rays (UHECRs) still remains a difficult puzzle for modern astrophysics. A major hurdle in the search for the sources is the fact that UHECRs are deflected by the Galactic magnetic field (GMF). Current knowledge of the GMF is limited, as most experimental measurements track line-of-sight--integrated quantities that are used to obtain best-fit parameters for global models including a large random component. The advent, however, of the Gaia era, with measurements of $\sim 10^9$ stellar parallaxes, in combination with upcoming large polarimetric surveys, make, for the first time, a 3D measurement of the GMF possible in principle. Such measurements can then be used to attempt a reconstruction of the trajectories of individual UHECRs through the Galaxy, in order to correct for their deflection. Motivated by these developments, in the present work, we study the limits of such a correction, by examining how its effectiveness depends on the uncertainty of any such future magnetic field measurements. To that end, we simulate attempts to reconstruct the trajectory of the cosmic ray by using hypothetical measurements of the GMF, based on values received from two recently updated GMF models. To simulate the uncertainty of a 3D measurement, random errors to these values are introduced separately for the plane-of-the-sky (POS) magnitude, the line-of-sight (LOS) magnitude and the POS direction. Our results highlight the conditions under which an effective correction is achievable. We find that the effectiveness is dependent on the particle rigidity and arrival direction, and can vary significantly depending on the GMF model used.


title: Generalized Eigenvectors of Isospectral Transformations,Spectral Equivalence and Reconstruction of Original Networks

abstract: Isospectral transformations (IT) of matrices and networks allow for compression of either object while keeping all the information about their eigenvalues and eigenvectors.We analyze here what happens to generalized eigenvectors under isospectral transformations and to what extent the initial network can be reconstructed from its compressed image under IT. We also generalize and essentially simplify the proof that eigenvectors are invariant under isospectral transformations and generalize and clarify the notion of spectral equivalence of networks.


title: Charting generalized supersoft supersymmetry

abstract: Without any shred of evidence for new physics from LHC, the last hiding spots of natural electroweak supersymmetry seem to lie either in compressed spectra or in spectra where scalars are suppressed with respect to the gauginos. While in the MSSM (or in any theory where supersymmetry is broken by the $F$-vev of a chiral spurion), a hierarchy between scalar and gaugino masses requires special constructions, it is automatic in scenarios where supersymmetry is broken by $D$-vev of a real spurion. In the latter framework, gaugino mediated contributions to scalar soft masses are finite (loop suppressed but not $\log$-enhanced), a feature often referred to as "supersoftness". Though phenomenologically attractive, pure supersoft models suffer from the $\mu$-problem, potential color-breaking minima, large $T$-parameter, etc. These problems can be overcome without sacrificing the model's virtues by departing from pure supersoftness and including $\mu$-type operators that use the same $D$-vev, a framework known as generalized supersoft supersymmetry. The main purpose of this paper is to point out that the new operators also solve the last remaining issue associated with supersoft spectra, namely that a right handed (RH) slepton is predicted to be the lightest superpartner, rendering the setup cosmologically unfeasible. In particular, we show that the $\mu$-operators in generalized supersoft generate a new source for scalar masses, which can raise the RH-slepton mass above bino due to corrections from renormalisation group evolutions (RGEs). In fact, a mild tuning can open up the bino--RH slepton coannihilation regime for a thermal dark matter. We derive the full set of RGEs required to determine the spectrum at low energies. Beginning with input conditions at a high scale, we show that completely viable spectra can be achieved.


title: Multi-mode four-wave mixing with a spatially-structured pump

abstract: We demonstrate a new four-wave mixing (4WM) geometry based on structured light. By utilizing near-field diffraction through a narrow slit, the pump beam is asymmetrically structured to modify the phase matching condition, generating multi-mode output in both the spatial and frequency domains. We show that the frequency parameter enables selection of various spatial-mode outputs, including a twin-beam geometry which preserves relative intensity squeezing shared between the two beams. The results suggest that the engineering of atomic states via structured light may provide a pathway to a diverse set of quantum resources based on multi-mode squeezed light.


title: Small zeros of Dirichlet $L$-functions of quadratic characters of prime modulus

abstract: In this paper, we investigate the distribution of the imaginary parts of zeros near the real axis of Dirichlet $L$-functions associated to the quadratic characters $\chi_{p}(\cdot)=(\cdot |p)$ with $p$ a prime number. Assuming the Generalized Riemann Hypothesis (GRH), we compute the one-level density for the zeros of this family of $L$-functions under the condition that the Fourier transform of the test function is supported on a closed subinterval of $(-1,1)$. We also write down the ratios conjecture for this family of $L$-functions a la Conrey, Farmer and Zirnbauer and derive a conjecture for the one-level density which is consistent with the main theorem of this paper and with the Katz-Sarnak prediction and includes lower order terms. Following the methods of \"Ozl\"uk and Snyder, we prove that GRH implies $L(\frac{1}{2},\chi_p)\neq 0$ for at least $75\%$ of the primes.


title: Co-formation of the Galactic disc and the stellar halo

abstract: Using a large sample of Main Sequence stars with 7-D measurements supplied by Gaia and SDSS, we study the kinematic properties of the local (within ~10 kpc from the Sun) stellar halo. We demonstrate that the halo's velocity ellipsoid evolves strongly with metallicity. At the low [Fe/H] end, the orbital anisotropy (the amount of motion in the radial direction compared to the tangential one) is mildly radial with 0.2<beta<0.4. However, for stars with [Fe/H]>-1.7 we measure extreme values of beta~0.9. Across the metallicity range considered, i.e. -3<[Fe/H]<-1, the stellar halo's spin is minimal, at the level of 20<v_theta (km/s) <30. Using a suite of cosmological zoom-in simulations of halo formation, we deduce that the observed acute anisotropy is inconsistent with the continuous accretion of dwarf satellites. Instead, we argue, the stellar debris in the inner halo were deposited in a major accretion event by a satellite with Mvir>10^10 Msun around the epoch of the Galactic disc formation, i.e. between 8 and 11 Gyr ago. The radical halo anisotropy is the result of the dramatic radialisation of the massive progenitor's orbit, amplified by the action of the growing disc.


title: Colored ray configurations

abstract: We study the cyclic color sequences induced at infinity by colored rays with apices being a given balanced finite bichromatic point set. We first study the case in which the rays are required to be pairwise disjoint. We derive a lower bound on the number of color sequences that can be realized from any such fixed point set and examine color sequences that can be realized regardless of the point set, exhibiting negative examples as well. We also provide a tight upper bound on the number of configurations that can be realized from a point set, and point sets for which there are asymptotically less configurations than that number. In addition, we provide algorithms to decide whether a color sequence is realizable from a given point set in a line or in general position. We address afterwards the variant of the problem where the rays are allowed to intersect. We prove that for some configurations and point sets, the number of ray crossings must be $\Theta(n^2)$ and study then configurations that can be realized by rays that pairwise cross. We show that there are point sets for which the number of configurations that can be realized by pairwise-crossing rays is asymptotically smaller than the number of configurations realizable by pairwise-disjoint rays. We provide also point sets from which any configuration can be realized by pairwise-crossing rays and show that there is no configuration that can be realized by pairwise-crossing rays from every point set.


title: Global stability of a distributed delayed viral model with general incidence rate

abstract: In this paper, we discussed an infinitely distributed delayed viral infection model with nonlinear immune response and general incidence rate. We proved the existence and uniqueness of the equilibria. By using the Lyapunov functional and LaSalle invariance principle, we obtained the conditions of global stabilities of the infection-free equilibrium, the immune-exhausted equilibrium and the endemic equilibrium. Numerical simulations are given to verify the analytical results.


title: Narrow Artificial Intelligence with Machine Learning for Real-Time Estimation of a Mobile Agents Location Using Hidden Markov Models

abstract: We propose to use a supervised machine learning technique to track the location of a mobile agent in real time. Hidden Markov Models are used to build artificial intelligence that estimates the unknown position of a mobile target moving in a defined environment. This narrow artificial intelligence performs two distinct tasks. First, it provides real-time estimation of the mobile agent's position using the forward algorithm. Second, it uses the Baum-Welch algorithm as a statistical learning tool to gain knowledge of the mobile target. Finally, an experimental environment is proposed, namely a video game that we use to test our artificial intelligence. We present statistical and graphical results to illustrate the efficiency of our method.


title: Predicting University Students' Academic Success and Major using Random Forests

abstract: In this article, a large data set containing every course taken by every undergraduate student in a major university in Canada over 10 years is analysed. Modern machine learning algorithms can use large data sets to build useful tools for the data provider, in this case, the university. In this article, two classifiers are constructed using random forests. To begin, the first two semesters of courses completed by a student are used to predict if they will obtain an undergraduate degree. Secondly, for the students that completed a program, their major is predicted using once again the first few courses they have registered to. A classification tree is an intuitive and powerful classifier and building a random forest of trees improves this classifier. Random forests also allow for reliable variable importance measurements. These measures explain what variables are useful to the classifiers and can be used to better understand what is statistically related to the students' situation. The results are two accurate classifiers and a variable importance analysis that provides useful information to university administrations.


title: A Correction for IUE UV Flux Distributions from Comparisons with CALSPEC

abstract: A collection of spectral energy distributions (SEDs) is available in the Hubble Space Telescope (HST) CALSPEC database that is based on calculated model atmospheres for pure hydrogen white dwarfs (WDs). A much larger set (~100,000) of UV SEDs covering the range (1150-3350A) with somewhat lower quality are available in the IUE database. IUE low-dispersion flux distributions are compared with CALSPEC to provide a correction that places IUE fluxes on the CALSPEC scale. While IUE observations are repeatable to only 4-10% in regions of good sensitivity, the average flux corrections have a precision of 2-3%. Our re-calibration places the IUE flux scale on the current UV reference standard and is relevant for any project based on IUE archival data, including our planned comparison of GALEX to the corrected IUE fluxes. IUE SEDs may be used to plan observations and cross-calibrate data from future missions, so the IUE flux calibration must be consistent with HST instrumental calibrations to the best possible precision.


abstract: Here we propose the Reweighted Autoencoded Variational Bayes for Enhanced Sampling (RAVE) method, a new iterative scheme that uses the deep learning framework of variational autoencoders to enhance sampling in molecular simulations. RAVE involves iterations between molecular simulations and deep learning in order to produce an increasingly accurate probability distribution along a low-dimensional latent space that captures the key features of the molecular simulation trajectory. Using the Kullback-Leibler divergence between this latent space distribution and the distribution of various trial reaction coordinates sampled from the molecular simulation, RAVE determines an optimum, yet nonetheless physically interpretable, reaction coordinate and optimum probability distribution. Both then directly serve as the biasing protocol for a new biased simulation, which is once again fed into the deep learning module with appropriate weights accounting for the bias, the procedure continuing until estimates of desirable thermodynamic observables are converged. Unlike recent methods using deep learning for enhanced sampling purposes, RAVE stands out in that (a) it naturally produces a physically interpretable reaction coordinate, (b) is independent of existing enhanced sampling protocols to enhance the fluctuations along the latent space identified via deep learning, and (c) it provides the ability to easily filter out spurious solutions learned by the deep learning procedure. The usefulness and reliability of RAVE is demonstrated by applying it to model potentials of increasing complexity, including computation of the binding free energy profile for a hydrophobic ligand-substrate system in explicit water with dissociation time of more than three minutes, in computer time at least twenty times less than that needed for umbrella sampling or metadynamics.

categories: physics.chem-ph cond-mat.soft cond-mat.stat-mech physics.comp-ph


title: Measurement of the Decays $\Lambda_c\to \Sigma\pi\pi$ at Belle

abstract: We report measurements of the branching fractions of the decays $\Lambda^+_c\to\Sigma^+\pi^-\pi^+$, $\Lambda^+_c\to\Sigma^0\pi^+\pi^0$ and $\Lambda^+_c\to\Sigma^+\pi^0\pi^0$ relative to the reference channel $\Lambda^+_c\to pK^-\pi^+$. The analysis is based on the full data sample collected at and close to $\Upsilon(4S)$ resonance by the Belle detector at the KEKB asymmetric-energy $e^+e^-$ collider corresponding to the integrated luminosity of 711 fb$^{-1}$. We measure ${\cal B}(\Lambda^+_c\rightarrow\Sigma^+\pi^-\pi^+)/{\cal B}(\Lambda^+_c\to pK^-\pi^+) = 0.706~\pm 0.003~\pm 0.029$, ${\cal B}(\Lambda^+_c\rightarrow\Sigma^0\pi^+\pi^0)/{\cal B}(\Lambda^+_c\to pK^-\pi^+) = 0.491~\pm 0.005~\pm 0.028$ and ${\cal B}(\Lambda^+_c\rightarrow\Sigma^+\pi^0\pi^0)/{\cal B}(\Lambda^+_c\to pK^-\pi^+) = 0.198~\pm 0.006~\pm 0.016$. The listed uncertainties are statistical and systematic, respectively.


title: State of the Practice for GIS Software

abstract: We present a reproducible method to analyze the state of software development practices in a given scientific domain and apply this method to Geographic Information Systems (GIS). The analysis is based on grading a set of 30 GIS products using a template of 56 questions based on 13 software qualities. The products range in scope and purpose from a complete desktop GIS systems, to stand-alone tools, to programming libraries/packages. The final ranking of the products is determined using the Analytic Hierarchy Process (AHP), a multicriteria decision making method that focuses on relative comparisons between products, rather than directly measuring qualities. The results reveal concerns regarding the correctness, maintainability, transparency and reproducibility of some GIS software. Three recommendations are presented as feedback to the GIS community: 1) Ensure each project has a requirements specification document, 2) Provide a wealth of support methods, such as an IRC (Internet Relay Chat) channel, a Stack Exchange tag for new questions, or opening the issue tracker for support requests, as well as the more traditional email-based methods, and, 3) Design product websites for maximum transparency (of the development process), for open source projects, provide a developer's guide.


title: The WISSH Quasars Project IV. BLR versus kpc-scale winds

abstract: We have undertaken a multi-band observing program aimed at obtaining a complete census of winds in a sample of WISE/SDSS selected hyper-luminous (WISSH) QSOs at z~2-4. We have analyzed the rest-frame optical (LBT/LUCI and VLT/SINFONI) and UV (SDSS) spectra of 18 randomly selected WISSH QSOs to measure the SMBH mass and study the properties of winds both in the NLR and BLR traced by blueshifted/skewed [OIII] and CIV emission lines, respectively. These WISSH QSOs are powered by SMBH with masses $\ge$10$^9$ Msun accreting at 0.4<$\lambda_{Edd}$<3.1. We have found the existence of two sub-populations characterized by the presence of outflows at different distances from the SMBH. One population ([OIII] sources) exhibits powerful [OIII] outflows, rest-frame EW (REW) of the CIV emission REW$_{CIV}\approx$20-40 A and modest CIV velocity shift (v$_{CIV}^{peak}$) with respect to the systemic redshift (<=2000 km/s). The second population (Weak [OIII] sources), representing ~70% of the analyzed WISSH QSOs, shows weak/absent [OIII] emission and an extremely large v$_{CIV}^{peak}$ (up to ~8000 km/s and REW$_{CIV}$<=20 A). We propose two explanations for the observed behavior of the strength of the [OIII] emission in terms of orientation effects of the line of sight and ionization cone. The dichotomy in the presence of BLR and NLR winds could be likely due to inclination effects considering a polar geometry scenario for the BLR winds. We find a strong correlation with L$_{Bol}$ and an anti-correlation with $\alpha_{ox}$, whereby the higher L$_{Bol}$, the steeper $\alpha_{ox}$ and the larger is the v$_{CIV}^{peak}$. Finally, the observed dependence v$_{CIV}^{peak}\propto L_{Bol}^{0.28\pm0.04}$ is consistent with radiatively driven winds scenario, where strong UV continuum is necessary to launch the wind and a weakness of the X-ray emission is fundamental to prevent overionization of the wind itself.


title: Cooling the Motion of a Silica Microsphere in a Magneto-Gravitational Trap in Ultra-High Vacuum

abstract: Levitated optomechanical systems, and particularly particles trapped in vacuum, provide unique platforms for studying the mechanical behavior of objects well-isolated from their environment. Ultimately, such systems may enable the study of fundamental questions in quantum mechanics, gravity, and other weak forces. While the optical trapping of nanoparticles has emerged as the prototypical levitated optomechanical system, it is not without problems due to the heating from the high optical intensity required, particularly when combined with a high vacuum environment. Here we investigate a magneto-gravitational trap in ultra-high vacuum. In contrast to optical trapping, we create an entirely passive trap for diamagnetic particles by utilizing the magnetic field generated by permanent magnets and the gravitational interaction. We demonstrate cooling the center of mass motion of a trapped silica microsphere from ambient temperature to an effective temperature near or below one milliKelvin in two degrees of freedom by optical feedback damping.


title: Asymptotic nonequivalence of density estimation and Gaussian white noise for small densities

abstract: It is well-known that density estimation on the unit interval is asymptotically equivalent to a Gaussian white noise experiment, provided the densities are sufficiently smooth and uniformly bounded away from zero. We show that a uniform lower bound, whose size we sharply characterize, is in general necessary for asymptotic equivalence to hold.


title: UMAP: Uniform Manifold Approximation and Projection for Dimension Reduction

abstract: UMAP (Uniform Manifold Approximation and Projection) is a novel manifold learning technique for dimension reduction. UMAP is constructed from a theoretical framework based in Riemannian geometry and algebraic topology. The result is a practical scalable algorithm that applies to real world data. The UMAP algorithm is competitive with t-SNE for visualization quality, and arguably preserves more of the global structure with superior run time performance. Furthermore, UMAP has no computational restrictions on embedding dimension, making it viable as a general purpose dimension reduction technique for machine learning.


title: Structural matrix algebras, generalized flags and gradings

abstract: We show that a structural matrix algebra $A$ is isomorphic to the endomorphism algebra of an algebraic-combinatorial object called a generalized flag. If the flag is equipped with a group grading, an algebra grading is induced on $A$. We classify the gradings obtained in this way as the orbits of the action of a double semidirect product on a certain set. Under some conditions on the associated graph, all good gradings on $A$ are of this type. As a bi-product, we obtain a new approach to compute the automorphism group of a structural matrix algebra.


title: Cover time for the frog model on trees

abstract: The frog model is a branching random walk on a graph in which particles branch only at unvisited sites. Consider an initial particle density of $\mu$ on the full $d$-ary tree of height $n$. If $\mu= \Omega( d^2)$, all of the vertices are visited in time $\Theta(n\log n)$ with high probability. Conversely, if $\mu = O(d)$ the cover time is $\exp(\Theta(\sqrt n))$ with high probability.


title: Hybrid Controller for Wind Turbine Generators to Ensure Adequate Frequency Response in Power Networks

abstract: Converter-interfaced power sources (CIPS) are hybrid control systems as they may switch between multiple operating modes. Due to increasing penetration, the hybrid behavior of CIPS, such as, wind turbine generators (WTG), may have significant impact on power system dynamics. In this paper, the frequency dynamics under inertia emulation and primary support from WTG is studied. A mode switching for WTG to ensure adequate frequency response is proposed. The switching instants are determined by our proposed concept of a region of safety (ROS), which is the initial set of safe trajectories. The barrier certificate methodology is employed to derive a new algorithm to obtain and enlarge the ROS for the given desired safe limits and the worst-case disturbance scenarios. Then critical switching instants and a safe recovery procedure are found. In addition, the emulated inertia and load-damping effect is derived in the time frame of inertia and primary frequency response, respectively. The theoretical results under critical cases are consistent with simulations and can be used as guidance for practical control design.


title: Coded Sparse Matrix Multiplication

abstract: In a large-scale and distributed matrix multiplication problem $C=A^{\intercal}B$, where $C\in\mathbb{R}^{r\times t}$, the coded computation plays an important role to effectively deal with "stragglers" (distributed computations that may get delayed due to few slow or faulty processors). However, existing coded schemes could destroy the significant sparsity that exists in large-scale machine learning problems, and could result in much higher computation overhead, i.e., $O(rt)$ decoding time. In this paper, we develop a new coded computation strategy, we call \emph{sparse code}, which achieves near \emph{optimal recovery threshold}, \emph{low computation overhead}, and \emph{linear decoding time} $O(nnz(C))$. We implement our scheme and demonstrate the advantage of the approach over both uncoded and current fastest coded strategies.


title: Extremal digraphs avoiding an orientation of $C_4$

abstract: Let $P_{2,2}$ be the orientation of $C_4$ which consists of two 2-paths with the same initial and terminal vertices. In this paper, we determine the maximum size of $P_{2,2}$-free digraphs of order $n$ as well as the extremal digraphs attaining the maximum size when $n\ge 13$.


title: Asymptotic analysis and energy quantization for the Lane-Emden problem in dimension two

abstract: We complete the study of the asymptotic behavior, as $p\rightarrow +\infty$, of the positive solutions to \[ \left\{\begin{array}{lr}-\Delta u= u^p & \mbox{in}\Omega\\ u=0 &\mbox{on}\partial \Omega \end{array}\right. \] when $\Omega$ is any smooth bounded domain in $\mathbb R^2$, started in [4]. In particular we show quantization of the energy to multiples of $8\pi e$ and prove convergence to $\sqrt{e}$ of the $L^{\infty}$-norm, thus confirming the conjecture made in [4].


title: GPU Accelerated Finite Element Assembly with Runtime Compilation

abstract: In recent years, high performance scientific computing on graphics processing units (GPUs) have gained widespread acceptance. These devices are designed to offer massively parallel threads for running code with general purpose. There are many researches focus on finite element method with GPUs. However, most of the works are specific to certain problems and applications. Some works propose methods for finite element assembly that is general for a wide range of finite element models. But the development of finite element code is dependent on the hardware architectures. It is usually complicated and error prone using the libraries provided by the hardware vendors. In this paper, we present architecture and implementation of finite element assembly for partial differential equations (PDEs) based on symbolic computation and runtime compilation technique on GPU. User friendly programming interface with symbolic computation is provided. At the same time, high computational efficiency is achieved by using runtime compilation technique. As far as we know, it is the first work using this technique to accelerate finite element assembly for solving PDEs. Experiments show that a one to two orders of speedup is achieved for the problems studied in the paper.


title: Superconducting properties of under- and over-doped Ba$_{x}$K$_{1-x}$BiO$_{3}$ perovskite oxide

abstract: In the present study, we investigate the thermodynamic properties of the Ba$_{x}$K$_{1-x}$BiO$_{3}$ (BKBO) superconductor in the under- ($x=0.5$) and over-doped ($x=0.7$) regime, within the framework of the Migdal-Eliashberg formalism. The analysis is conducted to verify that the electron-phonon pairing mechanism is responsible for the induction of the superconducting phase in the mentioned compound. In particular, we show that BKBO is characterized by the relatively high critical value of the Coulomb pseudopotential, which changes with doping level and does not follow the Morel-Anderson model. In what follows, the corresponding superconducting band gap size and related dimensionless ratio are estimated to increase with the doping, in agreement with the experimental predictions. Moreover the effective mass of electrons is found to take on high values in the entire doping and temperature region. Finally, the characteristic dimensionless ratios for the superconducting band gap, the critical magnetic field and the specific heat for the superconducting state are predicted to exceed the limits set within the Bardeen-Cooper-Schrieffer theory, suggesting pivotal role of the strong-coupling and retardation effects in the analyzed compound. Presented results supplement our previous investigations and account for the strong-coupling phonon-mediated character of the superconducting phase in BKBO at any doping level.


title: Mean-field Games for Bio-inspired Collective Decision-making in Dynamical Networks

abstract: Given a large number of homogeneous players that are distributed across three possible states, we consider the problem in which these players have to control their transition rates, while minimizing a cost. The optimal transition rates are based on the players' knowledge of their current state and of the distribution of all the other players, and this introduces mean-field terms in the running and the terminal cost. The first contribution involves a mean-field game model that brings together macroscopic and microscopic dynamics. We obtain the mean-field equilibrium associated with this model, by solving the corresponding initial-terminal value problem. We perform an asymptotic analysis to obtain a stationary equilibrium for the system. The second contribution involves the study of the microscopic dynamics of the system for a finite number of players that interact in a structured environment modeled by an interaction topology. The third contribution is the specialization of the model to describe honeybee swarms, virus propagation, and cascading failures in interconnected smart-grids. A numerical analysis is conducted which involves two types of cyber-attacks. We simulate in which ways failures propagate across the interconnected smart grids and the impact on the grids frequencies. We reframe our analysis within the context of Lyapunov's linearisation method and stability theory of nonlinear systems and Kuramoto coupled oscillators model.


title: The language (and series) of Hammersley-type processes

abstract: We study languages and formal power series associated to (variants of) Hammersley's process. We show that the ordinary Hammersley process yields a regular language and the Hammersley tree process yields deterministic context-free (but non-regular) languages. For the extension to intervals of the Hammersley process we show that there are two relevant formal languages. One of them leads to the same class of languages as the ordinary Hammersley tree process. The other one yields non-context-free languages. The results are motivated by the problem of studying the analog of the famous Ulam-Hammersley problem for heapable sequences. Towards this goal we also give an algorithm for computing formal power series associated to the variants of Hammersley's process. We employ these algorithms to settle the nature of the scaling constant, conjectured in previous work to be the golden ratio. Our results provide experimental support to this conjecture.


title: Integers represented by positive-definite quadratic forms and Petersson inner products

abstract: Let $Q$ be a positive-definite quaternary quadratic form with integer coefficients. We study the problem of giving bounds on the largest positive integer $n$ that is locally represented by $Q$ but not represented. Assuming that $n$ is relatively prime to $D(Q)$, the determinant of the Gram matrix of $Q$, we show that $n$ is represented provided that \[ n \gg \max \{ N(Q)^{3/2 + \epsilon} D(Q)^{5/4 + \epsilon}, N(Q)^{2 + \epsilon} D(Q)^{1 + \epsilon} \}. \] Here $N(Q)$ is the level of $Q$. We give three other bounds that hold under successively weaker local conditions on $n$. These results are proven by bounding the Petersson norm of the cuspidal part of the theta series, which is accomplished using an explicit formula for the Weil representation due to Scheithauer.


title: Generalized Master-Slave-Splitting Method and Application to Transmission-Distribution Coordinated Energy Management

abstract: Transmission-Distribution coordinated energy management (TDCEM) is recognized as a promising solution to the challenge of high DER penetration, but there is a lack of a distributed computation method that universally and effectively works for the TDCEM. To bridge this gap, a generalized mas-ter-slave-splitting (G-MSS) method is presented in this paper. This method is based on a general-purpose transmis-sion-distribution coordination model called G-TDCM, which thus enables the G-MSS to be applicable to most of the central functions of the TDCEM. In this G-MSS method, a basic heter-ogenous decomposition (HGD) algorithm is first derived from the HGD of the coupling constraints in the optimality conditions of the G-TDCM. Its optimality and convergence properties are then proved. Further, inspired by the conditions for conver-gence, a modified HGD algorithm that utilizes the subsystem's response function is developed and thus converges faster. The distributed G-MSS method is then demonstrated to successfully solve a series of central functions, e.g. power flow, contingency analysis, voltage stability assessment, economic dispatch and optimal power flow, of the TDCEM. The severe issues of over-voltage and erroneous assessment of the system security that are caused by DERs are thus resolved by the G-MSS method with modest computation cost.


title: A T-shaped double quantum dot system as a Fano interferometer: interplay of coherence and correlation upon spin currents

abstract: Based on Keldysh non-equilibrium Green function method, we have investigated spin current production in a hybrid T-shaped device, consisting of a central quantum dot connected to the leads and a side dot which only couples to the central dot. The topology of this structure allows for quantum interference of the different paths that go across the device, yielding Fano resonances in the spin dependent transport properties. Correlation effects are taken into account at the central dot and handled within a mean field approximation. Its interplay with the Fano effect is analyzed in the strong coupling regime. Non-vanishing spin currents are only obtained when the leads are ferromagnetic, the current being strongly dependent on the relative orientation of the lead polarizations. We calculate the conductance (spin and charge) by numerically differentiating the current, and a rich structure is obtained as a manifestation of quantum coherence and correlation effects. Increase of the Coulomb interaction produces localization of states at the side dot, largely suppressing Fano resonances. The interaction is also responsible for the negative values of the spin conductance in some regions of the voltage near resonances, effect which is the spin analog of the Esaki tunnel diode. We also analyze control of the currents via gate voltages applied to the dots, possibility which is interesting for practical operations.


title: Constraining Stellar Coronal Mass Ejections Through Multi-Wavelength Analysis of the Active M Dwarf EQ Peg

abstract: Stellar coronal mass ejections remain experimentally unconstrained, unlike their stellar flare counterparts which are observed ubiquitously across the electromagnetic spectrum. Low frequency radio bursts in the form of a type II burst offer the best means of identifying and constraining the rate and properties of stellar CMEs. CME properties can be further improved through the use of proposed solar-stellar scaling relations and multi-wavelength observations of CMEs through the use of type II bursts and the associated flare expected to occur occur alongside. We report on 20 hours of observation of the nearby, magnetically active, and well characterized M dwarf star EQ Peg. The observations are simultaneously observed with the Jansky Very Large Array at their P-band (230-470 MHz) and at the Apache Point observatory in the SDSS u' filter ($\lambda$ = 3557 \AA). Dynamic spectra of the P band data, constructed to search for signals in the frequency-time domains, did not reveal evidence for drifting radio bursts that could be ascribed to type II bursts. Given the sensitivity of our observations, we are able to place limits on the brightness temperature and source size of any bursts which may have occurred. Using solar scaling rations on four observed stellar flares, we predict CME parameters. Given the constraints on coronal density and photospheric field strength, our models suggest that the observed flares would have been insufficient to produce detectable type II bursts at our observed frequencies. We consider the implications of these results, and other recent findings, on stellar mass loss.


title: Locally Private Hypothesis Testing

abstract: We initiate the study of differentially private hypothesis testing in the local-model, under both the standard (symmetric) randomized-response mechanism (Warner, 1965, Kasiviswanathan et al, 2008) and the newer (non-symmetric) mechanisms (Bassily and Smith, 2015, Bassily et al, 2017). First, we study the general framework of mapping each user's type into a signal and show that the problem of finding the maximum-likelihood distribution over the signals is feasible. Then we discuss the randomized-response mechanism and show that, in essence, it maps the null- and alternative-hypotheses onto new sets, an affine translation of the original sets. We then give sample complexity bounds for identity and independence testing under randomized-response. We then move to the newer non-symmetric mechanisms and show that there too the problem of finding the maximum-likelihood distribution is feasible. Under the mechanism of Bassily et al (2007) we give identity and independence testers with better sample complexity than the testers in the symmetric case, and we also propose a $\chi^2$-based identity tester which we investigate empirically.


title: Battery Energy Storage Scheduling for Optimal Load Variance Minimization

abstract: Generation portfolio can be significantly altered due to the deployment of distributed energy resources (DER) in distribution networks and the concept of microgrid. Generally, distribution networks can operate in a more resilient and economic fashion through proper coordination of DER. However, due to the partially uncontrollable and stochastic nature of some DER, the variance of net load of distribution systems increases, which raises the operational cost and complicates operation for transmission companies. This motivates peak shaving and valley filling using energy storage units deployed in distribution systems. This paper aims at theoretical formulation of optimal load variance minimization, where the infinity norm of net load is minimized. Then, the problem is reformulated equivalently as a linear program. A case study is performed with capacity-limited battery energy storage model and the simplified power flow model of a radial distribution network. The influence of capacity limit and deployment location are studied.


title: On a transformation of Riordan moment sequences

abstract: We define a transformation that associates certain exponential moment sequences with ordinary moment sequences in a natural way. The ingredients of this transformation are series reversion, the Sumudu transform (a variant of the Laplace transform), and the inverting of generating functions. This transformation also has a simple interpretation in terms of continued fractions. It associates lattice path objects with permutation objects, and in particular it associates the Narayana triangle with the Eulerian triangle.


title: Using the existence of t-designs to prove Erd\H{o}s-Ko-Rado

abstract: In 1984, Wilson proved the Erd\H{o}s-Ko-Rado theorem for $t$-intersecting families of $k$-subsets of an $n$-set: he showed that if $n\ge(t+1)(k-t+1)$ and $\mathcal{F}$ is a family of $k$-subsets of an $n$-set such that any two members of $\mathcal{F}$ have at least $t$ elements in common, then $|\mathcal{F}|\le\binom{n-t}{k-t}$. His proof made essential use of a matrix whose origin is not obvious. In this paper we show that this matrix can be derived, in a sense, as a projection of $t$-$(n,k,1)$ design.


title: Difference equations related to Jacobi-type pencils

abstract: In this paper we study various difference equations related to Jacobi-type pencils. By a Jacobi-type pencil one means the following pencil: $J_5 - \lambda J_3$, where $J_3$ is a Jacobi matrix and $J_5$ is a semi-infinite real symmetric five-diagonal matrix with positive numbers on the second subdiagonal. The basic set of solutions for the corresponding $4$-th order difference equation is constructed. Spectral properties of the truncated pencil and some special matrix orthogonality relations are investigated. Classical type orthogonal polynomials satisfying a $4$-th order differential equation are constructed.


title: Pros and Cons of GAN Evaluation Measures

abstract: Generative models, in particular generative adversarial networks (GANs), have received significant attention recently. A number of GAN variants have been proposed and have been utilized in many applications. Despite large strides in terms of theoretical progress, evaluating and comparing GANs remains a daunting task. While several measures have been introduced, as of yet, there is no consensus as to which measure best captures strengths and limitations of models and should be used for fair model comparison. As in other areas of computer vision and machine learning, it is critical to settle on one or few good measures to steer the progress in this field. In this paper, I review and critically discuss more than 24 quantitative and 5 qualitative measures for evaluating generative models with a particular emphasis on GAN-derived models. I also provide a set of 7 desiderata followed by an evaluation of whether a given measure or a family of measures is compatible with them.


title: WIYN Open Cluster Study, LXXVI. Li Evolution Among Stars of Low/Intermediate Mass: The Metal-Deficient Open Cluster, NGC 2506

abstract: HYDRA spectra of 287 stars in the field of NGC 2506 from the turnoff through the giant branch are analyzed. With previous data, 22 are identified as probable binaries; 90 more are classified as potential non-members. Spectroscopic analyses of ~60 red giants and slowly rotating turnoff stars using line equivalent widths and a neural network approach lead to [Fe/H] = -0.27 +/- 0.07 (s.d.) and [Fe/H] = -0.27 +/- 0.06 (s.d.), respectively. Li abundances are derived for 145 probable single-star members, 44 being upper limits. Among turnoff stars outside the Li-dip, A(Li) = 3.04 +/- 0.16 (s.d.), with no trend with color, luminosity, or rotation speed. Evolving from the turnoff across the subgiant branch, there is a well-delineated decline to A(Li) ~1.25 at the giant branch base, coupled with the rotational spindown from between ~20 and 70 km/s to less than 20 km/s for stars entering the subgiant branch and beyond. A(Li) remains effectively constant from the giant branch base to the red giant clump level. A new member above the clump redefines the path of the first-ascent red giant branch; its Li is 0.6 dex below the first-ascent red giants. With one exception, all post-He-flash stars have upper limits to A(Li), at or below the level of the brightest first-ascent red giant. The patterns are in excellent qualitative agreement with the model predictions for low/intermediate-mass stars which undergo rotation-induced mixing at the turnoff and subgiant branch, first dredge-up, and thermohaline mixing beyond the red giant bump.


title: Skew braces and the Galois correspondence for Hopf Galois structures

abstract: Let $L/K$ be a Galois extension of fields with Galois group $\Gamma$, and suppose $L/K$ is also an $H$-Hopf Galois extension. Using the recently uncovered connection between Hopf Galois structures and skew left braces, we introduce a method to quantify the failure of surjectivity of the Galois correspondence from subHopf algebras of $H$ to intermediate subfields of $L/K$, given by the Fundamental Theorem of Hopf Galois Theory. Suppose $L \otimes_K H = LN$ where $N \cong (G, \star)$. Then there exists a skew left brace $(G, \star, \circ)$ where $(G, \circ) \cong \Gamma$. We show that there is a bijective correspondence between intermediate fields $E$ between $K$ and $L$ and certain sub-skew left braces of $G$, which we call the $\circ$-stable subgroups of $(G, \star)$. Counting these subgroups and comparing that number with the number of subgroups of $\Gamma \cong (G, \circ)$ describes how far the Galois correspondence for the $H$-Hopf Galois structure is from being surjective. The method is illustrated by a variety of examples.


title: Accuracy of the microcanonical Lanczos method to compute real-frequency dynamical spectral functions of quantum models at finite temperatures

abstract: We examine the accuracy of the microcanonical Lanczos method (MCLM) developed by Long, {\it et al.} [Phys. Rev. B {\bf 68}, 235106 (2003)] to compute dynamical spectral functions of interacting quantum models at finite temperatures. The MCLM is based on the microcanonical ensemble, which becomes exact in the thermodynamic limit. To apply the microcanonical ensemble at a fixed temperature, one has to find energy eigenstates with the energy eigenvalue corresponding to the internal energy in the canonical ensemble. Here, we propose to use thermal pure quantum state methods by Sugiura and Shimizu [Phys. Rev. Lett. {\bf 111}, 010401 (2013)] to obtain the internal energy. After obtaining the energy eigenstates using the Lanczos diagonalization method, dynamical quantities are computed via a continued fraction expansion, a standard procedure for Lanczos-based numerical methods. Using one-dimensional antiferromagnetic Heisenberg chains with $S=1/2$, we demonstrate that the proposed procedure is reasonably accurate even for relatively small systems.


title: Minimizing Latency to Support VR Social Interactions over Wireless Cellular Systems via Bandwidth Allocation

abstract: Immersive social interactions of mobile users are soon to be enabled within a virtual space, by means of virtual reality (VR) technologies and wireless cellular systems. In a VR mobile social network, the states of all interacting users should be updated synchronously and with low latency via two-way communications with edge computing servers. The resulting end-to-end latency depends on the relationship between the virtual and physical locations of the wireless VR users and of the edge servers. In this work, the problem of analyzing and optimizing the end-to-end latency is investigated for a simple network topology, yielding important insights into the interplay between physical and virtual geometries.


title: Estimating the Spectral Density of Large Implicit Matrices

abstract: Many important problems are characterized by the eigenvalues of a large matrix. For example, the difficulty of many optimization problems, such as those arising from the fitting of large models in statistics and machine learning, can be investigated via the spectrum of the Hessian of the empirical loss function. Network data can be understood via the eigenstructure of a graph Laplacian matrix using spectral graph theory. Quantum simulations and other many-body problems are often characterized via the eigenvalues of the solution space, as are various dynamic systems. However, naive eigenvalue estimation is computationally expensive even when the matrix can be represented; in many of these situations the matrix is so large as to only be available implicitly via products with vectors. Even worse, one may only have noisy estimates of such matrix vector products. In this work, we combine several different techniques for randomized estimation and show that it is possible to construct unbiased estimators to answer a broad class of questions about the spectra of such implicit matrices, even in the presence of noise. We validate these methods on large-scale problems in which graph theory and random matrix theory provide ground truth.


title: Learning Local Metrics and Influential Regions for Classification

abstract: The performance of distance-based classifiers heavily depends on the underlying distance metric, so it is valuable to learn a suitable metric from the data. To address the problem of multimodality, it is desirable to learn local metrics. In this short paper, we define a new intuitive distance with local metrics and influential regions, and subsequently propose a novel local metric learning method for distance-based classification. Our key intuition is to partition the metric space into influential regions and a background region, and then regulate the effectiveness of each local metric to be within the related influential regions. We learn local metrics and influential regions to reduce the empirical hinge loss, and regularize the parameters on the basis of a resultant learning bound. Encouraging experimental results are obtained from various public and popular data sets.


title: On variational solutions for whole brain serial-section histology using the computational anatomy random orbit model

abstract: This paper presents a variational framework for dense diffeomorphic atlas-mapping onto high-throughput histology stacks at the 20 um meso-scale. The observed sections are modelled as Gaussian random fields conditioned on a sequence of unknown section by section rigid motions and unknown diffeomorphic transformation of a three-dimensional atlas. To regularize over the high-dimensionality of our parameter space (which is a product space of the rigid motion dimensions and the diffeomorphism dimensions), the histology stacks are modelled as arising from a first order Sobolev space smoothness prior. We show that the joint maximum a-posteriori, penalized-likelihood estimator of our high dimensional parameter space emerges as a joint optimization interleaving rigid motion estimation for histology restacking and large deformation diffeomorphic metric mapping to atlas coordinates. We show that joint optimization in this parameter space solves the classical curvature non-identifiability of the histology stacking problem. The algorithms are demonstrated on a collection of whole-brain histological image stacks from the Mouse Brain Architecture Project.


title: Hybrid integration of silicon photonics circuits and InP lasers by photonic wire bonding

abstract: Efficient coupling of III-V light sources to silicon photonic circuits is one of the key challenges of integrated optics. Important requirements are low coupling losses, as well as small footprint and high yield of the overall assembly, along with the ability to use automated processes for large-scale production. In this paper, we demonstrate that photonic wire bonding addresses these challenges by exploiting direct-write two-photon lithography for in-situ fabrication of three-dimensional freeform waveguides between optical chips. In a series proof-of-concept experiments, we connect InP-based horizontal-cavity surface emitting lasers (HCSEL) to passive silicon photonic circuits with insertion losses down to 0.4 dB. To the best of our knowledge, this is the most efficient interface between an InP light source and a silicon photonic chip that has so far been demonstrated. Our experiments represent a key step in advancing photonic wire bonding to a universal integration platform for hybrid photonic multi-chip assemblies that combine known-good dies of different materials to high-performance hybrid multi-chip modules.


title: Don't Repeat Yourself: Seamless Execution and Analysis of Extensive Network Experiments

abstract: This paper presents MACI, the first bespoke framework for the management, the scalable execution, and the interactive analysis of a large number of network experiments. Driven by the desire to avoid repetitive implementation of just a few scripts for the execution and analysis of experiments, MACI emerged as a generic framework for network experiments that significantly increases efficiency and ensures reproducibility. To this end, MACI incorporates and integrates established simulators and analysis tools to foster rapid but systematic network experiments. We found MACI indispensable in all phases of the research and development process of various communication systems, such as i) an extensive DASH video streaming study, ii) the systematic development and improvement of Multipath TCP schedulers, and iii) research on a distributed topology graph pattern matching algorithm. With this work, we make MACI publicly available to the research community to advance efficient and reproducible network experiments.


title: Fabry-P\'erot Cavity Formed with Dielectric Metasurfaces in a Hollow-Core Fiber

abstract: We demonstrate a fiber-integrated Fabry-P\'erot cavity formed by attaching a pair of dielectric metasurfaces to the ends of a hollow-core photonic-crystal fiber segment. The metasurfaces consist of perforated membranes designed as photonic-crystal slabs that act as planar mirrors but can potentially allow injection of gases through their holes into the hollow core of the fiber. We have so far observed cavities with finesse of ~11 and Q factors of ~$4.5 \times 10^5$, but much higher values should be achievable with improved fabrication procedures. We expect this device to enable development of new fiber lasers, enhanced gas spectroscopy, and studies of fundamental light-matter interactions and nonlinear optics.


title: Bayesian Compressive Sensing with Circulant Matrix for Spectrum Sensing in Cognitive Radio Networks

abstract: For wideband spectrum sensing, compressive sensing has been proposed as a solution to speed up the high dimensional signals sensing and reduce the computational complexity. Compressive sensing consists of acquiring the essential information from a sparse signal and recovering it at the receiver based on an efficient sampling matrix and a reconstruction technique. In order to deal with the uncertainty, improve the signal acquisition performance, and reduce the randomness during the sensing and reconstruction processes, compressive sensing requires a robust sampling matrix and an efficient reconstruction technique. In this paper, we propose an approach that combines the advantages of a Circulant matrix with Bayesian models. This approach is implemented, extensively tested, and its results have been compared to those of l1 norm minimization with a Circulant or random matrix based on several metrics. These metrics are Mean Square Error, reconstruction error, correlation, recovery time, sampling time, and processing time. The results show that our technique is faster and more efficient.


title: Vector particle scattering on the lattice

abstract: In this work, we present an explicit form of the Luescher equation and consider the construction of the operators in different irreducible representations for the case of scattering of two vector particles. The formalism is applied to scalar QED in the Higgs Phase, where the $U(1)$ gauge boson acquires mass.


title: On-chip label-free plasmonic based imaging microscopy for microfluidics

abstract: In this work, we demonstrated on-chip label-free imaging microscopy using real and Fourier Plane (FP) microscopic dark field images of Surface Plasmons (SP), excited on engineered 1D and 2D low aspect ratio periodic plasmonic nanostructures. The periodic plasmonic nanostructures with period almost equal to resonance wavelength, were engineered to exhibit transmission resonances as transmission peaks in visible spectrum at normal incidence, without using extraordinary optical transmission phenomena. The plasmonic nanostructures exhibited a polarization rotation of 90 degree mediated by differential phase retardation in the SP mode due to Transverse Electric and Magnetic components. This was used to develop a dark field on-chip plasmonic polarization microscope for imaging SP excitation in real and Fourier planes. After successful integration of these plasmonic nanostructures with SU- 8 based microfluidic channels, a real-time monitoring of label free on chip sensing was demonstrated. Label-free on-chip imaging for interface of colorless miscible and immiscible analytes flowing on plasmonic nanostructures in the microfluidic channels were performed using color-selective filtering nature of plasmonic nanostructures. Since the imaging is realized on a chip and does not need any complicated and bulky arrangement, it will be well suited for on-chip point of care diagnostics.


title: Electrochemical and mechanical behaviors of dissimilar friction stir welding between 5086 and 6061 aluminum alloy

abstract: The electrochemical behavior and mechanical properties of friction stir welded AA5086 and AA6061 Al alloys were investigated. Micro-hardness measurements and tensile tests showed that the heat-affected zone (HAZ) in AA6061 had minimum hardness value (i.e., 88 HV) and served as failure site in the dissimilar weld. Corrosion testing revealed that the minimum value of Icorr appeared in the HAZ 5086 (0.54 uA/cm2) and HAZ 5086 was most resistant to corrosion. The AA 5086 side of the weld showed better corrosion resistance than the AA 6061 side.


title: Suppression of bend loss in writing of three-dimensional optical waveguides with femtosecond laser pulses

abstract: We provide a solution toward compact and low-loss three dimensional (3D) photonic circuits by femtosecond laser direct writing of 3D waveguides in fused silica. We suppress the bend loss by sandwiching the waveguide between a pair of walls formed by internal modification of glass. Our method allows to reduce the bend loss of a curved waveguide with a bending radius of 15 mm by more than one order of magnitude.


title: OAT: Attesting Operation Integrity of Embedded Devices

abstract: Due to the wide adoption of IoT/CPS systems, embedded devices(IoT frontends) become increasingly connected and mission-critical, which in turn has attracted advanced attacks (e.g., control-flow hijacks and data-only attacks). Unfortunately, IoT backends are unable to detect if such attacks have happened while receiving data, service requests, or operation status from IoT devices. As a result, currently, IoT backends are forced to blindly trust the IoT devices that they interact with. To fill this void, we first formulate a new security property for embedded devices, called "Operation Execution Integrity" or OEI. We then design and build a system, OAT, that enables remote OEI attestation for ARM-based bare-metal embedded devices. Our formulation of OEI captures the integrity of both control flow and critical data involved in an operation execution. Therefore, satisfying OEI entails that an operation execution is free of unexpected control and data manipulations, which existing attestation methods cannot check. Our design of OAT strikes a balance between prover's constraints (embedded devices' limited computing power and storage) and verifier's requirements(complete verifiability and forensic assistance). OAT uses a new control-flow measurement scheme, which enables light-weight and space-efficient collection of measurements (97% space reduction from the trace-based approach). OAT performs the remote control-flow verification through abstract execution, which is fast and deterministic. OAT also features lightweight integrity checking for critical data (74% fewer instrumentation needed than previous work). Our security analysis shows that OAT allows remote verifiers or IoT backends to detect both control-flow hijacks and data-only attacks that affect the execution of operations on IoT devices. In our evaluation using real embedded programs, OAT incurs a runtime overhead of 2.7%.


title: Using Intelligent Control to Improve PV Systems Efficiency

abstract: One of the important technologies in renewable energy is the Photovoltaic (PV), which is the direct conversion of light into electricity at the atomic level, and improving the efficiency of PV systems is one of the very important factors in getting the most out of this invaluable renewable resource of energy. While most research work we studied used conventional techniques to control two parameters at most, like power and change in power, or voltage and change in voltage, etc., we implemented unconventional techniques, namely intelligent control to control more than two parameters at a time, including change in temperature which had been ignored by many researchers for various reasons, as well as the use of probability theory to predict the location of power point and control how it would move before not after it did. Practically, we utilized available PV systems devices to test our controlled systems and we used simulations and compare our findings with previous work done by others in this area and our techniques showed good improvement in efficiency and we believe that it could open the door for other colleagues to add valuable work in this important field.


title: Metric Learning via Maximizing the Lipschitz Margin Ratio

abstract: In this paper, we propose the Lipschitz margin ratio and a new metric learning framework for classification through maximizing the ratio. This framework enables the integration of both the inter-class margin and the intra-class dispersion, as well as the enhancement of the generalization ability of a classifier. To introduce the Lipschitz margin ratio and its associated learning bound, we elaborate the relationship between metric learning and Lipschitz functions, as well as the representability and learnability of the Lipschitz functions. After proposing the new metric learning framework based on the introduced Lipschitz margin ratio, we also prove that some well known metric learning algorithms can be shown as special cases of the proposed framework. In addition, we illustrate the framework by implementing it for learning the squared Mahalanobis metric, and by demonstrating its encouraging results on eight popular datasets of machine learning.


title: Monopoles of Twelve Types in 3-Body Problems

abstract: We consider twelve different ways of modelling the 3-body problem in dimension $\geq 2$. These can be viewed as models of classical and quantum background independence. We show that a different type of monopole is realized in each's relational space: a type of reduced configuration space. 8 cases occur in 2-$d$, and 4 distinct ones in 3-$d$; these reflect counts of non-equivalent subgroup actions of $S_3 \times C_2$ and $S_3$ respectively. The $S_3$ acts on particle labels; the extra $C_2$ corresponds to the purely 2-$d$ option of whether or not to identify mirror images. The non-equivalent realization is due to a suite of subgroup, orbit space and stratification features. Our 2-$d$ monopoles include 4 known ones: a realization of Dirac's monopole in relational space rather than its more habitual setting of space, the 2-$d$ version of Iwai's monopole, and indistinguishable particle monopoles with and without mirror image identification. The 4 new ones are indistinguishable under a 2-particle label switch or under even permutations, in each case with optional mirror image identification. Our 4 3-$d$ monopoles are 2 known ones: the actual Iwai monopole and its already-announced indistinguishable-particles counterpart, and 2 new ones: the two-particle label switch and even permutation cases. All 4 3-$d$ cases are stratified. The three even-permutation cases are orbifolds, two with boundary, the 3-$d$ case's boundary constituting a separate stratum, giving a stratified orbifold. We document each of the 12 cases' underlying shape space and relational space, and each monopole's Hopf mathematics, global-section versus topological quantization dichotomy, Dirac string positioning, and Chern integral concordance with topological contributions form of Gauss--Bonnet Theorem.


title: Microscopic modeling of gas-surface scattering. I. A combined molecular dynamics-rate equation approach

abstract: A combination of first principle molecular dynamics (MD) simulations with a rate equation model (MD-RE approach) is presented to study the trapping and the scattering of rare gas atoms from metal surfaces. The temporal evolution of the atom fractions that are either adsorbed or scattered into the continuum is investigated in detail. We demonstrate that for this description one has to consider trapped, quasi-trapped and scattering states, and present an energetic definition of these states. The rate equations contain the transition probabilities between the states. We demonstrate how these rate equations can be derived from kinetic theory. Moreover, we present a rigorous way to determine the transition probabilities from a microscopic analysis of the particle trajectories generated by MD simulations. Once the system reaches quasi-equilibrium, the rates converge to stationary values, and the subsequent thermal adsorption/desorption dynamics is completely described by the rate equations without the need to perform further time-consuming MD simulations. As a proof of concept of our approach, MD simulations for argon atoms interacting with a platinum (111) surface are presented. A detailed deterministic trajectory analysis is performed, and the transition rates are constructed. The dependence of the rates on the incidence conditions and the lattice temperature is analyzed. Based on this example, we analyze the time scale of the gas-surface system to approach the quasi-stationary state. The MD-rate equation model has great relevance for the plasma-surface modeling as it makes an extension of accurate simulations to long, experimentally relevant time scales possible. Its application to the computation of atomic sticking probabilities is given in the second part (paper II).


title: The role of defects in the etching of graphene by intercalated oxygen

abstract: Graphene is one of the most promising 2D materials for various applications due to its unique electronic properties and high thermal stability. In previous studies, it was shown that when graphene is deposited onto some transition metal substrates, small molecules, such as O$_2$, intercalate between the graphene and the substrate and react to partially etch the graphene film when heated to desorb the intercalates. Here, carbon vacancy defects are intentionally formed on Gr/Ru(0001) and their effect on the intercalation of oxygen and etching of the graphene layer are investigated. 50 eV Ar$^+$ sputtering with a low fluence is used to create isolated single vacancy defects in the graphene overlayer and helium low energy ion scattering (LEIS) is employed for surface analysis. It is found that the defects both ease the intercalation of the oxygen molecules and improve the etching efficiency of the graphene during annealing.


title: Three-dimensional modeling of chiral nematic texture evolution under electric switching

abstract: Chiral nematic liquid crystals exhibit both a helical planar ground state with uniform twist and a metastable defect-rich focal conic texture, and can be switched between the two microstructures via application of transient voltage pulses. In this work, we model these electrically-induced texture transitions using finite difference methods to examine resulting microstructural evolution, the first time this transition has been modeled in three dimensions. We analyze the planar to focal conic, focal conic to planar, and planar to planar transitions depending on voltage pulse magnitude. We consider first the special case of chiral nematics with matched twist and bend elastic constants. Results show a variety of defect-rich morphologies in the disordered focal conic texture and demonstrate a fast recovery of the planar ground state on switching without formation of a transient planar state. We evaluate both texture microstructural evolution as well as cell capacitance. Beyond the single elastic constant approximation, we evaluate the planar to transient-planar as well as the planar to Helfrich-deformed transitions in simulations of a liquid crystal compound with different elastic constants. Our methods represent the evolving microstructure as a uniaxial director field, with relaxation dynamics calculated from a tensor representation so that half charge disclination defects are not suppressed. We discuss potential application of these computationally efficient three-dimensional modeling approaches for design and optimization of chiral nematic devices.


title: Cosmological Lithium Problems

abstract: We briefly describe the cosmological lithium problems followed by a summary of our recent theoretical work on the magnitude of the effects of electron screening, the possible existence of dark matter parallel universes and the use of non-extensive (Tsallis) statistics during big bang nucleosynthesis. Solutions within nuclear physics are also discussed and recent measurements of cross-sections based on indirect experimental techniques are summarized.


title: Non-Relativistic Fermion-Fermion Scattering in Higher Derivative Gravity

abstract: In this note, we examine the scattering of two identical fermions in theories where fermionic fields minimally coupled to higher derivative gravity. In particular, we consider the extension of general relativity with $R^2$ corrections or non-local terms. We expand the action of fermions around the flat space background and obtain two fermion-one graviton vertex. Then, by considering the scattering amplitude of two fermions, we calculate the non-relativistic limit and that obtain the potential for two fermion-fermion interaction which would be the usual Newtonian potential corrected with a Yukawa-like term. At the end, we briefly discuss the astronomical effects of such Yukawa-like potential by computing the gravitational pressure of a spherical star and use it for a white dwarf to obtain quantum corrections of Chandrasekhar radius.


title: Certified Robustness to Adversarial Examples with Differential Privacy

abstract: Adversarial examples that fool machine learning models, particularly deep neural networks, have been a topic of intense research interest, with attacks and defenses being developed in a tight back-and-forth. Most past defenses are best effort and have been shown to be vulnerable to sophisticated attacks. Recently a set of certified defenses have been introduced, which provide guarantees of robustness to norm-bounded attacks, but they either do not scale to large datasets or are limited in the types of models they can support. This paper presents the first certified defense that both scales to large networks and datasets (such as Google's Inception network for ImageNet) and applies broadly to arbitrary model types. Our defense, called PixelDP, is based on a novel connection between robustness against adversarial examples and differential privacy, a cryptographically-inspired formalism, that provides a rigorous, generic, and flexible foundation for defense.


title: Modeling of Teager Energy Operated Perceptual Wavelet Packet Coefficients with an Erlang-2 PDF for Real Time Enhancement of Noisy Speech

abstract: In this paper, for real time enhancement of noisy speech, a method of threshold determination based on modeling of Teager energy (TE) operated perceptual wavelet packet (PWP) coefficients of the noisy speech and noise by an Erlang-2 PDF is presented. The proposed method is computationally much faster than the existing wavelet packet based thresholding methods. A custom thresholding function based on a combination of mu-law and semisoft thresholding functions is designed and exploited to apply the statistically derived threshold upon the PWP coefficients. The proposed custom thresholding function works as a mu-law or a semisoft thresholding function or their combination based on the probability of speech presence and absence in a subband of the PWP transformed noisy speech. By using the speech files available in NOIZEUS database, a number of simulations are performed to evaluate the performance of the proposed method for speech signals in the presence of Gaussian white and street noises. The proposed method outperforms some of the state-of-the-art speech enhancement methods both at high and low levels of SNRs in terms of standard objective measures and subjective evaluations including formal listening tests.


title: Low-frequency pulse profile variation in PSR B2217+47: evidence for echoes from the interstellar medium

abstract: We have observed a complex and continuous change in the integrated pulse profile of PSR B2217+47, manifested as additional components trailing the main peak. These transient components are detected over 6 years at $150$ MHz using the LOw Frequency ARray (LOFAR), but they are not seen in contemporaneous Lovell observations at $1.5$ GHz. We argue that propagation effects in the ionized interstellar medium (IISM) are the most likely cause. The putative structures in the IISM causing the profile variation are roughly half-way between the pulsar and the Earth and have transverse radii $R \sim 30$ AU. We consider different models for the structures. Under the assumption of spherical symmetry, their implied average electron density is $\overline{n}_e \sim 100$ cm$^{-3}$. Since PSR B2217+47 is more than an order of magnitude brighter than the average pulsar population visible to LOFAR, similar profile variations would not have been identified in most pulsars, suggesting that subtle profile variations in low-frequency profiles might be more common than we have observed to date. Systematic studies of these variations at low frequencies can provide a new tool to investigate the proprieties of the IISM and the limits to the precision of pulsar timing.


title: Skyglow Changes Over Tucson, Arizona, Resulting From A Municipal LED Street Lighting Conversion

abstract: The transition from earlier lighting technologies to white light-emitting diodes (LEDs) is a significant change in the use of artificial light at night. LEDs emit considerably more short-wavelength light into the environment than earlier technologies on a per-lumen basis. Radiative transfer models predict increased skyglow over cities transitioning to LED unless the total lumen output of new lighting systems is reduced. The City of Tucson, Arizona (U.S.), recently converted its municipal street lighting system from a mixture of fully shielded high- and low-pressure sodium (HPS/LPS) luminaires to fully shielded 3000 K white LED luminaires. The lighting design intended to minimize increases to skyglow in order to protect the sites of nearby astronomical observatories without compromising public safety. This involved the migration of over 445 million fully shielded HPS/LPS lumens to roughly 142 million fully shielded 3000 K white LED lumens and an expected concomitant reduction in the amount of visual skyglow over Tucson. SkyGlow Simulator models predict skyglow decreases on the order of 10-20% depending on whether fully shielded or partly shielded lights are in use. We tested this prediction using visual night sky brightness estimates and luminance-calibrated, panchromatic all-sky imagery at 15 locations in and near the city. Data were obtained in 2014, before the LED conversion began, and in mid-2017 after approximately 95% of $\sim$18,000 luminaires was converted. Skyglow differed marginally, and in all cases with valid data changed by $<{\pm}$20%. Over the same period, the city's upward-directed optical radiance detected from Earth orbit decreased by approximately 7%. While these results are not conclusive, they suggest that LED conversions paired with dimming can reduce skyglow over cities.