title: Persistent Tensors and Multiqudit Entanglement Transformation

abstract: We construct a lower bound of the tensor rank for a new class of tensors, which we call persistent tensors. We present three specific families of persistent tensors, of which the lower bound is tight. We show that there is a chain of degenerations between these three families of minimal-rank persistent tensors that can be used to study the entanglement transformation between them. In addition, we show that these three families of persistent tensors are indeed different generalizations of multiqubit $\rm{W}$ states within multiqudit systems and are geometrically in the orbit closure of multiqudit $\rm{GHZ}$ states. Consequently, we show that one can obtain every one of the generalizations of $\rm{W}$ state from a multiqudit $\rm{GHZ}$ state via asymptotic Stochastic Local Operations and Classical Communication (SLOCC) with rate one. Finally, we extend the obtained lower bound of the tensor rank to direct sums with persistent summands and to even more general combinations of tensors, which we call block pyramidal tensors. As a result, we show that the tensor rank is multiplicative under the Kronecker and tensor products of minimal-rank persistent tensors with the $\rm{GHZ}$ tensor.


title: Trimer quantum spin liquid in a honeycomb array of Rydberg atoms

abstract: Quantum spin liquids are elusive but paradigmatic examples of strongly correlated quantum states that are characterized by long-range quantum entanglement. Recently, the direct signatures of a gapped topological $\mathbb{Z}_2$ spin liquid have been observed in a system of Rydberg atoms arrayed on the ruby lattice. Here, we illustrate the concrete realization of a fundamentally different class of spin liquids in a honeycomb array of Rydberg atoms. Exploring the quantum phase diagram of this system using both density-matrix renormalization group and exact diagonalization simulations, several density-wave-ordered phases are characterized and their origins explained. More interestingly, in the regime where third-nearest-neighbor atoms lie within the Rydberg blockade radius, we find a novel ground state -- with an emergent $\mathrm{U}(1)\times \mathrm{U}(1)$ local symmetry -- formed from superpositions of classical {\it trimer} configurations on the dual triangular lattice. The fidelity of this trimer spin liquid state can be enhanced via dynamical preparation, which we explain by a Rydberg-blockade-based projection mechanism associated with the smooth turnoff of the laser drive. Finally, we discuss the robustness of the trimer spin liquid phase under realistic experimental parameters and demonstrate that our proposal can be readily implemented in current Rydberg atom quantum simulators.

categories: cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el hep-lat quant-ph


title: Littlest Modular Seesaw

abstract: We present the first complete model of the Littlest Modular Seesaw, based on two right-handed neutrinos, within the framework of multiple modular symmetries, justifying the use of multiple moduli fields which take their values at 3 specific stabilizers of $\Gamma_4 \simeq S_4$, including a new phenomenological possibility. Using a semi-analytical approach, we perform a $\chi^2$ analysis of each case and show that good agreement with neutrino oscillation data is obtained, including predictive relations between the leptonic mixing angles and the ratio of light neutrino masses, which non-trivially agree with the experimental values. It is noteworthy that in this very predictive setup, the models fit the global fits of the experimental data remarkably well, both with and without the Super-Kamiokande atmospheric data, for both choices of stabilizers. By extending the model to include a weighton and the double cover group $\Gamma'_4 \simeq S'_4$, we are able to also account for the hierarchy of the charged leptons using modular symmetries, without altering the neutrino predictions.


title: Coherent error threshold for surface codes from Majorana delocalization

abstract: Statistical mechanics mappings provide key insights on quantum error correction. However, existing mappings assume incoherent noise, thus ignoring coherent errors due to, e.g., spurious gate rotations. We map the surface code with coherent errors, taken as $X$- or $Z$-rotations (replacing bit or phase flips), to a two-dimensional (2D) Ising model with complex couplings, and further to a 2D Majorana scattering network. Our mappings reveal both commonalities and qualitative differences in correcting coherent and incoherent errors. For both, the error-correcting phase maps, as we explicitly show by linking 2D networks to 1D fermions, to a $\mathbb{Z}_2$-nontrivial 2D insulator. However, beyond a rotation angle $\phi_\text{th}$, instead of a $\mathbb{Z}_2$-trivial insulator as for incoherent errors, coherent errors map to a Majorana metal. This $\phi_\text{th}$ is the theoretically achievable storage threshold. We numerically find $\phi_\text{th}\approx0.14\pi$. The corresponding bit-flip rate $\sin^2(\phi_\text{th})\approx 0.18$ exceeds the known incoherent threshold $p_\text{th}\approx0.11$.

categories: quant-ph cond-mat.dis-nn cond-mat.mes-hall cond-mat.stat-mech


title: Uncovering a population of gravitational lens galaxies with magnified standard candle SN Zwicky

abstract: Detecting gravitationally lensed supernovae is among the biggest challenges in astronomy. It involves a combination of two very rare phenomena: catching the transient signal of a stellar explosion in a distant galaxy and observing it through a nearly perfectly aligned foreground galaxy that deflects light towards the observer. High-cadence optical observations with the Zwicky Transient Facility, with an unparalleled large field of view, led to the detection of a multiply-imaged Type Ia supernova (SN Ia), ``SN Zwicky", a.k.a. SN 2022qmx. Magnified nearly twenty-five times, the system was found thanks to the ``standard candle" nature of SNe Ia. High-spatial resolution imaging with the Keck telescope resolved four images of the supernova with very small angular separation, corresponding to an Einstein radius of only $\theta_E =0.167"$ and almost identical arrival times. The small $\theta_E$ and faintness of the lensing galaxy is very unusual, highlighting the importance of supernovae to fully characterise the properties of galaxy-scale gravitational lenses, including the impact of galaxy substructures.


title: The Close AGN Reference Survey (CARS). A parsec scale multi-phase outflow in the super-Eddington NLS1 Mrk 1044

abstract: The interaction between Active Galactic Nuclei (AGN) and their host galaxies is scarcely resolved. Narrow-line Seyfert 1 (NLS1) galaxies are believed to represent AGN at early stages of their evolution and allow to observe AGN feeding and feedback processes at high accretion rates. We apply a spectroastrometric analysis to VLT MUSE NFM-AO observations of Mrk 1044, a nearby super-Eddington accreting NLS1. This allows us to map two ionised gas outflows traced by [O$\,$III] which have velocities of $-560\pm20\,{\rm km\:s}^{-1}$ and $-144 \pm 5 \,{\rm km\:s}^{-1}$. Both outflows are spatially unresolved and located close to the galaxy nucleus ($<1\,{\rm pc}$). They have gas densities higher than $10^5\,{\rm cm}^{-3}$, which implies that the BPT diagnostic cannot be used to constrain the underlying ionisation mechanism. We explore whether an expanding shell model can describe the velocity structure of Mrk 1044's unresolved multi-phase outflow. A kinematic analysis suggests that significant turbulence may be present in the ISM around the nucleus, which may lead to a condensation rain, potentially explaining the efficient feeding of Mrk 1044's AGN. We identify an additional ionised gas outflowing component that is spatially resolved. It has a velocity of $-211 \pm 22 \,{\rm km\:s}^{-1}$ and a projected size of $4.6 \pm 0.6 \,{\rm pc}$. Within the innermost 0.5" (160$\,{\rm pc}$) around the nucleus we detect modest star formation hidden by the beam-smeared emission from the outflow, which suggests that Mrk 1044's AGN phase set on recently. We estimate that the multi-phase outflow has been launched $<10^4 \,{\rm yrs}$ ago. It carries enough mass and energy to impact the host galaxy star formation on different spatial scales, highlighting the complexity of the AGN feeding and feedback cycle in its early stages.


title: Untwisting moir\'e physics: Almost ideal bands and fractional Chern insulators in periodically strained monolayer graphene

abstract: Moir\'e systems have emerged in recent years as a rich platform to study strong correlations. Here, we will discuss a simple, experimentally feasible setup based on periodically strained graphene that reproduces several key aspects of twisted moir\'e heterostructures -- but without introducing a twist. We consider a monolayer graphene sheet subject to a $C_2$-breaking periodic strain-induced psuedomagnetic field (PMF) with period $L_M \gg a$, along with a scalar potential of the same period. This system has {\it almost ideal} flat bands with valley-resolved Chern number $\pm 1$, where the deviation from ideal band geometry is analytically controlled and exponentially small in the dimensionless ratio $(L_M/l_B)^2$ where $l_B$ is the magnetic length corresponding to the maximum value of the PMF. Moreover, the scalar potential can tune the bandwidth far below the Coulomb scale, making this a very promising platform for strongly interacting topological phases. Using a combination of strong-coupling theory and self-consistent Hartree fock, we find quantum anomalous Hall states at integer fillings. At fractional filling, exact diagonaliztion reveals a fractional Chern insulator at parameters in the experimentally feasible range. Overall, we find that this system has larger interaction-induced gaps, smaller quasiparticle dispersion, and enhanced tunability compared to twisted graphene systems, even in their ideal limit.


title: The irresistible call of $\tau=i$

abstract: We analyze a large set of modular invariant models of lepton masses and mixing angles, pointing out that many of them prefer to live close to the self-dual point $\tau=i$. We show that in the vicinity of this point a universal behavior naturally emerges, independently from details of the theory such as the finite modular group acting on the lepton multiplets, the weights of the matter multiplets and even the form of the kinetic terms, which are not required to be neither minimal nor flavour blind. The neutrino mass spectrum is normally ordered and universal relations describe the scaling of the physical observables in terms of the parameter $|\tau-i|$.


title: ${\rm S{\scriptsize IM}BIG}$: Mock Challenge for a Forward Modeling Approach to Galaxy Clustering

abstract: Simulation-Based Inference of Galaxies (${\rm S{\scriptsize IM}BIG}$) is a forward modeling framework for analyzing galaxy clustering using simulation-based inference. In this work, we present the ${\rm S{\scriptsize IM}BIG}$ forward model, which is designed to match the observed SDSS-III BOSS CMASS galaxy sample. The forward model is based on high-resolution ${\rm Q{\scriptsize UIJOTE}}$ $N$-body simulations and a flexible halo occupation model. It includes full survey realism and models observational systematics such as angular masking and fiber collisions. We present the "mock challenge" for validating the accuracy of posteriors inferred from ${\rm S{\scriptsize IM}BIG}$ using a suite of 1,500 test simulations constructed using forward models with a different $N$-body simulation, halo finder, and halo occupation prescription. As a demonstration of ${\rm S{\scriptsize IM}BIG}$, we analyze the power spectrum multipoles out to $k_{\rm max} = 0.5\,h/{\rm Mpc}$ and infer the posterior of $\Lambda$CDM cosmological and halo occupation parameters. Based on the mock challenge, we find that our constraints on $\Omega_m$ and $\sigma_8$ are unbiased, but conservative. Hence, the mock challenge demonstrates that ${\rm S{\scriptsize IM}BIG}$ provides a robust framework for inferring cosmological parameters from galaxy clustering on non-linear scales and a complete framework for handling observational systematics. In subsequent work, we will use ${\rm S{\scriptsize IM}BIG}$ to analyze summary statistics beyond the power spectrum including the bispectrum, marked power spectrum, skew spectrum, wavelet statistics, and field-level statistics.


title: Combined X-ray and optical analysis to probe the origin of the plateau emission in $\gamma$-ray bursts afterglows

abstract: A large fraction of gamma-ray bursts (GRBs) shows a plateau phase during the X-ray afterglow emission, whose physical origin is still debated. In this work we define a sample of 30 GRBs with simultaneous X-ray and optical data during and after the plateau phase. Through a time-resolved spectral analysis of the X-ray plateaus, we test the consistency of the unabsorbed optical fluxes with those obtained via X-ray-to-optical spectral extrapolation by assuming a synchrotron spectrum. Combining X-ray with optical data, we find that 63% (19/30) GRBs are compatible with a single synchrotron spectrum thus suggesting that both the optical and X-ray radiations are produced from a single emitting region. For these GRBs we derive the temporal evolution of the break frequency and we compare it with the expectations predicted by several models. For 11/30 GRBs the optical emission is above the predicted range of values extrapolated from the X-rays in at least one temporal bin of the light curve. These GRBs may not be explained with a single zone emission, indicating the necessity of invoking two cooperating processes in order to explain the broad band spectral behaviour during X-ray plateaus. We discuss our findings in the framework of different scenarios invoked to explain the plateau feature, including the energy injection from a spinning-down magnetar and the high latitude emission from a structured jet.


title: Vacuum energy, the Casimir effect, and Newton's non-constant

abstract: We explore two hypotheses. First, the possibility that the quantum vacuum energy density of the Casimir effect contributes to a (local) gravitational vacuum energy density. Second, the possibility that a change in the gravitational coupling implies a change in the cosmological constant. We parametrize these two possibilities in a covariant framework and show that the next generation of Casimir experiments does have a surprisingly good chance of exploring this parameter space.


title: Asymptotic sign free in interacting fermion models

abstract: As an intrinsically-unbiased approach, quantum Monte Carlo (QMC) is of vital importance in understanding correlated phases of matter. Unfortunately, it often suffers notorious sign problem when simulating interacting fermion models. Here, we show for the first time that there exist interacting fermion models whose sign problem becomes less severe for larger system sizes and eventually disappears in the thermodynamic limit, which we dub as "asymptotic sign free". We demonstrate asymptotically-free sign in determinant QMC for various interacting models. Moreover, based on renormalization-group-like ideas we propose a heuristic understanding of the feature of asymptotic sign free. We believe that asymptotic sign free behavior could shed new lights to deepening our understanding of sign problem. More importantly, it can provide a promising way to decipher intriguing physics in correlated models which were conventionally thought not accessible by QMC.


title: Flavor-changing light bosons with accidental longevity

abstract: We consider a model with a complex scalar field that couples to $(e,\mu)$ or $(\mu,\tau)$ within the "longevity" window: $[|m_{l_1} - m_{l_2}|, m_{l_1} + m_{l_2}]$ in which $l_1$ and $l_2$ are the two different charged leptons. Within such a mass window, even a relatively large coupling (e.g. of the size commensurate with the current accuracy/discrepancy in the muon $g-2$ experiment) leads to long lifetimes and macroscopic propagation distance between production and decay points. We propose to exploit several existing neutrino experiments and one future experiment to probe the parameter space of this model. For the $\mu-e$ sector, we exploit the muonium decay branching ratio and the production and decay sequence at the LSND experiment, excluding the parametric region suggested by $g_\mu -2$ anomaly. For the $\tau-\mu$ sector, we analyze three main production mechanisms of scalars at beam dump experiments: the Drell-Yan process, the heavy meson decay, and the muon scattering. We explore the constraints from the past CHARM and NuTeV experiments, and evaluate sensitivity for the proposed beam dump experiment, SHiP. The latter can thoroughly probe the parameter space relevant for the $g_\mu - 2$ anomaly.


title: A Criterion for Strange Metallicity in the Lorenz Ratio

abstract: The Wiedemann-Franz (WF) law, stating that the Lorenz ratio $L = \kappa/(T\sigma)$ between the thermal and electrical conductivities in a metal approaches a universal constant $L_0=\pi^2 k_B^2/ (3 e^2)$ at low temperatures, is often interpreted as a signature of fermionic Landau quasi-particles. In contrast, we show that various models of weakly disordered non-Fermi liquids also obey the WF law at $T \to 0$. Instead, we propose using the leading low-temperature correction to the WF law, $L(T)-L_0$ (proportional to the inelastic scattering rate), to distinguish different types of strange metals. As an example, we demonstrate that in a solvable model of a marginal Fermi liquid, $L(T)-L_0\propto -T$. Using the quantum Boltzmann equation (QBE) approach, we find analogous behavior in a class of marginal- and non-Fermi liquids with a weakly momentum-dependent inelastic scattering. In contrast, in a Fermi liquid, $L(T)-L_0$ is proportional to $-T^2$. This holds even when the resistivity grows linearly with $T$, due to $T-$linear quasi-elastic scattering (as in the case of electron-phonon scattering at temperatures above the Debye frequency). Finally, by exploiting the QBE approach, we demonstrate that the transverse Lorenz ratio, $L_{xy} = \kappa_{xy}/(T\sigma_{xy})$, exhibits the same behavior.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.stat-mech


title: A $\sim$300 pc-sized core of Milky Way dark matter halo constrained from the OGLE micro-lensing sky map

abstract: We report the detection of a 282 $^{+34}_{-31}$ pc-sized core in the center of Milky Way dark matter halo at $68\%$ confidence level by using the micro-lensing event rate sky map data from the Optical Gravitational Lensing Experiment (OGLE) survey. We apply the spacial information of the micro-lensing sky map and model it with the detailed Milky Way dark matter halo Core/Cusp profile, and the fraction of dark matter in the form of Mini Dark Matter Structure (MDMS, $f_{\rm MDMS}=\Omega_{\rm MDMS}/\Omega_{\rm DM}$), e.g. primordial black hole, earth-mass subhalos, floating planets and so on. We find that this sky map can constrain both $f_{\rm MDMS}$ and the core size simultaneously without strong degeneracy while fully considering mass function of Milky Way stellar components from both the bulge and disk.


title: HD 191939 revisited: New and refined planet mass determinations, and a new planet in the habitable zone

abstract: HD 191939 (TOI-1339) is a nearby (d=54pc), bright (V=9mag), and inactive Sun-like star (G9 V) known to host a multi-planet transiting system. Ground-based spectroscopic observations confirmed the planetary nature of the three transiting sub-Neptunes (HD 191939 b, c, and d) originally detected by TESS and were used to measure the masses for planets b and c with 3$\sigma$ precision. These previous observations also reported the discovery of an additional Saturn-mass planet (HD 191939 e) and evidence for a further, very long-period companion (HD 191939 f). Here, we report the discovery of a new non-transiting planet in the system and a refined mass determination of HD 191939 d. The new planet, HD 191939 g, has a minimum mass of 13.5$\pm$2.0 M$_\oplus$ and a period of about 280 d. This period places the planet within the conservative habitable zone of the host star, and near a 1:3 resonance with HD 191939 e. The compilation of 362 radial velocity measurements with a baseline of 677 days from four different high-resolution spectrographs also allowed us to refine the properties of the previously known planets, including a 4.6$\sigma$ mass determination for planet d, for which only a 2$\sigma$ upper limit had been set until now. We confirm the previously suspected low density of HD 191939 d, which makes it an attractive target for attempting atmospheric characterisation. Overall, the planetary system consists of three sub-Neptunes interior to a Saturn-mass and a Uranus-mass planet plus a high-mass long-period companion. This particular configuration has no counterpart in the literature and makes HD 191939 an exceptional multi-planet transiting system with an unusual planet demographic worthy of future observation.


title: Dicke superradiance requires interactions beyond nearest-neighbors

abstract: Photon-mediated interactions within an excited ensemble of emitters can result in Dicke superradiance, where the emission rate is greatly enhanced, manifesting as a high-intensity burst at short times. The superradiant burst is most commonly observed in systems with long-range interactions between the emitters, although the minimal interaction range remains unknown. Here, we put forward a new theoretical method to bound the maximum emission rate by upper bounding the spectral radius of an auxiliary Hamiltonian. We harness this tool to prove that for an arbitrary ordered array with only nearest-neighbor interactions in all dimensions, a superradiant burst is not physically observable. We show that Dicke superradiance requires minimally the inclusion of next-nearest-neighbor interactions. For exponentially decaying interactions, the critical coupling is found to be asymptotically independent of the number of emitters in all dimensions, thereby defining the threshold interaction range where the collective enhancement balances out the decoherence effects. Our findings provide key physical insights to the understanding of collective decay in many-body quantum systems, and the designing of superradiant emission in physical systems for applications such as energy harvesting and quantum sensing.


title: Laser-driven first-order spin reorientation and Verwey phase transitions in the magnetite Fe$_3$O$_4$ beyond the range of thermodynamic equilibrium

abstract: Ultrafast photo-induced phase transitions occurring under the impact of femtosecond laser pulses provide versatile opport unities for switching solids between distinctly-different crystalline, electronic, and magnetic states and thus modify their functional properties in a significant way. In this paper, we report on the laser-induced spin reorientation and Verwey phase transitions in a single crystalline ferrimagnetic magnetite Fe$_3$O$_4$. Using femtosecond optical and magneto-optical pump-probe techniques, we define the range of the initial sample temperatures and laser fluences when partialor complete photo-induced phase transitions occur from a monoclinic insulating to a cubic metallic state with concomitant switching of magnetic anisotropy from the uniaxial to the cubic one. We thus reveal a connection between these phase transitions when driven by femtosecond laser pulses.Using transient linear and quadratic magneto-optical effects, we examine magnetization dynamics launched the switching of the magnetic anisotropy axis. We unveil the presence of the domains under going the laser-induced phase transitions even below the established threshold fluence for the transitions, as well as when the material is initially in the cubic phase. This is the manifestation of the first-order of these both laser-induced phase transitions beyond the range of thermodynamic equilibrium.


title: Stochastic thermodynamic bounds on logical circuit operation

abstract: Using a thermodynamically consistent, mesoscopic model for modern complementary metal-oxide-semiconductor transistors, we study an array of logical circuits and explore how their function is constrained by recent thermodynamic uncertainty relations when operating near thermal energies. For a single NOT gate, we find operating direction-dependent dynamics, and an optimal trade-off between dissipated heat and operation time certainty. For a memory storage device, we find an exponential relationship between the memory retention time and energy required to sustain that memory state. For a clock, we find that the certainty in the cycle time is maximized at biasing voltages near thermal energy, as is the trade-off between this certainty and the heat dissipated per cycle. We demonstrate that a simple control mechanism for the clock leads to a monotonic increase in cycle time certainty with biasing voltage alleviating its degradation at large biasing voltages. These results provide a framework for assessing thermodynamic costs of realistic computing devices, allowing for circuits to be designed and controlled for thermodynamically optimal operation.


title: Cover It Up! Bipartite Graphs Uncover Identifiability in Sparse Factor Analysis

abstract: Despite the popularity of factor models with sparse loading matrices, little attention has been given to formally address identifiability of these models beyond standard rotation-based identification such as the positive lower triangular constraint. To fill this gap, we present a counting rule on the number of nonzero factor loadings that is sufficient for achieving generic uniqueness of the variance decomposition in the factor representation. This is formalized in the framework of sparse matrix spaces and some classical elements from graph and network theory. Furthermore, we provide a computationally efficient tool for verifying the counting rule. Our methodology is illustrated for real data in the context of post-processing posterior draws in Bayesian sparse factor analysis.


title: A Systematic Comparison of Galaxy Cluster Temperatures Measured with NuSTAR and Chandra

abstract: Temperature measurements of galaxy clusters are used to determine their masses, which in turn are used to determine cosmological parameters. However, systematic differences between the temperatures measured by different telescopes imply a significant source of systematic uncertainty on such mass estimates. We perform the first systematic comparison between cluster temperatures measured with Chandra and NuSTAR. This provides a useful contribution to the effort of cross-calibrating cluster temperatures due to the harder response of NuSTAR compared with most other observatories. We measure average temperatures for 8 clusters observed with NuSTAR and Chandra. We fit the NuSTAR spectra in a hard (3-10 keV) energy band, and the Chandra spectra in both the hard and a broad (0.6-9 keV) band. We fit a power-law cross-calibration model to the resulting temperatures. At a Chandra temperature of 10 keV, the average NuSTAR temperature was $(10.5 \pm 3.7)\%$ and $(15.7 \pm 4.6)\%$ lower than Chandra for the broad and hard band fits respectively. We explored the impact of systematics from background modelling and multiphase temperature structure of the clusters, and found that these did not affect our results. Our sample are primarily merging clusters with complex thermal structures so are not ideal calibration targets. However, given the harder response of NuSTAR it would be expected to measure a higher average temperature than Chandra for a non-isothermal cluster, so we interpret our measurement as a lower limit on the difference in temperatures between NuSTAR and Chandra.


title: A configuration interaction correction on top of pair coupled cluster doubles

abstract: Numerous numerical studies have shown that geminal-based methods are a promising direction to model strongly correlated systems with low computational costs. Several strategies have been introduced to capture the missing dynamical correlation effects, which typically exploit \textit{a posteriori} corrections to account for correlation effects associated with broken-pair states or inter-geminal correlations. In this article, we scrutinize the accuracy of the pair coupled cluster doubles (pCCD) method extended by configuration interaction (CI) theory. Specifically, we benchmark various CI models, including, at most double excitations against selected CC corrections as well as conventional single-reference CC methods. A simple Davidson correction is also tested. The accuracy of the proposed pCCD-CI approaches is assessed for challenging small model systems such as the \ce{N2} and \ce{F2} dimers and various di- and triatomic actinide-containing compounds. In general, the proposed CI methods considerably improve spectroscopic constants compared to the conventional CCSD approach, provided a Davidson correction is included in the theoretical model. At the same time, their accuracy lies between the linearized frozen pCCD and frozen pCCD variants.


title: Measuring bar pattern speeds from single simulation snapshots

abstract: We describe methods to measure simultaneously the orientation angle $\psi$ and pattern speed $\Omega$ from single snapshots of simulated barred galaxies. Unlike previous attempts, our approach is unbiased, precise, and consistent in the sense that $\psi=\int\Omega\mathrm{d}t$. It can be extended to obtain the rate and axis of rotation, i.e. the vector $\boldsymbol{\Omega}$. We provide computer code implementing our method.


title: An Empirical Quantile Estimation Approach to Nonlinear Optimization Problems with Chance Constraints

abstract: We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. We demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.


title: Towards Inter-character Relationship-driven Story Generation

abstract: In this paper, we introduce the task of modeling interpersonal relationships for story generation. For addressing this task, we propose Relationships as Latent Variables for Story Generation, (ReLiSt). ReLiSt generates stories sentence by sentence and has two major components - a relationship selector and a story continuer. The relationship selector specifies a latent variable to pick the relationship to exhibit in the next sentence and the story continuer generates the next sentence while expressing the selected relationship in a coherent way. Our automatic and human evaluations demonstrate that ReLiSt is able to generate stories with relationships that are more faithful to desired relationships while maintaining the content quality. The relationship assignments to sentences during inference bring interpretability to ReLiSt.


title: Semi-Supervised Domain Adaptation for Cross-Survey Galaxy Morphology Classification and Anomaly Detection

abstract: In the era of big astronomical surveys, our ability to leverage artificial intelligence algorithms simultaneously for multiple datasets will open new avenues for scientific discovery. Unfortunately, simply training a deep neural network on images from one data domain often leads to very poor performance on any other dataset. Here we develop a Universal Domain Adaptation method DeepAstroUDA, capable of performing semi-supervised domain alignment that can be applied to datasets with different types of class overlap. Extra classes can be present in any of the two datasets, and the method can even be used in the presence of unknown classes. For the first time, we demonstrate the successful use of domain adaptation on two very different observational datasets (from SDSS and DECaLS). We show that our method is capable of bridging the gap between two astronomical surveys, and also performs well for anomaly detection and clustering of unknown data in the unlabeled dataset. We apply our model to two examples of galaxy morphology classification tasks with anomaly detection: 1) classifying spiral and elliptical galaxies with detection of merging galaxies (three classes including one unknown anomaly class); 2) a more granular problem where the classes describe more detailed morphological properties of galaxies, with the detection of gravitational lenses (ten classes including one unknown anomaly class).


title: Constructing endperiodic loxodromics of infinite-type arc graphs

abstract: We give general conditions to produce endperiodic homeomorphisms that act loxodromically on various arc graphs of infinite-type surfaces.


title: Quantum phase transitions in non-Hermitian $\mathcal{P}\mathcal{T}$-symmetric transverse-field Ising spin chains

abstract: We present a theoretical study of quantum phases and quantum phase transitions occurring in non-Hermitian $\mathcal{P}\mathcal{T}$-symmetric superconducting qubits chains described by a transverse-field Ising spin model. A non-Hermitian part of the Hamiltonian is implemented via imaginary staggered \textit{longitudinal } magnetic field, which corresponds to a local staggered gain and loss terms. By making use of a direct numerical diagonalization of the Hamiltonian for spin chains of a finite size $N$, we explore the dependencies of the energy spectrum, including the energy difference between the first excited and the ground states, the spatial correlation function of local polarization ($z$-component of local magnetization) on the adjacent spins interaction strength $J$ and the local gain (loss) parameter $\gamma$. A scaling procedure for the coherence length $\xi$ allows us to establish a complete quantum phase diagram of the system. We obtain two quantum phases for $J<0$, namely, $\mathcal{P}\mathcal{T}$-symmetry broken antiferromagnetic state and $\mathcal{P}\mathcal{T}$-symmetry preserved paramagnetic state, and the quantum phase transition line between them is the line of exception points. For $J>0$ the $\mathcal{P}\mathcal{T}$-symmetry of the ground state is retained in a whole region of parameter space of $J$ and $\gamma$, and a system shows \textit{two} intriguing quantum phase transitions between ferromagnetic and paramagnetic states for a fixed parameter $\gamma > 1$. We also provide the qualitative quantum phase diagram $\gamma-J$ derived in the framework of the Bethe-Peierls approximation that is in a good accord with numerically obtained results.

categories: quant-ph cond-mat.mes-hall cond-mat.stat-mech cond-mat.str-el cond-mat.supr-con


title: On the detection of synthetic images generated by diffusion models

abstract: Over the past decade, there has been tremendous progress in creating synthetic media, mainly thanks to the development of powerful methods based on generative adversarial networks (GAN). Very recently, methods based on diffusion models (DM) have been gaining the spotlight. In addition to providing an impressive level of photorealism, they enable the creation of text-based visual content, opening up new and exciting opportunities in many different application fields, from arts to video games. On the other hand, this property is an additional asset in the hands of malicious users, who can generate and distribute fake media perfectly adapted to their attacks, posing new challenges to the media forensic community. With this work, we seek to understand how difficult it is to distinguish synthetic images generated by diffusion models from pristine ones and whether current state-of-the-art detectors are suitable for the task. To this end, first we expose the forensics traces left by diffusion models, then study how current detectors, developed for GAN-generated images, perform on these new synthetic images, especially in challenging social-networks scenarios involving image compression and resizing. Datasets and code are available at github.com/grip-unina/DMimageDetection.


title: Thermal Simulation of Millimetre Wave Ablation of Geological Materials

abstract: This work is concerned with the numerical simulation of ablation of geological materials using a millimetre wave source. To this end, a new mathematical model is developed for a thermal approach to the problem, allowing for large scale simulations, whilst being able to include the strong temperature dependence of material parameters to ensure accurate modelling of power input into the rock. The model presented is implemented within an adaptive meshing framework, such that resolution can be placed where needed, for example at the borehole wall, to further improve the computational efficiency of large scale simulations. This approach allows for both the heating of the rock, and the removal of evaporated material, allowing rate of penetration and the shape of the resulting borehole to be quantified. The model is validated against experimental results, which indicates that the approach can accurately predict temperatures, and temperature gradients within the rock. The validated model is then exercised to obtain initial results demonstrating its capabilities for simulating the millimetre wave drilling process. The effects of the conditions at the surface of the rock are investigated, highlighting the importance of understanding the physical processes which occur between the wave guide and the rock. Additionally, the absorptivity of the rock, and the impact this has on the evaporation behaviour is considered. Simulations are carried out both for isotropic rock, and also for a multi-strata configuration. It is found that strata between similar rock types, such as granite and basalt, absorptive properties pose little problem for uniform drilling. However, larger variations in material parameters are shown to have strong implications on the evaporation behaviour of the wellbore, and hence the resulting structure.


title: Higher-Order Nodal Hinge States in Doped Superconducting Topological Insulator

abstract: Doped strong topological insulators are one of the most promising candidates to realize a fully gapped three-dimensional topological superconductor (TSC). In this letter, we revisit this system and reveal a possibility for higher-order topology which was previously missed. We find that over a finite-range of doping, the Fu-Berg superconducting pairing can give rise to both Majorana surface states, and nodal hinge states. Interestingly, we observe the coexistence of surface and hinge modes in the superconducting state only when there are both bulk and surface Fermi-surfaces in the normal state. Also, we find that the hinge modes can appear for normal states consisting of doped strong or weak topological insulators. In summary, this work may allow for the discovery of superconducting hinge modes in a well explored class of materials, i.e., doped strong or weak topological insulators.


title: Reduce, Reuse, Recycle: Improving Training Efficiency with Distillation

abstract: Methods for improving the efficiency of deep network training (i.e. the resources required to achieve a given level of model quality) are of immediate benefit to deep learning practitioners. Distillation is typically used to compress models or improve model quality, but it's unclear if distillation actually improves training efficiency. Can the quality improvements of distillation be converted into training speed-ups, or do they simply increase final model quality with no resource savings? We conducted a series of experiments to investigate whether and how distillation can be used to accelerate training using ResNet-50 trained on ImageNet and BERT trained on C4 with a masked language modeling objective and evaluated on GLUE, using common enterprise hardware (8x NVIDIA A100). We found that distillation can speed up training by up to 1.96x in ResNet-50 trained on ImageNet and up to 1.42x on BERT when evaluated on GLUE. Furthermore, distillation for BERT yields optimal results when it is only performed for the first 20-50% of training. We also observed that training with distillation is almost always more efficient than training without distillation, even when using the poorest-quality model as a teacher, in both ResNet-50 and BERT. Finally, we found that it's possible to gain the benefit of distilling from an ensemble of teacher models, which has O(n) runtime cost, by randomly sampling a single teacher from the pool of teacher models on each step, which only has a O(1) runtime cost. Taken together, these results show that distillation can substantially improve training efficiency in both image classification and language modeling, and that a few simple optimizations to distillation protocols can further enhance these efficiency improvements.


title: TOE: A Grid-Tagging Discontinuous NER Model Enhanced by Embedding Tag/Word Relations and More Fine-Grained Tags

abstract: So far, discontinuous named entity recognition (NER) has received increasing research attention and many related methods have surged such as hypergraph-based methods, span-based methods, and sequence-to-sequence (Seq2Seq) methods, etc. However, these methods more or less suffer from some problems such as decoding ambiguity and efficiency, which limit their performance. Recently, grid-tagging methods, which benefit from the flexible design of tagging systems and model architectures, have shown superiority to adapt for various information extraction tasks. In this paper, we follow the line of such methods and propose a competitive grid-tagging model for discontinuous NER. We call our model TOE because we incorporate two kinds of Tag-Oriented Enhancement mechanisms into a state-of-the-art (SOTA) grid-tagging model that casts the NER problem into word-word relationship prediction. First, we design a Tag Representation Embedding Module (TREM) to force our model to consider not only word-word relationships but also word-tag and tag-tag relationships. Concretely, we construct tag representations and embed them into TREM, so that TREM can treat tag and word representations as queries/keys/values and utilize self-attention to model their relationships. On the other hand, motivated by the Next-Neighboring-Word (NNW) and Tail-Head-Word (THW) tags in the SOTA model, we add two new symmetric tags, namely Previous-Neighboring-Word (PNW) and Head-Tail-Word (HTW), to model more fine-grained word-word relationships and alleviate error propagation from tag prediction. In the experiments of three benchmark datasets, namely CADEC, ShARe13 and ShARe14, our TOE model pushes the SOTA results by about 0.83%, 0.05% and 0.66% in F1, demonstrating its effectiveness.


title: A sufficient family of necessary inequalities for the compatibility of quantum marginals

abstract: The quantum marginal problem is concerned with characterizing which collections of quantum states on different subsystems are compatible in the sense that they are the marginals of some multipartite quantum state. Presented here is a countable family of inequalities, each of which is necessarily satisfied by any compatible collection of quantum states. Additionally, this family of inequalities is shown to be sufficient: every incompatible collection of quantum states will violate at least one inequality belonging to the family.


title: Constraints on Electron Acceleration in Gamma-Ray Bursts Afterglows from Radio Peaks

abstract: Studies of gamma-ray bursts (GRBs) and their multi-wavelength afterglows have led to insights in electron acceleration and emission properties from relativistic, high-energy astrophysical sources. Broadband modeling across the electromagnetic spectrum has been the primary means of investigating the physics behind these sources, although independent diagnostic tools have been developed to inform and corroborate assumptions made in particle acceleration simulations and broadband studies. We present a methodology to constrain three physical parameters related to electron acceleration in GRB blast waves: the fraction of shock energy in electrons, $\epsilon_e$; the fraction of electrons that gets accelerated into a power-law distribution of energies, $\xi_e$; and the minimum Lorentz factor of the accelerated electrons, $\gamma_m$. These parameters are constrained by observations of the peaks in radio afterglow light curves and spectral energy distributions. From a sample of 49 radio afterglows, we are able to find narrow distributions for these parameters, hinting at possible universality of the blast wave microphysics, although observational bias could play a role in this. Using radio peaks and considerations related to the prompt gamma-ray emission efficiency, we constrain the allowed parameter ranges for both $\epsilon_e$ and $\xi_e$ to within about one order of magnitude, $0.01\lesssim\epsilon_e\lesssim0.2$ and $0.1\lesssim\xi_e\lesssim1$. Such stringent constraints are inaccessible for $\xi_e$ from broadband studies due to model degeneracies.


title: Bounds in simple hexagonal lattice and classification of 11-stick knots

abstract: The stick number and the edge length of a knot type in the simple hexagonal lattice (sh-lattice) are the minimal numbers of sticks and edges required, respectively, to construct a knot of the given type in sh-lattice. By introducing a linear transformation between lattices, we prove that for any given knot both values in the sh-lattice are strictly less than the values in the cubic lattice. Finally, we show that the only non-trivial 11-stick knots in the sh-lattice are the trefoil knot ($3_1$) and the figure-eight knot ($4_1$).


title: Learning to Solve Voxel Building Embodied Tasks from Pixels and Natural Language Instructions

abstract: The adoption of pre-trained language models to generate action plans for embodied agents is a promising research strategy. However, execution of instructions in real or simulated environments requires verification of the feasibility of actions as well as their relevance to the completion of a goal. We propose a new method that combines a language model and reinforcement learning for the task of building objects in a Minecraft-like environment according to the natural language instructions. Our method first generates a set of consistently achievable sub-goals from the instructions and then completes associated sub-tasks with a pre-trained RL policy. The proposed method formed the RL baseline at the IGLU 2022 competition.


title: A control-oriented wind turbine dynamic simulation framework which resolves local atmospheric conditions

abstract: Wind turbines may experience local weather perturbation, which is not taken into account by the commonly-used wind turbine simulation packages. Without this information, it is extremely challenging to evaluate the controller performance with regard to the effect of the variation of local atmospheric conditions. On the other side, it is too late and costly to wait until field test time. To fill this gap, in this paper, we develop a control-oriented turbine dynamic simulation framework to evaluate the controller performance considering the perturbation of local atmospheric conditions. This goal is achieved by integrating an internal wind turbine (IWT) model in the Weather Research and Forecasting (WRF) simulation tool. The proposed framework is implemented on a 5MW reference wind turbine, where the effects of the local atmospheric conditions are illustrated. The controller performance results are compared with those derived from the Fatigue, Aerodynamics, Structures,and Turbulence (FAST) simulator as a validation. Simulation results show that the proposed WRF-IWT model is able to capture the turbine dynamics and controller performance regarding the perturbation of the complex local atmospheric conditions. The proposed framework can be leveraged to assist in designing controllers, which is more applicable to real-world wind conditions.


title: Proton and Helium Heating by Cascading Turbulence in a Low-beta Plasma

abstract: How ions are energized and heated is a fundamental problem in the study of energy dissipation in magnetized plasmas. In particular, the heating of heavy ions (including ${}^{4}\mathrm{He}^{2+}$, ${}^{3}\mathrm{He}^{2+}$ and others) has been a constant concern for understanding the microphysics of impulsive solar flares. In this article, via two-dimensional hybrid-kinetic Particle-in-Cell simulations, we study the heating of Helium ions (${}^{4}\mathrm{He}^{2+}$) by turbulence driven by cascading waves launched at large scales from the left-handed polarized Helium ion cyclotron wave branch of a multi-ion plasma composed of electrons, protons, and Helium ions. We find significant parallel (to the background magnetic field) heating for both Helium ions and protons due to the formation of beams and plateaus in their velocity distribution functions along the background magnetic field. The heating of Helium ions in the direction perpendicular to the magnetic field starts with a lower rate than that in the parallel direction, but overtakes the parallel heating after a few hundreds of the proton gyro-periods due to cyclotron resonances with mainly obliquely propagating waves induced by the cascade of injected Helium ion cyclotron waves at large scales. There is however little evidence for proton heating in the perpendicular direction due to the absence of left-handed polarized cyclotron waves near the proton cyclotron frequency. Our results are useful for understanding the preferential heating of ${}^{3}\mathrm{He}$ and other heavy ions in the ${}^{3}\mathrm{He}$-rich solar energetic particle events, in which Helium ions play a crucial role as a species of background ions regulating the kinetic plasma behavior.


title: Accelerated Design of Chalcogenide Glasses through Interpretable Machine Learning for Composition Property Relationships

abstract: Chalcogenide glasses possess several outstanding properties that enable several ground breaking applications, such as optical discs, infrared cameras, and thermal imaging systems. Despite the ubiquitous usage of these glasses, the composition property relationships in these materials remain poorly understood. Here, we use a large experimental dataset comprising approx 24000 glass compositions made of 51 distinct elements from the periodic table to develop machine learning models for predicting 12 properties, namely, annealing point, bulk modulus, density, Vickers hardness, Littleton point, Youngs modulus, shear modulus, softening point, thermal expansion coefficient, glass transition temperature, liquidus temperature, and refractive index. These models, by far, are the largest for chalcogenide glasses. Further, we use SHAP, a game theory based algorithm, to interpret the output of machine learning algorithms by analyzing the contributions of each element towards the models prediction of a property. This provides a powerful tool for experimentalists to interpret the models prediction and hence design new glass compositions with targeted properties. Finally, using the models, we develop several glass selection charts that can potentially aid in the rational design of novel chalcogenide glasses for various applications.


title: Towards Better Out-of-Distribution Generalization of Neural Algorithmic Reasoning Tasks

abstract: In this paper, we study the OOD generalization of neural algorithmic reasoning tasks, where the goal is to learn an algorithm (e.g., sorting, breadth-first search, and depth-first search) from input-output pairs using deep neural networks. First, we argue that OOD generalization in this setting is significantly different than common OOD settings. For example, some phenomena in OOD generalization of image classifications such as \emph{accuracy on the line} are not observed here, and techniques such as data augmentation methods do not help as assumptions underlying many augmentation techniques are often violated. Second, we analyze the main challenges (e.g., input distribution shift, non-representative data generation, and uninformative validation metrics) of the current leading benchmark, i.e., CLRS \citep{deepmind2021clrs}, which contains 30 algorithmic reasoning tasks. We propose several solutions, including a simple-yet-effective fix to the input distribution shift and improved data generation. Finally, we propose an attention-based 2WL-graph neural network (GNN) processor which complements message-passing GNNs so their combination outperforms the state-of-the-art model by a 3% margin averaged over all algorithms. Our code is available at: \url{https://github.com/smahdavi4/clrs}.


title: The effects of stellar rotation along the main sequence of the 100 Myr old massive cluster NGC 1850

abstract: Young star clusters enable us to study the effects of stellar rotation on an ensemble of stars of the same age and across a wide range in stellar mass and are therefore ideal targets for understanding the consequences of rotation on stellar evolution. We combine MUSE spectroscopy with HST photometry to measure the projected rotational velocities (Vsini) of 2,184 stars along the split main sequence and on the main sequence turn-off (MSTO) of the 100 Myr-old massive (10^5 M_sun) star cluster NGC 1850 in the Large Magellanic Cloud. At fixed magnitude, we observe a clear correlation between Vsini and colour, in the sense that fast rotators appear redder. The average Vsini values for stars on the blue and red branches of the split main sequence are ~100 km/s and ~200 km/s, respectively. The values correspond to about 25-30% and 50-60% of the critical rotation velocity and imply that rotation rates comparable to those observed in field stars of similar masses can explain the split main sequence. Our spectroscopic sample contains a rich population of ~200 fast rotating Be stars. The presence of shell features suggests that 23% of them are observed through their decretion disks, corresponding to a disk opening angle of 15 degrees. These shell stars can significantly alter the shape of the MSTO, hence care should be taken when interpreting this photometric feature. Overall, our findings impact our understanding of the evolution of young massive clusters and provide new observational constraints for testing stellar evolutionary models.


title: What we talk about when we talk about physics problem solving

abstract: I am a second-year cognitive science major, and as a student who has completed my physical science distributive requirements, I will likely never again come across Gauss's law. So why do I feel that the time and effort I devoted to solving Gauss's law problems was worth it? Partly, I inherently enjoy the learning process and the new perspective on the physical world I have acquired by understanding electromagnetism. But I was also inspired by the ways in which physics problems train the mind in effective problem-solving strategies. (Of course I was -- I am a cognitive science major!) Two themes emerged as I reflected on this realization. First, physics problems serve as useful toy models for more complex problems outside of physics, training us in broadly transferable problem-solving skills. Second, the physics problem-solving process invites us to reflect on our unique cognitive and affective processes. These themes are interconnected and complimentary. An improved metacognitive understanding of our minds facilitates solving progressively more complex problems, and the act of solving increasingly difficult problems provides further insight into our minds. In what follows, Professor Zosia Krusberg and I consider nine general lessons offered by the physics problem-solving process.


title: Thermals Forecast with Temps

abstract: This paper describes a simple method for estimating the strength of a thermal updraft from a Temp (i.e., a SkewT, Emagram, or similar) showing the temperature and dew point profile of the lower atmosphere. The data of the Temp can come from relevant weather models or from radiosondes.


title: Exploiting Kronecker structure in exponential integrators: fast approximation of the action of $\varphi$-functions of matrices via quadrature

abstract: In this article, we propose an algorithm for approximating the action of $\varphi-$functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the $\varphi-$functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive \emph{a priori} bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.


title: A Converse for Fault-tolerant Quantum Computation

abstract: As techniques for fault-tolerant quantum computation keep improving, it is natural to ask: what is the fundamental lower bound on redundancy? In this paper, we obtain a lower bound on the redundancy required for $\epsilon$-accurate implementation of a large class of operations that includes unitary operators. For the practically relevant case of sub-exponential depth and sub-linear gate size, our bound on redundancy is tighter than the known lower bounds. We obtain this bound by connecting fault-tolerant computation with a set of finite blocklength quantum communication problems whose accuracy requirements satisfy a joint constraint. The lower bound on redundancy obtained here leads to a strictly smaller upper bound on the noise threshold for non-degradable noise. Our bound directly extends to the case where noise at the outputs of a gate are non-i.i.d. but noise across gates are i.i.d.


title: Current mean values in the XYZ model

abstract: The XYZ model is an integrable spin chain which has an infinite set of conserved charges, but it lacks a global $U(1)$-symmetry. We consider the current operators, which describe the flow of the conserved quantities in this model. We derive an exact result for the current mean values, valid for any eigenstate in a finite volume with periodic boundary conditions. This result can serve as a basis for studying the transport properties of this model within Generalized Hydrodynamics.


title: A deletion-contraction long exact sequence for chromatic symmetric homology

abstract: Crew and Spirklt generalize Stanley's chromatic symmetric function to vertex-weighted graphs. One of the primary motivations for extending the chromatic symmetric function to vertex-weighted graphs is the existence of a deletion-contraction relation in this setting, which, as known, holds for the chromatic polynomial, but doesn't hold for the chromatic symmetric function. In this paper we find a categorification of their new invariant extending the definition of chromatic symmetric homology to vertex-weighted graphs. We prove the existence of a deletion-contraction long exact sequence for chromatic symmetric homology which lifts the deletion-contraction relation that holds for the extension of Crew and Spirklt. Moreover, the new categorification gives a useful computational tool and allow us to answer two questions left open by Chandler, Sazdanovic, Stella and Yip. In particular, we prove that, for a graph G with $n$ vertices, the maximal index with nonzero homology is not greater that $n$ - 1. Moreover, we show that the homology is non-trivial for all the indices between the minimum and the maximum with this property.


title: MHITNet: a minimize network with a hierarchical context-attentional filter for segmenting medical ct images

abstract: In the field of medical CT image processing, convolutional neural networks (CNNs) have been the dominant technique.Encoder-decoder CNNs utilise locality for efficiency, but they cannot simulate distant pixel interactions properly.Recent research indicates that self-attention or transformer layers can be stacked to efficiently learn long-range dependencies.By constructing and processing picture patches as embeddings, transformers have been applied to computer vision applications. However, transformer-based architectures lack global semantic information interaction and require a large-scale training dataset, making it challenging to train with small data samples. In order to solve these challenges, we present a hierarchical contextattention transformer network (MHITNet) that combines the multi-scale, transformer, and hierarchical context extraction modules in skip-connections. The multi-scale module captures deeper CT semantic information, enabling transformers to encode feature maps of tokenized picture patches from various CNN stages as input attention sequences more effectively. The hierarchical context attention module augments global data and reweights pixels to capture semantic context.Extensive trials on three datasets show that the proposed MHITNet beats current best practises


title: Photonic Dirac Waveguides

abstract: The Dirac equation is a paradigmatic model that describes a range of intriguing properties of relativistic spin-1/2 particles, from the existence of antiparticles to Klein tunneling. However, the Dirac-like equations have found application far beyond its original scope, and has been used to comprehend the properties of graphene and topological phases of matter. In the field of photonics, the opportunity to emulate Dirac physics has also enabled topological photonic insulators. In this paper, we demonstrate that judiciously engineered synthetic potentials in photonic Dirac systems can offer physical properties beyond both the elementary and quasi-particles, and topological realms. Specifically, we introduce a new class of optical Dirac waveguides, whose guided electromagnetic modes are endowed with pseudo-spin degree of freedom. Pseudo-spin coupled with the ability to engineer synthetic gauge potentials acting on it, enables control over the guided modes which is unattainable in conventional optical waveguides. In particular, we use a silicon nanophotonic metasurface that supports pseudo-spin degree of freedom as a testing platform to predict and experimentally confirm a spin-full nature of the Dirac waveguides. We also demonstrate that, for suitable trapping potentials, the guided modes exhibit spin-dependent field distributions, which gives rise to their distinct transport and radiative properties. Thereby, the Dirac waveguides manifest spin-dependent radiative lifetimes - the non-Hermitian spin-Hall effect - and open new avenues for spin-multiplexing, controlling characteristics of guided optical modes, and tuning light-matter interactions with photonic pseudo-spins.


title: Importance of Sample Selection in Exoplanet Atmosphere Population Studies

abstract: Understanding planet formation requires robust population studies, which are designed to reveal trends in planet properties. In this work, we aim to determine if different methods for selecting populations of exoplanets for atmospheric characterization with JWST could influence population-level inferences. We generate three hypothetical surveys of super-Earths/sub-Neptunes, each spanning a similar radius-insolation flux space. The survey samples are constructed based on three different selection criteria (evenly-spaced-by-eye, binned, and a quantitative selection function). Using an injection-recovery technique, we test how robustly individual-planet atmospheric parameters and population-level parameters can be retrieved. We find that all three survey designs result in equally suitable targets for individual atmospheric characterization, but not equally suitable targets for constraining population parameters. Only samples constructed with a quantitative method or that are sufficiently evenly-spaced-by-eye result in robust population parameter constraints. Furthermore, we find that the sample with the best targets for individual atmospheric study does not necessarily result in the best constrained population parameters. The method of sample selection must be considered. We also find that there may be large variability in population-level results with a sample that is small enough to fit in a single JWST cycle ($\sim$12 planets), suggesting that the most successful population-level analyses will be multi-cycle. Lastly, we infer that our exploration of sample selection is limited by the small number of transiting planets with measured masses around bright stars. Our results can guide future development of programs that aim to determine underlying trends in exoplanet atmospheric properties and, by extension, formation and evolution processes.


title: Refining the Masses and Radii of the Star Kepler-33 and its Five Transiting Planets

abstract: Kepler-33 hosts five validated transiting planets ranging in period from 5 to 41 days. The planets are in nearly co-planar orbits and exhibit remarkably similar (appropriately scaled) transit durations indicative of similar impact parameters. The outer three planets have radii of $3.5\lesssim R_{\rm p}/R_\oplus\lesssim4.7$ and are closely-packed dynamically, and thus transit timing variations can be observed. Photodynamical analysis of transit timing variations provide $2\sigma$ upper bounds on the eccentricity of the orbiting planets (ranging from $<0.02$ to $<0.2$) and the mean density of the host-star ($0.39_{-0.02}^{+0.01}\,{\rm g/cm^3}$). We combine \emph{Gaia} Early Data Release 3 parallax observations, the previously reported host-star effective temperature and metallicity, and our photodynamical model to refine properties of the host-star and the transiting planets. Our analysis yields well-constrained masses for Kepler-33~e ($6.6_{-1.0}^{+1.1}\,M_\oplus$) and f ($8.2_{-1.2}^{+1.6}\,M_\oplus$) along with $2\sigma$ upper limits for planets c ($<19\,M_\oplus$) and d ($<8.2\,M_\oplus$). We confirm the reported low bulk densities of planet d ($<0.4\,{\rm g/cm^3}$), e ($0.8\pm0.1\,{\rm g/cm^3}$), and f ($0.7\pm0.1\,{\rm g/cm^3}$). Based on comparisons with planetary evolution models, we find that Kepler-33~e and f exhibit relatively high envelope mass fractions of $f_{\rm env}=7.0_{-0.5}^{+0.6}\%$ and $f_{\rm env}=10.3\pm0.6\%$, respectively. Assuming a mass for planet d $\sim4\,M_\oplus$ suggests that it has $f_{\rm env}\gtrsim12\%$.


title: Magnetically Dominated Disks in Tidal Disruption Events and Quasi-Periodic Eruptions

abstract: The classical radiation pressure instability has been a persistent theoretical feature of thin, radiatively efficient accretion disks with accretion rates 1 to 100 per cent of the Eddington rate. But there is only limited evidence of its occurrence in nature: rapid heartbeat oscillations of a few X-ray binaries and now, perhaps, the new class of hourly X-ray transients called quasi-periodic eruptions (QPEs). The accretion disks formed in tidal disruption events (TDEs) have been observed to peacefully trespass through the range of unstable accretion rates without exhibiting any clear sign of the instability. We try to explain the occurrence or otherwise of this instability in these systems, by constructing steady state 1D models of thin magnetic accretion disks. The local magnetic pressure in the disk is assumed to be dominated by toroidal fields arising from a dynamo sourced by magneto-rotational instability (MRI). We choose a physically motivated criterion of MRI saturation, validated by recent magnetohydrodynamic simulations, to determine the strength of magnetic pressure in the disk. The resulting magnetic pressure support efficiently shrinks: (1) the parameter space of unstable mass accretion rates, explaining the absence of instability in systems such as TDEs and (2) the range of unstable radii in the inner accretion disk, which can shorten the quasi-periods of instability limit-cycles by more than three orders of magnitude, explaining the observed periods ( a few hrs) of QPEs. In addition to examining stability properties of strongly magnetized disks, we predict other observational signatures such as spectral hardening factors and jet luminosities to test the compatibility of our disk models with observations of apparently stable TDE disks.


title: Explicit-Implicit Domain Splitting for Two Phase Flows with Phase Transition

abstract: Two phase flows that include phase transition, especially phase creation, with a sharp interface remain a challenging task for numerics. We consider the isothermal Euler equations with phase transition between a liquid and a vapor phase. The phase interface is modeled as a sharp interface and the mass transfer across the phase boundary is modeled by a kinetic relation. Existence and uniqueness results were proven in Ref. \cite{Hantke2019a}. Using sharp interfaces for simulating nucleation and cavitation results in the grid containing tiny cells that are several orders of magnitude smaller than the remaining grid cells. This forces explicit time stepping schemes to take tiny time steps on these cells. As a remedy we suggest an explicit implicit domain splitting where the majority of the grid cells is treated explicitly and only the neighborhood of the tiny cells is treated implicitly. We use dual time stepping to solve the resulting small implicit systems. Our numerical results indicate that the new scheme is robust and provides significant speed-up compared to a fully explicit treatment.


title: Tree Representations of Brain Structural Connectivity via Persistent Homology

abstract: The brain structural connectome is generated by a collection of white matter fiber bundles constructed from diffusion weighted MRI (dMRI), acting as highways for neural activity. There has been abundant interest in studying how the structural connectome varies across individuals in relation to their traits, ranging from age and gender to neuropsychiatric outcomes. After applying tractography to dMRI to get white matter fiber bundles, a key question is how to represent the brain connectome to facilitate statistical analyses relating connectomes to traits. The current standard divides the brain into regions of interest (ROIs), and then relies on an adjacency matrix (AM) representation. Each cell in the AM is a measure of connectivity, e.g., number of fiber curves, between a pair of ROIs. Although the AM representation is intuitive, a disadvantage is the high-dimensionality due to the large number of cells in the matrix. This article proposes a simpler tree representation of the brain connectome, which is motivated by ideas in computational topology and takes topological and biological information on the cortical surface into consideration. We demonstrate that our tree representation preserves useful information and interpretability, while reducing dimensionality to improve statistical and computational efficiency. Applications to data from the Human Connectome Project (HCP) are considered and code is provided for reproducing our analyses.


title: Simplified Prophet Inequalities for Combinatorial Auctions

abstract: We consider prophet inequalities for XOS and MPH-$k$ combinatorial auctions and give a simplified proof for the existence of static and anonymous item prices which recover the state-of-the-art competitive ratios. Our proofs make use of a linear programming formulation which has a non-negative objective value if there are prices which admit a given competitive ratio $\alpha \geq 1$. Changing our perspective to dual space by an application of strong LP duality, we use an interpretation of the dual variables as probabilities to directly obtain our result. In contrast to previous work, our proofs do not require to argue about specific values of buyers for bundles, but only about the presence or absence of items. As a side remark, for any $k \geq 2$, this simplification also leads to a tiny improvement in the best competitive ratio for MPH-$k$ combinatorial auctions from $4k-2$ to $2k + 2 \sqrt{k(k-1)} -1$.


title: Inferring school district learning modalities during the COVID-19 pandemic with a hidden Markov model

abstract: In this study, learning modalities offered by public schools across the United States were investigated to track changes in the proportion of schools offering fully in-person, hybrid and fully remote learning over time. Learning modalities from 14,688 unique school districts from September 2020 to June 2021 were reported by Burbio, MCH Strategic Data, the American Enterprise Institute's Return to Learn Tracker and individual state dashboards. A model was needed to combine and deconflict these data to provide a more complete description of modalities nationwide. A hidden Markov model (HMM) was used to infer the most likely learning modality for each district on a weekly basis. This method yielded higher spatiotemporal coverage than any individual data source and higher agreement with three of the four data sources than any other single source. The model output revealed that the percentage of districts offering fully in-person learning rose from 40.3% in September 2020 to 54.7% in June of 2021 with increases across 45 states and in both urban and rural districts. This type of probabilistic model can serve as a tool for fusion of incomplete and contradictory data sources in support of public health surveillance and research efforts.


abstract: Event Detection, which aims to identify and classify mentions of event instances from unstructured articles, is an important task in Natural Language Processing (NLP). Existing techniques for event detection only use homogeneous one-hot vectors to represent the event type classes, ignoring the fact that the semantic meaning of the types is important to the task. Such an approach is inefficient and prone to overfitting. In this paper, we propose a Semantic Pivoting Model for Effective Event Detection (SPEED), which explicitly incorporates prior information during training and captures semantically meaningful correlations between input and events. Experimental results show that our proposed model achieves state-of-the-art performance and outperforms the baselines in multiple settings without using any external resources.


title: Universal location of Yang-Lee edge singularity in classic O(N) universality classes

abstract: Employing the functional renormalization group approach at next-to-leading order of the derivative expansion, we refine our earlier findings for the location of the Yang-Lee edge singularity in classic O(N) universality classes. For the universality classes of interest to QCD, in three dimensions, we found $|z_c|/R_\chi^{1/\gamma} = 1.612(9),\ 1.597(3)$ for $N=2$, $4$ correspondingly. We also established $|z_c| = 2.04(8),\ 1.69(3)$ for $N=2$, $4$ albeit with greater systematic error.


title: Alternative polynomial-time algorithm for Bipartite Matching

abstract: If $G$ is a bipartite graph, Hall's theorem \cite{H35} gives a condition for the existence of a matching of $G$ covering one side of the bipartition. This theorem admits a well-known algorithmic proof involving the repeated search of augmenting paths. We present here an alternative algorithm, using a game-theoretic formulation of the problem. We also show how to extend this formulation to the setting of balanced hypergraphs.


title: Detectability of Varied Hybridization Scenarios using Genome-Scale Hybrid Detection Methods

abstract: Hybridization events complicate the accurate reconstruction of phylogenies, as they lead to patterns of genetic heritability that are unexpected under traditional, bifurcating models of species trees. This has led to the development of methods to infer these varied hybridization events, both methods that reconstruct networks directly, and summary methods that predict individual hybridization events. However, a lack of empirical comparisons between methods - especially pertaining to large networks with varied hybridization scenarios - hinders their practical use. Here, we provide a comprehensive review of popular summary methods: TICR, MSCquartets, HyDe, Patterson's D-Statistic (ABBA-BABA), D3, and Dp. TICR and MSCquartets are based on quartet concordance factors gathered from gene tree topologies and Patterson's D-Statistic, D3, and Dp use site pattern frequencies to identify hybridization events. We then use simulated data to address questions of method accuracy and ideal use scenarios by testing methods against complex networks which depict gene flow events that differ in depth (timing), quantity (single vs. multiple, overlapping hybridizations), and rate of gene flow. We find that deeper or multiple hybridization events may introduce noise and weaken the signal of hybridization, leading to higher false negative rates across methods. Despite some forms of hybridization eluding quartet-based detection methods, MSCquartets displays high precision in most scenarios. While HyDe results in high false negative rates when tested on hybridizations involving ghost lineages, HyDe is the only method to be able to separate hybrid vs parent signals. Lastly, we test the methods on ultraconserved elements from the bee subfamily Nomiinae, finding the possibility of hybridization events between clades which correspond to regions of poor support in the species tree estimated in the original study.


title: MAgNET: A Graph U-Net Architecture for Mesh-Based Simulations

abstract: In many cutting-edge applications, high-fidelity computational models prove to be too slow for practical use and are therefore replaced by much faster surrogate models. Recently, deep learning techniques have increasingly been utilized to accelerate such predictions. To enable learning on large-dimensional and complex data, specific neural network architectures have been developed, including convolutional and graph neural networks. In this work, we present a novel encoder-decoder geometric deep learning framework called MAgNET, which extends the well-known convolutional neural networks to accommodate arbitrary graph-structured data. MAgNET consists of innovative Multichannel Aggregation (MAg) layers and graph pooling/unpooling layers, forming a graph U-Net architecture that is analogous to convolutional U-Nets. We demonstrate the predictive capabilities of MAgNET in surrogate modeling for non-linear finite element simulations in the mechanics of solids.


title: Construction of Precision sMDT Detector for ATLAS Muon Spectrometer Upgrade

abstract: This paper describes the small-diameter monitored drift-tube detector construction at the University of Michigan as a contribution to the ATLAS Muon Spectrometer upgrade for the high-luminosity Large Hadron Collider at CERN. Measurements of the first 30 chambers built at Michigan show that the drift tube wire position accuracy meets the specification of 20 microns. The positions of the platforms for alignment and magnetic field sensors are all installed well within the required precision. The cosmic ray test measurements show single wire tracking resolution of 100 +- 7 microns with an average detection efficiency above 99 %. The infrastructure, tooling, techniques, and procedures for chamber production are described in detail. The results from the chamber quality control tests of the first 30 constructed chambers are reported.


title: Tunable Dynamic Walking via Soft Twisted Beam Vibration

abstract: We propose a novel mechanism that propagates vibration through soft twisted beams, taking advantage of dynamically-coupled anisotropic stiffness to simplify the actuation of walking robots. Using dynamic simulation and experimental approaches, we show that the coupled stiffness of twisted beams with terrain contact can be controlled to generate a variety of complex trajectories by changing the frequency of the input signal. This work reveals how ground contact influences the system's dynamic behavior, supporting the design of walking robots inspired by this phenomenon. We also show that the proposed twisted beam produces a tunable walking gait from a single vibrational input.


title: Optimal Conservative Offline RL with General Function Approximation via Augmented Lagrangian

abstract: Offline reinforcement learning (RL), which refers to decision-making from a previously-collected dataset of interactions, has received significant attention over the past years. Much effort has focused on improving offline RL practicality by addressing the prevalent issue of partial data coverage through various forms of conservative policy learning. While the majority of algorithms do not have finite-sample guarantees, several provable conservative offline RL algorithms are designed and analyzed within the single-policy concentrability framework that handles partial coverage. Yet, in the nonlinear function approximation setting where confidence intervals are difficult to obtain, existing provable algorithms suffer from computational intractability, prohibitively strong assumptions, and suboptimal statistical rates. In this paper, we leverage the marginalized importance sampling (MIS) formulation of RL and present the first set of offline RL algorithms that are statistically optimal and practical under general function approximation and single-policy concentrability, bypassing the need for uncertainty quantification. We identify that the key to successfully solving the sample-based approximation of the MIS problem is ensuring that certain occupancy validity constraints are nearly satisfied. We enforce these constraints by a novel application of the augmented Lagrangian method and prove the following result: with the MIS formulation, augmented Lagrangian is enough for statistically optimal offline RL. In stark contrast to prior algorithms that induce additional conservatism through methods such as behavior regularization, our approach provably eliminates this need and reinterprets regularizers as "enforcers of occupancy validity" than "promoters of conservatism."

categories: cs.LG cs.AI math.OC math.ST stat.ML stat.TH


title: Using Unused: Non-Invasive Dynamic FaaS Infrastructure with HPC-Whisk

abstract: Modern HPC workload managers and their careful tuning contribute to the high utilization of HPC clusters. However, due to inevitable uncertainty it is impossible to completely avoid node idleness. Although such idle slots are usually too short for any HPC job, they are too long to ignore them. Function-as-a-Service (FaaS) paradigm promisingly fills this gap, and can be a good match, as typical FaaS functions last seconds, not hours. Here we show how to build a FaaS infrastructure on idle nodes in an HPC cluster in such a way that it does not affect the performance of the HPC jobs significantly. We dynamically adapt to a changing set of idle physical machines, by integrating open-source software Slurm and OpenWhisk. We designed and implemented a prototype solution that allowed us to cover up to 90\% of the idle time slots on a 50k-core cluster that runs production workloads.


title: SleepyWheels: An Ensemble Model for Drowsiness Detection leading to Accident Prevention

abstract: Around 40 percent of accidents related to driving on highways in India occur due to the driver falling asleep behind the steering wheel. Several types of research are ongoing to detect driver drowsiness but they suffer from the complexity and cost of the models. In this paper, SleepyWheels a revolutionary method that uses a lightweight neural network in conjunction with facial landmark identification is proposed to identify driver fatigue in real time. SleepyWheels is successful in a wide range of test scenarios, including the lack of facial characteristics while covering the eye or mouth, the drivers varying skin tones, camera placements, and observational angles. It can work well when emulated to real time systems. SleepyWheels utilized EfficientNetV2 and a facial landmark detector for identifying drowsiness detection. The model is trained on a specially created dataset on driver sleepiness and it achieves an accuracy of 97 percent. The model is lightweight hence it can be further deployed as a mobile application for various platforms.


title: A finite-dimensional approximation for partial differential equations on Wasserstein space

abstract: This paper presents a finite-dimensional approximation for a class of partial differential equations on the space of probability measures. These equations are satisfied in the sense of viscosity solutions. The main result states the convergence of the viscosity solutions of the finite-dimensional PDE to the viscosity solutions of the PDE on Wasserstein space, provided that uniqueness holds for the latter, and heavily relies on an adaptation of the Barles & Souganidis monotone scheme to our context, as well as on a key precompactness result for semimartingale measures. We illustrate this result with the example of the Hamilton-Jacobi-Bellman and Bellman-Isaacs equations arising in stochastic control and differential games, and propose an extension to the case of path-dependent PDEs.


title: Apple Silicon Performance in Scientific Computing

abstract: With the release of the Apple Silicon System-on-a-Chip processors, and the impressive performance shown in general use by both the M1 and M1 Ultra, the potential use for Apple Silicon processors in scientific computing is explored. Both the M1 and M1 Ultra are compared to current state-of-the-art data-center GPUs, including an NVIDIA V100 with PCIe, an NVIDIA V100 with NVLink, and an NVIDIA A100 with PCIe. The scientific performance is measured using the Scalable Heterogeneous Computing (SHOC) benchmark suite using OpenCL benchmarks. We find that both M1 processors outperform the GPUs in all benchmarks.


title: Moonraker -- Enceladus Multiple Flyby Mission

abstract: Enceladus, an icy moon of Saturn, possesses an internal water ocean and jets expelling ocean material into space. Cassini investigations indicated that the subsurface ocean could be a habitable environment having a complex interaction with the rocky core. Further investigation of the composition of the plume formed by the jets is necessary to fully understand the ocean, its potential habitability, and what it tells us about Enceladus' origin. Moonraker has been proposed as an ESA M-class mission designed to orbit Saturn and perform multiple flybys of Enceladus, focusing on traversals of the plume. The proposed Moonraker mission consists of an ESA-provided platform, with strong heritage from JUICE and Mars Sample Return, and carrying a suite of instruments dedicated to plume and surface analysis. The nominal Moonraker mission has a duration of 13.5 years. It includes a 23-flyby segment with 189 days allocated for the science phase, and can be expanded with additional segments if resources allow. The mission concept consists in investigating: i) the habitability conditions of present-day Enceladus and its internal ocean, ii) the mechanisms at play for the communication between the internal ocean and the surface of the South Polar Terrain, and iii) the formation conditions of the moon. Moonraker, thanks to state-of-the-art instruments representing a significant improvement over Cassini's payload, would quantify the abundance of key species in the plume, isotopic ratios, and physical parameters of the plume and the surface. Such a mission would pave the way for a possible future landed mission.


title: VIINTER: View Interpolation with Implicit Neural Representations of Images

abstract: We present VIINTER, a method for view interpolation by interpolating the implicit neural representation (INR) of the captured images. We leverage the learned code vector associated with each image and interpolate between these codes to achieve viewpoint transitions. We propose several techniques that significantly enhance the interpolation quality. VIINTER signifies a new way to achieve view interpolation without constructing 3D structure, estimating camera poses, or computing pixel correspondence. We validate the effectiveness of VIINTER on several multi-view scenes with different types of camera layout and scene composition. As the development of INR of images (as opposed to surface or volume) has centered around tasks like image fitting and super-resolution, with VIINTER, we show its capability for view interpolation and offer a promising outlook on using INR for image manipulation tasks.


title: ${\rm S{\scriptsize IM}BIG}$: A Forward Modeling Approach To Analyzing Galaxy Clustering

abstract: We present the first-ever cosmological constraints from a simulation-based inference (SBI) analysis of galaxy clustering from the new ${\rm S{\scriptsize IM}BIG}$ forward modeling framework. ${\rm S{\scriptsize IM}BIG}$ leverages the predictive power of high-fidelity simulations and provides an inference framework that can extract cosmological information on small non-linear scales, inaccessible with standard analyses. In this work, we apply ${\rm S{\scriptsize IM}BIG}$ to the BOSS CMASS galaxy sample and analyze the power spectrum, $P_\ell(k)$, to $k_{\rm max}=0.5\,h/{\rm Mpc}$. We construct 20,000 simulated galaxy samples using our forward model, which is based on high-resolution ${\rm Q{\scriptsize UIJOTE}}$ $N$-body simulations and includes detailed survey realism for a more complete treatment of observational systematics. We then conduct SBI by training normalizing flows using the simulated samples and infer the posterior distribution of $\Lambda$CDM cosmological parameters: $\Omega_m, \Omega_b, h, n_s, \sigma_8$. We derive significant constraints on $\Omega_m$ and $\sigma_8$, which are consistent with previous works. Our constraints on $\sigma_8$ are $27\%$ more precise than standard analyses. This improvement is equivalent to the statistical gain expected from analyzing a galaxy sample that is $\sim60\%$ larger than CMASS with standard methods. It results from additional cosmological information on non-linear scales beyond the limit of current analytic models, $k > 0.25\,h/{\rm Mpc}$. While we focus on $P_\ell$ in this work for validation and comparison to the literature, ${\rm S{\scriptsize IM}BIG}$ provides a framework for analyzing galaxy clustering using any summary statistic. We expect further improvements on cosmological constraints from subsequent ${\rm S{\scriptsize IM}BIG}$ analyses of summary statistics beyond $P_\ell$.


title: Privacy Induces Robustness: Information-Computation Gaps and Sparse Mean Estimation

abstract: We establish a simple connection between robust and differentially-private algorithms: private mechanisms which perform well with very high probability are automatically robust in the sense that they retain accuracy even if a constant fraction of the samples they receive are adversarially corrupted. Since optimal mechanisms typically achieve these high success probabilities, our results imply that optimal private mechanisms for many basic statistics problems are robust. We investigate the consequences of this observation for both algorithms and computational complexity across different statistical problems. Assuming the Brennan-Bresler secret-leakage planted clique conjecture, we demonstrate a fundamental tradeoff between computational efficiency, privacy leakage, and success probability for sparse mean estimation. Private algorithms which match this tradeoff are not yet known -- we achieve that (up to polylogarithmic factors) in a polynomially-large range of parameters via the Sum-of-Squares method. To establish an information-computation gap for private sparse mean estimation, we also design new (exponential-time) mechanisms using fewer samples than efficient algorithms must use. Finally, we give evidence for privacy-induced information-computation gaps for several other statistics and learning problems, including PAC learning parity functions and estimation of the mean of a multivariate Gaussian.


title: LARO: Learned Acquisition and Reconstruction Optimization to accelerate Quantitative Susceptibility Mapping

abstract: Quantitative susceptibility mapping (QSM) involves acquisition and reconstruction of a series of images at multi-echo time points to estimate tissue field, which prolongs scan time and requires specific reconstruction technique. In this paper, we present our new framework, called Learned Acquisition and Reconstruction Optimization (LARO), which aims to accelerate the multi-echo gradient echo (mGRE) pulse sequence for QSM. Our approach involves optimizing a Cartesian multi-echo k-space sampling pattern with a deep reconstruction network. Next, this optimized sampling pattern was implemented in an mGRE sequence using Cartesian fan-beam k-space segmenting and ordering for prospective scans. Furthermore, we propose to insert a recurrent temporal feature fusion module into the reconstruction network to capture signal redundancies along echo time. Our ablation studies show that both the optimized sampling pattern and proposed reconstruction strategy help improve the quality of the multi-echo image reconstructions. Generalization experiments show that LARO is robust on the test data with new pathologies and different sequence parameters. Our code is available at https://github.com/Jinwei1209/LARO.git.


title: Asymmetric transport for magnetic Dirac equations

abstract: This paper concerns the asymmetric transport associated with a low-energy interface Dirac model of graphene-type materials subject to external magnetic and electric fields. We show that the relevant physical observable, an interface conductivity, is quantized and robust to a large class of perturbations. These include defects that decay along or away from the interface, and sufficiently small or localized changes in the external fields. An explicit formula for the interface conductivity is given by a spectral flow.


title: Quantum Natural Language Generation on Near-Term Devices

abstract: The emergence of noisy medium-scale quantum devices has led to proof-of-concept applications for quantum computing in various domains. Examples include Natural Language Processing (NLP) where sentence classification experiments have been carried out, as well as procedural generation, where tasks such as geopolitical map creation, and image manipulation have been performed. We explore applications at the intersection of these two areas by designing a hybrid quantum-classical algorithm for sentence generation. Our algorithm is based on the well-known simulated annealing technique for combinatorial optimisation. An implementation is provided and used to demonstrate successful sentence generation on both simulated and real quantum hardware. A variant of our algorithm can also be used for music generation. This paper aims to be self-contained, introducing all the necessary background on NLP and quantum computing along the way.


title: Genuine multifractality in time series is due to temporal correlations

abstract: Based on the mathematical arguments formulated within the Multifractal Detrended Fluctuation Analysis (MFDFA) approach it is shown that in the uncorrelated time series from the Gaussian basin of attraction the effects resembling multifractality asymptotically disappear for positive moments when the length of time series increases. A hint is given that this applies to the negative moments as well and extends to the L\'evy stable regime of fluctuations. The related effects are also illustrated and confirmed by numerical simulations. This documents that the genuine multifractality in time series may only result from the long-range temporal correlations and the fatter distribution tails of fluctuations may broaden the width of singularity spectrum only when such correlations are present. The frequently asked question of what makes multifractality in time series - temporal correlations or broad distribution tails - is thus ill posed. In the absence of correlations only the bifractal or monofractal cases are possible. The former corresponds to the L\'evy stable regime of fluctuations while the latter to the ones belonging to the Gaussian basin of attraction in the sense of the Central Limit Theorem.


title: A Bayesian Framework on Asymmetric Mixture of Factor Analyser

abstract: Mixture of factor analyzer (MFA) model is an efficient model for the analysis of high dimensional data through which the factor-analyzer technique based on the covariance matrices reducing the number of free parameters. The model also provides an important methodology to determine latent groups in data. There are several pieces of research to extend the model based on the asymmetrical and/or with outlier datasets with some known computational limitations that have been examined in frequentist cases. In this paper, an MFA model with a rich and flexible class of skew normal (unrestricted) generalized hyperbolic (called SUNGH) distributions along with a Bayesian structure with several computational benefits have been introduced. The SUNGH family provides considerable flexibility to model skewness in different directions as well as allowing for heavy tailed data. There are several desirable properties in the structure of the SUNGH family, including, an analytically flexible density which leads to easing up the computation applied for the estimation of parameters. Considering factor analysis models, the SUNGH family also allows for skewness and heavy tails for both the error component and factor scores. In the present study, the advantages of using this family of distributions have been discussed and the suitable efficiency of the introduced MFA model using real data examples and simulation has been demonstrated.