title: A new 3D solar wind speed and density model based on interplanetary scintillation

abstract: The solar wind (SW) is an outflow of the solar coronal plasma, expanding supersonically throughout the heliosphere. SW particles interact by charge exchange with interstellar neutral atoms and on one hand, they modify the distribution of this gas in interplanetary space, and on the other hand they are seed population for heliospheric pickup ions and energetic neutral atoms (ENAs). The heliolatitudinal profiles of the SW speed and density evolve during the cycle of solar activity. A model of evolution of the SW speed and density is needed to interpret observations of ENAs, pickup ions, the heliospheric backscatter glow, etc. We derive the Warsaw Heliospheric Ionization Model 3DSW (WawHelIon 3DSW) based on interplanetary scintillation (IPS) tomography maps of the SW speed. We take the IPS tomography data from 1985 until 2020, compiled by \citet{tokumaru_etal:21a}. We derive a novel statistical method of filtering these data against outliers, we present a flexible analytic formula for the latitudinal profiles of the SW speed based on Legendre polynomials of varying order with additional restraining conditions at the poles, fit this formula to the yearly filtered data, and calculate the yearly SW density profiles using the latitudinally invariant SW energy flux, observed in the ecliptic plane. Despite application of refined IPS data set, a more sophisticated data filtering method, and a more flexible analytic model, the present results mostly agree with those obtained previously, which demonstrates the robustness of IPS studies of the SW structure.


title: A duality of scaffolds for translation association schemes

abstract: Scaffolds are certain tensors arising in the study of association schemes, and have been (implicitly) understood diagrammatically as digraphs with distinguished "root" nodes and with matrix edge weights, often taken from Bose-Mesner algebras. In this paper, we first present a slight modification of Martin's conjecture (2021) concerning a duality of scaffolds whose digraphs are embedded in a closed disk in the plane with root nodes all lying on the boundary circle, and then show that this modified conjecture holds true if we restrict ourselves to the class of translation association schemes, i.e., those association schemes that admit abelian regular automorphism groups.


title: Market Areas in General Equilibrium

abstract: This paper proposes a spatial model with a realistic geography where a continuous distribution of agents (e.g., farmers) engages in economic interactions with one location from a finite set (e.g., cities). The spatial structure of the equilibrium consists of a tessellation, i.e., a partition of space into a collection of mutually exclusive market areas. After proving the existence of a unique equilibrium, we characterize how the location of borders and, in the case with mobile labor, the set of inhabited cities change in response to economic shocks. To deal with a two-dimensional space, we draw on tools from computational geometry and from the theory of shape optimization. Finally, we provide an empirical application to illustrate the usefulness of the framework for applied work.


title: Reproducing size distributions of swarms of barchan dunes on Mars and Earth using a mean-field model

abstract: We apply a mean-field model of interactions between migrating barchan dunes, the CAFE model, which includes spontaneous calving, aggregation, fragmentation, and mass-exchange, yielding a steady-state size distribution that can be resolved for different choices of interaction parameters. The CAFE model is applied to empirically measured distributions of dune sizes in two barchan swarms in the north circumpolar region of Mars, three swarms in Morocco, and one in Mauritania, each containing more than 1000 bedforms. When the sizes of bedforms are rescaled by the mean size in each zone two attractor states appear, with the Tarfaya zones all displaying a common distribution and the Martian and Mauritanian zones sharing a different distribution. Comparison of these attractor states with the outputs of the CAFE model reveals that the Tarfaya-type distribution results from a preference for aggregation and fragmentation interactions whereas the Mars-Mauritania distribution is more likely a result of exchange-dominated interactions. We observe that there appears to be a greater number of collisions happening in Tarfaya than in the other areas which is consistent with a greater rate of aggregation-fragmentation processes as suggested by our model. Our comparison with the CAFE model also predicts a universal rule for the outputs of the different types of interactions with exchange interactions favouring the production of two dunes roughly equal in size. Fragmentation interactions often result in three bedforms with volumes approximately in the ratio 1:1:2. Finally, we find that spontaneous calving of dunes does not play an important role in shaping the size distributions in barchan swarms.

categories: physics.geo-ph astro-ph.EP physics.data-an stat.AP


title: Capacitive energy storage in single-file pores: Exactly-solvable models and simulations

abstract: Understanding charge storage in low-dimensional electrodes is crucial for developing novel ecologically friendly devices for capacitive energy storage and conversion, water desalination, etc. Exactly-solvable models allow in-depth analyses and essential physical insights into the charging mechanisms. So far, however, such analytical approaches have been mainly limited to lattice models. Herein, we develop a versatile, exactly-solvable, one-dimensional off-lattice model for charging single-file pores. Unlike the lattice model, this model shows an excellent quantitative agreement with three-dimensional Monte Carlo simulations. With analytical calculations and simulations, we show that the differential capacitance can be bell-shaped (one peak), camel-shaped (two peaks), or have four peaks. Transformations between these capacitance shapes can be induced by changing pore ionophilicity, cation-anion size asymmetry, or by adding solvent. We find that the camel-shaped capacitance, characteristic of dilute electrolytes, appears for strongly ionophilic pores with high ion densities, which we relate to charging mechanisms specific to narrow pores. We also derive a large-voltage asymptotic expression for the capacitance, showing that the capacitance decays to zero as the inverse square of the voltage, $C \sim u^{-2}$. This dependence follows from hard-core interactions and is not captured by the lattice model.


title: Measurement of the jet-particle $v_{2}$ in p--Pb and Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV with ALICE at the LHC

abstract: The elliptic flow of jet particles at midrapidity ($|\eta|$ < 0.8) is measured for the first time in p--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV and extended down to lower $\it{p}_{\rm T}$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. The jet particles are extracted with the two-particle correlation method using a two-dimensional fit, and their $v_{2}$ is calculated with the novel three-particle correlation technique. The $v_{2}$ of jet particles is found positive in high-multiplicity p--Pb collisions and exhibits no dependence on $\it{p}_{\rm T}$ within uncertainties. Comparisons with the jet-particle $v_{2}$ in Pb--Pb collisions and the inclusive charged-particle $v_{2}$ in both collision systems are also discussed, aiming to bring new insights into the understanding of the origin of the high-$\it{p}_{\rm T}$ azimuthal anisotropy observed in small collision systems.


title: TTDFT: A GPU accelerated Tucker tensor DFT code for large-scale Kohn-Sham DFT calculations

abstract: We present the Tucker tensor DFT (TTDFT) code which uses a tensor-structured algorithm with graphic processing unit (GPU) acceleration for conducting ground-state DFT calculations on large-scale systems. The Tucker tensor DFT algorithm uses a localized Tucker tensor basis computed from an additive separable approximation to the Kohn-Sham Hamiltonian. The discrete Kohn-Sham problem is solved using Chebyshev filtering subspace iteration method that relies on matrix-matrix multiplications of a sparse symmetric Hamiltonian matrix and a dense wavefunction matrix, expressed in the localized Tucker tensor basis. These matrix-matrix multiplication operations, which constitute the most computationally intensive step of the solution procedure, are GPU accelerated providing ~8-fold GPU-CPU speedup for these operations on the largest systems studied. The computational performance of the TTDFT code is presented using benchmark studies on aluminum nano-particles and silicon quantum dots with system sizes ranging up to ~7,000 atoms.


title: Dark matter dilution scenarios in the early universe

abstract: When the vacuum like energy of the Higgs potential within the standard model undergoes electroweak phase transition, an influx of entropy into the primordial plasma can lead to a significant dilution of frozen out dark matter density that was already present before the onset of the phase transition. The same effect can take place, if the early Universe was dominated by primordial black holes of small mass, evaporating before the period of Big Bang Nucleosynthesis. In this paper we calculate the dilution factor for the above mentioned scenarios.


title: Complexity, chaos and the moving D3-brane

abstract: We use the wave-function method developed in area of quantum information to investigate the quantum circuit complexity of the small quantum fluctuations around the probe $D_3$ brane moving in $AdS_5\times S^5$ bulk. In our consideration, the reference and target states are chosen as the vacuum state and the squeezed quantum state respectively. The evolution of parameters characterizing the squeezed quantum state are governed by the time-dependent $Schr\ddot{o}dinger$ equation, in which the Hamiltonian operator is derived from the perturbative action of $D_3$ brane. For a quantum chaotic system, some recent works indicate that the evolution of quantum circuit complexity could provide equivalent information like the out-of-time-order correlators. Basing on this inference, our results show that the quantum fluctuations around the non-BPS brane manifestly evolve into the chaotic regime at the late time, while the chaotic behavior is not easy to observe in case of BPS brane. In holographic viewpoint, it implies that the thermodynamic system consist of the $N=4$ supersymmetric particles in non-BPS states evolve into a chaotic system more easily than the one in BPS state.


title: The fate of pairing and spin-charge separation in the presence of long-range antiferromagnetism

abstract: We present a numerical study of competing orders in the 1D $t$-$J$ model with long-range RKKY-like staggered spin interactions. By circumventing the constraints imposed by Mermin-Wagner's theorem, this Hamiltonian can realize long-range N\'eel order at half-filling. We determine the full phase diagram as a function of the exchange and particle density using the density matrix renormalization group (DMRG) method. We show that pairing is disfavored and the AFM insulator and metallic phases are separated by a broad regime with phase segregation, before spin-charge separation re-emerges at low densities. Upon doping, interactions induce a confining potential that binds holons and spinons into full fledged fermionic quasi-particles in a range of parameters and densities. We numerically calculate the photoemission spectrum of the model, showing the appearance of a coherent quasi-particle band splitting away from the holon-spinon continuum with a width determined by $J$ that survives at finite doping. Comparison with analytical results using the self-consistent Born approximation (SCBA) and by solving the spinon-holon problem offer insight into the internal structure of the quasi-particles and help us explain the different features in the spectrum. We discuss how this simple toy-model can teach us about the phenomenology of its higher dimensional counterpart.


title: Optical leakage mitigation in ortho-mode transducer detectors for microwave applications

abstract: Planar ortho-mode transducers (OMTs) are a commonly used method of coupling optical signals between waveguides and on-chip circuitry and detectors. While the ideal OMT-waveguide coupling requires minimal disturbance to the waveguide, when used for mm-wave applications the waveguide is typically constructed from two sections to allow the OMT probes to be inserted into the waveguide. This break in the waveguide is a source of signal leakage and can lead to loss of performance and increased experimental systematic errors. Here we report on the development of new OMT-to-waveguide coupling structures with the goal of reducing leakage at the detector wafer interface. The pixel to pixel optical leakage due to the gap between the coupling waveguide and the backshort is reduced by means of a protrusion that passes through the OMT membrane and electrically connects the two waveguide sections on either side of the wafer. High frequency electromagnetic simulations indicate that these protrusions are an effective method to reduce optical leakage in the gap by ~80% percent, with a ~60% filling factor, relative to an standard OMT coupling architecture. Prototype devices have been designed to characterize the performance of the new design using a relative measurement with varying filling factors. We outline the simulation setup and results, and present a chip layout and sample box that will be used to perform the initial measurements.


title: The [OIII]$\lambda5007$ equivalent width distribution at z $\sim2$: The redshift evolution of the extreme emission line galaxies

abstract: We determine the [OIII]$\lambda5007$ equivalent width (EW) distribution of $1.700<\rm{z}<2.274$ rest-frame UV-selected (M$_{\rm{UV}}<-19$) star-forming galaxies in the GOODS North and South fields. We make use of deep HDUV broadband photometry catalogues for selection and 3D-HST WFC3/IR grism spectra for measurement of line properties. The [OIII]$\lambda5007$ EW distribution allows us to measure the abundance of extreme emission line galaxies (EELGs) within this population. We model a log-normal distribution to the [OIII]$\lambda5007$ rest-frame equivalent widths of galaxies in our sample, with location parameter $\mu=4.24\pm0.07$ and variance parameter $\sigma= 1.33\pm0.06$. This EW distribution has a mean [OIII]$\lambda5007$ EW of 168$\pm1\r{A}$. The fractions of $\rm{z}\sim2$ rest-UV-selected galaxies with [OIII]$\lambda5007$ EWs greater than $500, 750$ and $1000\r{A}$ are measured to be $6.8^{+1.0}_{-0.9}\%$, $3.6^{+0.7}_{-0.6}\%$, and $2.2^{+0.5}_{-0.4}\%$ respectively. The EELG fractions do not vary strongly with UV luminosity in the range ($-21.6<M_{\rm{UV}}<-19.0$) considered in this paper, consistent with findings at higher redshifts. We compare our results to $\rm{z}\sim5$ and $\rm{z}\sim7$ studies where candidate EELGs have been discovered through Spitzer/IRAC colours, and we identify rapid evolution with redshift in the fraction of star-forming galaxies observed in an extreme emission line phase (a rise by a factor $\sim10$ between $\rm{z}\sim2$ and $\rm{z}\sim7$). This evolution is consistent with an increased incidence of strong bursts in the galaxy population of the reionisation era. While this population makes a sub-dominant contribution of the ionising emissivity at $\rm{z}\simeq2$, EELGs are likely to dominate the ionising output in the reionisation era.


title: Uniting Parametric Uncertainty and Tipping Diagrams

abstract: Various subsystems of the Earth system may undergo critical transitions by passing a so-called tipping point, under sustained changes to forcing. For example, the Atlantic Meridional Overturning Circulation (AMOC) is of particular importance for North Atlantic heat transport and is thought to be potentially at risk of tipping. Given a model of such a subsystem that accurately includes the relevant physical processes, whether tipping occurs or not, will depend on model parameters that typically are uncertain. Reducing this parametric uncertainty is important to understand the likelihood of tipping behavior being present in the system and possible tipping locations. In this letter, we develop improved estimates for the parametric uncertainty by inferring probability distributions for the model parameters based on physical constraints and by using a Bayesian inversion technique. To visualize the impact of parametric uncertainty, we extend classical tipping diagrams by visualizing probabilistic bifurcation curves according to the inferred distribution of the model parameter. Furthermore, we highlight the uncertain locations of tipping points along the probabilistic bifurcation curves. We showcase our probabilistic visualizations of the tipping behavior using a simple box-model of the AMOC, the Stommel-Cessi model [5].


title: Tree-Cotree Decomposition of Isogeometric Mortared Spaces in H(curl) on Multi-Patch Domains

abstract: When applying isogeometric analysis to engineering problems, one often deals with multi-patch spline spaces that have incompatible discretisations, e.g. in the case of moving objects. In such cases mortaring has been shown to be advantageous. This contribution discusses the appropriate B-spline spaces needed for the solution of Maxwell's equations in the functions space H(curl) and the corresponding mortar spaces. The main contribution of this paper is to show that in formulations requiring gauging, as in the vector potential formulation of magnetostatic equations, one can remove the discrete kernel subspace from the mortared spaces by the graph-theoretical concept of a tree-cotree decomposition. The tree-cotree decomposition is done based on the control mesh, it works for non-contractible domains, and it can be straightforwardly applied independently of the degree of the B-spline bases. Finally, the simulation workflow is demonstrated using a realistic model of a rotating permanent magnet synchronous machine.


title: A useful lemma for calculating the Hausdorff dimension of certain sets in Engel expansions

abstract: Let $\{s_n\}$ and $\{t_n\}$ be two sequences of positive real numbers. Under some mild conditions on $\{s_n\}$ and $\{t_n\}$, we give the precise formula of the Hausdorff dimension of the set \[ \mathbb{E}(\{s_n\},\{t_n\}):=\Big\{x\in(0,1): s_{n}<d_{n}(x)\leq s_n+t_n, \forall n\geq1\Big\}, \] where $d_n(x)$ denotes the digit of the Engel expansion of $x$. This result improves the Lemma 2.6 of Shang and Wu (2021JNT), and is very useful for calculating the Hausdorff dimension of certain sets in Engel expansions.


title: Assessing present and future risk of water damage using building attributes, meteorology and topography

abstract: Weather-related risk makes the insurance industry inevitably concerned with climate and climate change. Buildings hit by pluvial flooding is a key manifestation of this risk, giving rise to compensations of the induced physical damages and business interruptions. In this work, we establish a nationwide, building-specific risk score for water damage associated with pluvial flooding in Norway. We fit a generalized additive model that relates the number of water damages to a wide range of explanatory variables that can be categorized into building attributes, climatological variables and topographical characteristics. The model assigns a risk score to every location in Norway, based on local topography and climate, which is not only useful for insurance companies, but also for city planning. Combining our model with an ensemble of climate projections allows us to project the (spatially varying) impacts of climate change on the risk of pluvial flooding towards the middle and end of the 21st century.


title: Shock-wave radio probing of solar wind sources in coronal magnetic fields

abstract: The Space Weather effects in the near-Earth environment as well as in atmospheres of other terrestrial planets arise by corpuscular radiation from the Sun, known as the solar wind. The solar magnetic fields govern the solar corona structure. Magnetic-field strength values in the solar wind sources - key information for modeling and forecasting the Space Weather climate - are derived from various solar space- and ground-based observations, but, so far not accounting for specific types of radio bursts. These are "fractured" type II radio bursts attributed to collisions of shock waves with coronal structures emitting the solar wind. Here, we report about radio observations of two "fractured" type II bursts to demonstrate a novel tool for probing of magnetic field variations in the solar wind sources. These results have direct impact on interpretations of this class of bursts and contribute to the current studies of the solar wind emitters.


title: Magnetism and Magnetotransport in the Kagome Antiferromagnet $\text{Mn}_3\text{Ge}$

abstract: We perform classical Monte Carlo and stochastic Landau-Lifshitz-Gilbert simulations to study temperature dependent magnetism of Kagome antiferromagnet (AFM) Weyl metal $\text{Mn}_3\text{Ge}$ and find that a long range chiral order sets in at a transition temperature well below the N{\'e}el temperature ($T_N$). Based on the crystalline symmetries, imposed by the chiral magnetic order, we argue for the presence of multiple iso-energetic Weyl nodes (nodes that are at same energy and with congruent Fermi surface around them) near chemical potential. Using the semi-classical Boltzmann equations, we show that the combined contribution to the net longitudinal magnetoconductance (LMC) and the planar Hall conductance (PHC) from tilted Weyl nodes can lead to signatures, qualitatively distinct from that of a single pair of Weyl nodes. In particular, we show that magnetic orders with different chiralities can give rise to different periods in LMC and PHC as a function of the in-plane magnetic field direction. This is ultimately related to differences in the symmetry-imposed constraints on the Weyl nodes.


title: Some recent advances in urban system science: models and data

abstract: Cities are characterized by the presence of a dense population with a high potential for interactions between individuals of diverse backgrounds. They appear in parallel to the Neolithic revolution a few millennia ago. The advantages brought in terms of agglomeration for economy, innovation, social and cultural advancements have kept them as a major landmark in recent human history. There are many different aspects to study in urban systems from a scientific point of view, just to name a few one can concentrate in demography and population evolution, mobility, economic output, land use and urban planning, home accessibility and real estate market, energy and water consumption, waste processing, health, education, integration of minorities, etc. In the last decade, the introduction of communication and information technologies have enormously facilitated the collection of datasets on these and other questions, making possible a more quantitative approach to city science. All these topics have been addressed in many works in the literature, and we do not intend to offer here a systematic review. Instead, we will only provide a brief taste of some of these above-mentioned aspects, which could serve as an introduction to a subsequent special number. Such a non-systematic view will lead us to leave outside many relevant papers, and for this we apologise.


title: Towards Comparative Physical Interpretation of Spatial Variability Aware Neural Networks: A Summary of Results

abstract: Given Spatial Variability Aware Neural Networks (SVANNs), the goal is to investigate mathematical (or computational) models for comparative physical interpretation towards their transparency (e.g., simulatibility, decomposability and algorithmic transparency). This problem is important due to important use-cases such as reusability, debugging, and explainability to a jury in a court of law. Challenges include a large number of model parameters, vacuous bounds on generalization performance of neural networks, risk of overfitting, sensitivity to noise, etc., which all detract from the ability to interpret the models. Related work on either model-specific or model-agnostic post-hoc interpretation is limited due to a lack of consideration of physical constraints (e.g., mass balance) and properties (e.g., second law of geography). This work investigates physical interpretation of SVANNs using novel comparative approaches based on geographically heterogeneous features. The proposed approach on feature-based physical interpretation is evaluated using a case-study on wetland mapping. The proposed physical interpretation improves the transparency of SVANN models and the analytical results highlight the trade-off between model transparency and model performance (e.g., F1-score). We also describe an interpretation based on geographically heterogeneous processes modeled as partial differential equations (PDEs).


title: A 400 Trillion-Grid Vlasov Simulation on Fugaku Supercomputer: Large-Scale Distribution of Cosmic Relic Neutrinos in a Six-dimensional Phase Space

abstract: We report a Vlasov simulation of cosmic relic neutrinos combined with N-body simulation of cold dark matter in the context of large-scale structure formation in the Universe performed on Fugaku supercomputer. Gravitational dynamics of the neutrinos is followed, for the first time, by directly integrating the Vlasov equation in a six-dimensional phase space. Our largest simulation combines the Vlasov simulation on 400 trillion grids and 330 billion-body calculations in a self-consistent manner, and reproduces accurately the nonlinear dynamics of neutrinos in the Universe. The novel high-order Vlasov solver is optimized by combining an array of state-of-the-art numerical schemes and fully utilizing the SIMD instructions on the A64FX processors. Time-To-Solution of our simulation is an order of magnitude shorter than the largest N-body simulations. The performance scales excellently with up to 147,456 nodes (7 million CPU cores) on Fugaku; the weak and strong scaling efficiencies are 82% - 96% and 82% - 93%, respectively.


title: First numerical evidence of Janssen-Oerding's prediction in a three-dimensional spin model far from equilibrium

abstract: Jansen and Oerding [H. K. Janssen, K. Oerding, J. Phys. A: Math. Gen. 27, 715 (1994)] predicted an interesting anomalous tricritical dynamic behavior in three-dimensional models via renormalization group theory. However, we verify a lack of literature about the computational verification of this universal behavior. Here, we used some tricks to capture the log corrections and the parameters predicted by these authors using the three-dimensional Blume-Capel model. In addition, we also performed a more detailed study of the dynamic localization of the phase diagram via power laws optimization. We quantify the crossover phenomena by computing the critical exponents near the tricritical point.


title: Trustworthy Pre-Processing of Sensor Data in Data On-chaining Workflows for Blockchain-based IoT Applications

abstract: Prior to provisioning sensor data to smart contracts, a pre-processing of the data on intermediate off-chain nodes is often necessary. When doing so, originally constructed cryptographic signatures cannot be verified on-chain anymore. This exposes an opportunity for undetected manipulation and presents a problem for applications in the Internet of Things where trustworthy sensor data is required on-chain. In this paper, we propose trustworthy pre-processing as enabler for end-to-end sensor data integrity in data on-chaining workflows. We define requirements for trustworthy pre-processing, present a model and common workflow for data on-chaining, select off-chain computation utilizing Zero-knowledge Proofs (ZKPs) and Trusted Execution Environments (TEEs) as promising solution approaches, and discuss both our proof-of-concept implementations and initial experimental, comparative evaluation results. The importance of trustworthy pre-processing and principle solution approaches are presented, addressing the major problem of end-to-end sensor data integrity in blockchain-based IoT applications.


title: Hybrid Quantum-Classical Algorithms for Loan Collection Optimization with Loan Loss Provisions

abstract: Banks are required to set aside funds in their income statement, known as a loan loss provision (LLP), to account for potential loan defaults and expenses. By treating the LLP as a global constraint, we propose a hybrid quantum-classical algorithm to solve Quadratic Constrained Binary Optimization (QCBO) models for loan collection optimization. The objective is to find a set of optimal loan collection actions that maximizes the expected net profit presented to the bank as well as the financial welfare in the financial network of loanees, while keeping the LLP at its minimum. Our algorithm consists of three parts: a classical divide-and-conquer algorithm to enable a large-scale optimization, a quantum alternating operator ansatz (QAOA) algorithm to maximize the objective function, and a classical sampling algorithm to handle the LLP. We apply the algorithm to a real-world data set with 600 loanees and 5 possible collection actions. The QAOA is performed using up to 35 qubits on a classical computer. We show that the presence of the QAOA can improve the expected net profit by approximately $70\%$, compared to when the QAOA is absent from the hybrid algorithm. Our work illustrates the use of near-term quantum devices to tackle real-world optimization problems.


title: From Theories on Styles to their Transfer in Text: Bridging the Gap with a Hierarchical Survey

abstract: Humans are naturally endowed with the ability to write in a particular style. They can, for instance, re-phrase a formal letter in an informal way, convey a literal message with the use of figures of speech or edit a novel by mimicking the style of some well-known authors. Automating this form of creativity constitutes the goal of style transfer. As a natural language generation task, style transfer aims at rewriting existing texts, and specifically, it creates paraphrases that exhibit some desired stylistic attributes. From a practical perspective, it envisions beneficial applications, like chatbots that modulate their communicative style to appear empathetic, or systems that automatically simplify technical articles for a non-expert audience. Several style-aware paraphrasing methods have attempted to tackle style transfer. A handful of surveys give a methodological overview of the field, but they do not support researchers to focus on specific styles. With this paper, we aim at providing a comprehensive discussion of the styles that have received attention in the transfer task. We organize them in a hierarchy, highlighting the challenges for the definition of each of them, and pointing out gaps in the current research landscape. The hierarchy comprises two main groups. One encompasses styles that people modulate arbitrarily, along the lines of registers and genres. The other group corresponds to unintentionally expressed styles, due to an author's personal characteristics. Hence, our review shows how these groups relate to one another, and where specific styles, including some that have not yet been explored, belong in the hierarchy. Moreover, we summarize the methods employed for different stylistic families, hinting researchers towards those that would be the most fitting for future research.


title: 2D-2FA: A New Dimension in Two-Factor Authentication

abstract: We propose a two-factor authentication (2FA) mechanism called 2D-2FA to address security and usability issues in existing methods. 2D-2FA has three distinguishing features: First, after a user enters a username and password on a login terminal, a unique $\textit{identifier}$ is displayed to her. She $\textit{inputs}$ the same identifier on her registered 2FA device, which ensures appropriate engagement in the authentication process. Second, a one-time PIN is computed on the device and $\textit{automatically}$ transferred to the server. Thus, the PIN can have very high entropy, making guessing attacks infeasible. Third, the identifier is also incorporated into the PIN computation, which renders $\textit{concurrent attacks}$ ineffective. Third-party services such as push-notification providers and 2FA service providers, do not need to be trusted for the security of the system. The choice of identifiers depends on the device form factor and the context. Users could choose to draw patterns, capture QR codes, etc. We provide a proof of concept implementation, and evaluate performance, accuracy, and usability of the system. We show that the system offers a lower error rate (about half) and better efficiency (2-3 times faster) compared to the commonly used PIN-2FA. Our study indicates a high level of usability with a SUS of 75, and a high perception of efficiency, security, accuracy, and adoptability.


title: A comparison of Cahn-Hilliard and Navier-Stokes-Cahn-Hilliard models on manifolds

abstract: We consider phase-field models with and without lateral flow for the numerical simulation of lateral phase separation and coarsening in lipid membranes. For the numerical solution of these models, we apply an unfitted finite element method that is flexible in handling complex and possibly evolving shapes in the absence of an explicit surface parametrization. Through several numerical tests, we investigate the effect of the presence of lateral flow on the evolution of phases. In particular, we focus on understanding how variable line tension, viscosity, membrane composition, and surface shape affect the pattern formation. Keywords: Lateral phase separation, surface Cahn-Hilliard equation, lateral flow, surface Navier-Stokes-Cahn-Hilliard system, TraceFEM


title: Higher differentiability for a class of problems under p,q subquadratic growth

abstract: We study the higher differentiability for nonlinear elliptic equation in divergence form $\mathcal{A}(x,Du)=b(x)$. The result covers the cases in which $\mathcal{A}(x, \xi)$ satisfies $p,q$ growth, with $1<p<2$ in $\xi$ and a Sobolev dependence of with respect to $x$. By means of an a-priori estimate we ensure the $W^{2,p}_{\mathrm{loc}}(\Omega)$-property for the solution of the boundary value problem.


title: Elasto-acoustic modelling and simulation for the seismic response of structures: The case of the Tahtal{\i} dam in the 2020 \.Izmir earthquake

abstract: As a mean to assess the risk dam structures are exposed to during earthquakes, we employ an abstract mathematical, three dimensional, elasto-acoustic coupled wave-propagation model taking into account (i) the dam structure itself, embedded into (ii) its surrounding topography, (iii) different material soil layers, (iv) the seismic source as well as (v) the reservoir lake filled with water treated as an acoustic medium. As a case study for extensive numerical simulations we consider the magnitude 7 seismic event of the 30$^{\rm th}$ of October 2020 taking place in the Icarian Sea (Greece) and the Tahtali dam around 30 km from there (Turkey). A challenging task is to resolve the multiple length scales that are present due to the huge differences in size between the dam building structure and the area of interest, considered for the propagation of the earthquake. Interfaces between structures and highly non-conforming meshes on different scales are resolved by means of a discontinuous Galerkin approach. The seismic source is modeled using inversion data about the real fault plane. Ultimately, we perform a real data driven, multi-scale, full source-to-site, physics based simulation based on the discontinuous Galerkin spectral element method, which allows to precisely validate the ground motion experienced along the Tahtali dam, comparing the synthetic seismograms against actually observed ones. A comparison with a more classical computational method, using a plane wave with data from a deconvolved seismogram reading as an input, is discussed.


title: An integrable semi-discretization of the modified Camassa-Holm equation with linear dispersion term

abstract: In the present paper, we are concerned with integrable discretization of a modified Camassa-Holm equation with linear dispersion term. The key of the construction is the semi-discrete analogue for a set of bilinear equations of the modified Camassa-Holm equation. Firstly, we show that these bilinear equations and their determinant solutions either in Gram-type or Casorati-type can be reduced from the discrete KP equation through Miwa transformation. Then, by scrutinizing the reduction process, we obtain a set of semi-discrete bilinear equations and their general soliton solution in Gram-type or Casorati-type determinant form. Finally, by defining dependent variables and discrete hodograph transformations, we are able to derive an integrable semi-discrete analogue of the modified Camassa-Holm equation. It is also shown that the semi-discrete modified Camassa-Holm equation converges to the continuous one in the continuum limit.


title: Quality control, data cleaning, imputation

abstract: This chapter addresses important steps during the quality assurance and control of RWD, with particular emphasis on the identification and handling of missing values. A gentle introduction is provided on common statistical and machine learning methods for imputation. We discuss the main strengths and weaknesses of each method, and compare their performance in a literature review. We motivate why the imputation of RWD may require additional efforts to avoid bias, and highlight recent advances that account for informative missingness and repeated observations. Finally, we introduce alternative methods to address incomplete data without the need for imputation.


title: Polariton interaction in one-dimensional arrays of atoms coupled to waveguides

abstract: Photons strongly coupled to material systems constitute a novel system for studying the dynamics of non-equilibrium quantum many-body systems. We give a fully analytical description of the dynamics of photons coupled to a one-dimensional array of two-level atoms. The description incorporates input, scattering inside the medium, and the ejection of photons from the array. We show that inelastic collisions, previously identified in small systems, also occur in infinite systems and elucidate the physics behind this effect. The developed theory constitutes an effective field theory for the dynamics, which can be used as a starting point for studies of many-body dynamics. We discuss different parameter regimes and touch upon possible applications and the limit of many excitations.


title: Upper and Lower Bounds for End-to-End Risks in Stochastic Robot Navigation

abstract: We present novel upper and lower bounds to estimate the collision probability of motion plans for autonomous agents with discrete-time linear Gaussian dynamics. Motion plans generated by planning algorithms cannot be perfectly executed by autonomous agents in reality due to the inherent uncertainties in the real world. Estimating collision probability is crucial to characterize the safety of trajectories and plan risk optimal trajectories. Our approach is an application of standard results in probability theory including the inequalities of Hunter, Kounias, Frechet, and Dawson. Using a ground robot navigation example, we numerically demonstrate that our method is considerably faster than the naive Monte Carlo sampling method and the proposed bounds are significantly less conservative than Boole's bound commonly used in the literature.


title: Internal Stark effect of single-molecule fluorescence

abstract: The optical properties of chromophores can be efficiently tuned by electrostatic fields generated in their close environment, a phenomenon that plays a central role for the optimization of complex functions within living organisms where it is known as internal Stark effect (ISE). Here, we realised an ISE experiment at the lowest possible scale, by monitoring the Stark shift generated by charges confined within a single chromophore on its emission energy. To this end, a scanning tunneling microscope (STM) functioning at cryogenic temperatures is used to sequentially remove the two central protons of a free-base phthalocyanine chromophore deposited on a NaCl-covered Ag(111) surface. STM-induced fluorescence measurements reveal spectral shifts that are associated to the electrostatic field generated by the internal charges remaining in the chromophores upon deprotonation.


title: Methane formation in cold regions from carbon atoms and molecular hydrogen

abstract: Methane is typically thought to be formed in the solid state on the surface of cold interstellar icy grain mantles via the successive atomic hydrogenation of a carbon atom. In the current work we investigate the potential role of molecular hydrogen in the CH$_4$ reaction network. We make use of an ultra-high vacuum cryogenic setup combining an atomic carbon atom beam and both atomic and/or molecular beams of hydrogen and deuterium on a H$_2$O ice. These experiments lead to the formation of methane isotopologues detected in situ through reflection absorption infrared spectroscopy. Most notably, CH$_4$ is formed in an experiment combining C atoms with H$_2$ on amorphous solid water, albeit slower than in experiments with H atoms present. Furthermore, CH$_2$D$_2$ is detected in an experiment of C atoms with H$_2$ and D$_2$ on H$_2$O ice. CD$_4$, however, is only formed when D atoms are present in the experiment. These findings have been rationalized by means of computational chemical insights. This leads to the following conclusions: a) the reaction C + H$_2$ -> CH$_2$ can take place, although not barrierless in the presence of water, b) the reaction CH + H$_2$ -> CH$_3$ is barrierless, but has not yet been included in astrochemical models, c) the reactions CH$_2$ + H$_2$ -> CH$_3$ + H and CH$_3$ + H$_2$ -> CH$_4$ + H can take place only via a tunneling mechanism and d) molecular hydrogen possibly plays a more important role in the solid-state formation of methane than assumed so far.


title: Computing the Invariant Circle and its Stable Manifolds for a 2-D Map by the Parameterization Method: Effective Algorithms and Rigorous Proofs of Convergence

abstract: We present and analyze rigorously a quadratically convergent algorithm to compute an invariant circle for 2-dimensional maps along with the corresponding foliation by stable manifolds. We prove that when the algorithm starts from an initial guess that satisfies the invariance equation very approximately (depending on some condition numbers, evaluated on the approximate solution), then the algorithm converges to a true solution which is close to the initial guess. The convergence is faster than exponential in smooth norms. The distance from the exact solution and the approximation is bounded by the initial error. This allows validating the numerical approximations (a-posteriori results). It also implies the usual persistence formulations since the exact solutions of the invariance equation for a model are approximate solutions for a similar model. The algorithm we present works irrespective of whether the dynamics on the invariant circle is a rotation or it is phase-locked. The condition numbers required do not involve any global qualitative properties of the map. They are obtained by evaluating derivatives of the initial guess, derivatives of the map in a neighborhood of the guess, performing algebraic operations, and taking suprema. The proof of the convergence is based on a general Nash-Moser implicit function theorem specially tailored for this problem. The Nash-Moser procedure has unusual properties. As it turns out, the regularity requirements are not very severe (only 2 derivatives suffice). We hope that this implicit function theorem may be of independent interest and have presented it in a self-contained appendix. The algorithm in this paper is very practical since it converges quadratically, and it requires moderate storage and operation count. Details of the implementation and results of the runs are described in a companion paper [YdlL21].


title: Universal fidelity reduction of quantum operations from weak dissipation

abstract: Quantum information processing is in real systems often limited by dissipation, stemming from remaining uncontrolled interaction with microscopic degrees of freedom. Given recent experimental progress, we consider weak dissipation, resulting in a small error probability per operation. Here, we find a simple formula for the fidelity reduction of any desired quantum operation. Interestingly, this reduction is independent of the specific operation; it depends only on the operation time and the dissipation. Using our formula, we investigate the situation where dissipation in different parts of the system have correlations, which is detrimental for the successful application of quantum error correction. Surprisingly, we find that a large class of correlations gives the same fidelity reduction as uncorrelated dissipation of similar strength.


title: Distributing Deep Learning Hyperparameter Tuning for 3D Medical Image Segmentation

abstract: Most research on novel techniques for 3D Medical Image Segmentation (MIS) is currently done using Deep Learning with GPU accelerators. The principal challenge of such technique is that a single input can easily cope computing resources, and require prohibitive amounts of time to be processed. Distribution of deep learning and scalability over computing devices is an actual need for progressing on such research field. Conventional distribution of neural networks consist in data parallelism, where data is scattered over resources (e.g., GPUs) to parallelize the training of the model. However, experiment parallelism is also an option, where different training processes are parallelized across resources. While the first option is much more common on 3D image segmentation, the second provides a pipeline design with less dependence among parallelized processes, allowing overhead reduction and more potential scalability. In this work we present a design for distributed deep learning training pipelines, focusing on multi-node and multi-GPU environments, where the two different distribution approaches are deployed and benchmarked. We take as proof of concept the 3D U-Net architecture, using the MSD Brain Tumor Segmentation dataset, a state-of-art problem in medical image segmentation with high computing and space requirements. Using the BSC MareNostrum supercomputer as benchmarking environment, we use TensorFlow and Ray as neural network training and experiment distribution platforms. We evaluate the experiment speed-up, showing the potential for scaling out on GPUs and nodes. Also comparing the different parallelism techniques, showing how experiment distribution leverages better such resources through scaling. Finally, we provide the implementation of the design open to the community, and the non-trivial steps and methodology for adapting and deploying a MIS case as the here presented.


title: Multiscale Finite Element Methods for an Elliptic Optimal Control Problem with Rough Coefficients

abstract: We investigate multiscale finite element methods for an elliptic distributed optimal control problem with rough coefficients. They are based on the (local) orthogonal decomposition methodology of M\aa lqvist and Peterseim.


title: A probabilistic view of latent space graphs and phase transitions

abstract: We study random graphs with latent geometric structure, where the probability of each edge depends on the underlying random positions corresponding to the two endpoints. We focus on the setting where this conditional probability is a general monotone increasing function of the inner product of two vectors; such a function can naturally be viewed as the cumulative distribution function of some independent random variable. We consider a one-parameter family of random graphs, characterized by the variance of this random variable, that smoothly interpolates between a random dot product graph and an Erd\H{o}s--R\'enyi random graph. We prove phase transitions of detecting geometry in these graphs, in terms of the dimension of the underlying geometric space and the variance parameter of the conditional probability. When the dimension is high or the variance is large, the graph is similar to an Erd\H{o}s--R\'enyi graph with the same edge density that does not possess geometry; in other parameter regimes, there is a computationally efficient signed triangle statistic that distinguishes them. The proofs make use of information-theoretic inequalities and concentration of measure phenomena.

categories: math.PR cs.IT cs.SI math.IT math.ST stat.TH


title: Carbon Atom Reactivity with Amorphous Solid Water: H$_2$O Catalyzed Formation of H$_2$CO

abstract: We report new computational and experimental evidence of an efficient and astrochemically relevant formation route to formaldehyde (H$_2$CO). This simplest carbonylic compound is central to the formation of complex organics in cold interstellar clouds, and is generally regarded to be formed by the hydrogenation of solid-state carbon monoxide. We demonstrate H$_2$CO formation via the reaction of carbon atoms with amorphous solid water. Crucial to our proposed mechanism is a concerted proton transfer catalyzed by the water hydrogen bonding network. Consequently, the reactions $^3$C + H$_2$O -> $^3$HCOH and $^1$HCOH -> $^1$H$_2$CO can take place with low or without barriers, contrary to the high-barrier traditional internal hydrogen migration. These low barriers or absence thereof explain the very small kinetic isotope effect in our experiments when comparing the formation of H$_2$CO to D$_2$CO. Our results reconcile the disagreement found in the literature on the reaction route: C + H$_2$O -> H$_2$CO.


title: Non-equilibrium quantum thermodynamics of a particle trapped in a controllable time-varying potential

abstract: Many advanced quantum techniques feature non-Gaussian dynamics, and the ability to manipulate the system in that domain is the next-stage in many experiments. One example of meaningful non-Gaussian dynamics is that of a double-well potential. Here we study the dynamics of a levitated nanoparticle undergoing the transition from an harmonic potential to a double-well in a realistic setting, subjecting to both thermalisation and localisation. We characterise the dynamics of the nanoparticle from a thermodynamic point-of-view, investigating the dynamics with the Wehrl entropy production and its rates. Furthermore, we investigate coupling regimes where the the quantum effect and thermal effect are of the same magnitude, and look at suitable squeezing of the initial state that provides the maximum coherence. The effects and the competitions of the unitary and the dissipative parts onto the system are demonstrated. We quantify the requirements to relate our results to a bonafide experiment with the presence of the environment, and discuss the experimental interpretations of our results in the end.


title: Balanced Excited Random Walk in Two Dimensions

abstract: We give non-trivial upper and lower bounds on the range of the so-called Balanced Excited Random Walk in two dimensions, and verify a conjecture of Benjamini, Kozma and Schapira. To the best of our knowledge these are the first non-trivial results for this 2-dimensional model


title: Laboratory observation and astronomical search of 1-cyano propargyl radical, HCCCHCN

abstract: The reaction between carbon atoms and vinyl cyanide, CH2CHCN, is a formation route to interstellar 3-cyano propargyl radical, CH2C3N, a species that has recently been discovered in space. The 1-cyano propargyl radical (HC3HCN), an isomer of CH2C3N, is predicted to be produced in the same reaction at least twice more effciently than CH2C3N. Hence, HC3HCN is a plausible candidate to be observed in space as well. We aim to generate the HC3HCN radical in the gas phase in order to investigate its rotational spectrum. The derived spectroscopic parameters for this species will be used to obtain reliable frequency predictions to support its detection in space.The HC3HCN radical was produced by an electric discharge, and its rotational spectrum was characterized using a Balle-Flygare narrowband-type Fourier-transform microwave spectrometer operating in the frequency region of 4-40 GHz. The spectral analysis was supported by high-level ab initio calculations. A total of 193 hyperfine components that originated from 12 rotational transitions, a- and b-type, were measured for the HC3HCN radical. The analysis allowed us to accurately determine 22 molecular constants, including rotational and centrifugal distortion constants as well as the fine and hyperfine constants. Transition frequency predictions were used to search for the HC3HCN radical in TMC-1 using the QUIJOTE survey between 30.1-50.4 GHz. We do not detect HC3HCN in TMC-1 and derive a 3sigma upper limit to its column density of 6.0e11 cm-2.


title: Friendly Cut Sparsifiers and Faster Gomory-Hu Trees

abstract: We devise new cut sparsifiers that are related to the classical sparsification of Nagamochi and Ibaraki [Algorithmica, 1992], which is an algorithm that, given an unweighted graph $G$ on $n$ nodes and a parameter $k$, computes a subgraph with $O(nk)$ edges that preserves all cuts of value up to $k$. We put forward the notion of a friendly cut sparsifier, which is a minor of $G$ that preserves all friendly cuts of value up to $k$, where a cut in $G$ is called friendly if every node has more edges connecting it to its own side of the cut than to the other side. We present an algorithm that, given a simple graph $G$, computes in almost-linear time a friendly cut sparsifier with $\tilde{O}(n \sqrt{k})$ edges. Using similar techniques, we also show how, given in addition a terminal set $T$, one can compute in almost-linear time a terminal sparsifier, which preserves the minimum $st$-cut between every pair of terminals, with $\tilde{O}(n \sqrt{k} + |T| k)$ edges. Plugging these sparsifiers into the recent $n^{2+o(1)}$-time algorithms for constructing a Gomory-Hu tree of simple graphs, along with a relatively simple procedure for handling the unfriendly minimum cuts, we improve the running time for moderately dense graphs (e.g., with $m=n^{1.75}$ edges). In particular, assuming a linear-time Max-Flow algorithm, the new state-of-the-art for Gomory-Hu tree is the minimum between our $(m+n^{1.75})^{1+o(1)}$ and the known $m n^{1/2+o(1)}$. We further investigate the limits of this approach and the possibility of better sparsification. Under the hypothesis that an $\tilde{O}(n)$-edge sparsifier that preserves all friendly minimum $st$-cuts can be computed efficiently, our upper bound improves to $\tilde{O}(m+n^{1.5})$ which is the best possible without breaking the cubic barrier for constructing Gomory-Hu trees in non-simple graphs.


title: Studies of high-field QED with the LUXE experiment at the European XFEL

abstract: The LUXE experiment aims at studying high-field QED in electron-laser and photon-laser interactions, with the 16.5 GeV electron beam of the European XFEL and a laser beam with power of up to 350 TW. The experiment will measure the spectra of electrons and photons in non-linear Compton scattering where production rates in excess of $10^9$ are expected per 1 Hz bunch crossing. At the same time positrons from pair creation in either the two-step trident process or the Breit--Wheeler process will be measured, where the expected rates range from $10^{-3}$ to $10^4$ per bunch crossing, depending on the laser power and focus. These measurements have to be performed in the presence of low-energy high radiation-background. To meet these challenges, for high-rate electron and photon fluxes, the experiment will use Cherenkov radiation detectors, scintillator screens, sapphire sensors as well as lead-glass monitors for back-scattering off the beam-dump. A four-layer silicon-pixel tracker and a compact sampling electromagnetic calorimeter will be used to measure the positron spectra. The layout of the experiment and the expected performance under the harsh radiation conditions will be presented.


title: Computing the Invariant Circle and the Foliation by Stable Manifolds for a 2-D Map by the Parameterization Method: Numerical Implementation and Results

abstract: We present and implement an algorithm for computing the invariant circle and the corresponding stable manifolds for 2-dimensional maps. The algorithm is based on the parameterization method, and it is backed up by an a-posteriori theorem established in [YdlL21]. The algorithm works irrespective of whether the internal dynamics in the invariant circle is a rotation or it is phase-locked. The algorithm converges quadratically and the number of operations and memory requirements for each step of the iteration is linear with respect to the size of the discretization. We also report on the result of running the implementation in some standard models to uncover new phenomena. In particular, we explored a bundle merging scenario in which the invariant circle loses hyperbolicity because the angle between the stable directions and the tangent becomes zero even if the rates of contraction are separated. We also discuss and implement a generalization of the algorithm to 3 dimensions, and implement it on the 3-dimensional Fattened Arnold Family (3D-FAF) map with non-resonant eigenvalues and present numerical results.


title: Sequential motion planning assisted by group actions

abstract: We study higher analogues of effective and effectual topological complexity of spaces equipped with a group action. These are $G$-homotopy invariant and are motivated by the (higher) motion planning problem of $G$-spaces for which their group action is thought of as an external system assisting the motion planning. Related to this interpretation we define what we call orbital topological complexity, which is also a $G$-homotopy invariant that provides an upper bound for the topological complexity of the quotient space by the group action. We apply these concepts to actions of the group of order two on orientable surfaces and spheres.


title: Application of 2-D Convolutional Neural Networks for Damage Detection in Steel Frame Structures

abstract: In this paper, we present an application of 2-D convolutional neural networks (2-D CNNs) designed to perform both feature extraction and classification stages as a single organism to solve the highlighted problems. The method uses a network of lighted CNNs instead of deep and takes raw acceleration signals as input. Using lighted CNNs, in which every one of them is optimized for a specific element, increases the accuracy and makes the network faster to perform. Also, a new framework is proposed for decreasing the data required in the training phase. We verified our method on Qatar University Grandstand Simulator (QUGS) benchmark data provided by Structural Dynamics Team. The results showed improved accuracy over other methods, and running time was adequate for real-time applications.


title: Stability and magnetic properties of grain boundaries in the inverse Heusler phase Fe$_2$CoGa and in bcc Fe

abstract: We investigate grain boundaries (GBs) in the cubic inverse Heusler phase Fe$_2$CoGa by means of first principles calculations based on density functional theory. Besides the energetic stability, the analysis focuses on the magnetic properties of a set of 16 GB structures in this intermetallic phase. We determine the integrated excess magnetization across the GB and analyze it in terms of the projected local magnetic moments of the atoms and their local Voronoi volumes. The results are systematically compared to those of corresponding GBs in body-centered cubic (bcc) Fe. The studied GBs in Fe$_2$CoGa may have a considerably increased magnetization at the GB, up to more than twice as much as in bcc Fe, depending on the GB type, while geometrical quantities like GB widening or local GB excess volume distributions are similar for both phases. We explain this difference by the higher flexibility of the ternary Fe$_2$CoGa phase in compensating the disturbance of a crystal defect by structural relaxation. The GB structures therefore have a lower energy accompanied with increased local magnetic moments of the Co and half of the Fe atoms within a distance of a few angstrom around the GB plane.


title: Charge Density Wave Activated Excitons in TiSe$_2$-MoSe$_2$ Heterostructures

abstract: Layered materials enable the assembly of a new class of heterostructures where lattice-matching is no longer a requirement. Interfaces in these heterostructures therefore become a fertile ground for unexplored physics as dissimilar phenomena can be coupled via proximity effects. In this article we identify an unexpected photoluminescence (PL) peak when MoSe$_2$ interacts with TiSe$_2$. A series of temperature-dependent and spatially-resolved PL measurements reveal this peak is unique to the TiSe$_2$-MoSe$_2$ interface, higher in energy compared to the neutral exciton, and exhibits exciton-like characteristics. The feature disappears at the TiSe$_2$ charge density wave transition, suggesting that the density wave plays an important role in the formation of this new exciton. We present several plausible scenarios regarding the origin of this peak that individually capture some aspects of our observations, but cannot fully explain this feature. These results therefore represent a fresh challenge for the theoretical community and provide a fascinating way to engineer excitons through interactions with charge density waves.


title: Contextuality, Fine-Tuning, and Teleological Explanation

abstract: I assess various proposals for the source of the intuition that there is something problematic about contextuality, ultimately concluding that contextuality is best thought of in terms of fine-tuning. I then argue that as with other fine-tuning problems in quantum mechanics, this behaviour can be understood as a manifestation of teleological features of physics. Finally I discuss several formal mathematical frameworks that have been used to analyse contextuality and consider how their results should be interpreted by scientific realists. In the course of this discussion I obtain several new mathematical results - I demonstrate that preparation contextuality is a form of fine-tuning, I show that measurement contextuality can be explained by appeal to a global constraint forbidding closed causal loops, and I demonstrate how negative probabilities can arise from a classical ontological model together with an epistemic restriction.


title: Dynamic atmospheres and winds of cool luminous giants, II. Gradual Fe enrichment of wind-driving silicate grains

abstract: The winds observed around AGB stars are generally attributed to radiation pressure on dust formed in the dynamical atmospheres of these long-period variables. The composition of wind-driving grains is affected by a feedback between their optical properties and the resulting heating due to stellar radiation. We explore the gradual Fe enrichment of wind-driving silicate grains in M-type AGB stars to derive typical values for Fe/Mg and test the effects on wind properties and synthetic spectra. We present new DARWIN models that allow for the growth of silicate grains with a variable Fe/Mg ratio and predict mass-loss rates, wind velocities, and grain properties. Synthetic spectra and other observables are computed with the COMA code. The self-regulating feedback between grain composition and radiative heating, in combination with quickly falling densities in the stellar wind, leads to low values of Fe/Mg, typically a few percent. Nevertheless, the models show distinct silicate features around 10 and 18 microns. Fe enrichment affects visual and near-IR photometry moderately, and the new DARWIN models agree well with observations in (J-K) vs. (V-K) and Spitzer color-color diagrams. The enrichment of the silicate dust with Fe is a secondary process, taking place in the stellar wind on the surface of large Fe-free grains that have initiated the outflow. Therefore, the mass-loss rates are basically unaffected, while the wind velocities tend to be slightly higher than in corresponding models with Fe-free silicate dust. The gradual Fe enrichment of silicate grains in the inner wind region should produce signatures observable in mid-IR spectro-interferometrical measurements. Mass-loss rates derived from existing DARWIN models, based on Fe-free silicates, can be applied to stellar evolution models since the mass-loss rates are not significantly affected by the inclusion of Fe in the silicate grains.


title: Hyperparameter Tuning is All You Need for LISTA

abstract: Learned Iterative Shrinkage-Thresholding Algorithm (LISTA) introduces the concept of unrolling an iterative algorithm and training it like a neural network. It has had great success on sparse recovery. In this paper, we show that adding momentum to intermediate variables in the LISTA network achieves a better convergence rate and, in particular, the network with instance-optimal parameters is superlinearly convergent. Moreover, our new theoretical results lead to a practical approach of automatically and adaptively calculating the parameters of a LISTA network layer based on its previous layers. Perhaps most surprisingly, such an adaptive-parameter procedure reduces the training of LISTA to tuning only three hyperparameters from data: a new record set in the context of the recent advances on trimming down LISTA complexity. We call this new ultra-light weight network HyperLISTA. Compared to state-of-the-art LISTA models, HyperLISTA achieves almost the same performance on seen data distributions and performs better when tested on unseen distributions (specifically, those with different sparsity levels and nonzero magnitudes). Code is available: https://github.com/VITA-Group/HyperLISTA.


title: A mobile ion in a Fermi sea

abstract: The remarkable single particle control of individual ions combined with the versatility of ultracold atomic gases makes hybrid ion-atom system an exciting new platform for quantum simulation of few- and many-body quantum physics. Here, we study theoretically the properties of a mobile ion immersed in a quantum degenerate gas of fermionic atoms. Using an effective low-energy atom-ion interaction together with a well established approach that includes exactly two-body correlations, we calculate the full spectral response of the ion and demonstrate the existence of several quasiparticle branches, which are charged analogues of the Fermi polaron observed in neutral atomic gases. Due to the long-range nature of the atom-ion interaction, these ionic Fermi polarons have several properties distinct from their neutral counterparts such as the simultaneous presence of several stable states and smooth transitions from repulsive to attractive polarons with increasing interaction strength. Surprisingly, the residue of the ionic polaron is shown to increase with the Fermi density for fixed interaction strength, which is in marked contrast to the neutral polaron. The properties of the ionic polaron approach that of the neutral polaron only in the low density limit where the average interparticle spacing is larger than the characteristic length of the atom-ion interaction. We finally analyse the effects of the Fermi gas on the molecular ions, which are bound atom-dimer states.


title: Generic countably infinite groups

abstract: Countably infinite groups (with a fixed underlying set) constitute a Polish space $G$ with a suitable metric, hence the Baire category theorem holds in $G$. We study isomorphism invariant subsets of $G$, which we call group properties. We say that the generic countably infinite group is of property $P$ if $P$ is comeager in $G$. We prove that every group property with the Baire property is either meager or comeager. We show that there is a comeager elementary equivalence class in $G$ but every isomorphism class is meager. We prove that the generic group is algebraically closed, simple, not finitely generated and not locally finite. We show that in the subspace of Abelian groups the generic group is isomorphic to the unique countable, divisible torsion group that contains every finite Abelian group. We sketch the model-theoretic setting in which many of our results can be generalized. We briefly discuss a connection with infinite games.


title: On the structure of the affine asymptotic Hecke algebras

abstract: According to a conjecture of Lusztig, the asymptotic affine Hecke algebra should admit a description in terms of the Grothedieck group of sheaves on the square of a finite set equivariant under the action of the centralizer of a nilpotent element in the reductive group. A weaker form of this statement, allowing for possible central extensions of stabilizers of that action, has been proved by the first named author with Ostrik. In the present paper we describe an example showing that nontrivial central extensions do arise, thus the above weaker statement is optimal. We also show that Lusztig's homomorphism from the affine Hecke algebra to the asymptotic affine Hecke algebra induces an isomorphism on cocenters and discuss the relation of the above central extensions to the structure of the cocenter.


title: Learning Co-segmentation by Segment Swapping for Retrieval and Discovery

abstract: The goal of this work is to efficiently identify visually similar patterns in images, e.g. identifying an artwork detail copied between an engraving and an oil painting, or recognizing parts of a night-time photograph visible in its daytime counterpart. Lack of training data is a key challenge for this co-segmentation task. We present a simple yet surprisingly effective approach to overcome this difficulty: we generate synthetic training pairs by selecting segments in an image and copy-pasting them into another image. We then learn to predict the repeated region masks. We find that it is crucial to predict the correspondences as an auxiliary task and to use Poisson blending and style transfer on the training pairs to generalize on real data. We analyse results with two deep architectures relevant to our joint image analysis task: a transformer-based architecture and Sparse Nc-Net, a recent network designed to predict coarse correspondences using 4D convolutions. We show our approach provides clear improvements for artwork details retrieval on the Brueghel dataset and achieves competitive performance on two place recognition benchmarks, Tokyo247 and Pitts30K. We also demonstrate the potential of our approach for unsupervised image collection analysis by introducing a spectral graph clustering approach to object discovery and demonstrating it on the object discovery dataset of \cite{rubinstein2013unsupervised} and the Brueghel dataset. Our code and data are available at http://imagine.enpc.fr/~shenx/SegSwap/.


title: Transformer Ensembles for Sexism Detection

abstract: This document presents in detail the work done for the sexism detection task at EXIST2021 workshop. Our methodology is built on ensembles of Transformer-based models which are trained on different background and corpora and fine-tuned on the provided dataset from the EXIST2021 workshop. We report accuracy of 0.767 for the binary classification task (task1), and f1 score 0.766, and for the multi-class task (task2) accuracy 0.623 and f1-score 0.535.


title: Feshbach resonances and molecule formation in ultracold mixtures of Rb and Yb($^3$P) atoms

abstract: We have investigated magnetically tunable Feshbach resonances in ultracold collisions of Rb with Yb in its metastable $^3$P$_2$ and $^3$P$_0$ states, using coupled-channel scattering and bound-state calculations. For the $^3$P$_2$ state, we find sharp resonances when both atoms are in their lowest Zeeman sublevels. However, these resonances are decayed by inelastic processes that produce Yb atoms in $^3$P$_1$ and $^3$P$_0$ states. The molecules that might be produced by magnetoassociation at the $^3$P$_2$ thresholds can decay by similar pathways and would have lifetimes no more than a few microseconds. For the $^3$P$_0$ state, by contrast, there are resonances that are promising for magnetoassociation. There are resonances due to both rotating and non-rotating molecular states that are significantly stronger than the analogous resonances for Yb($^1$S). The ones due to rotating states are denser in magnetic field; in contrast to Yb($^1$S), they exist even for bosonic isotopes of Yb($^3$P$_0$).


title: Learning to Be Cautious

abstract: A key challenge in the field of reinforcement learning is to develop agents that behave cautiously in novel situations. It is generally impossible to anticipate all situations that an autonomous system may face or what behavior would best avoid bad outcomes. An agent that could learn to be cautious would overcome this challenge by discovering for itself when and how to behave cautiously. In contrast, current approaches typically embed task-specific safety information or explicit cautious behaviors into the system, which is error-prone and imposes extra burdens on practitioners. In this paper, we present both a sequence of tasks where cautious behavior becomes increasingly non-obvious, as well as an algorithm to demonstrate that it is possible for a system to \emph{learn} to be cautious. The essential features of our algorithm are that it characterizes reward function uncertainty without task-specific safety information and uses this uncertainty to construct a robust policy. Specifically, we construct robust policies with a $k$-of-$N$ counterfactual regret minimization (CFR) subroutine given a learned reward function uncertainty represented by a neural network ensemble belief. These policies exhibit caution in each of our tasks without any task-specific safety tuning.


title: Geometry of Smooth Extremal Surfaces

abstract: We study the geometry of the smooth projective surfaces that are defined by Frobenius forms, a class of homogenous polynomials in prime characteristic recently shown to have minimal possible F-pure threshold among forms of the same degree. We call these surfaces $\textit{extremal surfaces}$, and show that their geometry is reminiscent of the geometry of smooth cubic surfaces, especially non-Frobenius split cubic surfaces of characteristic two, which are examples of extremal surfaces. For example, we show that an extremal surface $X$ contains $d^2(d^2-3d+3)$ lines where $d$ is the degree, which is notable since the number of lines on a complex surface is bounded above by a quadratic function in $d$. Whenever two of those lines meet, they determine a $d$-tangent plane to $X$ which consists of a union of $d$ lines meeting in one point; we count the precise number of such "star points" on $X$, showing that it is quintic in the degree, which recovers the fact that there are exactly 45 Eckardt points on an extremal cubic surface. Finally, we generalize the classical notion of a double six for cubic surfaces to a double $2d$ on an extremal surface of degree $d$. We show that, asymptotically in $d$, smooth extremal surfaces have at least $\frac{1}{16}d^{14}$ double $2d$'s. A key element of the proofs is using the large automorphism group of extremal surfaces which we show acts transitively on many sets, such as the set of (triples of skew) lines on the extremal surface. Extremal surfaces are closely related to finite Hermitian geometries, which we recover as the $\mathbb F_{q^2}$-rational points of special extremal surfaces defined by Hermitian forms over $\mathbb F_{q^2}$.


title: Contrastive prediction strategies for unsupervised segmentation and categorization of phonemes and words

abstract: We investigate the performance on phoneme categorization and phoneme and word segmentation of several self-supervised learning (SSL) methods based on Contrastive Predictive Coding (CPC). Our experiments show that with the existing algorithms there is a trade off between categorization and segmentation performance. We investigate the source of this conflict and conclude that the use of context building networks, albeit necessary for superior performance on categorization tasks, harms segmentation performance by causing a temporal shift on the learned representations. Aiming to bridge this gap, we take inspiration from the leading approach on segmentation, which simultaneously models the speech signal at the frame and phoneme level, and incorporate multi-level modelling into Aligned CPC (ACPC), a variation of CPC which exhibits the best performance on categorization tasks. Our multi-level ACPC (mACPC) improves in all categorization metrics and achieves state-of-the-art performance in word segmentation.


title: Context-independent mapping and free choice are equivalent

abstract: Free choice (or statistical independence) assumption in a hidden variable model (HVM) means that the settings chosen by experimenters do not depend on the values of the hidden variable. The assumption of context-independent (CI) mapping in an HVM means that the results of a measurement do not depend on settings for other measurements. If the measurements are spacelike separated, this assumption is known as local causality. Both free choice and CI mapping assumptions are considered necessary for derivation of the Bell-type criteria of contextuality/nonlocality. It is known, however, for a variety of special cases, that the two assumptions are not logically independent. We show here, in complete generality, for any system of random variables with or without disturbance/signaling, that an HVM that postulates CI mapping is equivalent to an HVM that postulates free choice. If one denies the possibility that a given empirical scenario can be described by an HVM in which measurements depend on other measurements' settings, free choice violations should be denied too, and vice versa. KEYWORDS: Contextuality; context-independent mapping; free choice; local causality; nonlocality.


title: Physics-informed linear regression is competitive with two Machine Learning methods in residential building MPC

abstract: Because physics-based building models are difficult to obtain as each building is individual, there is an increasing interest in generating models suitable for building MPC directly from measurement data. Machine learning methods have been widely applied to this problem and validated mostly in simulation; there are, however, few studies on a direct comparison of different models or validation in real buildings to be found in the literature. Methods that are indeed validated in application often lead to computationally complex non-convex optimization problems. Here we compare physics-informed Autoregressive-Moving-Average with Exogenous Inputs (ARMAX) models to Machine Learning models based on Random Forests and Input Convex Neural Networks and the resulting convex MPC schemes in experiments on a practical building application with the goal of minimizing energy consumption while maintaining occupant comfort, and in a numerical case study. We demonstrate that Predictive Control in general leads to savings between 26% and 49% of heating and cooling energy, compared to the building's baseline hysteresis controller. Moreover, we show that all model types lead to satisfactory control performance in terms of constraint satisfaction and energy reduction. However, we also see that the physics-informed ARMAX models have a lower computational burden, and a superior sample efficiency compared to the Machine Learning based models. Moreover, even if abundant training data is available, the ARMAX models have a significantly lower prediction error than the Machine Learning models, which indicates that the encoded physics-based prior of the former cannot independently be found by the latter.


title: On the use of uncertainty in classifying Aedes Albopictus mosquitoes

abstract: The re-emergence of mosquito-borne diseases (MBDs), which kill hundreds of thousands of people each year, has been attributed to increased human population, migration, and environmental changes. Convolutional neural networks (CNNs) have been used by several studies to recognise mosquitoes in images provided by projects such as Mosquito Alert to assist entomologists in identifying, monitoring, and managing MBD. Nonetheless, utilising CNNs to automatically label input samples could involve incorrect predictions, which may mislead future epidemiological studies. Furthermore, CNNs require large numbers of manually annotated data. In order to address the mentioned issues, this paper proposes using the Monte Carlo Dropout method to estimate the uncertainty scores in order to rank the classified samples to reduce the need for human supervision in recognising Aedes albopictus mosquitoes. The estimated uncertainty was also used in an active learning framework, where just a portion of the data from large training sets was manually labelled. The experimental results show that the proposed classification method with rejection outperforms the competing methods by improving overall performance and reducing entomologist annotation workload. We also provide explainable visualisations of the different regions that contribute to a set of samples' uncertainty assessment.


title: The Jacobi operator and its Donoghue $m$-functions

abstract: In this paper we construct Donoghue $m$-functions for the Jacobi differential operator in $L^2\big((-1,1); (1-x)^{\alpha} (1+x)^{\beta} dx\big)$, associated to the differential expression \begin{align*} \begin{split} \tau_{\alpha,\beta} = - (1-x)^{-\alpha} (1+x)^{-\beta}(d/dx) \big((1-x)^{\alpha + 1}(1+x)^{\beta + 1}\big) (d/dx),& \\ x \in (-1,1), \; \alpha, \beta \in \mathbb{R}, \end{split} \end{align*} whenever at least one endpoint, $x=\pm 1$, is in the limit circle case. In doing so, we provide a full treatment of the Jacobi operator's $m$-functions corresponding to coupled boundary conditions whenever both endpoints are in the limit circle case, a topic not covered in the literature.


title: Improving the quality of generative models through Smirnov transformation

abstract: Solving the convergence issues of Generative Adversarial Networks (GANs) is one of the most outstanding problems in generative models. In this work, we propose a novel activation function to be used as output of the generator agent. This activation function is based on the Smirnov probabilistic transformation and it is specifically designed to improve the quality of the generated data. In sharp contrast with previous works, our activation function provides a more general approach that deals not only with the replication of categorical variables but with any type of data distribution (continuous or discrete). Moreover, our activation function is derivable and therefore, it can be seamlessly integrated in the backpropagation computations during the GAN training processes. To validate this approach, we evaluate our proposal against two different data sets: a) an artificially rendered data set containing a mixture of discrete and continuous variables, and b) a real data set of flow-based network traffic data containing both normal connections and cryptomining attacks. To evaluate the fidelity of the generated data, we analyze both their results in terms of quality measures of statistical nature and also regarding the use of these synthetic data to feed a nested machine learning-based classifier. The experimental results evince a clear outperformance of the GAN network tuned with this new activation function with respect to both a na\"ive mean-based generator and a standard GAN. The quality of the data is so high that the generated data can fully substitute real data for training the nested classifier without a fall in the obtained accuracy. This result encourages the use of GANs to produce high-quality synthetic data that are applicable in scenarios in which data privacy must be guaranteed.


title: Frequency conversion of abruptly autofocusing waves

abstract: Abruptly autofocusing waves and associated ring-Airy (RA) beams are attracting increasing interest owing to their fascinating properties such as their ability of abruptly autofocusing to small F-number. Optical frequency conversion via nonlinear interactions can further expand their applications to new area, yet are rarely studied. In this work, we report the frequency conversion of RA beams via sum-frequency generation using perfect flattop and common Gauss beams as the pump beams. The nonlinear transformation of the spatial complex amplitude of the signal and associated influences on autofocusing behavior, under different conditions of interaction location (i.e., original, autofocusing, and Fourier planes) and pump structure, were systematically studied and experimentally investigated. This proof-of principle demonstration provides a general guideline to build the frequency interface for abruptly autofocusing waves and a reference for relevant studies involving nonlinear transformation of abruptly autofocusing waves.


title: Spinons and damped phonons in spin-1/2 quantum-liquid Ba$_{4}$Ir${}_3$O${}_{10}$ observed by Raman scattering

abstract: In spin-1/2 Mott insulators, non-magnetic quantum liquid phases are often argued to arise when the system shows no magnetic ordering, but identifying positive signatures of these phases or related spinon quasiparticles can be elusive. Here we use Raman scattering to provide three signatures for spinons in a possible spin-orbit quantum liquid material Ba${}_4$Ir${}_3$O${}_{10}$: (1) A broad hump, which we show can arise from Luttinger Liquid spinons in Raman with parallel photon polarizations normal to 1D chains; (2) Strong phonon damping from phonon-spin coupling via the spin-orbit interaction; and (3) the absence of (1) and (2) in the magnetically ordered phase that is produced when 2% of Ba is substituted by Sr ((Ba${}_{0.98}$Sr${}_{0.02}$)${}_4$Ir${}_3$O${}_{10}$). The phonon damping via itinerant spinons seen in this quantum-liquid insulator suggests a new mechanism for enhancing thermoelectricity in strongly correlated conductors, through a neutral quantum liquid that need not affect electronic transport.


title: Non-equilibrium fluctuations and nonlinear response of an active bath

abstract: We analyze the dynamics of a passive colloidal probe immersed in an active bath using an optical trap to study three physical processes: (1) the non-equilibrium fluctuations transferred to the probe by the active bath, (2) the friction experienced by the probe as it is driven through the active bath, and (3) the force relaxation of the probe returning to its equilibrium position. We measure the local force dynamics where all of the following characteristics are of $\mathcal{O}(1)$: the size of the probe colloid relative to active bath particle; the size of the probe colloid relative to the characteristic run-length of an active particle; and the timescale of probe movement to the persistence time of an active particle. We find at P\'{e}clet (Pe) $\ll 1$ the active suspension exhibits shear thinning down to the solvent viscosity (but not below); at $0.85 <$ Pe $\leq 5.1$ the active bath shear thickens; and at Pe $\geq 8.5$ the effective viscosity of the active bath shows a decreased effect of thickening and plateaus. These results are in agreement with recent modeling and simulations of the nonlinear rheology of an isotropic active bath, providing experimental verification, and suggesting the model predictions extends to moderately dense suspensions. Further, we observe that the distribution of force fluctuations depends on Pe, unlike in passive equilibrium baths. Lastly, we measure the energy transfer rate from the active bath to the probe to be $\langle J \rangle \approx 10^3$ $k_B T/$s, which leads to an increase in the effective diffusion of the probe by a factor of $\sim 2$.


title: Takagi factorization of matrices depending on parameters and locating degeneracies of singular values

abstract: In this work we consider the Takagi factorization of a matrix valued function depending on parameters. We give smoothness and genericity results and pay particular attention to the concerns caused by having either a singular value equal to $0$ or multiple singular values. For these phenomena, we give theoretical results showing that their co-dimension is $2$, and we further develop and test numerical methods to locate in parameter space values where these occurrences take place. Numerical study of the density of these occurrences is performed.


title: Performance Analysis of Dual-Hop THz Wireless Transmission for Backhaul Applications

abstract: THz transmissions suffer from pointing errors due to antenna misalignment and incur higher path loss from the molecular absorption in addition to the channel fading. In this paper, we employ an amplify-and-forward (AF) dual-hop relaying to mitigate the effect of pointing errors and extend the range of the THz wireless system for backhaul connectivity. We provide statistical analysis on the performance of the considered system by deriving analytical expressions for the outage probability, average bit-error-rate (BER), average signal-to-noise ratio (SNR), and a lower bound on the ergodic capacity over independent and identical (i.i.d) $\alpha$-$\mu$ fading combined with the statistical effect of pointing errors. Using computer simulations, we validate the derived analysis of the relay-assisted system. We also demonstrate the effect of the system parameters on outage probability and average BER with the help of diversity order. We show that data rates up to several \mbox{Gbps} can be achieved using THz transmissions, which is desirable for next-generation wireless systems, especially for backhaul applications.


title: Entropy Stable Discontinuous Galerkin Methods for Balance Laws in Non-Conservative Form: Applications to the Euler Equations with Gravity

abstract: In this work a non-conservative balance law formulation is considered that encompasses the rotating, compressible Euler equations for dry atmospheric flows. We develop a semi-discretely entropy stable discontinuous Galerkin method on curvilinear meshes using a generalization of flux differencing for numerical fluxes in fluctuation form. The method uses the skew-hybridized formulation of the element operators to ensure that, even in the presence of under-integration on curvilinear meshes, the resulting discretization is entropy stable. Several atmospheric flow test cases in one, two, and three dimensions confirm the theoretical entropy stability results as well as show the high-order accuracy and robustness of the method.


title: Good labeling property of simple nested fractals

abstract: We show various criteria to verify if a given nested fractal has a good labeling property, inter alia we present a characterization of GLP for fractals with an odd number of essential fixed points. We show a convenient reduction of area to be investigated in verification of GLP and give examples that further reduction is impossible. We prove that if a number of essential fixed points is a power of two, then a fractal must have GLP and that there are no values other than primes or powers of two guaranteeing GLP. For all other numbers of essential fixed points we are able to construct examples having and other not having GLP.


title: High Precision Ringdown Modeling: Multimode Fits and BMS Frames

abstract: Quasi-normal mode (QNM) modeling is an invaluable tool for characterizing remnant black holes, studying strong gravity, and testing GR. Only recently have QNM studies begun to focus on multimode fitting to numerical relativity (NR) strain waveforms. As GW observatories become even more sensitive they will be able to resolve higher-order modes. Consequently, multimode QNM fits will be critically important, and in turn require a more thorough treatment of the asymptotic frame at $\mathscr{I}^+$. The first main result of this work is a method for systematically fitting a QNM model containing many modes to a numerical waveform produced using Cauchy-characteristic extraction (CCE), an extraction technique which is known to resolve memory effects. We choose the modes to model based on their power contribution to the residual between numerical and model waveforms. We show that the all-mode strain mismatch improves by a factor of $\sim10^5$ when using multimode fitting as opposed to only fitting the $(2,\pm2,n)$ modes. Our most significant result addresses a critical point that has been overlooked in the QNM literature: the importance of matching the Bondi-van der Burg-Metzner-Sachs (BMS) frame of the numerical waveform to that of the QNM model. We show that by mapping the numerical waveforms$-$which exhibit the memory effect$-$to a BMS frame known as the super rest frame, there is an improvement of $\sim10^5$ in the all-mode strain mismatch compared to using a strain waveform whose BMS frame is not fixed. Furthermore, we find that by mapping CCE waveforms to the super rest frame, we can obtain all-mode mismatches that are, on average, a factor of $\sim4$ better than using the publicly-available extrapolated waveforms. We illustrate the effectiveness of these modeling enhancements by applying them to families of waveforms produced by NR and comparing our results to previous QNM studies.


title: Efficient Representation of Interaction Patterns with Hyperbolic Hierarchical Clustering for Classification of Users on Twitter

abstract: Social media platforms play an important role in democratic processes. During the 2019 General Elections of India, political parties and politicians widely used Twitter to share their ideals, advocate their agenda and gain popularity. Twitter served as a ground for journalists, politicians and voters to interact. The organic nature of these interactions can be upended by malicious accounts on Twitter, which end up being suspended or deleted from the platform. Such accounts aim to modify the reach of content by inorganically interacting with particular handles. These interactions are a threat to the integrity of the platform, as such activity has the potential to affect entire results of democratic processes. In this work, we design a feature extraction framework which compactly captures potentially insidious interaction patterns. Our proposed features are designed to bring out communities amongst the users that work to boost the content of particular accounts. We use Hyperbolic Hierarchical Clustering (HypHC) which represents the features in the hyperbolic manifold to further separate such communities. HypHC gives the added benefit of representing these features in a lower dimensional space -- thus serving as a dimensionality reduction technique. We use these features to distinguish between different classes of users that emerged in the aftermath of the 2019 General Elections of India. Amongst the users active on Twitter during the elections, 2.8% of the users participating were suspended and 1% of the users were deleted from the platform. We demonstrate the effectiveness of our proposed features in differentiating between regular users (users who were neither suspended nor deleted), suspended users and deleted users. By leveraging HypHC in our pipeline, we obtain F1 scores of upto 93%.


title: Seismic response of cylinder assemblies in axial flow

abstract: Earthquakes are a great challenge for the safety of nuclear reactors. To address this challenge, we need to better understand how the reactor core responds to seismic forcing. The reactor core is made of fuel assemblies, which are themselves composed of flexible fuel rods immersed in a strong axial flow. This gives rise to strongly-coupled fluid-structure interactions whose accurate modelling generally requires high computational costs. In this paper, we introduce a new model able to capture the mechanical response of the reactor core subjected to seismic forcing with low computational costs. This model is based on potential flow theory for the fluid part and Euler-Bernoulli beam theory for the structural part allowing us to predict the response to seismic forcing in presence of axial flow.. The linear equations are solved in the Fourier space to decrease computational time. For validation purposes, we first use the proposed model to compute the response of a single cylinder in axial flow. We then implement a multiple cylinder geometry made of 4 fuel assemblies, each made of 8x8 cylinders, corresponding to an experimental facility available at CEA. The comparison between numerical results and experiments show good agreement. The model can correctly predict the added mass. It can also qualitatively capture the coupling between assemblies and the effect of confinement. This shows that a potential flow approach can give insight into the complex fluid-structure interactions within a nuclear reactor and, in particular, be used to predict the response to seismic forcing at low computational cost.


title: Finite Callan-Symanzik renormalisation for multiple scalar fields

abstract: We study a finite, divergence free approach to renormalisation originally proposed in the early '70s by Blaer and Young, and Callan. It is based on equations similar to the Callan-Symanzik equations, and introduced in the context of $\lambda \phi^4$ theory. We generalise this method to the case of two interacting scalar fields, with obvious generalisation to an arbitrary number of fields.


title: Delayed Propagation Transformer: A Universal Computation Engine towards Practical Control in Cyber-Physical Systems

abstract: Multi-agent control is a central theme in the Cyber-Physical Systems (CPS). However, current control methods either receive non-Markovian states due to insufficient sensing and decentralized design, or suffer from poor convergence. This paper presents the Delayed Propagation Transformer (DePT), a new transformer-based model that specializes in the global modeling of CPS while taking into account the immutable constraints from the physical world. DePT induces a cone-shaped spatial-temporal attention prior, which injects the information propagation and aggregation principles and enables a global view. With physical constraint inductive bias baked into its design, our DePT is ready to plug and play for a broad class of multi-agent systems. The experimental results on one of the most challenging CPS -- network-scale traffic signal control system in the open world -- show that our model outperformed the state-of-the-art expert methods on synthetic and real-world datasets. Our codes are released at: https://github.com/VITA-Group/DePT.