title: On the Optimality and Convergence Properties of the Iterative Learning Model Predictive Controller

abstract: In this technical note we analyse the performance improvement and optimality properties of the Learning Model Predictive Control (LMPC) strategy for linear deterministic systems. The LMPC framework is a policy iteration scheme where closed-loop trajectories are used to update the control policy for the next execution of the control task. We show that, when a Linear Independence Constraint Qualification (LICQ) condition holds, the LMPC scheme guarantees strict iterative performance improvement and optimality, meaning that the closed-loop cost evaluated over the entire task converges asymptotically to the optimal cost of the infinite-horizon control problem. Compared to previous works this sufficient LICQ condition can be easily checked, it holds for a larger class of systems and it can be used to adaptively select the prediction horizon of the controller, as demonstrated by a numerical example.


title: Spin-valley collective modes of the electron liquid in graphene

abstract: We develop the theory of collective modes supported by a Fermi liquid of electrons in pristine graphene. Under reasonable assumptions regarding the electron-electron interaction, all the modes but the plasmon are over-damped. In addition to the $SU(2)$ symmetric spin mode, these include also the valley imbalance modes obeying a $U(1)$ symmetry, and a $U(2)$ symmetric valley spin imbalance mode. We derive the interactions and diffusion constants characterizing the over-damped modes. The corresponding relaxation rates set fundamental constraints on graphene valley- and spintronics applications.


title: Kernel Aggregated Fast Multipole Method: Efficient summation of Laplace and Stokes kernel functions

abstract: Many different simulation methods for Stokes flow problems involve a common computationally intense task -- the summation of a kernel function over $O(N^2)$ pairs of points. One popular technique is the Kernel Independent Fast Multipole Method (KIFMM), which constructs a spatial adaptive octree for all points and places a small number of equivalent multipole and local equivalent points around each octree box, and completes the kernel sum with $O(N)$ cost, using these equivalent points. Simpler kernels can be used between these equivalent points to improve the efficiency of KIFMM. Here we present further extensions and applications to this idea, to enable efficient summations and flexible boundary conditions for various kernels. We call our method the Kernel Aggregated Fast Multipole Method (KAFMM), because it uses different kernel functions at different stages of octree traversal. We have implemented our method as an open-source software library STKFMM based on the high performance library PVFMM, with support for Laplace kernels, the Stokeslet, regularized Stokeslet, Rotne-Prager-Yamakawa (RPY) tensor, and the Stokes double-layer and traction operators. Open and periodic boundary conditions are supported for all kernels, and the no-slip wall boundary condition is supported for the Stokeslet and RPY tensor. The package is designed to be ready-to-use as well as being readily extensible to additional kernels.

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


title: Diagnostic data integration using deep neural networks for real-time plasma analysis

abstract: Recent advances in acquisition equipment is providing experiments with growing amounts of precise yet affordable sensors. At the same time an improved computational power, coming from new hardware resources (GPU, FPGA, ACAP), has been made available at relatively low costs. This led us to explore the possibility of completely renewing the chain of acquisition for a fusion experiment, where many high-rate sources of data, coming from different diagnostics, can be combined in a wide framework of algorithms. If on one hand adding new data sources with different diagnostics enriches our knowledge about physical aspects, on the other hand the dimensions of the overall model grow, making relations among variables more and more opaque. A new approach for the integration of such heterogeneous diagnostics, based on composition of deep variational autoencoders, could ease this problem, acting as a structural sparse regularizer. This has been applied to RFX-mod experiment data, integrating the soft X-ray linear images of plasma temperature with the magnetic state. However to ensure a real-time signal analysis, those algorithmic techniques must be adapted to run in well suited hardware. In particular it is shown that, attempting a quantization of neurons transfer functions, such models can be modified to create an embedded firmware. This firmware, approximating the deep inference model to a set of simple operations, fits well with the simple logic units that are largely abundant in FPGAs. This is the key factor that permits the use of affordable hardware with complex deep neural topology and operates them in real-time.


title: Panoster: End-to-end Panoptic Segmentation of LiDAR Point Clouds

abstract: Panoptic segmentation has recently unified semantic and instance segmentation, previously addressed separately, thus taking a step further towards creating more comprehensive and efficient perception systems. In this paper, we present Panoster, a novel proposal-free panoptic segmentation method for LiDAR point clouds. Unlike previous approaches relying on several steps to group pixels or points into objects, Panoster proposes a simplified framework incorporating a learning-based clustering solution to identify instances. At inference time, this acts as a class-agnostic segmentation, allowing Panoster to be fast, while outperforming prior methods in terms of accuracy. Without any post-processing, Panoster reached state-of-the-art results among published approaches on the challenging SemanticKITTI benchmark, and further increased its lead by exploiting heuristic techniques. Additionally, we showcase how our method can be flexibly and effectively applied on diverse existing semantic architectures to deliver panoptic predictions.


title: CNN Profiler on Polar Coordinate Images for Tropical Cyclone Structure Analysis

abstract: Convolutional neural networks (CNN) have achieved great success in analyzing tropical cyclones (TC) with satellite images in several tasks, such as TC intensity estimation. In contrast, TC structure, which is conventionally described by a few parameters estimated subjectively by meteorology specialists, is still hard to be profiled objectively and routinely. This study applies CNN on satellite images to create the entire TC structure profiles, covering all the structural parameters. By utilizing the meteorological domain knowledge to construct TC wind profiles based on historical structure parameters, we provide valuable labels for training in our newly released benchmark dataset. With such a dataset, we hope to attract more attention to this crucial issue among data scientists. Meanwhile, a baseline is established with a specialized convolutional model operating on polar-coordinates. We discovered that it is more feasible and physically reasonable to extract structural information on polar-coordinates, instead of Cartesian coordinates, according to a TC's rotational and spiral natures. Experimental results on the released benchmark dataset verified the robustness of the proposed model and demonstrated the potential for applying deep learning techniques for this barely developed yet important topic.


title: Non-abelian Thouless pumping in a photonic lattice

abstract: Non-abelian gauge fields emerge naturally in the description of adiabatically evolving quantum systems having degenerate levels. Here we show that they also play a role in Thouless pumping in the presence of degenerate bands. To this end we consider a photonic Lieb lattice having two degenerate non-dispersive modes and we show that, when the lattice parameters are slowly modulated, the propagation of the photons bear the fingerprints of the underlying non-abelian gauge structure. The non-dispersive character of the bands enables a high degree of control on photon propagation. Our work paves the way to the generation and detection of non-abelian gauge fields in photonic and optical lattices.


title: On $BT_1$ group schemes and Fermat Jacobians

abstract: Let $p$ be a prime number and let $k$ be an algebraically closed field of characteristic $p$. A $BT_1$ group scheme over $k$ is a finite commutative group scheme which arises as the kernel of $p$ on a $p$-divisible (Barsotti--Tate) group. We compare three classifications of $BT_1$ group schemes, due in large part to Kraft, Ekedahl, and Oort, and defined using words, canonical filtrations, and permutations. Using this comparison, we determine the Ekedahl--Oort types of Fermat quotient curves and we compute four invariants of the $p$-torsion group schemes of these curves.


title: Efimov effect at the Kardar-Parisi-Zhang roughening transition

abstract: Surface growth governed by the Kardar-Parisi-Zhang (KPZ) equation in dimensions higher than two undergoes a roughening transition from smooth to rough phases with increasing the nonlinearity. It is also known that the KPZ equation can be mapped onto quantum mechanics of attractive bosons with a contact interaction, where the roughening transition corresponds to a binding transition of two bosons with increasing the attraction. Such critical bosons in three dimensions actually exhibit the Efimov effect, where a three-boson coupling turns out to be relevant under the renormalization group so as to break the scale invariance down to a discrete one. On the basis of these facts linking the two distinct subjects in physics, we predict that the KPZ roughening transition in three dimensions shows either the discrete scale invariance or no intrinsic scale invariance.


title: Sizeless: Predicting the optimal size of serverless functions

abstract: Serverless functions are a cloud computing paradigm where the provider takes care of resource management tasks such as resource provisioning, deployment, and auto-scaling. The only resource management task that developers are still in charge of is selecting how much resources are allocated to each worker instance. However, selecting the optimal size of serverless functions is quite challenging, so developers often neglect it despite its significant cost and performance benefits. Existing approaches aiming to automate serverless functions resource sizing require dedicated performance tests, which are time-consuming to implement and maintain. In this paper, we introduce an approach to predict the optimal resource size of a serverless function using monitoring data from a single resource size. As our approach does not require dedicated performance tests, it enables cloud providers to implement resource sizing on a platform level and automate the last resource management task associated with serverless functions. We evaluate our approach on three different serverless applications, where it selects the optimal memory size for 71.7% of the serverless functions and the second-best memory size for 22.3% of the serverless functions, which results in an average speedup of 43.6% while simultaneously decreasing average costs by 10.2%.


title: Invariant chains in algebra and discrete geometry

abstract: We relate finite generation of cones, monoids, and ideals in increasing chains (the local situation) to equivariant finite generation of the corresponding limit objects (the global situation). For cones and monoids there is no analog of Noetherianity as in the case of ideals and we demonstrate this in examples. As a remedy, we find local-global correspondences for finite generation. These results are derived from a more general framework that relates finite generation under closure operations to equivariant finite generation under general families of maps. We also give a new proof that non-saturated Inc-invariant chains of ideals stabilize, closing a gap in the literature.


title: Science Extraction from TESS Observations of Known Exoplanet Hosts

abstract: The transit method of exoplanet discovery and characterization has enabled numerous breakthroughs in exoplanetary science. These include measurements of planetary radii, mass-radius relationships, stellar obliquities, bulk density constraints on interior models, and transmission spectroscopy as a means to study planetary atmospheres. The Transiting Exoplanet Survey Satellite (TESS) has added to the exoplanet inventory by observing a significant fraction of the celestial sphere, including many stars already known to host exoplanets. Here we describe the science extraction from TESS observations of known exoplanet hosts during the primary mission. These include transit detection of known exoplanets, discovery of additional exoplanets, detection of phase signatures and secondary eclipses, transit ephemeris refinement, and asteroseismology as a means to improve stellar and planetary parameters. We provide the statistics of TESS known host observations during Cycle 1 & 2, and present several examples of TESS photometry for known host stars observed with a long baseline. We outline the major discoveries from observations of known hosts during the primary mission. Finally, we describe the case for further observations of known exoplanet hosts during the TESS extended mission and the expected science yield.


title: The public debt multiplier

abstract: We study the effects on economic activity of a pure temporary change in government debt and the relationship between the debt multiplier and the level of debt in an overlapping generations framework. The debt multiplier is positive but quite small during normal times while it is much larger during crises. Moreover, it increases with the steady state level of debt. Hence, the call for fiscal consolidation during recessions seems ill-advised. Finally, a rise in the steady state debt-to-GDP level increases the steady state real interest rate providing more room for manoeuvre to monetary policy to fight deflationary shocks.


title: Polymer Informatics with Multi-Task Learning

abstract: Modern data-driven tools are transforming application-specific polymer development cycles. Surrogate models that can be trained to predict the properties of new polymers are becoming commonplace. Nevertheless, these models do not utilize the full breadth of the knowledge available in datasets, which are oftentimes sparse; inherent correlations between different property datasets are disregarded. Here, we demonstrate the potency of multi-task learning approaches that exploit such inherent correlations effectively, particularly when some property dataset sizes are small. Data pertaining to 36 different properties of over $13, 000$ polymers (corresponding to over $23,000$ data points) are coalesced and supplied to deep-learning multi-task architectures. Compared to conventional single-task learning models (that are trained on individual property datasets independently), the multi-task approach is accurate, efficient, scalable, and amenable to transfer learning as more data on the same or different properties become available. Moreover, these models are interpretable. Chemical rules, that explain how certain features control trends in specific property values, emerge from the present work, paving the way for the rational design of application specific polymers meeting desired property or performance objectives.


title: Observation of an accidental bound state in the continuum in a chain of dielectric disks

abstract: Being a general wave phenomenon, bound states in the continuum (BICs) appear in acoustic, hydrodynamic, and photonic systems of various dimensionalities. Here, we report the first experimental observation of an accidental electromagnetic BIC in a one-dimensional periodic chain of coaxial ceramic disks. We show that the accidental BIC manifests itself as a narrow peak in the transmission spectra of the chain placed between two loop antennas. We demonstrate a linear growth of the radiative quality factor of the BICs with the number of disks that is well-described with a tight-binding model. We estimate the number of the disks when the radiation losses become negligible in comparison to material absorption and, therefore, the chain can be considered practically as infinite. The presented analysis is supported by near-field measurements of the BIC profile. The obtained results provide useful guidelines for practical implementations of structures with BICs opening new horizons for the development of radio-frequency and optical metadevices.


title: Unconstraint minimization of continuous convex functions. Application to LP

abstract: Our contribution in this paper is two folded. We consider first the case of linear programming with real coefficients and give a method which allows the computation of a new upper bound on the distance from the origin to a feasible point. Next we present an application of the ellipsoid method to form a novel algorithm which allows one to search for the continuous convex functions minimum without the prior knowledge on the search radius. If stopped early, the proposed algorithm can give proofs that the optimal value has not been reached yet, hence the user can opt for more iterations. However, if prior guarantees exist on the existence of an optimum point in a given ball, the algorithm is guaranteed to find it. For such a case we prove a polynomial upper bound on the number of flops required to obtain the solution. The presented algorithm is then applied to linear programming. We further develop our algorithm and provide a method to answer the feasibility question of linear programs with real coefficients in a number of flops bounded above by a polynomial in the number of variables and the number of constraints. However, in case of feasibility, our method does not generate an actual feasible point. For obtaining such a point, we have to make some assumptions on the subgradients of a certain function.


title: Semi-Grant-Free NOMA: Ergodic Rates Analysis with Random Deployed Users

abstract: Semi-grant-free (Semi-GF) non-orthogonal multiple access (NOMA) enables grant-free (GF) and grant-based (GB) users to share the same resource blocks, thereby balancing the connectivity and stability of communications. This letter analyzes ergodic rates of Semi-GF NOMA systems. First, this paper exploits a Semi-GF protocol, denoted as dynamic protocol, for selecting GF users into the occupied GB channels via the GB user's instantaneous received power. Under this protocol, the closed-form analytical and approximated expressions for ergodic rates are derived. The numerical results illustrate that the GF user (weak NOMA user) has a performance upper limit, while the ergodic rate of the GB user (strong NOMA user) increases linearly versus the transmit signal-to-noise ratio.


title: Semi-abelian analogues of Schanuel Conjecture and applications

abstract: In this article we study Semi-abelian analogues of Schanuel conjecture. As showed by the first author, Schanuel Conjecture is equivalent to the Generalized Period Conjecture applied to 1-motives without abelian part. Extending her methods, the second, the third and the fourth authors have introduced the Abelian analogue of Schanuel Conjecture as the Generalized Period Conjecture applied to 1-motives without toric part. As a first result of this paper, we define the Semi-abelian analogue of Schanuel Conjecture as the Generalized Period Conjecture applied to 1-motives. C. Cheng et al. proved that Schanuel conjecture implies the algebraic independence of the values of the iterated exponential and the values of the iterated logarithm, answering a question of M. Waldschmidt. The second, the third and the fourth authors have investigated a similar question in the setup of abelian varieties: the Weak Abelian Schanuel conjecture implies the algebraic independence of the values of the iterated abelian exponential and the values of an iterated generalized abelian logarithm. The main result of this paper is that a Relative Semi-abelian conjecture implies the algebraic independence of the values of the iterated semi-abelian exponential and the values of an iterated generalized semi-abelian logarithm.


title: Slicing a single wireless collision channel among throughput- and timeliness-sensitive services

abstract: The fifth generation (5G) wireless system has a platform-driven approach, aiming to support heterogeneous connections with very diverse requirements. The shared wireless resources should be sliced in a way that each user perceives that its requirement has been met. Heterogeneity challenges the traditional notion of resource efficiency, as the resource usage has cater for, e.g. rate maximization for one user and timeliness requirement for another user. This paper treats a model for radio access network (RAN) uplink, where a throughput-demanding broadband user shares wireless resources with an intermittently active user that wants to optimize the timeliness, expressed in terms of latency-reliability or Age of Information (AoI). We evaluate the trade-offs between throughput and timeliness for Orthogonal Multiple Access (OMA) as well as Non-Orthogonal Multiple Access (NOMA) with successive interference cancellation (SIC). We observe that NOMA with SIC, in a conservative scenario with destructive collisions, is just slightly inferior to that of OMA, which indicates that it may offer significant benefits in practical deployments where the capture effect is frequently encountered. On the other hand, finding the optimal configuration of NOMA with SIC depends on the activity pattern of the intermittent user, to which OMA is insensitive.


title: Effect of two parallel intruders on net work during granular penetrations

abstract: The impact of single passive intruders into granular particles has been studied in detail. However, the impact force produced by multiple intruders separated at a distance from one another, and hence the effect of their presence in close proximity to one another, is largely unexplored. Here, we use numerical simulations and laboratory experiments to study the force response of two parallel rods intruding vertically into granular media while varying the gap spacing between them. We also explored the effect of variations in friction, intruder size, and particle size on the force response. The net work ($W$) of the two rods over the depth of intrusion was measured, and the instantaneous velocities of particles over the duration of intrusion were calculated by simulations. We found that the work done by the intruders changes with distance between them. We observed a peak in $W$ at a gap spacing of $\sim$3 particle diameters, which was up to 25% greater than $W$ at large separation (\textgreater 11 particle diameters), beyond which the net work plateaued. This peak was likely due to less particle flow between intruders as we found a larger number of strong forces---identified as force chains---in the particle domain at gaps surrounding the peak force. Although higher friction caused greater force generation during intrusion, the gap spacing between the intruders at which the peak work was generated remained unchanged. Larger intruder sizes resulted in greater net work with the peak in $W$ occurring at slightly larger intruder separations. Taken together, our results show that peak work done by two parallel intruders remained within a narrow range, remaining robust to most other tested parameters.


title: Predictions for Strong Lens Detections with the Nancy Grace Roman Space Telescope

abstract: Strong gravitational lensing is an ideal tool for mapping the distribution of dark matter and for testing the values of cosmological parameters. Forthcoming surveys such as Euclid and the Rubin Observatory LSST should increase the number of known strong lensing systems significantly, with for example over 100,000 such systems predicted in the Euclid wide survey. In this short research note we predict that approximately 17,000 strong gravitational lenses will also be detectible in the Nancy Grace Roman Space Telescope 2000 square degree survey, using the LensPop gravitational lensing model. We present predicted distributions in source and deflector redshifts, magnitudes, and magnifications. Although this survey is not primarily designed as a strong lensing detection experiment, it will still provide a large complementary catalogue to shallower, wider-area forthcoming lensing discovery projects.


title: Improving Perceptual Quality by Phone-Fortified Perceptual Loss using Wasserstein Distance for Speech Enhancement

abstract: Speech enhancement (SE) aims to improve speech quality and intelligibility, which are both related to a smooth transition in speech segments that may carry linguistic information, e.g. phones and syllables. In this study, we propose a novel phone-fortified perceptual loss (PFPL) that takes phonetic information into account for training SE models. To effectively incorporate the phonetic information, the PFPL is computed based on latent representations of the wav2vec model, a powerful self-supervised encoder that renders rich phonetic information. To more accurately measure the distribution distances of the latent representations, the PFPL adopts the Wasserstein distance as the distance measure. Our experimental results first reveal that the PFPL is more correlated with the perceptual evaluation metrics, as compared to signal-level losses. Moreover, the results showed that the PFPL can enable a deep complex U-Net SE model to achieve highly competitive performance in terms of standardized quality and intelligibility evaluations on the Voice Bank-DEMAND dataset.


title: Self-assembling kinetics: Accessing a new design space via differentiable statistical-physics models

abstract: The inverse problem of designing component interactions to target emergent structure is fundamental to numerous applications in biotechnology, materials science, and statistical physics. Equally important is the inverse problem of designing emergent kinetics, but this has received considerably less attention. Using recent advances in automatic differentiation, we show how kinetic pathways can be precisely designed by directly differentiating through statistical-physics models, namely free energy calculations and molecular dynamics simulations. We consider two systems that are crucial to our understanding of structural self-assembly: bulk crystallization and small nanoclusters. In each case we are able to assemble precise dynamical features. Using gradient information, we manipulate interactions among constituent particles to tune the rate at which these systems yield specific structures of interest. Moreover, we use this approach to learn non-trivial features about the high-dimensional design space, allowing us to accurately predict when multiple kinetic features can be simultaneously and independently controlled. These results provide a concrete and generalizable foundation for studying non-structural self-assembly, including kinetic properties as well as other complex emergent properties, in a vast array of systems.


title: Probing jet launching in neutron star X-ray binaries: the variable and polarized jet of SAX J1808.4-3658

abstract: We report on an optical photometric and polarimetric campaign on the accreting millisecond X-ray pulsar (AMXP) SAX J1808.4-3658 during its 2019 outburst. The emergence of a low-frequency excess in the spectral energy distribution in the form of a red excess above the disc spectrum (seen most prominently in z, i and R-bands) is observed as the outburst evolves. This is indicative of optically thin synchrotron emission due to a jet, as seen previously in this source and in other AMXPs during outburst. At the end of the outburst decay, the source entered a reflaring state. The low-frequency excess is still observed during the reflares. Our optical (BVRI) polarimetric campaign shows variable linear polarization (LP) throughout the outburst. We show that this is intrinsic to the source, with low-level but significant detections (0.2-2%) in all bands. The LP spectrum is red during both the main outburst and the reflaring state, favoring a jet origin for this variable polarization over other interpretations, such as Thomson scattering with free electrons from the disc or the propelled matter. During the reflaring state, a few episodes with stronger LP level (1-2 %) are observed. The low-level, variable LP is suggestive of strongly tangled magnetic fields near the base of the jet. These results clearly demonstrate how polarimetry is a powerful tool for probing the magnetic field structure in X-ray binary jets, similar to AGN jets.


title: Transferability of machine learning potentials: Protonated water neural network potential applied to the protonated water hexamer

abstract: A previously published neural network potential for the description of protonated water clusters up to the protonated water tetramer, H$^+$(H$_2$O)$_4$, at essentially converged coupled cluster accuracy (J. Chem. Theory Comput. 16, 88 (2020)) is applied to the protonated water hexamer, H$^+$(H$_2$O)$_6$ -- a system that the neural network has never seen before. Although being in the extrapolation regime, it is shown that the potential not only allows for quantum simulations from ultra-low temperatures $\sim$ 1 K up to 300 K, but that it is able to describe the new system very accurately compared to explicit coupled cluster calculations. This transferability of the model is rationalized by the similarity of the atomic environments encountered for the larger cluster compared to the environments in the training set of the model. Compared to the interpolation regime the quality of the model is reduced by roughly one order of magnitude, but most of the difference to the coupled cluster reference comes from global shifts of the potential energy surface, while local energy fluctuations are well recovered. These results suggest that the application of neural network potentials in extrapolation regimes can provide useful results and might be more general than usually thought.


title: On the algebraic definition of total rotation in RSA

abstract: Total rotation is a quantity that has been used for years in RSA. However, its definition has no mathematical sense, since the Euler angles do not form a vector space, since angles cannot define a multiplication group. With this work I tried to give a mathematical definition of the total rotation connecting the Euler description of the rotations with the helical axis. The approximation for small angles was used to connect Euler's angles and helical angle. With this approximation Euler angles acquire the properties of a vector space and it is possible to justify the meaning of this parameter. Validation test showed that total rotation has an approximation error between 5\% and 7\% for angles in the range $\left[-\frac{\pi}{6}, \frac{\pi}{6} \right]$. Since usually RSA uses smaller angle ranges, the approximation is perfectly suitable for use in RSA.


title: Invariant Functions On Cluster Ensembles

abstract: We define the notion of an invariant function on a cluster ensemble with respect to an action of the cluster modular group on its associated function fields. We realize many examples of previously studied functions as elements of this type of invariant ring and give many new examples. We show that these invariants have geometric and number theoretic interpretations, and classify them for ensembles associated to affine Dynkin diagrams.


title: On the dynamics of thin layers of viscous flows inside another viscous fluid

abstract: In this work we will study the dynamics of a thin layer of a viscous fluid which is embedded in the interior of another viscous fluid. The resulting flow can be approximated by means of the solutions of a free boundary problem for the Stokes equation in which one of the unknowns is the shape of a curve which approximates the geometry of the thin layer of fluid. We also derive the equation yielding the thickness of this fluid. This model, that will be termed as the Geometric Free Boundary Problem, will be derived using matched asymptotic expansions. We will prove that the Geometric Free Boundary Problem is well posed and the solutions of the thickness equation are well defined (in particular they do not yield breaking of fluid layers) as long as the solutions of the Geometric Free Boundary Problem exist.


title: Ensemble sampler for infinite-dimensional inverse problems

abstract: We introduce a new Markov chain Monte Carlo (MCMC) sampler for infinite-dimensional inverse problems. Our new sampler is based on the affine invariant ensemble sampler, which uses interacting walkers to adapt to the covariance structure of the target distribution. We extend this ensemble sampler for the first time to infinite-dimensional function spaces, yielding a highly efficient gradient-free MCMC algorithm. Because our new ensemble sampler does not require gradients or posterior covariance estimates, it is simple to implement and broadly applicable.


title: Latent Fibrations: Fibrations for Categories of Partial Maps

abstract: Latent fibrations are an adaptation, appropriate for categories of partial maps (as presented by restriction categories), of the usual notion of fibration. The paper initiates the development of the basic theory of latent fibrations and explores some key examples. Latent fibrations cover a wide variety of examples, some of which are partial versions of standard fibrations, and some of which are particular to partial map categories (particularly those that arise in computational settings). Latent fibrations with various special properties are identified: hyperconnected latent fibrations, in particular, are shown to support the construction of a fibrational dual; this is important to reverse differential programming and, more generally, in the theory of lenses.


title: Rise of the Titans: Gas Excitation and Feedback in a Binary Hyper-Luminous Dusty Starburst Galaxy at z~6

abstract: We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^11 (alpha_CO/1.0) Msun. We also find that the excitation of the star-forming gas is overall moderate for a z>5 dusty starburst, which is consistent with its moderate dust temperature. A high density, high kinetic temperature gas component embedded in the gas reservoir is required to fully explain the CO line ladder. This component is likely associated with the "maximum starburst" nuclei in the two merging galaxies, which are separated by only (140+/-13) km/s along the line of sight and 9.0 kpc in projection. The kinematic structure of both components is consistent with galaxy disks, but this interpretation remains limited by the spatial resolution of the current data. The OH+ features are only detected towards the northern component, which is also the one that is more enshrouded in dust and thus remains undetected up to 1.6 um even in our sensitive new HST/WFC3 imaging. The absorption component of the OH+ line is blueshifted and peaks near the CO and continuum emission peak while the emission is redshifted and peaks offset by 1.7 kpc from the CO and continuum emission peak, suggesting that the gas is associated with a massive molecular outflow from the intensely star-forming nucleus that supplies 125 Msun/yr of enriched gas to its halo.


title: Nanoscale plasticity and neuromorphic dynamics in silicon suboxide RRAM

abstract: Resistive random-access memories, also known as memristors, whose resistance can be modulated by the electrically driven formation and disruption of conductive filaments within an insulator, are promising candidates for neuromorphic applications due to their scalability, low-power operation and diverse functional behaviours. However, understanding the dynamics of individual filaments, and the surrounding material, is challenging, owing to the typically very large cross-sectional areas of test devices relative to the nanometre scale of individual filaments. In the present work, conductive atomic force microscopy is used to study the evolution of conductivity at the nanoscale in a fully CMOS-compatible silicon suboxide thin film. Distinct filamentary plasticity and background conductivity enhancement are reported, suggesting that device behaviour might be best described by composite core (filament) and shell (background conductivity) dynamics. Furthermore, constant current measurements demonstrate an interplay between filament formation and rupture, resulting in current-controlled voltage spiking in nanoscale regions, with an estimated optimal energy consumption of 25 attojoules per spike. This is very promising for extremely low-power neuromorphic computation and suggests that the dynamic behaviour observed in larger devices should persist and improve as dimensions are scaled down.


title: Novel imaging technique for $\alpha$-particles using a fast optical camera

abstract: A new imaging technique for $\alpha$-particles using a fast optical camera focused on a thin scintillator is presented. As $\alpha$-particles interact in a thin layer of LYSO fast scintillator, they produce a localized flash of light. The light is collected with a lens to an intensified optical camera, Tpx3Cam, with single photon sensitivity and excellent spatial & temporal resolutions. The interactions of photons with the camera is reconstructed by means of a custom algorithm, capable of discriminating single photons using time and spatial information.


title: Quasinormal modes and their overtones at the common horizon in a binary black hole merger

abstract: It is expected that all astrophysical black holes in equilibrium are well described by the Kerr solution. Moreover, any black hole far away from equilibrium, such as one initially formed in a compact binary merger or by the collapse of a massive star, will eventually reach a final equilibrium Kerr state. At sufficiently late times in this process of reaching equilibrium, we expect that the black hole is modeled as a perturbation around the final state. The emitted gravitational waves will then be damped sinusoids with frequencies and damping times given by the quasinormal mode spectrum of the final Kerr black hole. An observational test of this scenario, often referred to as black hole spectroscopy, is one of the major goals of gravitational wave astronomy. It was recently suggested that the quasinormal mode description including the higher overtones might hold even right after the remnant black hole is first formed. At these times, the black hole is expected to be highly dynamical and nonlinear effects are likely to be important. In this paper we investigate this remarkable scenario in terms of the horizon dynamics. Working with high accuracy simulations of a simple configuration, namely the head-on collision of two nonspinning black holes with unequal masses, we study the dynamics of the final common horizon in terms of its shear and its multipole moments. We show that they are indeed well described by a superposition of ringdown modes as long as a sufficiently large number of higher overtones are included. This description holds even for the highly dynamical final black hole shortly after its formation. We discuss the implications and caveats of this result for black hole spectroscopy and for our understanding of the approach to equilibrium.


title: A Study on Efficiency in Continual Learning Inspired by Human Learning

abstract: Humans are efficient continual learning systems; we continually learn new skills from birth with finite cells and resources. Our learning is highly optimized both in terms of capacity and time while not suffering from catastrophic forgetting. In this work we study the efficiency of continual learning systems, taking inspiration from human learning. In particular, inspired by the mechanisms of sleep, we evaluate popular pruning-based continual learning algorithms, using PackNet as a case study. First, we identify that weight freezing, which is used in continual learning without biological justification, can result in over $2\times$ as many weights being used for a given level of performance. Secondly, we note the similarity in human day and night time behaviors to the training and pruning phases respectively of PackNet. We study a setting where the pruning phase is given a time budget, and identify connections between iterative pruning and multiple sleep cycles in humans. We show there exists an optimal choice of iteration v.s. epochs given different tasks.


title: The statistical reliability of 267 GHz JCMT observations of Venus: No significant evidence for phosphine absorption

abstract: In the light of the recent announcement of the discovery of the potential biosignature phosphine in the atmosphere of Venus I present an independent reanalysis of the original JCMT data to assess the statistical reliability of the detection. Two line detection methods are explored, low order polynomial fits and higher order multiple polynomial fits. A non-parametric bootstrap analysis reveals that neither line detection method is able to recover a statistically significant detection. Similar to the results of other reanalyses of ALMA Venus spectra, the polynomial fitting process results in false positive detections in the JCMT spectrum. There is thus no significant evidence for phosphine absorption in the JCMT Venus spectra.


title: Wave-particle energy transfer directly observed in an ion cyclotron wave

abstract: Context. The first studies with Parker Solar Probe (PSP) data have made significant progress toward the understanding of the fundamental properties of ion cyclotron waves in the inner heliosphere. The survey mode particle measurements of PSP, however, did not make it possible to measure the coupling between electromagnetic fields and particles on the time scale of the wave periods. Aims. We present a novel approach to study wave-particle energy exchange with PSP. Methods. We use the Flux Angle operation mode of the Solar Probe Cup in conjunction with the electric field measurements and present a case study when the Flux Angle mode measured the direct interaction of the proton velocity distribution with an ion cyclotron wave. Results. Our results suggest that the energy transfer from fields to particles on the timescale of a cyclotron period is equal to approximately 3-6% of the electromagnetic energy flux. This rate is consistent with the hypothesis that the ion cyclotron wave was locally generated in the solar wind.


title: Testing gravity with the Milky Way: Yukawa potential

abstract: We test a Yukawa correction to the Newtonian potential, making use of our own Galaxy - the Milky Way - as a testbed. We include as free parameter the Yukawa strength and range and the dark matter profile parameters, and compare several morphologies for the bulge, gas, and disk components, also using Bayesian model selection criteria. We employ up-to-date datasets for both the visible (baryonic) component of the Milky Way, and for the tracers of the gravitational potential (the Rotation Curve). We find that the data are consistent with the Newtonian potential, and constrain the Yukawa coupling $\beta$ to be negative and $\lambda$ to range along the curve $\lambda = a|\beta|^{c}$ with $a = (0.77 \pm 0.06)$ kpc and $c = -0.503\substack{+0.016 \\ -0.019}$.


title: Chronic, cortex-wide imaging of specific cell populations during behavior

abstract: Measurements of neuronal activity across brain areas are important for understanding the neural correlates of cognitive and motor processes like attention, decision-making, and action selection. However, techniques that allow cellular resolution measurements are expensive and require a high degree of technical expertise, which limits their broad use. Widefield imaging of genetically encoded indicators is a high throughput, cost effective, and flexible approach to measure activity of specific cell populations with high temporal resolution and a cortex-wide field of view. Here we outline our protocol for assembling a widefield setup, a surgical preparation to image through the intact skull, and imaging neural activity chronically in behaving, transgenic mice that express a calcium indicator in specific subpopulations of cortical neurons. Further, we highlight a processing pipeline that leverages novel, cloud-based methods to analyze large-scale imaging datasets. The protocol targets labs that are seeking to build macroscopes, optimize surgical procedures for long-term chronic imaging, and/or analyze cortex-wide neuronal recordings.


title: Exact Self-Consistent Effective Hamiltonian Theory

abstract: We propose a general variational fermionic many-body wavefunction that generates an effective Hamiltonian in a quadratic form, which can then be exactly solved. The theory can be constructed within the density functional theory framework, and a self-consistent scheme is proposed for solving the exact density functional theory. We apply the theory to structurally-disordered systems, symmetric and asymmetric Hubbard dimers, and the corresponding lattice models. The single fermion excitation spectra show a persistent gap due to the fermionic-entanglement-induced pairing condensate. For disordered systems, the density of states at the edge of the gap diverges in the thermodynamic limit, suggesting a topologically ordered phase. A sharp resonance is predicted as the gap is not dependent on the temperature of the system. For the symmetric Hubbard model, the gap for both half-filling and doped case suggests that the quantum phase transition between the antiferromagnetic and superconducting phases is continuous.


title: Explicit stabilized multirate method for stiff stochastic differential equations

abstract: Stabilized explicit methods are particularly efficient for large systems of stiff stochastic differential equations (SDEs) due to their extended stability domain. However, they loose their efficiency when a severe stiffness is induced by very few "fast" degrees of freedom, as the stiff and nonstiff terms are evaluated concurrently. Therefore, inspired by [A. Abdulle, M. J. Grote, and G. Rosilho de Souza, Preprint (2020), arXiv:2006.00744] we introduce a stochastic modified equation whose stiffness depends solely on the "slow" terms. By integrating this modified equation with a stabilized explicit scheme we devise a multirate method which overcomes the bottleneck caused by a few severely stiff terms and recovers the efficiency of stabilized schemes for large systems of nonlinear SDEs. The scheme is not based on any scale separation assumption of the SDE and therefore it is employable for problems stemming from the spatial discretization of stochastic parabolic partial differential equations on locally refined grids. The multirate scheme has strong order 1/2, weak order 1 and its stability is proved on a model problem. Numerical experiments confirm the efficiency and accuracy of the scheme.


title: Additive processes on the unit circle and Loewner chains

abstract: This paper defines the notion of generators for a class of decreasing radial Loewner chains which are only continuous with respect to time. For this purpose, "Loewner's integral equation" which generalizes Loewner's differential equation is defined and analyzed. The definition of generators is motivated by the L\'evy-Khintchine representation for additive processes on the unit circle. Actually, we can and do introduce a homeomorphism between the above class of Loewner chains and the set of the distributions of increments of additive processes equipped with suitable topologies. On the other hand, from the viewpoint of non-commutative probability theory, the above generators also induce bijections with some other objects: in particular, monotone convolution hemigroups and free convolution hemigroups. Finally, the generators of Loewner chains constructed from free convolution hemigroups via subordination are computed.


title: Reinforcement Learning for Sparse-Reward Object-Interaction Tasks in a First-person Simulated 3D Environment

abstract: First-person object-interaction tasks in high-fidelity, 3D, simulated environments such as the AI2Thor virtual home-environment pose significant sample-efficiency challenges for reinforcement learning (RL) agents learning from sparse task rewards. To alleviate these challenges, prior work has provided extensive supervision via a combination of reward-shaping, ground-truth object-information, and expert demonstrations. In this work, we show that one can learn object-interaction tasks from scratch without supervision by learning an attentive object-model as an auxiliary task during task learning with an object-centric relational RL agent. Our key insight is that learning an object-model that incorporates object-attention into forward prediction provides a dense learning signal for unsupervised representation learning of both objects and their relationships. This, in turn, enables faster policy learning for an object-centric relational RL agent. We demonstrate our agent by introducing a set of challenging object-interaction tasks in the AI2Thor environment where learning with our attentive object-model is key to strong performance. Specifically, we compare our agent and relational RL agents with alternative auxiliary tasks to a relational RL agent equipped with ground-truth object-information, and show that learning with our object-model best closes the performance gap in terms of both learning speed and maximum success rate. Additionally, we find that incorporating object-attention into an object-model's forward predictions is key to learning representations which capture object-category and object-state.


title: A fast and scalable computational framework for large-scale and high-dimensional Bayesian optimal experimental design

abstract: We develop a fast and scalable computational framework to solve large-scale and high-dimensional Bayesian optimal experimental design problems. In particular, we consider the problem of optimal observation sensor placement for Bayesian inference of high-dimensional parameters governed by partial differential equations (PDEs), which is formulated as an optimization problem that seeks to maximize an expected information gain (EIG). Such optimization problems are particularly challenging due to the curse of dimensionality for high-dimensional parameters and the expensive solution of large-scale PDEs. To address these challenges, we exploit two essential properties of such problems: the low-rank structure of the Jacobian of the parameter-to-observable map to extract the intrinsically low-dimensional data-informed subspace, and the high correlation of the approximate EIGs by a series of approximations to reduce the number of PDE solves. We propose an efficient offline-online decomposition for the optimization problem: an offline stage of computing all the quantities that require a limited number of PDE solves independent of parameter and data dimensions, and an online stage of optimizing sensor placement that does not require any PDE solve. For the online optimization, we propose a swapping greedy algorithm that first construct an initial set of sensors using leverage scores and then swap the chosen sensors with other candidates until certain convergence criteria are met. We demonstrate the efficiency and scalability of the proposed computational framework by a linear inverse problem of inferring the initial condition for an advection-diffusion equation, and a nonlinear inverse problem of inferring the diffusion coefficient of a log-normal diffusion equation, with both the parameter and data dimensions ranging from a few tens to a few thousands.


title: Hopf actions of some quantum groups on path algebras

abstract: Our first collection of results parametrize (filtered) actions of a quantum Borel $U_q(\mathfrak{b}) \subset U_q(\mathfrak{sl}_2)$ on the path algebra of an arbitrary (finite) quiver. When $q$ is a root of unity, we give necessary and sufficient conditions for these actions to factor through corresponding finite-dimensional quotients, generalized Taft algebras $T(r,n)$ and small quantum groups $U_q(\mathfrak{sl}_2)$. In the second part of the paper, we shift to the language of tensor categories. Here we consider a quiver path algebra equipped with an action of a Hopf algebra $H$ to be a tensor algebra in the tensor category of representations $H$. Such a tensor algebra is generated by an algebra and bimodule in this tensor category. Our second collection of results describe the corresponding bimodule categories via an equivalence with categories of representations of certain explicitly described quivers with relations.


title: Room Temperature Quantum Coherent Revival in an Ensemble of Artificial Atoms

abstract: We report a demonstration of the hallmark concept of quantum optics: periodic collapse and revival of quantum coherence (QCR) in a room temperature ensemble of quantum dots (QD). Control over quantum states, inherent to QCR, together with the dynamical QD properties present an opportunity for practical room temperature building blocks of quantum information processing. The amplitude decay of QCR is dictated by the QD homogeneous linewidth, thus, enabling its extraction in a double-pulse Ramsey-type experiment. The more common photon echo technique was also invoked and yielded the same linewidth. Measured electrical bias and temperature dependencies of the transverse relaxation times enable to determine the two main decoherence mechanisms: carrier-carrier and carrier-phonon scatterings.


title: Orientation-dependent stabilization of MgCr$_2$O$_4$ spinel thin films

abstract: AB$_2$O$_4$ normal spinels with a magnetic B site can host a variety of magnetic and orbital frustrations leading to spin-liquid phases and field-induced phase transitions. Here we report the first epitaxial growth of (111)-oriented MgCr$_2$O$_4$ thin films. By characterizing the structural and electronic properties of films grown along (001) and (111) directions, the influence of growth orientation has been studied. Despite distinctly different growth modes observed during deposition, the comprehensive characterization reveals no measurable disorder in the cation distribution nor multivalency issue for Cr ions in either orientation. Contrary to a naive expectation, the (111) stabilized films exhibit a smoother surface and a higher degree of crystallinity than (001)-oriented films. The preference in growth orientation is explained within the framework of heteroepitaxial stabilization in connection to a significantly lower (111) surface energy. These findings open broad opportunities in the fabrication of 2D kagome-triangular heterostructures with emergent magnetic behavior inaccessible in bulk crystals.


title: Extending Gravitational Wave Extraction Using Weyl Characteristic Fields

abstract: We present a detailed methodology for extracting the full set of Newman-Penrose Weyl scalars from numerically generated spacetimes without requiring a tetrad that is completely orthonormal or perfectly aligned to the principal null directions. We also describe how to implement an extrapolation technique for computing the Weyl scalars' contribution at asymptotic null infinity in postprocessing. These methods have been used to produce $\Psi_4$ and $h$ waveforms for the Simulating eXtreme Spacetimes (SXS) waveform catalog and now have been expanded to produce the entire set of Weyl scalars. These new waveform quantities are critical for the future of gravitational wave astronomy in order to understand the finite-amplitude gauge differences that can occur in numerical waveforms. We also present a new analysis of the accuracy of waveforms produced by the Spectral Einstein Code. While ultimately we expect Cauchy characteristic extraction to yield more accurate waveforms, the extraction techniques described here are far easier to implement and have already proven to be a viable way to produce production-level waveforms that can meet the demands of current gravitational-wave detectors.


title: Forecasting Hamiltonian dynamics without canonical coordinates

abstract: Conventional neural networks are universal function approximators, but because they are unaware of underlying symmetries or physical laws, they may need impractically many training data to approximate nonlinear dynamics. Recently introduced Hamiltonian neural networks can efficiently learn and forecast dynamical systems that conserve energy, but they require special inputs called canonical coordinates, which may be hard to infer from data. Here we significantly expand the scope of such networks by demonstrating a simple way to train them with any set of generalised coordinates, including easily observable ones.


title: Gwbench: a novel Fisher information package for gravitational-wave benchmarking

abstract: We present a new Python package, gwbench, implementing the well-established Fisher information formalism as a fast and straightforward tool for the purpose of gravitational-wave benchmarking, i.e. the estimation of signal-to-noise ratios and measurement errors of gravitational waves observed by a network of detectors. Such an infrastructure is necessary due to the high computational cost of Bayesian parameter estimation methods which renders them less effective for the scientific assessment of gravitational waveforms, detectors, and networks of detectors, especially when determining their effects on large populations of gravitational-wave sources spread throughout the universe. gwbench further gives quick access to detector locations and sensitivities, while including the effects of Earth's rotation on the latter, as well as waveform models and their derivatives, while giving access to the host of waveforms available in the LSC Algorithm Library. With the provided functionality, gwbench is relevant for a wide variety of applications in gravitational-wave astronomy such as waveform modeling, detector development, cosmology, and tests of general relativity.


title: Micromobility in Smart Cities: A Closer Look at Shared Dockless E-Scooters via Big Social Data

abstract: The micromobility is shaping first- and last-mile travels in urban areas. Recently, shared dockless electric scooters (e-scooters) have emerged as a daily alternative to driving for short-distance commuters in large cities due to the affordability, easy accessibility via an app, and zero emissions. Meanwhile, e-scooters come with challenges in city management, such as traffic rules, public safety, parking regulations, and liability issues. In this paper, we collected and investigated 5.8 million scooter-tagged tweets and 144,197 images, generated by 2.7 million users from October 2018 to March 2020, to take a closer look at shared e-scooters via crowdsourcing data analytics. We profiled e-scooter usages from spatial-temporal perspectives, explored different business roles (i.e., riders, gig workers, and ridesharing companies), examined operation patterns (e.g., injury types, and parking behaviors), and conducted sentiment analysis. To our best knowledge, this paper is the first large-scale systematic study on shared e-scooters using big social data.


title: Shortest closed curve to inspect a sphere

abstract: We show that in Euclidean 3-space any closed curve which lies outside the unit sphere and contains the sphere within its convex hull has length at least $4\pi$. Equality holds only when the curve is composed of $4$ semicircles of length $\pi$, arranged in the shape of a baseball seam, as conjectured by V. A. Zalgaller in 1996.


title: WASP-186 and WASP-187: two hot Jupiters discovered by SuperWASP and SOPHIE with additional observations by TESS

abstract: We present the discovery of two new hot Jupiters identified from the WASP survey, WASP-186b and WASP-187b (TOI-1494.01 and TOI-1493.01). Their planetary nature was established from SOPHIE spectroscopic observations, and additional photometry was obtained from TESS. Stellar parameters for the host stars are derived from spectral line, IRFM, and isochrone placement analyses. These parameters are combined with the photometric and radial velocity data in an MCMC method to determine the planetary properties. WASP-186b is a massive Jupiter (4.22 +/- 0.18 M_J, 1.11 +/-0.03 R_J) orbiting a mid-F star on a 5.03 day eccentric (e=0.327 +/- 0.008) orbit. WASP-187b is a low density (0.80 +/- 0.09 M_J, 1.64 +/- 0.05 R_J) planet in a 5.15 day circular orbit around a slightly evolved early F-type star.


title: Rosella: A Self-Driving Distributed Scheduler for Heterogeneous Clusters

abstract: Large-scale interactive web services and advanced AI applications make sophisticated decisions in real-time, based on executing a massive amount of computation tasks on thousands of servers. Task schedulers, which often operate in heterogeneous and volatile environments, require high throughput, i.e., scheduling millions of tasks per second, and low latency, i.e., incurring minimal scheduling delays for millisecond-level tasks. Scheduling is further complicated by other users' workloads in a shared system, other background activities, and the diverse hardware configurations inside datacenters. We present Rosella, a new self-driving, distributed approach for task scheduling in heterogeneous clusters. Rosella automatically learns the compute environment and adjusts its scheduling policy in real-time. The solution provides high throughput and low latency simultaneously because it runs in parallel on multiple machines with minimum coordination and only performs simple operations for each scheduling decision. Our learning module monitors total system load and uses the information to dynamically determine optimal estimation strategy for the backends' compute-power. Rosella generalizes power-of-two-choice algorithms to handle heterogeneous workers, reducing the max queue length of O(log n) obtained by prior algorithms to O(log log n). We evaluate Rosella with a variety of workloads on a 32-node AWS cluster. Experimental results show that Rosella significantly reduces task response time, and adapts to environment changes quickly.


title: Space-Time Covid-19 Bayesian SIR modeling in South Carolina

abstract: The Covid-19 pandemic has spread across the world since the beginning of 2020. Many regions have experienced its effects. The state of South Carolina in the USA has seen cases since early March 2020 and a primary peak in early April 2020. A lockdown was imposed on April 6th but lifting of restrictions started on April 24th. The daily case and death data as reported by NCHS (deaths) via the New York Times GitHUB repository have been analyzed and approaches to modeling of the data are presented. Prediction is also considered and the role of asymptomatic transmission is assessed as a latent unobserved effect. Two different time periods are examined and one step prediction is provided.


title: Identifying Entangled Physics Relationships through Sparse Matrix Decomposition to Inform Plasma Fusion Design

abstract: A sustainable burn platform through inertial confinement fusion (ICF) has been an ongoing challenge for over 50 years. Mitigating engineering limitations and improving the current design involves an understanding of the complex coupling of physical processes. While sophisticated simulations codes are used to model ICF implosions, these tools contain necessary numerical approximation but miss physical processes that limit predictive capability. Identification of relationships between controllable design inputs to ICF experiments and measurable outcomes (e.g. yield, shape) from performed experiments can help guide the future design of experiments and development of simulation codes, to potentially improve the accuracy of the computational models used to simulate ICF experiments. We use sparse matrix decomposition methods to identify clusters of a few related design variables. Sparse principal component analysis (SPCA) identifies groupings that are related to the physical origin of the variables (laser, hohlraum, and capsule). A variable importance analysis finds that in addition to variables highly correlated with neutron yield such as picket power and laser energy, variables that represent a dramatic change of the ICF design such as number of pulse steps are also very important. The obtained sparse components are then used to train a random forest (RF) surrogate for predicting total yield. The RF performance on the training and testing data compares with the performance of the RF surrogate trained using all design variables considered. This work is intended to inform design changes in future ICF experiments by augmenting the expert intuition and simulations results.


title: Ground Roll Suppression using Convolutional Neural Networks

abstract: Seismic data processing plays a major role in seismic exploration as it conditions much of the seismic interpretation performance. In this context, generating reliable post-stack seismic data depends also on disposing of an efficient pre-stack noise attenuation tool. Here we tackle ground roll noise, one of the most challenging and common noises observed in pre-stack seismic data. Since ground roll is characterized by relative low frequencies and high amplitudes, most commonly used approaches for its suppression are based on frequency-amplitude filters for ground roll characteristic bands. However, when signal and noise share the same frequency ranges, these methods usually deliver also signal suppression or residual noise. In this paper we take advantage of the highly non-linear features of convolutional neural networks, and propose to use different architectures to detect ground roll in shot gathers and ultimately to suppress them using conditional generative adversarial networks. Additionally, we propose metrics to evaluate ground roll suppression, and report strong results compared to expert filtering. Finally, we discuss generalization of trained models for similar and different geologies to better understand the feasibility of our proposal in real applications.


title: On Linearizability and the Termination of Randomized Algorithms

abstract: We study the question of whether the "termination with probability 1" property of a randomized algorithm is preserved when one replaces the atomic registers that the algorithm uses with linearizable (implementations of) registers. We show that in general this is not so: roughly speaking, every randomized algorithm A has a corresponding algorithm A' that solves the same problem if the registers that it uses are atomic or strongly-linearizable, but does not terminate if these registers are replaced with "merely" linearizable ones. Together with a previous result shown in [15], this implies that one cannot use the well-known ABD implementation of registers in message-passing systems to automatically transform any randomized algorithm that works in shared-memory systems into a randomized algorithm that works in message-passing systems: with a strong adversary the resulting algorithm may not terminate.


title: Safety-Aware Cascade Controller Tuning Using Constrained Bayesian Optimization

abstract: This paper presents an automated, model-free, data-driven method for the safe tuning of PID cascade controller gains based on Bayesian optimization. The optimization objective is composed of data-driven performance metrics and modeled using Gaussian processes. We further introduce a data-driven constraint that captures the stability requirements from system data. Numerical evaluation shows that the proposed approach outperforms relay feedback autotuning and quickly converges to the global optimum, thanks to a tailored stopping criterion. We demonstrate the performance of the method in simulations and experiments. For experimental implementation, in addition to the introduced safety constraint, we integrate a method for automatic detection of the critical gains and extend the optimization objective with a penalty depending on the proximity of the current candidate points to the critical gains. The resulting automated tuning method optimizes system performance while ensuring stability and standardization


title: Credit Valuation Adjustment in Credit Risk with Simultaneous Defaults Possibility

abstract: In a series of recent papers, Damiano Brigo, Andrea Pallavicini, and co-authors have shown that the value of a contract in a Credit Valuation Adjustment (CVA) setting, being the sum of the cash flows, can be represented as a solution of a decoupled forward-backward stochastic differential equation (FBSDE). CVA is the difference between the risk-free portfolio value and the true portfolio value that takes into account the possibility of a counter party's default. In other words, CVA is the market value of counter party credit risk. This has achieved noteworthy importance after the 2008 financial debacle, where counter party risk played an under-modeled but huge risk. In their analysis, Brigo et al make the classical assumption of conditional independence of the default times, given the risk free market filtration. This does not allow for the increasingly likely case of simultaneous defaults. We weaken their assumption, replacing it with a martingale orthogonality condition. This in turn changes the form of the BSDE that arises from the model.


title: SDSS-IV MaNGA: When is morphology imprinted on galaxies?

abstract: It remains an open question as to how long ago the morphology that we see in a present-day galaxy was typically imprinted. Studies of galaxy populations at different redshifts reveal that the balance of morphologies has changed over time, but such snapshots cannot uncover the typical timescales over which individual galaxies undergo morphological transformation, nor which are the progenitors of today's galaxies of different types. However, these studies also show a strong link between morphology and star-formation rate over a large range in redshift, which offers an alternative probe of morphological transformation. We therefore derive the evolution in star-formation rate and stellar mass of a sample of 4342 galaxies in the SDSS-IV MaNGA survey through a stellar population "fossil record" approach, and show that the average evolution of the population shows good agreement with known behaviour from previous studies. Although the correlation between a galaxy's contemporaneous morphology and star-formation rate is strong over a large range of lookback times, we find that a galaxy's present-day morphology only correlates with its relatively recent (~2 Gyr) star-formation history. We therefore find strong evidence that morphological transitions to galaxies' current appearance occurred on timescales as short as a few billion years.


title: Inference of ventricular activation properties from non-invasive electrocardiography

abstract: The realisation of precision cardiology requires novel techniques for the non-invasive characterisation of individual patients' cardiac function to inform therapeutic and diagnostic decision-making. The electrocardiogram (ECG) is the most widely used clinical tool for cardiac diagnosis. Its interpretation is, however, confounded by functional and anatomical variability in heart and torso. In this study, we develop new computational techniques to estimate key ventricular activation properties for individual subjects by exploiting the synergy between non-invasive electrocardiography and image-based torso-biventricular modelling and simulation. More precisely, we present an efficient sequential Monte Carlo approximate Bayesian computation-based inference method, integrated with Eikonal simulations and torso-biventricular models constructed based on clinical cardiac magnetic resonance (CMR) imaging. The method also includes a novel strategy to treat combined continuous (conduction speeds) and discrete (earliest activation sites) parameter spaces, and an efficient dynamic time warping-based ECG comparison algorithm. We demonstrate results from our inference method on a cohort of twenty virtual subjects with cardiac volumes ranging from 74 cm3 to 171 cm3 and considering low versus high resolution for the endocardial discretisation (which determines possible locations of the earliest activation sites). Results show that our method can successfully infer the ventricular activation properties from non-invasive data, with higher accuracy for earliest activation sites, endocardial speed, and sheet (transmural) speed in sinus rhythm, rather than the fibre or sheet-normal speeds.


title: Decimation and Interleaving Operations in One-Sided Symbolic Dynamics

abstract: This paper studies subsets of one-sided shift spaces on a finite alphabet. Such subsets arise in symbolic dynamics, in fractal constructions, and in number theory. We study a family of decimation operations, which extract subsequences of symbol sequences in infinite arithmetic progressions, and show they are closed under composition. We also study a family of $n$-ary interleaving operations, one for each $n \ge 1$. Given subsets $X_0, X_1, ..., X_{n-1}$ of the shift space, the $n$-ary interleaving operator produces a set whose elements combine individual elements ${\bf x}_i$, one from each $X_i$, by interleaving their symbol sequences cyclically in arithmetic progressions $(\bmod\,n)$. We determine algebraic relations between decimation and interleaving operators and the shift operator. We study set-theoretic $n$-fold closure operations $X \mapsto X^{[n]}$, which interleave decimations of $X$ of modulus level $n$. A set is $n$-factorizable if $X=X^{[n]}$. The $n$-fold interleaving operators are closed under composition and are idempotent. To each $X$ we assign the set $\mathcal{N}(X)$ of all values $n \ge 1$ for which $X= X^{[n]}$. We characterize the possible sets $\mathcal{N}(X)$ as nonempty sets of positive integers that form a distributive lattice under the divisibility partial order and are downward closed under divisibility. We show that all sets of this type occur. We introduce a class of weakly shift-stable sets and show that this class is closed under all decimation, interleaving, and shift operations. This class includes all shift-invariant sets. We study two notions of entropy for subsets of the full one-sided shift and show that they coincide for weakly shift-stable $X$, but can be different in general. We give a formula for entropy of interleavings of weakly shift-stable sets in terms of individual entropies.


title: Fractional parts of non-integer powers of primes

abstract: Let $\alpha > 0$ be any fixed non-integer, $I$ be any subinterval of $[0; 1)$. In the paper, we prove an analogue of Bombieri-Vinogradov theorem for the set of primes $p$ satisfying the condition $\{ p^{\alpha} \} \in I$. This strengthens the previous result of Gritsenko and Zinchenko.


title: Away from Trolley Problems and Toward Risk Management

abstract: As automated vehicles receive more attention from the media, there has been an equivalent increase in the coverage of the ethical choices a vehicle may be forced to make in certain crash situations with no clear safe outcome. Much of this coverage has focused on a philosophical thought experiment known as the "trolley problem," and substituting an automated vehicle for the trolley and the car's software for the bystander. While this is a stark and straightforward example of ethical decision making for an automated vehicle, it risks marginalizing the entire field if it is to become the only ethical problem in the public's mind. In this chapter, I discuss the shortcomings of the trolley problem, and introduce more nuanced examples that involve crash risk and uncertainty. Risk management is introduced as an alternative approach, and its ethical dimensions are discussed.


title: StencilFlow: Mapping Large Stencil Programs to Distributed Spatial Computing Systems

abstract: Spatial computing devices have been shown to significantly accelerate stencil computations, but have so far relied on unrolling the iterative dimension of a single stencil operation to increase temporal locality. This work considers the general case of mapping directed acyclic graphs of heterogeneous stencil computations to spatial computing systems, assuming large input programs without an iterative component. StencilFlow maximizes temporal locality and ensures deadlock freedom in this setting, providing end-to-end analysis and mapping from a high-level program description to distributed hardware. We evaluate our generated architectures on a Stratix 10 FPGA testbed, yielding 1.31 TOp/s and 4.18 TOp/s on single-device and multi-device, respectively, demonstrating the highest performance recorded for stencil programs on FPGAs to date. We then leverage the framework to study a complex stencil program from a production weather simulation application. Our work enables productively targeting distributed spatial computing systems with large stencil programs, and offers insight into architecture characteristics required for their efficient execution in practice.


title: Analysis of high-order velocity moments in a strained channel flow

abstract: In the current study, model expressions for fifth-order velocity moments obtained from the truncated Gram-Charlier series expansions model for a turbulent flow field probability density function are validated using data from direct numerical simulation (DNS) of a planar turbulent flow in a strained channel. Simplicity of the model expressions, the lack of unknown coefficients, and their applicability to non-Gaussian turbulent flows make this approach attractive to use for closing turbulent models based on the Reynolds-averaged Navier-Stokes equations. The study confirms validity of the model expressions. It also shows that the imposed flow deformation improves an agreement between the model and DNS profiles for the fifth-order moments in the flow buffer zone including when the flow reverses its direction. The study reveals sensitivity of particularly odd velocity moments to the grid resolution. A new length scale is proposed as a criterion for the grid generation near the wall and in the other flow areas dominated by high mean velocity gradients when higher-order statistics have to be collected from DNS.


title: Symmetry properties of superconducting order parameter in Sr$_2$RuO$_4$

abstract: Soon after the discovery of superconductivity in Sr2RuO4 (SRO) in 1994, it was conjectured that its order parameter (OP) has a form similar to that realized in the superfluid He-3, namely, odd parity and spin triplet. The chiral p-wave pairing realized in He-3 A-phase was favored by several early experiments, in particular, the muon spin rotation and the Knight shift measurements published in 1998. However, the original Knight shift result was called into question in early 2019, raising the question as to whether the chiral p-wave, or even the spin-triplet pairing itself, is indeed realized in SRO. In this brief pedagogical review, we will address this question by counterposing the currently accepted results of Knight shift, polarized neutron scattering (PNS), spin counterflow half-quantum vortex (HQV), and Josephson experiments with predictions made both by standard BCS theory and by more general arguments based only on (1) the symmetry of the Hamiltonian including the spin-orbital terms, (2) thermodynamics, and (3) the qualitative experimental features of SRO. We also introduce a notation for triplet states alternative to the more popular d-vector one which we believe well suited to SRO. We conclude that the most recent Knight shift and PNS experiments do not exclude in the bulk the odd-parity, spin-triplet "helical" states allowed by the symmetry group of SRO but do exclude the chiral p-wave state. On the other hand, the Josephson and the HQV experiments showed that the pairing in SRO cannot be of the even-parity, spin-singlet type, and in the surface region or in samples of mesoscopic size, the d-vector must have a c-axis component, thus excluding apparently all bulk p-wave states except the chiral p-wave. Possible resolution of this rather glaring prima facie contradiction is discussed, taking into account other important experiments on SRO touched upon only briefly in this article.


title: Geometric Sampling of Networks

abstract: Motivated by the methods and results of manifold sampling based on Ricci curvature, we propose a similar approach for networks. To this end we make appeal to three types of discrete curvature, namely the graph Forman-, full Forman- and Haantjes-Ricci curvatures for edge-based and node-based sampling. We present the results of experiments on real life networks, as well as for square grids arising in Image Processing. Moreover, we consider fitting Ricci flows and we employ them for the detection of networks' backbone. We also develop embedding kernels related to the Forman-Ricci curvatures and employ them for the detection of the coarse structure of networks, as well as for network visualization with applications to SVM. The relation between the Ricci curvature of the original manifold and that of a Ricci curvature driven discretization is also studied.


title: Exploring complex networks with the ICON R package

abstract: We introduce ICON, an R package that contains 1075 complex network datasets in a standard edgelist format. All provided datasets have associated citations and have been indexed by the Colorado Index of Complex Networks - also referred to as ICON. In addition to supplying a large and diverse corpus of useful real-world networks, ICON also implements an S3 generic to work with the network and ggnetwork R packages for network analysis and visualization, respectively. Sample code in this report also demonstrates how ICON can be used in conjunction with the igraph package. Currently, the Comprehensive R Archive Network hosts ICON v0.4.0. We hope that ICON will serve as a standard corpus for complex network research and prevent redundant work that would be otherwise necessary by individual research groups. The open source code for ICON and for this reproducible report can be found at https://github.com/rrrlw/ICON.


title: Design Diversity for Improving Efficiency and Reducing Risk in Oil and Gas Well Stimulation under Uncertain Reservoir Conditions

abstract: Hydraulic fracturing stimulates fracture swarm in reservoir formation though pressurized injection fluid. However restricted by the availability of formation data, the variability embraced by reservoir keeps uncertain, driving unstable gas recovery along with low resource efficiency, being responsible for resource scarcity, contaminated water, and injection-induced earthquake. Resource efficiency is qualified though new determined energy efficiency, a scale of recovery and associated environmental footprint. To maximize energy efficiency while minimize its' variation, we issue picked designs at reservoir conditions dependent optimal probabilities, assembling high efficiency portfolios and low risk portfolios for portfolio combination, which balance the variation and efficiency at optimal by adjusting the proportion of each portfolio. Relative to regular design for one well, the optimal portfolio combination applied in multiple wells receive remarkable variation reduction meanwhile substantial energy efficiency increase, in response to the call of more recovery per unit investment and less environment cost per unit nature gas extracted.


title: Projection on particle number and angular momentum: Example of triaxial Bogoliubov quasiparticle states

abstract: Many quantal many-body methods that aim at the description of self-bound nuclear or mesoscopic electronic systems make use of auxiliary wave functions that break one or several of the symmetries of the Hamiltonian in order to include correlations associated with the geometrical arrangement of the system's constituents. Such reference states have been used already for a long time within self-consistent methods that are either based on effective valence-space Hamiltonians or energy density functionals, and they are presently also gaining popularity in the design of novel ab-initio methods. A fully quantal treatment of a self-bound many-body system, however, requires the restoration of the broken symmetries through the projection of the many-body wave functions of interest onto good quantum numbers. The goal of this work is three-fold. First, we want to give a general presentation of the formalism of the projection method starting from the underlying principles of group representation theory. Second, we want to investigate formal and practical aspects of the numerical implementation of particle-number and angular-momentum projection of Bogoliubov quasiparticle vacua, in particular with regard of obtaining accurate results at minimal computational cost. Third, we want to analyze the numerical, computational and physical consequences of intrinsic symmetries of the symmetry-breaking states when projecting them.


title: A Visuospatial Dataset for Naturalistic Verb Learning

abstract: We introduce a new dataset for training and evaluating grounded language models. Our data is collected within a virtual reality environment and is designed to emulate the quality of language data to which a pre-verbal child is likely to have access: That is, naturalistic, spontaneous speech paired with richly grounded visuospatial context. We use the collected data to compare several distributional semantics models for verb learning. We evaluate neural models based on 2D (pixel) features as well as feature-engineered models based on 3D (symbolic, spatial) features, and show that neither modeling approach achieves satisfactory performance. Our results are consistent with evidence from child language acquisition that emphasizes the difficulty of learning verbs from naive distributional data. We discuss avenues for future work on cognitively-inspired grounded language learning, and release our corpus with the intent of facilitating research on the topic.


title: Uncertainty Relations in Hydrodynamics

abstract: The uncertainty relations in hydrodynamics are numerically studied. We first give a review for the formulation of the generalized uncertainty relations in the stochastic variational method (SVM), following the paper by two of the present authors [Phys. Lett. A382, 1472 (2018)]. In this approach, the origin of the finite minimum value of uncertainty is attributed to the non-differentiable (virtual) trajectory of a quantum particle and then both of the Kennard and Robertson-Schr\"{o}dinger inequalities in quantum mechanics are reproduced. The same non-differentiable trajectory is applied to the motion of fluid elements in hydrodynamics. By introducing the standard deviations of position and momentum for fluid elements, the uncertainty relations in hydrodynamics are derived. These are applicable even to the Gross-Pitaevskii equation and then the field-theoretical uncertainty relation is reproduced. We further investigate numerically the derived relations and find that the behaviors of the uncertainty relations for liquid and gas are qualitatively different. This suggests that the uncertainty relations in hydrodynamics are used as a criterion to classify liquid and gas in fluid.


title: Simulation of detection and scattering of sound waves by the lateral line of a fish

abstract: An explicitly solvable model of sound sources detection and sound wave scattering by the lateral line of a fish or amphibian based on the theory of self-adjoint perturbations of the Laplace operator by arrays of additional boundary conditions at isolated points is proposed and discussed.


title: Frequency comb generation in a pulse-pumped normal dispersion Kerr mini-resonator

abstract: Kerr microresonators driven in the normal dispersion regime typically require the presence of localized dispersion perturbations, such as those induced by avoided mode crossings, to initiate the formation of optical frequency combs. In this work, we experimentally demonstrate that this requirement can be lifted by driving the resonator with a pulsed pump source. We also show that controlling the desynchronization between the pump repetition rate and the cavity free-spectral range (FSR) provides a simple mechanism to tune the center frequency of the output comb. Using a fiber mini-resonator with a radius of only 6 cm we experimentally present spectrally flat combs with a bandwidth of 3 THz whose center frequency can be tuned by more than 2 THz. By driving the cavity at harmonics of its 0.54 GHz FSR, we are able to generate combs with line spacings selectable between 0.54 and 10.8 GHz. The ability to tune both the center frequency and frequency spacing of the output comb highlights the flexibility of this platform. Additionally, we demonstrate that under conditions of large pump-cavity desynchronization, the same cavity also supports a new form of Raman-assisted anomalous dispersion cavity soliton.


title: Speech-Based Emotion Recognition using Neural Networks and Information Visualization

abstract: Emotions recognition is commonly employed for health assessment. However, the typical metric for evaluation in therapy is based on patient-doctor appraisal. This process can fall into the issue of subjectivity, while also requiring healthcare professionals to deal with copious amounts of information. Thus, machine learning algorithms can be a useful tool for the classification of emotions. While several models have been developed in this domain, there is a lack of userfriendly representations of the emotion classification systems for therapy. We propose a tool which enables users to take speech samples and identify a range of emotions (happy, sad, angry, surprised, neutral, clam, disgust, and fear) from audio elements through a machine learning model. The dashboard is designed based on local therapists' needs for intuitive representations of speech data in order to gain insights and informative analyses of their sessions with their patients.


title: Maximal green sequences for string algebras

abstract: Maximal green sequences are important objects in representation theory, cluster algebras, and string theory. The two fundamental questions about maximal green sequences are whether a given algebra admits such sequences and, if so, does it admit only finitely many. We study maximal green sequences in the case of string algebras and give sufficient conditions on the algebra that ensure an affirmative answer to these questions.


title: Measure, dimension, and complexity of the transient motion in Hamiltonian systems

abstract: Hamiltonian systems that are either open, leaking, or contain holes in the phase space possess solutions that eventually escape the system's domain. The motion described by such escape orbits before crossing the escape threshold can be understood as a transient behavior. In this work, we introduce a numerical method to visually illustrate and quantify the transient motion in Hamiltonian systems based on the transient measure, a finite-time version of the natural measure. We apply this method to two physical systems: the single-null divertor tokamak, described by a symplectic map; and the Earth-Moon system, as modeled by the planar circular restricted three-body problem. Our results portray how different locations for the ensemble of initial conditions may lead to different transient dynamical scenarios in both systems. We show that these scenarios can be properly quantified from a geometrical aspect, the transient correlation dimension, and a dynamical aspect, the transient complexity coefficient.


title: Construction Payment Automation Using Blockchain-Enabled Smart Contracts and Reality Capture Technologies

abstract: This paper presents a smart contract-based solution for autonomous administration of construction progress payments. It bridges the gap between payments (cash flow) and the progress assessments at job sites (product flow) enabled by reality capture technologies and building information modeling (BIM). The approach eliminates the reliance on the centralized and heavily intermediated mechanisms of existing payment applications. The construction progress is stored in a distributed manner using content addressable file sharing; it is broadcasted to a smart contract which automates the on-chain payment settlements and the transfer of lien rights. The method was successfully used for processing payments to 7 subcontractors in two commercial construction projects where progress monitoring was performed using a camera-equipped unmanned aerial vehicle (UAV) and an unmanned ground vehicle (UGV) equipped with a laser scanner. The results show promise for the method's potential for increasing the frequency, granularity, and transparency of payments. The paper is concluded with a discussion of implications for project management, introducing a new model of project as a singleton state machine.


title: Accurate Prostate Cancer Detection and Segmentation on Biparametric MRI using Non-local Mask R-CNN with Histopathological Ground Truth

abstract: Purpose: We aimed to develop deep machine learning (DL) models to improve the detection and segmentation of intraprostatic lesions (IL) on bp-MRI by using whole amount prostatectomy specimen-based delineations. We also aimed to investigate whether transfer learning and self-training would improve results with small amount labelled data. Methods: 158 patients had suspicious lesions delineated on MRI based on bp-MRI, 64 patients had ILs delineated on MRI based on whole mount prostatectomy specimen sections, 40 patients were unlabelled. A non-local Mask R-CNN was proposed to improve the segmentation accuracy. Transfer learning was investigated by fine-tuning a model trained using MRI-based delineations with prostatectomy-based delineations. Two label selection strategies were investigated in self-training. The performance of models was evaluated by 3D detection rate, dice similarity coefficient (DSC), 95 percentile Hausdrauff (95 HD, mm) and true positive ratio (TPR). Results: With prostatectomy-based delineations, the non-local Mask R-CNN with fine-tuning and self-training significantly improved all evaluation metrics. For the model with the highest detection rate and DSC, 80.5% (33/41) of lesions in all Gleason Grade Groups (GGG) were detected with DSC of 0.548[0.165], 95 HD of 5.72[3.17] and TPR of 0.613[0.193]. Among them, 94.7% (18/19) of lesions with GGG > 2 were detected with DSC of 0.604[0.135], 95 HD of 6.26[3.44] and TPR of 0.580[0.190]. Conclusion: DL models can achieve high prostate cancer detection and segmentation accuracy on bp-MRI based on annotations from histologic images. To further improve the performance, more data with annotations of both MRI and whole amount prostatectomy specimens are required.


title: Linear Regression Games: Convergence Guarantees to Approximate Out-of-Distribution Solutions

abstract: Recently, invariant risk minimization (IRM) (Arjovsky et al.) was proposed as a promising solution to address out-of-distribution (OOD) generalization. In Ahuja et al., it was shown that solving for the Nash equilibria of a new class of "ensemble-games" is equivalent to solving IRM. In this work, we extend the framework in Ahuja et al. for linear regressions by projecting the ensemble-game on an $\ell_{\infty}$ ball. We show that such projections help achieve non-trivial OOD guarantees despite not achieving perfect invariance. For linear models with confounders, we prove that Nash equilibria of these games are closer to the ideal OOD solutions than the standard empirical risk minimization (ERM) and we also provide learning algorithms that provably converge to these Nash Equilibria. Empirical comparisons of the proposed approach with the state-of-the-art show consistent gains in achieving OOD solutions in several settings involving anti-causal variables and confounders.


title: An environmental dependence of the physical and structural properties in the Hydra Cluster galaxies

abstract: The nearby Hydra Cluster ($\sim$50 Mpc) is an ideal laboratory to understand, in detail, the influence of the environment on the morphology and quenching of galaxies in dense environments. We study the Hydra cluster galaxies in the inner regions ($1R_{200}$) of the cluster using data from the Southern Photometric Local Universe Survey (S-PLUS), which uses 12 narrow and broad band filters in the visible region of the spectrum. We analyse structural (S\'ersic index, effective radius) and physical (colours, stellar masses and star formation rates) properties. Based on this analysis, we find that $\sim$88 percent of the Hydra cluster galaxies are quenched. Using the Dressler-Schectman test approach, we also find that the cluster shows possible substructures. Our analysis of the phase-space diagram together with DBSCAN algorithm indicates that Hydra shows an additional substructure that appears to be in front of the cluster centre, which is still falling into it. Our results, thus, suggest that the Hydra Cluster might not be relaxed. We analyse the median S\'ersic index as a function of wavelength and find that for red ($(u-r)\geq$2.3) and early-type galaxies it displays a slight increase towards redder filters (13 and 18 percent, for red and early-type respectively) whereas for blue+green ($(u-r)$<2.3) galaxies it remains constant. Late-type galaxies show a small decrease of the median S\'ersic index toward redder filters. Also, the S\'ersic index of galaxies, and thus their structural properties, do not significantly vary as a function of clustercentric distance and density within the cluster; and this is the case regardless of the filter.