title: Provably accurate simulation of gauge theories and bosonic systems

abstract: Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit $\Lambda$ on each local quantum number, and if the initial state has low local quantum numbers, then an error at most $\epsilon$ can be achieved by choosing $\Lambda$ to scale polylogarithmically with $\epsilon^{-1}$, an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on $\Lambda$ that achieves accuracy $\epsilon$, obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy $\epsilon$ by truncating local quantum numbers at $\Lambda=\textrm{polylog}(\epsilon^{-1})$.


title: $6$D Anomaly-Free Matter Spectrum in F-theory on Singular Spaces

abstract: In this paper we study the 6d localized charged matter spectrum of F-theory directly on a singular elliptic Calabi-Yau 3-fold, i.e. without smoothing via resolution or deformation of the entire fibration. Given only the base surface, discriminant locus, and the $SL(2,\mathbb{Z})$ local system, we propose a general prescription for determining the charged matter spectrum localized at intersections of seven-branes, using the technology of string junctions. More precisely, at each codimension-$2$ collision of seven-branes, we determine the local seven-brane content and compute the number of massless string junctions modulo the action of the $SL(2,\mathbb{Z})$ monodromy. We find agreement with the predicted results from $6$d anomaly cancellation in all cases considered. Examples include a generic Weierstrass model with arbitrary Kodaira fiber intersecting an $I_1$, as well as cases with jointly charged matter localized at intersections of non-abelian seven-branes.


title: Is the high-energy neutrino event IceCube-200530A associated with a hydrogen-rich superluminous supernova?

abstract: The Zwicky Transient Facility (ZTF) follow-up campaign of alerts released by the IceCube Neutrino Observatory has led to the likely identification of the transient AT2019fdr as the source of the neutrino event IC200530A. AT2019fdr was initially suggested to be a tidal disruption event in a Narrow-Line Seyfert 1 galaxy. However, the combination of its spectral properties, color evolution, and feature-rich light curve suggests that AT2019fdr may be a Type IIn superluminous supernova. In the latter scenario, IC200530A may have been produced via inelastic proton-proton collisions between the relativistic protons accelerated at the forward shock and the cold protons of the circumstellar medium. Here, we investigate this possibility and find that at most $4.6\times 10^{-2}$ muon neutrino and antineutrino events are expected to be detected by the IceCube Neutrino Observatory within $394$ days of discovery in the case of excellent discrimination of the atmospheric background. After correcting for the Eddington bias, which occurs when a single cosmic neutrino event is adopted to infer the neutrino emission at the source, we conclude that IC200530A may originate from the hydrogen-rich superluminous supernova AT2019fdr.


title: Lyman Continuum Galaxy Candidates in COSMOS

abstract: Star-forming galaxies are the sources likely to have reionized the universe. As we cannot observe them directly due to the opacity of the intergalactic medium at $z\gtrsim5$, we study $z\sim3\text{--}5$ galaxies as proxies to place observational constraints on cosmic reionization. Using new deep \textit{Hubble Space Telescope} rest-frame UV F336W and F435W imaging (30-orbit, $\sim40$~arcmin$^2$, $\sim29\text{--}30$~mag depth at 5$\sigma$), we attempt to identify a sample of Lyman continuum galaxies (LCGs). These are individual sources that emit ionizing flux below the Lyman break ($<912~\text{\AA}$). This population would allow us to constrain cosmic reionization parameters such as the number density and escape fraction ($f_{\rm esc}$) of ionizing sources. We compile a comprehensive parent sample that does not rely on the Lyman-break technique for redshifts. We present three new spectroscopic candidates at $z\sim3.7\text{--}4.4$, and 32 new photometric candidates. The high-resolution multi-band HST imaging and new Keck/Low Resolution Imaging Spectrometer (LRIS) redshifts make these promising spectroscopic LCG candidates. Using both a traditional and probabilistic approach, we find the most likely $f_{\rm esc}$ values for the three spectroscopic LCG candidates are $>100\%$, and therefore not physical. We are unable to confirm the true nature of these sources with the best available imaging and direct blue Keck/LRIS spectroscopy. More spectra, especially from the new class of 30 m telescopes, will be required to build a statistical sample of LCGs to place firm observational constraints on cosmic reionization.


title: Scaling relations and baryonic cycling in local star-forming galaxies. III. Outflows, effective yields and metal loading factors

abstract: Gas accretion and stellar feedback processes link metal content, star formation, and gas and stellar mass (and the potential depth) in star-forming galaxies. Constraining this hypersurface has been challenging because of the need for measurements of HI and HII gas masses spanning a broad parameter space. A recent step forward has been achieved through the "Metallicity And Gas for Mass Assembly" (MAGMA) sample of local star-forming galaxies, which consists of homogeneously-determined parameters and a significant quantity of dwarf galaxies, with stellar masses as low as $\sim 10^5 - 10^{6}\, M_{\odot}$. Here, we adopt a "standard" galactic chemical evolution model, with which we can quantify stellar-driven outflows. In particular, we constrain the difference between the mass-loading in accretion and outflows and the wind metal-loading factor. The resulting model reproduces very well the local mass-metallicity relation, and the observed trends of metallicity with gas fraction. Although the difference in mass loading between accreted and expelled gas is extremely difficult to constrain, we find indications that, on average, the amount of gas acquired through accretion is roughly the same as the gas lost through bulk stellar outflows, corresponding to a "gas equilibrium" scenario. In agreement with previous work, the wind metal-loading factor shows a steep increase toward lower mass and circular velocity, indicating that low-mass galaxies are more efficient at expelling metals, thus shaping the mass-metallicity relation. Effective yields are found to increase with mass up to an inflection mass threshold, with a mild decline at larger masses and circular velocities. A comparison of our results for metal loading in outflows with the expectations for their mass loading favors momentum-driven winds at low masses, rather than energy-driven ones. (abridged)


title: KiDS & Euclid: Cosmological implications of a pseudo angular power spectrum analysis of KiDS-1000 cosmic shear tomography

abstract: We present a tomographic weak lensing analysis of the Kilo Degree Survey Data Release 4 (KiDS-1000), using a new pseudo angular power spectrum estimator (pseudo-$C_{\ell}$) under development for the ESA Euclid mission. Over 21 million galaxies with shape information are divided into five tomographic redshift bins, ranging from 0.1 to 1.2 in photometric redshift. We measured pseudo-$C_{\ell}$ using eight bands in the multipole range $76<\ell<1500$ for auto- and cross-power spectra between the tomographic bins. A series of tests were carried out to check for systematic contamination from a variety of observational sources including stellar number density, variations in survey depth, and point spread function properties. While some marginal correlations with these systematic tracers were observed, there is no evidence of bias in the cosmological inference. B-mode power spectra are consistent with zero signal, with no significant residual contamination from E/B-mode leakage. We performed a Bayesian analysis of the pseudo-$C_{\ell}$ estimates by forward modelling the effects of the mask. Assuming a spatially flat $\Lambda$CDM cosmology, we constrained the structure growth parameter $S_8 = \sigma_8(\Omega_{\rm m}/0.3)^{1/2} = 0.754_{-0.029}^{+0.027}$. When combining cosmic shear from KiDS-1000 with baryon acoustic oscillation and redshift space distortion data from recent Sloan Digital Sky Survey (SDSS) measurements of luminous red galaxies, as well as the Lyman-$\alpha$ forest and its cross-correlation with quasars, we tightened these constraints to $S_8 = 0.771^{+0.006}_{-0.032}$. These results are in very good agreement with previous KiDS-1000 and SDSS analyses and confirm a $\sim 3\sigma$ tension with early-Universe constraints from cosmic microwave background experiments.


title: Challenges for Unsupervised Anomaly Detection in Particle Physics

abstract: Anomaly detection relies on designing a score to determine whether a particular event is uncharacteristic of a given background distribution. One way to define a score is to use autoencoders, which rely on the ability to reconstruct certain types of data (background) but not others (signals). In this paper, we study some challenges associated with variational autoencoders, such as the dependence on hyperparameters and the metric used, in the context of anomalous signal (top and $W$) jets in a QCD background. We find that the hyperparameter choices strongly affect the network performance and that the optimal parameters for one signal are non-optimal for another. In exploring the networks, we uncover a connection between the latent space of a variational autoencoder trained using mean-squared-error and the optimal transport distances within the dataset. We then show that optimal transport distances to representative events in the background dataset can be used directly for anomaly detection, with performance comparable to the autoencoders. Whether using autoencoders or optimal transport distances for anomaly detection, we find that the choices that best represent the background are not necessarily best for signal identification. These challenges with unsupervised anomaly detection bolster the case for additional exploration of semi-supervised or alternative approaches.


title: Euclidean D-branes in Type IIB string theory on Calabi-Yau threefolds

abstract: We compute the contribution of Euclidean D-branes in type IIB string theory on Calabi-Yau threefolds to the metric on the hypermultiplet moduli space in the large volume, weak coupling limit. Our results are in perfect agreement with the predictions based on S-duality, mirror symmetry and supersymmetry.


title: $S^5$: The Orbital and Chemical Properties of One Dozen Stellar Streams

abstract: We report the kinematic, orbital, and chemical properties of 12 stellar streams with no evident progenitors, using line-of-sight velocities and metallicities from the Southern Stellar Stream Spectroscopic Survey ($S^5$), proper motions from $Gaia$ EDR3, and distances derived from distance tracers or the literature. This data set provides the largest homogeneously analyzed set of streams with full 6D kinematics and metallicities. All streams have heliocentric distances between ${\sim}10-50$ kpc. The velocity and metallicity dispersions show that half of the stream progenitors were disrupted dwarf galaxies (DGs), while the other half originated from disrupted globular clusters (GCs), hereafter referred to as DG and GC streams. Based on the mean metallicities of the streams and the mass-metallicity relation, the luminosities of the progenitors of the DG streams range between Carina and Ursa Major I ($-9.5\lesssim M_V\lesssim-5.5$). Four of the six GC streams have mean metallicities of [Fe/H]$< -2$, more metal-poor than typical Milky Way (MW) GCs at similar distances. Interestingly, the 300S and Jet GC streams are the only streams on retrograde orbits in our dozen stream sample. Finally, we compare the orbital properties of the streams with known DGs and GCs in the MW, finding several possible associations. Some streams appear to have been accreted with the recently discovered Gaia-Enceladus-Sausage system, and others suggest that GCs were formed in and accreted together with the progenitors of DG streams whose stellar masses are similar to Draco to Carina ($\sim10^5-10^6M_\odot$).


title: Absence of superdiffusion in certain random spin models

abstract: The dynamics of spin at finite temperature in the spin-1/2 Heisenberg chain was found to be superdiffusive in numerous recent numerical and experimental studies. Theoretical approaches to this problem have emphasized the role of nonabelian SU(2) symmetry as well as integrability but the associated methods cannot be readily applied when integrability is broken. We examine spin transport in a spin-1/2 chain in which the exchange couplings fluctuate in space and time around a nonzero mean $J$, a model introduced by De Nardis et al. [Phys. Rev. Lett. 127, 057201 (2021)]. We show that operator dynamics in the strong noise limit at infinite temperature can be analyzed using conventional perturbation theory as an expansion in $J$. We find that regular diffusion persists at long times, albeit with an enhanced diffusion constant. The finite time spin dynamics is analyzed and compared with matrix product operator simulations.


title: The impact of metallicity-dependent dust destruction on the dust-to-metals ratio in galaxies

abstract: The ratio of the mass of interstellar dust to the total mass of metals (the dust-to-metals/DTM ratio) tends to increase with metallicity. This can be explained by the increasing efficiency of grain growth in the interstellar medium (ISM) at higher metallicities, with a corollary being that the low DTM ratios seen at low metallicities are due to inefficient stellar dust production. This interpretation assumes that the efficiency of dust destruction in the ISM is constant, whereas it might be expected to increase at low metallicity; the decreased cooling efficiency of low-metallicity gas should result in more post-shock dust destruction via thermal sputtering. We show that incorporating a sufficiently strong metallicity dependence into models of galaxy evolution removes the need for low stellar dust yields. The contribution of stellar sources to the overall dust budget may be significantly underestimated, and that of grain growth overestimated, by models assuming a constant destruction efficiency.


title: Top-quark mass effects in H+jet and H+2 jets production

abstract: We present calculations of Higgs boson production via gluon-gluon fusion in association with one or two additional jets at next-to-leading order in QCD. The calculation of H+jet is exact in the treatment of the top-quark mass, whereas for the H+2 jets calculation the two-loop virtual amplitudes are approximated via a reweighting with leading-order mass effects, while keeping all top-quark mass effects in the real radiation contributions. For H+jet production, this study extends a previous calculation, revealing an error in the previous results. For total and differential cross sections, we present new results and compare the QCD corrections with the infinite top-mass limit, for which we find a strikingly good agreement if all amplitudes are rescaled by the leading-order mass dependence.


title: Characterizing quantum instruments: from non-demolition measurements to quantum error correction

abstract: In quantum information processing quantum operations are often processed alongside measurements which result in classical data. Due to the information gain of classical measurement outputs non-unitary dynamical processes can take place on the system, for which common quantum channel descriptions fail to describe the time evolution. Quantum measurements are correctly treated by means of so-called quantum instruments capturing both classical outputs and post-measurement quantum states. Here we present a general recipe to characterize quantum instruments alongside its experimental implementation and analysis. Thereby, the full dynamics of a quantum instrument can be captured, exhibiting details of the quantum dynamics that would be overlooked with common tomography techniques. For illustration, we apply our characterization technique to a quantum instrument used for the detection of qubit loss and leakage, which was recently implemented as a building block in a quantum error correction (QEC) experiment (Nature 585, 207-210 (2020)). Our analysis reveals unexpected and in-depth information about the failure modes of the implementation of the quantum instrument. We then numerically study the implications of these experimental failure modes on QEC performance, when the instrument is employed as a building block in QEC protocols on a logical qubit. Our results highlight the importance of careful characterization and modelling of failure modes in quantum instruments, as compared to simplistic hardware-agnostic phenomenological noise models, which fail to predict the undesired behavior of faulty quantum instruments. The presented methods and results are directly applicable to generic quantum instruments.


title: Deprojection of external barred galaxies from photometry

abstract: The observations of external galaxies are projected to the 2D sky plane. Reconstructing the 3D intrinsic density distribution of a galaxy from the 2D image is challenging, especially for barred galaxies, but is a critical step for constructing galactic dynamical models. Here we present a method for deprojecting barred galaxies and we validate the method by testing against mock images created from an N-body simulation with a peanut-shaped bar. We decompose a galaxy image into a bulge (including a bar) and a disk. By subtracting the disk from the original image a barred bulge remains. We perform multi-Gaussian expansion (MGE) fit to each component, then we deproject them separately by considering the barred bulge is triaxial while the disk is axisymmetric. We restrict the barred bulge to be aligned in the disk plane and has a similar thickness to the disk in the outer regions. The 3D density distribution is thus constructed by combining the barred bulge and the disk. Our model can generally recover the 3D density distribution of disk and inner barred bulge regions, although not a perfect match to the peanut-shaped structure. By using the same initial conditions, we integrate the orbits in our model-inferred potential and the true potential by freezing the N-body simulation. We find that 85% of all these orbits have similar morphologies in these two potentials, and our model supports the orbits that generate a boxy/peanut-shaped structure and an elongated bar similar to these in the true potential.


title: Considering user agreement in learning to predict the aesthetic quality

abstract: How to robustly rank the aesthetic quality of given images has been a long-standing ill-posed topic. Such challenge stems mainly from the diverse subjective opinions of different observers about the varied types of content. There is a growing interest in estimating the user agreement by considering the standard deviation of the scores, instead of only predicting the mean aesthetic opinion score. Nevertheless, when comparing a pair of contents, few studies consider how confident are we regarding the difference in the aesthetic scores. In this paper, we thus propose (1) a re-adapted multi-task attention network to predict both the mean opinion score and the standard deviation in an end-to-end manner; (2) a brand-new confidence interval ranking loss that encourages the model to focus on image-pairs that are less certain about the difference of their aesthetic scores. With such loss, the model is encouraged to learn the uncertainty of the content that is relevant to the diversity of observers' opinions, i.e., user disagreement. Extensive experiments have demonstrated that the proposed multi-task aesthetic model achieves state-of-the-art performance on two different types of aesthetic datasets, i.e., AVA and TMGA.


title: Constraints on the maximum densities of neutron stars from postmerger gravitational waves with third-generation observations

abstract: Using data from 289 numerical relativity simulations of merging binary neutron stars, we identify, for the first time, a robust quasi-universal relation connecting the postmerger peak gravitational-wave frequency and the value of the density at the center of the maximum mass nonrotating neutron star. This relation offers a new possibility for precision equation-of-state constraints with next-generation ground-based gravitational-wave interferometers. Mock Einstein Telescope observations of fiducial events indicate that Bayesian inferences can constrain the maximum density to ${\sim}15\%$ ($90\%$ confidence level) for a single signal at the minimum sensitivity threshold for a detection. If the postmerger signal is included in a full-spectrum (inspiral-merger-postmerger) analysis of such signal, the pressure-density function can be tightly constrained up to the maximum density, and the maximum neutron star mass can be measured with an accuracy better than $12\%$ ($90\%$ confidence level).


title: The impact of cosmic-ray attenuation on the carbon cycle emission in molecular clouds

abstract: Observations of the emission of the carbon cycle species (C, C+ CO) are commonly used to diagnose gas properties in the interstellar medium but are significantly sensitive to the cosmic-ray ionization rate. The carbon-cycle chemistry is known to be quite sensitive to the cosmic-ray ionization rate, $\zeta$, controlled by the flux of low-energy cosmic rays which get attenuated through molecular clouds. However, astrochemical models commonly assume a constant cosmic-ray ionization rate in the clouds. We investigate the effect of cosmic-ray attenuation on the emission of carbon cycle species from molecular clouds, in particular the [CII] 158 $\mu$m, [CI] 609 $\mu$m and CO (J = 1 - 0) 115.27 GHz lines. We use a post-processed chemical model of diffuse and dense simulated molecular clouds and quantify the variation in both column densities and velocity integrated line emission of the carbon cycle with different cosmic-ray ionization rate models. We find that the abundances and column densities of carbon cycle species is significantly impacted by the chosen cosmic-ray ionization rate model: no single constant ionization rate can reproduce the abundances modelled with an attenuated cosmic-ray model. Further, we show that constant ionization rate models fail to simultaneously reproduce the integrated emission of the lines we consider, and their deviations from a physically derived cosmic-ray attenuation model is too complex to be simply corrected. We demonstrate that the two clouds we model exhibit a similar average $A_{\rm V, eff}$ -- $n_{\rm H}$ relationship, resulting in an average relation between the cosmic-ray ionization rate and density $\zeta(n_{\rm H})$. We conclude by providing a number of implementation recommendations for CRs in astrochemical models, but emphasize the necessity for column-dependent cosmic-ray ionization rate prescriptions.


title: Supercurrent-induced topological phase transitions

abstract: We show that finite current in superconductors can induce topological phase transitions, as a result of the deformation of the quasiparticle spectrum by a finite center-of-mass (COM) momentum of the Cooper pairs. To show the wide applicability of this mechanism, we examine the topological properties of three prototypical systems, the Kitaev chain, $s$-wave superconductors, and $d$-wave superconductors. We introduce a finite COM momentum as an external field corresponding to supercurrent and show that all the models exhibit current-induced topological phase transitions. We also discuss the possibility of observing the phase transitions in experiments and the relation to the other finite COM momentum pairing states.

categories: cond-mat.supr-con cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el


title: Spin-lasing in bimodal quantum dot micropillar cavities

abstract: Spin-controlled lasers are highly interesting photonic devices and have been shown to provide ultra-fast polarization dynamics in excess of 200 GHz. In contrast to conventional semiconductor lasers their temporal properties are not limited by the intensity dynamics, but are governed primarily by the interaction of the spin dynamics with the birefringent mode splitting that determines the polarization oscillation frequency. Another class of modern semiconductor lasers are high-beta emitters which benefit from enhanced light-matter interaction due to strong mode confinement in low-mode-volume microcavities. In such structures, the emission properties can be tailored by the resonator geometry to realize for instance bimodal emission behavior in slightly elliptical micropillar cavities. We utilize this attractive feature to demonstrate and explore spin-lasing effects in bimodal high-beta quantum dot micropillar lasers. The studied microlasers with a beta-factor of 4% show spin-laser effects with experimental polarization oscillation frequencies up to 15 GHz and predicted frequencies up to about 100 GHz which are controlled by the ellipticity of the resonator. Our results reveal appealing prospects for very compact, ultra-fast and energy-efficient spin-lasers and can pave the way for future purely electrically injected spin-lasers enabled by short injection path lengths.


title: Language Modelling via Learning to Rank

abstract: We consider language modelling (LM) as a multi-label structured prediction task by re-framing training from solely predicting a single ground-truth word to ranking a set of words which could continue a given context. To avoid annotating top-$k$ ranks, we generate them using pre-trained LMs: GPT-2, BERT, and Born-Again models. This leads to a rank-based form of knowledge distillation (KD). We also develop a method using $N$-grams to create a non-probabilistic teacher which generates the ranks without the need of a pre-trained LM. We confirm the hypotheses that we can treat LMing as a ranking task and that we can do so without the use of a pre-trained LM. We show that rank-based KD generally improves perplexity (PPL), often with statistical significance, when compared to Kullback-Leibler-based KD. Surprisingly, given the simplicity of the method, $N$-grams act as competitive teachers and achieve similar performance as using either BERT or a Born-Again model teachers. GPT-2 always acts as the best teacher, though, and using it and a Transformer-XL student on Wiki-02, rank-based KD reduces a cross-entropy baseline from 65.27 to 55.94 and against a KL-based KD of 56.70.


title: Open-Domain Question-Answering for COVID-19 and Other Emergent Domains

abstract: Since late 2019, COVID-19 has quickly emerged as the newest biomedical domain, resulting in a surge of new information. As with other emergent domains, the discussion surrounding the topic has been rapidly changing, leading to the spread of misinformation. This has created the need for a public space for users to ask questions and receive credible, scientific answers. To fulfill this need, we turn to the task of open-domain question-answering, which we can use to efficiently find answers to free-text questions from a large set of documents. In this work, we present such a system for the emergent domain of COVID-19. Despite the small data size available, we are able to successfully train the system to retrieve answers from a large-scale corpus of published COVID-19 scientific papers. Furthermore, we incorporate effective re-ranking and question-answering techniques, such as document diversity and multiple answer spans. Our open-domain question-answering system can further act as a model for the quick development of similar systems that can be adapted and modified for other developing emergent domains.


title: Finite-time teleportation phase transition in random quantum circuits

abstract: How long does it take to entangle two distant qubits in a quantum circuit evolved by generic unitary dynamics? We show that if the time evolution is followed by measurements of all but two infinitely separated test qubits, then the entanglement between them can undergo a phase transition and become nonzero at a finite critical time $t_c$. The fidelity of teleporting a quantum state from an input qubit to an infinitely distant output qubit shows the same critical onset. Specifically, these finite-time transitions occur in short-range interacting two-dimensional random unitary circuits and in sufficiently long-range interacting one-dimensional circuits. The phase transition is understood by mapping the random continuous-time evolution to a finite-temperature thermal state of an effective spin Hamiltonian, where the inverse temperature equals the evolution time in the circuit. In this framework, the entanglement between two distant qubits at times $t>t_c$ corresponds to the emergence of long-range ferromagnetic spin correlations below the critical temperature. We verify these predictions using numerical simulation of Clifford circuits and propose potential realizations in existing platforms for quantum simulation.


title: The Complexity of Bipartite Gaussian Boson Sampling

abstract: Gaussian boson sampling is a model of photonic quantum computing that has attracted attention as a platform for building quantum devices capable of performing tasks that are out of reach for classical devices. There is therefore significant interest, from the perspective of computational complexity theory, in solidifying the mathematical foundation for the hardness of simulating these devices. We show that, under the standard Anti-Concentration and Permanent-of-Gaussians conjectures, there is no efficient classical algorithm to sample from ideal Gaussian boson sampling distributions (even approximately) unless the polynomial hierarchy collapses. The hardness proof holds in the regime where the number of modes scales quadratically with the number of photons, a setting in which hardness was widely believed to hold but that nevertheless had no definitive proof. Crucial to the proof is a new method for programming a Gaussian boson sampling device so that the output probabilities are proportional to the permanents of submatrices of an arbitrary matrix. This technique is a generalization of Scattershot BosonSampling that we call BipartiteGBS. We also make progress towards the goal of proving hardness in the regime where there are fewer than quadratically more modes than photons (i.e., the high-collision regime) by showing that the ability to approximate permanents of matrices with repeated rows/columns confers the ability to approximate permanents of matrices with no repetitions. The reduction suffices to prove that GBS is hard in the constant-collision regime.


title: Tutorial on laser linewidths

abstract: In this tutorial, the physical origins and mathematical analyses of laser linewidths are reviewed. The semi-classical model is based on an equation for the light-mode amplitude that includes random source terms, one term for each process that affects the amplitude (stimulated and spontaneous emission, stimulated absorption, and facet and material loss). Although the source terms are classical, their assigned strengths are consistent with the laws of quantum optics. Analysis of this equation shows that the laser linewidth is proportional to the sum of the (positive) source strengths for all gain and loss processes. Three-level and semiconductor lasers have broader linewidths than comparable four-level lasers, because stimulated absorption and the stimulated emission that compensates it both contribute to the linewidth.


title: Tutorial on stochastic systems

abstract: In this tutorial, three examples of stochastic systems are considered: A strongly-damped oscillator, a weakly-damped oscillator and an undamped oscillator (integrator) driven by noise. The evolution of these systems is characterized by the temporal correlation functions and spectral densities of their displacements, which are determined and discussed. Damped oscillators reach steady stochastic states. Their correlations are decreasing functions of the difference between the sample times and their spectra have peaks near their resonance frequencies. An undamped oscillator never reaches a steady state. Its energy increases with time and its spectrum is sharply peaked at low frequencies. The required mathematical methods and physical concepts are explained on a just-in-time basis, and some theoretical pitfalls are mentioned. The insights one gains from studies of oscillators can be applied to a wide variety of physical systems, such as atom and semiconductor lasers, which will be discussed in a subsequent tutorial.


title: Radiative Generation of dS from AdS

abstract: The large volume scenario (LVS) of type IIB string compactifications has a robust supersymmetry breaking minimum with a negative cosmological constant (CC). We argue that radiative corrections below the Kaluza-Klein (KK) scale can result in a positive CC, though the string theory generated CC is negative, if some mild conditions on the spectrum of low energy fluctuations are satisfied. This would make the so-called deSitter swampland conjecture (even if true at a high scale) physically irrelevant.


title: Interpretable AI forecasting for numerical relativity waveforms of quasi-circular, spinning, non-precessing binary black hole mergers

abstract: We present a deep-learning artificial intelligence model that is capable of learning and forecasting the late-inspiral, merger and ringdown of numerical relativity waveforms that describe quasi-circular, spinning, non-precessing binary black hole mergers. We used the NRHybSur3dq8 surrogate model to produce train, validation and test sets of $\ell=|m|=2$ waveforms that cover the parameter space of binary black hole mergers with mass-ratios $q\leq8$ and individual spins $|s^z_{\{1,2\}}| \leq 0.8$. These waveforms cover the time range $t\in[-5000\textrm{M}, 130\textrm{M}]$, where $t=0M$ marks the merger event, defined as the maximum value of the waveform amplitude. We harnessed the ThetaGPU supercomputer at the Argonne Leadership Computing Facility to train our AI model using a training set of 1.5 million waveforms. We used 16 NVIDIA DGX A100 nodes, each consisting of 8 NVIDIA A100 Tensor Core GPUs and 2 AMD Rome CPUs, to fully train our model within 3.5 hours. Our findings show that artificial intelligence can accurately forecast the dynamical evolution of numerical relativity waveforms in the time range $t\in[-100\textrm{M}, 130\textrm{M}]$. Sampling a test set of 190,000 waveforms, we find that the average overlap between target and predicted waveforms is $\gtrsim99\%$ over the entire parameter space under consideration. We also combined scientific visualization and accelerated computing to identify what components of our model take in knowledge from the early and late-time waveform evolution to accurately forecast the latter part of numerical relativity waveforms. This work aims to accelerate the creation of scalable, computationally efficient and interpretable artificial intelligence models for gravitational wave astrophysics.


title: Full-shape cosmology analysis of SDSS-III BOSS galaxy power spectrum using emulator-based halo model: a $5\%$ determination of $\sigma_8$

abstract: We present the results obtained from the full-shape cosmology analysis of the redshift-space power spectra for 4 galaxy samples of the SDSS-III BOSS DR12 galaxy catalog over $0.2 < z < 0.75$. For the theoretical template, we use an emulator that was built from an ensemble set of $N$-body simulations, which enables fast and accurate computation of the redshift-space power spectrum of halos. Combining with the halo occupation distribution to model the halo-galaxy connection, we can compute the redshift-space power spectrum of BOSS-like galaxies in the flat $\Lambda$CDM cosmology. In our cosmology inference, we use the power spectrum monopole, quadrupole and hexadecapole and include 7 nuisance parameters to model uncertainties in the halo-galaxy connection for each galaxy sample, but do not use any information on the abundance of galaxies. We demonstrate a validation of our analysis pipeline using the mock catalogs of BOSS-like galaxies, generated using different recipes of the halo-galaxy connection and including the assembly bias effect. Assuming weak priors on cosmological parameters, except for $\Omega_{\rm b}h^2$ and $n_{\rm s}$, we show that our model well reproduces the BOSS power spectra. Including the power spectrum information up to $k_{\rm max}=0.25\,h{\rm Mpc}^{-1}$, we find $\Omega_\mathrm{m}=0.301^{+0.012}_{-0.011}$, $H_0=68.2 \pm 1.4~\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$, and $\sigma_8=0.786^{+0.036}_{-0.037}$, for the mode and 68\% credible interval, after marginalization over nuisance parameters. We find little improvement in the cosmological parameters beyond a maximum wavelength $k_{\rm max}\simeq 0.2\,h\,{\rm Mpc}^{-1}$ due to the shot noise domination and marginalization of the halo-galaxy connection parameters. Our results are consistent with the Planck CMB results within $1\sigma$ statistical uncertainties.


title: Collective phenomena in pp interactions with high multiplicity

abstract: The latest results on the search for collective phenomena in pp interactions at the U-70 accelerator at IHEP (Protvino) are presented. This the long-term experiment has been carried out at the SVD-2 setup. We received the evidence of the pion (Bose-Einstein) condensate formation. Two noticeable peaks in the angular distribution of charged pions are interpreted as Cherenkov radiation of gluons by quarks. The leading effect disappears and the secondary system becomes isotropic in all directions at high multiplicity region.


title: AstraeusIV : Quantifying the star formation histories of galaxies in the Epoch of Reionization

abstract: We use the \textsc{astraeus} framework, that couples an N-body simulation with a semi-analytic model for galaxy formation and a semi-numerical model for reionization, to quantify the star formation histories (SFHs) of galaxies in the first billion years. Exploring four models of radiative feedback, we fit the SFH of each galaxy at $z>5$ as $\mathrm{log}(\mathrm{SFR}(z))=-\alpha(1 + z)+\beta$; star formation is deemed stochastic if it deviates from this fit by more than $\Delta_\mathrm{SFR}=0.6\,$dex. Our key findings are: (i) The fraction of stellar mass formed and time spent in the stochastic phase decrease with increasing stellar mass and redshift $z$. While galaxies with stellar masses of $M_\star\sim10^7M_\odot$ at $z\sim5~(10)$ form $\sim70\%~(20\%)$ of their stellar mass in the stochastic phase, this reduces to $<10\%$ at all redshifts for galaxies with $M_\star > 10^{10}M_\odot$; (ii) the fractional mass assembled and lifetime spent in the stochastic phase do not significantly change with the radiative feedback model used; (iii) at all redshifts, $\alpha$ increases (decreases for the strongest radiative feedback model) with stellar mass for galaxies with $M_\star\lesssim 10^{8.5}M_\odot$ and converges to $\sim0.18$ for more massive galaxies; $\beta$ always increases with stellar mass. Our proposed fits can reliably recover the stellar masses and mass-to-light ratios for galaxies with $M_\star\sim10^{8-10.5}M_\odot$ and $M_{UV}\sim-17~{\rm to}~-23$ at $z\sim 5-9$. This physical model can therefore be used to derive the SFHs for galaxies observed by a number of forthcoming instruments.


title: Block Contextual MDPs for Continual Learning

abstract: In reinforcement learning (RL), when defining a Markov Decision Process (MDP), the environment dynamics is implicitly assumed to be stationary. This assumption of stationarity, while simplifying, can be unrealistic in many scenarios. In the continual reinforcement learning scenario, the sequence of tasks is another source of nonstationarity. In this work, we propose to examine this continual reinforcement learning setting through the block contextual MDP (BC-MDP) framework, which enables us to relax the assumption of stationarity. This framework challenges RL algorithms to handle both nonstationarity and rich observation settings and, by additionally leveraging smoothness properties, enables us to study generalization bounds for this setting. Finally, we take inspiration from adaptive control to propose a novel algorithm that addresses the challenges introduced by this more realistic BC-MDP setting, allows for zero-shot adaptation at evaluation time, and achieves strong performance on several nonstationary environments.


title: Bounds on entries in Bianchi group generators

abstract: Upper and lower bounds are given for the maximum Euclidean curvature among faces in Bianchi's fundamental polyhedron for $PSL_2(O)$ in the upper-half space model of hyperbolic space, where $O$ is an imaginary quadratic ring of integers with discriminant $\Delta$. We prove these bounds are asymptotically within $(\log |\Delta|)^{8.54}$ of one another. This improves on the previous best upper-bound, which is roughly off by a factor between $\Delta^2$ and $|\Delta|^{5/2}$ depending on the smallest prime dividing $\Delta$. The gap between our upper and lower bounds is determined by an analog of Jacobsthal's function, introduced here for imaginary quadratic fields.


title: "Mysteries" of Modern Physics and the Fundamental Constants $c$, $h$, and $G$

abstract: We review how the kinematic structures of special relativity and quantum mechanics both stem from the relativity principle, i.e., "no preferred reference frame" (NPRF). Essentially, NPRF applied to the measurement of the speed of light $c$ gives the light postulate and leads to the geometry of Minkowski spacetime, while NPRF applied to the measurement of Planck's constant $h$ gives "average-only" projection and leads to the denumerable-dimensional Hilbert space of quantum mechanics. These kinematic structures contain the counterintuitive aspects ("mysteries") of time dilation, length contraction, and quantum entanglement. In this essay, we extend the application of NPRF to the gravitational constant $G$ and show that it leads to the "mystery" of the contextuality of mass in general relativity. Thus, we see an underlying coherence and integrity in modern physics via its "mysteries" and the fundamental constants $c$, $h$, and $G$.


title: Guided Policy Search using Sequential Convex Programming for Initialization of Trajectory Optimization Algorithms

abstract: Nonlinear trajectory optimization algorithms have been developed to handle optimal control problems with nonlinear dynamics and nonconvex constraints in trajectory planning. The performance and computational efficiency of many trajectory optimization methods are sensitive to the initial guess, i.e., the trajectory guess needed by the recursive trajectory optimization algorithm. Motivated by this observation, we tackle the initialization problem for trajectory optimization via policy optimization. To optimize a policy, we propose a guided policy search method that has two key components: i) Trajectory update; ii) Policy update. The trajectory update involves offline solutions of a large number of trajectory optimization problems from different initial states via Sequential Convex Programming (SCP). Here we take a single SCP step to generate the trajectory iterate for each problem. In conjunction with these iterates, we also generate additional trajectories around each iterate via a feedback control law. Then all these trajectories are used by a stochastic gradient descent algorithm to update the neural network policy, i.e., the policy update step. As a result, the trained policy makes it possible to generate trajectory candidates that are close to the optimality and feasibility and that provide excellent initial guesses for the trajectory optimization methods. We validate the proposed method via a real-world 6-degree-of-freedom powered descent guidance problem for a reusable rocket.


title: Representational Continuity for Unsupervised Continual Learning

abstract: Continual learning (CL) aims to learn a sequence of tasks without forgetting the previously acquired knowledge. However, recent CL advances are restricted to supervised continual learning (SCL) scenarios. Consequently, they are not scalable to real-world applications where the data distribution is often biased and unannotated. In this work, we focus on unsupervised continual learning (UCL), where we learn the feature representations on an unlabelled sequence of tasks and show that reliance on annotated data is not necessary for continual learning. We conduct a systematic study analyzing the learned feature representations and show that unsupervised visual representations are surprisingly more robust to catastrophic forgetting, consistently achieve better performance, and generalize better to out-of-distribution tasks than SCL. Furthermore, we find that UCL achieves a smoother loss landscape through qualitative analysis of the learned representations and learns meaningful feature representations. Additionally, we propose Lifelong Unsupervised Mixup (LUMP), a simple yet effective technique that interpolates between the current task and previous tasks' instances to alleviate catastrophic forgetting for unsupervised representations.


title: Cloud-Assisted Collaborative Road Information Discovery with Gaussian Process: Application to Road Profile Estimation

abstract: There is an increasing popularity in exploiting modern vehicles as mobile sensors to obtain important road information such as potholes, black ice and road profile. Availability of such information has been identified as a key enabler for next-generation vehicles with enhanced safety, efficiency, and comfort. However, existing road information discovery approaches have been predominately performed in a single-vehicle setting, which is inevitably susceptible to vehicle model uncertainty and measurement errors. To overcome these limitations, this paper presents a novel cloud-assisted collaborative estimation framework that can utilize multiple heterogeneous vehicles to iteratively enhance estimation performance. Specifically, each vehicle combines its onboard measurements with a cloud-based Gaussian process (GP), crowdsourced from prior participating vehicles as "pseudo-measurements", into a local estimator to refine the estimation. The resultant local onboard estimation is then sent back to the cloud to update the GP, where we utilize a noisy input GP (NIGP) method to explicitly handle uncertain GPS measurements. We employ the proposed framework to the application of collaborative road profile estimation. Promising results on extensive simulations and hardware-in-the-loop experiments show that the proposed collaborative estimation can significantly enhance estimation and iteratively improve the performance from vehicle to vehicle, despite vehicle heterogeneity, model uncertainty, and measurement noises.


title: WAFFLE: Weighted Averaging for Personalized Federated Learning

abstract: In federated learning, model personalization can be a very effective strategy to deal with heterogeneous training data across clients. We introduce WAFFLE (Weighted Averaging For Federated LEarning), a personalized collaborative machine learning algorithm that leverages stochastic control variates for faster convergence. WAFFLE uses the Euclidean distance between clients' updates to weigh their individual contributions and thus minimize the personalized model loss on the specific agent of interest. Through a series of experiments, we compare our new approach to two recent personalized federated learning methods--Weight Erosion and APFL--as well as two general FL methods--Federated Averaging and SCAFFOLD. Performance is evaluated using two categories of non-identical client data distributions--concept shift and label skew--on two image data sets (MNIST and CIFAR10). Our experiments demonstrate the comparative effectiveness of WAFFLE, as it achieves or improves accuracy with faster convergence.


title: Hyperspace Selections Avoiding Points

abstract: In this paper, we deal with a hyperspace selection problem in the setting of connected spaces. We present two solutions of this problem illustrating the difference between selections for the nonempty closed sets, and those for the at most two-point sets. In the first case, we obtain a characterisation of compact orderable spaces. In the latter case -- that of selections for at most two-point sets, the same selection property is equivalent to the existence of a ternary relation on the space, known as a cyclic order, and gives a characterisation of the so called weakly cyclically orderable spaces.


title: Output Space Entropy Search Framework for Multi-Objective Bayesian Optimization

abstract: We consider the problem of black-box multi-objective optimization (MOO) using expensive function evaluations (also referred to as experiments), where the goal is to approximate the true Pareto set of solutions by minimizing the total resource cost of experiments. For example, in hardware design optimization, we need to find the designs that trade-off performance, energy, and area overhead using expensive computational simulations. The key challenge is to select the sequence of experiments to uncover high-quality solutions using minimal resources. In this paper, we propose a general framework for solving MOO problems based on the principle of output space entropy (OSE) search: select the experiment that maximizes the information gained per unit resource cost about the true Pareto front. We appropriately instantiate the principle of OSE search to derive efficient algorithms for the following four MOO problem settings: 1) The most basic em single-fidelity setting, where experiments are expensive and accurate; 2) Handling em black-box constraints} which cannot be evaluated without performing experiments; 3) The discrete multi-fidelity setting, where experiments can vary in the amount of resources consumed and their evaluation accuracy; and 4) The em continuous-fidelity setting, where continuous function approximations result in a huge space of experiments. Experiments on diverse synthetic and real-world benchmarks show that our OSE search based algorithms improve over state-of-the-art methods in terms of both computational-efficiency and accuracy of MOO solutions.


title: FlexiTerm: A more efficient implementation of flexible multi-word term recognition

abstract: Terms are linguistic signifiers of domain-specific concepts. Automated recognition of multi-word terms (MWT) in free text is a sequence labelling problem, which is commonly addressed using supervised machine learning methods. Their need for manual annotation of training data makes it difficult to port such methods across domains. FlexiTerm, on the other hand, is a fully unsupervised method for MWT recognition from domain-specific corpora. Originally implemented in Java as a proof of concept, it did not scale well, thus offering little practical value in the context of big data. In this paper, we describe its re-implementation in Python and compare the performance of these two implementations. The results demonstrated major improvements in terms of efficiency, which allow FlexiTerm to transition from the proof of concept to the production-grade application.


title: Masking Effects in Combined Hardness and Stiffness Rendering Using an Encountered-Type Haptic Display

abstract: Rendering stable hard surfaces is an important problem in haptics for many tasks, including training simulators for orthopedic surgery or dentistry. Current impedance devices cannot provide enough force and stiffness to render a wall, and the high friction and inertia of admittance devices make it difficult to render free space. We propose to address these limitations by combining haptic augmented reality, untethered haptic interaction, and an encountered-type haptic display. We attach a plate with the desired hardness on the kinesthetic device's end-effector, which the user interacts with using an untethered stylus. This method allows us to directly change the hardness of the end-effector based on the rendered object. In this paper, we evaluate how changing the hardness of the end-effector can mask the device's stiffness and affect the user's perception. The results of our human subject experiment indicate that when the end-effector is made of a hard material, it is difficult for users to perceive when the underlying stiffness being rendered by the device is changed, but this stiffness change is easy to distinguish while the end-effector is made of a soft material. These results show promise for our approach in avoiding the limitations of haptic devices when rendering hard surfaces.


title: Bundle Networks: Fiber Bundles, Local Trivializations, and a Generative Approach to Exploring Many-to-one Maps

abstract: Many-to-one maps are ubiquitous in machine learning, from the image recognition model that assigns a multitude of distinct images to the concept of "cat" to the time series forecasting model which assigns a range of distinct time-series to a single scalar regression value. While the primary use of such models is naturally to associate correct output to each input, in many problems it is also useful to be able to explore, understand, and sample from a model's fibers, which are the set of input values $x$ such that $f(x) = y$, for fixed $y$ in the output space. In this paper we show that popular generative architectures are ill-suited to such tasks. Motivated by this we introduce a novel generative architecture, a Bundle Network, based on the concept of a fiber bundle from (differential) topology. BundleNets exploit the idea of a local trivialization wherein a space can be locally decomposed into a product space that cleanly encodes the many-to-one nature of the map. By enforcing this decomposition in BundleNets and by utilizing state-of-the-art invertible components, investigating a network's fibers becomes natural.


title: Cut norm discontinuity of triangular truncation of graphons

abstract: The space of $L^p$ graphons, symmetric measurable functions $w: [0,1]^2 \to \mathbb{R}$ with finite $p$-norm, features heavily in the study of sparse graph limit theory. We show that the triangular cut operator $M_{\chi}$ acting on this space is not continuous with respect to the cut norm. This is achieved by showing that as $n\to \infty$, the norm of the triangular truncation operator $\mathcal{T}_n$ on symmetric matrices equipped with the cut norm grows to infinity as well. Due to the density of symmetric matrices in the space of $L^p$ graphons, the norm growth of $\mathcal{T}_n$ generalizes to the unboundedness of $M_{\chi}$. We also show that the norm of $\mathcal{T}_n$ grows to infinity on symmetric matrices equipped with the operator norm.


title: Electroweak Symmetry Breaking and WIMP-FIMP Dark Matter

abstract: Electroweak Symmetry Breaking (EWSB) is known to produce a massive universe that we live in. However, it may also provide an important boundary for freeze-in or freeze-out of dark matter (DM) connected to Standard Model via Higgs portal as processes contributing to DM relic differ across the boundary. We explore such possibilities in a two-component DM framework, where a massive $U(1)_X$ gauge boson DM freezes-in and a scalar singlet DM freezes-out, that inherits the effect of EWSB for both the cases in a correlated way. Amongst different possibilities, we study two sample cases; first when one DM component freezes in and the other freezes out from thermal bath both necessarily $before$ EWSB and the second, when both freeze-in and freeze-out occur $after$ EWSB. We find some prominent distinctive features in the available parameter space of the model for these two cases, after addressing relic density and the recent most direct search constraints from XENON1T, some of which can be borrowed in a model independent way.


title: ADMM-DAD net: a deep unfolding network for analysis compressed sensing

abstract: In this paper, we propose a new deep unfolding neural network based on the ADMM algorithm for analysis Compressed Sensing. The proposed network jointly learns a redundant analysis operator for sparsification and reconstructs the signal of interest. We compare our proposed network with a state-of-the-art unfolded ISTA decoder, that also learns an orthogonal sparsifier. Moreover, we consider not only image, but also speech datasets as test examples. Computational experiments demonstrate that our proposed network outperforms the state-of-the-art deep unfolding network, consistently for both real-world image and speech datasets.


title: Scattering for the defocusing, cubic nonlinear Schr{\"o}dinger equation with initial data in a critical space

abstract: In this note we prove scattering for a defocusing nonlinear Schr{\"o}dinger equation with initial data lying in a critical Besov space. In addition, we obtain polynomial bounds on the scattering size as a function of the critical Besov norm.


title: Spectral Convergence of Symmetrized Graph Laplacian on manifolds with boundary

abstract: We study the spectral convergence of a symmetrized Graph Laplacian matrix induced by a Gaussian kernel evaluated on pairs of embedded data, sampled from a manifold with boundary, a sub-manifold of $\mathbb{R}^m$. Specifically, we deduce the convergence rates for eigenpairs of the discrete Graph-Laplacian matrix to the eigensolutions of the Laplace-Beltrami operator that are well-defined on manifolds with boundary, including the homogeneous Neumann and Dirichlet boundary conditions. For the Dirichlet problem, we deduce the convergence of the \emph{truncated Graph Laplacian}, which is recently numerically observed in applications, and provide a detailed numerical investigation on simple manifolds. Our method of proof relies on the min-max argument over a compact and symmetric integral operator, leveraging the RKHS theory for spectral convergence of integral operator and a recent pointwise asymptotic result of a Gaussian kernel integral operator on manifolds with boundary.


title: Pion and Kaon Distribution Amplitudes up to twist-3 in the QCD Instanton Vacuum

abstract: We discuss the pion and kaon distribution amplitudes up to twist-3 in the context of the random instanton vacuum (RIV). We construct explicitly the pertinent quasi-pion and quasi-kaon distributions in the RIV, and analyze them in leading order in the diluteness factor, at a resolution fixed by the inverse instanton size. The distribution amplitudes (DA) follow from the large momentum limit. The results at higher resolution are discussed using QCD evolution, and compared to their asymptotic limits and some lattice results.


title: Scaling Laws for the Few-Shot Adaptation of Pre-trained Image Classifiers

abstract: Empirical science of neural scaling laws is a rapidly growing area of significant importance to the future of machine learning, particularly in the light of recent breakthroughs achieved by large-scale pre-trained models such as GPT-3, CLIP and DALL-e. Accurately predicting the neural network performance with increasing resources such as data, compute and model size provides a more comprehensive evaluation of different approaches across multiple scales, as opposed to traditional point-wise comparisons of fixed-size models on fixed-size benchmarks, and, most importantly, allows for focus on the best-scaling, and thus most promising in the future, approaches. In this work, we consider a challenging problem of few-shot learning in image classification, especially when the target data distribution in the few-shot phase is different from the source, training, data distribution, in a sense that it includes new image classes not encountered during training. Our current main goal is to investigate how the amount of pre-training data affects the few-shot generalization performance of standard image classifiers. Our key observations are that (1) such performance improvements are well-approximated by power laws (linear log-log plots) as the training set size increases, (2) this applies to both cases of target data coming from either the same or from a different domain (i.e., new classes) as the training data, and (3) few-shot performance on new classes converges at a faster rate than the standard classification performance on previously seen classes. Our findings shed new light on the relationship between scale and generalization.


title: Scalable Graph Embedding LearningOn A Single GPU

abstract: Graph embedding techniques have attracted growing interest since they convert the graph data into continuous and low-dimensional space. Effective graph analytic provides users a deeper understanding of what is behind the data and thus can benefit a variety of machine learning tasks. With the current scale of real-world applications, most graph analytic methods suffer high computation and space costs. These methods and systems can process a network with thousands to a few million nodes. However, scaling to large-scale networks remains a challenge. The complexity of training graph embedding system requires the use of existing accelerators such as GPU. In this paper, we introduce a hybrid CPU-GPU framework that addresses the challenges of learning embedding of large-scale graphs. The performance of our method is compared qualitatively and quantitatively with the existing embedding systems on common benchmarks. We also show that our system can scale training to datasets with an order of magnitude greater than a single machine's total memory capacity. The effectiveness of the learned embedding is evaluated within multiple downstream applications. The experimental results indicate the effectiveness of the learned embedding in terms of performance and accuracy.


title: Convergence Rates of Decentralized Gradient Methods over Cluster Networks

abstract: We present an analysis for the performance of decentralized consensus-based gradient (DCG) methods for solving optimization problems over a cluster network of nodes. This type of network is composed of a number of densely connected clusters with a sparse connection between them. Decentralized algorithms over cluster networks have been observed to constitute two-time-scale dynamics, where information within any cluster is mixed much faster than the one across clusters. Based on this observation, we present a novel analysis to study the convergence of the DCG methods over cluster networks. In particular, we show that these methods converge at a rate $\ln(T)/T$ and only scale with the number of clusters, which is relatively small to the size of the network. Our result improves the existing analysis, where these methods are shown to scale with the size of the network. The key technique in our analysis is to consider a novel Lyapunov function that captures the impact of multiple time-scale dynamics on the convergence of this method. We also illustrate our theoretical results by a number of numerical simulations using DCG methods over different cluster networks.


title: Speed limits on classical chaos

abstract: Uncertainty in the initial conditions of dynamical systems can cause exponentially fast divergence of trajectories, a signature of deterministic chaos. Here, we derive a classical uncertainty relation that sets a speed limit on the rates of local observables underlying this behavior. For systems with a time-invariant stability matrix, this general speed limit simplifies to classical analogues of the Mandelstam-Tamm versions of the time-energy uncertainty relation. This classical bound derives from our definition of Fisher information in terms of Lyapunov vectors on tangent space, analogous to the quantum Fisher information defined in terms of wavevectors on Hilbert space. This information measures fluctuations in local stability of the state space and sets a lower bound on the time of classical, dynamical systems to evolve between two distinguishable states. The bounds it sets apply to systems that are open or closed, conservative or dissipative, actively driven or passively evolving, and directly connect the geometries of phase space and information.


title: Approximation of the Set of Integrable Trajectories of the Control Systems with Limited Control Resources

abstract: In this paper an approximation of the set of multivariable and $L_2$ integrable trajectories of the control system described by Urysohn type integral equation is considered. It is assumed that the system is affine with respect to the control vector. The admissible control functions are chosen from the closed ball of the space $L_2$, centered at the origin with radius $\rho$. Step by step way, the set of admissible control functions is replaced by the set of controls, consisting of a finite number of piecewise-constant control functions. It is proved that under appropriate choosing of the discretization parameters, the set of trajectories generated by a finite number of piecewise-constant control functions is an internal approximation of the set of trajectories.


title: Deep Chandra observations of diffuse hot plasma in M83

abstract: It is widely believed that galaxy formation and evolution is regulated by stellar mechanical feedback in forms of fast stellar winds and supernova explosions. However, the coupling of this feedback with the interstellar medium remains poorly understood. We examine how the coupling may be traced by diffuse soft X-ray emission in M83 -- a nearby face-on spiral galaxy undergoing active star formation, based chiefly on 729~ks Chandra observations. Our main findings are 1) the X-ray emission is enhanced not only along the galaxy's grand spiral arms, but also clearly in their downstreams; 2) the spectrum of the emission can be well characterized by a super-solar metallicity plasma with a lognormal temperature distribution, plus an X-ray absorption of a lognormal column density distribution; 3) the intensity of the emission is strongly anti-correlated with the dust obscuration seen in optical images of the galaxy. These findings suggest A) the morphology of the X-ray emission is likely due to the convolution of the feedback heating of the plasma with its thermal and dynamical evolution; B) the X-ray emission, accounting for ~10% of the feedback energy input rate, probably traces only the high-energy tail of the radiation from the plasma; C) a good fraction of the recent star forming regions seems sufficiently energetic to produce multi-phased outflows, likely responsible for much of the dust obscuration and X-ray absorption. Direct confrontation of the findings with theories/simulations could help to understand the underlying astrophysics of the coupling and how the hot plasma shapes the interstellar medium.


title: A Novel Clustering-Based Algorithm for Continuous and Non-invasive Cuff-Less Blood Pressure Estimation

abstract: Extensive research has been performed on continuous, non-invasive, cuffless blood pressure (BP) measurement using artificial intelligence algorithms. This approach involves extracting certain features from physiological signals like ECG, PPG, ICG, BCG, etc. as independent variables and extracting features from Arterial Blood Pressure (ABP) signals as dependent variables, and then using machine learning algorithms to develop a blood pressure estimation model based on these data. The greatest challenge of this field is the insufficient accuracy of estimation models. This paper proposes a novel blood pressure estimation method with a clustering step for accuracy improvement. The proposed method involves extracting Pulse Transit Time (PTT), PPG Intensity Ratio (PIR), and Heart Rate (HR) features from Electrocardiogram (ECG) and Photoplethysmogram (PPG) signals as the inputs of clustering and regression, extracting Systolic Blood Pressure (SBP) and Diastolic Blood Pressure (DBP) features from ABP signals as dependent variables, and finally developing regression models by applying Gradient Boosting Regression (GBR), Random Forest Regression (RFR), and Multilayer Perceptron Regression (MLP) on each cluster. The method was implemented using the MIMICII dataset with the silhouette criterion used to determine the optimal number of clusters. The results showed that because of the inconsistency, high dispersion, and multi-trend behavior of the extracted features vectors, the accuracy can be significantly improved by running a clustering algorithm and then developing a regression model on each cluster, and finally weighted averaging of the results based on the error of each cluster. When implemented with 5 clusters and GBR, this approach yielded an MAE of 2.56 for SBP estimates and 2.23 for DBP estimates, which were significantly better than the best results without clustering (DBP: 6.27, SBP: 6.36).


title: Bandits Don't Follow Rules: Balancing Multi-Facet Machine Translation with Multi-Armed Bandits

abstract: Training data for machine translation (MT) is often sourced from a multitude of large corpora that are multi-faceted in nature, e.g. containing contents from multiple domains or different levels of quality or complexity. Naturally, these facets do not occur with equal frequency, nor are they equally important for the test scenario at hand. In this work, we propose to optimize this balance jointly with MT model parameters to relieve system developers from manual schedule design. A multi-armed bandit is trained to dynamically choose between facets in a way that is most beneficial for the MT system. We evaluate it on three different multi-facet applications: balancing translationese and natural training data, or data from multiple domains or multiple language pairs. We find that bandit learning leads to competitive MT systems across tasks, and our analysis provides insights into its learned strategies and the underlying data sets.


title: Refining bridge-block decompositions through two-stage and recursive tree partitioning

abstract: In transmission networks power flows and network topology are deeply intertwined due to power flow physics. Recent literature shows that specific network substructures named bridge-blocks prevent line failures from propagating globally. A two-stage and recursive tree partitioning approach have been proposed to create more bridge-blocks in transmission networks, improving their robustness against cascading line failures. In this paper we consider the problem of refining the bridge-block decomposition of a given power network with minimal impact on the maximum congestion. We propose two new solution methods, depending on the preferred power flow model. More specifically, (i) we introduce a novel MILP-based approach that uses the DC approximation to solve more efficiently the second-stage optimization problem of the two-stage approach and (ii) we show how the existing recursive approach can be extended to work with AC power flows, drastically improving the running times when compared to the pre-existing AC-based two-stage method.


title: Study of positional encoding approaches for Audio Spectrogram Transformers

abstract: Transformers have revolutionized the world of deep learning, specially in the field of natural language processing. Recently, the Audio Spectrogram Transformer (AST) was proposed for audio classification, leading to state of the art results in several datasets. However, in order for ASTs to outperform CNNs, pretraining with ImageNet is needed. In this paper, we study one component of the AST, the positional encoding, and propose several variants to improve the performance of ASTs trained from scratch, without ImageNet pretraining. Our best model, which incorporates conditional positional encodings, significantly improves performance on Audioset and ESC-50 compared to the original AST.


title: Mixed-integer linear programming approaches for tree partitioning of power networks

abstract: In transmission networks, power flows and network topology are deeply intertwined due to power flow physics. Recent literature shows that a specific more hierarchical network structure can effectively inhibit the propagation of line failures across the entire system. In particular, a novel approach named tree partitioning has been proposed, which seeks to bolster the robustness of power networks through strategic alterations in network topology, accomplished via targeted line switching actions. Several tree partitioning problem formulations have been proposed by considering different objectives, among which power flow disruption and network congestion. Furthermore, various heuristic methods based on a two-stage and recursive approach have been proposed. The present work provides a general framework for tree partitioning problems based on mixed-integer linear programming (MILP). In particular, we present a novel MILP formulation to optimally solve tree partitioning problems and also propose a two-stage heuristic based on MILP. We perform extensive numerical experiments to solve two tree partitioning problem variants, demonstrating the excellent performance of our solution methods. Lastly, through exhaustive cascading failure simulations, we compare the effectiveness of various tree partitioning strategies and show that, on average, they can achieve a substantial reduction in lost load compared to the original topologies.


title: Higher theta series for unitary groups over function fields

abstract: In previous work, we defined certain virtual fundamental classes for special cycles on the moduli stack of Hermitian shtukas, and related them to the higher derivatives of non-singular Fourier coefficients of Siegel-Eisenstein series. In the present article, we construct virtual fundamental classes in greater generality, including those expected to relate to the higher derivatives of singular Fourier coefficients. We assemble these classes into "higher" theta series, which we conjecture to be modular. Two types of evidence are presented: structural properties affirming that the cycle classes behave as conjectured under certain natural operations such as intersection products, and verification of modularity in several special situations. One innovation underlying these results is a new approach to special cycles in terms of derived algebraic geometry.


title: Bag-of-Vectors Autoencoders for Unsupervised Conditional Text Generation

abstract: Text autoencoders are often used for unsupervised conditional text generation by applying mappings in the latent space to change attributes to the desired values. Recently, Mai et al. (2020) proposed Emb2Emb, a method to learn these mappings in the embedding space of an autoencoder. However, their method is restricted to autoencoders with a single-vector embedding, which limits how much information can be retained. We address this issue by extending their method to Bag-of-Vectors Autoencoders (BoV-AEs), which encode the text into a variable-size bag of vectors that grows with the size of the text, as in attention-based models. This allows to encode and reconstruct much longer texts than standard autoencoders. Analogous to conventional autoencoders, we propose regularization techniques that facilitate learning meaningful operations in the latent space. Finally, we adapt Emb2Emb for a training scheme that learns to map an input bag to an output bag, including a novel loss function and neural architecture. Our empirical evaluations on unsupervised sentiment transfer show that our method performs substantially better than a standard autoencoder.


title: Sustainability Through Cognition Aware Safety Systems -- Next Level Human-Machine-Interaction

abstract: Industrial Safety deals with the physical integrity of humans, machines and the environment when they interact during production scenarios. Industrial Safety is subject to a rigorous certification process that leads to inflexible settings, in which all changes are forbidden. With the progressing introduction of smart robotics and smart machinery to the factory floor, combined with an increasing shortage of skilled workers, it becomes imperative that safety scenarios incorporate a flexible handling of the boundary between humans, machines and the environment. In order to increase the well-being of workers, reduce accidents, and compensate for different skill sets, the configuration of machines and the factory floor should be dynamically adapted, while still enforcing functional safety requirements. The contribution of this paper is as follows: (1) We present a set of three scenarios, and discuss how industrial safety mechanisms could be augmented through dynamic changes to the work environment in order to decrease potential accidents, and thus increase productivity. (2) We introduce the concept of a Cognition Aware Safety System (CASS) and its architecture. The idea behind CASS is to integrate AI based reasoning about human load, stress, and attention with AI based selection of actions to avoid the triggering of safety stops. (3) And finally, we will describe the required performance measurement dimensions for a quantitative performance measurement model to enable a comprehensive (triple bottom line) impact assessment of CASS. Additionally we introduce a detailed guideline for expert interviews to explore the feasibility of the approach for given scenarios.


title: On the Convergence Theory for Hessian-Free Bilevel Algorithms

abstract: Bilevel optimization has arisen as a powerful tool in modern machine learning. However, due to the nested structure of bilevel optimization, even gradient-based methods require second-order derivative approximations via Jacobian- or/and Hessian-vector computations, which can be costly and unscalable in practice. Recently, Hessian-free bilevel schemes have been proposed to resolve this issue, where the general idea is to use zeroth- or first-order methods to approximate the full hypergradient of the bilevel problem. However, we empirically observe that such approximation can lead to large variance and unstable training, but estimating only the response Jacobian matrix as a partial component of the hypergradient turns out to be extremely effective. To this end, we propose a new Hessian-free method, which adopts the zeroth-order-like method to approximate the response Jacobian matrix via taking difference between two optimization paths. Theoretically, we provide the convergence rate analysis for the proposed algorithms, where our key challenge is to characterize the approximation and smoothness properties of the trajectory-dependent estimator, which can be of independent interest. This is the first known convergence rate result for this type of Hessian-free bilevel algorithms. Experimentally, we demonstrate that the proposed algorithms outperform baseline bilevel optimizers on various bilevel problems. Particularly, in our experiment on few-shot meta-learning with ResNet-12 network over the miniImageNet dataset, we show that our algorithm outperforms baseline meta-learning algorithms, while other baseline bilevel optimizers do not solve such meta-learning problems within a comparable time frame.


title: Proper orientations and proper chromatic number

abstract: The proper chromatic number $\Vec{\chi}(G)$ of a graph $G$ is the minimum $k$ such that there exists an orientation of the edges of $G$ with all vertex-outdegrees at most $k$ and such that for any adjacent vertices, the outdegrees are different. Two major conjectures about the proper chromatic number are resolved. First it is shown, that $\Vec{\chi}(G)$ of any planar graph $G$ is bounded (in fact, it is at most 14). Secondly, it is shown that for every graph, $\Vec{\chi}(G)$ is at most $O(\frac{r\log r}{\log\log r})+\tfrac{1}{2}\MAD(G)$, where $r=\chi(G)$ is the usual chromatic number of the graph, and $\MAD(G)$ is the maximum average degree taken over all subgraphs of $G$. Several other related results are derived. Our proofs are based on a novel notion of fractional orientations.


title: Estimating the effects of a California gun control program with Multitask Gaussian Processes

abstract: Gun violence is a critical public safety concern in the United States. In 2006 California implemented a unique firearm monitoring program, the Armed and Prohibited Persons System (APPS), to address gun violence in the state. The APPS program first identifies those firearm owners who become prohibited from owning one due to federal or state law, then confiscates their firearms. Our goal is to assess the effect of APPS on California murder rates using annual, state-level crime data across the US for the years before and after the introduction of the program. To do so, we adapt a non-parametric Bayesian approach, multitask Gaussian Processes (MTGPs), to the panel data setting. MTGPs allow for flexible and parsimonious panel data models that nest many existing approaches and allow for direct control over both dependence across time and dependence across units, as well as natural uncertainty quantification. We extend this approach to incorporate non-Normal outcomes, auxiliary covariates, and multiple outcome series, which are all important in our application. We also show that this approach has attractive Frequentist properties, including a representation as a weighting estimator with separate weights over units and time periods. Applying this approach, we find that the increased monitoring and enforcement from the APPS program substantially decreased homicides in California. We also find that the effect on murder is driven entirely by declines in gun-related murder with no measurable effect on non-gun murder. Estimated cost per murder avoided are substantially lower than conventional estimates of the value of a statistical life, suggesting a very high benefit-cost ratio for this enforcement effort.


title: Out-of-Distribution Robustness in Deep Learning Compression

abstract: In recent years, deep neural network (DNN) compression systems have proved to be highly effective for designing source codes for many natural sources. However, like many other machine learning systems, these compressors suffer from vulnerabilities to distribution shifts as well as out-of-distribution (OOD) data, which reduces their real-world applications. In this paper, we initiate the study of OOD robust compression. Considering robustness to two types of ambiguity sets (Wasserstein balls and group shifts), we propose algorithmic and architectural frameworks built on two principled methods: one that trains DNN compressors using distributionally-robust optimization (DRO), and the other which uses a structured latent code. Our results demonstrate that both methods enforce robustness compared to a standard DNN compressor, and that using a structured code can be superior to the DRO compressor. We observe tradeoffs between robustness and distortion and corroborate these findings theoretically for a specific class of sources.


title: Conductivity of Concentrated Electrolytes

abstract: The conductivity of ionic solutions is arguably their most important trait, being widely used in electrochemical, biochemical, and environmental applications. The Debye-H\"uckel-Onsager theory successfully predicts the conductivity at very low ionic concentrations of up to a few millimolars, but there is no well-established theory applicable at higher concentrations. We study the conductivity of ionic solutions using a stochastic density functional theory, paired with a modified Coulomb interaction that accounts for the hard-core repulsion between the ions. The modified potential suppresses unphysical, short-range electrostatic interactions, which are present in the Debye-H\"uckel-Onsager theory. Our results for the conductivity show very good agreement with experimental data up to 3 molars, without any fit parameters. We provide a compact expression for the conductivity, accompanied by a simple analytical approximation.


title: Leveraging Generative Models for Covert Messaging: Challenges and Tradeoffs for "Dead-Drop" Deployments

abstract: State of the art generative models of human-produced content are the focus of many recent papers that explore their use for steganographic communication. In particular, generative models of natural language text. Loosely, these works (invertibly) encode message-carrying bits into a sequence of samples from the model, ultimately yielding a plausible natural language covertext. By focusing on this narrow steganographic piece, prior work has largely ignored the significant algorithmic challenges, and performance-security tradeoffs, that arise when one actually tries to build a messaging pipeline around it. We make these challenges concrete, by considering the natural application of such a pipeline: namely, "dead-drop" covert messaging over large, public internet platforms (e.g. social media sites). We explicate the challenges and describe approaches to overcome them, surfacing in the process important performance and security tradeoffs that must be carefully tuned. We implement a system around this model-based format-transforming encryption pipeline, and give an empirical analysis of its performance and (heuristic) security.


title: Distributed and Localized Model Predictive Control. Part I: Synthesis and Implementation

abstract: The increasing presence of large-scale distributed systems highlights the need for scalable control strategies where only local communication is required. Moreover, in safety-critical systems it is imperative that such control strategies handle constraints in the presence of disturbances. In response to this need, we present the Distributed and Localized Model Predictive Control (DLMPC) algorithm for large-scale linear systems. DLMPC is a distributed closed-loop model predictive control (MPC) scheme wherein only local state and model information needs to be exchanged between subsystems for the computation and implementation of control actions. We use the System Level Synthesis (SLS) framework to reformulate the centralized MPC problem, and show that this allows us to naturally impose localized communication constraints between sub-controllers. The structure of the resulting problem can be exploited to develop an Alternating Direction Method of Multipliers (ADMM) based algorithm that allows for distributed and localized computation of closed-loop control policies. We demonstrate that computational complexity of the subproblems solved by each subsystem in DLMPC is independent of the size of the global system. To the best of our knowledge, DLMPC is the first MPC algorithm that allows for the scalable distributed computation as well as implementation of distributed closed-loop control policies, and seemingly deals with additive disturbances. In our companion paper, we show that this approach enjoys recursive feasibility and asymptotic stability.


title: Random Walks with Variable Restarts for Negative-Example-Informed Label Propagation

abstract: Label propagation is frequently encountered in machine learning and data mining applications on graphs, either as a standalone problem or as part of node classification. Many label propagation algorithms utilize random walks (or network propagation), which provide limited ability to take into account negatively-labeled nodes (i.e., nodes that are known to be not associated with the label of interest). Specialized algorithms to incorporate negatively labeled samples generally focus on learning or readjusting the edge weights to drive walks away from negatively-labeled nodes and toward positively-labeled nodes. This approach has several disadvantages, as it increases the number of parameters to be learned, and does not necessarily drive the walk away from regions of the network that are rich in negatively-labeled nodes. We reformulate random walk with restarts and network propagation to enable "variable restarts", that is the increased likelihood of restarting at a positively-labeled node when a negatively-labeled node is encountered. Based on this reformulation, we develop CusTaRd, an algorithm that effectively combines variable restart probabilities and edge re-weighting to avoid negatively-labeled nodes. In addition to allowing variable restarts, CusTaRd samples negatively-labeled nodes from neighbors of positively-labeled nodes to better characterize the difference between positively and negatively labeled nodes. To assess the performance of CusTaRd, we perform comprehensive experiments on four network datasets commonly used in benchmarking label propagation and node classification algorithms. Our results show that CusTaRd consistently outperforms competing algorithms that learn/readjust edge weights, and sampling of negatives from the close neighborhood of positives further improves predictive accuracy.


title: Structure of Stellar Remnants with Coupling to a Light Scalar

abstract: In this paper we study how a Yukawa coupling of the Standard Model fermions to a light scalar field effects the stellar structure of cold stellar remnants such as neutron stars. We elucidate the stellar structure phenomenology using a simple model of a massive scalar coupled to a single dominant fermion with no other interactions. For a broad scalar mass range ($10^{-10}\,\mathrm{eV}\ll m_\phi\ll10^3\,\mathrm{eV}$ for neutron stars) we show that the equation-of-state and stellar structure depends only the effective coupling $\mathfrak{g}=\frac{g_f\,m_f}{m_\phi}$, where $g_f$ is the Yukawa coupling, $m_f$ the fermion mass and $m_\phi$ is the scalar kinematic mass at nuclear densities. If $\mathfrak{g}>\mathcal{O}(1)$ the Yukawa coupled matter exhibits various anomalous behaviors including hydrodynamic instability, negative pressure, distinct phases (soft and hard) of matter with sharp phase boundaries between them and with vacuum. These anomalies can lead to stars consisting of only soft, only hard or hybrid of soft and hard matter. These stars can have varying signs of the slope of the mass-radius relation, anomalously large and small masses, gaps in allowed radii, multiple radii for the same mass, very thin crusts and radiate anomalously large amounts of energy when they form (in the form of neutrinos for neutron stars). To the extent that these anomalies have not and/or will not be observed limits the effective coupling to $\mathfrak{g}<\mathcal{O}(1)$. We argue this phenomenology is generic to realistic models of stars with Yukawa coupled matter.


title: Acousto-optically driven single-shot ultrafast optical imaging

abstract: Driven by many applications in a wide span of scientific fields, a myriad of advanced ultrafastimaging techniques have emerged in the last decade, featuring record-high imaging speeds abovea trillion-frame-per-second with long sequence depths. Although bringing remarkable insights invarious ultrafast phenomena, their application out of a laboratory environment is however limitedin most cases, either by the cost, complexity of operation or by an heavy data processing. We thenreport a flexible single-shot imaging technique combining sequentially-timed all-optical mappingphotography (STAMP) with acousto-optics programmable dispersive filtering. The full controlover the acquisition parameters is enabled via the spectro-temporal tailoring of the imaging pulsesin an electrically-driven spectral phase and amplitude shaper in which the pulse shaping in boththe temporal and spectral domains is controlled through the interaction of the light field with anacoustic wave. Here, contrary to most single-shot techniques, the frame rate, exposure time and frameintensities can be independently adjusted in a wide range of pulse durations and chirp values, makingthe system remarkably versatile and user-friendly. The imaging speed of the system as well as itsflexibility are validated by visualizing ultrashort events on both the picosecond and nanosecond timescales. With the perspective of real-world applications and to achieve the highest technical simplicity,we eventually demonstrate its lensless operation based on digital in-line holography. The virtues andlimitations as well as the potential improvements of this on-demand ultrafast imaging method arecritically discussed.


title: Impedance response of ionic liquids in long slit pores

abstract: We study the dynamics of ionic liquids in a thin slit pore geometry. Beginning with the classical and dynamic density functional theories for systems of charged hard spheres, an asymptotic procedure leads to a simplified model which incorporates both the accurate resolution of the ion layering (perpendicular to the slit pore wall) and the ion transport in the pore length. This reduced-order model enables qualitative comparisons between different ionic liquids and electrode pore sizes at low numerical expense. We derive semi-analytical expressions for the impedance response of the reduced-order model involving numerically computable sensitivities, and obtain effective finite-space Warburg elements valid in the high and low frequency limits. Additionally, we perform time-dependent numerical simulations to recover the impedance response as a cross-validation step. We investigate the dependence of the impedance response on system parameters and the choice of density functional theory used. The inclusion of electrostatic effects beyond mean-field qualitatively changes the dependence of the characteristic response time on the pore width. We observe peaks in the response time as a function of pore width, with height and location depending on the potential difference imposed. We discuss how the calculated dynamic properties can be used together with equilibrium results to optimise ionic liquid supercapacitors for a given application.


title: Weak and Strong Convergence of Generalized Proximal Point Algorithms with Relaxed Parameters

abstract: In this work, we propose and study a framework of generalized proximal point algorithms associated with a maximally monotone operator. We indicate sufficient conditions on the regularization and relaxation parameters of generalized proximal point algorithms for the equivalence of the boundedness of the sequence of iterations generated by this algorithm and the non-emptiness of the zero set of the maximally monotone operator, and for the weak and strong convergence of the algorithm. Our results cover or improve many results on generalized proximal point algorithms in our references. Improvements of our results are illustrated by comparing our results with related known ones.


title: Light, delayed: The Shapiro Effect and the Newtonian Limit

abstract: The Shapiro effect, also known as the gravitational time delay, is close kin to the gravitational deflection of light that was the central topic of our Summer School. It is also an interesting test bed for exploring a topic that provides the foundations for most of the calculations we have done in this school, yet is highly complex when treated more rigorously: the question of the Newtonian limit, and of the post-Newtonian corrections that must be applied to include the leading-order effects of general relativity. This contribution discusses simplified derivations for the gravitational redshift and the Shapiro effect, as well as astrophysical situations in which the Shapiro effect can be measured.


title: Unconditional uniqueness for the Benjamin-Ono equation

abstract: We study the unconditional uniqueness of solutions to the Benjamin-Ono equation with initial data in $H^{s}$, both on the real line and on the torus. We use the gauge transformation of Tao and two iterations of normal form reductions via integration by parts in time. By employing a refined Strichartz estimate we establish the result below the regularity threshold $s=1/6$. As a by-product of our proof, we also obtain a nonlinear smoothing property on the gauge variable at the same level of regularity.


title: Algebraic Reasoning of Quantum Programs via Non-idempotent Kleene Algebra

abstract: We investigate the algebraic reasoning of quantum programs inspired by the success of classical program analysis based on Kleene algebra. One prominent example of such is the famous Kleene Algebra with Tests (KAT), which has furnished both theoretical insights and practical tools. The succinctness of algebraic reasoning would be especially desirable for scalable analysis of quantum programs, given the involvement of exponential-size matrices in most of the existing methods. A few key features of KAT including the idempotent law and the nice properties of classical tests, however, fail to hold in the context of quantum programs due to their unique quantum features, especially in branching. We propose Non-idempotent Kleene Algebra (NKA) as a natural alternative and identify complete and sound semantic models for NKA as well as their quantum interpretations. In light of applications of KAT, we demonstrate algebraic proofs in NKA of quantum compiler optimization and the normal form of quantum while-programs. Moreover, we extend NKA with Tests (i.e., NKAT), where tests model quantum predicates following effect algebra, and illustrate how to encode propositional quantum Hoare logic as NKAT theorems.