title: Dynamic State Estimation for Radial Microgrid Protection

abstract: Microgrids are localized electrical grids with control capability that are able to disconnect from the traditional grid to operate autonomously. They strengthen grid resilience, help mitigate grid disturbances, and support a flexible grid by enabling the integration of distributed energy resources. Given the likely presence of critical loads, the proper protection of microgrids is of vital importance; however, this is complicated in the case of inverter-interfaced microgrids where low fault currents preclude the use of conventional time-overcurrent protection. This paper introduces and investigates the application of dynamic state estimation, a generalization of differential protection, for the protection of radial portions of microgrids (or distribution networks); both phasor-based and dynamic approaches are investigated for protection. It is demonstrated through experiments on three case-study systems that dynamic state estimation is capable of correctly identifying model parameters for both normal and faulted operation.


title: Multiwavelength constraints on the unidentified Galactic TeV sources HESS J1427$-$608, HESS J1458$-$608, and new VHE $\gamma$-ray source candidates

abstract: The H.E.S.S. Galactic Plane Survey (HGPS) revealed 78 TeV sources among which 47 are not clearly associated with a known object. We present a multiwavelength approach to constrain the origin of the emission from unidentified HGPS sources. We present a generic pipeline that explores a large database of multiwavelength archival data toward any region in the Galactic plane. Along with a visual inspection of the retrieved multiwavelength observations to search for faint and uncataloged counterparts, we derive a radio spectral index that helps disentangle thermal from nonthermal emission and a mean magnetic field through X-ray and TeV data in case of a leptonic scenario. We also search for a spectral connection between the GeV and the TeV regimes with the Fermi-LAT cataloged sources that may be associated with the unidentified HGPS source. We complete the association procedure with catalogs of known objects and with the source catalogs from instruments whose data are retrieved. The method is applied on two unidentified sources, namely HESS J1427$-$608 and HESS J1458$-$608, for which the multiwavelength constraints favor the pulsar wind nebula (PWN) scenario. We model their broadband nonthermal spectra in a leptonic scenario with a magnetic field $B \lesssim 10$ $\mu$G, which is consistent with that obtained from ancient PWNe. We place both sources within the context of the TeV PWN population to estimate the spin-down power and the characteristic age of the putative pulsar. We also shed light on two possibly significant $\gamma$-ray excesses in the HGPS: the first is located in the south of the unidentified source HESS J1632$-$478 and the second is spatially coincident with the synchrotron-emitting supernova remnant G28.6$-$0.1. The multiwavelength counterparts found toward both $\gamma$-ray excesses make these promising candidates for being new very-high energy $\gamma$-ray sources.


title: Testing Simultaneous Diagonalizability

abstract: This paper proposes novel methods to test for simultaneous diagonalization of possibly asymmetric matrices. Motivated by various applications, a two-sample test as well as a generalization for multiple matrices are proposed. A partial version of the test is also studied to check whether a partial set of eigenvectors is shared across samples. Additionally, a novel algorithm for the considered testing methods is introduced. Simulation studies demonstrate favorable performance for all designs. Finally, the theoretical results are utilized to decouple vector autoregression models into multiple univariate time series, and to test for the same stationary distribution in recurrent Markov chains. These applications are demonstrated using macroeconomic indices of 8 countries and streamflow data, respectively.


title: The Grothendieck Construction in Categorical Network Theory

abstract: In this thesis, we present a flexible framework for specifying and constructing operads which are suited to reasoning about network construction. The data used to present these operads is called a \emph{network model}, a monoidal variant of Joyal's combinatorial species. The construction of the operad required that we develop a monoidal lift of the Grothendieck construction. We then demonstrate how concepts like priority and dependency can be represented in this framework. For the former, we generalize Green's graph products of groups to the context of universal algebra. For the latter, we examine the emergence of monoidal fibrations from the presence of catalysts in Petri nets.


title: The geometry and DSZ quantization of four-dimensional supergravity

abstract: We implement the Dirac-Schwinger-Zwanziger integrality condition on four-dimensional classical ungauged supergravity and use it to obtain its duality-covariant, gauge-theoretic, differential-geometric model on an oriented four-manifold $M$ of arbitrary topology. Classical bosonic supergravity is completely determined by a submersion $\pi$ over $M$ equipped with a complete Ehresmann connection, a vertical euclidean metric, and a vertically-polarized flat symplectic vector bundle $\Xi$. Building on these structures, we implement the Dirac-Schwinger-Zwanziger integrality condition through the choice of an element in the degree-two sheaf cohomology group with coefficients in a locally constant sheaf $\mathcal{L}\subset \Xi$ valued in the groupoid of integral symplectic spaces. We show that this data determines a Siegel principal bundle $P_{\mathfrak{t}}$ of fixed type $\mathfrak{t}\in \mathbb{Z}^{n_v}$ whose connections provide the global geometric description of the local electromagnetic gauge potentials of the theory. Furthermore, we prove that the Maxwell gauge equations of the theory reduce to the polarized self-duality condition determined by $\Xi$ on the connections of $P_{\mathfrak{t}}$. In addition, we investigate the continuous and discrete U-duality groups of the theory, characterizing them through short exact sequences and realizing the latter through the gauge group of $P_{\mathfrak{t}}$ acting on its adjoint bundle. This elucidates the geometric origin of U-duality, which we explore in several examples, illustrating its dependence on the topology of the fiber bundles $\pi$ and $P_{\mathfrak{t}}$ as well as on the isomorphism type of $\mathcal{L}$.


title: Collaborative Experience between Scientific Software Projects using Agile Scrum Development

abstract: Developing sustainable software for the scientific community requires expertise in software engineering and domain science. This can be challenging due to the unique needs of scientific software, the insufficient resources for software engineering practices in the scientific community, and the complexity of developing for evolving scientific contexts. While open-source software can partially address these concerns, it can introduce complicating dependencies and delay development. These issues can be reduced if scientists and software developers collaborate. We present a case study wherein scientists from the SuperNova Early Warning System collaborated with software developers from the Scalable Cyberinfrastructure for Multi-Messenger Astrophysics project. The collaboration addressed the difficulties of open-source software development, but presented additional risks to each team. For the scientists, there was a concern of relying on external systems and lacking control in the development process. For the developers, there was a risk in supporting a user-group while maintaining core development. These issues were mitigated by creating a second Agile Scrum framework in parallel with the developers' ongoing Agile Scrum process. This Agile collaboration promoted communication, ensured that the scientists had an active role in development, and allowed the developers to evaluate and implement the scientists' software requirements. The collaboration provided benefits for each group: the scientists actuated their development by using an existing platform, and the developers utilized the scientists' use-case to improve their systems. This case study suggests that scientists and software developers can avoid scientific computing issues by collaborating and that Agile Scrum methods can address emergent concerns.


title: Multiscale prediction of microstructure length scales in metallic alloy casting

abstract: In this article, we combine casting experiments and quantitative simulations to present a novel multiscale modeling approach to predict local primary dendritic spacings in metallic alloys solidified in conditions relevant to industrial casting processes. To this end, primary dendritic spacings were measured in instrumented casting experiments in Al-Cu alloys containing 1\,wt\% and 4\,wt\% of Cu, and they were compared to spacing stability ranges and average spacings in dendritic arrays simulated using phase-field (PF) and dendritic needle network (DNN) models. It is first shown that PF and DNN lead to similar results for the Al-1\,wt\%Cu alloy, using a dendrite tip selection constant calculated with PF in the DNN simulations. PF simulations cannot achieve quantitative predictions for the Al-4\,wt\%Cu alloy because they are too computationally demanding due to the large separation of scale between tip radius and diffusion length, a characteristic feature of non-dilute alloys. Nevertheless, the results of DNN simulations for non-dilute Al-Cu alloys are in overall good agreement with our experimental results as well as with those of an extensive literature review. Simulations consistently suggest a widening of the PDAS stability range with a decrease of the temperature gradient as the microstructure goes from cellular-dendrites to well-developed hierarchical dendrites.


title: Minimax Off-Policy Evaluation for Multi-Armed Bandits

abstract: We study the problem of off-policy evaluation in the multi-armed bandit model with bounded rewards, and develop minimax rate-optimal procedures under three settings. First, when the behavior policy is known, we show that the Switch estimator, a method that alternates between the plug-in and importance sampling estimators, is minimax rate-optimal for all sample sizes. Second, when the behavior policy is unknown, we analyze performance in terms of the competitive ratio, thereby revealing a fundamental gap between the settings of known and unknown behavior policies. When the behavior policy is unknown, any estimator must have mean-squared error larger -- relative to the oracle estimator equipped with the knowledge of the behavior policy -- by a multiplicative factor proportional to the support size of the target policy. Moreover, we demonstrate that the plug-in approach achieves this worst-case competitive ratio up to a logarithmic factor. Third, we initiate the study of the partial knowledge setting in which it is assumed that the minimum probability taken by the behavior policy is known. We show that the plug-in estimator is optimal for relatively large values of the minimum probability, but is sub-optimal when the minimum probability is low. In order to remedy this gap, we propose a new estimator based on approximation by Chebyshev polynomials that provably achieves the optimal estimation error. Numerical experiments on both simulated and real data corroborate our theoretical findings.


title: A nonabelian Brunn-Minkowski inequality

abstract: Henstock and Macbeath asked in 1953 whether the Brunn-Minkowski inequality can be generalized to nonabelian locally compact groups; questions along the same line were also asked by Hrushovski, McCrudden, and Tao. We obtain here such an inequality and prove that it is sharp for helix-free locally compact groups, which includes real linear algebraic groups, Nash groups, semisimple Lie groups with finite center, solvable Lie groups, etc. The proof follows an induction on dimension strategy; new ingredients include an understanding of the role played by maximal compact subgroups of Lie groups, a necessary modified form of the inequality which is also applicable to nonunimodular locally compact groups, and a proportionated averaging trick.


title: Implementing Admittance Relaying for Microgrid Protection

abstract: The rapid increase of distributed energy resources has led to the widespread deployment of microgrids. These flexible and efficient local energy grids are able to operate in both grid-connected mode and islanded mode; they are interfaced to the main power system by a fast semiconductor switch and commonly make use of inverter-interfaced generation. This paper focuses on inverter interfaced microgrids, which present a challenge for protection as they do not provide the high short-circuit current necessary for conventional time-overcurrent protection. The application of admittance relaying for the protection of inverter-interfaced microgrids is investigated as a potential solution. The comparison of analytical and simulated results of performed four experiments prove the suitability of admittance relaying for microgrids protection.


title: Parametrically driven Kerr cavity solitons

abstract: Temporal cavity solitons are optical pulses that propagate indefinitely in nonlinear resonators. They are currently attracting a lot of attention, both for their many potential applications and for their connection to other fields of science. Cavity solitons are phase locked to a driving laser. This is what distinguishes them from laser dissipative solitons and the main reason why they are excellent candidates for precision applications such as optical atomic clocks. To date, the focus has been on driving Kerr solitons close to their carrier frequency, in which case a single stable localised solution exists for fixed parameters. Here we experimentally demonstrate, for the first time, Kerr cavity solitons excitation around twice their carrier frequency. In that configuration, called parametric driving, two solitons of opposite phase may coexist. We use a fibre resonator that incorporates a quadratically nonlinear section and excite stable solitons by scanning the driving frequency. Our experimental results are in excellent agreement with a seminal amplitude equation, highlighting connections to hydrodynamic and mechanical systems, amongst others. Furthermore, we experimentally confirm that two different phase-locked solitons may be simultaneously excited and harness this multiplicity to generate a string of random bits, thereby extending the pool of applications of Kerr resonators to random number generators and Ising machines.


title: Hamiltonian Control of Magnetic Field Lines: Computer Assisted Results Proving the Existence of KAM Barriers

abstract: We reconsider a control theory for Hamiltonian systems, that was introduced on the basis of KAM theory and applied to a model of magnetic field in previous articles. By a combination of Frequency Analysis and of a rigorous (Computer Assisted) KAM algorithm we prove that in the phase space of the magnetic field, due to the control term, a set of invariant tori appear, and it acts as a transport barrier. Our analysis, which is common (but often also limited) to Celestial Mechanics, is based on a normal form approach; it is also quite general and can be applied to quasi-integrable Hamiltonian systems satisfying a few additional mild assumptions. As a novelty with respect to the works that in the last two decades applied Computer Assisted Proofs into the framework of KAM theory, we provide all the codes allowing to produce our results. They are collected in a software package that is publicly available from the {\it Mendeley Data} repository. All these codes are designed in such a way to be easy-to-use, also for what concerns eventual adaptations for applications to similar problems.


title: Deterministic photonic quantum computation in a synthetic time dimension

abstract: Photonics offers unique advantages as a substrate for quantum information processing, but imposes fundamental scalability challenges. Nondeterministic schemes impose massive resource overheads, while deterministic schemes require prohibitively many identical quantum emitters to realize sizeable quantum circuits. Here we propose a scalable architecture for a photonic quantum computer which needs minimal quantum resources to implement any quantum circuit: a single coherently controlled atom. Optical switches endow a photonic quantum state with a synthetic time dimension by modulating photon-atom couplings. Quantum operations applied to the atomic qubit can be teleported onto the photonic qubits via projective measurement, and arbitrary quantum circuits can be compiled into a sequence of these teleported operators. This design negates the need for many identical quantum emitters to be integrated into a photonic circuit and allows effective all-to-all connectivity between photonic qubits. The proposed device has a machine size which is independent of quantum circuit depth, does not require single-photon detectors, operates deterministically, and is robust to experimental imperfections.


title: Contrast-enhanced spectral mammography with a photon-counting detector

abstract: Purpose: Spectral imaging is a method in medical x-ray imaging to extract information about the object constituents by the material-specific energy dependence of x-ray attenuation. The authors have investigated a photon-counting spectral imaging system with two energy bins for contrast-enhanced mammography. System optimization and the potential benefit compared to conventional non-energy-resolved absorption imaging was studied. Methods: A framework for system characterization was set up that included quantum and anatomical noise and a theoretical model of the system was benchmarked to phantom measurements. Results: Optimal combination of the energy-resolved images corresponded approximately to minimization of the anatomical noise, which is commonly referred to as energy subtraction. In that case, an ideal-observer detectability index could be improved close to 50% compared to absorption imaging in the phantom study. Optimization with respect to the signal-to-quantum-noise ratio, commonly referred to as energy weighting, yielded only a minute improvement. In a simulation of a clinically more realistic case, spectral imaging was predicted to perform approximately 30% better than absorption imaging for an average glandularity breast with an average level of anatomical noise. For dense breast tissue and a high level of anatomical noise, however, a rise in detectability by a factor of 6 was predicted. Another ~70-90% improvement was found to be within reach for an optimized system. Conclusions: Contrast-enhanced spectral mammography is feasible and beneficial with the current system, and there is room for additional improvements. Inclusion of anatomical noise is essential for optimizing spectral imaging systems.


title: An efficient pre-object collimator based on an x-ray lens

abstract: A multi-prism lens (MPL) is a refractive x-ray lens with one-dimensional focusing properties. If used as a pre-object collimator in a scanning system for medical x-ray imaging, it reduces the divergence of the radiation and improves on photon economy compared to a slit collimator. Potential advantages include shorter acquisition times, a reduced tube loading, or improved resolution. We present the first images acquired with an MPL in a prototype for a scanning mammography system. The lens showed a gain of flux of 1.32 compared to a slit collimator at equal resolution, or a gain in resolution of 1.31-1.44 at equal flux. We expect the gain of flux in a clinical set-up with an optimized MPL and a custom-made absorption filter to reach 1.67, or 1.45-1.54 gain in resolution.


title: Energy resolution of a photon-counting silicon strip detector

abstract: A photon-counting silicon strip detector with two energy thresholds was investigated for spectral X-ray imaging in a mammography system. Preliminary studies already indicate clinical benefit of the detector, and the purpose of the present study is optimization with respect to energy resolution. Factors relevant for the energy response were measured, simulated, or gathered from previous studies, and used as input parameters to a cascaded detector model. Threshold scans over several X-ray spectra were used to calibrate threshold levels to energy, and to validate the model. The energy resolution of the detector assembly was assessed to range over DeltaE/E = 0.12-0.26 in the mammography region. Electronic noise dominated the peak broadening, followed by charge sharing between adjacent detector strips, and a channel-to-channel threshold spread. The energy resolution may be improved substantially if these effects are reduced to a minimum. Anti-coincidence logic mitigated double counting from charge sharing, but erased the energy resolution of all detected events, and optimization of the logic is desirable. Pile-up was found to be of minor importance at typical mammography rates.


title: The variational method, backreactions, and the absorption probability in Wald type problems

abstract: We argue that the variational method in Wald type thought experiments, involves order of magnitude problems when one imposes the fact that $\delta M$ is inherently a first order quantity itself. One observes that the contribution of the second order perturbations is actually of the fourth order. Therefore backreactions have to be explicitly calculated. Here, we re-consider the overspinning problem for Kerr-Newman black holes interacting with test fields. We calculate the backreaction effects due to the induced increase in the angular velocity of the event horizon, which brings a partial solution to the overspinning problem. To bring an ultimate solution, we argue that the absorption probability should be taken into account in Wald type problems where black holes interact with test fields. This fundamentally alters the course of the analysis of the thought experiments. Due to the fact that a small fraction of the challenging modes is absorbed by the black holes, overspinning is prevented for both nearly extremal and extremal cases. Some extreme cases are easily fixed by backreaction effects. The arguments do not apply to the generic overspinning by fermionic fields for which the absorption probability is positive definite.


title: $R$-process enhancements of Gaia-Enceladus in GALAH DR3

abstract: The dominant site of production of $r$-process elements remains unclear despite recent observations of a neutron star merger. Observational constraints on the properties of the sites can be obtained by comparing $r$-process abundances in different environments. The recent Gaia data releases and large samples from high-resolution optical spectroscopic surveys are enabling us to compare $r$-process element abundances between stars formed in an accreted dwarf galaxy, Gaia-Enceladus, and those formed in the Milky Way. We aim to understand the origin of $r$-process elements in Gaia-Enceladus. We first construct a sample of stars to study Eu abundances without being affected by the detection limit. We then kinematically select 71 Gaia-Enceladus stars and 93 in-situ stars from the Galactic Archaeology with HERMES (GALAH) DR3, of which 50 and 75 stars can be used to study Eu reliably. Gaia-Enceladus stars clearly show higher ratios of [{Eu}/{Mg}] than in-situ stars. High [{Eu}/{Mg}] along with low [{Mg}/{Fe}] are also seen in relatively massive satellite galaxies such as the LMC, Fornax, and Sagittarius dwarfs. On the other hand, unlike these galaxies, Gaia-Enceladus does not show enhanced [{Ba}/{Eu}] or [{La}/{Eu}] ratios suggesting a lack of significant $s$-process contribution. From comparisons with simple chemical evolution models, we show that the high [{Eu}/{Mg}] of Gaia-Enceladus can naturally be explained by considering $r$-process enrichment by neutron-star mergers with delay time distribution that follows a similar power-law as type~Ia supernovae but with a shorter minimum delay time.


title: Rare Earth Ions Doped Mixed Crystals for Fast Quantum Computers with Optical Frequency Qubits

abstract: The possibility of using mixed crystals highly doped with rare earth ions (REIs) as physical systems for creating fast quantum computers with a sampling time of nanoseconds is discussed. The electronic 4f states of rare earth ions with small values of the diagonal elements of the Judd-Ofelt matrix U(2) are proposed as optical frequency qubit levels. CNOT and other conditional gate operations are performed by exciting the rare earth ion into the 4f state with a large diagonal element of U(2), causing a Stark blockade. It is found that the main interaction responsible for this blockade is the quadrupole-quadrupole interaction. The large inhomogeneous broadening of the frequencies of the electronic transitions in mixed crystals and the weak interaction of 4f electrons with phonons make it possible to achieve a high computation rate and a long decoherence time of the qubits. An ensemble of closest REIs is described that can act as an OQC instance; the frequencies of the corresponding qubits can be found using the spectral hole burning method.


title: The Observed Rate of Binary Black Hole Mergers can be Entirely Explained by Globular Clusters

abstract: Since the first signal in 2015, the gravitational-wave detections of merging binary black holes (BBHs) by the LIGO and Virgo collaborations (LVC) have completely transformed our understanding of the lives and deaths of compact object binaries, and have motivated an enormous amount of theoretical work on the astrophysical origin of these objects. We show that the phenomenological fit to the redshift-dependent merger rate of BBHs from Abbott et al. (2020) is consistent with a purely dynamical origin for these objects, and that the current merger rate of BBHs from the LVC could be explained entirely with globular clusters alone. While this does not prove that globular clusters are the dominant formation channel, we emphasize that many formation scenarios could contribute a significant fraction of the current LVC rate, and that any analysis that assumes a single (or dominant) mechanism for producing BBH mergers is implicitly using a specious astrophysical prior.


title: On Monte-Carlo methods in convex stochastic optimization

abstract: We develop a novel procedure for estimating the optimizer of general convex stochastic optimization problems of the form $\min_{x\in\mathcal{X}} \mathbb{E}[F(x,\xi)]$, when the given data is a finite independent sample selected according to $\xi$. The procedure is based on a median-of-means tournament, and is the first procedure that exhibits the optimal statistical performance in heavy tailed situations: we recover the asymptotic rates dictated by the central limit theorem in a non-asymptotic manner once the sample size exceeds some explicitly computable threshold. Additionally, our results apply in the high-dimensional setup, as the threshold sample size exhibits the optimal dependence on the dimension (up to a logarithmic factor). The general setting allows us to recover recent results on multivariate mean estimation and linear regression in heavy-tailed situations and to prove the first sharp, non-asymptotic results for the portfolio optimization problem.

categories: math.ST math.OC math.PR q-fin.MF stat.TH


title: On the derivation of the Khmaladze transforms

abstract: Some 40 years ago Khmaladze introduced a transform which greatly facilitated the distribution free goodness of fit testing of statistical hypotheses. In the last decade, he has published a related transform, broadly offering an alternative means to the same end. The aim of this paper is to derive these transforms using relatively elementary means, making some simplifications, but losing little in the way of generality. In this way it is hoped to make these transforms more accessible and more widely used in statistical practice. We also propose a change of name of the second transform to the Khmaladze rotation, in order to better reflect its nature.


title: A note on an open conjecture in rational dynamical system

abstract: This note is an attempt with the open conjecture $8$ proposed by the authors of\cite{G. LADAS} which states: Assume $\alpha,\beta, \lambda \in [0,\infty)$. Then every positive solution of the difference equation : \begin{equation*} z_{n+1}=\frac{\alpha+z_{n}\beta +z_{n-1}\lambda}{z_{n-2}},\quad n=0,1,\ldots \end{equation*} is bounded if and only if $\beta=\lambda$. We will use a construction of sub-energy function and properties of Todd's difference equation to disprove that conjecture in general.


title: Multiple giant eruptions and X-ray emission in the recoiling AGN/LBV candidate SDSS1133

abstract: We present a comprehensive analysis of 20 years worth of multi-color photometric light curves, multi-epoch optical spectra, and X-ray data of an off-nuclear variable object SDSS1133 in Mrk 177 at $z=0.0079$. The UV-optical light curves reveal that SDSS1133 experienced four outbursts in 2001, 2014, 2019, and 2021. The persistent UV-optical luminosity in the non-outbursting state is $\sim 10^{41}$ erg/s with small-scale flux variations, and peak luminosities during the outbursts reach $\sim 10^{42}$ erg/s. The optical spectra exhibit enduring broad hydrogen Balmer P-Cygni profiles with the absorption minimum at $\sim -2,000$ km/s, indicating the presence of fast moving ejecta. {\it Chandra} detected weak X-ray emission at a $0.3-10$ keV luminosity of $L_{X} = 4 \times 10^{38}$ erg/s after the 2019 outburst. These lines of evidence suggests that SDSS1133 is an extreme luminous blue variable (LBV) star experiencing multiple giant eruptions with interactions of the ejected shell with different shells and/or circumstellar medium (CSM), and disfavors the recoiling Active Galactic Nuclei (AGN) scenario suggested in the literature. We suggest that pulsational pair-instability may provide a viable explanation for the multiple energetic eruptions in SDSS1133. If the current activity of SDSS1133 is a precursor of a supernova explosion, we may be able to observe a few additional giant eruptions and then the terminal supernova explosion or collapse to a massive black hole in future observations.


title: Connection among environment, cloud-cloud collision speed, and star formation activity in the strongly barred galaxy NGC1300

abstract: Cloud-cloud collision (CCC) has been suggested as a mechanism to induce massive star formation. Recent simulations suggest that a CCC speed is different among galactic-scale environments, which is responsible for observed differences in star formation activity. In particular, a high-speed CCC is proposed as a cause of star formation suppression in the bar regions in barred spiral galaxies. Focusing on the strongly barred galaxy NGC1300, we investigate the CCC speed. We find the CCC speed in the bar and bar-end tend to be higher than that in the arm. The estimated CCC speed is $\sim20~\rm km~s^{-1}$, $\sim16~\rm km~s^{-1}$, and $\sim11~\rm km~s^{-1}$ in the bar, bar-end, and arm, respectively. Although the star formation activity is different in the bar and bar-end, the CCC speed and the number density of high-speed CCC with $> 20~\rm km~s^{-1}$ are high in both regions, implying the existence of other parameters that control the star formation. The difference in molecular gas mass (average density) of the giant molecular clouds (GMCs) between the bar (lower mass and lower density) and bar-end (higher mass and higher density) may be cause for the different star formation activity. Combining with our previous study (Maeda et al.), the leading candidates of causes for the star formation suppression in the bar in NGC1300 are the presence of a large amount of diffuse molecular gases and high-speed CCCs between low mass GMCs.


title: Low Surface Brightness Galaxies in z > 1 Galaxy Clusters: HST approaches the Progenitors of Local Ultra Diffuse Galaxies

abstract: Ultra Diffuse Galaxies (UDGs), a type of large Low Surface Brightness (LSB) galaxies with particularly large effective radii (r_eff > 1.5 kpc), are now routinely studied in the local (z<0.1) universe. While they are found to be abundant in clusters, groups, and in the field, their formation mechanisms remain elusive and an active topic of debate. New insights may be found by studying their counterparts at higher redshifts (z>1.0), even though cosmological surface brightness dimming makes them particularly diffcult to detect and study there. This work uses the deepest Hubble Space Telescope (HST) imaging stacks of z > 1 clusters, namely: SPT-CL J2106-5844 and MOO J1014+0038. These two clusters, at z=1.13 and z=1.23, were monitored as part of the HST See-Change program. Compared to the Hubble Extreme Deep Field (XDF) as reference field, we find statistical over-densities of large LSB galaxies in both clusters. Based on stellar population modelling and assuming no size evolution, we find that the faintest sources we can detect are about as bright as expected for the progenitors of the brightest local UDGs. We find that the LSBs we detect in SPT-CL J2106-5844 and MOO J1014-5844 already have old stellar populations that place them on the red sequence. Correcting for incompleteness, and based on an extrapolation of local scaling relations, we estimate that distant UDGs are relatively under-abundant compared to local UDGs by a factor ~ 3. A plausible explanation for the implied increase with time would be a significant size growth of these galaxies in the last ~ 8 Gyr, as also suggested by hydrodynamical simulations.


title: How baryons can significantly bias cluster count cosmology

abstract: We quantify two main pathways through which baryonic physics biases cluster count cosmology. We create mock cluster samples that reproduce the baryon content inferred from X-ray observations. We link clusters to their counterparts in a dark matter-only universe, whose abundances can be predicted robustly, by assuming the dark matter density profile is not significantly affected by baryons. We derive weak lensing halo masses and infer the best-fitting cosmological parameters $\Omega_\mathrm{m}$, $S_8=\sigma_8(\Omega_\mathrm{m}/0.3)^{0.2}$, and $w_0$ from the mock cluster sample. We find that because of the need to accommodate the change in the density profile due to the ejection of baryons, weak lensing mass calibrations are only unbiased if the concentration is left free when fitting the reduced shear with NFW profiles. However, even unbiased total mass estimates give rise to biased cosmological parameters if the measured mass functions are compared with predictions from dark matter-only simulations. This bias dominates for haloes with $m_\mathrm{500c} < 10^{14.5} \, h^{-1} \, \mathrm{M_\odot}$. For a stage IV-like cluster survey without mass estimation uncertainties, an area $\approx 15000 \, \mathrm{deg^2}$ and a constant mass cut of $m_\mathrm{200m,min} = 10^{14} \, h^{-1} \, \mathrm{M_\odot}$, the biases are $-11 \pm 1 \, \%$ in $\Omega_\mathrm{m}$, $-3.29 \pm 0.04 \, \%$ in $S_8$, and $9 \pm 1.5 \, \%$ in $w_0$. The statistical significance of the baryonic bias depends on how accurately the actual uncertainty on individual cluster mass estimates is known. We suggest that rather than the total halo mass, the (re-scaled) dark matter mass inferred from the combination of weak lensing and observations of the hot gas, should be used for cluster count cosmology.


title: Sub-arcsecond Imaging of the Complex Organic Chemistry in Massive Star-forming Region G10.6-0.4

abstract: Massive star-forming regions exhibit an extremely rich and diverse chemistry, which in principle provides a wealth of molecular probes, as well as laboratories for interstellar prebiotic chemistry. Since the chemical structure of these sources displays substantial spatial variation among species on small scales (${\lesssim}10^4$ au), high angular resolution observations are needed to connect chemical structures to local environments and inform astrochemical models of massive star formation. To address this, we present ALMA 1.3 mm observations toward OB cluster-forming region G10.6-0.4 (hereafter "G10.6") at a resolution of 0.14$^{\prime\prime}$ (700 au). We find highly-structured emission from complex organic molecules (COMs) throughout the central 20,000 au, including two hot molecular cores and several shells or filaments. We present spatially-resolved rotational temperature and column density maps for a large sample of COMs and warm gas tracers. These maps reveal a range of gas substructure in both O- and N-bearing species. We identify several spatial correlations that can be explained by existing models of COM formation, including NH$_2$CHO/HNCO and CH$_3$OCHO/CH$_3$OCH$_3$, but also observe unexpected distributions and correlations which suggest that our current understanding of COM formation is far from complete. Importantly, complex chemistry is observed throughout G10.6, rather than being confined to hot cores. The COM composition appears to be different in the cores compared to the more extended structures, which illustrates the importance of high spatial resolution observations of molecular gas in elucidating the physical and chemical processes associated with massive star formation.


title: Know the enemy: 2D Fermi liquids

abstract: We describe an analytical theory investigating the regime of validity of the Fermi liquid theory in interacting, via the long-range Coulomb coupling, two-dimensional Fermi systems comparing it with with the corresponding 3D systems. We find that the 2D Fermi liquid theory and 2D quasiparticles are robust up to high energies and temperatures of the order of Fermi energy above the Fermi surface, very similar to the corresponding three-dimensional situation. We calculate the phase diagram in the frequency-temperature space separating the collisionless ballistic regime and the collision-dominated hydrodynamic regime for 2D and 3D interacting electron systems. We also provide the temperature corrections up to third order for the renormalized effective mass, and comment on the validity of 2D Wiedemann-Franz law and 2D Kadawoki-Woods relation.


title: Charming ALPs

abstract: Axion-like particles (ALPs) are ubiquitous in models of new physics explaining some of the most pressing puzzles of the Standard Model. However, until relatively recently, little attention has been paid to its interplay with flavour. In this work, we study in detail the phenomenology of ALPs that exclusively interact with up-type quarks at the tree-level, which arise in some well-motivated ultra-violet completions such as QCD-like dark sectors or Froggatt-Nielsen type models of flavour. Our study is performed in the low-energy effective theory to highlight the key features of these scenarios in a model independent way. We derive all the existing constraints on these models and demonstrate how upcoming experiments at fixed-target facilities and the LHC can probe a vast region of the parameter space, which is currently not excluded by cosmological and astrophysical bounds. We also emphasize how a future measurement of the currently unavailable meson decay $D \to \pi + \rm{invisible}$ could complement these upcoming searches and help to probe a large unexplored region of their parameter space.


title: Disorder-induced topology in quench dynamics

abstract: We study the effect of strong disorder on topology and entanglement in quench dynamics. Although disorder-induced topological phases have been well studied in equilibrium, the disorder-induced topology in quench dynamics has not been explored. In this work, we predict a disorder-induced topology of post-quench states characterized by the quantized dynamical Chern number and the crossings in the entanglement spectrum in $(1+1)$ dimensions. The dynamical Chern number undergoes transitions from zero to unity, and back to zero when increasing the disorder strength. The boundaries between different dynamical Chern numbers are determined by delocalized critical points in the post-quench Hamiltonian with the strong disorder. An experimental realization in quantum walks is discussed.


title: Stiefel Liquids: Possible Non-Lagrangian Quantum Criticality from Intertwined Orders

abstract: We propose a new type of quantum liquids, dubbed Stiefel liquids, based on $2+1$ dimensional nonlinear sigma models on target space $SO(N)/SO(4)$, supplemented with Wess-Zumino-Witten terms. We argue that the Stiefel liquids form a class of critical quantum liquids with extraordinary properties, such as large emergent symmetries, a cascade structure, and nontrivial quantum anomalies. We show that the well known deconfined quantum critical point and $U(1)$ Dirac spin liquid are unified as two special examples of Stiefel liquids, with $N=5$ and $N=6$, respectively. Furthermore, we conjecture that Stiefel liquids with $N>6$ are non-Lagrangian, in the sense that under renormalization group they flow to infrared (conformally invariant) fixed points that cannot be described by any renormalizable continuum Lagrangian. Such non-Lagrangian states are beyond the paradigm of parton gauge mean-field theory familiar in the study of exotic quantum liquids in condensed matter physics. The intrinsic absence of (conventional or parton-like) mean-field construction also means that, within the traditional approaches, it will be difficult to decide whether a non-Lagrangian state can actually emerge from a specific UV system (such as a lattice spin system). For this purpose we hypothesize that a quantum state is emergible from a lattice system if its quantum anomalies match with the constraints from the (generalized) Lieb-Schultz-Mattis theorems. Based on this hypothesis, we find that some of the non-Lagrangian Stiefel liquids can indeed be realized in frustrated quantum spin systems, for example, on triangular or Kagome lattice, through the intertwinement between non-coplanar magnetic orders and valence-bond-solid orders.

categories: cond-mat.str-el cond-mat.quant-gas hep-th math-ph math.MP


title: Non-Halo Structures and their Effects on Gravitational Lensing

abstract: Anomalies in the flux-ratios of the images of quadruply-lensed quasars have been used to constrain the nature of dark matter. Assuming these lensing perturbations are caused by dark matter haloes, it is currently possible to constrain the mass of a hypothetical Warm Dark Matter (WDM) particle to be $m_\chi > 5.2$ keV. However, the assumption that perturbations are only caused by DM haloes might not be correct as other structures, such as filaments and pancakes, exist and make up a significant fraction of the mass in the universe, ranging between 5$\%$ -- 50$\%$ depending on the dark matter model. Using novel fragmentation-free simulations of 1 and 3keV WDM cosmologies we study these "non-halo" structures and estimate their impact on flux-ratio observations. We find that these structures display sharp density gradients with short correlation lengths, and can contribute more to the lensing signal than all haloes up to the half-mode mass combined, thus reducing the differences expected among WDM models. We estimate that non-halo structures can be the dominant cause of line-of-sight flux-ratio anomalies in very warm, but already excluded, $m_x \sim 1 \rm{keV}$ scenarios. For colder cases $m_x \gtrsim 3 \rm{keV}$, we estimate that non-haloes can contribute about $5 - 10\%$ of the total flux-ratio signal.


title: Continuum approach to real time dynamics of 1+1D gauge field theory: out of horizon correlations of the Schwinger model

abstract: We develop a truncated Hamiltonian method to study nonequilibrium real time dynamics in the Schwinger model - the quantum electrodynamics in D=1+1. This is a purely continuum method that captures reliably the invariance under local and global gauge transformations and does not require a discretisation of space-time. We use it to study a phenomenon that is expected not to be tractable using lattice methods: we show that the 1+1D quantum electrodynamics admits the dynamical horizon violation effect which was recently discovered in the case of the sine-Gordon model. Following a quench of the model, oscillatory long-range correlations develop, manifestly violating the horizon bound. We find that the oscillation frequencies of the out-of-horizon correlations correspond to twice the masses of the mesons of the model suggesting that the effect is mediated through correlated meson pairs. We also report on the cluster violation in the massive version of the model, previously known in the massless Schwinger model. The results presented here reveal a novel nonequilibrium phenomenon in 1+1D quantum electrodynamics and make a first step towards establishing that the horizon violation effect is present in gauge field theory.

categories: hep-th cond-mat.stat-mech gr-qc quant-ph


title: Randomizing multi-product formulas for Hamiltonian simulation

abstract: Quantum simulation, the simulation of quantum processes on quantum computers, suggests a path forward for the efficient simulation of problems in condensed-matter physics, quantum chemistry, and materials science. While the majority of quantum simulation algorithms are deterministic, a recent surge of ideas has shown that randomization can greatly benefit algorithmic performance. In this work, we introduce a scheme for quantum simulation that unites the advantages of randomized compiling on the one hand and higher-order multi-product formulas, as they are used for example in linear-combination-of-unitaries (LCU) algorithms or quantum error mitigation, on the other hand. In doing so, we propose a framework of randomized sampling that is expected to be useful for programmable quantum simulators and present two new multi-product formula algorithms tailored to it. Our framework reduces the circuit depth by circumventing the need for oblivious amplitude amplification required by the implementation of multi-product formulas using standard LCU methods, rendering it especially useful for early quantum computers used to estimate the dynamics of quantum systems instead of performing full-fledged quantum phase estimation. Our algorithms achieve a simulation error that shrinks exponentially with the circuit depth. To corroborate their functioning, we prove rigorous performance bounds as well as the concentration of the randomized sampling procedure. We demonstrate the functioning of the approach for several physically meaningful examples of Hamiltonians, including fermionic systems and the Sachdev-Ye-Kitaev model, for which the method provides a favorable scaling in the effort.


title: The Gas Content and Stripping of Local Group Dwarf Galaxies

abstract: The gas content of the complete compilation of Local Group dwarf galaxies (119 within 2 Mpc) is presented using HI survey data. Within the virial radius of the Milky Way (224 kpc here), 53 of 55 dwarf galaxies are devoid of gas to limits of M$_{\rm HI}<10^4$ M$_\odot$. Within the virial radius of M31 (266 kpc), 27 of 30 dwarf galaxies are devoid of gas (with limits typically $<10^5$ M$_\odot$). Beyond the virial radii of the Milky Way and M31, the majority of the dwarf galaxies have detected HI gas and have HI masses higher than the limits. When the relationship between gas content and distance is investigated using a Local Group virial radius, more of the non-detected dwarf galaxies are within this radius (85$\pm1$ of the 93 non-detected dwarf galaxies) than within the virial radii of the Milky Way and M31. Using the Gaia proper motion measurements available for 38 dwarf galaxies, the minimum gas density required to completely strip them of gas is calculated. Halo densities between $10^{-5}$ and $5 \times 10^{-4}$ cm$^{-3}$ are typically required for instantaneous stripping at perigalacticon. When compared to halo density with radius expectations from simulations and observations, 80% of the dwarf galaxies with proper motions are consistent with being stripped by ram pressure at Milky Way pericenter. The results suggest a diffuse gaseous galactic halo medium is important in quenching dwarf galaxies, and that a Local Group medium also potentially plays a role.


title: Dark Matter Constraints from a Unified Analysis of Strong Gravitational Lenses and Milky Way Satellite Galaxies

abstract: Joint analyses of small-scale cosmological structure probes are relatively unexplored and promise to advance measurements of microphysical dark matter properties using heterogeneous data. Here, we present a multidimensional analysis of dark matter substructure using strong gravitational lenses and the Milky Way (MW) satellite galaxy population, accounting for degeneracies in model predictions and using covariances in the constraining power of these individual probes for the first time. We simultaneously infer the projected subhalo number density and the half-mode mass describing the suppression of the subhalo mass function in thermal relic warm dark matter (WDM), $M_{\mathrm{hm}}$, using the semianalytic model $\mathrm{\texttt{Galacticus}}$ to connect the subhalo population inferred from MW satellite observations to the strong lensing host halo mass and redshift regime. Combining MW satellite and strong lensing posteriors in this parameter space yields $M_{\mathrm{hm}}<10^{7.0}\ M_{\mathrm{\odot}}$ (WDM particle mass $m_{\mathrm{WDM}}>9.7\ \mathrm{keV}$) at $95\%$ confidence and disfavors $M_{\mathrm{hm}}=10^{7.4}\ M_{\mathrm{\odot}}$ ($m_{\mathrm{WDM}}=7.4\ \mathrm{keV}$) with a 20:1 marginal likelihood ratio, improving limits on $m_{\mathrm{WDM}}$ set by the two methods independently by $\sim 30\%$. These results are marginalized over the line-of-sight contribution to the strong lensing signal, the mass of the MW host halo, and the efficiency of subhalo disruption due to baryons and are robust to differences in the disruption efficiency between the MW and strong lensing regimes at the $\sim 10\%$ level. This work paves the way for unified analyses of next-generation small-scale structure measurements covering a wide range of scales and redshifts.


title: Combined Constraints on First Generation Leptoquarks

abstract: In this article we perform a combined analysis of low energy precision constraints and LHC searches for leptoquarks which couple to first generation fermions. Considering all ten leptoquark representations, five scalar and five vector ones, we study at the precision frontier the constraints from $K\to\pi\nu\nu$, $K\to\pi e^+e^-$, $K^0-\bar K^0$ and $D^0-\bar D^0$ mixing, as well as from experiments searching for parity violation (APV and QWEAK). We include LHC searches for $s$-channel single resonant production, pair production and Drell-Yan-like signatures of leptoquarks. Interestingly, we find that the recent non-resonant di-lepton analysis of ATLAS provides stronger bounds than the resonant searches recasted so far to constrain $t$-channel production of leptoquarks. Taking into account all these bounds, we observe that none of the leptoquark representations can address the so-called "Cabibbo angle anomaly" via a direct contribution to super-allowed beta decays.


title: Primordial Black Holes in the Excursion Set Theory

abstract: We study primordial black holes (PBHs) formation in the excursion set theory (EST) in a vast range of PBHs masses with and without confirmed constraints on their abundance. In this work, we introduce a new concept of the first touch in the context of EST for PBHs formation. This new framework takes into account the earlier horizon reentry of smaller masses. Our study shows that in the EST, it is possible to produce PBHs in different mass range, with enhanced power spectrum, which could make up all dark matter. We also show that in a broad blue-tilted power spectrum, the production of PBHs is dominated by smaller masses. Our analysis put an upper limit $\sim\,$0.1 on the amplitude of the curvature power spectrum at length scales relevant for PBHs formation.


title: Solving Quadratic Unconstrained Binary Optimization with divide-and-conquer and quantum algorithms

abstract: Quadratic Unconstrained Binary Optimization (QUBO) is a broad class of optimization problems with many practical applications. To solve its hard instances in an exact way, known classical algorithms require exponential time and several approximate methods have been devised to reduce such cost. With the growing maturity of quantum computing, quantum algorithms have been proposed to speed up the solution by using either quantum annealers or universal quantum computers. Here we apply the divide-and-conquer approach to reduce the original problem to a collection of smaller problems whose solutions can be assembled to form a single Polynomial Binary Unconstrained Optimization instance with fewer variables. This technique can be applied to any QUBO instance and leads to either an all-classical or a hybrid quantum-classical approach. When quantum heuristics like the Quantum Approximate Optimization Algorithm (QAOA) are used, our proposal leads to a double advantage: a substantial reduction of quantum resources, specifically an average of ~42% fewer qubits to solve MaxCut on random 3-regular graphs, together with an improvement in the quality of the approximate solutions reached.


title: Symmetry-resolved dynamical purification in synthetic quantum matter

abstract: When a quantum system initialized in a product state is subjected to either coherent or incoherent dynamics, the entropy of any of its connected partitions generically increases as a function of time, signalling the inevitable spreading of (quantum) information throughout the system. Here, we show that, in the presence of continuous symmetries and under ubiquitous experimental conditions, symmetry-resolved information spreading is inhibited due to the competition of coherent and incoherent dynamics: in given quantum number sectors, entropy decreases as a function of time, signalling dynamical purification. Such dynamical purification bridges between two distinct short and intermediate time regimes, characterized by a log-volume and log-area entropy law, respectively. It is generic to symmetric quantum evolution, and as such occurs for different partition geometry and topology, and classes of (local) Liouville dynamics. We then develop a protocol to measure symmetry-resolved entropies and negativities in synthetic quantum systems based on the random unitary toolbox, and demonstrate the generality of dynamical purification using experimental data from trapped ion experiments [Brydges et al., Science 364, 260 (2019)]. Our work shows that symmetry plays a key role as a magnifying glass to characterize many-body dynamics in open quantum systems, and, in particular, in noisy-intermediate scale quantum devices.


title: Presence versus absence of Two-Dimensional Fermi Surface Anomalies

abstract: We theoretically consider Fermi surface anomalies manifesting in the temperature dependent quasiparticle properties of two-dimensional (2D) interacting electron systems, comparing and contrasting with the corresponding 3D Fermi liquid situation. In particular, employing microscopic many body perturbative techniques, we obtain analytically the leading-order and the next-to-leading-order interaction corrections to the renormalized effective mass for three distinct physical interaction models: electron-phonon, electron-paramagnon, and electron-electron Coulomb coupling. We find that the 2D renormalized effective mass does not develop any Fermi surface anomaly due to electron-phonon interaction, manifesting $\mathcal{O}(T^2)$ temperature correction and thus remaining consistent with the Sommerfeld expansion of the non-interacting Fermi function, in contrast to the corresponding 3D situation where the temperature correction to the renormalized effective mass has the anomalous $T^2 \log T$ behavior. By contrast, both electron-paramagnon and electron-electron interactions lead to the anomalous $\mathcal{O}(T)$ corrections to the 2D effective mass renormalization in contrast to $T^2 \log T$ behavior in the corresponding 3D interacting systems. We provide detailed analytical results, and comment on the conditions under which a $T^2 \log T$ term could possibly arise in the 2D quasiparticle effective mass from electron-phonon interactions. We also compare results for the temperature dependent specific heat in the interacting 2D and 3D Fermi systems, using the close connection between the effective mass and specific heat.


title: Internet of Predictable Things (IoPT) Framework to Increase Cyber-Physical System Resiliency

abstract: During the last two decades, distributed energy systems, especially renewable energy sources (RES), have become more economically viable with increasing market share and penetration levels on power systems. In addition to decarbonization and decentralization of energy systems, digitalization has also become very important. The use of artificial intelligence (AI), advanced optimization algorithms, Industrial Internet of Things (IIoT), and other digitalization frameworks makes modern power system assets more intelligent, while vulnerable to cybersecurity risks. This paper proposes the concept of the Internet of Predictable Things (IoPT) that incorporates advanced data analytics and machine learning methods to increase the resiliency of cyber-physical systems against cybersecurity risks. The proposed concept is demonstrated using a cyber-physical system testbed under a variety of cyber attack scenarios as a proof of concept (PoC).


title: Merger-Ringdown Consistency: A New Test of Strong Gravity using Deep Learning

abstract: The gravitational waves emitted during the coalescence of binary black holes are an excellent probe to test the behaviour of strong gravity. In this paper, we propose a new test called the `merger-ringdown consistency test` that focuses on probing the imprints of the dynamics in strong-gravity around the black-holes during the plunge-merger and ringdown phase. Furthermore, we present a scheme that allows us to efficiently combine information across multiple ringdown observations to perform a statistical null test of GR using the detected BH population. We present a proof-of-concept study for this test using simulated binary black hole ringdowns embedded in the next-generation ground-based detector noise. We demonstrate the feasibility of our test using a deep learning framework, setting a precedence for performing precision tests of gravity with neural networks.


title: Simultaneous supply and demand constraints in input-output networks: The case of Covid-19 in Germany, Italy, and Spain

abstract: Natural and anthropogenic disasters frequently affect both the supply and demand side of an economy. A striking recent example is the Covid-19 pandemic which has created severe disruptions to economic output in most countries. These direct shocks to supply and demand will propagate downstream and upstream through production networks. Given the exogenous shocks, we derive a lower bound on total shock propagation. We find that even in this best case scenario network effects substantially amplify the initial shocks. To obtain more realistic model predictions, we study the propagation of shocks bottom-up by imposing different rationing rules on industries if they are not able to satisfy incoming demand. Our results show that economic impacts depend strongly on the emergence of input bottlenecks, making the rationing assumption a key variable in predicting adverse economic impacts. We further establish that the magnitude of initial shocks and network density heavily influence model predictions.


title: Tangent infinity-categories and Goodwillie calculus

abstract: We make precise the analogy between Goodwillie's calculus of functors in homotopy theory and the differential calculus of smooth manifolds by introducing a higher-categorical framework of which both theories are examples. That framework is an extension to infinity-categories of the tangent categories of Cockett and Cruttwell (introduced originally by Rosick\'y). The basic data of a tangent infinity-category consist of an endofunctor, that plays the role of the tangent bundle construction, together with various natural transformations that mimic structure possessed by the ordinary tangent bundles of smooth manifolds. The role of the tangent bundle functor in Goodwillie calculus is played by Lurie's tangent bundle for infinity-categories, introduced to generalize the cotangent complexes of Andr\'e, Quillen and Illusie. We show that Lurie's construction admits the additional structure maps and satisfies the conditions needed to form a tangent infinity-category which we refer to as the Goodwillie tangent structure. Cockett and Cruttwell (and others) have started to develop various aspects of differential geometry in the abstract context of tangent categories, and we begin to apply those ideas to Goodwillie calculus. For example, we show that the role of Euclidean spaces in the calculus of manifolds is played in Goodwillie calculus by the stable infinity-categories. We also show that Goodwillie's n-excisive functors are the direct analogues of n-jets of smooth maps between manifolds; to state that connection precisely, we develop a notion of tangent (infinity,2)-category and show that Goodwillie calculus is best understood in that context.


title: Policy choices can help keep 4G and 5G universal broadband affordable

abstract: The United Nations Broadband Commission has committed the international community to accelerate universal broadband. However, the cost of meeting this objective, and the feasibility of doing so on a commercially viable basis, are not well understood. Using scenario analysis, this paper compares the global cost-effectiveness of different infrastructure strategies for the developing world to achieve universal 4G or 5G mobile broadband. Utilizing remote sensing and demand forecasting, least-cost network designs are developed for eight representative low and middle-income countries (Malawi, Uganda, Kenya, Senegal, Pakistan, Albania, Peru and Mexico), the results from which form the basis for aggregation to the global level. The cost of meeting a minimum 10 Mbps per user is estimated at USD 1.7 trillion using 5G Non-Standalone, approximately 0.6% of annual GDP for the developing world over the next decade. However, by creating a favorable regulatory environment, governments can bring down these costs by as much as three quarters, to USD 0.5 trillion (approximately 0.2% of annual GDP), and avoid the need for public subsidy. Providing governments make judicious choices, adopting fiscal and regulatory regimes conducive to lowering costs, universal broadband may be within reach of most developing countries over the next decade.


title: Sub-damped Lyman alpha systems in the XQ-100 survey II -- Chemical evolution at 2.4<z<4.3

abstract: We present the measured gas-phase metal column densities in 155 sub-damped Lyman alpha systems (subDLAs) with the aim to investigate the contribution of subDLAs to the chemical evolution of the Universe. The sample was identified within the absorber-blind XQ-100 quasar spectroscopic survey over the redshift range 2.4<=z<=4.3. Using all available column densities of the ionic species investigated (mainly CIV, SiII, MgII, SiIV, AlII, FeII, CII, and OI; in order of decreasing detection frequency), we estimate the ionization-corrected gas-phase metallicity of each system using Markov Chain Monte Carlo techniques to explore a large grid of Cloudy ionization models. Without accounting for ionization and dust depletion effects, we find that the HI-weighted gas-phase metallicity evolution of subDLAs are consistent with damped Lyman alpha systems (DLAs). When ionization corrections are included, subDLAs are systematically more metal-poor than DLAs (between ~0.5 sigma and ~3 sigma significance) by up to ~1.0 dex over the redshift range 3<=z<=4.3. The correlation of gas-phase [Si/Fe] with metallicity in subDLAs appears to be consistent with that of DLAs, suggesting that the two classes of absorbers have a similar relative dust depletion pattern. As previously seen for Lyman limit systems, the gas-phase [C/O] in subDLAs remains constantly solar for all metallicities indicating that both subDLAs and Lyman limit systems could trace carbon-rich ejecta, potentially in circumgalactic environments.


title: A Census of the Extended Neutral Hydrogen Around 18 MHONGOOSE Galaxies

abstract: We present the analysis of the diffuse, low column density HI environment of 18 MHONGOOSE galaxies. We obtained deep observations with the Robert C. Byrd Green Bank Telescope, and reached down to a 3sigma column density detection limit of NHI=6.3x10^{17} cm^{-2} over a 20 km/s linewidth. We analyze the environment around these galaxies, with a focus on HI gas that reaches column densities below NHI=10^{19} cm^{-2}. We calculate the total amount of HI gas in and around the galaxies revealing that nearly all of these galaxies contained excess HI outside of their disks. We quantify the amount of diffuse gas in the maps of each galaxy, defined by HI gas with column densities below 10^{19} cm^{-2}, and find a large spread in percentages of diffuse gas. However, by binning the percentage of diffuse HI into quarters, we find that the bin with the largest number of galaxies is the lowest quartile (0-25\% diffuse HI). We identified several galaxies which may be undergoing gas accretion onto the galaxy disk using multiple methods of analysis, including azimuthally averaging column densities beyond the disk, and identifying structure within our integrated intensity (Moment 0) maps. We measured HI mass outside the disks of most of our galaxies, with rising cumulative flux even at large radii. We also find a strong correlation between the fraction of diffuse gas in a galaxy and its baryonic mass, and test this correlation using both Spearman and Pearson correlation coefficients. We see evidence of a dark matter halo mass threshold of M_{halo}~10^{11.1} \msun{} in which galaxies with high fractions of diffuse HI all reside below. It is in this regime in which cold-mode accretion should dominate. Finally, we suggest a rotation velocity of v_{rot}~80 km\s as an upper threshold to find diffuse gas-dominated galaxies.


title: The role of the nominal iron content in the structural, compositional and physical properties of BaFe$_{2+{\delta}}$S$_3$

abstract: BaFe$_2$S$_3$ is a quasi-one-dimensional antiferromagnetic insulator that becomes superconducting under hydrostatic pressure. The magnetic ordering temperature, $T_N$, as well as the presence of superconductivity have been found to be sample dependent. It has been argued that the Fe content may play a decisive role, with the use of 5%mol excess Fe being reportedly required during the synthesis to optimize the magnetic ordering temperature and the superconducting properties. However, it is yet unclear whether an Fe off-stoichiometry is actually present in the samples, and how it affects the structural, magnetic and transport properties. Here, we present a systematic study of compositional, structural and physical properties of BaFe$_{2+\delta}$S$_3$ as a function of the nominal Fe excess $\delta$. As $\delta$ increases, we observe the presence of an increasing fraction of secondary phases but no systematic change in the composition or crystal structure of the main phase. Magnetic susceptibility curves are influenced by the presence of magnetic secondary phases. The previously reported maximum of $T_N$ at $\delta$=0.1 was not confirmed. Samples with nominal $\delta$=0 present the lowest $T_N$ and the resistivity anomaly at the highest temperature $T^*$ while, for $\delta \geq 0.05$, both quantities and the transport gap are seemingly $\delta$-independent. Finally, we show that crystals free of ferromagnetic spurious phases can be obtained by remelting samples with nominal $\delta$=0.05 in a Bridgman process.


title: Machine learning applications for COVID-19: A state-of-the-art review

abstract: The COVID-19 pandemic has galvanized the machine learning community to create new solutions that can help in the fight against the virus. The body of literature related to applications of machine learning and artificial intelligence to COVID-19 is constantly growing. The goal of this article is to present the latest advances in machine learning research applied to COVID-19. We cover four major areas of research: forecasting, medical diagnostics, drug development, and contact tracing. We review and analyze the most successful state of the art studies. In contrast to other existing surveys on the subject, our article presents a high level overview of the current research that is sufficiently detailed to provide an informed insight.


title: Safe and Efficient Model-free Adaptive Control via Bayesian Optimization

abstract: Adaptive control approaches yield high-performance controllers when a precise system model or suitable parametrizations of the controller are available. Existing data-driven approaches for adaptive control mostly augment standard model-based methods with additional information about uncertainties in the dynamics or about disturbances. In this work, we propose a purely data-driven, model-free approach for adaptive control. Tuning low-level controllers based solely on system data raises concerns on the underlying algorithm safety and computational performance. Thus, our approach builds on GoOSE, an algorithm for safe and sample-efficient Bayesian optimization. We introduce several computational and algorithmic modifications in GoOSE that enable its practical use on a rotational motion system. We numerically demonstrate for several types of disturbances that our approach is sample efficient, outperforms constrained Bayesian optimization in terms of safety, and achieves the performance optima computed by grid evaluation. We further demonstrate the proposed adaptive control approach experimentally on a rotational motion system.


title: Evolution of circumstellar discs in young star-forming regions

abstract: The evolution of circumstellar discs is influenced by their surroundings. The relevant processes include external photoevaporation due to nearby stars, and dynamical truncations. The impact of these processes on disc populations depends on the star-formation history and on the dynamical evolution of the region. Since star formation history and the phase-space characteristics of the stars are important for the evolution of the discs, we start simulating the evolution of the star cluster with the results of molecular cloud collapse simulations. In the simulation we form stars with circumstellar discs, which can be affected by different processes. Our models account for the viscous evolution of the discs, internal and external photoevaporation of gas, external photoevaporation of dust, and dynamical truncations. All these processes are resolved together with the dynamical evolution of the cluster, and the evolution of the stars. An extended period of star formation, lasting for at least 2 Myr, results in some discs being formed late. These late formed discs have a better chance of survival because the cluster gradually expands with time, and a lower local stellar density reduces the effects of photoevaporation and dynamical truncation. Late formed discs can then be present in regions of high UV radiation, solving the proplyd lifetime problem. We also find a considerable fraction of discs that lose their gas content, but remain sufficiently rich in solids to be able to form a rocky planetary system.


title: Source Counts Spanning Eight Decades of Flux Density at 1.4 GHz

abstract: Brightness-weighted differential source counts $S^2 n(S)$ spanning the eight decades of flux density between $0.25\,\mu\mathrm{Jy}$ and 25 Jy at 1.4 GHz were measured from (1) the confusion brightness distribution in the MeerKAT DEEP2 image below $10\,\mu\mathrm{Jy}$, (2) counts of DEEP2 sources between $10\,\mu\mathrm{Jy}$ and $2.5\,\mathrm{mJy}$, and (3) counts of NVSS sources stronger than $2.5\,\mathrm{mJy}$. We present our DEEP2 catalog of $1.7 \times 10^4$ discrete sources complete above $S = 10\,\mu\mathrm{Jy}$ over $\Omega = 1.04\,\mathrm{deg}^2$. The brightness-weighted counts converge as $S^2 n(S) \propto S^{1/2}$ below $S = 10\,\mu\mathrm{Jy}$, so $>99\%$ of the $\Delta T_\mathrm{b} \sim 0.06\,\mathrm{K}$ sky brightness produced by active galactic nuclei and $\approx96\%$ of the $\Delta T_\mathrm{b} \sim 0.04\,\mathrm{K}$ added by star-forming galaxies has been resolved into sources with $S \geq 0.25\,\mu\mathrm{Jy}$. The $\Delta T_\mathrm{b} \approx 0.4\,\mathrm{K}$ excess brightness measured by ARCADE 2 cannot be produced by faint sources smaller than $\approx 50\,\mathrm{kpc}$ if they cluster like galaxies.


title: Factorization method for some inhomogeneous Lienard equations

abstract: We obtain closed-form solutions of several inhomogeneous Lienard equations by the factorization method. The two factorization conditions involved in the method are turned into a system of first-order differential equations containing the forcing term. In this way, one can find the forcing terms that lead to integrable cases. Because of the reduction of order feature of factorization, the solutions are simultaneously solutions of first-order differential equations with polynomial nonlinearities. The illustrative examples of Lienard solutions obtained in this way generically have rational parts, and consequently display singularities.


title: Spicules and downflows in the solar chromosphere

abstract: High-speed downflows have been observed in the solar transition region (TR) and lower corona for many decades. Despite their abundance, it has been hard to find signatures of such downflows in the solar chromosphere. In this work, we target an enhanced network region that shows ample occurrences of rapid spicular downflows in the \halpha\ spectral line that could potentially be linked to high-speed TR downflowing counterparts. We used the $k$-means algorithm to classify the spectral profiles of on-disk spicules in \halpha{} and \cak{} data observed from the Swedish 1-m Solar Telescope (SST) and employed an automated detection method based on advanced morphological image processing operations to detect such downflowing features, in conjunction with rapid blue-shifted and red-shifted excursions (RBEs and RREs). We report the existence of a new category of RREs (termed as downflowing RRE) for the first time that, contrary to earlier interpretation, are associated with chromospheric field-aligned downflows moving towards the strong magnetic field regions. Statistical analysis performed on nearly 20,000 RBEs and 15,000 RREs (including the downflowing counterparts), detected in our 97~min long dataset, shows that the downflowing RREs are very similar to RBEs and RREs except for their oppositely directed plane-of-sky motion. Furthermore, we also find that RBEs, RREs and downflowing RREs can be represented by a wide range of spectral profiles with varying Doppler offsets, and \halpha{} line core widths, both along and perpendicular to the spicule axis, that causes them to be associated with multiple substructures that evolve together. We speculate that these rapid plasma downflows could well be the chromospheric counterparts of the commonly observed TR downflows.


title: Semi-device-independent framework based on restricted distrust in prepare-and-measure experiments

abstract: A semi-device-independent framework for prepare-and-measure experiments is introduced in which an experimenter can tune the degree of distrust in the performance of the quantum devices. In this framework, a receiver operates an uncharacterised measurement device and a sender operates a preparation device that emits states with a bounded fidelity with respect to a set of target states. No assumption on Hilbert space dimension is required. The set of quantum correlations is investigated and bounded from both the interior and the exterior. Furthermore, the optimal performance of quantum state discrimination with bounded distrust is derived and applied to certification of detection efficiency. Quantum-over-classical advantages are demonstrated and the magnitude of distrust compatible with such advantages is explored. Finally, efficient schemes for semi-device-independent random number generation are developed.


title: Multi-Task Network Pruning and Embedded Optimization for Real-time Deployment in ADAS

abstract: Camera-based Deep Learning algorithms are increasingly needed for perception in Automated Driving systems. However, constraints from the automotive industry challenge the deployment of CNNs by imposing embedded systems with limited computational resources. In this paper, we propose an approach to embed a multi-task CNN network under such conditions on a commercial prototype platform, i.e. a low power System on Chip (SoC) processing four surround-view fisheye cameras at 10 FPS. The first focus is on designing an efficient and compact multi-task network architecture. Secondly, a pruning method is applied to compress the CNN, helping to reduce the runtime and memory usage by a factor of 2 without lowering the performances significantly. Finally, several embedded optimization techniques such as mixed-quantization format usage and efficient data transfers between different memory areas are proposed to ensure real-time execution and avoid bandwidth bottlenecks. The approach is evaluated on the hardware platform, considering embedded detection performances, runtime and memory bandwidth. Unlike most works from the literature that focus on classification task, we aim here to study the effect of pruning and quantization on a compact multi-task network with object detection, semantic segmentation and soiling detection tasks.


title: The Devils in the Point Clouds: Studying the Robustness of Point Cloud Convolutions

abstract: Recently, there has been a significant interest in performing convolution over irregularly sampled point clouds. Since point clouds are very different from regular raster images, it is imperative to study the generalization of the convolution networks more closely, especially their robustness under variations in scale and rotations of the input data. This paper investigates different variants of PointConv, a convolution network on point clouds, to examine their robustness to input scale and rotation changes. Of the variants we explored, two are novel and generated significant improvements. The first is replacing the multilayer perceptron based weight function with much simpler third degree polynomials, together with a Sobolev norm regularization. Secondly, for 3D datasets, we derive a novel viewpoint-invariant descriptor by utilizing 3D geometric properties as the input to PointConv, in addition to the regular 3D coordinates. We have also explored choices of activation functions, neighborhood, and subsampling methods. Experiments are conducted on the 2D MNIST & CIFAR-10 datasets as well as the 3D SemanticKITTI & ScanNet datasets. Results reveal that on 2D, using third degree polynomials greatly improves PointConv's robustness to scale changes and rotations, even surpassing traditional 2D CNNs for the MNIST dataset. On 3D datasets, the novel viewpoint-invariant descriptor significantly improves the performance as well as robustness of PointConv. We achieve the state-of-the-art semantic segmentation performance on the SemanticKITTI dataset, as well as comparable performance with the current highest framework on the ScanNet dataset among point-based approaches.


title: Implicit Bias of Linear RNNs

abstract: Contemporary wisdom based on empirical studies suggests that standard recurrent neural networks (RNNs) do not perform well on tasks requiring long-term memory. However, precise reasoning for this behavior is still unknown. This paper provides a rigorous explanation of this property in the special case of linear RNNs. Although this work is limited to linear RNNs, even these systems have traditionally been difficult to analyze due to their non-linear parameterization. Using recently-developed kernel regime analysis, our main result shows that linear RNNs learned from random initializations are functionally equivalent to a certain weighted 1D-convolutional network. Importantly, the weightings in the equivalent model cause an implicit bias to elements with smaller time lags in the convolution and hence, shorter memory. The degree of this bias depends on the variance of the transition kernel matrix at initialization and is related to the classic exploding and vanishing gradients problem. The theory is validated in both synthetic and real data experiments.


title: The isomorphism problem for tensor algebras of multivariable dynamical systems

abstract: We resolve the isomorphism problem for tensor algebras of unital multivariable dynamical systems. Specifically we show that unitary equivalence after a conjugation for multivariable dynamical systems is a complete invariant for complete isometric isomorphisms between their tensor algebras. In particular, this settles a conjecture of Davidson and Kakariadis relating to work of Arveson from the sixties, and extends related work of Kakariadis and Katsoulis.


title: An improvement of a saddle point theorem and some of its applications

abstract: In this paper, we establish an improved version of a saddle point theorem ([4]) removing a weak lower semicontinuity assumption at all. We then revisit some of the applications of that theorem in the light of such an improvement. For instance, we obtain the following very general result of local nature: Let $(H,\langle\cdot,\cdot\rangle)$ be a real Hilbert space and $\Phi:B_{\rho}\to H$ a $C^{1,1}$ function, with $\Phi(0)\neq 0$. Then, for each $r>0$ small enough, there exist only two points points $x^*, u^*\in S_r$, such that $$\max\{\langle \Phi(x^*),x^*-x\rangle, \langle \Phi(x),x^*-x\rangle\}< 0\ ,$$ for all $x\in B_r\setminus \{x^*\}$, $$\|\Phi(u^*)-u^*\|=dist(\Phi(u^*),B_r)$$ and $$\|\Phi(x)-u^*\|<\|\Phi(x)-x\|$$ for all $x\in B_r\setminus \{u^*\}$, where $$B_r=\{x\in H : \|x\|\leq r\}$$ and $$S_r=\{x\in H : \|x\|=r\}\ .$$


title: Effective field theory for closed strings near the Hagedorn temperature

abstract: We discuss interacting, closed, bosonic and superstrings in thermal equilibrium at temperatures close to the Hagedorn temperature in flat space. We calculate S-matrix elements of the strings at the Hagedorn temperature and use them to construct a low-energy effective action for interacting strings near the Hagedorn temperature. We show, in particular, that the four-point amplitude of massless winding modes leads to a positive quartic interaction. Furthermore, the effective field theory has a generalized conformal structure, namely, it is conformally invariant when the temperature is assigned an appropriate scaling dimension. Then, we show that the equations of motion resulting from the effective action possess a winding-mode-condensate background solution above the Hagedorn temperature and present a worldsheet conformal field theory, similar to a Sine-Gordon theory, that corresponds to this solution. We find that the Hagedorn phase transition in our setup is second order, in contrast to a first-order transition that was found previously in different setups.


title: Estimation of the electrostatic effects in the LISA-Pathfinder critical test mass dynamics via the method of moments

abstract: LISA-Pathfinder is an ESA space mission flown between 2015 and 2017 to demonstrate a technological maturity sufficient for building a gravitational waves telescope in space, such as the Laser Interferometer Space Antenna (LISA). A pair of cubic test masses is hosted inside the LISA-Pathfinder spacecraft and shielded from any force other than the interplanetary gravitational field. The purity of the shielding gives the performance of the mission. There are a number of aspects that had to be confirmed in-flight. One of them is the transition phase from the launch configuration, when the test masses are locked, to the science free-falling configuration. Each test mass is initially released from the mechanical constraints via a dedicated mechanism and then captured by an electrostatic control system. In fact, each test mass is surrounded by a set of electrodes for actuation and sensing purposes. The performance criterion of the release is the final velocity of the test mass relative to the spacecraft, with an upper threshold set to 5 $\mu m/s$. The LISA-Pathfinder first in-flight release velocities highlighted an unexpected dynamics with large linear and angular velocities. The electrostatic control was successful, but only relying on a manual procedure that cannot be considered as baseline for LISA. This paper helps investigating the in-flight non-compliance by dealing with the modeling of the electrostatic environment around each test mass and its contribution to the release and capture dynamics. The electrostatic model is based on the method of moments, a boundary element numerical technique suitable for estimating forces and capacitances between conductors. We also provide a short overview of the method, which can be used for the analysis of other phenomena within LISA and for the design of future gravitational waves telescopes and space projects.


title: On abelian subgroups of finite groups

abstract: We consider abelain subgroups of small index in finite groups. More generally, we consider subgroups such that the product of their index by the index of their centralizer is small.


title: The induced permittivity increment of electrorheological fluids in an applied electric field in association with chain formation: A Brownian Dynamics simulation study

abstract: We report Brownian Dynamics simulation results for the relative permittivity of electrorheological (ER) fluids in an applied electric field. The relative permittivity of an ER fluid can be calculated from the Clausius-Mosotti (CM) equation in the small applied field limit. When a strong field is applied, however, the ER spheres are organized into chains and assemblies of chains in which case the ER spheres are polarized not only by the external field but by each other. This manifests itself in an enhanced dielectric response, e.g., in an increase in the relative permittivity. The correction to the relative permittivity and the time dependence of this correction is simulated on the basis of a model in which the ER particles are represented as polarizable spheres. In this model, the spheres are also polarized by each other in addition to the applied field. Our results are qualitatively similar to those obtained by Horv\'ath and Szalai experimentally (\textit{Phys. Rev. E}, \textbf{86}, 061403, 2012). We report characteristic time constants obtained from bi-exponential fits that can be associated with formation of pairs and short chains as well as with aggregation of chains. The electric field dependence of the induced dielectric increment reveals the same qualitative behavior that experiments did: three regions with different slopes corresponding to different aggregation processes are identified.


title: A New Weak Choice Principle

abstract: For every natural number $n$ we introduce a new weak choice principle $\mathrm{nRC_{fin}}$: Given any infinite set $x$, there is an infinite subset $y\subseteq x$ and a selection function $f$ that chooses an $n$-element subset from every finite $z\subseteq y$ containing at least $n$ elements. By constructing new permutation models built on a set of atoms obtained as Fra\"iss\'e limits, we will study the relation of $\mathrm{nRC_{fin}}$ to the weak choice principles $\mathrm{RC_m}$ (that has already been studied by Montenegro, Halbeisen and Tachtsis): Given any infinite set $x$, there is an infinite subset $y\subseteq x$ with a choice function $f$ on the family of all $m$-element subsets of $y$. Moreover, we prove a stronger analogue of Montenegros results when we study the relation between $\mathrm{nRC_{fin}}$ and $\mathrm{kC_{fin}^-}$ which is defined by: Given any infinite family $\mathcal{F}$ of finite sets of cardinality greater than $k$, there is an infinite subfamily $\mathcal{A}\subseteq \mathcal{F}$ with a selection function $f$ that chooses a $k$-element subset from each $A\in\mathcal{A}$.


title: Porcupine: A Synthesizing Compiler for Vectorized Homomorphic Encryption

abstract: Homomorphic encryption (HE) is a privacy-preserving technique that enables computation directly on encrypted data. Despite its promise, HE has seen limited use due to performance overheads and compilation challenges. Recent work has made significant advances to address the performance overheads but automatic compilation of efficient HE kernels remains relatively unexplored. This paper presents Porcupine, an optimizing compiler, and HE DSL named Quill to automatically generate HE code using program synthesis. HE poses three major compilation challenges: it only supports a limited set of SIMD-like operators, it uses long-vector operands, and decryption can fail if ciphertext noise growth is not managed properly. Quill captures the underlying HE operator behavior that enables Porcupine to reason about the complex trade-offs imposed by the challenges and generate optimized, verified HE kernels. To improve synthesis time, we propose a series of optimizations including a sketch design tailored to HE and instruction restriction to narrow the program search space. We evaluate Procupine using a set of kernels and show speedups of up to 51% (11% geometric mean) compared to heuristic-driven hand-optimized kernels. Analysis of Porcupine's synthesized code reveals that optimal solutions are not always intuitive, underscoring the utility of automated reasoning in this domain.


title: An Analytical Expression for the Effective Area of the Step-Index Single-Mode Optical Fiber

abstract: The nonlinear coefficient {\gamma} of the step-index single-mode fiber (SMF) is inversely proportional to the effective area of that fiber. An analytical expression for the effective area of the hybrid HE11-mode of the SMF is derived for the first time, using its exact electromagnetic field expressions. The effective area is found in terms of generalized hypergeometric functions. A simpler expression is also deduced under the weakly-guided fiber approximation, and is found to be in agreement with established theory.


title: Tasting Nuclear Pasta Made with Classical Molecular Dynamics Simulations

abstract: Nuclear clusters or voids in the inner crust of neutron stars were predicted to have various shapes collectively nicknamed nuclear pasta. The recent review in Ref. \cite{Lopez1} by L\'opez, Dorso and Frank summarized their systematic investigations into properties especially the morphological and thermodynamical phase transitions of the nuclear pasta within a Classical Molecular Dynamics model, providing further stimuli to find more observational evidences of the predicted nuclear pasta in neutron stars.


title: Scalable Optimization for Wind Farm Control using Coordination Graphs

abstract: Wind farms are a crucial driver toward the generation of ecological and renewable energy. Due to their rapid increase in capacity, contemporary wind farms need to adhere to strict constraints on power output to ensure stability of the electricity grid. Specifically, a wind farm controller is required to match the farm's power production with a power demand imposed by the grid operator. This is a non-trivial optimization problem, as complex dependencies exist between the wind turbines. State-of-the-art wind farm control typically relies on physics-based heuristics that fail to capture the full load spectrum that defines a turbine's health status. When this is not taken into account, the long-term viability of the farm's turbines is put at risk. Given the complex dependencies that determine a turbine's lifetime, learning a flexible and optimal control strategy requires a data-driven approach. However, as wind farms are large-scale multi-agent systems, optimizing control strategies over the full joint action space is intractable. We propose a new learning method for wind farm control that leverages the sparse wind farm structure to factorize the optimization problem. Using a Bayesian approach, based on multi-agent Thompson sampling, we explore the factored joint action space for configurations that match the demand, while considering the lifetime of turbines. We apply our method to a grid-like wind farm layout, and evaluate configurations using a state-of-the-art wind flow simulator. Our results are competitive with a physics-based heuristic approach in terms of demand error, while, contrary to the heuristic, our method prolongs the lifetime of high-risk turbines.


title: Higher-order Fabry-P\'erot Interferometer from Topological Hinge States

abstract: We propose an intrinsic 3D Fabry-P\'erot type interferometer, coined "higher-order interferometer", that utilizes the chiral hinge states of second-order topological insulators and cannot be equivalently mapped to 2D space because of higher-order topology. Quantum interference patterns in the two-terminal conductance of this interferometer are controllable not only by tuning the strength but also, particularly, by rotating the direction of the magnetic field applied perpendicularly to the transport direction. Remarkably, the conductance exhibits a characteristic beating pattern with multiple frequencies with respect to field strength or direction. Our novel interferometer provides feasible and robust magneto-transport signatures to probe the particular hinge states of higher-order topological insulators.


title: Parallel-in-time high-order multiderivative IMEX solvers

abstract: In this work, we present a novel class of parallelizable high-order time integration schemes for the approximate solution of additive ODEs. The methods achieve high order through a combination of a suitable quadrature formula involving multiple derivatives of the ODE's right-hand side and a predictor-corrector ansatz. The latter approach is designed in such a way that parallelism in time is made possible. We present thorough analysis as well as numerical results that showcase scaling opportunities of methods from this class of solvers.


title: The Complexity of Monitoring Hyperproperties

abstract: We study the runtime verification of hyperproperties, expressed in the temporal logic HyperLTL, as a means to inspect a system with respect to security polices. Runtime monitors for hyperproperties analyze trace logs that are organized by common prefixes in the form of a tree-shaped Kripke structure, or are organized both by common prefixes and by common suffixes in the form of an acyclic Kripke structure. Unlike runtime verification techniques for trace properties, where the monitor tracks the state of the specification but usually does not need to store traces, a monitor for hyperproperties repeatedly model checks the growing Kripke structure. This calls for a rigorous complexity analysis of the model checking problem over tree-shaped and acyclic Kripke structures. We show that for trees, the complexity in the size of the Kripke structure is L-complete independently of the number of quantifier alternations in the HyperLTL formula. For acyclic Kripke structures, the complexity is PSPACE-complete (in the level of the polynomial hierarchy that corresponds to the number of quantifier alternations). The combined complexity in the size of the Kripke structure and the length of the HyperLTL formula is PSPACE-complete for both trees and acyclic Kripke structures, and is as low as NC for the relevant case of trees and alternation-free HyperLTL formulas. Thus, the size and shape of both the Kripke structure and the formula have significant impact on the complexity of the model checking problem.


title: DyLoc: Dynamic Localization for Massive MIMO Using Predictive Recurrent Neural Networks

abstract: This paper presents a data-driven localization framework with high precision in time-varying complex multipath environments, such as dense urban areas and indoors, where GPS and model-based localization techniques come short. We consider the angle-delay profile (ADP), a linear transformation of channel state information (CSI), in massive MIMO systems and show that ADPs preserve users' motion when stacked temporally. We discuss that given a static environment, future frames of ADP time-series are predictable employing a video frame prediction algorithm. We express that a deep convolutional neural network (DCNN) can be employed to learn the background static scattering environment. To detect foreground changes in the environment, corresponding to path blockage or addition, we introduce an algorithm taking advantage of the trained DCNN. Furthermore, we present DyLoc, a data-driven framework to recover distorted ADPs due to foreground changes and to obtain precise location estimations. We evaluate the performance of DyLoc in several dynamic scenarios employing DeepMIMO dataset to generate geo-tagged CSI datasets for indoor and outdoor environments. We show that previous DCNN-based techniques fail to perform with desirable accuracy in dynamic environments, while DyLoc pursues localization precisely. Moreover, simulations show that as the environment gets richer in terms of the number of multipath, DyLoc gets more robust to foreground changes.


title: Pregeometry and euclidean quantum gravity

abstract: Einstein's general relativity can emerge from pregeometry, with the metric composed of more fundamental fields. We formulate euclidean pregeometry as a $SO(4)$ - Yang-Mills theory. In addition to the gauge fields we include a vector field in the vector representation of the gauge group. The gauge - and diffeomorphism - invariant kinetic terms for these fields permit a well-defined euclidean functional integral, in contrast to metric gravity with the Einstein-Hilbert action. The propagators of all fields are well behaved at short distances, without tachyonic or ghost modes. The long distance behavior is governed by the composite metric and corresponds to general relativity. In particular, the graviton propagator is free of ghost or tachyonic poles despite the presence of higher order terms in a momentum expansion of the inverse propagator. This pregeometry seems to be a valid candidate for euclidean quantum gravity, without obstructions for analytic continuation to a Minkowski signature of the metric.