title: Homogenization and Convergence Rates for Periodic Parabolic Equations with Highly Oscillating Potentials

abstract: This paper considers a family of second-order periodic parabolic equations with highly oscillating potentials, which have been considered many times for the time-varying potentials in stochastic homogenization. Following a standard two-scale expansions illusion, we can guess and succeed in determining the homogenized equation in different cases that the potentials satisfy the corresponding assumptions, based on suitable uniform estimates of the $L^2(0,T;H^1(\Omega))$-norm for the solutions. To handle the more singular case and obtain the convergence rates in $L^\infty(0,T;L^2(\Omega))$, we need to estimate the Hessian term as well as the t-derivative term more exactly, which may be depend on $\varepsilon$. The difficulty is to find suitable uniform estimates for the $L^2(0,T;H^1(\Omega))$-norm and suitable estimates for the higher order derivative terms.


title: Stochastic Composition Optimization of Functions without Lipschitz Continuous Gradient

abstract: In this paper, we study stochastic optimization of two-level composition of functions without Lipschitz continuous gradient. The smoothness property is generalized by the notion of relative smoothness which provokes the Bregman gradient method. We propose three Stochastic Composition Bregman Gradient algorithms for the three possible relatively smooth compositional scenarios and provide their sample complexities to achieve an $\epsilon$-approximate stationary point. For the smooth of relatively smooth composition, the first algorithm requires $O(\epsilon^{-2})$ calls to the stochastic oracles of the inner function value and gradient as well as the outer function gradient. When both functions are relatively smooth, the second algorithm requires $O(\epsilon^{-3})$ calls to the inner function value stochastic oracle and $O(\epsilon^{-2})$ calls to the inner and outer functions gradients stochastic oracles. We further improve the second algorithm by variance reduction for the setting where just the inner function is smooth. The resulting algorithm requires $O(\epsilon^{-5/2})$ calls to the inner function value stochastic oracle, $O(\epsilon^{-3/2})$ calls to the inner function gradient and $O(\epsilon^{-2})$ calls to the outer function gradient stochastic oracles. Finally, we numerically evaluate the performance of these three algorithms over two different examples.


title: Heavy meson thresholds in Born-Oppenheimer Effective field theory

abstract: We consider heavy meson-antimeson pairs and their coupling to quarkonium in the context of nonrelativistic EFTs incorporating the adiabatic expansion. We work out all the leading order couplings of quarkonium to heavy meson-antimeson pairs and obtain their contributions to the masses and widths of quarkonia. We match the new potentials terms to NRQCD. Using the available lattice data for the coupled system of quarkonium and the lowest laying heavy meson-antimeson pair, we extract the mixing potential and use it to compute numerically the contributions of $D\bar{D}(B\bar{B})$ and $D_s\bar{D}_s(B_s\bar{B}_s)$ to the masses and widths of the charmonium (bottomonium) states for $l=0,1,2$ and up to $n=6$ covering the states in threshold region. When a quarkonium state and a heavy meson-antimeson pair are separated by small energy gaps, their interactions can be described by a threshold EFT with contact interactions. We work out the matching between the two EFTs obtaining the couplings of the threshold EFT in terms of the mixing potential and quarkonium wave functions.


title: Constraints on Nonrelativistic-QCD Long-Distance Matrix Elements from $J/\psi$ Plus $W$/$Z$ Production at the LHC

abstract: We study the associated production of prompt $J/\psi$ mesons and $W$ or $Z$ bosons within the factorization approach of nonrelativistic QCD (NRQCD) at next-to-leading order in $\alpha_s$, via intermediate color singlet ${^3}S_1^{[1]}$ and ${^3}P_J^{[1]}$ and color octet ${^1S}_0^{[8]}$, ${^3S}_1^{[8]}$ and ${^3P}_J^{[8]}$ states. Requiring for our predictions to be compatible with recent ATLAS measurements yields stringent new constraints on charmonium long-distance matrix elements (LDMEs) being nonperturbative, process-independent input parameters. Considering four popular LDME sets fitted to data of single $J/\psi$ inclusive production, we find that one is marginally compatible with the data, with central predictions typically falling short by a factor of three, one is unfavored, the factor of shortfall being about one order of magnitude, and two violate cross section positivity for direct $J/\psi+W/Z$ production. The large rate of prompt $J/\psi$ plus $W$ production observed by ATLAS provides strong evidence for the color octet mechanism inherent to NRQCD factorization, the leading color singlet contribution entering only at $\mathcal{O}(G_F\alpha_s^4)$, beyond the order considered here.


title: Cycle Encoding of a StyleGAN Encoder for Improved Reconstruction and Editability

abstract: GAN inversion aims to invert an input image into the latent space of a pre-trained GAN. Despite the recent advances in GAN inversion, there remain challenges to mitigate the tradeoff between distortion and editability, i.e. reconstructing the input image accurately and editing the inverted image with a small visual quality drop. The recently proposed pivotal tuning model makes significant progress towards reconstruction and editability, by using a two-step approach that first inverts the input image into a latent code, called pivot code, and then alters the generator so that the input image can be accurately mapped into the pivot code. Here, we show that both reconstruction and editability can be improved by a proper design of the pivot code. We present a simple yet effective method, named cycle encoding, for a high-quality pivot code. The key idea of our method is to progressively train an encoder in varying spaces according to a cycle scheme: W->W+->W. This training methodology preserves the properties of both W and W+ spaces, i.e. high editability of W and low distortion of W+. To further decrease the distortion, we also propose to refine the pivot code with an optimization-based method, where a regularization term is introduced to reduce the degradation in editability. Qualitative and quantitative comparisons to several state-of-the-art methods demonstrate the superiority of our approach.


title: Characterising Determinism in MBQCs involving Pauli Measurements

abstract: We introduce a new characterisation of determinism in measurement-based quantum computing. The one-way model of computation consists in performing local measurements over a large entangled state represented by a graph. The ability to perform an overall deterministic computation requires a correction strategy because of the non-determinism of each measurement. The existence of such correction strategy depends on the underlying graph and the basis of the performed measurements. GFlow is a well-known graphical characterisation of robust determinism in MBQC when every measurement is performed in some specific planes of the Bloch sphere. While Pauli measurements are ubiquitous in MBQC, GFlow fails to be necessary for determinism when a measurement-based quantum computation involves Pauli measurements. As a consequence, Pauli Flow was designed more than 15 years ago as a generalisation of GFlow to handle MBQC with Pauli measurements: Pauli flow guarantees robust determinism, however it has been shown more recently that it fails to be a necessary condition. We introduce a further extension called Extended Pauli Flow that we prove necessary and sufficient for robust determinism.


title: Discovery of fulvenallene in TMC-1 with the QUIJOTE line survey

abstract: We report the detection of fulvenallene ($c$-C$_5$H$_4$CCH$_2$) in the direction of TMC-1 with the QUIJOTE line survey. Thirty rotational transitions with $K_a$=0,1,2,3 and $J$=9-15 were detected. The best rotational temperature fitting of the data is 9\,K and a derived column density is (2.7$\pm$0.3)$\times$10$^{12}$ cm$^{-2}$, which is only a factor of 4.4 below that of its potential precursor cyclopentadiene ($c$-C$_5$H$_6$), and 1.4--1.9 times higher than that of the ethynyl derivatives of cyclopentadiene. We searched for fulvene ($c$-C$_5$H$_4$CH$_2$), a CH$_2$ derivative of cyclopentadiene, for which we derive a 3$\sigma$ upper limit to its column density of (3.5$\pm$0.5)$\times$10$^{12}$ cm$^{-2}$. Upper limits were also obtained for toluene (C$_6$H$_5$CH$_3$) and styrene (C$_6$H$_5$C$_2$H$_3$), the methyl and vinyl derivatives of benzene. Fulvenallene and ethynyl cyclopentadiene are likely formed in the reaction between cyclopentadiene ($c$-C$_5$H$_6$) and the ehtynyl radical (CCH). However, the bottom-up gas-phase synthesis of cycles in TMC-1 underestimates the abundance of cyclopentadiene by two orders of magnitude, which strengthens the need to study all possible chemical pathways to cyclisation in cold dark cloud environments, such as TMC-1. However, the inclusion of the reaction between C$_3$H$_3^+$ and C$_2$H$_4$ produces a good agreement between model and observed abundances.


title: Data-Centric Epidemic Forecasting: A Survey

abstract: The COVID-19 pandemic has brought forth the importance of epidemic forecasting for decision makers in multiple domains, ranging from public health to the economy as a whole. While forecasting epidemic progression is frequently conceptualized as being analogous to weather forecasting, however it has some key differences and remains a non-trivial task. The spread of diseases is subject to multiple confounding factors spanning human behavior, pathogen dynamics, weather and environmental conditions. Research interest has been fueled by the increased availability of rich data sources capturing previously unobservable facets and also due to initiatives from government public health and funding agencies. This has resulted, in particular, in a spate of work on 'data-centered' solutions which have shown potential in enhancing our forecasting capabilities by leveraging non-traditional data sources as well as recent innovations in AI and machine learning. This survey delves into various data-driven methodological and practical advancements and introduces a conceptual framework to navigate through them. First, we enumerate the large number of epidemiological datasets and novel data streams that are relevant to epidemic forecasting, capturing various factors like symptomatic online surveys, retail and commerce, mobility, genomics data and more. Next, we discuss methods and modeling paradigms focusing on the recent data-driven statistical and deep-learning based methods as well as on the novel class of hybrid models that combine domain knowledge of mechanistic models with the effectiveness and flexibility of statistical approaches. We also discuss experiences and challenges that arise in real-world deployment of these forecasting systems including decision-making informed by forecasts. Finally, we highlight some challenges and open problems found across the forecasting pipeline.


title: Polynomial Threshold Functions for Decision Lists

abstract: For $S \subseteq \{0,1\}^n$ a Boolean function $f \colon S \to \{-1,1\}$ is a polynomial threshold function (PTF) of degree $d$ and weight $W$ if there is a polynomial $p$ with integer coefficients of degree $d$ and with sum of absolute coefficients $W$ such that $f(x) = \text{sign}(p(x))$ for all $x \in S$. We study a representation of decision lists as PTFs over Boolean cubes $\{0,1\}^n$ and over Hamming balls $\{0,1\}^{n}_{\leq k}$. As our first result, we show that for all $d = O\left( \left( \frac{n}{\log n}\right)^{1/3}\right)$ any decision list over $\{0,1\}^n$ can be represented by a PTF of degree $d$ and weight $2^{O(n/d^2)}$. This improves the result by Klivans and Servedio [Klivans, Servedio, 2006] by a $\log^2 d$ factor in the exponent of the weight. Our bound is tight for all $d = O\left( \left( \frac{n}{\log n}\right)^{1/3}\right)$ due to the matching lower bound by Beigel [Beigel, 1994]. For decision lists over a Hamming ball $\{0,1\}^n_{\leq k}$ we show that the upper bound on weight above can be drastically improved to $n^{O(\sqrt{k})}$ for $d = \Theta(\sqrt{k})$. We also show that similar improvement is not possible for smaller degrees by proving the lower bound $W = 2^{\Omega(n/d^2)}$ for all $d = O(\sqrt{k})$. \end{abstract}


title: On Decentralizing Federated Reinforcement Learning in Multi-Robot Scenarios

abstract: Federated Learning (FL) allows for collaboratively aggregating learned information across several computing devices and sharing the same amongst them, thereby tackling issues of privacy and the need of huge bandwidth. FL techniques generally use a central server or cloud for aggregating the models received from the devices. Such centralized FL techniques suffer from inherent problems such as failure of the central node and bottlenecks in channel bandwidth. When FL is used in conjunction with connected robots serving as devices, a failure of the central controlling entity can lead to a chaotic situation. This paper describes a mobile agent based paradigm to decentralize FL in multi-robot scenarios. Using Webots, a popular free open-source robot simulator, and Tartarus, a mobile agent platform, we present a methodology to decentralize federated learning in a set of connected robots. With Webots running on different connected computing systems, we show how mobile agents can perform the task of Decentralized Federated Reinforcement Learning (dFRL). Results obtained from experiments carried out using Q-learning and SARSA by aggregating their corresponding Q-tables, show the viability of using decentralized FL in the domain of robotics. Since the proposed work can be used in conjunction with other learning algorithms and also real robots, it can act as a vital tool for the study of decentralized FL using heterogeneous learning algorithms concurrently in multi-robot scenarios.


title: Multi-Task Learning Framework for Emotion Recognition in-the-wild

abstract: This paper presents our system for the Multi-Task Learning (MTL) Challenge in the 4th Affective Behavior Analysis in-the-wild (ABAW) competition. We explore the research problems of this challenge from three aspects: 1) For obtaining efficient and robust visual feature representations, we propose MAE-based unsupervised representation learning and IResNet/DenseNet-based supervised representation learning methods; 2) Considering the importance of temporal information in videos, we explore three types of sequential encoders to capture the temporal information, including the encoder based on transformer, the encoder based on LSTM, and the encoder based on GRU; 3) For modeling the correlation between these different tasks (i.e., valence, arousal, expression, and AU) for multi-task affective analysis, we first explore the dependency between these different tasks and propose three multi-task learning frameworks to model the correlations effectively. Our system achieves the performance of $1.7607$ on the validation dataset and $1.4361$ on the test dataset, ranking first in the MTL Challenge. The code is available at https://github.com/AIM3-RUC/ABAW4.


title: Alterfactual Explanations -- The Relevance of Irrelevance for Explaining AI Systems

abstract: Explanation mechanisms from the field of Counterfactual Thinking are a widely-used paradigm for Explainable Artificial Intelligence (XAI), as they follow a natural way of reasoning that humans are familiar with. However, all common approaches from this field are based on communicating information about features or characteristics that are especially important for an AI's decision. We argue that in order to fully understand a decision, not only knowledge about relevant features is needed, but that the awareness of irrelevant information also highly contributes to the creation of a user's mental model of an AI system. Therefore, we introduce a new way of explaining AI systems. Our approach, which we call Alterfactual Explanations, is based on showing an alternative reality where irrelevant features of an AI's input are altered. By doing so, the user directly sees which characteristics of the input data can change arbitrarily without influencing the AI's decision. We evaluate our approach in an extensive user study, revealing that it is able to significantly contribute to the participants' understanding of an AI. We show that alterfactual explanations are suited to convey an understanding of different aspects of the AI's reasoning than established counterfactual explanation methods.


title: Latest observations on the low energy excess in CRESST-III

abstract: The CRESST experiment observes an unexplained excess of events at low energies. In the current CRESST-III data-taking campaign we are operating detector modules with different designs to narrow down the possible explanations. In this work, we show first observations of the ongoing measurement, focusing on the comparison of time, energy and temperature dependence of the excess in several detectors. These exclude dark matter, radioactive backgrounds and intrinsic sources related to the crystal bulk as a major contribution.


title: Synchronization in pairs of opto-thermally driven mechanically coupled micro-oscillators

abstract: We study the phenomenon of synchronization in pairs of doubly clamped, mechanically coupled silicon micro-oscillators. A continuous-wave laser beam is used to drive the micro-beams into limit cycle oscillations and to detect the oscillations using interferometry. Devices of different dimensions are used to introduce frequency detuning, and short silicon bridges connecting the micro-beams are used as mechanical coupling between the oscillators. The region of synchronization is plotted for the MEMS system in the detuning vs. coupling parameter space and compared with the numerical analysis of a corresponding, lumped-parameter model. Three states of oscillations are observed i.e. the drift state, quasi-periodic state, and the synchronized state. The numerical model also distinguishes between in-phase and out-of-phase synchronization where out-of-phase synchronization is observed at low coupling strengths and low frequency detuning.


title: Solutions of the Schr\"{o}dinger equation for anisotropic dipole-dipole interaction plus isotropic van der Waals interaction

abstract: By generalizing Bo Gao's approach [Phys. Rev. A 58, 1728 (1998)] for solving the Schr\"{o}dinger equation for an isotropic van der Waals (vdW) potential to the systems with a multi-scale anisotropic long-range interaction, we derive the solutions for the Schr\"{o}dinger equation for an anisotropic dipole-dipole interaction plus an isotropic attractive vdW potential, i.e., ${C_d(1-3\cos^2\theta)}/{r^3}-{C_6}/{r^6}$, which is projected to the subspace with angular momentum $l\leq l_{\rm cut}$, with $l_{\rm cut}$ being an arbitrary angular-momentum cutoff. Here $\theta$ is the polar angle of the coordinate $\boldsymbol{r}$ and $r=|\boldsymbol{r}|$. The asymptotic behaviors of these solutions for $r\rightarrow 0$ and $r\rightarrow \infty$ are obtained. These results can be used in the research of collisions and chemical reactions between ultra-cold polar molecules in a static electric field. Our approach to derive the solutions can be applied to the systems with a general long-range potential $\sum_{\lambda= 2}^{\lambda_{\rm max}} {V_\lambda(\theta,\varphi)}/{r^\lambda}$, with $\varphi$ being the azimuthal angle of $\boldsymbol{r}$, and thus can be used in various problems on molecule-molecule interaction.

categories: cond-mat.quant-gas physics.atom-ph physics.chem-ph quant-ph


title: P4TE: PISA Switch Based Traffic Engineering in Fat-Tree Data Center Networks

abstract: This work presents P4TE, an in-band traffic monitoring, load-aware packet forwarding, and flow rate controlling mechanism for traffic engineering in fat-tree topology-based data center networks using PISA switches. It achieves sub-RTT reaction time to change in network conditions, improved flow completion time, and balanced link utilization. Unlike the classical probe-based monitoring approach, P4TE uses an in-band monitoring approach to identify traffic events in the data plane. Based on these events, it re-adjusts the priorities of the paths. It uses a heuristic-based load-aware forwarding path selection mechanism to respond to changing network conditions and control the flow rate by sending feedback to the end hosts. It is implementable on emerging v1model.p4 architecture-based programmable switches and capable of maintaining the line-rate performance. Our evaluation shows that P4TE uses a small amount of resources in the PISA pipeline and achieves an improved flow completion time than ECMP and HULA.


title: To what extent can we analyze Kotlin programs using existing Java taint analysis tools? (Extended Version)

abstract: As an alternative to Java, Kotlin has gained rapid popularity since its introduction and has become the default choice for developing Android apps. However, due to its interoperability with Java, Kotlin programs may contain almost the same security vulnerabilities as their Java counterparts. Hence, we question: to what extent can one use an existing Java static taint analysis on Kotlin code? In this paper, we investigate the challenges in implementing a taint analysis for Kotlin compared to Java. To answer this question, we performed an exploratory study where each Kotlin construct was examined and compared to its Java equivalent. We identified 18 engineering challenges that static-analysis writers need to handle differently due to Kotlin's unique constructs or the differences in the generated bytecode between the Kotlin and Java compilers. For eight of them, we provide a conceptual solution, while six of those we implemented as part of SecuCheck-Kotlin, an extension to the existing Java taint analysis SecuCheck.


title: Broken axial symmetry as essential feature for a consistent modelling of various observables in heavy nuclei

abstract: Although most nuclear spectroscopy as well as atomic hyperfine structure data do not deliver accurate information on nuclear axiality the ad-hoc assumption of symmetry about one axis found widespread use in nuclear model calculations. In the theoretical interpretation of nuclear properties as well as in the analysis of experimental data triaxiality was considered - if at all - only for some, often exotic, nuclides. A breaking of axial symmetry combined to a spin-independent moment of inertia results in a surprisingly simple heuristic triaxial parametrization of the yrast sequence in all heavy nuclei, including well deformed ones. No additional fit parameters are needed in detailed studies of the mass and charge dependence of the electric dipole strength in the range of and outside of giant dipole resonances. Allowing triaxiality also avoids the introduction of an arbitrary level density parameter \=a to fit the accurate values observed in n-capture experiments and \=a can be taken from nuclear matter studies. A combination of this value to the yrast energies no longer based on axiality and the related I(I+1) rule results in agreement to data independent of spin. And predictions for radiative neutron capture as derived on the basis of non-axiality are improved as well. The experimentally favoured broken axial symmetry is in accord to HFB and MC-shell model calculations already for nuclei in the valley of stability.


title: Evidence of dual Shapiro steps in a Josephson junctions array

abstract: The modern primary voltage standard is based on the AC Josephson effect and the ensuing Shapiro steps, where a microwave tone applied to a Josephson junction yields a constant voltage $hf/2e$ ($h$ is Planck's constant and $e$ the electron charge) determined by only the microwave frequency $f$ and fundamental constants. Duality arguments for current and voltage have long suggested the possibility of dual Shapiro steps -- that a Josephson junction device could produce current steps with heights determined only on the applied frequency. In this report, we embed an ultrasmall Josephson junction in a high impedance array of larger junctions to reveal dual Shapiro steps. For multiple frequencies, we detect that the AC response of the circuit is synchronised with the microwave tone at frequency $f$, and the corresponding emergence of flat steps in the DC response with current $2ef$, equal to the tunnelling of a Cooper pair per tone period. This work sheds new light on phase-charge duality, omnipresent in condensed matter physics, and extends it to Josephson circuits. Looking forward, it opens a broad range of possibilities for new experiments in the field of circuit quantum electrodynamics and is an important step towards the long-sought closure of the quantum metrology electrical triangle.


title: Inference for high-dimensional split-plot designs with different dimensions between groups

abstract: In repeated Measure Designs with multiple groups, the primary purpose is to compare different groups in various aspects. For several reasons, the number of measurements and therefore the dimension of the observation vectors can depend on the group, making the usage of existing approaches impossible. We develop an approach which can be used not only for a possibly increasing number of groups $a$, but also for group-depending dimension $d_i$, which is allowed to go to infinity. This is a unique high-dimensional asymptotic framework impressing through its variety and do without usual conditions on the relation between sample size and dimension. It especially includes settings with fixed dimensions in some groups and increasing dimensions in other ones, which can be seen as semi-high-dimensional. To find a appropriate statistic test new and innovative estimators are developed, which can be used under these diverse settings on $a,d_i$ and $n_i$ without any adjustments. We investigated the asymptotic distribution of a quadratic-form-based test statistic and developed an asymptotic correct test. Finally, an extensive simulation study is conducted to investigate the role of the single group's dimension.


title: A subwavelength atomic array switched by a single Rydberg atom

abstract: Enhancing light-matter coupling at the level of single quanta is essential for numerous applications in quantum science. The cooperative optical response of subwavelength atomic arrays has been found to open new pathways for such strong light-matter couplings, while simultaneously offering access to multiple spatial modes of the light field. Efficient single-mode free-space coupling to such arrays has been reported, but the spatial control over the modes of outgoing light fields has remained elusive. Here, we demonstrate such spatial control over the optical response of an atomically thin mirror formed by a subwavelength array of atoms in free space using a single controlled ancilla atom excited to a Rydberg state. The switching behavior is controlled by the admixture of a small Rydberg fraction to the atomic mirror, and consequently strong dipolar Rydberg interactions with the ancilla. Driving Rabi oscillations on the ancilla atom, we demonstrate coherent control of the transmission and reflection of the array. These results represent a step towards the realization of quantum coherent metasurfaces, the demonstration of controlled atom-photon entanglement and deterministic engineering of quantum states of light.


title: Scalable Spatio-Temporal Smoothing via Hierarchical Sparse Cholesky Decomposition

abstract: We propose an approximation to the forward-filter-backward-sampler (FFBS) algorithm for large-scale spatio-temporal smoothing. FFBS is commonly used in Bayesian statistics when working with linear Gaussian state-space models, but it requires inverting covariance matrices which have the size of the latent state vector. The computational burden associated with this operation effectively prohibits its applications in high-dimensional settings. We propose a scalable spatio-temporal FFBS approach based on the hierarchical Vecchia approximation of Gaussian processes, which has been previously successfully used in spatial statistics. On simulated and real data, our approach outperformed a low-rank FFBS approximation.


title: A Physical-Constraint-Preserving Finite Volume WENO Method for Special Relativistic Hydrodynamics on Unstructured Meshes

abstract: This paper presents a highly robust third-order accurate finite volume weighted essentially non-oscillatory (WENO) method for special relativistic hydrodynamics on unstructured triangular meshes. We rigorously prove that the proposed method is physical-constraint-preserving (PCP), namely, always preserves the positivity of the pressure and the rest-mass density as well as the subluminal constraint on the fluid velocity. The method is built on a highly efficient compact WENO reconstruction on unstructured meshes, a simple PCP limiter, the provably PCP property of the Harten--Lax--van Leer flux, and third-order strong-stability-preserving time discretization. Due to the relativistic effects, the primitive variables (namely, the rest-mass density, velocity, and pressure) are highly nonlinear implicit functions in terms of the conservative variables, making the design and analysis of our method nontrivial. To address the difficulties arising from the strong nonlinearity, we adopt a novel quasilinear technique for the theoretical proof of the PCP property. Three provable convergence-guaranteed iterative algorithms are also introduced for the robust recovery of primitive quantities from admissible conservative variables. We also propose a slight modification to an existing WENO reconstruction to ensure the scaling invariance of the nonlinear weights and thus to accommodate the homogeneity of the evolution operator, leading to the advantages of the modified WENO reconstruction in resolving multi-scale wave structures. Extensive numerical examples are presented to demonstrate the robustness, expected accuracy, and high resolution of the proposed method.

categories: math.NA astro-ph.IM cs.NA physics.comp-ph physics.flu-dyn


title: Inertia-gravity-wave diffusion by geostrophic turbulence: the impact of flow time dependence

abstract: The scattering of three-dimensional inertia-gravity waves by a turbulent geostrophic flow leads to the redistribution of their action through what is approximately a diffusion process in wavevector space. The corresponding diffusivity tensor was obtained by Kafiabad, Savva & Vanneste (2019, J. Fluid Mech., 869, R7) under the assumption of a time-independent geostrophic flow. We relax this assumption to examine how the weak diffusion of wave action across constant-frequency cones that results from the slow time dependence of the geostrophic flow affects the distribution of wave energy. We find that the stationary wave-energy spectrum that arises from a single-frequency wave forcing is localised within a thin boundary layer around the constant-frequency cone, with a thickness controlled by the acceleration spectrum of the geostrophic flow. We obtain an explicit analytic formula for the wave-energy spectrum which shows good agreement with the results of a high-resolution simulation of the Boussinesq equations.


title: Green, Quantized Federated Learning over Wireless Networks: An Energy-Efficient Design

abstract: In this paper, a green-quantized FL framework, which represents data with a finite precision level in both local training and uplink transmission, is proposed. Here, the finite precision level is captured through the use of quantized neural networks (QNNs) that quantize weights and activations in fixed-precision format. In the considered FL model, each device trains its QNN and transmits a quantized training result to the base station. Energy models for the local training and the transmission with quantization are rigorously derived. To minimize the energy consumption and the number of communication rounds simultaneously, a multi-objective optimization problem is formulated with respect to the number of local iterations, the number of selected devices, and the precision levels for both local training and transmission while ensuring convergence under a target accuracy constraint. To solve this problem, the convergence rate of the proposed FL system is analytically derived with respect to the system control variables. Then, the Pareto boundary of the problem is characterized to provide efficient solutions using the normal boundary inspection method. Design insights on balancing the tradeoff between the two objectives while achieving a target accuracy are drawn from using the Nash bargaining solution and analyzing the derived convergence rate. Simulation results show that the proposed FL framework can reduce energy consumption until convergence by up to 70\% compared to a baseline FL algorithm that represents data with full precision without damaging the convergence rate.


title: Hybrid photon-phonon blockade

abstract: We describe a novel type of blockade in a hybrid mode generated by linear coupling of photonic and phononic modes. We refer to this effect as hybrid photon-phonon blockade and show how it can be generated and detected in a driven nonlinear optomechanical superconducting system. Thus, we study boson-number correlations in the photon, phonon, and hybrid modes in linearly coupled microwave and mechanical resonators with a superconducting qubit inserted in one of them. We find such system parameters for which we observe eight types of different combinations of either blockade or tunnelling effects (defined via the sub- and super-Poissonian statistics, respectively) for photons, phonons, and hybrid bosons. In particular, we find that the hybrid photon-phonon blockade can be generated by mixing the photonic and phononic modes which do not exhibit blockade.


title: Image Synthesis with Disentangled Attributes for Chest X-Ray Nodule Augmentation and Detection

abstract: Lung nodule detection in chest X-ray (CXR) images is common to early screening of lung cancers. Deep-learning-based Computer-Assisted Diagnosis (CAD) systems can support radiologists for nodule screening in CXR. However, it requires large-scale and diverse medical data with high-quality annotations to train such robust and accurate CADs. To alleviate the limited availability of such datasets, lung nodule synthesis methods are proposed for the sake of data augmentation. Nevertheless, previous methods lack the ability to generate nodules that are realistic with the size attribute desired by the detector. To address this issue, we introduce a novel lung nodule synthesis framework in this paper, which decomposes nodule attributes into three main aspects including shape, size, and texture, respectively. A GAN-based Shape Generator firstly models nodule shapes by generating diverse shape masks. The following Size Modulation then enables quantitative control on the diameters of the generated nodule shapes in pixel-level granularity. A coarse-to-fine gated convolutional Texture Generator finally synthesizes visually plausible nodule textures conditioned on the modulated shape masks. Moreover, we propose to synthesize nodule CXR images by controlling the disentangled nodule attributes for data augmentation, in order to better compensate for the nodules that are easily missed in the detection task. Our experiments demonstrate the enhanced image quality, diversity, and controllability of the proposed lung nodule synthesis framework. We also validate the effectiveness of our data augmentation on greatly improving nodule detection performance.


title: Neural Greedy Pursuit for Feature Selection

abstract: We propose a greedy algorithm to select $N$ important features among $P$ input features for a non-linear prediction problem. The features are selected one by one sequentially, in an iterative loss minimization procedure. We use neural networks as predictors in the algorithm to compute the loss and hence, we refer to our method as neural greedy pursuit (NGP). NGP is efficient in selecting $N$ features when $N \ll P$, and it provides a notion of feature importance in a descending order following the sequential selection procedure. We experimentally show that NGP provides better performance than several feature selection methods such as DeepLIFT and Drop-one-out loss. In addition, we experimentally show a phase transition behavior in which perfect selection of all $N$ features without false positives is possible when the training data size exceeds a threshold.


title: A Near-Linear Time Sampler for the Ising Model with External Field

abstract: We give a near-linear time sampler for the Gibbs distribution of the ferromagnetic Ising models with edge activities $\boldsymbol{\beta} > 1$ and external fields $\boldsymbol{\lambda}<1$ (or symmetrically, $\boldsymbol{\lambda}>1$) on general graphs with bounded or unbounded maximum degree. Our algorithm is based on the field dynamics given in [CFYZ21]. We prove the correctness and efficiency of our algorithm by establishing spectral independence of distribution of the random cluster model and the rapid mixing of Glauber dynamics on the random cluster model in a low-temperature regime, which may be of independent interest.


title: Electronic materials with nanoscale curved geometries

abstract: Research into electronic nanomaterials has recently seen a growing focus into the synthesis of structures with unconventional curved geometries including bent wires in planar systems and three-dimensional architectures obtained by rolling up nanomembranes. The inclusion of these geometries has led to the prediction and observation of a series of novel effects that either result from shape-driven modifications of the electronic motion or from an intrinsic change of electronic and magnetic properties due to peculiar confinement effects. Moreover, local strains often generated by curvature also trigger the appearance of new phenomena due to the essential role played by electromechanical coupling in solids. Here we review the recent developments in the discovery of these shape-, confinement- and strain-induced curvature effects at the nanoscale, and discuss their potential use in electronic and spintronic devices.


title: Gravitational potential energy of a multi-component galactic disk

abstract: We calculate ab initio the gravitational potential energy per unit area for a gravitationally coupled multi-component galactic disk of stars and gas, which is given as the integration over vertical density distribution, vertical gravitational force, and vertical distance. This is based on the method proposed by Camm for a single-component disk, which we extend here for a multi-component disk by deriving the expression of the energy explicitly at any galactocentric radius R. For a self-consistent distribution, the density and force are obtained by jointly solving the equation of vertical hydrostatic equilibrium and the Poisson equation. Substituting the numerical values for the density distribution and force obtained for the coupled system, in the derived expression of the energy, we find that the energy of each component remains unchanged compared to the energy for the corresponding single-component case. We explain this surprising result by simplifying the above expression for the energy of a component analytically, which turns out to be equal to the surface density times the squared vertical velocity dispersion of the component. However, the energy required to raise a unit test mass to a certain height z from the mid-plane is higher in the coupled case. The system is therefore more tightly bound closer to the mid-plane, and hence it is harder to disturb it due to an external tidal encounter.


title: Mimetic Models: Ethical Implications of AI that Acts Like You

abstract: An emerging theme in artificial intelligence research is the creation of models to simulate the decisions and behavior of specific people, in domains including game-playing, text generation, and artistic expression. These models go beyond earlier approaches in the way they are tailored to individuals, and the way they are designed for interaction rather than simply the reproduction of fixed, pre-computed behaviors. We refer to these as mimetic models, and in this paper we develop a framework for characterizing the ethical and social issues raised by their growing availability. Our framework includes a number of distinct scenarios for the use of such models, and considers the impacts on a range of different participants, including the target being modeled, the operator who deploys the model, and the entities that interact with it.


title: Incentive Design for Large Congestion Games: Publicness-Specific Bayes Correlated Wardrop Equilibrium (Extended Abstract)

abstract: The travel costs of the players (travelers) in anonymous congestion games depend on their choices of routes and also on the states of the transportation network such as incidents, weather, and road work. In this extended abstract, we consider an incomplete-information environment in which the realizations of the states are unobserved by the travelers. We study how a planner can incentivize the travelers to behave in her favor by strategically designing what and how the travelers get informed about the realizations of the states.


title: Parametric models and information geometry on W*-algebras

abstract: We introduce the notion of smooth parametric model of normal positive linear functionals on possibly infinite-dimensional W*-algebras generalizing the notions of parametric models used in classical and quantum information geometry. We then use the Jordan product naturally available in this context in order to define a Riemannian metric tensor on parametric models satsfying suitable regularity conditions. This Riemannian metric tensor reduces to the Fisher-Rao metric tensor, or to the Fubini-Study metric tensor, or to the Bures-Helstrom metric tensor when suitable choices for the W*-algebra and the models are made.


title: Composition Theorems for Interactive Differential Privacy

abstract: An interactive mechanism is an algorithm that stores a data set and answers adaptively chosen queries to it. The mechanism is called differentially private, if any adversary cannot distinguish whether a specific individual is in the data set by interacting with the mechanism. We study composition properties of differential privacy in concurrent compositions. In this setting, an adversary interacts with k interactive mechanisms in parallel and can interleave its queries to the mechanisms arbitrarily. Previously, Vadhan and Wang [2021] proved an optimal concurrent composition theorem for pure-differential privacy. We significantly generalize and extend their results. Namely, we prove optimal parallel composition properties for several major notions of differential privacy in the literature, including approximate DP, R\'enyi DP, and zero-concentrated DP. Our results demonstrate that the adversary gains no advantage by interleaving its queries to independently running mechanisms. Hence, interactivity is a feature that differential privacy grants us for free. Concurrently and independently of our work, Vadhan and Zhang [2022] proved an optimal concurrent composition theorem for f-DP [Dong et al., 2022], which implies our result for the approximate DP case.


title: Positivity-Preserving Well-Balanced Central Discontinuous Galerkin Schemes for the Euler Equations under Gravitational Fields

abstract: This paper designs and analyzes positivity-preserving well-balanced (WB) central discontinuous Galerkin (CDG) schemes for the Euler equations with gravity. A distinctive feature of these schemes is that they not only are WB for a general known stationary hydrostatic solution, but also can preserve the positivity of the fluid density and pressure. The standard CDG method does not possess this feature, while directly applying some existing WB techniques to the CDG framework may not accommodate the positivity and keep other important properties at the same time. In order to obtain the WB and positivity-preserving properties simultaneously while also maintaining the conservativeness and stability of the schemes, a novel spatial discretization is devised in the CDG framework based on suitable modifications to the numerical dissipation term and the source term approximation. The modifications are based on a crucial projection operator for the stationary hydrostatic solution, which is proposed for the first time in this work. This novel projection has the same order of accuracy as the standard $L^2$-projection, can be explicitly calculated, and is easy to implement without solving any optimization problems. More importantly, it ensures that the projected stationary solution has the same cell averages on both the primal and dual meshes, which is a key to achieve the desired properties of our schemes. Based on some convex decomposition techniques, rigorous positivity-preserving analyses for the resulting WB CDG schemes are carried out. Several one- and two-dimensional numerical examples are performed to illustrate the desired properties of these schemes, including the high-order accuracy, the WB property, the robustness for simulations involving the low pressure or density, high resolution for the discontinuous solutions and the small perturbations around the equilibrium state.


title: RCLane: Relay Chain Prediction for Lane Detection

abstract: Lane detection is an important component of many real-world autonomous systems. Despite a wide variety of lane detection approaches have been proposed, reporting steady benchmark improvements over time, lane detection remains a largely unsolved problem. This is because most of the existing lane detection methods either treat the lane detection as a dense prediction or a detection task, few of them consider the unique topologies (Y-shape, Fork-shape, nearly horizontal lane) of the lane markers, which leads to sub-optimal solution. In this paper, we present a new method for lane detection based on relay chain prediction. Specifically, our model predicts a segmentation map to classify the foreground and background region. For each pixel point in the foreground region, we go through the forward branch and backward branch to recover the whole lane. Each branch decodes a transfer map and a distance map to produce the direction moving to the next point, and how many steps to progressively predict a relay station (next point). As such, our model is able to capture the keypoints along the lanes. Despite its simplicity, our strategy allows us to establish new state-of-the-art on four major benchmarks including TuSimple, CULane, CurveLanes and LLAMAS.


title: Derivative of a hypergraph as a tool for linguistic pattern analysis

abstract: The search for linguistic patterns, stylometry and forensic linguistics have in the theory of complex networks, their structures and associated mathematical tools, allies with which to model and analyze texts. In this paper we present a new model supported by several mathematical structures such as the hypergraphs or the concept of derivative graph to introduce a new methodology able to analyze the mesoscopic relationships between sentences, paragraphs, chapters and texts, focusing not only in a quantitative index but also in a new mathematical structure that will be of singular help to both: detecting the style of an author and determining the language level of a text. In addition, these new mathematical structures may be useful to detect similarity and dissimilarity in texts and, eventually, even plagiarism.


title: Properties of the gradient squared of the discrete Gaussian free field

abstract: In this paper we study the properties of the centered (norm of the) gradient squared of the discrete Gaussian free field in $U_{\epsilon}=U/\epsilon\cap \mathbb{Z}^d$, $U\subset \mathbb{R}^d$ and $d\geq 2$. The covariance structure of the field is a function of the transfer current matrix and this relates the model to a class of systems (e.g. height-one field of the Abelian sandpile model or pattern fields in dimer models) that have a Gaussian limit due to the rapid decay of the transfer current. Indeed, we prove that the properly rescaled field converges to white noise in an appropriate local Besov-H\"older space. Moreover, under a different rescaling, we determine the $k$-point correlation function and cumulants on $U_{\epsilon}$ and in the continuum limit as $\epsilon\to 0$. This result is related to the analogue limit for the height-one field of the Abelian sandpile (\citet{durre}), with the same conformally covariant property in $d=2$.


title: Residence time distributions for in-line chaotic mixers

abstract: We investigate the distributions of residence time for in-line chaotic mixers; in particular, we consider the Kenics, the F-mixer and the Multi-level laminating mixer, and also a synthetic model that mimics their behavior and allows exact mathematical calculations. We show that whatever the number of elements of mixer involved, the distribution possesses a $t^{-3}$ tail, so that its shape is always far from Gaussian. This $t^{-3}$ tail also invalidates the use of second-order moment and variance. As a measure for the width of the distribution, we consider the mean absolute deviation and show that, unlike the standard deviation, it converges in the limit of large sample size. Finally, we analyze the performances of the different in-line mixers from the residence-time point of view when varying the number of elements and the shape of the cross-section.


title: A General Wasserstein Framework for Data-driven Distributionally Robust Optimization: Tractability and Applications

abstract: Data-driven distributionally robust optimization is a recently emerging paradigm aimed at finding a solution that is driven by sample data but is protected against sampling errors. An increasingly popular approach, known as Wasserstein distributionally robust optimization (DRO), achieves this by applying the Wasserstein metric to construct a ball centred at the empirical distribution and finding a solution that performs well against the most adversarial distribution from the ball. In this paper, we present a general framework for studying different choices of a Wasserstein metric and point out the limitation of the existing choices. In particular, while choosing a Wasserstein metric of a higher order is desirable from a data-driven perspective, given its less conservative nature, such a choice comes with a high price from a robustness perspective - it is no longer applicable to many heavy-tailed distributions of practical concern. We show that this seemingly inevitable trade-off can be resolved by our framework, where a new class of Wasserstein metrics, called coherent Wasserstein metrics, is introduced. Like Wasserstein DRO, distributionally robust optimization using the coherent Wasserstein metrics, termed generalized Wasserstein distributionally robust optimization (GW-DRO), has all the desirable performance guarantees: finite-sample guarantee, asymptotic consistency, and computational tractability. The worst-case expectation problem in GW-DRO is in general a nonconvex optimization problem, yet we provide new analysis to prove its tractability without relying on the common duality scheme. Our framework, as shown in this paper, offers a fruitful opportunity to design novel Wasserstein DRO models that can be applied in various contexts such as operations management, finance, and machine learning.


title: Exploring shock-capturing schemes for Particles on Demand simulation of compressible flows

abstract: In this exploratory study, we apply shock-capturing schemes within the framework of the Particles on Demand kinetic model to simulate compressible flows with mild and strong shock waves and discontinuities. The model is based on the semi-Lagrangian method where the information propagates along the characteristics while a set of shock-capturing concepts such as the total variation diminishing and weighted essentially non-oscillatory schemes are employed to capture the discontinuities and the shock-waves. The results show that the reconstruction schemes are able to remove the oscillations at the location of the shock waves and together with the Galilean invariance nature of the Particles on Demand model, stable simulations of mild to extreme compressible benchmarks can be carried out. Moreover, the essential numerical properties of the reconstruction schemes such as their spectral analysis and order of accuracy are discussed.


title: Bayesian Generational Population-Based Training

abstract: Reinforcement learning (RL) offers the potential for training generally capable agents that can interact autonomously in the real world. However, one key limitation is the brittleness of RL algorithms to core hyperparameters and network architecture choice. Furthermore, non-stationarities such as evolving training data and increased agent complexity mean that different hyperparameters and architectures may be optimal at different points of training. This motivates AutoRL, a class of methods seeking to automate these design choices. One prominent class of AutoRL methods is Population-Based Training (PBT), which have led to impressive performance in several large scale settings. In this paper, we introduce two new innovations in PBT-style methods. First, we employ trust-region based Bayesian Optimization, enabling full coverage of the high-dimensional mixed hyperparameter search space. Second, we show that using a generational approach, we can also learn both architectures and hyperparameters jointly on-the-fly in a single training run. Leveraging the new highly parallelizable Brax physics engine, we show that these innovations lead to large performance gains, significantly outperforming the tuned baseline while learning entire configurations on the fly. Code is available at https://github.com/xingchenwan/bgpbt.


title: A Level Set Kalman Filter Approach to Estimate the Circadian Phase and its Uncertainty from Wearable Data

abstract: The circadian clock is an internal timer that coordinates the daily rhythms of behavior and physiology, including sleep and hormone secretion. Accurately tracking the state of the circadian clock, or circadian phase, holds immense potential for precision medicine. Wearable devices present an opportunity to estimate the circadian phase in the real world, as they can non-invasively monitor various physiological outputs influenced by the circadian clock. However, accurately estimating circadian phase from wearable data remains challenging, primarily due to the lack of methods that integrate minute-by-minute wearable data with prior knowledge of the circadian phase. To address this issue, we propose a framework that integrates multi-time scale physiological data to estimate the circadian phase, along with an efficient implementation algorithm based on Bayesian inference and a new state space estimation method called the level set Kalman filter. Our numerical experiments indicate that our approach outperforms previous methods for circadian phase estimation consistently. Furthermore, our method enables us to examine the contribution of noise from different sources to the estimation, which was not feasible with prior methods. We found that internal noise unrelated to external stimuli is a crucial factor in determining estimation results. Lastly, we developed a user-friendly computational package and applied it to real-world data to demonstrate the potential value of our approach. Our results provide a foundation for systematically understanding the real-world dynamics of the circadian clock.


title: High harmonic generation from a subwavelength dielectric resonator

abstract: Higher-order optical harmonics entered the realm of nanostructured solids being observed recently in optical gratings and metasurfaces with a subwavelength thickness. Structuring materials at the subwavelength scale allows for resonant enhancing of the efficiency of nonlinear processes and reducing the size of high-harmonic sources. Here we report the observation of up to a seventh harmonic generated from a single subwavelength resonator made of AlGaAs material. This process is enabled by careful engineering of the resonator geometry for supporting optical modes associated with a quasi-bound state in the continuum in the mid-infrared spectral range at around {\lambda}=3.7 {\mu}m pump wavelength. The resonator volume measures ~ 0.1 {\lambda}^3. The resonant modes are excited with an azimuthally polarized tightly focused beam. We evaluate the contributions of perturbative and non-perturbative nonlinearities to the harmonic generation process. Our work proves the possibility to miniaturize solid-state sources of high harmonics to the subwavelength volumes.


title: Bounding generalization error with input compression: An empirical study with infinite-width networks

abstract: Estimating the Generalization Error (GE) of Deep Neural Networks (DNNs) is an important task that often relies on availability of held-out data. The ability to better predict GE based on a single training set may yield overarching DNN design principles to reduce a reliance on trial-and-error, along with other performance assessment advantages. In search of a quantity relevant to GE, we investigate the Mutual Information (MI) between the input and final layer representations, using the infinite-width DNN limit to bound MI. An existing input compression-based GE bound is used to link MI and GE. To the best of our knowledge, this represents the first empirical study of this bound. In our attempt to empirically falsify the theoretical bound, we find that it is often tight for best-performing models. Furthermore, it detects randomization of training labels in many cases, reflects test-time perturbation robustness, and works well given only few training samples. These results are promising given that input compression is broadly applicable where MI can be estimated with confidence.


title: Tur\'an Inequalities for Infinite Product Generating Functions

abstract: In the $1970$s, Nicolas proved that the partition function $p(n)$ is log-concave for $ n > 25$. In \cite{HNT21}, a precise conjecture on the log-concavity for the plane partition function $\func{pp}(n)$ for $n >11$ was stated. This was recently proven by Ono, Pujahari, and Rolen. In this paper, we provide a general picture. We associate to double sequences $\{g_d(n)\}_{d,n}$ with $g_d(1)=1$ and $$0 \leq g_{d}\left( n\right) - n^{d}\leq g_{1}\left( n\right) \left( n-1\right) ^{d-1}$$ polynomials $\{P_n^{g_d}(x)\}_{d,n}$ given by \begin{equation*} \sum_{n=0}^{\infty} P_n^{g_d}(x) \, q^n := \func{exp}\left( x \sum_{n=1}^{\infty} g_d(n) \frac{q^n}{n} \right) =\prod_{n=1}^{\infty} \left( 1 - q^n \right)^{-x f_d(n)}. \end{equation*} We recover $ p(n)= P_n^{\sigma_1}(1)$ and $\func{pp}\left( n\right) = P_n^{\sigma_2}(1)$, where $\sigma_d (n):= \sum_{\ell \mid n} \ell^d$ and $f_d(n)= n^{d-1}$. Let $n \geq 6$. Then the sequence $\{P_n^{\sigma_d}(1)\}_d$ is log-concave for almost all $d$ if and only if $n$ is divisible by $3$. Let $\func{id}(n)=n$. Then $P_n^{\func{id}}(x) = \frac{x}{n} L_{n-1}^{(1)}(-x)$, where $L_{n}^{\left( \alpha \right) }\left( x\right) $ denotes the $\alpha$-associated Laguerre polynomial. In this paper, we invest in Tur\'an inequalities \begin{equation*} \Delta_{n}^{g_d}(x) := \left( P_n^{g_d}(x) \right)^2 - P_{n-1}^{g_d}(x) \, P_{n+1}^{g_d}(x) \geq 0. \end{equation*} Let $n \geq 6$ and $0 \leq x < 2 - \frac{12}{n+4}$. Then $n$ is divisible by $3$ if and only if $\Delta_{n}^{g_d}(x) \geq 0$ for almost all $d$. Let $n \geq 6$ and $n \not\equiv 2 \pmod{3}$. Then the condition on $x$ can be reduced to $x \geq 0$. We determine explicit bounds. As an analogue to Nicolas' result, we have for $g_1= \func{id}$ that $\Delta_{n}^{\func{id}}(x) \geq 0$ for all $x \geq 0 $ and all $n$.


title: An Elementary Proof of a Theorem of Hardy and Ramanujan

abstract: Let $Q(n)$ denote the number of integers $1 \leq q \leq n$ whose prime factorization $q= \prod^{t}_{i=1}p^{a_i}_i$ satisfies $a_1\geq a_2\geq \ldots \geq a_t$. Hardy and Ramanujan proved that $$ \log Q(n) \sim \frac{2\pi}{\sqrt{3}} \sqrt{\frac{\log(n)}{\log\log(n)}}\;. $$ Before proving the above precise asymptotic formula, they studied in great detail what can be obtained concerning $Q(n)$ using purely elementary methods, and were only able to obtain much cruder lower and upper bounds using such methods. In this paper we show that it is in fact possible to obtain a purely elementary (and much shorter) proof of the Hardy--Ramanujan Theorem. Towards this goal, we first give a simple combinatorial argument, showing that $Q(n)$ satisfies a (pseudo) recurrence relation. This enables us to replace almost all the hard analytic part of the original proof with a short inductive argument.


title: Influence of grain growth on CO2 ice spectroscopic profiles : Modelling for dense cores and disks

abstract: Interstellar dust grain growth in dense clouds and protoplanetary disks, even moderate, affects the observed interstellar ice profiles as soon as a significant fraction of dust grains is in the size range close to the wave vector at the considered wavelength. The continuum baseline correction made prior to analysing ice profiles influences the subsequent analysis and hence the estimated ice composition, typically obtained by band fitting using thin film ice mixture spectra. We model the effect of grain growth on ice mantle spectroscopic profiles, focusing on CO2 to see how it can affect interstellar ice mantle spectral analysis and interpretation. Using the Discrete Dipole Approximation for Scattering and Absorption of Light, the mass absorption coefficients of distributions of grains composed of ellipsoidal silicate cores with water and carbon dioxide ice mantles are calculated. A few other ice mantle compositions are also calculated. We explore the size distribution evolution from dense clouds to simulate the first steps of grain growth up to three microns in size. The results are injected into RADMC-3D full scattering radiative transfer models of spherical clouds and protoplanetary disk templates to retrieve observable spectral energy distributions. We focus on calculated profile of the CO2 antisymmetric stretching mode ice band profile at 4.27 microns, a meaningful indicator of grain growth. The observed profiles toward dense cores with the Infrared space observatory and Akari satellites already showed profiles possibly indicative of moderate grain growth.The observation of protoplanetary disks at high inclination with the JWST should present distorted profiles that will put constraints on the extent of dust growth. The more evolved the dust size distribution, the more the extraction of the ice mantle composition will require both understanding and taking into account grain growth.


title: Det6D: A Ground-Aware Full-Pose 3D Object Detector for Improving Terrain Robustness

abstract: Accurate 3D object detection with LiDAR is critical for autonomous driving. Existing research is all based on the flat-world assumption. However, the actual road can be complex with steep sections, which breaks the premise. Current methods suffer from performance degradation in this case due to difficulty correctly detecting objects on sloped terrain. In this work, we propose Det6D, the first full-degree-of-freedom 3D object detector without spatial and postural limitations, to improve terrain robustness. We choose the point-based framework by founding their capability of detecting objects in the entire spatial range. To predict full-degree poses, including pitch and roll, we design a ground-aware orientation branch that leverages the local ground constraints. Given the difficulty of long-tail non-flat scene data collection and 6D pose annotation, we present Slope-Aug, a data augmentation method for synthesizing non-flat terrain from existing datasets recorded in flat scenes. Experiments on various datasets demonstrate the effectiveness and robustness of our method in different terrains. We further conducted an extended experiment to explore how the network predicts the two extra poses. The proposed modules are plug-and-play for existing point-based frameworks. The code is available at https://github.com/HITSZ-NRSL/De6D.


title: SphereFed: Hyperspherical Federated Learning

abstract: Federated Learning aims at training a global model from multiple decentralized devices (i.e. clients) without exchanging their private local data. A key challenge is the handling of non-i.i.d. (independent identically distributed) data across multiple clients that may induce disparities of their local features. We introduce the Hyperspherical Federated Learning (SphereFed) framework to address the non-i.i.d. issue by constraining learned representations of data points to be on a unit hypersphere shared by clients. Specifically, all clients learn their local representations by minimizing the loss with respect to a fixed classifier whose weights span the unit hypersphere. After federated training in improving the global model, this classifier is further calibrated with a closed-form solution by minimizing a mean squared loss. We show that the calibration solution can be computed efficiently and distributedly without direct access of local data. Extensive experiments indicate that our SphereFed approach is able to improve the accuracy of multiple existing federated learning algorithms by a considerable margin (up to 6% on challenging datasets) with enhanced computation and communication efficiency across datasets and model architectures.


title: Batchelor, Saffman, and Kazantsev spectra in galactic small-scale dynamos

abstract: The magnetic fields in galaxy clusters and probably also in the interstellar medium are believed to be generated by a small-scale dynamo. Theoretically, during its kinematic stage, it is characterized by a Kazantsev spectrum, which peaks at the resistive scale. It is only slightly shallower than the Saffman spectrum that is expected for random and causally connected magnetic fields. Causally disconnected fields have the even steeper Batchelor spectrum. Here we show that all three spectra are present in the small-scale dynamo. During the kinematic stage, the Batchelor spectrum occurs on scales larger than the energy-carrying scale of the turbulence, and the Kazantsev spectrum on smaller scales within the inertial range of the turbulence -- even for a magnetic Prandtl number of unity. In the saturated state, the dynamo develops a Saffman spectrum on large scales, suggestive of the build-up of long-range correlations. At large magnetic Prandtl numbers, elongated structures are seen in synthetic synchrotron emission maps showing the parity-even E polarization. We also observe a significant excess in the E polarization over the parity-odd B polarization at subresistive scales, and a deficiency at larger scales. This finding is at odds with the observed excess in the Galactic microwave foreground emission, which is believed to be associated with larger scales. The E and B polarizations may be highly non-Gaussian and skewed in the kinematic regime of the dynamo. For dust emission, however, the polarized emission is always nearly Gaussian, and the excess in the E polarization is much weaker.


title: LOTUS: A (non-)LTE Optimization Tool for Uniform derivation of Stellar atmospheric parameters

abstract: Precise fundamental atmospheric stellar parameters and abundance determination of individual elements in stars are important for all stellar population studies. Non-Local Thermodynamic Equilibrium (Non-LTE; hereafter NLTE) models are often important for such high precision, however, can be computationally complex and expensive, which renders the models less utilized in spectroscopic analyses. To alleviate the computational burden of such models, we developed a robust 1D, LTE and NLTE fundamental atmospheric stellar parameter derivation tool, $\texttt{LOTUS}$, to determine the effective temperature $T_{\mathrm{eff}}$, surface gravity $\log g$, metallicity $\mbox{[Fe/H]}$ and microturbulent velocity $v_{\mathrm{mic}}$ for FGK type stars, from equivalent width (EW) measurements of Fe I and Fe II lines. We utilize a generalized curve of growth method to take into account the EW dependencies of each Fe I and Fe II line on the corresponding atmospheric stellar parameters. A global differential evolution optimization algorithm is then used to derive the optimized fundamental parameters. Additionally, $\texttt{LOTUS}$ can determine precise uncertainties for each stellar parameter using a Markov Chain Monte Carlo (MCMC) algorithm. We test and apply $\texttt{LOTUS}$ on a sample of benchmark stars, as well as stars with available asteroseismic surface gravities from the K2 survey, and metal-poor stars from $R$-process Alliance (RPA) survey. We find very good agreement between our NLTE-derived parameters in $\texttt{LOTUS}$ to non-spectroscopic values within $T_{\mathrm{eff}}=\pm 30$ K and $\log g=\pm 0.10$ dex for benchmark stars. We provide open access of our code, as well as of the interpolated pre-computed NLTE EW grids available on Github, and documentation with working examples on Readthedocs.


title: New high-precision strong lensing modeling of Abell 2744. Preparing for JWST observations

abstract: We present a new strong lensing (SL) model of the Hubble Frontier Fields galaxy cluster Abell 2744, at z=0.3072, by exploiting archival Hubble Space Telescope (HST) multi-band imaging and Multi Unit Spectroscopic Explorer (MUSE) follow-up spectroscopy. The lens model considers 90 spectroscopically confirmed multiple images (from 30 background sources), which represents the largest secure sample for this cluster field prior to the recently acquired James Webb Space Telescope observations. The inclusion of the sub-structures within several extended sources as model constraints allows us to accurately characterize the inner total mass distribution of the cluster and the position of the cluster critical lines. We include the lensing contribution of 225 cluster members, 202 of which are spectroscopically confirmed. We also measure the internal velocity dispersion of 85 cluster galaxies to independently estimate the role of the subhalo mass component in the lens model. We investigate the effect of the cluster environment on the total mass reconstruction of the cluster core with two different mass parameterizations. We consider the mass contribution from three external clumps, either based on previous weak-lensing studies, or extended HST imaging of luminous members around the cluster core. In the latter case, the observed positions of the multiple images are better reproduced, with a remarkable accuracy of 0.37", a factor of $\sim2$ smaller than previous lens models. We develop and make publicly available a Strong Lensing Online Tool (SLOT) to exploit the predictive power and the full statistical information of this and future models, through a simple graphical interface. We plan to apply our high-precision SL model to the first analysis of the GLASS-JWST-ERS program, specifically to measure the intrinsic physical properties of high-$z$ galaxies from robust magnification maps.


title: Multiple solutions and profile description for a nonlinear Schr\"odinger-Bopp-Podolsky-Proca system on a manifold

abstract: We prove a multiplicity result for \begin{equation*} \begin{cases} -\varepsilon^{2}\Delta_g u+\omega u+q^{2}\phi u=|u|^{p-2}u\\[1mm] -\Delta_g \phi +a^{2}\Delta_g^{2} \phi + m^2 \phi =4\pi u^{2} \end{cases} \text{ in }M, \end{equation*} where $(M,g)$ is a smooth and compact $3$-dimensional Riemannian manifold without boundary, $p\in(4,6)$, $a,m,q\neq 0$, $\varepsilon>0$ small enough. The proof of this result relies on Lusternik-Schnirellman category. We also provide a profile description for low energy solutions.


title: Unrolled algorithms for group synchronization

abstract: The group synchronization problem involves estimating a collection of group elements from noisy measurements of their pairwise ratios. This task is a key component in many computational problems, including the molecular reconstruction problem in single-particle cryo-electron microscopy (cryo-EM). The standard methods to estimate the group elements are based on iteratively applying linear and non-linear operators, and are not necessarily optimal. Motivated by the structural similarity to deep neural networks, we adopt the concept of algorithm unrolling, where training data is used to optimize the algorithm. We design unrolled algorithms for several group synchronization instances, including synchronization over the group of 3-D rotations: the synchronization problem in cryo-EM. We also apply a similar approach to the multi-reference alignment problem. We show by numerical experiments that the unrolling strategy outperforms existing synchronization algorithms in a wide variety of scenarios.


title: Fusion cross section for the reaction $^{18}$O+$^{12}$C at 16.7 MeV/nucleon

abstract: The absolute cross section for the complete fusion and the total reaction cross section as a function of the charge of the heavy residue have been indirectly estimated for the reaction $^{18}$O+$^{12}$C at 16.7 MeV/nucleon. Data have been collected with the GARFIELD+RCo setup at LNL (Italy). The mbarn/count normalization factor has been extracted using the statistical model HF$\ell$ (Hauser-Feshbach light).


title: Feasibility study for precisely measuring the EIC ${}^3$He beam polarization with the Polarized Atomic Hydrogen Gas Jet Target polarimeter at RHIC

abstract: The Polarized Atomic Hydrogen Gas Jet Target polarimeter (HJET) is used to measure the absolute proton beam polarization, $\sigma_P^\text{syst}/P\!\lesssim\!0.5\%$, at the Relativistic Heavy Ion Collider. Here I consider the possibility of employing HJET to measure the ${}^3\text{He}$ ($h$) beam polarization at the Electron-Ion Collider (EIC). The dominant contribution to the ratio of the $h^\uparrow{p}$ and $p^\uparrow{h}$ analyzing powers, which is needed for such measurements, can be easily calculated using well-known values of the proton and helion magnetic moments, but some corrections should be applied to achieve the required accuracy. It was found that corrections due to absorption and ${}^3\text{He}$ breakup effectively cancel in the ratio and a correction due to hadronic spin-flip amplitudes can be derived from the proton beam measurements. As a result, the anticipated systematic uncertainty in the measured ${}^3\text{He}$ beam polarization can satisfy the EIC requirement $\sigma_P^\text{syst}/P\!\lesssim\!1\%$.


title: A Modest View of the Black Hole Information Paradox

abstract: Thirty years ago, John Preskill concluded "that the information loss paradox may well presage a revolution in fundamental physics" and mused that "Conceivably, the puzzle of black hole evaporation portends a scientific revolution as sweeping as that that led to the formulation of quantum theory in the early 20th century." Many still agree with this assessment. On the other hand, it seems to me the "paradox" has little to do with the physical world but rather, at best, simply points out the possible inconsistency of two, already disparate, theories (mathematical models) of nature, general relativity and quantum mechanics, with virtually no conceivable observational consequences. The information paradox hinges on the concepts of a pure quantum state, the unitarity of quantum mechanics, and Hawking's semi-classical calculation of black hole evaporation. I used the qualifier "at best" above because, for me, the concept of a quantum state is far more restrictive than required by the paradox while unitarity is not a property of nature but rather of a mathematical model and is already violated by the process of making a measurement. Furthermore, the semi-classical calculation of Hawking is surely of limited applicability.


title: Machine Learning Assisted Resistive Force Theory for Helical Structures at Low Reynolds Number

abstract: The hydrodynamic forces on a slender rod in a fluid medium at low Reynolds number can be modeled using resistive force theories (RFTs) or slender body theories (SBTs). The former represent the forces by local drag coefficients and are computationally cheap; however, they are physically inaccurate when long-range hydrodynamic interaction is involved. The later are physically accurate but require solving integral equations and, therefore, are computationally expensive. This paper investigates RFTs in comparison with state-of-the art SBT methods. During the process, a neural network-based hydrodynamic model that -- similar to RFTs -- relies on local drag coefficients for computational efficiency was developed. However, the network is trained using data from an SBT (regularized stokeslet segments method). The $R^2$ value of the trained coefficients were $\sim 0.99$ with mean absolute error of $1.6\times10^{-2}$. The machine learning resistive force theory (MLRFT) accounts for local hydrodynamic forces distribution, the dependence on rotational and translational speeds and directions, and geometric parameters of the slender object. We show that, when classical RFT fails to accurately predict the forces, torques, and drags on slender rods under low Reynolds number flows, MLRFT exhibits good agreement with physically accurate SBT simulations. In terms of computational speed, MLRFT forgoes the need of solving an inverse problem and, therefore, requires negligible computation time in comparison with SBT. MLRFT presents a computationally inexpensive hydrodynamic model for flagellar propulsion can be used in the design and optimization of biomimetic flagellated robots and analysis of bacterial locomotion.


title: Exploring the configurational space of amorphous graphene with machine-learned atomic energies

abstract: Two-dimensionally extended amorphous carbon ("amorphous graphene") is a prototype system for disorder in 2D, showing a rich and complex configurational space that is yet to be fully understood. Here we explore the nature of amorphous graphene with an atomistic machine-learning (ML) model. We create structural models by introducing defects into ordered graphene through Monte-Carlo bond switching, defining acceptance criteria using the machine-learned local, atomic energies associated with a defect, as well as the nearest-neighbor (NN) environments. We find that physically meaningful structural models arise from ML atomic energies in this way, ranging from continuous random networks to paracrystalline structures. Our results show that ML atomic energies can be used to guide Monte-Carlo structural searches in principle, and that their predictions of local stability can be linked to short- and medium-range order in amorphous graphene. We expect that the former point will be relevant more generally to the study of amorphous materials, and that the latter has wider implications for the interpretation of ML potential models.


title: Momentum Injection via Dielectric Barrier Discharge Actuators in Low-Speed External Flow

abstract: Dielectric barrier discharge (DBD) plasma actuators can generate a wall jet without moving parts through interaction between ionized and neutral molecules in an electric field. The coupling between electro-hydrodynamic, turbulence, and viscous effects in the flow boundary layer remains unclear and deserves careful investigation. We present an experimental investigation of momentum injection by DBD actuators in a U_external = 5 m/s and U_external = 11 m/s co-flow and counter-flow configuration over a range of VAC = 12 kV - 19.5 kV peak-to-peak at a frequency of 2 kHz. In the co-flow configuration, the DBD actuator adds momentum to the boundary layer, similar to an electrohydrodynamic (EHD) jet in quiescent conditions. In the counter-flow configuration, flow separation is observed at free stream velocity U_external = 5 m/s. The momentum displacement in the counter-flow configuration is ~ 6x greater than EHD jet momentum in a quiescent environment. Both co-flow and counter-flow momentum injections show diminishing effects with increased external flow speed. This work highlights that the resulting flow pattern is not a simple superposition of the EHD jet and the free stream but is determined by a balance between the inertial, viscous and Coulombic forces of the EHD and the external flow. The velocity profiles and momentum characteristics can be used to validate numerical models and inform the design of DBD actuators for active flow control.


title: Geometric Features Informed Multi-person Human-object Interaction Recognition in Videos

abstract: Human-Object Interaction (HOI) recognition in videos is important for analyzing human activity. Most existing work focusing on visual features usually suffer from occlusion in the real-world scenarios. Such a problem will be further complicated when multiple people and objects are involved in HOIs. Consider that geometric features such as human pose and object position provide meaningful information to understand HOIs, we argue to combine the benefits of both visual and geometric features in HOI recognition, and propose a novel Two-level Geometric feature-informed Graph Convolutional Network (2G-GCN). The geometric-level graph models the interdependency between geometric features of humans and objects, while the fusion-level graph further fuses them with visual features of humans and objects. To demonstrate the novelty and effectiveness of our method in challenging scenarios, we propose a new multi-person HOI dataset (MPHOI-72). Extensive experiments on MPHOI-72 (multi-person HOI), CAD-120 (single-human HOI) and Bimanual Actions (two-hand HOI) datasets demonstrate our superior performance compared to state-of-the-arts.


title: Interstellar Complex Organic Molecules in SiO-traced massive outflows

abstract: The interstellar medium contains dust and gas, in which at high densities and cold conditions molecules can proliferate. Interstellar Complex Organic Molecules (iCOMs) are C-bearing species that contain at least six atoms. As they are detected in young stellar objects, iCOMs are expected to inhabit the early stages of the star formation evolution. In this study, we try to determine which iCOMs are present in the outflow component of massive protostars. To do this, we analyzed the morphological extension of blue- and red-shifted iCOMs emission in a sample of eleven massive protostars employing mapping observations at one mm within a ~1 GHz bandwidth for both IRAM-30m and APEX telescopes. We modeled the iCOMs emission of the central pointing spectra of our objects using the XCLASS LTE radiative transfer code. We detected the presence of several iCOMs such as: CH3OH, 13CH3OH, CH3OCHO, C2H5C15N, and c-C3H2)CH2. In G034.41+0.24, G327.29-0.58, G328.81+0.63, G333.13-0.43, G340.97-1.02, G351.45+0.66 and G351.77-0.54, the iCOMs lines show a faint broad line profile. Due to the offset peak positions of the blue- and red-shifted emission, covering from ~0.1 to 0.5 pc, these wings are possibly related to movements external to the compact core, such as large-scale low-velocity outflows. We have also established a correlation between the parent iCOM molecule CH3OH and the shock tracer SiO, reinforcing the hypothesis that shock environments provide the conditions to boost the formation of iCOMs via gas-phase reactions.


title: Common graphs with arbitrary connectivity and chromatic number

abstract: A graph $H$ is common if the number of monochromatic copies of $H$ in a 2-edge-colouring of the complete graph $K_n$ is asymptotically minimised by the random colouring. We prove that, given $k,r>0$, there exists a $k$-connected common graph with chromatic number at least $r$. The result is built upon the recent breakthrough of Kr\'a\v{l}, Volec, and Wei who obtained common graphs with arbitrarily large chromatic number and answers a question of theirs.


title: Panic! At the Disks: First Rest-frame Optical Observations of Galaxy Structure at $z > 3$ with JWST in the SMACS 0723 Field

abstract: We present early results regarding the morphological and structural properties of galaxies seen with the James Webb Space Telescope at $z > 3$ in the Early Release Observations of SMACS 0723, a galaxy cluster at $z=0.39$. We investigate, for the first time, the optical morphologies of a significant number of $z > 3$ galaxies with accurate photometric redshifts in this field to determine the form of galaxy structure in the relatively early universe. We use visual morphologies and \textsc{Morfometryka} measures to perform quantitative morphology measurements, both parametric with light profile fitting (S\'ersic indices) and non-parametric (CAS values). Using these, we measure the relative fraction of disk, spheroidal, and peculiar galaxies at $3 < z < 8$. We discover the surprising result that at $z > 1.5$ disk galaxies dominate the overall fraction of morphologies, with a factor of $\sim 10$ relative higher number of disk galaxies than seen by the Hubble Space Telescope at these redshifts. Our visual morphological estimates of galaxies align closely with their locations in CAS parameter space and their S\'ersic indices.


title: Prior-Independent Auctions for Heterogeneous Bidders

abstract: We study the design of prior-independent auctions in a setting with heterogeneous bidders. In particular, we consider the setting of selling to $n$ bidders whose values are drawn from $n$ independent but not necessarily identical distributions. We work in the robust auction design regime, where we assume the seller has no knowledge of the bidders' value distributions and must design a mechanism that is prior-independent. While there have been many strong results on prior-independent auction design in the i.i.d. setting, not much is known for the heterogeneous setting, even though the latter is of significant practical importance. Unfortunately, no prior-independent mechanism can hope to always guarantee any approximation to Myerson's revenue in the heterogeneous setting; similarly, no prior-independent mechanism can consistently do better than the second-price auction. In light of this, we design a family of (parametrized) randomized auctions which approximates at least one of these benchmarks: For heterogeneous bidders with regular value distributions, our mechanisms either achieve a good approximation of the expected revenue of an optimal mechanism (which knows the bidders' distributions) or exceeds that of the second-price auction by a certain multiplicative factor. The factor in the latter case naturally trades off with the approximation ratio of the former case. We show that our mechanism is optimal for such a trade-off between the two cases by establishing a matching lower bound. Our result extends to selling $k$ identical items to heterogeneous bidders with an additional $O\big(\ln^2 k\big)$-factor in our trade-off between the two cases.


title: Maximally entangled proton and charged hadron multiplicity in Deep Inelastic Scattering

abstract: We study the proposal by Kharzeev-Levin to determine entanglement entropy in Deep Inelastic Scattering (DIS) from parton distribution functions (PDFs) and to relate the former to the entropy of final state hadrons. We find several uncertainties in the current comparison to data, in particular the overall normalization, the relation between charged versus total hadron multiplicity in the comparison to experimental results as well as different methods to determine the number of partons in Deep Inelastic Scattering. We further provide a comparison to data based on leading order HERA PDF as well as PDFs obtained from an unintegrated gluon distribution subject to next-to-leading order Balitsky-Fadin-Kuraev-Lipatov and Balitsky-Kovchegov evolution. Within uncertainties we find good agreement with H1 data. We provide also predictions for entropy at lower photon virtualities, where non-linear QCD dynamics is expected to become relevant.


title: Drag coefficient of a rising bubble in a shear-thinning fluid using the power-law scheme coupled with a Cahn-Hilliard equation with a variable mobility: A lattice Boltzmann study and comparison with experiment

abstract: This study aims to investigate the behavior of multicomponent fluid flows consisting of Newtonian and non-Newtonian components, especially terminal velocity of a rising bubble in a power-law fluid. A recent lattice Boltzmann (LB) model is extended using power-law scheme to be able to simulate both Newtonian and non-Newtonian fluid flows at high density and viscosity ratios. Also, a variable mobility is introduced in this study to minimize the unphysical error around small bubbles in the domain. A three-component fluid flow system is examined using a constant and variable mobility. It is shown that each component has more stability using variable mobility while constant mobility causes interface dissipation, leading to mass loss gradually. In addition, two test cases including power-law fluid flows driven between two parallel plates are conducted to show the accuracy and capability of the model. To find a grid-independent computational domain, a grid independency test is carried out to show that a 200*400 domain size is suitable for our computations. Then, terminal velocity of a rising bubble is compared to an existing correlation in the literature, indicating that the results are in good agreement with existing study so that average relative error in six different cases is 5.66 %. Also, the simulated examples show good conformity to experimental results over a range of the Reynolds and Eotvos numbers.


title: Deep equilibrium networks are sensitive to initialization statistics

abstract: Deep equilibrium networks (DEQs) are a promising way to construct models which trade off memory for compute. However, theoretical understanding of these models is still lacking compared to traditional networks, in part because of the repeated application of a single set of weights. We show that DEQs are sensitive to the higher order statistics of the matrix families from which they are initialized. In particular, initializing with orthogonal or symmetric matrices allows for greater stability in training. This gives us a practical prescription for initializations which allow for training with a broader range of initial weight scales.


title: Band alignment and interlayer hybridisation in transition metal dichalcogenide/hexagonal boron nitride heterostructures

abstract: In van der Waals heterostructures, the relative alignment of bands between layers, and the resulting band hybridisation, are key factors in determining a range of electronic properties. This work examines these effects for heterostructures of transition metal dichalcogenides (TMDs) and hexagonal boron nitride (hBN), an ubiquitous combination given the role of hBN as an encapsulating material. By comparing results of density functional calculations with experimental angle-resolved photoemission spectroscopy (ARPES) results, we explore the hybridisation between the valence states of the TMD and hBN layers, and show that it introduces avoided crossings between the TMD and hBN bands, with umklapp processes opening `ghost' avoided crossings in individual bands. Comparison between DFT and ARPES spectra for the MoSe$_2$/hBN heterostructure shows that the valence bands of MoSe$_2$ and hBN are significantly further separated in energy in experiment as compared to DFT. We then show that a novel scissor operator can be applied to the hBN valence states in the DFT calculations, to correct the band alignment and enable quantitative comparison to ARPES, explaining avoided crossings and other features of band visibility in the ARPES spectra.


title: Two Remarkably Luminous Galaxy Candidates at $z\approx10-12$ Revealed by JWST

abstract: The first few hundred Myrs at $z>10$ mark the last major uncharted epoch in the history of the Universe, where only a single galaxy (GNz11 at $z\approx11$) is currently spectroscopically confirmed. Here we present a search for luminous $z>10$ galaxies with $JWST$/NIRCam photometry spanning $\approx1-5\mu$m and covering 49 arcmin$^{2}$ from the public Early Release Science programs (CEERS and GLASS). Our most secure candidates are two $M_{\rm{UV}}\approx-21$ systems: GLASS-z12 and GLASS-z10. These galaxies display abrupt $\gtrsim1.8$ mag breaks in their spectral energy distributions, consistent with complete absorption of flux bluewards of Lyman-$\alpha$ that is redshifted to $z=12.4^{+0.1}_{-0.3}$ and $z=10.4^{+0.4}_{-0.5}$. Lower redshift interlopers such as quiescent galaxies with strong Balmer breaks would be comfortably detected at $>5\sigma$ in multiple bands where instead we find no flux. From SED modeling we infer that these galaxies have already built up $\sim 10^9$ solar masses in stars over the $\lesssim300-400$ Myrs after the Big Bang. The brightness of these sources enable morphological constraints. Tantalizingly, GLASS-z10 shows a clearly extended exponential light profile, potentially consistent with a disk galaxy of $r_{\rm{50}}\approx0.7$ kpc. These sources, if confirmed, join GNz11 in defying number density forecasts for luminous galaxies based on Schechter UV luminosity functions, which require a survey area $>10\times$ larger than we have studied here to find such luminous sources at such high redshifts. They extend evidence from lower redshifts for little or no evolution in the bright end of the UV luminosity function into the cosmic dawn epoch, with implications for just how early these galaxies began forming. This, in turn, suggests that future deep $JWST$ observations may identify relatively bright galaxies to much earlier epochs than might have been anticipated.


title: Regret Minimization with Noisy Observations

abstract: In a typical optimization problem, the task is to pick one of a number of options with the lowest cost or the highest value. In practice, these cost/value quantities often come through processes such as measurement or machine learning, which are noisy, with quantifiable noise distributions. To take these noise distributions into account, one approach is to assume a prior for the values, use it to build a posterior, and then apply standard stochastic optimization to pick a solution. However, in many practical applications, such prior distributions may not be available. In this paper, we study such scenarios using a regret minimization model. In our model, the task is to pick the highest one out of $n$ values. The values are unknown and chosen by an adversary, but can be observed through noisy channels, where additive noises are stochastically drawn from known distributions. The goal is to minimize the regret of our selection, defined as the expected difference between the highest and the selected value on the worst-case choices of values. We show that the na\"ive algorithm of picking the highest observed value has regret arbitrarily worse than the optimum, even when $n = 2$ and the noises are unbiased in expectation. On the other hand, we propose an algorithm which gives a constant-approximation to the optimal regret for any $n$. Our algorithm is conceptually simple, computationally efficient, and requires only minimal knowledge of the noise distributions.


title: Early results from GLASS-JWST. III: Galaxy candidates at z$\sim$9-15

abstract: We present the results of a first search for galaxy candidates at z$\sim$9--15 on deep seven-bands NIRCam imaging acquired as part of the GLASS-JWST Early Release Science Program on a flanking field of the Frontier Fields cluster A2744. Candidates are selected via two different renditions of the Lyman-break technique, isolating objects at z$\sim$9-11, and z$\sim$9-15, respectively, supplemented by photometric redshifts obtained with two independent codes. We find five color-selected candidates at z$>$9, plus one additional candidate with photometric redshift z$_{phot}\geq$9. In particular, we identify two bright candidates at $M_{UV}\simeq -21$ that are unambiguously placed at $z\simeq 10.6$ and $z\simeq 12.2$, respectively. The total number of galaxies discovered at $z>9$ is in line with the predictions of a non-evolving LF. The two bright ones at $z>10$ are unexpected given the survey volume, although cosmic variance and small number statistics limits general conclusions. This first search demonstrates the unique power of JWST to discover galaxies at the high redshift frontier. The candidates are ideal targets for spectroscopic follow-up in Cycle$-2$.


title: Design and Simulated Performance of Calorimetry Systems for the ECCE Detector at the Electron Ion Collider

abstract: We describe the design and performance the calorimeter systems used in the ECCE detector design to achieve the overall performance specifications cost-effectively with careful consideration of appropriate technical and schedule risks. The calorimeter systems consist of three electromagnetic calorimeters, covering the combined pseudorapdity range from -3.7 to 3.8 and two hadronic calorimeters. Key calorimeter performances which include energy and position resolutions, reconstruction efficiency, and particle identification will be presented.


title: Simulating Heisenberg Interactions in the Ising Model with Strong Drive Fields

abstract: The time-evolution of an Ising model with large driving fields over discrete time intervals is shown to be reproduced by an effective XXZ-Heisenberg model at leading order in the inverse field strength. For specific orientations of the drive field, the dynamics of the XXX-Heisenberg model is reproduced. These approximate equivalences, valid above a critical driving field strength set by dynamical phase transitions in the Ising model, are expected to enable quantum devices that natively evolve qubits according to the Ising model to simulate more complex systems.


title: All Paths Lead to Rome

abstract: All roads lead to Rome is the core idea of the puzzle game Roma. It is played on an $n \times n$ grid consisting of quadratic cells. Those cells are grouped into boxes of at most four neighboring cells and are either filled, or to be filled, with arrows pointing in cardinal directions. The goal of the game is to fill the empty cells with arrows such that each box contains at most one arrow of each direction and regardless where we start, if we follow the arrows in the cells, we will always end up in the special Roma-cell. In this work, we study the computational complexity of the puzzle game Roma and show that completing a Roma board according to the rules is an \NP-complete task, counting the number of valid completions is #Ptime-complete, and determining the number of preset arrows needed to make the instance \emph{uniquely} solvable is $\Sigma_2^P$-complete. We further show that the problem of completing a given Roma instance on an $n\times n$ board cannot be solved in time $\mathcal{O}\left(2^{{o}(n)}\right)$ under ETH and give a matching dynamic programming algorithm based on the idea of Catalan structures.