title: Quantization of rationally deformed Morse potentials by Wronskian transforms of Romanovski-Bessel polynomials

abstract: The paper advances the suggestion by Odake and Sasaki to re-write eigenfunctions of rationally deformed Morse potentials in terms of Wronskians of Laguerre polynomials in the reciprocal argument. It is shown that the constructed quasi-rational seed solutions of the Schrodinger equation with the Morse potential are formed by generalized Bessel polynomials with degree-independent indexes. As a new achievement we can point to the construction of the Darboux-Crum net of isospectral rational potentials using Wronskians of generalized Bessel polynomials with no positive zeros. It was then proved that any solvable rational Darboux-Crum transform of the Morse potential can be re-expressed in terms of Wronskians of the latter polynomials accompanied by juxtaposed pairs of Romanovski-Bessel polynomials.


title: A note on Penrose limits

abstract: As a footnote to arXiv:1909.07756, I show that, given a (time-like) umbilic 3-surface $\Sigma$ in a 4-dimensional space-time $M$, the Penrose limit taken along any null geodesic $\Gamma$ which lies in $\Sigma$ is a diagonalisable plane-wave.


title: Characterizing and contrasting structural proton transport mechanisms in azole hydrogen bond networks using ab initio molecular dynamics

abstract: Imidazole and 1,2,3-triazole are promising hydrogen-bonded heterocycles that conduct protons via a structural mechanism and whose derivatives are present in systems ranging from biological proton channels to proton exchange membrane fuel cells. Here, we leverage multiple time-stepping to perform ab initio molecular dynamics of imidazole and 1,2,3-triazole at the nanosecond timescale. We show that despite the close structural similarities of these compounds, their proton diffusion constants vary by over an order of magnitude. Our simulations reveal the reasons for these differences in diffusion constants, which range from the degree of hydrogen-bonded chain linearity to the effect of the central nitrogen atom in 1,2,3-triazole on proton transport. In particular, we uncover evidence of two "blocking" mechanisms in 1,2,3-triazole, where covalent and hydrogen bonds formed by the central nitrogen atom limit the mobility of protons. Our simulations thus provide insights into the origins of the experimentally observed 10-fold difference in proton conductivity.


title: Axion Echos from the Supernova Graveyard

abstract: Stimulated decays of axion dark matter, triggered by a source in the sky, could produce a photon flux along the continuation of the line of sight, pointing backward to the source. The strength of this so-called axion "echo" signal depends on the entire history of the source and could still be strong from sources that are dim today but had a large flux density in the past, such as supernova remnants (SNRs). This echo signal turns out to be most observable in the radio band. We study the sensitivity of radio telescopes such as the Square Kilometer Array (SKA) to echo signals generated by SNRs that have already been observed, and show that SKA could reach axion-photon couplings of order $g_{a\gamma\gamma} \sim \mathcal{O}(10^{-11}) \,\mathrm{GeV}^{-1}$ for axion masses $m_a \lesssim 10^{-5}\;\mathrm{eV}$. In addition, we show projections of the detection reach for signals coming from old SNRs and from newly born supernovae that could be detected in the future. Intriguingly, an observable echo signal could come from old "ghost" SNRs which were very bright in the past but are now so dim that they haven't been observed.

categories: hep-ph astro-ph.CO astro-ph.GA astro-ph.HE


title: Measuring the Density Fields around Bright Quasars at $z\sim 6$ with XQR-30 Spectra

abstract: Measuring the density of the intergalactic medium using quasar sightlines in the epoch of reionization is challenging due to the saturation of Lyman-$\alpha$ absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Lyman-$\alpha$ absorption is not saturated. In this study, we use $10$ high-resolution ($R\gtrsim 10,000$) $z\sim 6$ quasar spectra from the extended XQR-30 sample to measure the density field in the quasar proximity zones. We find a variety of environments within $3$ pMpc distance from the quasars. We compare the observed density cumulative distribution function (CDF) with models from the $\textit{Cosmic Reionization on Computers}$ simulation, and find a good agreement between $1.5$ to $3$pMpc from the quasar. This region is far away from the quasar hosts and hence approaching the mean density of the universe, which allows us to use the CDF to set constraints on the cosmological parameter $\sigma_8=0.6\pm0.3$. The uncertainty is mainly due to the limited number of high-quality quasar sightlines currently available. Utilizing the more than $>200$ known quasars at $z\gtrsim 6$, this method will allow us in the future to tighten the constraint on $\sigma_8$ to the percent level. In the region closer to the quasar within $1.5$ pMpc, we find the density is higher than predicted in the simulation by $1.23 \pm 0.17$, suggesting the typical host dark matter halo mass of a bright quasar ($M_{\rm 1450}<-26.5$) at $z\sim 6$ is $\log_{\rm 10} (M/M_\odot)=12.5^{+0.4}_{-0.7}$.


title: Magneto-Thermal Instability In Galaxy Clusters I: Theory and Two-Dimensional Simulations

abstract: Determining the origin of turbulence in galaxy clusters, and quantifying its transport of heat, is an outstanding problem, with implications for our understanding of their thermodynamic history and structure. As the dilute plasma of the intracluster medium (ICM) is magnetized, heat and momentum travel preferentially along magnetic field lines. This anisotropy triggers a class of buoyancy instabilities that destabilize the ICM, and whose turbulent motions can augment or impede heat transport. We focus on the magneto-thermal instability (MTI), which may be active in the periphery of galaxy clusters. We aim to take a fresh look at the problem and construct a general theory that explains the MTI saturation mechanism and provides scalings and estimates for the turbulent kinetic energy, magnetic energy, and heat flux. We simulate MTI turbulence with a Boussinesq code, SNOOPY, which, in contrast to previous work, allows us to perform an extensive sampling of the parameter space. In two dimensions the saturation mechanism involves an inverse cascade carrying kinetic energy from the short MTI injection scales to larger scales, where it is arrested by the stable entropy stratification; at a characteristic "buoyancy scale", the energy is dumped into large-scale g-modes, which subsequently dissipate. Consequently, the entropy stratification sets an upper limit on the size and strength of turbulent eddies. Meanwhile, the MTI conveys a substantial fraction of heat, despite the tangled geometry of the magnetic field. In a companion paper, these results are extended to three-dimensional flows, and compared to real cluster observations.


title: The Impact of the First Galaxies on Cosmic Dawn and Reionization

abstract: The formation of the first galaxies during cosmic dawn and reionization (at redshifts $z=5-30$), triggered the last major phase transition of our universe, as hydrogen evolved from cold and neutral to hot and ionized. The 21-cm line of neutral hydrogen will soon allow us to map these cosmic milestones and study the galaxies that drove them. To aid in interpreting these observations, we upgrade the publicly available code {\tt 21cmFAST}. We introduce a new, flexible parametrization of the additive feedback from: an inhomogeneous, $H_2$-dissociating (Lyman-Werner; LW) background; and dark matter -- baryon relative velocities; which recovers results from recent, small-scale hydrodynamical simulations with both effects. We perform a large, "best-guess" simulation as the 2021 installment of the Evolution of 21-cm Structure (EOS) project. This improves the previous release with a galaxy model that reproduces the observed UV luminosity functions (UVLFs), and by including a population of molecular-cooling galaxies. The resulting 21-cm global signal and power spectrum are significantly weaker, primarily due to a more rapid evolution of the star-formation rate density required to match the UVLFs. Nevertheless, we forecast high signal-to-noise detections for both HERA and the SKA. We demonstrate how the stellar-to-halo mass relation of the unseen, first galaxies can be inferred from the 21-cm evolution. Finally, we show that the spatial modulation of X-ray heating due to relative velocities provides a unique acoustic signature that is detectable at $z \approx 10-15$ in our fiducial model. Ours are the first public simulations with joint inhomogeneous LW and relative-velocity feedback across the entire cosmic dawn and reionization.


title: Axion dark matter-induced echo of supernova remnants

abstract: Axions are a theoretically promising dark matter (DM) candidate. In the presence of radiation from bright astrophysical sources at radio frequencies, nonrelativistic DM axions can undergo stimulated decay to two nearly back-to-back photons, meaning that bright sources of radio waves will have a counterimage ("gegenschein") in nearly the exact opposite spatial direction. The counterimage will be spectrally distinct from backgrounds, taking the form of a narrow radio line centered at $\nu = m_a/4\pi$ with a width determined by Doppler broadening in the DM halo, $\Delta \nu/\nu \sim 10^{-3}$. In this work, we show that the axion decay-induced echoes of supernova remnants may be bright enough to be detectable. Their non-detection may be able to set the strongest limits to date on axion DM in the $\sim 1-10 \, \mu$eV mass range where there are gaps in coverage from existing experiments.


title: Wigner crystallization at large fine structure constant

abstract: We consider the fate of the Wigner crystal state in a two dimensional system of massive Dirac electrons as the effective fine structure constant $\alpha$ is increased. In a Dirac system, larger $\alpha$ naively corresponds to stronger electron-electron interactions, but it also implies a stronger interband dielectric response that effectively renormalizes the electron charge. We calculate the critical density and critical temperature associated with quantum and thermal melting of the Wigner crystal state using two independent approaches. We show that at $\alpha \gg 1$, the Wigner crystal state is best understood in terms of logarithmically-interacting electrons, and that both the critical density and the melting temperature approach a universal, $\alpha$-independent value. We discuss our results in the context of recent experiments in twisted bilayer graphene near the magic angle.


title: Regulation of star formation by large scale gravito-turbulence

abstract: A simple model for star formation based on supernova (SN) feedback and gravitational heating via the collapse of perturbations in gravitationally unstable disks reproduces the Schmidt-Kennicutt relation between the star formation rate (SFR) per unit area, $\Sigma_{SFR}$, and the gas surface density, $\Sigma_g$, remarkably well. The gas velocity dispersion, $\sigma_g$, is derived self-consistently in conjunction with $\Sigma_{SFR}$ and is found to match the observations. Gravitational instability triggers {"gravito-turbulence"} at the scale of the least stable perturbation mode, boosting $\sigma_g$ at $\Sigma_g> \, \Sigma_g^\textrm{thr}=50\, {\rm M}_\odot\, {\rm pc}^{-2}$, and contributing to the pressure needed to carry the disk weight vertically. $\Sigma_{SFR}$ is reduced to the observed level at $ \Sigma_g > \Sigma_g^\textrm{thr}$, whereas at lower surface densities, SN feedback is the prevailing energy source. Our proposed star formation recipes require efficiencies of order 1\%, and the Toomre parameter, $Q$, for the joint gaseous and stellar disk is predicted to be close to the critical value for marginal stability for $\Sigma_g< \, \Sigma_g^\textrm{thr}$, spreading to lower values and larger gas velocity dispersion at higher $\Sigma_g$.


title: COSMOS2020: A panchromatic view of the Universe to $z\sim10$ from two complementary catalogs

abstract: The Cosmic Evolution Survey (COSMOS) has become a cornerstone of extragalactic astronomy. Since the last public catalog in 2015, a wealth of new imaging and spectroscopic data has been collected in the COSMOS field. This paper describes the collection, processing, and analysis of this new imaging data to produce a new reference photometric redshift catalog. Source detection and multi-wavelength photometry is performed for 1.7 million sources across the $2\,\mathrm{deg}^{2}$ of the COSMOS field, $\sim$966,000 of which are measured with all available broad-band data using both traditional aperture photometric methods and a new profile-fitting photometric extraction tool, The Farmer, which we have developed. A detailed comparison of the two resulting photometric catalogs is presented. Photometric redshifts are computed for all sources in each catalog utilizing two independent photometric redshift codes. Finally, a comparison is made between the performance of the photometric methodologies and of the redshift codes to demonstrate an exceptional degree of self-consistency in the resulting photometric redshifts. The $i<21$ sources have sub-percent photometric redshift accuracy and even the faintest sources at $25<i<27$ reach a precision of $5\,\%$. Finally, these results are discussed in the context of previous, current, and future surveys in the COSMOS field. Compared to COSMOS2015, reaches the same photometric redshift precision at almost one magnitude deeper. Both photometric catalogs and their photometric redshift solutions and physical parameters will be made available through the usual astronomical archive systems (ESO Phase 3, IPAC IRSA, and CDS).


title: The PLATO field selection process I. Identification and content of the long-pointing fields

abstract: PLATO (PLAnetary Transits and Oscillations of stars) is an ESA M-class satellite planned for launch by end 2026 and dedicated to the wide-field search of transiting planets around bright and nearby stars, with a strong focus on discovering habitable rocky planets hosted by solar-like stars. The choice of the fields to be pointed at is a crucial task since it has a direct impact on the scientific return of the mission. In this paper we describe and discuss the formal requirements and the key scientific prioritization criteria that have to be taken into account in the Long-duration Observation Phase (LOP) field selection, and apply a quantitative metric to guide us in this complex optimization process. We identify two provisional LOP fields, one for each hemisphere (LOPS1, LOPN1), and discuss their properties and stellar content. While additional fine-tuning shall be applied to LOP selection before the definitive choice (to be made two years before launch), we expect their position will not move by more than a few degrees with respect to what is proposed in this paper.


title: Fractionalization and anyonic statistics in the integer quantum Hall collider

abstract: One remarkable feature of strongly correlated systems is the phenomenon of fractionalization where quasiparticles carry only a fraction of the charge or spin of the elementary constituents. Such quasiparticles often present anyonic statistics in two dimensions and lie at the heart of the fractional quantum Hall effect. We discuss the observation of fractionalization and anyonic statistics already in the integer quantum Hall effect coupled to a metallic island. A continuous fractional emitter is proposed, which sends dilute beams of non-integer charges, and its full counting statistics is obtained. The fractional charge is governed solely by the number of ballistic channels covered by the island and it is one half of the electron charge for a single ballistic channel. We further characterize the mixing of two such fractional beams through a quantum point contact beam splitter. We predict negative cross-correlations, in strong contrast with free electrons, that depend on the double exchange phase between electrons and the fractional charges emulating anyons. The result is similar to a genuine fractional edge state as recently measured at filling $\nu = 1/3$. We revisit the physical interpretation of this experiment and point towards a direct braiding measurement rather than a deviation from fermionic antibunching.


title: Supercool Composite Dark Matter Beyond 100 TeV

abstract: Dark Matter could be a composite state of a confining sector with an approximate scale symmetry. We consider the case where the associated pseudo-Goldstone boson, the dilaton, mediates its interactions with the Standard Model. When the confining phase transition in the early universe is supercooled, its dynamics allows for Dark Matter masses up to $10^6$ TeV. We derive the precise parameter space compatible with all experimental constraints, finding that this scenario can be tested partly by telescopes and entirely by gravitational waves.


title: Thermal keV dark matter in a gauged B-L model with cosmic inflation

abstract: We investigate the possibility of keV scale thermal dark matter (DM) in a gauged $B - L$ extension of the standard model with three right-handed neutrinos (RHN) and one vector like fermion in the context of cosmic inflation. The complex singlet scalar field responsible for the spontaneous breaking of $B-L$ gauge symmetry is non-minimally coupled to gravity and serves the role of inflaton. The keV scale vector like fermion DM gives rise to the possibility of warm dark matter, but it gets overproduced thermally. The subsequent entropy dilution due to one of the RHN decay can bring the thermal abundance of DM within the observed limit. The dynamics of both the DM and the diluter are regulated by the $B-L$ model parameters which are also restricted by the requirement of successful inflationary dynamics.We constrain the model parameter space from the requirement of producing sufficient entropy dilution to obtain correct DM relic, inflationary observables along with other phenomenological constraints. Interestingly, we obtain unique predictions for the order of lightest active neutrino mass ($m_{\nu_l}\lesssim 10^{-14}$ eV) as function of $B-L$ gauge coupling for keV scale DM. The proposed framework also explains the origin of observed baryon asymmetry from the decay of other two heavier RHN by overcoming the entropy dilution effect.


title: Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer observations of the Euclid deep fields and calibration fields

abstract: We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.


title: Mapping physical parameters in Orion KL at high spatial resolution

abstract: The Orion Kleinmann-Low nebula (Orion KL) is notoriously complex and exhibits a range of physical and chemical components. We conducted high angular resolution (sub-arcsecond) observations of $^13$CH$_{3}$OH $\nu=0$ ($\sim$0.3$^{\prime\prime}$ and $\sim$0.7$^{\prime\prime}$) and CH$_3$CN $\nu_8=1$ ($\sim$0.2$^{\prime\prime}$ and $\sim$0.9$^{\prime\prime}$) line emission with the Atacama Large Millimeter/submillimeter Array (ALMA) to investigate Orion KL's structure on small spatial scales (${\le}350$ au). Gas kinematics, excitation temperatures, and column densities were derived from the molecular emission via a pixel-by-pixel spectral line fitting of the image cubes, enabling us to examine the small-scale variation of these parameters. Sub-regions of the Hot Core have a higher excitation temperature in a 0.2$^{\prime\prime}$ beam than a 0.9$^{\prime\prime}$ beam, indicative of possible internal sources of heating. Furthermore, the velocity field includes a bipolar ${\sim}7{-}8$ km s$^{-1}$ feature with a southeast-northwest orientation against the surrounding ${\sim}4{-}5$ km s$^{-1}$ velocity field, which may be due to an outflow. We also find evidence of a possible source of internal heating toward the Northwest Clump, since the excitation temperature there is higher in a smaller beam versus a larger beam. Finally, the region southwest of the Hot Core (Hot Core-SW) presents itself as a particularly heterogeneous region bridging the Hot Core and Compact Ridge. Additional studies to identify the (hidden) sources of luminosity and heating within Orion KL are necessary to better understand the nebula and its chemistry.


title: The two-loop bispectrum of large-scale structure

abstract: The bispectrum is the leading non-Gaussian statistic in large-scale structure, carrying valuable information on cosmology that is complementary to the power spectrum. To access this information, we need to model the bispectrum in the weakly non-linear regime. In this work we present the first two-loop, i.e., next-to-next-to-leading order perturbative description of the bispectrum within an effective field theory (EFT) framework. Using an analytic expansion of the perturbative kernels up to $F_6$ we derive a renormalized bispectrum that is demonstrated to be independent of the UV cutoff. We show that the EFT parameters associated with the four independent second-order EFT operators known from the one-loop bispectrum are sufficient to absorb the UV sensitivity of the two-loop contributions in the double-hard region. In addition, we employ a simplified treatment of the single-hard region, introducing one extra EFT parameter at two-loop order. We compare our results to N-body simulations using the realization-based grid-PT method and find good agreement within the expected range, as well as consistent values for the EFT parameters. The two-loop terms start to become relevant at $k\approx 0.07h~\mathrm{Mpc}^{-1}$. The range of wavenumbers with percent-level agreement, independently of the shape, extends from $0.08h~\mathrm{Mpc}^{-1}$ to $0.15h~\mathrm{Mpc}^{-1}$ when going from one to two loops at $z=0$. In addition, we quantify the impact of using exact instead of Einstein-de-Sitter kernels for the one-loop bispectrum, and discuss in how far their impact can be absorbed into a shift of the EFT parameters.


title: The unusually high dark matter concentration of the galaxy group NGC 1600

abstract: We investigate the properties of the dark matter (DM) halo surrounding the nearby galaxy group NGC 1600. Through the use of deep (252 ks) Chandra observations and 64.3 ks of XMM-Newton observations, we construct surface brightness profiles in multiple energy bands in order to perform hydrostatic equilibrium analysis of the hot plasma within NGC 1600. Regardless of the DM model profile assumed, we measure a halo concentration (c$_{200}$) that is an extreme, positive outlier of the $\Lambda$CDM c$_{200}$-M$_{200}$ relation. For a typical NFW DM profile, we measure c$_{200}\!=\!26.7\pm1.4$ and M$_{200}\!=\!(2.0\pm0.2)\times10^{13}$ M$_\odot$; assuming a similar halo mass, the average concentration expected is c$_{200}=6-7$ for the theoretical $\Lambda$CDM c-M relation. Such a high concentration is similar to that of well-known fossil groups MRK 1216 and NGC 6482. While NGC 1600 exhibits some properties of a fossil group, it fails to meet the X-ray luminosity threshold of L$_X>5\times10^{41}$ erg s$^{-1}$. Whether or not it is considered a fossil group, the high concentration value makes it part of a select group of galaxy groups.


title: Measurement of the Relativistic Sunyaev-Zeldovich Corrections in RX J1347.5-1145

abstract: We present a measurement of the relativistic corrections to the thermal Sunyaev-Zel'dovich (SZ) effect spectrum, the rSZ effect, toward the massive galaxy cluster RX J1347.5-1145 by combining sub-mm images from Herschel-SPIRE with mm-wave Bolocam maps. Our analysis simultaneously models the SZ effect signal, the population of cosmic infrared background (CIB) galaxies, and galactic cirrus dust emission in a manner that fully accounts for their spatial and frequency-dependent correlations. Gravitational lensing of background galaxies by RX J1347.5-1145 is included in our methodology based on a mass model derived from HST observations. Utilizing a set of realistic mock observations, we employ a forward modelling approach that accounts for the non-Gaussian covariances between observed astrophysical components to determine the posterior distribution of SZ effect brightness values consistent with the observed data. We determine a maximum a posteriori (MAP) value of the average Comptonization parameter of the intra-cluster medium (ICM) within R$_{2500}$ to be $\langle y \rangle_{2500} = 1.56 \times 10^{-4}$, with corresponding 68~per cent credible interval $[1.42,1.63] \times 10^{-4}$, and a MAP ICM electron temperature of $\langle \textrm{T}_{\textrm{sz}} \rangle_{2500} = 22.4$~keV with 68~per cent credible interval spanning $[10.4,33.0]$~keV. This is in good agreement with the pressure-weighted temperature obtained from {\it Chandra} X-ray observations, $\langle \textrm{T}_{\textrm{x,pw}}\rangle_{2500} = 17.4 \pm 2.3$~keV. We aim to apply this methodology to comparable existing data for a sample of 39 galaxy clusters, with an estimated uncertainty on the ensemble mean $\langle \textrm{T}_{\textrm{sz}} \rangle_{2500}$ at the $\simeq 1$~keV level, sufficiently precise to probe ICM physics and to inform X-ray temperature calibration.


title: The Cosmic Ultraviolet Baryon Survey (CUBS) IV: The Complex Multiphase Circumgalactic Medium as Revealed by Partial Lyman Limit Systems

abstract: We present a detailed study of two partial Lyman limit systems (pLLSs) of neutral hydrogen column density $N_\mathrm{H\,I}\approx(1-3)\times10^{16}\,\mathrm{cm}^{-2}$ discovered at $z=0.5$ in the Cosmic Ultraviolet Baryon Survey (CUBS). Available far-ultraviolet spectra from the Hubble Space Telescope Cosmic Origins Spectrograph and optical echelle spectra from MIKE on the Magellan Telescopes enable a comprehensive ionization analysis of diffuse circumgalactic gas based on resolved kinematics and abundance ratios of atomic species spanning five different ionization stages. These data provide unambiguous evidence of kinematically aligned multi-phase gas that masquerades as a single-phase structure and can only be resolved by simultaneous accounting of the full range of observed ionic species. Both systems are resolved into multiple components with inferred $\alpha$-element abundance varying from $[\alpha/\text{H}]\approx\!{-0.8}$ to near solar and densities spanning over two decades from $\log n_\mathrm{H}\mathrm{cm}^{-3}\approx\!-2.2$ to $<-4.3$. Available deep galaxy survey data from the CUBS program taken with VLT/MUSE, Magellan/LDSS3-C and Magellan/IMACS reveal that the $z=0.47$ system is located 55 kpc from a star-forming galaxy with prominent Balmer absorption of stellar mass $M_\star\approx2\times10^{10}M_\odot$, while the $z=0.54$ system resides in an over-dense environment of 11 galaxies within 750 kpc in projected distance, with the most massive being a luminous red galaxy of $M_\star\approx2\times10^{11}M_\odot$ at 375 kpc. The study of these two pLLSs adds to an emerging picture of the complex, multiphase circumgalactic gas that varies in chemical abundances and density on small spatial scales in diverse galaxy environments. The inhomogeneous nature of metal enrichment and density revealed in observations must be taken into account in theoretical models of diffuse halo gas.


title: Near-UV Reddening Observed in the Reflectance Spectrum of High Inclination Centaur 2012 DR30

abstract: Centaurs with high orbital inclinations and perihelia (i > 60 degrees; q >= 15 au) are a small group of poorly understood minor planets that are predicted to enter the giant planet region of the Solar System from the inner Oort Cloud. As such they are one of the few samples of relatively unaltered Oort Cloud material that can currently be directly observed. Here we present two new reflectance spectra of one of the largest of these objects, 2012 DR30, in order to constrain its color and surface composition. Contrary to reports that 2012 DR30 has variable optical color, we find that consistent measurements of its spectral gradient from most new and published datasets at 0.55-0.8 micron agree with a spectral gradient of S ~ 10+/-1 %/0.1 micron within their uncertainties. The spectral variability of 2012 DR30 at Near-UV/blue and Near-Infrared wavelengths, however, is still relatively unconstrained; self-consistent rotationally resolved followup observations are needed to characterise any spectral variation in those regions. We tentatively confirm previous detections of water ice on the surface of 2012 DR30 , and also consistently observe a steady steepening of the gradient of its spectrum from wavelengths of around 0.6 micron towards Near-UV wavelengths. Plausible surface materials responsible for the observed reddening may include ferric oxides contained within phyllosilicates, and aromatic refractory organics.


title: Frequency Centric Defense Mechanisms against Adversarial Examples

abstract: Adversarial example (AE) aims at fooling a Convolution Neural Network by introducing small perturbations in the input image.The proposed work uses the magnitude and phase of the Fourier Spectrum and the entropy of the image to defend against AE. We demonstrate the defense in two ways: by training an adversarial detector and denoising the adversarial effect. Experiments were conducted on the low-resolution CIFAR-10 and high-resolution ImageNet datasets. The adversarial detector has 99% accuracy for FGSM and PGD attacks on the CIFAR-10 dataset. However, the detection accuracy falls to 50% for sophisticated DeepFool and Carlini & Wagner attacks on ImageNet. We overcome the limitation by using autoencoder and show that 70% of AEs are correctly classified after denoising.


title: Charged particle dynamics in the vicinity of black hole from vector-tensor theory of gravity immersed in an external magnetic field

abstract: The dynamics of charged particle in the vicinity of event horizon of a rotating charged black hole (BH) with gravitomagnetic charge from a class of vector-tensor theories of modified gravity immersed in an axially symmetric magnetic field (Bfield) has been studied. The presence of Bfield reduces the radii of innermost stable circular orbit (ISCO) of charged particle and the opposite phenomenon occurs due the parameter $\beta$ which describes the deviation of modified gravity considered here from the usual Einstein - Maxwell gravity also acts to increase the radius of ISCO. Angular momentum also plays the role regarding the motion of particle which implies that larger the angular momentum easier for a particle to escape. We study the stability of orbits using Lyapunov characteristic exponent which implies more stability in the presence of Bfield. Generally, in the presence of Bfield it is easier for a particle to escape than its absence. We find twist in the trajectories of charged particles close to event horizon due to spin of BH in the presence of Bfield and there is no twist in the absence of it. Finally, we comment on possible utilization of our findings for the relativistic jets and magnetohydrodynamical out flows.


title: Provably Robust Model-Centric Explanations for Critical Decision-Making

abstract: We recommend using a model-centric, Boolean Satisfiability (SAT) formalism to obtain useful explanations of trained model behavior, different and complementary to what can be gleaned from LIME and SHAP, popular data-centric explanation tools in Artificial Intelligence (AI). We compare and contrast these methods, and show that data-centric methods may yield brittle explanations of limited practical utility. The model-centric framework, however, can offer actionable insights into risks of using AI models in practice. For critical applications of AI, split-second decision making is best informed by robust explanations that are invariant to properties of data, the capability offered by model-centric frameworks.


title: A Kpc Scale Molecular Wave in the Inner Galaxy: Feather of the Milky Way?

abstract: We report the discovery of a velocity coherent, kpc-scale molecular structure towards the Galactic center region with an angular extent of 30deg and an aspect ratio of 60:1. The kinematic distance of the CO structure ranges between 4.4 to 6.5 kpc. Analysis of the velocity data and comparison with the existing spiral arm models support that a major portion of this structure is either a sub-branch of the Norma arm or an inter-arm giant molecular filament, likely to be a kpc-scale feather (or spur) of the Milky Way, similar to those observed in nearby spiral galaxies. The filamentary cloud is at least 2.0 kpc in extent, considering the uncertainties in the kinematic distances, and it could be as long as 4 kpc. The vertical distribution of this highly elongated structure reveals a pattern similar to that of a sinusoidal wave. The exact mechanisms responsible for the origin of such a kpc-scale filament and its wavy morphology remains unclear. The distinct wave-like shape and its peculiar orientation makes this cloud, named as the Gangotri wave, one of the largest and most intriguing structures identified in the Milky Way.


title: CausalAF: Causal Autoregressive Flow for Safety-Critical Driving Scenario Generation

abstract: Generating safety-critical scenarios, which are crucial yet difficult to collect, provides an effective way to evaluate the robustness of autonomous driving systems. However, the diversity of scenarios and efficiency of generation methods are heavily restricted by the rareness and structure of safety-critical scenarios. Therefore, existing generative models that only estimate distributions from observational data are not satisfying to solve this problem. In this paper, we integrate causality as a prior into the scenario generation and propose a flow-based generative framework, Causal Autoregressive Flow (CausalAF). CausalAF encourages the generative model to uncover and follow the causal relationship among generated objects via novel causal masking operations instead of searching the sample only from observational data. By learning the cause-and-effect mechanism of how the generated scenario causes risk situations rather than just learning correlations from data, CausalAF significantly improves learning efficiency. Extensive experiments on three heterogeneous traffic scenarios illustrate that CausalAF requires much fewer optimization resources to effectively generate safety-critical scenarios. We also show that using generated scenarios as additional training samples empirically improves the robustness of autonomous driving algorithms.


title: Quantum field theory in a de Sitter universe transiting to the radiation stage

abstract: We study some physical aspects of quantum field theory in a two stage universe starting from the inflationary de Sitter and transiting into the radiation dominated stage. We look into the time evolution of the primordial vacuum states, associated with the (i) comoving and (ii) Bunch-Davies modes. We show how the power spectrum for a comoving observer, obtained from the excitation of the aforementioned states defined in the de Sitter stage, changes as the universe transits into the radiation stage. In addition, we also develop a methodology to transfer the well known result of particle creation in the static de Sitter frame, originating from the aforementioned vacuum states, while the universe makes a transition to the next (radiation dominated) stage.


title: Rapid IoT Device Identification at the Edge

abstract: Consumer Internet of Things (IoT) devices are increasingly common in everyday homes, from smart speakers to security cameras. Along with their benefits come potential privacy and security threats. To limit these threats we must implement solutions to filter IoT traffic at the edge. To this end the identification of the IoT device is the first natural step. In this paper we demonstrate a novel method of rapid IoT device identification that uses neural networks trained on device DNS traffic that can be captured from a DNS server on the local network. The method identifies devices by fitting a model to the first seconds of DNS second-level-domain traffic following their first connection. Since security and privacy threat detection often operate at a device specific level, rapid identification allows these strategies to be implemented immediately. Through a total of 51,000 rigorous automated experiments, we classify 30 consumer IoT devices from 27 different manufacturers with 82% and 93% accuracy for product type and device manufacturers respectively.


title: Vacuum laser acceleration of super-ponderomotive electrons using relativistic transparency injection

abstract: Intense lasers can accelerate electrons to very high energy over a short distance. Such compact accelerators have several potential applications including fast ignition, high energy physics, and radiography. Among the various schemes of laser-based electron acceleration, vacuum laser acceleration has the merits of super-high acceleration gradient and great simplicity. Yet its realization has been difficult because injecting free electrons into the fast-oscillating laser field is not trivial. Here we demonstrate free-electron injection and subsequent vacuum laser acceleration of electrons up to 20 MeV using the relativistic transparency effect. When a high-contrast intense laser drives a thin solid foil, electrons from the dense opaque plasma are first accelerated to near-light speed by the standing laser wave in front of the solid foil and subsequently injected into the transmitted laser field as the opaque plasma becomes relativistically transparent. It is possible to further optimize the electron injection/acceleration by manipulating the laser polarization, incident angle, and temporal pulse shaping. Our result also sheds new light on the fundamental relativistic transparency process, crucial for producing secondary particle and light sources.


title: Precise distances from OGLE-IV member RR Lyrae stars in six bulge globular clusters

abstract: Context. RR Lyrae stars are useful standard candles allowing one to derive accurate distances for old star clusters. Based on the recent catalogues from OGLE-IV and Gaia Early Data Release 3 (EDR3), the distances can be improved for a few bulge globular clusters. Aims. The aim of this work is to derive an accurate distance for the following six moderately metal-poor, relatively high-reddening bulge globular clusters: NGC 6266, NGC 6441, NGC 6626, NGC 6638, NGC 6642, and NGC 6717. Methods. We combined newly available OGLE-IV catalogues of variable stars containing mean I magnitudes, with Clement's previous catalogues containing mean V magnitudes, and with precise proper motions from Gaia EDR3. Astrometric membership probabilities were computed for each RR Lyrae, in order to select those compatible with the cluster proper motions. Applying luminosity-metallicity relations derived from BaSTI $\alpha$-enhanced models (He-enhanced for NGC 6441 and canonical He for the other clusters), we updated the distances with relatively low uncertainties. Results. Distances were derived with the I and V bands, with a $5-8\%$ precision. We obtained 6.6 kpc, 13.1 kpc, 5.6 kpc, 9.6 kpc, 8.2 kpc, and 7.3 kpc for NGC 6266, NGC 6441, NGC 6626, NGC 6638, NGC 6642, and NGC 6717, respectively. The results are in excellent agreement with the literature for all sample clusters, considering the uncertainties. Conclusions. The present method of distance derivation, based on recent data of member RR Lyrae stars, updated BaSTI models, and robust statistical methods, proved to be consistent. A larger sample of clusters will be investigated in a future work.


title: Stochastic Low Frequency Variability in 3-Dimensional Radiation Hydrodynamical Models of Massive Star Envelopes

abstract: Increasing main sequence stellar luminosity with stellar mass leads to the eventual dominance of radiation pressure in stellar envelope hydrostatic balance. As the luminosity approaches the Eddington limit, additional instabilities (beyond conventional convection) can occur. These instabilities readily manifest in the outer envelopes of OB stars, where the opacity increase associated with iron yields density and gas pressure inversions in 1D models. Additionally, recent photometric surveys (e.g. TESS) have detected excess broadband low frequency variability in power spectra of OB star lightcurves, called stochastic low frequency variability (SLFV). This motivates our novel 3D Athena++ radiation hydrodynamical (RHD) simulations of two 35$\,$M$_\odot$ star envelopes (the outer $\approx$15$\%$ of the stellar radial extent), one on the zero-age main sequence and the other in the middle of the main sequence. Both models exhibit turbulent motion far above and below the conventional iron opacity peak convection zone (FeCZ), obliterating any ``quiet" part of the near-surface region and leading to velocities at the photosphere of 10-100$\,$km$\,$s$^{-1}$, directly agreeing with spectroscopic data. Surface turbulence also produces SLFV in model lightcurves with amplitudes and power-law slopes that are strikingly similar to those of observed stars. The characteristic frequencies associated with SLFV in our models are comparable to the thermal time in the FeCZ ($\approx$3-7$\,$days$^{-1}$). These simulations, which have no free parameters, are directly validated by observations and, though more models are needed, we remain optimistic that 3D RHD models of main sequence O star envelopes exhibit SLFV originating from the FeCZ.


title: On a range of exponents for absence of Lavrentiev phenomenon for double phase functionals

abstract: For a class of functionals having the $(p,q)$-growth, we establish an improved range of exponents $p$, $q$ for which the Lavrentiev phenomenon does not occur. The proof is based on a standard mollification argument and Young convolution inequality. Our contribution is two-fold. First, we observe that it is sufficient to regularise only bounded functions. Second, we exploit the $L^{\infty}$ bound on the function rather than the $L^p$ estimate on the gradient. Our proof does not rely on the properties of minimizers to variational problems but it is rather a consequence of the underlying Musielak-Orlicz function spaces. Moreover, our method works for unbounded boundary data, the variable exponent functionals and vectorial problems. In addition, the result seems to be optimal for $p\le d$.


title: Mixed-state TQFTs

abstract: In this short note, we propose a generalization of Atiyah type TQFTs from pure states to mixed states in the sense that the Hilbert space of pure states associated to a space manifold is replaced by a quantum coherent space related to density matrices. Atiyah type TQFT is a symmetric monoidal functor from the Bord category of manifolds to the category Vec of finite dimensional vector spaces. In this paper, we define mixed-state TQFTs by replacing the target category Vec by QCS--the category of quantum coherent spaces, then a mixed-state TQFT is simply a symmetric monoidal functor from Bord to QCS. We also discuss how to construct interesting examples from subsystem quantum error correction codes beyond the trivial ones--all Atiyah type TQFTs.


title: Collaborative Uncertainty in Multi-Agent Trajectory Forecasting

abstract: Uncertainty modeling is critical in trajectory forecasting systems for both interpretation and safety reasons. To better predict the future trajectories of multiple agents, recent works have introduced interaction modules to capture interactions among agents. This approach leads to correlations among the predicted trajectories. However, the uncertainty brought by such correlations is neglected. To fill this gap, we propose a novel concept, collaborative uncertainty(CU), which models the uncertainty resulting from the interaction module. We build a general CU-based framework to make a prediction model to learn the future trajectory and the corresponding uncertainty. The CU-based framework is integrated as a plugin module to current state-of-the-art (SOTA) systems and deployed in two special cases based on multivariate Gaussian and Laplace distributions. In each case, we conduct extensive experiments on two synthetic datasets and two public, large-scale benchmarks of trajectory forecasting. The results are promising: 1) The results of synthetic datasets show that CU-based framework allows the model to appropriately approximate the ground-truth distribution. 2) The results of trajectory forecasting benchmarks demonstrate that the CU-based framework steadily helps SOTA systems improve their performances. Especially, the proposed CU-based framework helps VectorNet improve by 57cm regarding Final Displacement Error on nuScenes dataset. 3) The visualization results of CU illustrate that the value of CU is highly related to the amount of the interactive information among agents.


title: Boosted CVaR Classification

abstract: Many modern machine learning tasks require models with high tail performance, i.e. high performance over the worst-off samples in the dataset. This problem has been widely studied in fields such as algorithmic fairness, class imbalance, and risk-sensitive decision making. A popular approach to maximize the model's tail performance is to minimize the CVaR (Conditional Value at Risk) loss, which computes the average risk over the tails of the loss. However, for classification tasks where models are evaluated by the zero-one loss, we show that if the classifiers are deterministic, then the minimizer of the average zero-one loss also minimizes the CVaR zero-one loss, suggesting that CVaR loss minimization is not helpful without additional assumptions. We circumvent this negative result by minimizing the CVaR loss over randomized classifiers, for which the minimizers of the average zero-one loss and the CVaR zero-one loss are no longer the same, so minimizing the latter can lead to better tail performance. To learn such randomized classifiers, we propose the Boosted CVaR Classification framework which is motivated by a direct relationship between CVaR and a classical boosting algorithm called LPBoost. Based on this framework, we design an algorithm called $\alpha$-AdaLPBoost. We empirically evaluate our proposed algorithm on four benchmark datasets and show that it achieves higher tail performance than deterministic model training methods.


title: Deletion-contraction for a unified Laplacian and applications

abstract: We define a graph Laplacian with vertex weights in addition to the more classical edge weights, which unifies the combinatorial Laplacian and the normalised Laplacian. Moreover, we give a combinatorial interpretation for the coefficients of the weighted Laplacian characteristic polynomial in terms of weighted spanning forests and use this to prove a deletion-contraction relation. We prove various interlacing theorems relating to deletion and contraction, as well as to rectangular tilings, drawing on the work of Brooks, Smith, Stone and Tutte on square tilings. Additionally, we show that the weighted Laplacian also satisfies a vertex analogue of deletion-contraction. We give applications of weighted Laplacian eigenvalues to sparse cuts, independent sets and graph colouring, and establish new cases of a conjecture of Stanley on distinguishing nonisomorphic trees.


title: Can't Fool Me: Adversarially Robust Transformer for Video Understanding

abstract: Deep neural networks have been shown to perform poorly on adversarial examples. To address this, several techniques have been proposed to increase robustness of a model for image classification tasks. However, in video understanding tasks, developing adversarially robust models is still unexplored. In this paper, we aim to bridge this gap. We first show that simple extensions of image based adversarially robust models slightly improve the worst-case performance. Further, we propose a temporal attention regularization scheme in Transformer to improve the robustness of attention modules to adversarial examples. We illustrate using a large-scale video data set YouTube-8M that the final model (A-ART) achieves close to non-adversarial performance on its adversarial example set. We achieve 91% GAP on adversarial examples, whereas baseline Transformer and simple adversarial extensions achieve 72.9% and 82% respectively, showing significant improvement in robustness over the state-of-the-art.


title: Open microscopy in the life sciences: Quo Vadis?

abstract: Light microscopy allows observing cellular features and objects with sub-micrometer resolution. As such, light microscopy has been playing a fundamental role in the life sciences for more than a hundred years. Fueled by the availability of mass-produced electronics and hardware, publicly shared documentation and building instructions, open-source software, wide access to rapid prototyping and 3D printing, and the enthusiasm of contributors and users involved, the concept of open microscopy has been gaining incredible momentum, bringing new sophisticated tools to an expanding user base. Here, we will first discuss the ideas behind open science and open microscopy before highlighting recent projects and developments in open microscopy. We argue that the availability of well-designed open hardware and software solutions targeting broad user groups or even non-experts, will increasingly be relevant to cope with the increasing complexity of cutting-edge imaging technologies. We will then extensively discuss the current and future challenges of open microscopy.

categories: q-bio.OT physics.bio-ph physics.data-an physics.ins-det physics.optics


title: Changes in the near-surface shear layer of the Sun

abstract: We use helioseismic data obtained over two solar cycles to determine whether there are changes in the near-surface shear layer (NSSL). We examine this by determining the radial gradient of the solar rotation rate. The radial gradient itself shows a solar-cycle dependence, and the changes are more pronounced in the active latitudes than at adjoining higher latitudes; results at the highest latitudes (greater than about70 degrees) are unreliable. The pattern changes with depth, even within the NSSL. We find that the near-surface shear layer is deeper at lower latitudes than at high latitudes and that the extent of the layer also shows a small solar-cycle related change.


title: On sensitivity of meta-learning to support data

abstract: Meta-learning algorithms are widely used for few-shot learning. For example, image recognition systems that readily adapt to unseen classes after seeing only a few labeled examples. Despite their success, we show that modern meta-learning algorithms are extremely sensitive to the data used for adaptation, i.e. support data. In particular, we demonstrate the existence of (unaltered, in-distribution, natural) images that, when used for adaptation, yield accuracy as low as 4\% or as high as 95\% on standard few-shot image classification benchmarks. We explain our empirical findings in terms of class margins, which in turn suggests that robust and safe meta-learning requires larger margins than supervised learning.


title: Haefliger's approach for spherical knots modulo immersions

abstract: We show that for the spaces of spherical embeddings modulo immersions $\bar{Emb}(S^n,S^{n+q})$ and long embeddings modulo immersions $\bar{Emb}_{\partial}(D^n,D^{n+q})$, the set of connected components is isomorphic to $\pi_{n+1}(SG,SG_q)$ for $q\geq 3$. As a consequence, we show that all the terms of the long exact sequence of the triad $(SG;SO,SG_q)$ have a geometric meaning relating to spherical embeddings and immersions.


title: Energetics of the Charge Generation in Organic Donor-Acceptor Interfaces

abstract: Non-fullerene acceptor (NFA) materials have posed new paradigms for the design of organic solar cells (OSC), whereby efficient carrier generation is obtained with small driving forces, in order to maximize the open-circuit voltage. In this paper we use a coarse-grained mixed quantum-classical method, that combines Ehrenfest and Redfield theories, to shed light on charge generation process in small energy offset interfaces. We have investigated the influence of the energetic driving force as well as the vibronic effects on the charge generation and photovoltaic energy conversion. By analyzing the effects of the Holstein and Peierls vibrational couplings, we find that vibrational couplings produce an overall effect of improving the charge generation. However, the two vibronic mechanisms play different roles: the Holstein relaxation mechanism decreases the charge generation whereas the Peierls mechanism always assists the charge generation. Moreover, by examining the electron-hole binding energy as a function of time, we evince two distinct regimes for the charge separation: the temperature independent excitonic spread on a sub-100 fs timescale and the complete dissociation of the charge-transfer state that occurs on the timescale of tens to hundreds of picoseconds, depending on the temperature. The quantum dynamics of the system exhibits the three regimes of the Marcus electron transfer kinetics as the energy offset of the interface is varied.


title: Modelling of spatial infection spread through heterogeneous population: from lattice to PDE models

abstract: We present a simple model for the spread of an infection that incorporates spatial variability in population density. Starting from first principle considerations, we explore how a novel PDE with state-dependent diffusion can be obtained. This model exhibits higher infection rates in the areas of higher population density, a feature that we argue to be consistent with epidemiological observations. The model also exhibits an infection wave whose speed varies with population density. In addition, we demonstrate the possibility that an infection can ``jump'' (i.e., tunnel) across areas of low population density towards the areas of high population density. We briefly touch upon the data reported for coronavirus spread in the Canadian province of Nova Scotia as a case example with a number of qualitatively similar features as our model. Lastly, we propose a number of generalizations of the model towards future studies.


title: Unbiased Graph Embedding with Biased Graph Observations

abstract: Graph embedding techniques are pivotal in real-world machine learning tasks that operate on graph-structured data, such as social recommendation and protein structure modeling. Embeddings are mostly performed on the node level for learning representations of each node. Since the formation of a graph is inevitably affected by certain sensitive node attributes, the node embeddings can inherit such sensitive information and introduce undesirable biases in downstream tasks. Most existing works impose ad-hoc constraints on the node embeddings to restrict their distributions for unbiasedness/fairness, which however compromise the utility of the resulting embeddings. In this paper, we propose a principled new way for unbiased graph embedding by learning node embeddings from an underlying bias-free graph, which is not influenced by sensitive node attributes. Motivated by this new perspective, we propose two complementary methods for uncovering such an underlying graph, with the goal of introducing minimum impact on the utility of the embeddings. Both our theoretical justification and extensive experimental comparisons against state-of-the-art solutions demonstrate the effectiveness of our proposed methods.


title: High-Energy Neutrino Production in Clusters of Galaxies

abstract: In this work, we compute the contribution from clusters of galaxies to the diffuse neutrino background. Clusters of galaxies can potentially produce cosmic rays (CRs) up to very-high energies via large-scale shocks and turbulent acceleration. Due to their unique magnetic-field configuration, CRs with energy $\leq 10^{17}$ eV can be trapped within these structures over cosmological time scales, and generate secondary particles, including neutrinos and gamma rays, through interactions with the background gas and photons. We employ three-dimensional cosmological magnetohydrodynamical simulations of structure formation to model the turbulent intergalactic medium. We use the distribution of clusters within this cosmological volume to extract the properties of this population. We propagate CRs in this environment using multi-dimensional Monte Carlo simulations across different redshifts (from $z \sim 5 \; \text{to} \; z =0$), considering all relevant photohadronic, photonuclear, and hadronuclear interactions. We also include the cosmological evolution of the CR sources. We find that, for CRs injected with a spectral index $1.5 - 2.7$ and cutoff energy $E_{max} = 10^{16} - 10^{17}$ eV, clusters contribute to a substantial fraction to the diffuse flux observed by the IceCube Neutrino Observatory, and most of the contribution comes from clusters with $M > 10^{14} \; M_{\odot}$ and redshift $z < 0.3$.


title: The Integrated Sachs Wolfe effect: unWISE and Planck constraints on Dynamical Dark Energy

abstract: CMB photons redshift and blueshift as they move through gravitational potentials $\Phi$ while propagating across the Universe. If the potential is not constant in time, the photons will pick up a net redshift or blueshift, known as the Integrated Sachs-Wolfe (ISW) effect. In the $z \ll 1000$ universe, $\dot{\Phi}$ is nonzero on large scales when the Universe transitions from matter to dark energy domination. This effect is only detectable in cross-correlation with large-scale structure at $z \sim 1$. In this paper we present a 3.2$\sigma$ detection of the ISW effect using cross-correlations between unWISE infrared galaxies and Planck CMB temperature maps. We use 3 tomographic galaxy samples spanning $0 < z < 2$, allowing us to fully probe the dark energy domination era and the transition into matter domination. This measurement is consistent with $\Lambda$CDM ($A_{\rm ISW} = 0.96 \pm 0.30$). We study constraints on a particular class of dynamical dark energy models (where the dark energy equation of state is different in matter and dark energy domination), finding that unWISE-ISW improves constraints from type Ia supernovae due to improved constraints on the time evolution of dark energy. When combining with BAO measurements, we obtain the tightest constraints on specific dynamical dark energy models. In the context of a phenomenological model for freezing quintessence, the Mocker model, we constrain the dark energy density within 10% at $z < 2$ using ISW, BAO and supernovae. Moreover, the ISW measurement itself provides an important independent check when relaxing assumptions about the theory of gravity, as it is sensitive to the gravitational potential rather than the expansion history.


title: Probing hadronic interactions with measurements from the Pierre Auger Observatory

abstract: The Pierre Auger Observatory is the largest facility in the world to study ultra-high-energy cosmic rays. It has a hybrid detection technique that combines the observation of the longitudinal development of extensive air showers and the measurement of their particles at the ground. This capability has opened the possibility to probe hadronic interactions taking place at energies well beyond those accessible by human-made accelerators. In this report, we present a selection of the latest results on hadronic interactions with measurements from the Pierre Auger Observatory. These data span over three decades in energy, showing the tension between data from the muon component of air showers and predictions based on the most updated hadronic interaction models.


title: Wire-Cell 3D Pattern Recognition Techniques for Neutrino Event Reconstruction in Large LArTPCs: Algorithm Description and Quantitative Evaluation with MicroBooNE Simulation

abstract: Wire-Cell is a 3D event reconstruction package for liquid argon time projection chambers. Through geometry, time, and drifted charge from multiple readout wire planes, 3D space points with associated charge are reconstructed prior to the pattern recognition stage. Pattern recognition techniques, including track trajectory and $dQ/dx$ (ionization charge per unit length) fitting, 3D neutrino vertex fitting, track and shower separation, particle-level clustering, and particle identification are then applied on these 3D space points as well as the original 2D projection measurements. A deep neural network is developed to enhance the reconstruction of the neutrino interaction vertex. Compared to traditional algorithms, the deep neural network boosts the vertex efficiency by a relative 30\% for charged-current $\nu_e$ interactions. This pattern recognition achieves 80-90\% reconstruction efficiencies for primary leptons, after a 65.8\% (72.9\%) vertex efficiency for charged-current $\nu_e$ ($\nu_\mu$) interactions. Based on the resulting reconstructed particles and their kinematics, we also achieve 15-20\% energy reconstruction resolutions for charged-current neutrino interactions.


title: The development of a split-tail heliosphere and the role of non-ideal processes: a comparison of the BU and Moscow models

abstract: Global models of the heliosphere are critical tools used in the interpretation of heliospheric observations. There are several three-dimensional magnetohydrodynamic (MHD) heliospheric models that rely on different strategies and assumptions. Until now only one paper has compared global heliosphere models, but without magnetic field effects. We compare the results of two different MHD models, the BU and Moscow models. Both models use identical boundary conditions to compare how different numerical approaches and physical assumptions contribute to the heliospheric solution. Based on the different numerical treatments of discontinuities, the BU model allows for the presence of magnetic reconnection, while the Moscow model does not. Both models predict collimation of the solar outflow in the heliosheath by the solar magnetic field and produce a split-tail where the solar magnetic field confines the charged solar particles into distinct north and south columns that become lobes. In the BU model, the ISM flows between the two lobes at large distances due to MHD instabilities and reconnection. Reconnection in the BU model at the port flank affects the draping of the interstellar magnetic field in the immediate vicinity of the heliopause. Different draping in the models cause different ISM pressures, yielding different heliosheath thicknesses and boundary locations, with the largest effects at high latitudes. The BU model heliosheath is 15% thinner and the heliopause is 7% more inwards at the north pole relative to the Moscow model. These differences in the two plasma solutions may manifest themselves in energetic neutral atom measurements of the heliosphere.


title: Poitou-Tate Duality for totally positive Galois cohomology

abstract: In this paper, we establish a Poitou-Tate's global duality for totally positive Galois cohomology. We illustrate this result in the case of the twisted module "\`a la Tate" $\mathbb{Z}_{2}(i)$, $i$ integer.


title: Internally-driven warps in disc galaxies

abstract: Any perturbation to a disc galaxy that creates a misalignment between the planes of the inner and outer disc, will excite a slowly evolving bending wave in the outer disc. The torque from the stiff inner disc drives a retrograde, leading-spiral bending wave that grows in amplitude as it propagates outward over a period of several Gyr. This behaviour creates warps that obey the rules established from observations, and operates no matter what the original cause of the misalignment between the inner and outer disc. The part of the disc left behind by the outwardly propagating wave is brought into alignment with the inner disc. Here we confirm that mild warps in simulations of disc galaxies can be excited by shot noise in the halo, as was recently reported. We show that the quadrupole component of the noise creates disc distortions most effectively. Bending waves caused by shot noise in carefully constructed equilibrium simulations of isolated galaxies are far too mild to be observable, but perturbations from halo substructure and galaxy assembly must excite larger amplitude bending waves in real galaxies.


title: Signature of a heliotail organized by the solar magnetic field and the role of non-ideal processes in modeled IBEX ENA maps: a comparison of the BU and Moscow MHD models

abstract: Energetic neutral atom (ENA) models typically require post-processing routines to convert the distributions of plasma and H atoms into ENA maps. Here we investigate how two different kinetic-MHD models of the heliosphere (the BU and Moscow models) manifest in modeled ENA maps using the same prescription and how they compare with Interstellar Boundary Explorer (IBEX) observations. Both MHD models treat the solar wind as a single-ion plasma for protons, which include thermal solar wind ions, pick-up ions (PUIs), and electrons. Our ENA prescription partitions the plasma into three distinct ion populations (thermal solar wind, PUIs transmitted and ones energized at the termination shock) and models the populations with Maxwellian distributions. Both kinetic-MHD heliospheric models produce a heliotail with heliosheath plasma organized by the solar magnetic field into two distinct north and south columns that become lobes of high mass flux flowing down the heliotail, though in the BU model the ISM flows between the two lobes at distances in the heliotail larger than 300 AU. While our prescription produces similar ENA maps for the two different plasma and H atom solutions at the IBEX-Hi energy range (0.5 - 6 keV), the modeled ENA maps require a scaling factor of ~2 to be in agreement with the data. This problem is present in other ENA models with the Maxwellian approximation of multiple ion species and indicates that a higher neutral density or some acceleration of PUIs in the heliosheath is required.


title: Fisheries Management in Congested Waters: A Game-Theoretic Assessment of the East China Sea

abstract: Fisheries in the East China Sea (ECS) face multiple concerning trends. Aside from depleted stocks caused by overfishing, illegal encroachments by fishermen from one nation into another's legal waters are a common occurrence. This behavior presumably could be stopped via strong monitoring, controls, and surveillance (MCS), but MCS is routinely rated below standards for nations bordering the ECS. This paper generalizes the ECS to a model of a congested maritime environment, defined as an environment where multiple nations can fish in the same waters with equivalent operating costs, and uses game-theoretic analysis to explain why the observed behavior persists in the ECS. The paper finds that nations in congested environments are incentivized to issue excessive quotas, which in turn tacitly encourages illegal fishing and extracts illegal rent from another's legal waters. This behavior couldn't persist in the face of strong MCS measures, and states are thus likewise incentivized to use poor MCS. A bargaining problem is analyzed to complement the noncooperative game, and a key finding is the nation with lower nonoperating costs has great leverage during the bargain.


title: Evaluating Serverless Architecture for Big Data Enterprise Applications

abstract: In this paper, we investigate serverless computing for performing large scale data processing with cloudnative primitives.


title: On the Effects of Artificial Data Modification

abstract: Data distortion is commonly applied in vision models during both training (e.g methods like MixUp and CutMix) and evaluation (e.g. shape-texture bias and robustness). This data modification can introduce artificial information. It is often assumed that the resulting artefacts are detrimental to training, whilst being negligible when analysing models. We investigate these assumptions and conclude that in some cases they are unfounded and lead to incorrect results. Specifically, we show current shape bias identification methods and occlusion robustness measures are biased and propose a fairer alternative for the latter. Subsequently, through a series of experiments we seek to correct and strengthen the community's perception of how augmenting affects learning of vision models. Based on our empirical results we argue that the impact of the artefacts must be understood and exploited rather than eliminated.


title: Nonparametric Matrix Estimation with One-Sided Covariates

abstract: Consider the task of matrix estimation in which a dataset $X \in \mathbb{R}^{n\times m}$ is observed with sparsity $p$, and we would like to estimate $\mathbb{E}[X]$, where $\mathbb{E}[X_{ui}] = f(\alpha_u, \beta_i)$ for some Holder smooth function $f$. We consider the setting where the row covariates $\alpha$ are unobserved yet the column covariates $\beta$ are observed. We provide an algorithm and accompanying analysis which shows that our algorithm improves upon naively estimating each row separately when the number of rows is not too small. Furthermore when the matrix is moderately proportioned, our algorithm achieves the minimax optimal nonparametric rate of an oracle algorithm that knows the row covariates. In simulated experiments we show our algorithm outperforms other baselines in low data regimes.


title: Rademacher Random Projections with Tensor Networks

abstract: Random projection (RP) have recently emerged as popular techniques in the machine learning community for their ability in reducing the dimension of very high-dimensional tensors. Following the work in [30], we consider a tensorized random projection relying on Tensor Train (TT) decomposition where each element of the core tensors is drawn from a Rademacher distribution. Our theoretical results reveal that the Gaussian low-rank tensor represented in compressed form in TT format in [30] can be replaced by a TT tensor with core elements drawn from a Rademacher distribution with the same embedding size. Experiments on synthetic data demonstrate that tensorized Rademacher RP can outperform the tensorized Gaussian RP studied in [30]. In addition, we show both theoretically and experimentally, that the tensorized RP in the Matrix Product Operator (MPO) format is not a Johnson-Lindenstrauss transform (JLT) and therefore not a well-suited random projection map


title: Diachronic Text Mining Investigation of Therapeutic Candidates for COVID-19

abstract: Diachronic text mining has frequently been applied to long-term linguistic surveys of word meaning and usage shifts over time. In this paper we apply short-term diachronic text mining to a rapidly growing corpus of scientific publications on COVID-19 captured in the CORD-19 dataset in order to identify co-occurrences and analyze the behavior of potential candidate treatments. We used a data set associated with a COVID-19 drug re-purposing study from Oak Ridge National Laboratory. This study identified existing candidate coronavirus treatments, including drugs and approved compounds, which had been analyzed and ranked according to their potential for blocking the ability of the SARS-COV-2 virus to invade human cells. We investigated the occurrence of these candidates in temporal instances of the CORD-19 corpus. We found that at least 25% of the identified terms occurred in temporal instances of the corpus to the extent that their frequency and contextual dynamics could be evaluated. We identified three classes of behaviors: those where frequency and contextual shifts were small and positively correlated; those where there was no correlation between frequency and contextual changes; and those where there was a negative correlation between frequency and contextual shift. We speculate that the latter two patterns are indicative that a target candidate therapeutics is undergoing active evaluation. The patterns we detected demonstrate the potential benefits of using diachronic text mining techniques with a large dynamic text corpus to track drug-repurposing activities across international clinical and laboratory settings.


title: Video-based fully automatic assessment of open surgery suturing skills

abstract: The goal of this study was to develop new reliable open surgery suturing simulation system for training medical students in situation where resources are limited or in the domestic setup. Namely, we developed an algorithm for tools and hands localization as well as identifying the interactions between them based on simple webcam video data, calculating motion metrics for assessment of surgical skill. Twenty-five participants performed multiple suturing tasks using our simulator. The YOLO network has been modified to a multi-task network, for the purpose of tool localization and tool-hand interaction detection. This was accomplished by splitting the YOLO detection heads so that they supported both tasks with minimal addition to computer run-time. Furthermore, based on the outcome of the system, motion metrics were calculated. These metrics included traditional metrics such as time and path length as well as new metrics assessing the technique participants use for holding the tools. The dual-task network performance was similar to that of two networks, while computational load was only slightly bigger than one network. In addition, the motion metrics showed significant differences between experts and novices. While video capture is an essential part of minimally invasive surgery, it is not an integral component of open surgery. Thus, new algorithms, focusing on the unique challenges open surgery videos present, are required. In this study, a dual-task network was developed to solve both a localization task and a hand-tool interaction task. The dual network may be easily expanded to a multi-task network, which may be useful for images with multiple layers and for evaluating the interaction between these different layers.


title: The Value of Information When Deciding What to Learn

abstract: All sequential decision-making agents explore so as to acquire knowledge about a particular target. It is often the responsibility of the agent designer to construct this target which, in rich and complex environments, constitutes a onerous burden; without full knowledge of the environment itself, a designer may forge a sub-optimal learning target that poorly balances the amount of information an agent must acquire to identify the target against the target's associated performance shortfall. While recent work has developed a connection between learning targets and rate-distortion theory to address this challenge and empower agents that decide what to learn in an automated fashion, the proposed algorithm does not optimally tackle the equally important challenge of efficient information acquisition. In this work, building upon the seminal design principle of information-directed sampling (Russo & Van Roy, 2014), we address this shortcoming directly to couple optimal information acquisition with the optimal design of learning targets. Along the way, we offer new insights into learning targets from the literature on rate-distortion theory before turning to empirical results that confirm the value of information when deciding what to learn.


title: Global sensitivity analysis of rare event probabilities

abstract: By their very nature, rare event probabilities are expensive to compute; they are also delicate to estimate as their value strongly depends on distributional assumptions on the model parameters. Hence, understanding the sensitivity of the computed rare event probabilities to the hyper-parameters that define the distribution law of the model parameters is crucial. We show that by (i) accelerating the calculation of rare event probabilities through subset simulation and (ii) approximating the resulting probabilities through a polynomial chaos expansion, the global sensitivity of such problems can be analyzed through a double-loop sampling approach. The resulting method is conceptually simple and computationally efficient; its performance is illustrated on a subsurface flow application and on an analytical example.


title: Multistationarity in cyclic sequestration-transmutation networks

abstract: We consider a natural class of reaction networks which consist of reactions where either two species can inactivate each other (i.e., sequestration), or some species can be transformed into another (i.e., transmutation), in a way that gives rise to a feedback cycle. We completely characterize the capacity of multistationarity of these networks. This is especially interesting because such networks provide simple examples of "atoms of multistationarity", i.e., minimal networks that can give rise to multiple positive steady states


title: Biological learning in key-value memory networks

abstract: In neuroscience, classical Hopfield networks are the standard biologically plausible model of long-term memory, relying on Hebbian plasticity for storage and attractor dynamics for recall. In contrast, memory-augmented neural networks in machine learning commonly use a key-value mechanism to store and read out memories in a single step. Such augmented networks achieve impressive feats of memory compared to traditional variants, yet their biological relevance is unclear. We propose an implementation of basic key-value memory that stores inputs using a combination of biologically plausible three-factor plasticity rules. The same rules are recovered when network parameters are meta-learned. Our network performs on par with classical Hopfield networks on autoassociative memory tasks and can be naturally extended to continual recall, heteroassociative memory, and sequence learning. Our results suggest a compelling alternative to the classical Hopfield network as a model of biological long-term memory.


title: On the $L_\infty$ formulation of Chern-Simons theories

abstract: $L_{\infty}$ algebras have been recently studied as algebraic frameworks in the formulation of gauge theories in which the gauge symmetries and the dynamics of the interacting theories are contained in a set of products acting on a graded vector space. On the other hand, FDAs are differential algebras that generalize Lie algebras by including higher-degree differential forms on their differential equations. In this article, we review the dual relation between FDAs and $L_{\infty}$ algebras. We study the formulation of standard Chern--Simons theories in terms of $L_{\infty}$ algebras and extend the results to FDA-based gauge theories. We focus on two cases, namely a flat (or zero-curvature) theory and a generalized Chern--Simons theory, both including high-degree differential forms as fundamental fields.


title: Search for an anomalous excess of inclusive charged-current $\nu_e$ interactions in the MicroBooNE experiment using Wire-Cell reconstruction

abstract: We report a search for an anomalous excess of inclusive charged-current (CC) $\nu_e$ interactions using the Wire-Cell event reconstruction package in the MicroBooNE experiment, which is motivated by the previous observation of a low-energy excess (LEE) of electromagnetic events from the MiniBooNE experiment. With a single liquid argon time projection chamber detector, the measurements of $\nu_{\mu}$ CC interactions as well as $\pi^0$ interactions are used to constrain signal and background predictions of $\nu_e$ CC interactions. A data set collected from February 2016 to July 2018 corresponding to an exposure of 6.369 $\times$ 10$^{20}$ protons on target from the Booster Neutrino Beam at FNAL is analyzed. With $x$ representing an overall normalization factor and referred to as the LEE strength parameter, we select 56 fully contained $\nu_e$ CC candidates while expecting 69.6 $\pm$ 8.0 (stat.) $\pm$ 5.0 (sys.) and 103.8 $\pm$ 9.0 (stat.) $\pm$ 7.4 (sys.) candidates after constraints for the absence (eLEE$_{x=0}$) of the median signal strength derived from the MiniBooNE observation and the presence (eLEE$_{x=1}$) of that signal strength, respectively. Under a nested hypothesis test using both rate and shape information in all available channels, the best-fit $x$ is determined to be 0 (eLEE$_{x=0}$) with a 95.5% confidence level upper limit of $x$ at 0.502. Under a simple-vs-simple hypotheses test, the eLEE$_{x=1}$ hypothesis is rejected at 3.75$\sigma$, while the eLEE$_{x=0}$ hypothesis is shown to be consistent with the observation at 0.45$\sigma$. In the context of the eLEE model, the estimated 68.3% confidence interval of the $\nu_e$ hypothesis to explain the LEE observed in the MiniBooNE experiment is disfavored at a significance level of more than 2.6$\sigma$ (3.0$\sigma$) considering MiniBooNE's full (statistical) uncertainties.


title: Fast Asymptotic-Numerical Method For Coarse Mesh Particle Simulation In Channels Of Arbitrary Cross Section

abstract: Particles traveling through inertial microfluidic devices migrate to focusing streamlines. We present a numerical method that calculates migration velocities of particles in inertial microfluidic channels of arbitrary cross section by representing particles by singularities. Refinements to asymptotic analysis are given that improve the regularity of the singularity approximation, making finite element approximations of flow and pressure fields more effective. Sample results demonstrate that the method is computationally efficient and able to capture bifurcations in particle focusing positions due to changes in channel shape and Reynolds number.


title: Adversarial Attacks and Defenses for Social Network Text Processing Applications: Techniques, Challenges and Future Research Directions

abstract: The growing use of social media has led to the development of several Machine Learning (ML) and Natural Language Processing(NLP) tools to process the unprecedented amount of social media content to make actionable decisions. However, these MLand NLP algorithms have been widely shown to be vulnerable to adversarial attacks. These vulnerabilities allow adversaries to launch a diversified set of adversarial attacks on these algorithms in different applications of social media text processing. In this paper, we provide a comprehensive review of the main approaches for adversarial attacks and defenses in the context of social media applications with a particular focus on key challenges and future research directions. In detail, we cover literature on six key applications, namely (i) rumors detection, (ii) satires detection, (iii) clickbait & spams identification, (iv) hate speech detection, (v)misinformation detection, and (vi) sentiment analysis. We then highlight the concurrent and anticipated future research questions and provide recommendations and directions for future work.


title: CHIP: CHannel Independence-based Pruning for Compact Neural Networks

abstract: Filter pruning has been widely used for neural network compression because of its enabled practical acceleration. To date, most of the existing filter pruning works explore the importance of filters via using intra-channel information. In this paper, starting from an inter-channel perspective, we propose to perform efficient filter pruning using Channel Independence, a metric that measures the correlations among different feature maps. The less independent feature map is interpreted as containing less useful information$/$knowledge, and hence its corresponding filter can be pruned without affecting model capacity. We systematically investigate the quantification metric, measuring scheme and sensitiveness$/$reliability of channel independence in the context of filter pruning. Our evaluation results for different models on various datasets show the superior performance of our approach. Notably, on CIFAR-10 dataset our solution can bring $0.90\%$ and $0.94\%$ accuracy increase over baseline ResNet-56 and ResNet-110 models, respectively, and meanwhile the model size and FLOPs are reduced by $42.8\%$ and $47.4\%$ (for ResNet-56) and $48.3\%$ and $52.1\%$ (for ResNet-110), respectively. On ImageNet dataset, our approach can achieve $40.8\%$ and $44.8\%$ storage and computation reductions, respectively, with $0.15\%$ accuracy increase over the baseline ResNet-50 model. The code is available at https://github.com/Eclipsess/CHIP_NeurIPS2021.


title: Global well-posedness for a system of quasilinear wave equations on a product space

abstract: We consider a system of quasilinear wave equations on the product space $\mathbb{R}^{1+3}\times \mathbb{S}^1$, which we want to see as a toy model for Einstein equations with additional compact dimensions. We show global existence for small and regular initial data with polynomial decay at infinity. The method combines energy estimates on hyperboloids inside the light cone and weighted energy estimates outside the light cone.


title: Steering Majorana braiding via skyrmion-vortex pairs: a scalable platform

abstract: Majorana zero modes are quasiparticles that hold promise as building blocks for topological quantum computing. However, the litmus test for their detection, the observation of exotic non-abelian statistics revealed by braiding, has so far eluded experimental efforts. Here we take advantage of the fact that skyrmion-vortex pairs in superconductor-ferromagnet heterostructures harboring Majorana zero modes can be easily manipulated in two spatial dimensions. We adiabatically braid the hybrid topological structures and explicitly confirm the non-abelian statistics of the Majorana zero modes numerically using a self-consistent calculation of the superconducting order parameter. Our proposal of controlling skyrmion-vortex pairs provides the necessary leeway toward a scalable topological quantum computing platform.


title: A direct test of Auger cascade induced nucleation from heavy element contamination in C$_3$F$_8$ bubble chambers

abstract: Understanding and quantifying the gamma-induced bubble nucleation background in clean nuclear recoil detection bubble chambers is of utmost importance to bubble chamber based dark matter searches. We present data confirming the hypothesis that large Auger cascades from high-Z elements such as iodine and xenon dramatically increase the response of C$_3$F$_8$ bubble chambers to gamma rays. These tests, performed with a small calibration bubble chamber filled with C$_3$F$_8$+$\mathcal{O}$(10) ppm xenon, show that the probability of bubble nucleation scales with the rate of xenon inner-shell vacancies, reaching values >$10\%$ per K-Shell vacancy for Seitz thresholds of interest to future dark matter searches in bubble chambers. We also place an upper limit on bubble nucleation probability for argon Auger events, relevant to large future bubble chambers which may contain some residual atmospheric argon after the active fluid fill.


title: Combining Recurrent, Convolutional, and Continuous-time Models with Linear State-Space Layers

abstract: Recurrent neural networks (RNNs), temporal convolutions, and neural differential equations (NDEs) are popular families of deep learning models for time-series data, each with unique strengths and tradeoffs in modeling power and computational efficiency. We introduce a simple sequence model inspired by control systems that generalizes these approaches while addressing their shortcomings. The Linear State-Space Layer (LSSL) maps a sequence $u \mapsto y$ by simply simulating a linear continuous-time state-space representation $\dot{x} = Ax + Bu, y = Cx + Du$. Theoretically, we show that LSSL models are closely related to the three aforementioned families of models and inherit their strengths. For example, they generalize convolutions to continuous-time, explain common RNN heuristics, and share features of NDEs such as time-scale adaptation. We then incorporate and generalize recent theory on continuous-time memorization to introduce a trainable subset of structured matrices $A$ that endow LSSLs with long-range memory. Empirically, stacking LSSL layers into a simple deep neural network obtains state-of-the-art results across time series benchmarks for long dependencies in sequential image classification, real-world healthcare regression tasks, and speech. On a difficult speech classification task with length-16000 sequences, LSSL outperforms prior approaches by 24 accuracy points, and even outperforms baselines that use hand-crafted features on 100x shorter sequences.


title: Fair Sequential Selection Using Supervised Learning Models

abstract: We consider a selection problem where sequentially arrived applicants apply for a limited number of positions/jobs. At each time step, a decision maker accepts or rejects the given applicant using a pre-trained supervised learning model until all the vacant positions are filled. In this paper, we discuss whether the fairness notions (e.g., equal opportunity, statistical parity, etc.) that are commonly used in classification problems are suitable for the sequential selection problems. In particular, we show that even with a pre-trained model that satisfies the common fairness notions, the selection outcomes may still be biased against certain demographic groups. This observation implies that the fairness notions used in classification problems are not suitable for a selection problem where the applicants compete for a limited number of positions. We introduce a new fairness notion, ``Equal Selection (ES),'' suitable for sequential selection problems and propose a post-processing approach to satisfy the ES fairness notion. We also consider a setting where the applicants have privacy concerns, and the decision maker only has access to the noisy version of sensitive attributes. In this setting, we can show that the perfect ES fairness can still be attained under certain conditions.


title: Learning Collaborative Policies to Solve NP-hard Routing Problems

abstract: Recently, deep reinforcement learning (DRL) frameworks have shown potential for solving NP-hard routing problems such as the traveling salesman problem (TSP) without problem-specific expert knowledge. Although DRL can be used to solve complex problems, DRL frameworks still struggle to compete with state-of-the-art heuristics showing a substantial performance gap. This paper proposes a novel hierarchical problem-solving strategy, termed learning collaborative policies (LCP), which can effectively find the near-optimum solution using two iterative DRL policies: the seeder and reviser. The seeder generates as diversified candidate solutions as possible (seeds) while being dedicated to exploring over the full combinatorial action space (i.e., sequence of assignment action). To this end, we train the seeder's policy using a simple yet effective entropy regularization reward to encourage the seeder to find diverse solutions. On the other hand, the reviser modifies each candidate solution generated by the seeder; it partitions the full trajectory into sub-tours and simultaneously revises each sub-tour to minimize its traveling distance. Thus, the reviser is trained to improve the candidate solution's quality, focusing on the reduced solution space (which is beneficial for exploitation). Extensive experiments demonstrate that the proposed two-policies collaboration scheme improves over single-policy DRL framework on various NP-hard routing problems, including TSP, prize collecting TSP (PCTSP), and capacitated vehicle routing problem (CVRP).


title: First ${\beta}$-decay spectroscopy of $^{135}$In and new ${\beta}$-decay branches of $^{134}$In

abstract: The $\beta$ decay of the neutron-rich $^{134}$In and $^{135}$In was investigated experimentally in order to provide new insights into the nuclear structure of the tin isotopes with magic proton number $Z=50$ above the $N=82$ shell. The $\beta$-delayed $\gamma$-ray spectroscopy measurement was performed at the ISOLDE facility at CERN, where indium isotopes were selectively laser-ionized and on-line mass separated. Three $\beta$-decay branches of $^{134}$In were established, two of which were observed for the first time. Population of neutron-unbound states decaying via $\gamma$ rays was identified in the two daughter nuclei of $^{134}$In, $^{134}$Sn and $^{133}$Sn, at excitation energies exceeding the neutron separation energy by 1 MeV. The $\beta$-delayed one- and two-neutron emission branching ratios of $^{134}$In were determined and compared with theoretical calculations. The $\beta$-delayed one-neutron decay was observed to be dominant $\beta$-decay branch of $^{134}$In even though the Gamow-Teller resonance is located substantially above the two-neutron separation energy of $^{134}$Sn. Transitions following the $\beta$ decay of $^{135}$In are reported for the first time, including $\gamma$ rays tentatively attributed to $^{135}$Sn. In total, six new levels were identified in $^{134}$Sn on the basis of the $\beta \gamma \gamma$ coincidences observed in the $^{134}$In and $^{135}$In $\beta$ decays. A transition that might be a candidate for deexciting the missing neutron single-particle $13/2^+$ state in $^{133}$Sn was observed in both $\beta$ decays and its assignment is discussed. Experimental level schemes of $^{134}$Sn and $^{135}$Sn are compared with shell-model predictions. Using the fast timing technique, half-lives of the $2^+$, $4^+$ and $6^+$ levels in $^{134}$Sn were determined.


title: Revisiting Batch Norm Initialization

abstract: Batch normalization (BN) is comprised of a normalization component followed by an affine transformation and has become essential for training deep neural networks. Standard initialization of each BN in a network sets the affine transformation scale and shift to 1 and 0, respectively. However, after training we have observed that these parameters do not alter much from their initialization. Furthermore, we have noticed that the normalization process can still yield overly large values, which is undesirable for training. We revisit the BN formulation and present a new initialization method and update approach for BN to address the aforementioned issues. Experiments are designed to emphasize and demonstrate the positive influence of proper BN scale initialization on performance, and use rigorous statistical significance tests for evaluation. The approach can be used with existing implementations at no additional computational cost. Source code is available at https://github.com/osu-cvl/revisiting-bn-init.


title: Gravitational distortion on photon state at the vicinity of the Earth

abstract: As a photon propagates along a null geodesic, the space-time curvature around the geodesic distorts its wave function. We utilise the Fermi coordinates adapted to a general null geodesic, and derive the equation for interaction between the Riemann tensor and the photon wave function. The equation is solved by being mapped to a time-dependent Schr\"odinger equation in $(2+1)$ dimensions. The results show that as a Gaussian time-bin wavepacket with a narrow bandwidth travels over a null geodesic, it gains an extra phase that is a function of the Riemann tensor evaluated and integrated over the propagation trajectory. This extra phase is calculated for communication between satellites around the Earth, and is shown to be measurable by current technology.


title: Generic character of charge and spin density waves in superconducting cuprates

abstract: Understanding the nature of charge density waves (CDW) in cuprate superconductors has been complicated by material specific differences. A striking example is the opposite doping dependence of the CDW ordering wavevector in La-based and Y-based compounds, the two families where charge ordering is strongest and best characterized. Here we report a combined resonant soft X-ray scattering (RSXS) and neutron scattering study of charge and spin density waves in isotopically enriched La$_{1.8-x}$ Eu$_{0.2}$ Sr$_{x}$ CuO$_{4}$ over a range of doping $0.07 \leq x \leq 0.20$. For all dopings studied by RSXS, we find that the CDW amplitude is approximately temperature-independent and develops well above experimentally accessible temperatures. Surprisingly, the CDW ordering wavevector shows a non-monotonic temperature dependence, with a sudden change occurring at temperatures near the SDW onset temperature. We describe this behavior with a Landau-Ginzburg theory for an incommensurate CDW in a metallic system with a finite charge compressibility and CDW-SDW coupling. Our Landau-Ginzburg analysis suggests that the ordering wavevector at high temperatures decreases with increased doping. This behavior is opposite to the trend at low temperatures and highly reminiscent of the doping dependence seen in YBa$_2$ Cu$_3$ O$_{6+\delta}$ , suggesting a common origin of the CDW in hole-doped cuprate superconductors.