title: Propagating Conditions and the Time of ICMEs Arrival: A Comparison of the Effective Acceleration Model with ENLIL and DBEM Models

abstract: The Effective Acceleration Model (EAM) predicts the Time-of-Arrival (ToA) of the Coronal Mass Ejection (CME) driven shock and the average speed within the sheath at 1 AU. The model is based on the assumption that the ambient solar wind interacts with the interplanetary CME (ICME) resulting in constant acceleration or deceleration. The upgraded version of the model (EAMv3), presented here, incorporates two basic improvements: (a) a new technique for the calculation of the acceleration (or deceleration) of the ICME from the Sun to 1 AU and (b) a correction for the CME plane-of-sky speed. A validation of the upgraded EAM model is performed via comparisons to predictions from the ensemble version of the Drag-Based model (DBEM) and the WSA-ENLIL+Cone ensemble model. A common sample of 16 CMEs/ICMEs, in 2013-2014, is used for the comparison. Basic performance metrics such as the mean absolute error (MAE), mean error (ME) and root mean squared error (RMSE) between observed and predicted values of ToA are presented. MAE for EAM model was 8.7$\pm$1.6 hours while for DBEM and ENLIL was 14.3$\pm$2.2 and 12.8$\pm$1.7 hours, respectively. ME for EAM was -1.4$\pm$2.7 hours in contrast with -9.7$\pm$3.4 and -6.1$\pm$3.3 hours from DBEM and ENLIL. We also study the hypothesis of stronger deceleration in the interplanetary (IP) space utilizing the EAMv3 and DBEM models. In particularly, the DBEM model perform better when a greater value of drag parameter, of order of a factor of 3, is used in contrast to previous studies. EAMv3 model shows a deceleration of ICMEs at greater distances, with a mean value of 0.72 AU.


title: SPAA: Stealthy Projector-based Adversarial Attacks on Deep Image Classifiers

abstract: Light-based adversarial attacks use spatial augmented reality (SAR) techniques to fool image classifiers by altering the physical light condition with a controllable light source, e.g., a projector. Compared with physical attacks that place hand-crafted adversarial objects, projector-based ones obviate modifying the physical entities, and can be performed transiently and dynamically by altering the projection pattern. However, subtle light perturbations are insufficient to fool image classifiers, due to the complex environment and project-and-capture process. Thus, existing approaches focus on projecting clearly perceptible adversarial patterns, while the more interesting yet challenging goal, stealthy projector-based attack, remains open. In this paper, for the first time, we formulate this problem as an end-to-end differentiable process and propose a Stealthy Projector-based Adversarial Attack (SPAA) solution. In SPAA, we approximate the real Project-and-Capture process using a deep neural network named PCNet, then we include PCNet in the optimization of projector-based attacks such that the generated adversarial projection is physically plausible. Finally, to generate both robust and stealthy adversarial projections, we propose an algorithm that uses minimum perturbation and adversarial confidence thresholds to alternate between the adversarial loss and stealthiness loss optimization. Our experimental evaluations show that SPAA clearly outperforms other methods by achieving higher attack success rates and meanwhile being stealthier, for both targeted and untargeted attacks.


title: Evolution of Digital Advertising Strategies during the 2020 US Presidential Primary

abstract: Political advertising on digital platforms has grown dramatically in recent years as campaigns embrace new ways of targeting supporters and potential voters. Previous scholarship shows that digital advertising has both positive effects on democratic politics through increased voter knowledge and participation, and negative effects through user manipulation, opinion echo-chambers, and diminished privacy. However, research on election campaign strategies has focused primarily on traditional media, such as television. Here, we examine how political campaign dynamics have evolved in response to the growth of digital media by analyzing the advertising strategies of US presidential election campaigns during the 2020 primary cycle. To identify geographic and temporal trends, we employ regression analyses of campaign spending across nearly 600,000 advertisements published on Facebook. We show that campaigns heavily target voters in candidates' home states during the "invisible primary" stage before shifting to states with early primaries.


title: GNN-XML: Graph Neural Networks for Extreme Multi-label Text Classification

abstract: Extreme multi-label text classification (XMTC) aims to tag a text instance with the most relevant subset of labels from an extremely large label set. XMTC has attracted much recent attention due to massive label sets yielded by modern applications, such as news annotation and product recommendation. The main challenges of XMTC are the data scalability and sparsity, thereby leading to two issues: i) the intractability to scale to the extreme label setting, ii) the presence of long-tailed label distribution, implying that a large fraction of labels have few positive training instances. To overcome these problems, we propose GNN-XML, a scalable graph neural network framework tailored for XMTC problems. Specifically, we exploit label correlations via mining their co-occurrence patterns and build a label graph based on the correlation matrix. We then conduct the attributed graph clustering by performing graph convolution with a low-pass graph filter to jointly model label dependencies and label features, which induces semantic label clusters. We further propose a bilateral-branch graph isomorphism network to decouple representation learning and classifier learning for better modeling tail labels. Experimental results on multiple benchmark datasets show that GNN-XML significantly outperforms state-of-the-art methods while maintaining comparable prediction efficiency and model size.


title: The evolution of the low-density HI intergalactic medium from z=3.6 to 0: Data, transmitted flux and HI column density

abstract: We present a new, uniform analysis of the HI transmitted flux (F) and HI column density (N(HI)) distribution in the low-density IGM as a function of redshift z for 0<z<3.6 using 55 HST/COS FUV (Delta(z)=7.2 at z<0.5), five HST/STIS+COS NUV (Delta(z)=1.3 at z~1) and 24 VLT/UVES and Keck/HIRES (Delta(z)=11.6 at 1.7<z<3.6) AGN spectra. We performed a consistent, uniform Voigt profile analysis to combine spectra taken with different instruments, to reduce systematics and to remove metal-line contamination. We confirm previously known conclusions on firmer quantitative grounds in particular by improving the measurements at z~1. Two flux statistics at 0<F<1, the mean HI flux and the flux probability distribution function (PDF), show that considerable evolution occurs from z=3.6 to z=1.5, after which it slows down to become effectively stable for z<0.5. However, there are large sightline variations. For the HI column density distribution function (CDDF, f proptional to N(HI)^(-beta)) at log (N(HI)/1cm^-2)=[13.5, 16.0], beta increases as z decreases from beta~1.60 at z~3.4 to beta~1.82 at z~0.1. The CDDF shape at lower redshifts can be reproduced by a small amount of clockwise rotation of a higher-z CDDF with a slightly larger CDDF normalisation. The absorption line number per z (dn/dz) shows a similar evolutionary break at z~1.5 as seen in the flux statistics. High-N(HI) absorbers evolve more rapidly than low-N(HI) absorbers to decrease in number or cross-section with time. The individual dn/dz shows a large scatter at a given z. The scatter increases toward lower z, possibly caused by a stronger clustering at lower z.


title: Understanding Learned Reward Functions

abstract: In many real-world tasks, it is not possible to procedurally specify an RL agent's reward function. In such cases, a reward function must instead be learned from interacting with and observing humans. However, current techniques for reward learning may fail to produce reward functions which accurately reflect user preferences. Absent significant advances in reward learning, it is thus important to be able to audit learned reward functions to verify whether they truly capture user preferences. In this paper, we investigate techniques for interpreting learned reward functions. In particular, we apply saliency methods to identify failure modes and predict the robustness of reward functions. We find that learned reward functions often implement surprising algorithms that rely on contingent aspects of the environment. We also discover that existing interpretability techniques often attend to irrelevant changes in reward output, suggesting that reward interpretability may need significantly different methods from policy interpretability.


title: A Decision Tree Lifted Domain for Analyzing Program Families with Numerical Features (Extended Version)

abstract: Lifted (family-based) static analysis by abstract interpretation is capable of analyzing all variants of a program family simultaneously, in a single run without generating any of the variants explicitly. The elements of the underlying lifted analysis domain are tuples, which maintain one property per variant. Still, explicit property enumeration in tuples, one by one for all variants, immediately yields combinatorial explosion. This is particularly apparent in the case of program families that, apart from Boolean features, contain also numerical features with big domains, thus admitting astronomic configuration spaces. The key for an efficient lifted analysis is proper handling of variability-specific constructs of the language (e.g., feature-based runtime tests and #if directives). In this work, we introduce a new symbolic representation of the lifted abstract domain that can efficiently analyze program families with numerical features. This makes sharing between property elements corresponding to different variants explicitly possible. The elements of the new lifted domain are constraint-based decision trees, where decision nodes are labeled with linear constraints defined over numerical features and the leaf nodes belong to an existing single-program analysis domain. To illustrate the potential of this representation, we have implemented an experimental lifted static analyzer, called SPLNUM^2Analyzer, for inferring invariants of C programs. It uses existing numerical domains (e.g., intervals, octagons, polyhedra) from the APRON library as parameters. An empirical evaluation on benchmarks from SV-COMP and BusyBox yields promising preliminary results indicating that our decision trees-based approach is effective and outperforms the tuple-based approach, which is used as a baseline lifted analysis based on abstract interpretation.


title: The preservability of the curvature-adaptedness along the mean curvature flow

abstract: In this paper, we investigate the preservability of the curvature-adaptedness along the mean curvature flow starting from a compact curvature-adapted hypersurface in locally symmetric spaces, where the curvature-adaptedness means that the shape operator and the normal Jacobi operator of the hypersurface commute.


title: Performance of Dense Coding and Teleportation for Random States --Augmentation via Pre-processing

abstract: In order to understand the resourcefulness of a natural quantum system in quantum communication tasks, we study the dense coding capacity (DCC) and teleportation fidelity (TF) of Haar uniformly generated random multipartite states of various ranks. We prove that when a rank-2 two-qubit state, a Werner state, and a pure state possess the same amount of entanglement, the DCC of a rank-2 state belongs to the envelope made by pure and Werner states. In a similar way, we obtain an upper bound via the generalized Greenberger-Horne-Zeilinger state for rank-2 three-qubit states when the dense coding with two senders and a single receiver is performed and entanglement is measured in the senders:receiver bipartition. The normalized frequency distribution of DCC for randomly generated two-, three- and four-qubit density matrices with global as well as local decodings at the receiver's end are reported. The estimation of mean DCC for two-qubit states is found to be in good agreement with the numerical simulations. Universally, we observe that the performance of random states for dense coding as well as teleportation decreases with the increase of the rank of states which we have shown to be surmounted by the local pre-processing operations performed on the shared states before starting the protocols, irrespective of the rank of the states. The local pre-processing employed here is based on positive operator valued measurements along with classical communication and we show that unlike dense coding with two-qubit random states, senders' operations are always helpful to probabilistically enhance the capabilities of implementing dense coding as well as teleportation.


title: Diversity of phase transitions and phase separations in active fluids

abstract: Active matter is not only indispensable to our understanding of diverse biological processes, but also provides a fertile ground for discovering novel physics. Many emergent properties impossible for equilibrium systems have been demonstrated in active systems. These emergent features include motility-induced phase separation, long-ranged ordered (collective motion) phase in two dimensions, and order-disorder phase co-existences (banding and reverse-banding regimes). Here, we unify these diverse phase transitions and phase co-existences into a single formulation based on generic hydrodynamic equations for active fluids. We also reveal a novel co-moving co-existence phase and a putative novel critical point.


title: Virtual Classrooms and Real Harms: Remote Learning at U.S. Universities

abstract: Universities have been forced to rely on remote educational technology to facilitate the rapid shift to online learning. In doing so, they acquire new risks of security vulnerabilities and privacy violations. To help universities navigate this landscape, we develop a model that describes the actors, incentives, and risks, informed by surveying 49 educators and 14 administrators at U.S. universities. Next, we develop a methodology for administrators to assess security and privacy risks of these products. We then conduct a privacy and security analysis of 23 popular platforms using a combination of sociological analyses of privacy policies and 129 state laws, alongside a technical assessment of platform software. Based on our findings, we develop recommendations for universities to mitigate the risks to their stakeholders.


title: MuPix10: First Results from the Final Design

abstract: Many years of research and development of High Voltage Monolithic Active Pixel Sensors (HVMAPS) have culminated in the final design for the Mu3e pixel sensor. MuPix10 is a fully monolithic sensor with an active pixel matrix size of $20\times20\,\mathrm{mm}^2$ produced in the $180\,\mathrm{nm}$ HV-CMOS process at TSI Semiconductors. The pixel size is $80\times80\,\mathrm{\mu m}^2$. Hits are read out using a column-drain architecture and sent over up to four serial links with up to $1.6\,\left.\mathrm{Gbit}\middle/\mathrm{s}\right.$ each. By means of DC/DC converters and exclusive usage of on-chip biasing, MuPix10 is fully operable with a minimal set of electrical connections. This is an integral requirement by the Mu3e experiment since it enables the construction of ultra-thin pixel modules with $0.1\,$% of a radiation length per layer. First results from lab characterisation and testbeam campaigns are presented.


title: A Linear-Time Algorithm for Minimum $k$-Hop Dominating Set of a Cactus Graph

abstract: Given a graph $G=(V,E)$ and an integer $k \ge 1$, a $k$-hop dominating set $D$ of $G$ is a subset of $V$, such that, for every vertex $v \in V$, there exists a node $u \in D$ whose hop-distance from $v$ is at most $k$. A $k$-hop dominating set of minimum cardinality is called a minimum $k$-hop dominating set. In this paper, we present linear-time algorithms that find a minimum $k$-hop dominating set in unicyclic and cactus graphs. To achieve this, we show that the $k$-dominating set problem on unicycle graph reduces to the piercing circular arcs problem, and show a linear-time algorithm for piercing sorted circular arcs, which improves the best known $O(n\log n)$-time algorithm.


title: On Vacuum Fluctuations in Quantum Gravity and Interferometer Arm Fluctuations

abstract: We propose a simple model of spacetime vacuum fluctuations motivated by AdS/CFT, where the vacuum is described by a thermal density matrix, $\rho = \frac{e^{-K}}{\mbox{Tr}(e^{-K})}$ with $K$ the modular Hamiltonian. In AdS/CFT, both the expectation value of $K$ and its fluctuations $\langle \Delta K^2 \rangle$ have been calculated; both obey an area law identical to the Bekenstein-Hawking area law of black hole mechanics: $\langle K \rangle = \langle \Delta K^2 \rangle = \frac{A}{4 G_N}$, where $A$ is the area of an (extremal) entangling surface. It has also been shown that $\Delta K$ gravitates in AdS, and hence generates metric fluctuations. These theoretical results are intriguing, but it is not known how to precisely extend such ideas about holographic quantum gravity to ordinary flat space. We take the approach of considering whether experimental signatures in metric fluctuations could determine properties of the vacuum of quantum gravity in flat space. In particular, we propose a theoretical model motived by the AdS/CFT calculations that reproduces the most important features of modular Hamiltonian fluctuations; the model consists of a high occupation number bosonic degree of freedom. We show that if this theory couples through ordinary gravitational couplings to the mirrors in an interferometer with strain sensitivity similar to what will be available for gravitational waves, vacuum fluctuations could be observable.


title: Extreme learning machine collocation for the numerical solution of elliptic PDEs with sharp gradients

abstract: We introduce a new numerical method based on machine learning to approximate the solution of elliptic partial differential equations with collocation using a set of sigmoidal functions. We show that a feedforward neural network with a single hidden layer with sigmoidal functions and fixed, random, internal weights and biases can be used to compute accurately a collocation solution. The choice to fix internal weights and bias leads to the so-called Extreme Learning Machine network. We discuss how to determine the range for both internal weights and biases in order to obtain a good underlining approximating space, and we explore the required number of collocation points. We demonstrate the efficiency of the proposed method with several one-dimensional diffusion-advection-reaction problems that exhibit steep behaviors, such as boundary layers. The boundary conditions are imposed directly as collocation equations. We point out that there is no need of training the network, as the proposed numerical approach results to a linear problem that can be easily solved using least-squares. Numerical results show that the proposed method achieves a good accuracy. Finally, we compare the proposed method with finite differences and point out the significant improvements in terms of computational cost, thus avoiding the time-consuming training phase.


title: Domain walls and CP violation with left right supersymmetry:implications for leptogenesis and electron EDM

abstract: Low scale leptogenesis scenarios are difficult to verify due to our inability to relate the parameters involved in the early universe processes with the low energy or collider observables. Here we show that one can in principle relate the parameters giving rise to the transient $CP$ violating phase involved in leptogenesis with those that can be deduced from the observation of electric dipole moment (EDM) of the electron. We work out the details of this in the context of the left right symmetric supersymmetric model (LRSUSY) which provides a strong connection between such parameters. In particular, we show that baryon asymmetry requirements imply the scale $M_{B-L}$ of $U(1)_{B-L}$ symmetry breaking to be larger than $10^{4.5}~\mathrm{GeV}$. Moreover the scale $M_R$ of $SU(2)_R$ symmetry breaking is tightly constrained to lie in a narrow band significantly below $M_{B-L}^2 / M_{EW}$. These are the most stringent constraints on the parameter space of LRSUSY model being considered.


title: Search for double $\beta$-decay modes of $^{64}$Zn using purified zinc

abstract: The production of ultra-pure raw material is a crucial step to ensure the required background level in rare event searches. In this work, we establish an innovative technique developed to produce high-purity (99.999%) granular zinc. We demonstrate the effectiveness of the refining procedure by measuring the internal contaminations of the purified zinc with a high-purity germanium detector at the Laboratori Nazionali del Gran Sasso. The total activity of cosmogenic activated nuclides is measured at the level of a few mBq/kg, as well as limits on naturally occurring radionuclides are set to less than mBq/kg. The excellent radiopurity of the zinc sample allows us to search for electron capture with positron emission and neutrinoless double electron capture of $^{64}$Zn, setting the currently most stringent lower limits on their half-lives, $T_{1/2}^{\varepsilon\beta^+} > 2.7 \times 10^{21}\text{yr}$ (90% C.I.), and $T_{1/2}^{0\nu2\varepsilon}> 2.6 \times 10^{21}\text{yr}$ (90% C.I.), respectively.


title: Hindsight and Sequential Rationality of Correlated Play

abstract: Driven by recent successes in two-player, zero-sum game solving and playing, artificial intelligence work on games has increasingly focused on algorithms that produce equilibrium-based strategies. However, this approach has been less effective at producing competent players in general-sum games or those with more than two players than in two-player, zero-sum games. An appealing alternative is to consider adaptive algorithms that ensure strong performance in hindsight relative to what could have been achieved with modified behavior. This approach also leads to a game-theoretic analysis, but in the correlated play that arises from joint learning dynamics rather than factored agent behavior at equilibrium. We develop and advocate for this hindsight rationality framing of learning in general sequential decision-making settings. To this end, we re-examine mediated equilibrium and deviation types in extensive-form games, thereby gaining a more complete understanding and resolving past misconceptions. We present a set of examples illustrating the distinct strengths and weaknesses of each type of equilibrium in the literature, and prove that no tractable concept subsumes all others. This line of inquiry culminates in the definition of the deviation and equilibrium classes that correspond to algorithms in the counterfactual regret minimization (CFR) family, relating them to all others in the literature. Examining CFR in greater detail further leads to a new recursive definition of rationality in correlated play that extends sequential rationality in a way that naturally applies to hindsight evaluation.


title: Three-body capture, ejection, and the demographics of bound objects in binary systems

abstract: We study the capture of light objects of arbitrary velocity by binary systems. Extending results for the capture of comets in the solar system, we develop a simple geometric characterization of the capture cross section, leading directly to the distribution of orbital parameters of captured objects. We use the same framework to study the lifetimes of these bound orbits prior to ejection, and find that a simplified version of the \"Opik--Arnold approach readily yields a closed-form estimate for the ejection rate that agrees well with numerical experiments. Without any detailed-balance assumptions, our results make manifest the characteristics of close encounters leading to capture and ejection. As an application of our results, we demonstrate the estimation of the equilibrium population of captured dark matter particles in a binary system.


title: Neurosymbolic AI: The 3rd Wave

abstract: Current advances in Artificial Intelligence (AI) and Machine Learning (ML) have achieved unprecedented impact across research communities and industry. Nevertheless, concerns about trust, safety, interpretability and accountability of AI were raised by influential thinkers. Many have identified the need for well-founded knowledge representation and reasoning to be integrated with deep learning and for sound explainability. Neural-symbolic computing has been an active area of research for many years seeking to bring together robust learning in neural networks with reasoning and explainability via symbolic representations for network models. In this paper, we relate recent and early research results in neurosymbolic AI with the objective of identifying the key ingredients of the next wave of AI systems. We focus on research that integrates in a principled way neural network-based learning with symbolic knowledge representation and logical reasoning. The insights provided by 20 years of neural-symbolic computing are shown to shed new light onto the increasingly prominent role of trust, safety, interpretability and accountability of AI. We also identify promising directions and challenges for the next decade of AI research from the perspective of neural-symbolic systems.


title: INeRF: Inverting Neural Radiance Fields for Pose Estimation

abstract: We present iNeRF, a framework that performs mesh-free pose estimation by "inverting" a Neural RadianceField (NeRF). NeRFs have been shown to be remarkably effective for the task of view synthesis - synthesizing photorealistic novel views of real-world scenes or objects. In this work, we investigate whether we can apply analysis-by-synthesis via NeRF for mesh-free, RGB-only 6DoF pose estimation - given an image, find the translation and rotation of a camera relative to a 3D object or scene. Our method assumes that no object mesh models are available during either training or test time. Starting from an initial pose estimate, we use gradient descent to minimize the residual between pixels rendered from a NeRF and pixels in an observed image. In our experiments, we first study 1) how to sample rays during pose refinement for iNeRF to collect informative gradients and 2) how different batch sizes of rays affect iNeRF on a synthetic dataset. We then show that for complex real-world scenes from the LLFF dataset, iNeRF can improve NeRF by estimating the camera poses of novel images and using these images as additional training data for NeRF. Finally, we show iNeRF can perform category-level object pose estimation, including object instances not seen during training, with RGB images by inverting a NeRF model inferred from a single view.


title: Asymptotic stability of small ground states for NLS under random perturbations

abstract: We consider the cubic Schr\"odinger equation on the euclidean space perturbed by a short-range potential $V$. The presence of a negative simple eigenvalue for $-\Delta+V$ gives rise to a curve of small and localized nonlinear ground states that yield some time-periodic solutions known to be asymptotically stable in the energy space. We study the persistence of these coherent states under rough perturbations. We shall construct a large measure set of small scaling-supercritical solutions below the energy space that display some asymptotic stability behavior. The main difficulty is the need to handle the interactions of localized and dispersive terms in the modulation equations. To do so, we use a critical-weighted strategy to combine probabilistic nonlinear estimates in critical spaces based on $U^p, V^q$ (controlling higher order terms) with some local energy decay estimates (controlling lower order terms). We also revisit in the perturbed setting the analysis of B\'enyi, Oh and Pocovnicu on the probabilistic global well-posedness and scattering for small supercritical initial data. We use a distorted Fourier transform and semiclassical functional calculus to generalize probabilistic and bilinear Strichartz estimates.


title: Automatic Standardization of Colloquial Persian

abstract: The Iranian Persian language has two varieties: standard and colloquial. Most natural language processing tools for Persian assume that the text is in standard form: this assumption is wrong in many real applications especially web content. This paper describes a simple and effective standardization approach based on sequence-to-sequence translation. We design an algorithm for generating artificial parallel colloquial-to-standard data for learning a sequence-to-sequence model. Moreover, we annotate a publicly available evaluation data consisting of 1912 sentences from a diverse set of domains. Our intrinsic evaluation shows a higher BLEU score of 62.8 versus 61.7 compared to an off-the-shelf rule-based standardization model in which the original text has a BLEU score of 46.4. We also show that our model improves English-to-Persian machine translation in scenarios for which the training data is from colloquial Persian with 1.4 absolute BLEU score difference in the development data, and 0.8 in the test data.


title: The moving frame method for iterated-integrals: orthogonal invariants

abstract: Geometric features, robust to noise, of curves in Euclidean space are of great interest for various applications such as machine learning and image analysis. We apply the Fels-Olver's moving frame method (for geometric features) paired with the log-signature transform (for robust features) to construct a set of integral invariants under rigid motions for curves in $\mathbb{R}^d$ from the iterated-integral signature. In particular we show that one can algorithmically construct a set of invariants that characterize the equivalence class of the truncated iterated-integrals signature under orthogonal transformations which yields a characterization of a curve in $\mathbb{R}^d$ under rigid motions (and tree-like extensions) and an explicit method to compare curves up to these transformations.


title: Uso di un software di Geometria Dinamica nella formazione dei futuri insegnanti

abstract: In this paper we make some considerations about using a dynamic geometry software for teaching history of mathematica at university level. After a short introduction to the software GeoGebra, we discuss four activities. The first is an exploration of the menu of the tools menu following the first part of the first book of the elements of Euclid. The second is an activity about pentagons and pentagrams which aims at finding a more natural way of discovering their properties than that suggested by Euclid as an illustration of the method of synthesis and analysis. The third is an activity on circle inversion and its connection with polarity and the construction of the fourth harmonic point after three given point on a line. The activity continues with the reading of Poncelet's ideas about ideal secants of conics and how to use GeoGebra to make sense of them. The fourth activity shows how to use the possibility to move objects in GeoGebra to make sense of the geometric intuition behind Cremona and Hirst letters about the Bellavitis-Hirst inversion of a circle. We discuss the reasons why we believe that these kind of activities may help students to read and understand original writings in projective geometry and why we believe that they may be very effective in the training of prospective teachers of mathematics.


title: Differential Projective Modules over Differential Rings, II

abstract: Differential modules over a commutative differential ring R which are finitely generated projective as ring modules, with differential homomorphisms, form an additive category, so their isomorphism classes form a monoid. We study the quotient monoid of this monoid by the submonoid of isomorphism classes of free modules with component wise derivation. This quotient monoid has the reduced K group of R (ignoring the derivation) as an image and contains the reduced K group of the constants of R as its subgroup of units. This monoid provides a description of the isomorphism classes of differential projective R modules up to an equivalence.


title: The fate of $V_1$ vector leptoquarks: the impact of future flavour data

abstract: Motivated by the recent experimental progress on the $B$-meson decay anomalies (in particular the angular observables in $B\to K^\ast\mu\mu$), we rely on a simplified-model approach to study the prospects of vector leptoquarks in what concerns numerous flavour observables, identifying several promising decay modes which would allow to (indirectly) probe such an extension. Our findings suggest that the confirmation of the $B$-meson decay anomalies, in parallel with positive signals (at Belle II or LHCb) for $\tau\to \phi \mu$, $B_{(s)}$-meson decays to $\tau^+ \tau^-$ and $\tau^+ \mu^-$ ($\tau^+ e^-$) final states, as well as an observation of certain charged lepton flavour violation decays (at COMET or Mu2e), would contribute to strengthen the case for this scenario. We also illustrate how the evolution of the experimental determination of $R_{D^{(*)}}$ could be instrumental in falsifying an explanation of the anomalous $B$-meson decay data via a vector $V_1$ leptoquark.


title: Absence of a Dirac gap in ferromagnetic Cr$_x$(Bi$_{0.1}$Sb$_{0.9}$)$_{2-x}$Te$_3$

abstract: Magnetism breaks the time reversal symmetry expected to open a Dirac gap in 3D topological insulators that consequently leads to quantum anomalous Hall effect. The most common approach of inducing ferromagnetic state is by doping magnetic 3$d$ elements into bulk of 3D topological insulators. In Cr$_{0.15}$(Bi$_{0.1}$Sb$_{0.9}$)$_{1.85}$Te$_3$, the material where the quantum anomalous Hall effect was initially discovered at temperatures much lower than the ferromagnetic transition, $T_C$, the scanning tunneling microscopy studies have reported a large Dirac gap $\sim20-100$ meV. The discrepancy between the low temperature of quantum anomalous Hall effect ($\ll T_C$) and large spectroscopic Dirac gaps ($\gg T_C$) found in magnetic topological insulators remains puzzling. Here, we used angle-resolved photoemission spectroscopy to study the surface electronic structure of pristine and potassium doped surface of Cr$_{0.15}$(Bi$_{0.1}$Sb$_{0.9}$)$_{1.85}$Te$_3$. Upon potassium deposition, the $p$-type surface state of pristine sample was turned into an $n$-type, allowing spectroscopic observation of Dirac point. We find a gapless surface state, with no evidence of a large Dirac gap reported in tunneling studies.


title: Relaxation in an Extended Bosonic Josephson Junction

abstract: We present a detailed analysis of the relaxation dynamics in an extended bosonic Josephson junction. We show that stochastic classical field simulations using Gross-Pitaevskii equations in three spatial dimensions reproduce the main experimental findings of M. Pigneur et al., Phys. Rev. Lett. 120, 173601 (2018). We give an analytic solution describing the short time evolution through multimode dephasing. For longer times, the observed relaxation to a phase locked state is caused by nonlinear dynamics beyond the sine-Gordon model, induced by the longitudinal confinement potential and persisting even at zero temperature. Finally, we analyze different experimentally relevant trapping geometries to mitigate these effects. Our results provide the basis for future experimental implementations aiming to study nonlinear and quantum effects of the relaxation in extended bosonic Josephson junctions.


title: Absolute determination of the single-photon optomechanical coupling rate via a Hopf bifurcation

abstract: We establish a method for the determination of the single-photon optomechanical coupling rate, which characterizes the radiation pressure interaction in an optomechanical system. The estimation of the rate with which a mechanical oscillator, initially in a thermal state, undergoes a Hopf bifurcation, and reaches a limit cycle, allows us to determine the single-photon optomechanical coupling rate in a simple and consistent way. Most importantly, and in contrast to other methods, our method does not rely on knowledge of the system's bath temperature and on a calibration of the signal. We provide the theoretical framework, and experimentally validate this method, providing a procedure for the full characterization of an optomechanical system, which could be extended outside cavity optomechanics, whenever a resonator is driven into a limit cycle by the appropriate interaction with another degree of freedom.


title: Quantum Control in the Unitary Sphere: Lambda-S1 and its Categorical Model

abstract: In a recent paper, a realizability technique has been used to give a semantics of a quantum lambda calculus. Such a technique gives rise to an infinite number of valid typing rules, without giving preference to any subset of those. In this paper, we introduce a valid subset of typing rules, defining an expressive enough quantum calculus. Then, we propose a categorical semantics for it. Such a semantics consists of an adjunction between the category of distributive-action spaces of value distributions (that is, linear combinations of values in the lambda calculus), and the category of sets of value distributions.


title: Stable Big Bang formation for Einstein's equations: The complete sub-critical regime

abstract: For $(t,x) \in (0,\infty)\times\mathbb{T}^D$, the generalized Kasner solutions are a family of explicit solutions to various Einstein-matter systems that start out smooth but then develop a Big Bang singularity as $t \downarrow 0$, i.e., curvature blowup along a spacelike hypersurface. The family is parameterized by the Kasner exponents $\widetilde{q}_1,\cdots,\widetilde{q}_D \in \mathbb{R}$, which satisfy two algebraic constraints. There are heuristics in the mathematical physics literature, going back more than 50 years, suggesting that the Big Bang formation should be stable under perturbations of the Kasner initial data, given say at $\lbrace t = 1 \rbrace$, as long as the exponents are "sub-critical" in the following sense: $\mathop{\max_{I,J,B=1,\cdots,D}}_{I < J} \{\widetilde{q}_I+\widetilde{q}_J-\widetilde{q}_B\}<1$. Previous works have shown the stability of the singularity under stronger assumptions: 1) the Einstein-scalar field system with $D = 3$ and $\widetilde{q}_1 \approx \widetilde{q}_2 \approx \widetilde{q}_3 \approx 1/3$ or 2) the Einstein-vacuum equations for $D \geq 39$ with $\max_{I=1,\cdots,D} |\widetilde{q}_I| < 1/6$. We prove that the Kasner singularity is dynamically stable for \emph{all} sub-critical Kasner exponents, thereby justifying the heuristics in the full regime where stable monotonic-type curvature blowup is expected. We treat the $1+D$-dimensional Einstein-scalar field system for $D \geq3$ and the $1+D$ dimensional Einstein-vacuum equations for $D \geq 10$. Moreover, for the Einstein-vacuum equations in $1+3$ dimensions, where instabilities are in general expected, we prove that all singular Kasner solutions have stable Big Bangs under polarized $U(1)$-symmetric perturbations of their initial data. Our results hold for open sets of initial data in Sobolev spaces without symmetry, apart from our work on polarized $U(1)$-symmetric solutions.


title: Constraining the magnetic field properties of Bok globule B335 using SOFIA/HAWC+

abstract: Thanks to their well-defined shape and mostly isolated locations, Bok globules are suitable objects for studying the physics of low-mass star formation. To study the magnetic field of the prototypical Bok globule B335, we obtained a spatially resolved polarization map with SOFIA/HAWC+ at a wavelength of 214$\,\mu$m. For the first time, these observations reveal that polarization holes in Bok globules, that is, the decrease in polarization degree towards their dense centers, also occur in the far-infrared wavelength regime. The observed polarization pattern is uniform with a mean polarization angle of 48$^\circ\pm $26$^\circ$ and a magnetic field strength of $\sim$ 142$\,\mu$G. Moreover, we use complementary polarimetic data for B335 obtained at near-infrared to millimeter wavelengths to analyze and constrain the magnetic field across different scales. By applying the 3D Monte-Carlo radiative transfer code POLARIS (Reissl et al. 2016), we developed a model for the density and magnetic field structure as well as for the dust properties of this globule. We conclude that the column density towards the center of B335 is too low to cause the observed polarization hole in B335 via dichroic absorption (Brauer et al. 2016). Furthermore, we conclude that the effect of self-scattering has no significant impact on the observed polarization. Adopting dust-grain alignment via the radiative torque mechanism, a combination of the interstellar radiation field and the central star as radiation sources is consistent with the decrease in polarization degree at the outer regions of B335 ($\approx\,$10$^4\,$au from the core). However, the model fails to explain the low polarization degree within the inner 5000 au.


title: Galactic Bar Resonances Inferred from Kinematically Hot Stars in Gaia EDR3

abstract: Using a numerical simulation of an isolated barred disc galaxy, we first demonstrate that the resonances of the inner bar structure induce more prominent features in the action space distribution for the kinematically hotter stars, which are less sensitive to the local perturbation, such as the transient spiral arms. Then, we analyse the action distribution for the kinematically hotter stars selected from the Gaia EDR3 data as the stars with higher values of radial and vertical actions. We find several resonance features, including two new features, in the angular momentum distribution similar to what are seen in our numerical simulations. We show that the bar pattern speeds of about $\Omega_{\rm bar}\sim34$~km~s$^{-1}$~kpc$^{-1}$ and 42~km~s$^{-1}$~kpc$^{-1}$ explain all these features equally well. The resonance features we find correspond to the inner 4:1, co-rotation, outer 4:1, outer Lindblad and outer 4:3 (co-rotation, outer 4:1, outer Lindblad, outer 4:3 and outer 1:1) resonances, when $\Omega_{\rm bar}\sim34$ (42) km~s$^{-1}$~kpc$^{-1}$ is assumed.


title: GeV Scale Inelastic Dark Matter with Dark Photon Mediator via Direct Detection and Cosmological/Laboratory Constraints

abstract: We propose a new candidate of GeV scale inelastic dark matter (DM). Our construction has an anomaly-free $U(1)_X^{}$ gauge group with dark photon mediator, and can realize either scalar or fermionic inelastic DM. It is highly predictive and testable. We study the scattering rate of light inelastic DM with electrons in the XENON1T experiment and with nuclei in the XENON1T, CRESST-III, CDEX-1B and DarkSide-50 experiments. We resolve the recent XENON1T anomaly via electron recoil detection. Combining the XENON1T constraints from both electron recoils and nuclear recoils (including Migdal effect), we predict the inelastic DM mass $\lesssim 1.5$GeV. We further analyze the bounds by the DM relic abundance, the lifetime of the heavier DM component, and laboratory constraints, from which we identify the viable parameter space for the future probe. This provides an important benchmark for the theories and experimental tests of GeV scale inelastic DM.


title: Planck spectroscopy

abstract: All spectrometers rely on some mechanism to achieve spectral selectivity; common examples include gratings, prisms, and interferometers with moving mirrors. We experimentally demonstrated and validated a spectroscopic technique -- here dubbed Planck spectroscopy -- that measures the spectral emissivity of a surface using only a temperature-controlled stage and a detector, without any wavelength-selective optical components. Planck spectroscopy involves the measurement of temperature-dependent thermally emitted power, where the spectral selectivity is realized via the temperature- and wavelength dependence of Planck's law. We experimentally demonstrated and validated Planck spectroscopy in the mid infrared, for wavelengths from 3 to 13 um -- limited primarily by the bandwidth of our detector -- with resolution of approximately 1 um. The minimalistic setup of Planck spectroscopy can be implemented using infrared cameras to achieve low-cost infrared hyperspectral imaging and imaging ellipsometry.


title: Flatland-RL : Multi-Agent Reinforcement Learning on Trains

abstract: Efficient automated scheduling of trains remains a major challenge for modern railway systems. The underlying vehicle rescheduling problem (VRSP) has been a major focus of Operations Research (OR) since decades. Traditional approaches use complex simulators to study VRSP, where experimenting with a broad range of novel ideas is time consuming and has a huge computational overhead. In this paper, we introduce a two-dimensional simplified grid environment called "Flatland" that allows for faster experimentation. Flatland does not only reduce the complexity of the full physical simulation, but also provides an easy-to-use interface to test novel approaches for the VRSP, such as Reinforcement Learning (RL) and Imitation Learning (IL). In order to probe the potential of Machine Learning (ML) research on Flatland, we (1) ran a first series of RL and IL experiments and (2) design and executed a public Benchmark at NeurIPS 2020 to engage a large community of researchers to work on this problem. Our own experimental results, on the one hand, demonstrate that ML has potential in solving the VRSP on Flatland. On the other hand, we identify key topics that need further research. Overall, the Flatland environment has proven to be a robust and valuable framework to investigate the VRSP for railway networks. Our experiments provide a good starting point for further research and for the participants of the NeurIPS 2020 Flatland Benchmark. All of these efforts together have the potential to have a substantial impact on shaping the mobility of the future.


title: AutoSelect: Automatic and Dynamic Detection Selection for 3D Multi-Object Tracking

abstract: 3D multi-object tracking is an important component in robotic perception systems such as self-driving vehicles. Recent work follows a tracking-by-detection pipeline, which aims to match past tracklets with detections in the current frame. To avoid matching with false positive detections, prior work filters out detections with low confidence scores via a threshold. However, finding a proper threshold is non-trivial, which requires extensive manual search via ablation study. Also, this threshold is sensitive to many factors such as target object category so we need to re-search the threshold if these factors change. To ease this process, we propose to automatically select high-quality detections and remove the efforts needed for manual threshold search. Also, prior work often uses a single threshold per data sequence, which is sub-optimal in particular frames or for certain objects. Instead, we dynamically search threshold per frame or per object to further boost performance. Through experiments on KITTI and nuScenes, our method can filter out $45.7\%$ false positives while maintaining the recall, achieving new S.O.T.A. performance and removing the need for manually threshold tuning.


title: Probing Few-Shot Generalization with Attributes

abstract: Despite impressive progress in deep learning, generalizing far beyond the training distribution is an important open challenge. In this work, we consider few-shot classification, and aim to shed light on what makes some novel classes easier to learn than others, and what types of learned representations generalize better. To this end, we define a new paradigm in terms of attributes -- simple building blocks of which concepts are formed -- as a means of quantifying the degree of relatedness of different concepts. Our empirical analysis reveals that supervised learning generalizes poorly to new attributes, but a combination of self-supervised pretraining with supervised finetuning leads to stronger generalization. The benefit of self-supervised pretraining and supervised finetuning is further investigated through controlled experiments using random splits of the attribute space, and we find that predictability of test attributes provides an informative estimate of a model's generalization ability.


title: Encoding Classical Information in Gauge Subsystems of Quantum Codes

abstract: We show how to construct hybrid quantum-classical codes from subsystem codes by encoding the classical information into the gauge qudits using gauge fixing. Unlike previous work on hybrid codes, we allow for two separate minimum distances, one for the quantum information and one for the classical information. We give an explicit construction of hybrid codes from two classical linear codes using Bacon-Casaccino subsystem codes, as well as several new examples of good hybrid code.


title: Self-Supervised Learning of Lidar Segmentation for Autonomous Indoor Navigation

abstract: We present a self-supervised learning approach for the semantic segmentation of lidar frames. Our method is used to train a deep point cloud segmentation architecture without any human annotation. The annotation process is automated with the combination of simultaneous localization and mapping (SLAM) and ray-tracing algorithms. By performing multiple navigation sessions in the same environment, we are able to identify permanent structures, such as walls, and disentangle short-term and long-term movable objects, such as people and tables, respectively. New sessions can then be performed using a network trained to predict these semantic labels. We demonstrate the ability of our approach to improve itself over time, from one session to the next. With semantically filtered point clouds, our robot can navigate through more complex scenarios, which, when added to the training pool, help to improve our network predictions. We provide insights into our network predictions and show that our approach can also improve the performances of common localization techniques.


title: Anisotropic proximity-induced superconductivity and edge supercurrent in Kagome metal, K1-xV3Sb5

abstract: Materials with transition metals in triangular lattices are of great interest for their potential combination of strong correlation, exotic magnetism and electronic topology. Kagome nets are of particular importance since the discovery of geometrically frustrated magnetism and topological band structures in crystals like Herbertsmithite and Fe3Sn2, respectively. KV3Sb5 was discovered to be a layered topological metal with a Kagome net of vanadium. Here, we fabricated Josephson Junctions (JJ) of K1-xV3Sb5 and induced superconductivity over long junction lengths. Through magnetoresistance and current vs. phase measurements, we observed magnetic field sweeping direction dependent magnetoresistance, and an anisotropic interference pattern with a Fraunhofer pattern for in-plane magnetic field, but a suppression of critical current for out-of-plane magnetic field. These results indicate an anisotropic internal magnetic field in K1-xV3Sb5 which influences the superconducting coupling in the junction, possibly giving rise to spin-triplet superconductivity. In addition, the observation of long-lived fast oscillations shows evidence of spatially localized conducting channels arising from edge states. These observations pave the way for studying unconventional superconductivity and Josephson device based on Kagome metals with electron correlation and topology.


title: Are Fewer Labels Possible for Few-shot Learning?

abstract: Few-shot learning is challenging due to its very limited data and labels. Recent studies in big transfer (BiT) show that few-shot learning can greatly benefit from pretraining on large scale labeled dataset in a different domain. This paper asks a more challenging question: "can we use as few as possible labels for few-shot learning in both pretraining (with no labels) and fine-tuning (with fewer labels)?". Our key insight is that the clustering of target samples in the feature space is all we need for few-shot finetuning. It explains why the vanilla unsupervised pretraining (poor clustering) is worse than the supervised one. In this paper, we propose transductive unsupervised pretraining that achieves a better clustering by involving target data even though its amount is very limited. The improved clustering result is of great value for identifying the most representative samples ("eigen-samples") for users to label, and in return, continued finetuning with the labeled eigen-samples further improves the clustering. Thus, we propose eigen-finetuning to enable fewer shot learning by leveraging the co-evolution of clustering and eigen-samples in the finetuning. We conduct experiments on 10 different few-shot target datasets, and our average few-shot performance outperforms both vanilla inductive unsupervised transfer and supervised transfer by a large margin. For instance, when each target category only has 10 labeled samples, the mean accuracy gain over the above two baselines is 9.2% and 3.42 respectively.


title: Binning is Sinning (Supernova Version): The Impact of Self-Calibration in Cosmological Analyses with Type Ia Supernovae

abstract: Recent cosmological analyses (e.g., JLA, Pantheon) of Type Ia Supernova (SNIa) have propagated systematic uncertainties into a covariance matrix and either binned or smoothed the systematic vectors in redshift space. We demonstrate that systematic error budgets of these analyses can be improved by a factor of $\sim1.5\times$ with the use of unbinned and unsmoothed covariance matrices. To understand this, we employ a separate approach that simultaneously fits for cosmological parameters and additional self-calibrating scale parameters that constrain the size of each systematic. We show that the covariance-matrix approach and scale-parameter approach yield equivalent results, implying that in both cases the data can self-calibrate certain systematic uncertainties, but that this ability is hindered when information is binned or smoothed in redshift space. We review the top systematic uncertainties in current analyses and find that the reduction of systematic uncertainties in the unbinned case depends on whether a systematic is consistent with varying the cosmological model and whether or not the systematic can be described by additional correlations between SN properties and luminosity. Furthermore, we show that the power of self-calibration increases with the size of the dataset, which presents a tremendous opportunity for upcoming analyses of photometrically classified samples, like those of Legacy Survey of Space and Time (LSST) and the Nancy Grace Roman Telescope (NGRST). However, to take advantage of self-calibration in large, photometrically-classified samples, we must first address the issue that binning is required in currently-used photometric methodologies.


title: Robust Consistent Video Depth Estimation

abstract: We present an algorithm for estimating consistent dense depth maps and camera poses from a monocular video. We integrate a learning-based depth prior, in the form of a convolutional neural network trained for single-image depth estimation, with geometric optimization, to estimate a smooth camera trajectory as well as detailed and stable depth reconstruction. Our algorithm combines two complementary techniques: (1) flexible deformation-splines for low-frequency large-scale alignment and (2) geometry-aware depth filtering for high-frequency alignment of fine depth details. In contrast to prior approaches, our method does not require camera poses as input and achieves robust reconstruction for challenging hand-held cell phone captures containing a significant amount of noise, shake, motion blur, and rolling shutter deformations. Our method quantitatively outperforms state-of-the-arts on the Sintel benchmark for both depth and pose estimations and attains favorable qualitative results across diverse wild datasets.


title: Ultrafast reorientation of the N\'eel vector in antiferromagnetic Dirac semimetals

abstract: Antiferromagnets exhibit distinctive characteristics such as ultrafast dynamics and robustness against perturbative fields, thereby attracting considerable interest in fundamental physics and technological applications. Recently, it was revealed that the N\'eel vector can be switched by a current-induced staggered (N\'eel) spin-orbit torque in antiferromagnets with the parity-time symmetry, and furthermore, a nonsymmorphic symmetry enables the control of Dirac fermions. However, the real-time dynamics of the magnetic and electronic structures remain largely unexplored. Here, we propose a theory of the ultrafast dynamics in antiferromagnetic Dirac semimetals and show that the N\'eel vector is rotated in the picosecond timescale by the terahertz-pulse-induced N\'eel spin-orbit torque and other torques originating from magnetic anisotropies. This reorientation accompanies the modulation of the mass of Dirac fermions and can be observed in real time by the magneto-optical effects. Our results provide a theoretical basis for emerging ultrafast antiferromagnetic spintronics combined with the topological aspects of materials.


title: Portrait Neural Radiance Fields from a Single Image

abstract: We present a method for estimating Neural Radiance Fields (NeRF) from a single headshot portrait. While NeRF has demonstrated high-quality view synthesis, it requires multiple images of static scenes and thus impractical for casual captures and moving subjects. In this work, we propose to pretrain the weights of a multilayer perceptron (MLP), which implicitly models the volumetric density and colors, with a meta-learning framework using a light stage portrait dataset. To improve the generalization to unseen faces, we train the MLP in the canonical coordinate space approximated by 3D face morphable models. We quantitatively evaluate the method using controlled captures and demonstrate the generalization to real portrait images, showing favorable results against state-of-the-arts.


title: The product on $\W$-spaces of rational forms

abstract: We explore the notion of the spaces $\W_{z_1, \ldots, z_n}$ of rational differential forms with complex formal parameteres $(z_1, \ldots, z_n)$ for $n \ge 0$, and define a product between theire elements. Let $V$ be a quasi-conformal grading-restricted vertex algebra, $W$ be its module, $\overline{W}$ be the algebraic completion of $W$, and $\W_{z_1, \ldots, z_n}$ be the space of rational differential forms in $(z_1, \ldots, z_n)$. Using geometric interpretation in terms of sewing two Riemann spheres with a number of marked points, we introduce a product of elements of two spaces $\W_{x_1, \ldots, x_k}$ and $\W_{y_1, \ldots, y_n}$, and study its properties. The product takes values in $\W_{x_1, \ldots, x_k; y_1, \ldots, y_n}$. We prove that the product is defined by an absolutely convergent series. In applications, for two spaces $C^k_m(V, \W)$ and $C^n_{m'}(V, \W)$ (introduced in \cite{Huang}) of chain-cochain double complex associated to a grading-restricted vertex algerba $V$ (which provides an example of $\W$ introduced in \cite{Huang}) we define a product between them coherent with the differential of the complex. We prove that the product brings about a map to the space $C^{k+n}_{m+m'}(V, \W)$, and satisfy an analogue of Leibniz formula.


title: Synergistic Integration of Optical and Microwave Satellite Data for Crop Yield Estimation

abstract: Developing accurate models of crop stress, phenology and productivity is of paramount importance, given the increasing need of food. Earth observation remote sensing data provides a unique source of information to monitor crops in a temporally resolved and spatially explicit way. In this study, we propose the combination of multisensor (optical and microwave) remote sensing data for crop yield estimation and forecasting using two novel approaches. We first propose the lag between Enhanced Vegetation Index derived from MODIS and Vegetation Optical Depth derived from SMAP as a new joint metric combining the information from the two satellite sensors in a unique feature or descriptor. Our second approach avoids summarizing statistics and uses machine learning to combine full time series of EVI and VOD. This study considers two statistical methods, a regularized linear regression and its nonlinear extension called kernel ridge regression to directly estimate the county-level surveyed total production, as well as individual yields of the major crops grown in the region: corn, soybean and wheat. The study area includes the US Corn Belt, and we use agricultural survey data from the National Agricultural Statistics Service (USDA-NASS) for year 2015 for quantitative assessment.


title: A Sentiment Analysis Approach to the Prediction of Market Volatility

abstract: Prediction and quantification of future volatility and returns play an important role in financial modelling, both in portfolio optimization and risk management. Natural language processing today allows to process news and social media comments to detect signals of investors' confidence. We have explored the relationship between sentiment extracted from financial news and tweets and FTSE100 movements. We investigated the strength of the correlation between sentiment measures on a given day and market volatility and returns observed the next day. The findings suggest that there is evidence of correlation between sentiment and stock market movements: the sentiment captured from news headlines could be used as a signal to predict market returns; the same does not apply for volatility. Also, in a surprising finding, for the sentiment found in Twitter comments we obtained a correlation coefficient of -0.7, and p-value below 0.05, which indicates a strong negative correlation between positive sentiment captured from the tweets on a given day and the volatility observed the next day. We developed an accurate classifier for the prediction of market volatility in response to the arrival of new information by deploying topic modelling, based on Latent Dirichlet Allocation, to extract feature vectors from a collection of tweets and financial news. The obtained features were used as additional input to the classifier. Thanks to the combination of sentiment and topic modelling our classifier achieved a directional prediction accuracy for volatility of 63%.


title: Data-driven Method for Estimating Aircraft Mass from Quick Access Recorder using Aircraft Dynamics and Multilayer Perceptron Neural Network

abstract: Accurate aircraft-mass estimation is critical to airlines from the safety-management and performance-optimization viewpoints. Overloading an aircraft with passengers and baggage might result in a safety hazard. In contrast, not fully utilizing an aircraft's payload-carrying capacity undermines its operational efficiency and airline profitability. However, accurate determination of the aircraft mass for each operating flight is not feasible because it is impractical to weigh each aircraft component, including the payload. The existing methods for aircraft-mass estimation are dependent on the aircraft- and engine-performance parameters, which are usually considered proprietary information. Moreover, the values of these parameters vary under different operating conditions while those of others might be subject to large estimation errors. This paper presents a data-driven method involving use of the quick access recorder (QAR)-a digital flight-data recorder-installed on all aircrafts to record the initial aircraft climb mass during each flight. The method requires users to select appropriate parameters among several thousand others recorded by the QAR using physical models. The selected data are subsequently processed and provided as input to a multilayer perceptron neural network for building the model for initial-climb aircraft-mass prediction. Thus, the proposed method offers the advantages of both the model-based and data-driven approaches for aircraft-mass estimation. Because this method does not explicitly rely on any aircraft or engine parameter, it is universally applicable to all aircraft types. In this study, the proposed method was applied to a set of Boeing 777-300ER aircrafts, the results of which demonstrated reasonable accuracy. Airlines can use this tool to better utilize aircraft's payload.


title: Ensemble of Discriminators for Domain Adaptation in Multiple Sound Source 2D Localization

abstract: This paper introduces an ensemble of discriminators that improves the accuracy of a domain adaptation technique for the localization of multiple sound sources. Recently, deep neural networks have led to promising results for this task, yet they require a large amount of labeled data for training. Recording and labeling such datasets is very costly, especially because data needs to be diverse enough to cover different acoustic conditions. In this paper, we leverage acoustic simulators to inexpensively generate labeled training samples. However, models trained on synthetic data tend to perform poorly with real-world recordings due to the domain mismatch. For this, we explore two domain adaptation methods using adversarial learning for sound source localization which use labeled synthetic data and unlabeled real data. We propose a novel ensemble approach that combines discriminators applied at different feature levels of the localization model. Experiments show that our ensemble discrimination method significantly improves the localization performance without requiring any label from the real data.


title: Blending MPC & Value Function Approximation for Efficient Reinforcement Learning

abstract: Model-Predictive Control (MPC) is a powerful tool for controlling complex, real-world systems that uses a model to make predictions about future behavior. For each state encountered, MPC solves an online optimization problem to choose a control action that will minimize future cost. This is a surprisingly effective strategy, but real-time performance requirements warrant the use of simple models. If the model is not sufficiently accurate, then the resulting controller can be biased, limiting performance. We present a framework for improving on MPC with model-free reinforcement learning (RL). The key insight is to view MPC as constructing a series of local Q-function approximations. We show that by using a parameter $\lambda$, similar to the trace decay parameter in TD($\lambda$), we can systematically trade-off learned value estimates against the local Q-function approximations. We present a theoretical analysis that shows how error from inaccurate models in MPC and value function estimation in RL can be balanced. We further propose an algorithm that changes $\lambda$ over time to reduce the dependence on MPC as our estimates of the value function improve, and test the efficacy our approach on challenging high-dimensional manipulation tasks with biased models in simulation. We demonstrate that our approach can obtain performance comparable with MPC with access to true dynamics even under severe model bias and is more sample efficient as compared to model-free RL.


title: Von Neumann entropy in a dispersive cavity

abstract: We study the von Neumann entropy of the partial trace of a system of two two-level atoms (qubits) in a dispersive cavity where the atoms are interacting collectively with a single mode electromagnetic field in the cavity. We make a contrast of this entanglement entropy with the spin squeezing property of the system. We find a close relationship between this von Neumann entropy and the spin squeezing of the system. We observe that this entanglement entropy attains its maximum value when the spin squeezing in the system goes to its maximum and is zero when the spin squeezing in the system vanishes.


title: The Zeno and anti-Zeno effects: studying modified decay rates for spin-boson models with both strong and weak system-environment couplings

abstract: In this paper, we look into what happens to a quantum system under repeated measurements if system evolution is removed before each measurement is performed. Beginning with investigating a single two-level system coupled to two independent baths of harmonic oscillators, we move to replacing it with a large collection of such systems, thereby invoking the large spin-boson model. Whereas each of our two-level systems interacts strongly with one of the aforementioned baths, it interacts weakly with the other. A polaron transformation is used to make it possible for the problem in the strong coupling regime to be treated with perturbation theory. We find that the case involving a single two-level system exhibits qualitative and quantitative differences from the case involving a collection of them; however, the general effects of strong and weak couplings turn out to be the same as those in the presence of system evolution, something which allows us to establish that system evolution has no practical bearing on any of these effects.


title: Gap Filling of Biophysical Parameter Time Series with Multi-Output Gaussian Processes

abstract: In this work we evaluate multi-output (MO) Gaussian Process (GP) models based on the linear model of coregionalization (LMC) for estimation of biophysical parameter variables under a gap filling setup. In particular, we focus on LAI and fAPAR over rice areas. We show how this problem cannot be solved with standard single-output (SO) GP models, and how the proposed MO-GP models are able to successfully predict these variables even in high missing data regimes, by implicitly performing an across-domain information transfer.


title: A dearth of young and bright massive stars in the Small Magellanic Cloud

abstract: Massive star evolution at low metallicity is closely connected to many fields in high-redshift astrophysics, but poorly understood. The Small Magellanic Cloud (SMC) is a unique laboratory to study this because of its metallicity of 0.2 Zsol, its proximity, and because it is currently forming stars. We used a spectral type catalog in combination with GAIA magnitudes to calculate temperatures and luminosities of bright SMC stars. By comparing these with literature studies, we tested the validity of our method, and using GAIA data, we estimated the completeness of stars in the catalog as a function of luminosity. This allowed us to obtain a nearly complete view of the most luminous stars in the SMC. When then compared with stellar evolution predictions. We also calculated the extinction distribution, the ionizing photon production rate, and the star formation rate. Our results imply that the SMS hosts only 30 very luminous main-sequence stars (M > 40 Msol; L > 10^5 Lsol), which are far fewer than expected from the number of stars in the luminosity range 3*10^4 < L/Lsol < 3*10^5 and from the typically quoted star formation rate in the SMC. Even more striking, we find that for masses above M > 20 Msol, stars in the first half of their hydrogen-burning phase are almost absent. This mirrors a qualitatively similar peculiarity that is known for the Milky Way and Large Magellanic Cloud. This amounts to a lack of hydrogen-burning counterparts of helium-burning stars, which is more pronounced for higher luminosities. We argue that a declining star formation rate or a steep initial mass function are unlikely to be the sole explanations for the dearth of young bright stars. Instead, many of these stars might be embedded in their birth clouds, although observational evidence for this is weak. We discuss implications for cosmic reionization and the top end of the initial mass function.


title: The link between star formation and gas in nearby galaxies

abstract: Observations of the interstellar medium are key to deciphering the physical processes regulating star formation in galaxies. However, observational uncertainties and detection limits can bias the interpretation unless carefully modeled. Here I re-analyze star formation rates and gas masses of a representative sample of nearby galaxies with the help of multi-dimensional Bayesian modeling. Typical star forming galaxies are found to lie in a 'star forming plane' largely independent of their stellar mass. Their star formation activity is tightly correlated with the molecular and total gas content, while variations of the molecular-gas-to-star conversion efficiency are shown to be significantly smaller than previously reported. These data-driven findings suggest that physical processes that modify the overall galactic gas content, such as gas accretion and outflows, regulate the star formation activity in typical nearby galaxies, while a change in efficiency triggered by, e.g., galaxy mergers or gas instabilities, may boost the activity of starbursts.


title: Symmetry constraints on superconductivity in twisted bilayer graphene: Fractional vortices, $4e$ condensates or non-unitary pairing

abstract: When two graphene sheets are twisted relative to each other by a small angle, enhanced correlations lead to superconductivity whose origin remains under debate. Here, we derive some general constraints on superconductivity in twisted bilayer graphene (TBG), independent of its underlying mechanism. Neglecting weak coupling between valleys, the global symmetry group of TBG consists of independent spin rotations in each valley in addition to valley charge rotations, $ {\rm SU}(2) \times {\rm SU}(2) \times {\rm U}_V(1) $. This symmetry is further enhanced to a full ${\rm SU}(4)$ in the idealized chiral limit. In both cases, we show that any charge $2e$ pairing must break the global symmetry. Additionally, if the pairing is unitary the resulting superconductor admits fractional vortices. This leads to the following general statement: Any superconducting condensate in either symmetry class has to satisfy one of three possibilities: (i) the superconducting pairing is non-unitary, (ii) the superconducting condensate has charge $2e$ but admits at least half quantum vortices which carry a flux of $h/4e$, or (iii) the superconducting condensate has charge $2me$, $m>1$, with vortices carrying $h/2me$ flux. The latter possibility can be realized by a symmetric charge $4e$ superconductor ($m=2$). Non-unitary pairing (i) is expected in TBG for superconductors observed in the vicinity of flavor polarized states. On the other hand, in the absence of flavor polarization, e.g. in the vicinity of charge neutrality, one of the two exotic possibilities (ii) and (iii) is expected. We sketch how all three scenarios can be realized in different limits within a strong coupling theory of superconductivity based on skyrmions. Finally we discuss the effect of symmetry lowering anisotropies and experimental implications of these scenarios.


title: Crossover between Strongly-coupled and Weakly-coupled Exciton Superfluids

abstract: In fermionic systems, superconductivity and superfluidity are enabled through the condensation of fermion pairs. The nature of this condensate can be tuned by varying the pairing strength, with weak coupling yielding a BCS-like condensate and strong coupling resulting in a BEC-like process. However, demonstration of this cross-over has remained elusive in electronic systems. Here we study graphene double-layers separated by an atomically thin insulator. Under applied magnetic field, electrons and holes couple across the barrier to form bound magneto-excitons whose pairing strength can be continuously tuned by varying the effective layer separation. Using temperature-dependent Coulomb drag and counter-flow current measurements, we demonstrate the capability to tune the magneto-exciton condensate through the entire weak-coupling to strong-coupling phase diagram. Our results establish magneto-exciton condensates in graphene as a model platform to study the crossover between two Bosonic quantum condensate phases in a solid state system.


title: Quantum Spin Hall Effect in Ta$_2$M$_3$Te$_5$ (M = Pd, Ni)

abstract: Quantum spin Hall (QSH) effect with great promise for the potential application in spintronics and quantum computing has attracted extensive research interest from both theoretical and experimental researchers. Here, we predict monolayer Ta$_2$Pd$_3$Te$_5$ can be a QSH insulator based on first-principles calculations. The interlayer binding energy in the layered van der Waals compound Ta$_2$Pd$_3$Te$_5$ is 19.6 meV/A$^2$; thus, its monolayer/thin-film structures could be readily obtained by exfoliation. The band inversion near the Fermi level ($E_F$) is an intrinsic characteristic, which happens between Ta-$5d$ and Pd-$4d$ orbitals without spin-orbit coupling (SOC). The SOC effect opens a global gap and makes the system a QSH insulator. With the $d$-$d$ band-inverted feature, the nontrivial topology in monolayer Ta$_2$Pd$_3$Te$_5$ is characterized by the time-reversal topological invariant $\mathbb Z_2=1$, which is computed by the one-dimensional (1D) Wilson loop method as implemented in our first-principles calculations. The helical edge modes are also obtained using surface Green's function method. Our calculations show that the QSH state in Ta$_2M_3$Te$_5$ ($M=$ Pd, Ni) can be tuned by external strain. These monolayers and thin films provide feasible platforms for realizing QSH effect as well as related devices.


title: Mixed EW-QCD two-loop amplitudes for $q\bar{q} \to \ell^+\ell^-$ and $\gamma_5$ scheme independence of multi-loop corrections

abstract: We perform a dedicated study of the $q \bar{q}$-initiated two-loop electroweak-QCD Drell-Yan scattering amplitude in dimensional regularization schemes for vanishing light quark and lepton masses. For the relative order $\alpha$ and $\alpha_s$ one-loop Standard Model corrections, details of our comparison to the original literature are given. The infrared pole terms of the mixed two-loop amplitude are governed by a known generalization of the dipole formula and we show explicitly that exactly the same two-loop polarized hard scattering functions are obtained in both the standard 't Hooft-Veltman-Breitenlohner-Maison $\gamma_5$ scheme and Kreimer's anticommuting $\gamma_5$ scheme.


title: Simulations of the star-forming molecular gas in an interacting M51-like galaxy: cloud population statistics

abstract: To investigate how molecular clouds react to different environmental conditions at a galactic scale, we present a catalogue of giant molecular clouds resolved down to masses of $\sim 10$~M$_{\odot}$ from a simulation of the entire disc of an interacting M51-like galaxy and a comparable isolated galaxy. Our model includes time-dependent gas chemistry, sink particles for star formation and supernova feedback, meaning we are not reliant on star formation recipes based on threshold densities and can follow the physics of the cold molecular phase. We extract giant molecular clouds at a given timestep of the simulations and analyse their properties. In the disc of our simulated galaxies, spiral arms seem to act merely as snowplows, gathering gas and clouds without dramatically affecting their properties. In the centre of the galaxy, on the other hand, environmental conditions lead to larger, more massive clouds. While the galaxy interaction has little effect on cloud masses and sizes, it does promote the formation of counter-rotating clouds. We find that the identified clouds seem to be largely gravitationally unbound at first glance, but a closer analysis of the hierarchical structure of the molecular interstellar medium shows that there is a large range of virial parameters with a smooth transition from unbound to mostly bound for the densest structures. The common observation that clouds appear to be virialised entities may therefore be due to CO bright emission highlighting a specific level in this hierarchical binding sequence. The small fraction of gravitationally bound structures found suggests that low galactic star formation efficiencies may be set by the process of cloud formation and initial collapse.


title: Pulse envelope effect on nonlinear Compton scattering in electron-laser collisions

abstract: Nonlinear Compton scattering is calculated for the collision of an electron with a plane wave pulse. A mid infra-red (IR) peak arises in the photon spectrum due to long-range interference associated with the pulse envelope. The case of a flat-top pulse is studied as a toy model for pulse envelope effects and reduced to two final-state momentum integrations; the case of a sine-squared pulse is studied numerically. A perturbative expansion in charge-field coupling reveals that already at intermediate intensities, many orders are required to correctly capture the structure of the mid-IR peak. By regularising the classical result, it is shown that the mid-IR peak is due to plane-wave ponderomotive effects on the pulse envelope. Finally, it is shown that the mid-IR peak can be isolated using energy, angle and polarisation filters.


title: Pion crystals hosting topologically stable baryons

abstract: We construct analytic (3+1)-dimensional inhomogeneous and topologically non-trivial pion systems using chiral perturbation theory. We discuss the effect of isospin asymmetry with vanishing electromagnetic interactions as well as some particular configurations with non-vanishing electromagnetic interactions. The inhomogeneous configurations of the pion fields are characterized by a non-vanishing topological charge that can be identified with baryons surrounded by a cloud of pions. This system supports a topologically protected persistent superflow. When the electromagnetic field is turned on the superflow corresponds to an electromagnetic supercurrent.


title: Stellar Winds Drive Strong Variations in Exoplanet Evaporative Outflows and Transit Absorption Signatures

abstract: Stellar wind and photon radiation interactions with a planet can cause atmospheric depletion, which may have a potentially catastrophic impact on a planet's habitability. While the implications of photoevaporation on atmospheric erosion have been researched to some degree, studies of the influence of the stellar wind on atmospheric loss are in their infancy. Here, we use three-dimensional magnetohydrodynamic simulations to model the effect of the stellar wind on the magnetosphere and outflow of a hypothetical planet, modeled to have an H-rich evaporating envelope with a pre-defined mass loss rate, orbiting in the habitable zone close to a low-mass M dwarf. We take the TRAPPIST-1 system as a prototype, with our simulated planet situated at the orbit of TRAPPIST-1e. We show that the atmospheric outflow is dragged and accelerated upon interaction with the wind, resulting in a diverse range of planetary magnetosphere morphologies and plasma distributions as local stellar wind conditions change. We consider the implications of the wind-outflow interaction on potential hydrogen Lyman-alpha (Lya) observations of the planetary atmosphere during transits. The Lya observational signatures depend strongly on the local wind conditions at the time of the observation and can be subject to considerable variation on timescales as short as an hour. Our results indicate that observed variations in exoplanet Lya transit signatures could be explained by wind-outflow interaction.


title: Transmon platform for quantum computing challenged by chaotic fluctuations

abstract: From the perspective of many body physics, the transmon qubit architectures currently developed for quantum computing are systems of coupled nonlinear quantum resonators. A significant amount of intentional frequency detuning (disorder) is required to protect individual qubit states against the destabilizing effects of nonlinear resonator coupling. Here we investigate the stability of this variant of a many-body localized (MBL) phase for system parameters relevant to current quantum processors of two different types, those using untunable qubits (IBM type) and those using tunable qubits (Delft/Google type). Applying three independent diagnostics of localization theory -- a Kullback-Leibler analysis of spectral statistics, statistics of many-body wave functions (inverse participation ratios), and a Walsh transform of the many-body spectrum -- we find that these computing platforms are dangerously close to a phase of uncontrollable chaotic fluctuations.


title: Exact analytical treatment of multiqubit noisy dynamics in exchange-coupled semiconductor spin qubits

abstract: Charge noise remains the primary obstacle to the development of quantum information technologies with semiconductor spin qubits. We use an exact analytical calculation to determine the effects of quasistatic charge noise on a ring of three equally spaced exchange-coupled quantum dots. We calculate the disorder-averaged return probability from a specific initial state, and use it to determine the coherence time T2* and show that it depends only on the disorder strength and not the mean interaction strength. We also use a perturbative approach to investigate other arrangements of three or four qubits, finding that the return probability contains multiple oscillation frequencies. These oscillations decay in a Gaussian manner, determined by differences in energy levels of the Hamiltonian. We give quantitative values for gate times resulting in several target fidelities. We find that the decoherence time decreases with increasing number of qubits. Our work provides useful analytical insight into the charge noise dynamics of coupled spin qubits.


title: Testing $f(R)$ Gravity With Scale Dependent Cosmic Void Velocity Profiles

abstract: We study the impact of cosmological scale modifications to General Relativity on the dynamics of halos within voids by comparing N-body simulations incorporating Hu-Sawicki $f(R)$ gravity, with $|f_{R0}|=10^{-6}$ and $10^{-5}$, to those of $\Lambda$CDM. By examining the radial velocity statistics within voids classified based on their size and density-profile, as "rising" ($R$-type) or "shell" ($S$-type), we find that halo motions in small $R$-type voids, with effective radius $<15 Mpc/h$, reveal distinctive differences between $f(R)$ and $\Lambda$CDM cosmologies. To understand this observed effect, we study the linear and nonlinear fifth forces, and develop an iterative algorithm to accurately solve the non-linear fifth force equation. We use this to characterize the Chameleon screening mechanism in voids and contrast the behavior with that observed in gravitationally collapsed objects. The force analysis underscores how smaller $R$-type voids exhibit the highest ratios of fifth force to Newtonian force, which source distinguishable differences in the velocity profiles and thereby provide rich environments in which to constrain gravity.


title: Redshift space three-point correlation function of IGM at $z<0.48$

abstract: The Ly$\alpha$ forest decomposed into Voigt profile components allow us to study clustering properties of the intergalactic medium and its dependence on various physical quantities. Here, we report the first detections of probability excess of low-z (i.e z<0.48) Ly$\alpha$ absorber triplets over a scale of $r_\parallel\leq8$pMpc with a maximum amplitude of $8.76^{+1.96}_{-1.65}$ at a longitudinal separation of 1-2pMpc. We measure non-zero three-point correlation ($\zeta=4.76^{+1.98}_{-1.67}$) only at this scale with reduced three-point correlation value of Q=$0.95^{+0.39}_{-0.38}$. The measured $\zeta$ shows an increasing trend with increasing HI column density ($N_{HI}$) while Q does not show any $N_{HI}$ dependence. About 88% of the triplets contributing to $\zeta$ (at z$\le0.2$) have nearby galaxies (whose distribution is known to be complete for 0.1L$_*$ at z<0.1 and for L$_*$ at z~0.25 and within 20' to the quasar sightlines) within a velocity separation of 500$kms^{-1}$ and a median impact parameter of 405pkpc. The measured impact parameters are consistent with appreciable number of triplets at z$\le0.2$ not originating from individual galaxies but tracing the underlying galaxy distribution. Frequency of occurrence of high-b absorbers in triplets (~85%) is a factor~3 higher than that found among the full sample (~32%). Using four different cosmological simulations, we quantify the effect of peculiar velocities, feedback effects and show that most of the observed trends are broadly reproduced. However, $\zeta$ at small scales ($r_\parallel<1$pMpc) and b-dependence of $\zeta$ in simulations are found inconsistent with the observations. This could either be related to the fact that none of these simulations reproduce the observed b-distribution and $N_{HI}$ distribution for $N_{HI}>10^{14}$cm$^{-2}$ self-consistently or to the widespread of signal-to-noise ratio in the observed data.


title: The TW Hya Rosetta Stone Project III: Resolving the Gaseous Thermal Profile of the Disk

abstract: The thermal structure of protoplanetary disks is a fundamental characteristic of the system that has wide reaching effects on disk evolution and planet formation. In this study, we constrain the 2D thermal structure of the protoplanetary disk TW Hya structure utilizing images of seven CO lines. This includes new ALMA observations of 12CO J=2-1 and C18O J=2-1 as well as archival ALMA observations of 12CO J=3-2, 13CO J=3-2, 6-5, C18O J= 3-2, 6-5. Additionally, we reproduce a Herschel observation of the HD J=1-0 line flux, the spectral energy distribution, and utilize a recent quantification of CO radial depletion in TW Hya. These observations were modeled using the thermochemical code RAC2D, and our best fit model reproduces all spatially resolved CO surface brightness profiles. The resulting thermal profile finds a disk mass of 0.025 Msun and a thin upper layer of gas depleted of small dust with a thickness of approx 1.2% of the corresponding radius. Using our final thermal structure, we find that CO alone is not a viable mass tracer as its abundance is degenerate with the total H2 surface density. Different mass models can readily match the spatially resolved CO line profiles with disparate abundance assumptions. Mass determination requires additional knowledge and, in this work, HD provides the additional constraint to derive the gas mass and supports the inference of CO depletion in the TW Hya disk. Our final thermal structure confirms the use of HD as a powerful probe of protoplanetary disk mass. Additionally, the method laid out in this paper is an employable strategy for extraction of disk temperatures and masses in the future.


title: A machine learning approach to galaxy properties: joint redshift-stellar mass probability distributions with Random Forest

abstract: We demonstrate that highly accurate joint redshift-stellar mass probability distribution functions (PDFs) can be obtained using the Random Forest (RF) machine learning (ML) algorithm, even with few photometric bands available. As an example, we use the Dark Energy Survey (DES), combined with the COSMOS2015 catalogue for redshifts and stellar masses. We build two ML models: one containing deep photometry in the $griz$ bands, and the second reflecting the photometric scatter present in the main DES survey, with carefully constructed representative training data in each case. We validate our joint PDFs for $10,699$ test galaxies by utilizing the copula probability integral transform and the Kendall distribution function, and their univariate counterparts to validate the marginals. Benchmarked against a basic set-up of the template-fitting code BAGPIPES, our ML-based method outperforms template fitting on all of our predefined performance metrics. In addition to accuracy, the RF is extremely fast, able to compute joint PDFs for a million galaxies in just under $6$ min with consumer computer hardware. Such speed enables PDFs to be derived in real time within analysis codes, solving potential storage issues. As part of this work we have developed GALPRO, a highly intuitive and efficient Python package to rapidly generate multivariate PDFs on-the-fly. GALPRO is documented and available for researchers to use in their cosmology and galaxy evolution studies.

categories: astro-ph.GA astro-ph.CO astro-ph.IM cs.LG


title: An Algorithm for the Separation-Preserving Transition of Clusterings

abstract: The separability of clusters is one of the most desired properties in clustering. There is a wide range of settings in which different clusterings of the same data set appear. We are interested in applications where there is a need for an explicit, gradual transition of one separable clustering into another one. This transition should be a sequence of simple, natural steps that upholds separability of the clusters throughout. We design an algorithm for such a transition. We exploit the intimate connection of separability and linear programming over bounded-shape partition and transportation polytopes: separable clusterings lie on the boundary of partition polytopes, form a subset of the vertices of the corresponding transportation polytopes, and circuits of both polytopes are readily interpreted as sequential or cyclical exchanges of items between clusters. This allows for a natural approach to achieve the desired transition through a combination of two walks: an edge walk between two so-called radial clusterings in a transportation polytope, computed through an adaptation of classical tools of sensitivity analysis and parametric programming; and a walk from a separable clustering to a corresponding radial clustering, computed through a tailored, iterative routine updating cluster sizes and re-optimizing the cluster assignment of items.


title: Colorful points in the XY regime of XXZ quantum magnets

abstract: In the XY regime of the XXZ Heisenberg model phase diagram, we demonstrate that the origin of magnetically ordered phases is influenced by the presence of solvable points with exact quantum coloring ground states featuring a quantum-classical correspondence. Using exact diagonalization and density matrix renormalization group calculations, for both the square and the triangular lattice magnets, we show that the ordered physics of the solvable points in the extreme XY regime, at $\frac{J_z}{J_\perp}=-1$ and $\frac{J_z}{J_\perp}=-\frac{1}{2}$ respectively with $J_\perp > 0$, adiabatically extends to the more isotropic regime $\frac{J_z}{J_\perp} \sim 1$. We highlight the projective structure of the coloring ground states to compute the correlators in fixed magnetization sectors which enables an understanding of the features in the static spin structure factors and correlation ratios. These findings are contrasted with an anisotropic generalization of the celebrated one-dimensional Majumdar-Ghosh model, which is also found to be (ground state) solvable. For this model, both exact dimer and three-coloring ground states exist at $\frac{J_z}{J_\perp}=-\frac{1}{2}$ but only the two dimer ground states survive for any $\frac{J_z}{J_\perp} > -\frac{1}{2}$.


title: The Cepheid Distance to the Narrow-Line Seyfert 1 Galaxy NGC 4051

abstract: We derive a distance of $D = 16.6 \pm 0.3$~Mpc ($\mu=31.10\pm0.04$~mag) to the archetypal narrow-line Seyfert 1 galaxy NGC 4051 based on Cepheid Period--Luminosity relations and new Hubble Space Telescope multiband imaging. We identify 419 Cepheid candidates and estimate the distance at both optical and near-infrared wavelengths using subsamples of precisely-photometered variables (123 and 47 in the optical and near-infrared subsamples, respectively). We compare our independent photometric procedures and distance-estimation methods to those used by the SH0ES team and find agreement to 0.01~mag. The distance we obtain suggests an Eddington ratio $\dot{m} \approx 0.2$ for NGC 4051, typical of narrow-line Seyfert 1 galaxies, unlike the seemingly-odd value implied by previous distance estimates. We derive a peculiar velocity of $-490\pm34$~km~s$^{-1}$ for NGC 4051, consistent with the overall motion of the Ursa Major Cluster in which it resides. We also revisit the energetics of the NGC 4051 nucleus, including its outflow and mass accretion rates.


title: Investigating $\gamma$-ray halos around three HAWC bright sources in Fermi-LAT data

abstract: Numerous extended sources around Galactic pulsars have shown significant $\gamma$-ray emission from GeV to TeV energies, revealing hundreds of TeV energy electrons scattering off of the underlying photon fields through inverse Compton scattering (ICS). HAWC TeV gamma-ray observations of few-degree extended emission around the pulsars Geminga and Monogem, and LAT GeV emission around Geminga, suggest that systems older than 100 kyr have multi-TeV $e^\pm$ propagating beyond the SNR-PWN system into the interstellar medium. Following the discovery of few $\gamma$-ray sources by HAWC at energies E$>100$~TeV, we investigate the presence of an extended $\gamma$-ray emission in Fermi-LAT data around the three brightest sources detected by HAWC up to 100 TeV. We find an extended emission of $\theta_{68} = 1.00^{+0.05}_{-0.07}$ deg around eHWC J1825-134 and $\theta_{68} = 0.71\pm0.10$ deg eHWC J1907+063. The analysis with ICS templates on Fermi-LAT data point to diffusion coefficient values which are significantly lower than the average Galactic one. When studied along with HAWC data, the $\gamma$-ray Fermi-LAT data provide invaluable insight into the very high-energy electron and positron parent populations.