title: A systematic search for galaxy protocluster cores at the transition epoch of their star formation activity

abstract: The redshift of $z\sim1.5$ is the transition epoch of protoclusters (PCs) from the star-forming phase into the quenching phase, and hence an appropriate era to investigate the build up of the quenched population. We define a `core' as the most massive halo in a given PC, where environmental effects are likely to work most effectively, and search for cores at $1<z<1.5$. We use a photometric redshift catalogue of a wide (effective area of $\sim22.2\,\mathrm{deg}^{2}$) and deep ($i\sim26.8\,\mathrm{mag}$) optical survey with Subaru Hyper-Suprime Cam. Regarding galaxies with $\log(M_{*}/M_{\odot})>11.3$ as the central galaxies of PC cores, we estimate their average halo mass by clustering analysis and find it to be $\log(M_\mathrm{h}/M_{\odot})\sim13.7$. An expected mass growth by the IllustrisTNG simulation and the observed overdensities around them suggest that the PC cores we find are progenitors of present-day clusters. Classifying our galaxy sample into red and blue galaxies, we calculate the stellar mass function (SMF) and the red galaxy fraction. The SMFs in the PC cores are more-top heavy than field, implying early high-mass galaxy formation and disruption of low-mass galaxies. We also find that the red fraction increases with stellar mass, consistent with stellar-mass dependent environmental quenching recently found at $z>1$. Interestingly, although the cores with red and blue centrals have similar halo masses, only those with red centrals show a significant red fraction excess compared to the field, suggesting a conformity effect. Some observational features of PC cores may imply that the conformity is caused by assembly bias.


title: The non-$\ell$-part of the number of spanning trees in abelian $\ell$-towers of multigraphs

abstract: Let $\ell$ and $p$ be two distinct primes. We study the $p$-adic valuation of the number of spanning trees in an abelian $\ell$-tower of connected multigraphs. This is analogous to the classical theorem of Washington--Sinnott on the growth of the $p$-part of the class group in a cyclotomic $\mathbb{Z}_\ell$-extension of abelian extensions of $\mathbb{Q}$. Furthermore, we show that under certain hypotheses, the number of primes dividing the number of spanning trees is unbounded in such a tower.


title: Online Slice Reconfiguration for End-to-End QoE in 6G Applications

abstract: End-to-end (E2E) quality of experience (QoE) for 6G applications depends on the synchronous allocation of networking and computing resources, also known as slicing. However, the relationship between the resources and the E2E QoE outcomes is typically stochastic and non-stationary. Existing works consider known resource demands for slicing and formulate optimization problems for slice reconfiguration. In this work, we create and manage slices by learning the relationship between E2E QoE and resources. We develop a gradient-based online slice reconfiguration algorithm (OSRA) to reconfigure and manage slices in resource-constrained scenarios for radio access networks (RAN). We observe that our methodology meets the QoE requirements with high accuracy compared to existing approaches. It improves upon the existing approaches by approximately 98\% for bursty traffic variations. Our algorithm has fast convergence and achieves low E2E delay violations for lower priority slices.


title: Charge density wave order in kagome metal AV$_3$Sb$_5$ (A= Cs, Rb, K)

abstract: We employ polarization-resolved electronic Raman spectroscopy and density functional theory to study the primary and secondary order parameters, as well as their interplay, in the charge density wave (CDW) state of the kagome metal AV3Sb5. Previous x-ray diffraction data at 15K established that the CDW order in CsV3Sb5 comprises of a 2x2x4 structure: one layer of inverse-star-of-David and three consecutive layers of star-of-David pattern. We analyze the lattice distortions based the 2x2x4 structure at 15K, and find that U lattice distortion is the primary order parameter while M and L distortions are secondary order parameters for Vanadium displacements. This conclusion is confirmed by the calculation of bare susceptibility that shows a broad peak at around qz=0.25 along the hexagonal Brillouin zone face central line (U-line). We also identify several phonon modes emerging in the CDW state, which are lattice vibration modes related to V and Sb atoms as well as alkali atoms. The detailed temperature evolution of these modes' frequencies, HWHM, and integrated intensities support a phase diagram with two successive structural phase transitions in CsV3Sb5: the first one with a primary order parameter appearing at TS=94K and the second isostructural one appearing at around 70K.


title: Roots of the identity operator and proximal mappings: (classical and phantom) cycles and gap vectors

abstract: Recently, Simons provided a lemma for a support function of a closed convex set in a general Hilbert space and used it to prove the geometry conjecture on cycles of projections. In this paper, we extend Simons's lemma to closed convex functions, show its connections to Attouch-Thera duality, and use it to characterize (classical and phantom) cycles and gap vectors of proximal mappings.


title: Persistent Extension and Analogous Bars: Data-Induced Relations Between Persistence Barcodes

abstract: A central challenge in topological data analysis is the interpretation of barcodes. The classical algebraic-topological approach to interpreting homology classes is to build maps to spaces whose homology carries semantics we understand and then to appeal to functoriality. However, we often lack such maps in real data; instead, we must rely on a cross-dissimilarity measure between our observations of a system and a reference. In this paper, we develop a pair of computational homological algebra approaches for relating persistent homology classes and barcodes: persistent extension, which enumerates potential relations between cycles from two complexes built on the same vertex set, and the method of analogous bars, which utilizes persistent extension and the witness complex built from a cross-dissimilarity measure to provide relations across systems. We provide an implementation of these methods and demonstrate their use in comparing cycles between two samples from the same metric space and determining whether topology is maintained or destroyed under clustering and dimensionality reduction.


title: Multi Order Coverage data structure to plan multi-messenger observations

abstract: We describe the use of Multi Order Coverage (MOC) maps as a practical way to manage complex regions of the sky for the planning of multi-messenger observations. MOC maps are a data structure that provides a multi-resolution representation of irregularly shaped and fragmentary regions over the sky based on the HEALPix (Hierarchical Equal Area isoLatitude Pixelization) tessellation. We present a new application of MOC, in combination with the \texttt{astroplan} observation planning package, to enable the efficient computation of sky regions and the visibility of these regions from a specific location on the Earth at a particular time. Using the example of the low-latency gravitational-wave alerts, and a simulated observational campaign with three observatories, we show that the use of MOC maps allows a high level of interoperability to support observing schedule plans. Gravitational-wave detections have an associated credible region localization on the sky. We demonstrate that these localizations can be encoded as MOC maps, and how they can be used in visualisation tools, and processed (filtered, combined) and also their utility for access to Virtual Observatory services which can be queried 'by MOC' for data within the region of interest. The ease of generating the MOC maps and the fast access to data means that the whole system can be very efficient, so that any updates on the gravitational-wave sky localization can be quickly taken into account and the corresponding adjustments to observing schedule plans can be rapidly implemented. We provide example python code as a practical example of these methods. In addition, a video demonstration of the entire workflow is available.


title: Heterodyne spectrometer sensitivity limit for quantum networking

abstract: Optical heterodyne detection-based spectrometers are attractive due to their relatively simple construction and ultra-high resolution. Here we demonstrate a proof-of-principle single-mode optical-fiber-based heterodyne spectrometer which has picometer resolution and quantum-limited sensitivity around 1550 nm. Moreover, we report a generalized quantum limit of detecting broadband multi-spectral-temporal-mode light using heterodyne detection, which provides a sensitivity limit on a heterodyne detection-based optical spectrometer. We then compare this sensitivity limit to several spectrometer types and dim light sources of interest, such as, spontaneous parametric downconversion, Raman scattering, and spontaneous four-wave mixing. We calculate the heterodyne spectrometer is significantly less sensitive than a single-photon detector and unable to detect these dim light sources, except for the brightest and narrowest-bandwidth examples.


title: `Next Generation' Reservoir Computing: an Empirical Data-Driven Expression of Dynamical Equations in Time-Stepping Form

abstract: Next generation reservoir computing based on nonlinear vector autoregression (NVAR) is applied to emulate simple dynamical system models and compared to numerical integration schemes such as Euler and the $2^\text{nd}$ order Runge-Kutta. It is shown that the NVAR emulator can be interpreted as a data-driven method used to recover the numerical integration scheme that produced the data. It is also shown that the approach can be extended to produce high-order numerical schemes directly from data. The impacts of the presence of noise and temporal sparsity in the training set is further examined to gauge the potential use of this method for more realistic applications.


title: Reverse-Engineering Information Presentations: Recovering Hierarchical Grouping from Layouts of Visual Elements

abstract: Visual elements in an information presentation are often spatially and semantically grouped hierarchically for effective message delivery. Studying the hierarchical grouping information can help researchers and designers better explore layout structures and understand design demographics. However, recovering hierarchical grouping is challenging due to a large number of possibilities for compositing visual elements into a single-page design. This paper introduces an automatic approach that takes the layout of visual elements as input and returns the hierarchical grouping as output. To understand information presentations, we first contribute a dataset of 23,072 information presentations with diverse layouts to the community. Next, we propose our technique with a Transformer-based model to predict relatedness between visual elements and a bottom-up algorithm to produce the hierarchical grouping. Finally, we evaluate our technique through a technical experiment and a user study with 30 designers. The results show that the proposed technique is promising.


title: Structure Formation in Dark Matter Particle Production Cosmology

abstract: We investigate a cosmological scenario in which the dark matter particles can be created during the evolution of the Universe. By regarding the Universe as an open thermodynamic system and using non-equilibrium thermodynamics, we examine the mechanism of gravitational particle production. In this setup, we study the large-scale structure (LSS) formation of the Universe in the Newtonian regime of perturbations and derive the equations governing the evolution of the dark matter overdensities. Then, we implement the cosmological data from Planck 2018 CMB measurements, SNe Ia and BAO observations, as well as the Riess et al. (2019) local measurement for $H_0$ to provide some cosmological constraints for the parameters of our model. We see that the best case of our scenario ($\chi_{{\rm tot}}^{2}=3834.40$) fits the observational data better than the baseline $\Lambda$CDM model ($\chi_{{\rm tot}}^{2} = 3838.00$) at the background level. We moreover estimate the growth factor of linear perturbations and show that the best case of our model ($\chi_{f\sigma_{8}}^{2}=39.85$) fits the LSS data significantly better than the $\Lambda$CDM model ($\chi_{f\sigma_{8}}^{2}=45.29$). Consequently, our model also makes a better performance at the level of the linear perturbations compared to the standard cosmological model.


title: Abstract nonlinear evolution inclusions of second order with applications in visco-elasto-plasticity

abstract: Existence of strong solutions of an abstract Cauchy problem for a class of doubly nonlinear evolution inclusion of second order is established via a semi-implicit time discretization method. The principal parts of the operators acting on $u$ and $u'$ are multi-valued subdifferential operators and are discretized implicitly. A non-variational and non-monotone perturbation acting nonlineary on $u$ and $u'$ is allowed and discretized explicitly in time. The convergence of a variational approximation scheme is established using methods from convex analysis. In addition, it is proven that the solution satisfies an energy-dissipation equality. Applications of the abstract theory to various examples, e.g., a model in visco-elastic-plasticity, are provided.


title: Aitchison's Compositional Data Analysis 40 Years On: A Reappraisal

abstract: The development of John Aitchison's approach to compositional data analysis is followed since his paper read to the Royal Statistical Society in 1982. Aitchison's logratio approach, which was proposed to solve the problematic aspects of working with data with a fixed sum constraint, is summarized and reappraised. It is maintained that the properties on which this approach was originally built, the main one being subcompositional coherence, are not required to be satisfied exactly -- quasi-coherence is sufficient, that is near enough to being coherent for all practical purposes. This opens up the field to using simpler data transformations, such as power transformations, that permit zero values in the data. The additional property of exact isometry, which was subsequently introduced and not in Aitchison's original conception, imposed the use of isometric logratio transformations, but these are complicated and problematic to interpret, involving ratios of geometric means. If this property is regarded as important in certain analytical contexts, for example unsupervised learning, it can be relaxed by showing that regular pairwise logratios, as well as the alternative quasi-coherent transformations, can also be quasi-isometric, meaning they are close enough to exact isometry for all practical purposes. It is concluded that the isometric and related logratio transformations such as pivot logratios are not a prerequisite for good practice, although many authors insist on their obligatory use. This conclusion is fully supported here by case studies in geochemistry and in genomics, where the good performance is demonstrated of pairwise logratios, as originally proposed by Aitchison, or Box-Cox power transforms of the original compositions where no zero replacements are necessary.


title: Multiple Genome Analytics Framework: The Case of All SARS-CoV-2 Complete Variants

abstract: Pattern detection and string matching are fundamental problems in computer science and the accelerated expansion of bioinformatics and computational biology have made them a core topic for both disciplines. The SARS-CoV-2 pandemic has made such problems more demanding with hundreds or thousands of new genome variants discovered every week, because of constant mutations, and there is a desperate need for fast and accurate analyses. The requirement for computational tools for genomic analyses, such as sequence alignment, is very important, although, in most cases the resources and computational power required are enormous. The presented Multiple Genome Analytics Framework combines data structures and algorithms, specifically built for text mining and pattern detection, that can help to efficiently address several computational biology and bioinformatics problems concurrently with minimal resources. A single execution of advanced algorithms, with space and time complexity O(nlogn), is enough to acquire knowledge on all repeated patterns that exist in multiple genome sequences and this information can be used from other meta-algorithms for further meta-analyses. The potential of the proposed framework is demonstrated with the analysis of more than 300,000 SARS-CoV-2 genome sequences and the detection of all repeated patterns with length up to 60 nucleotides in these sequences. These results have been used to provide answers to questions such as common patterns among all variants, sequence alignment, palindromes and tandem repeats detection, different organism genome comparisons, polymerase chain reaction primers detection, etc.


title: D-Box: DMA-enabled Compartmentalization for Embedded Applications

abstract: Embedded and Internet-of-Things (IoT) devices have seen an increase in adoption in many domains. The security of these devices is of great importance as they are often used to control critical infrastructure, medical devices, and vehicles. Existing solutions to isolate microcontroller (MCU) resources in order to increase their security face significant challenges such as specific hardware unavailability, Memory Protection Unit (MPU) limitations and a significant lack of Direct Memory Access (DMA) support. Nevertheless, DMA is fundamental for the power and performance requirements of embedded applications. In this paper, we present D-Box, a systematic approach to enable secure DMA operations for compartmentalization solutions of embedded applications using real-time operating systems (RTOS). D-Box defines a reference architecture and a workflow to protect DMA operations holistically. It provides practical methods to harden the kernel and define capability-based security policies for easy definition of DMA operations with strong security properties. We implemented a D-Box prototype for the Cortex-M3/M4 on top of the popular FreeRTOS-MPU (F-MPU). The D-Box procedures and a stricter security model enabled DMA operations, yet it exposed 41 times less ROP (return-orienting-programming) gadgets when compared with the standard F-MPU. D-Box adds only a 2% processor overhead while reducing the power consumption of peripheral operation benchmarks by 18.2%. The security properties and performance of D-Box were tested and confirmed on a real-world case study of a Programmable Logic Controller (PLC) application.


title: Anti-Holomorphic Modes in Vortex Lattices

abstract: A continuum theory of linearized Helmholtz-Kirchoff point vortex dynamics about a steadily rotating lattice state is developed by two separate methods: firstly by a direct procedure, secondly by taking the long-wavelength limit of Tkachenko's exact solution for a triangular vortex lattice. Solutions to the continuum theory are found, described by arbitrary anti-holomorphic functions, and give power-law localized edge modes. Numerical results for finite lattices show excellent agreement to the theory.

categories: physics.flu-dyn cond-mat.quant-gas cond-mat.soft math-ph math.MP


title: A Tight Reverse Minkowski Inequality for the Epstein Zeta Function

abstract: We prove that if $\mathcal{L} \subset \mathbb{R}^n$ is a lattice such that $\det(\mathcal{L}') \geq 1$ for all sublattices $\mathcal{L}' \subseteq \mathcal{L}$, then \[ \sum_{\substack{\mathbf{y}\in\mathcal{L}\\\mathbf{y}\neq\mathbf0}} (\|\mathbf{y}\|^2+q)^{-s} \leq \sum_{\substack{\mathbf{z} \in \mathbb{Z}^n\\\mathbf{z}\neq\mathbf{0}}} (\|\mathbf{z}\|^2+q)^{-s} \] for all $s > n/2$ and all $0 \leq q \leq (2s-n)/(n+2)$, with equality if and only if $\mathcal{L}$ is isomorphic to $\mathbb{Z}^n$.


title: First-order continuation method for steady-state variably saturated groundwater flow modeling

abstract: Recently, the nonlinearity continuation method has been used to numerically solve boundary value problems for steady-state Richards equation. The method can be considered as a predictor-corrector procedure with the simplest form which has been applied to date having a trivial, zeroth-order predictor. In this article, effect of a more sophisticated predictor technique is examined. Numerical experiments are performed with finite volume and mimetic finite difference discretizations on various problems, including real-life examples.


title: An Intelligent System for Multi-topic Social Spam Detection in Microblogging

abstract: The communication revolution has perpetually reshaped the means through which people send and receive information. Social media is an important pillar of this revolution and has brought profound changes to various aspects of our lives. However, the open environment and popularity of these platforms inaugurate windows of opportunities for various cyber threats, thus social networks have become a fertile venue for spammers and other illegitimate users to execute their malicious activities. These activities include phishing hot and trendy topics and posting a wide range of contents in many topics. Hence, it is crucial to continuously introduce new techniques and approaches to detect and stop this category of users. This paper proposes a novel and effective approach to detect social spammers. An investigation into several attributes to measure topic-dependent and topic-independent users' behaviours on Twitter is carried out. The experiments of this study are undertaken on various machine learning classifiers. The performance of these classifiers are compared and their effectiveness is measured via a number of robust evaluation measures. Further, the proposed approach is benchmarked against state-of-the-art social spam and anomalous detection techniques. These experiments report the effectiveness and utility of the proposed approach and embedded modules.


title: Sharp estimates on random hyperplane tessellations

abstract: We study the problem of generating a hyperplane tessellation of an arbitrary set $T$ in $\mathbb{R}^n$, ensuring that the Euclidean distance between any two points corresponds to the fraction of hyperplanes separating them up to a pre-specified error $\delta$. We focus on random gaussian tessellations with uniformly distributed shifts and derive sharp bounds on the number of hyperplanes $m$ that are required. Surprisingly, our lower estimates falsify the conjecture that $m\sim \ell_*^2(T)/\delta^2$, where $\ell_*^2(T)$ is the gaussian width of $T$, is optimal.


title: Azimuthal Asymmetries in Unpolarised Semi-Inclusive DIS at COMPASS

abstract: In semi-inclusive deep inelastic scattering (SIDIS) the non-zero transverse momentum of partons induces azimuthal dependence of the cross-section. For unpolarised nucleon, three azimuthal modulations that can be related to different combinations of twist-two or higher twist transverse momentum dependent PDFs and fragmentation functions arise: the so-called Cahn effect reflected in $\cos\phi_\mathrm{h}$ modulation, the $\cos2\phi_\mathrm{h}$ term related to the Boer--Mulders PDF and $\sin\phi_\mathrm{h}$ twist-three effect known as beam-spin asymmetry. In 2016 and 2017, the COMPASS experiment at CERN collected a large sample of SIDIS events using a longitudinally polarised 160~GeV/$c$ muon beam scattering on a liquid hydrogen target. Amplitudes of the aforementioned azimuthal modulations have been extracted from part of the data. A new procedure has been developed to subtract a background coming from the decay of diffractively produced vector mesons. The results presented in this talk qualitatively agree with earlier COMPASS results obtained with an isoscalar target.


title: Reproducible, incremental representation learning with Rosetta VAE

abstract: Variational autoencoders are among the most popular methods for distilling low-dimensional structure from high-dimensional data, making them increasingly valuable as tools for data exploration and scientific discovery. However, unlike typical machine learning problems in which a single model is trained once on a single large dataset, scientific workflows privilege learned features that are reproducible, portable across labs, and capable of incrementally adding new data. Ideally, methods used by different research groups should produce comparable results, even without sharing fully trained models or entire data sets. Here, we address this challenge by introducing the Rosetta VAE (R-VAE), a method of distilling previously learned representations and retraining new models to reproduce and build on prior results. The R-VAE uses post hoc clustering over the latent space of a fully-trained model to identify a small number of Rosetta Points (input, latent pairs) to serve as anchors for training future models. An adjustable hyperparameter, $\rho$, balances fidelity to the previously learned latent space against accommodation of new data. We demonstrate that the R-VAE reconstructs data as well as the VAE and $\beta$-VAE, outperforms both methods in recovery of a target latent space in a sequential training setting, and dramatically increases consistency of the learned representation across training runs.


title: First-principles simulation of light-ion microscopy of graphene

abstract: The extreme sensitivity of 2D materials to defects and nanostructure requires precise imaging techniques to verify presence of desirable and absence of undesirable features in the atomic geometry. Helium-ion beams have emerged as a promising materials imaging tool, achieving up to 20 times higher resolution and 10 times larger depth-of-field than conventional or environmental scanning electron microscopes. Here, we offer first-principles theoretical insights to advance ion-beam imaging of atomically thin materials by performing real-time time-dependent density functional theory simulations of single impacts of 10-200 keV light ions in free-standing graphene. We predict that detecting electrons emitted from the back of the material (the side from which the ion exits) would result in up to 3 times higher signal and up to 5 times higher contrast images, making 2D materials especially compelling targets for ion-beam microscopy. We also find that the charge induced in the graphene equilibrates on a sub-fs time scale, leading to only slight disturbances in the carbon lattice that are unlikely to damage the atomic structure for any of the beam parameters investigated here.


title: The Pad\'e matrix pencil method with spurious pole information assimilation

abstract: We present a novel method for calculating Pad\'e approximants that is capable of eliminating spurious poles placed at the point of development and of identifying and eliminating spurious poles created by precision limitations and/or noisy coefficients. Information contained in in the eliminated poles is assimilated producing a reduced order Pad\'e approximant (PA). While the [m+k/m] conformation produced by the algorithm is flexible, the m value of the rational approximant produced by the algorithm reported here is determined by the number of spurious poles eliminated. Spurious poles due to coefficient noise/precision limitations are identified using an evidence-based filter parameter applied to the singular values of a matrix comprised of the series coefficients. The rational function poles are found directly by solving a generalized eigenvalue problem defined by a matrix pencil. Spurious poles place at the point of development, responsible in some algorithms for degeneracy, are identified by their magnitudes. Residues are found by solving an overdetermined linear matrix equation. The method is compared with the so-called Robust Pad\'e Approximation (RPA) method and shown to be competitive on the problems studied. By eliminating spurious poles, particularly in functions with branch points, such as those encountered solving the power-flow problem, solution of these complex-valued problems is made more reliable.


title: New interior model of neutron stars

abstract: The Tolman VII solution is considered by some as one of the few analytical solutions to Einstein's equations, which describes approximately well the interior of neutron stars (NSs). This solution is characterized by the mass $M$, radius $R$, and an energy density that varies quadratically with the radial coordinate $r$. Recently, Jiang and Yagi proposed a modification of this solution, the so-called modified Tolman VII (MTVII) solution, by introducing an additional quartic term to the energy density radial profile. The MTVII solution is an approximate solution to Einstein's equation, which includes a new parameter $\alpha$ that allows the solution to have a better agreement with the energy density profiles for realistic NSs. Here we consider the MTVII solution, showing that for certain values of the parameter $\alpha$ and compactness $\mathcal{C}$ this solution manifests a region of negative pressure near the surface which leads to negative values of the tidal Love number. To alleviate these drawbacks, we introduce an exact version of the MTVII solution obtained by solving numerically Einstein's equations for the MTVII energy density profile. As an application of our new exact MTVII (EMTVII) solution, we calculate the tidal Love number and tidal deformability, as a function of $\mathcal{C}$, for different values of the parameter $\alpha$. We find that the EMTVII solution predicts a positive tidal Love number for the whole range of allowed values of parameters $(\mathcal{C},\alpha)$, in agreement with previous results for realistic NSs.


title: Semilinear wave equations on accelerated expanding FLRW spacetimes

abstract: We identify a large class of systems of semilinear wave equations, on fixed accelerated expanding FLRW spacetimes, with nearly at spatial slices, for which we prove small data future global well-posedness. The family of systems we consider is large in the sense that, among other examples, it includes general wave maps, as well as natural generalizations of some of Fritz John's "blow up" equations (whose future blow up disappears, in our setting, as a consequence of the spacetime expansion). We also establish decay upper bounds, which are sharp within the family of systems under analysis.


title: Coexistence of Multiple Stacking Charge Density Waves in Kagome Superconductor ${\mathrm{CsV}}_3{\mathrm{Sb}}_5$

abstract: The recently discovered Kagome family ${\mathrm{AV}}_3{\mathrm{Sb}}_5$ (A = K, Rb, Cs) exhibits rich physical phenomena, including non-trivial topological electronic structure, giant anomalous Hall effect, charge density waves (CDW) and superconductivity. Notably, CDW in ${\mathrm{AV}}_3{\mathrm{Sb}}_5$ is evidenced to intertwine with its superconductivity and topology, but its nature remains elusive. Here, we combine x-ray scattering experiments and density-functional theory calculations to investigate the CDWs in ${\mathrm{CsV}}_3{\mathrm{Sb}}_5$ and demonstrate the coexistence of 2 $\times$ 2 $\times$ 2 and 2 $\times$ 2 $\times$ 4 CDW stacking phases. Competition between these CDW phases is revealed by tracking the temperature evolution of CDW intensities, which also manifests in different transition temperatures during warming- and cooling measurements. We also identify a meta-stable quenched state of ${\mathrm{CsV}}_3{\mathrm{Sb}}_5$ after fast-cooling process. Our study demonstrates the coexistence of competing CDW stacking in ${\mathrm{CsV}}_3{\mathrm{Sb}}_5$, offering new insights in understanding the novel properties of this system.


title: Probabilistic design of optimal sequential decision-making algorithms in learning and control

abstract: This survey is focused on certain sequential decision-making problems that involve optimizing over probability functions. We discuss the relevance of these problems for learning and control. The survey is organized around a framework that combines a problem formulation and a set of resolution methods. The formulation consists of an infinite-dimensional optimization problem. The methods come from approaches to search optimal solutions in the space of probability functions. Through the lenses of this overarching framework we revisit popular learning and control algorithms, showing that these naturally arise from suitable variations on the formulation mixed with different resolution methods. A running example, for which we make the code available, complements the survey. Finally, a number of challenges arising from the survey are also outlined.


title: Parallel Neural Local Lossless Compression

abstract: The recently proposed Neural Local Lossless Compression (NeLLoC), which is based on a local autoregressive model, has achieved state-of-the-art (SOTA) out-of-distribution (OOD) generalization performance in the image compression task. In addition to the encouragement of OOD generalization, the local model also allows parallel inference in the decoding stage. In this paper, we propose two parallelization schemes for local autoregressive models. We discuss the practicalities of implementing the schemes and provide experimental evidence of significant gains in compression runtime compared to the previous, non-parallel implementation.


title: How I learned to stop worrying and love the curse of dimensionality: an appraisal of cluster validation in high-dimensional spaces

abstract: The failure of the Euclidean norm to reliably distinguish between nearby and distant points in high dimensional space is well-known. This phenomenon of distance concentration manifests in a variety of data distributions, with iid or correlated features, including centrally-distributed and clustered data. Unsupervised learning based on Euclidean nearest-neighbors and more general proximity-oriented data mining tasks like clustering, might therefore be adversely affected by distance concentration for high-dimensional applications. While considerable work has been done developing clustering algorithms with reliable high-dimensional performance, the problem of cluster validation--of determining the natural number of clusters in a dataset--has not been carefully examined in high-dimensional problems. In this work we investigate how the sensitivities of common Euclidean norm-based cluster validity indices scale with dimension for a variety of synthetic data schemes, including well-separated and noisy clusters, and find that the overwhelming majority of indices have improved or stable sensitivity in high dimensions. The curse of dimensionality is therefore dispelled for this class of fairly generic data schemes.


title: Depth optimization of CZ, CNOT, and Clifford circuits

abstract: We seek to develop better upper bound guarantees on the depth of quantum CZ gate, CNOT gate, and Clifford circuits than those reported previously. We focus on the number of qubits $n\,{\leq}\,$1,345,000 [1], which represents the most practical use case. Our upper bound on the depth of CZ circuits is $\lfloor n/2 + 0.4993{\cdot}\log^2(n) + 3.0191{\cdot}\log(n) - 10.9139\rfloor$, improving best known depth by a factor of roughly 2. We extend the constructions used to prove this upper bound to obtain depth upper bound of $\lfloor n + 1.9496{\cdot}\log^2(n) + 3.5075{\cdot}\log(n) - 23.4269 \rfloor$ for CNOT gate circuits, offering an improvement by a factor of roughly $4/3$ over state of the art, and depth upper bound of $\lfloor 2n + 2.9487{\cdot}\log^2(n) + 8.4909{\cdot}\log(n) - 44.4798\rfloor$ for Clifford circuits, offering an improvement by a factor of roughly $5/3$.


title: The Fairness Field Guide: Perspectives from Social and Formal Sciences

abstract: Over the past several years, a slew of different methods to measure the fairness of a machine learning model have been proposed. However, despite the growing number of publications and implementations, there is still a critical lack of literature that explains the interplay of fair machine learning with the social sciences of philosophy, sociology, and law. We hope to remedy this issue by accumulating and expounding upon the thoughts and discussions of fair machine learning produced by both social and formal (specifically machine learning and statistics) sciences in this field guide. Specifically, in addition to giving the mathematical and algorithmic backgrounds of several popular statistical and causal-based fair machine learning methods, we explain the underlying philosophical and legal thoughts that support them. Further, we explore several criticisms of the current approaches to fair machine learning from sociological and philosophical viewpoints. It is our hope that this field guide will help fair machine learning practitioners better understand how their algorithms align with important humanistic values (such as fairness) and how we can, as a field, design methods and metrics to better serve oppressed and marginalized populaces.


title: Learning Enhancement of CNNs via Separation Index Maximizing at the First Convolutional Layer

abstract: In this paper, a straightforward enhancement learning algorithm based on Separation Index (SI) concept is proposed for Convolutional Neural Networks (CNNs). At first, the SI as a supervised complexity measure is explained its usage in better learning of CNNs for classification problems illustrate. Then, a learning strategy proposes through which the first layer of a CNN is optimized by maximizing the SI, and the further layers are trained through the backpropagation algorithm to learn further layers. In order to maximize the SI at the first layer, A variant of ranking loss is optimized by using the quasi least square error technique. Applying such a learning strategy to some known CNNs and datasets, its enhancement impact in almost all cases is demonstrated.


title: The Ideal Membership Problem and Abelian Groups

abstract: Given polynomials $f_0,\dots, f_k$ the Ideal Membership Problem, IMP for short, asks if $f_0$ belongs to the ideal generated by $f_1,\dots, f_k$. In the search version of this problem the task is to find a proof of this fact. The IMP is a well-known fundamental problem with numerous applications, for instance, it underlies many proof systems based on polynomials such as Nullstellensatz, Polynomial Calculus, and Sum-of-Squares. Although the IMP is in general intractable, in many important cases it can be efficiently solved. Mastrolilli [SODA'19] initiated a systematic study of IMPs for ideals arising from Constraint Satisfaction Problems (CSPs), parameterized by constraint languages, denoted IMP($\Gamma$). The ultimate goal of this line of research is to classify all such IMPs accordingly to their complexity. Mastrolilli achieved this goal for IMPs arising from CSP($\Gamma$) where $\Gamma$ is a Boolean constraint language, while Bulatov and Rafiey [ArXiv'21] advanced these results to several cases of CSPs over finite domains. In this paper we consider IMPs arising from CSPs over `affine' constraint languages, in which constraints are subgroups (or their cosets) of direct products of Abelian groups. This kind of CSPs include systems of linear equations and are considered one of the most important types of tractable CSPs. Some special cases of the problem have been considered before by Bharathi and Mastrolilli [MFCS'21] for linear equation modulo 2, and by Bulatov and Rafiey [ArXiv'21] to systems of linear equations over $GF(p)$, $p$ prime. Here we prove that if $\Gamma$ is an affine constraint language then IMP($\Gamma$) is solvable in polynomial time assuming the input polynomial has bounded degree.


title: Continuous limits of large plant-pollinator random networks and some applications

abstract: We study a stochastic individual-based model of interacting plant and pollinator species through a bipartite graph: each species is a node of the graph, an edge representing interactions between a pair of species. The dynamics of the system depends on the between- and within-species interactions: pollination by insects increases plant reproduction rate but has a cost which can increase plant death rate, depending on pollinators density. Pollinators reproduction is increased by the resources harvested on plants. Each species is characterized by a trait corresponding to its degree of generalism. This trait determines the structure of the interactions graph and the quantity of resources exchanged between species. Our model includes in particular nested or modular networks. Deterministic approximations of the stochastic measure-valued process by systems of ordinary differential equations or integro-differential equations are established and studied, when the population is large or when the graph is dense and can be replaced with a graphon. The long-time behaviors of these limits are studied and central limit theorems are established to quantify the difference between the discrete stochastic individual-based model and the deterministic approximations. Finally, studying the continuous limits of the interaction network and the resulting PDEs, we show that nested plant-pollinator communities are expected to collapse towards a coexistence between a single pair of species of plants and pollinators.


title: Implementation of a generalized precession parameter in the RIFT parameter estimation algorithm

abstract: Since the initial discovery of gravitational waves in 2015, significant developments have been made towards waveform interpretation and estimation of compact binary source parameters. We present herein an implementation of the generalized precession parameter $\langle \chi_p \rangle$, which averages over all angular variations on the precession timescale, within the RIFT parameter estimation framework. Relative to the originally-proposed precession parameter $\chi_p$, which characterizes the single largest dynamical spin in a binary, this new parameter $\langle \chi_p \rangle$ has a unique domain $1 < \langle \chi_p \rangle < 2$, which is exclusive to binaries with two precessing spins. After reviewing the physical differences between these two parameters, we describe how $\langle \chi_p \rangle$ was implemented in RIFT and apply it to all 36 events from the second half of the Advanced LIGO and Advanced Virgo third operating run (O3b). In O3b, ten events show significant amounts of precession $\langle \chi_p \rangle > 0.5$. Of particular interest is GW191109_010717; we show it has a $\sim28\%$ probability that the originating system necessarily contains two misaligned spins.


title: Using Survey Data to Obtain More Representative Site Samples for Impact Studies

abstract: To improve the generalizability of impact evaluations, recent research has examined statistical methods for selecting representative samples of sites. However, these methods rely on having rich data on impact moderators for all sites in the target population. This paper offers a new approach to selecting sites for impact studies when rich data on impact moderators are available, but only from a survey based on a representative sample of the impact study's target population. Survey data are used to (1) estimate the proportion of sites in the population with certain characteristics, and (2) set limits on the number of sites with different characteristics that the sample can include. The Principal Investigator enforces the limits to ensure that certain types of sites are not overrepresented in the final sample. These limits can be layered on top of site selection and recruitment approaches to improve the representativeness of the sample.


title: Assemble Foundation Models for Automatic Code Summarization

abstract: Automatic code summarization is beneficial to daily software development since it could help reduce the requirement of manual writing. Currently, artificial intelligence is undergoing a paradigm shift. The foundation models pretrained on massive data and finetuned to downstream tasks surpass specially customized models. This trend inspired us to consider reusing foundation models instead of learning from scratch. Thereby, we propose a flexible and robust approach for automatic code summarization, based on neural models. We assemble available foundation models, such as CodeBERT and GPT-2, into a single neural model named AdaMo. Moreover, we utilize Gaussian noise as the simulation of contextual information to optimize the latent representation. Furthermore, we introduce two adaptive schemes from the perspective of knowledge transfer, namely continuous pretraining and intermediate finetuning, and design intermediate stage tasks for general sequence-to-sequence learning. Finally, we evaluate AdaMo against a benchmark dataset for code summarization, by comparing it with state-of-the-art models.


title: Time reversal of spinal processes for linear and non-linear branching processes near stationarity

abstract: We consider a stochastic individual-based population model with competition, trait-structure affecting reproduction and survival, and changing environment. The changes of traits are described by jump processes, and the dynamics can be approximated in large population by a non-linear PDE with a non-local mutation operator. Using the fact that this PDE admits a non-trivial stationary solution, we can approximate the non-linear stochastic population process by a linear birth-death process where the interactions are frozen, as long as the population remains close to this equilibrium. This allows us to derive, when the population is large, the equation satisfied by the ancestral lineage of an individual uniformly sampled at a fixed time $T$, which is the path constituted of the traits of the ancestors of this individual in past times $t\leq T$. This process is a time inhomogeneous Markov process, but we show that the time reversal of this process possesses a very simple structure (e.g. time-homogeneous and independent of $T$). This extends recent results where the authors studied a similar model with a Laplacian operator but where the methods essentially relied on the Gaussian nature of the mutations.


title: The FACET Project: Forward Aperture CMS ExTension to search for new Long-Lived Particles

abstract: FACET is a proposed new subsysem for CMS to search for portals such as dark photons, dark higgs, heavy neutral leptons and axion-like particles in the very forward direction at the High Luminosity LHC. Such particles can penetrate up to 50 m of iron and then decay inside a 14 m$^3$ vacuum pipe made by enlarging an 18 m long section of the LHC pipe to a radius of 50 cm.


title: Learning-Based MIMO Channel Estimation under Spectrum Efficient Pilot Allocation and Feedback

abstract: Wireless links using massive MIMO transceivers are vital for next generation wireless communications networks networks. Precoding in Massive MIMO transmission requires accurate downlink channel state information (CSI). Many recent works have effectively applied deep learning (DL) to jointly train UE-side compression networks for delay domain CSI and a BS-side decoding scheme. Vitally, these works assume that the full delay domain CSI is available at the UE, but in reality, the UE must estimate the delay domain based on a limited number of frequency domain pilots. In this work, we propose a linear pilot-to-delay (P2D) estimator that transforms sparse frequency pilots to the truncated delay CSI. We show that the P2D estimator is accurate under frequency downsampling, and we demonstrate that the P2D estimate can be effectively utilized with existing autoencoder-based CSI estimation networks. In addition to accounting for pilot-based estimates of downlink CSI, we apply unrolled optimization networks to emulate iterative solutions to compressed sensing (CS), and we demonstrate better estimation performance than prior autoencoder-based DL networks. Finally, we investigate the efficacy of trainable CS networks for in a differential encoding network for time-varying CSI estimation, and we propose a new network, MarkovNet-ISTA-ENet, comprised of both a CS network for initial CSI estimation and multiple autoencoders to estimate the error terms. We demonstrate that this heterogeneous network has better asymptotic performance than networks comprised of only one type of network.


title: Towards a Data Privacy-Predictive Performance Trade-off

abstract: Machine learning is increasingly used in the most diverse applications and domains, whether in healthcare, to predict pathologies, or in the financial sector to detect fraud. One of the linchpins for efficiency and accuracy in machine learning is data utility. However, when it contains personal information, full access may be restricted due to laws and regulations aiming to protect individuals' privacy. Therefore, data owners must ensure that any data shared guarantees such privacy. Removal or transformation of private information (de-identification) are among the most common techniques. Intuitively, one can anticipate that reducing detail or distorting information would result in losses for model predictive performance. However, previous work concerning classification tasks using de-identified data generally demonstrates that predictive performance can be preserved in specific applications. In this paper, we aim to evaluate the existence of a trade-off between data privacy and predictive performance in classification tasks. We leverage a large set of privacy-preserving techniques and learning algorithms to provide an assessment of re-identification ability and the impact of transformed variants on predictive performance. Unlike previous literature, we confirm that the higher the level of privacy (lower re-identification risk), the higher the impact on predictive performance, pointing towards clear evidence of a trade-off.


title: Dark Energy Survey Year 3 Results: Three-Point Shear Correlations and Mass Aperture Moments

abstract: We present high signal-to-noise measurements of three-point shear correlations and the third moment of the mass aperture statistic using the first 3 years of data from the Dark Energy Survey. We additionally obtain the first measurements of the configuration and scale dependence of the four three-point shear correlations which carry cosmological information. With the third-order mass aperture statistic, we present tomographic measurements over angular scales of 4 to 60 arcminutes with a combined statistical significance of 15.0$\sigma$. Using the tomographic information and measuring also the second-order mass aperture, we additionally obtain a skewness parameter and its redshift evolution. We find that the amplitudes and scale-dependence of these shear 3pt functions are in qualitative agreement with measurements in a mock galaxy catalog based on N-body simulations, indicating promise for including them in future cosmological analyses. We validate our measurements by showing that B-modes, parity-violating contributions and PSF modeling uncertainties are negligible, and determine that the measured signals are likely to be of astrophysical and gravitational origin.


title: Stellar dating using chemical clocks and Bayesian inference

abstract: Dating stars is a major challenge with a deep impact on many astrophysical fields. One of the most promising techniques for this is using chemical abundances. Recent space- and ground-based facilities have improved the quantity of stars with accurate observations. This has opened the door for using Bayesian inference tools to maximise the information we can extract from them. Our aim is to present accurate and reliable stellar age estimates of FGK stars using chemical abundances and stellar parameters. We used one of the most flexible Bayesian inference techniques (hierarchical Bayesian models) to exceed current possibilities in the use of chemical abundances for stellar dating. Our model is a data-driven model. We used a training set that has been presented in the literature with ages estimated with isochrones and accurate stellar abundances and general characteristics. The core of the model is a prescription of certain abundance ratios as linear combinations of stellar properties including age. We gathered four different testing sets to assess the accuracy, precision, and limits of our model. We also trained a model using chemical abundances alone. We found that our age estimates and those coming from asteroseismology, other accurate sources, and also with ten Gaia benchmark stars agree well. The mean absolute difference of our estimates compared with those used as reference is 0.9 Ga, with a mean difference of 0.01 Ga. When using open clusters, we reached a very good agreement for Hyades, NGC 2632, Ruprecht 147, and IC4651. We also found outliers that are a reflection of chemical peculiarities and/or stars at the limit of the validity ranges of the training set. The model that only uses chemical abundances shows slightly worse mean absolute difference (1.18 Ga) and mean difference (-0.12 Ga).


title: Examining and Mitigating the Impact of Crossbar Non-idealities for Accurate Implementation of Sparse Deep Neural Networks

abstract: Recently several structured pruning techniques have been introduced for energy-efficient implementation of Deep Neural Networks (DNNs) with lesser number of crossbars. Although, these techniques have claimed to preserve the accuracy of the sparse DNNs on crossbars, none have studied the impact of the inexorable crossbar non-idealities on the actual performance of the pruned networks. To this end, we perform a comprehensive study to show how highly sparse DNNs, that result in significant crossbar-compression-rate, can lead to severe accuracy losses compared to unpruned DNNs mapped onto non-ideal crossbars. We perform experiments with multiple structured-pruning approaches (such as, C/F pruning, XCS and XRS) on VGG11 and VGG16 DNNs with benchmark datasets (CIFAR10 and CIFAR100). We propose two mitigation approaches - Crossbar column rearrangement and Weight-Constrained-Training (WCT) - that can be integrated with the crossbar-mapping of the sparse DNNs to minimize accuracy losses incurred by the pruned models. These help in mitigating non-idealities by increasing the proportion of low conductance synapses on crossbars, thereby improving their computational accuracies.


title: NLP in Human Rights Research -- Extracting Knowledge Graphs About Police and Army Units and Their Commanders

abstract: In this working paper we explore the use of an NLP system to assist the work of Security Force Monitor (SFM). SFM creates data about the organizational structure, command personnel and operations of police, army and other security forces, which assists human rights researchers, journalists and litigators in their work to help identify and bring to account specific units and personnel alleged to have committed abuses of human rights and international criminal law. This working paper presents an NLP system that extracts from English language news reports the names of security force units and the biographical details of their personnel, and infers the formal relationship between them. Published alongside this working paper are the system's code and training dataset. We find that the experimental NLP system performs the task at a fair to good level. Its performance is sufficient to justify further development into a live workflow that will give insight into whether its performance translates into savings in time and resource that would make it an effective technical intervention.


title: Contextual Bandits for Advertising Campaigns: A Diffusion-Model Independent Approach (Extended Version)

abstract: Motivated by scenarios of information diffusion and advertising in social media, we study an influence maximization problem in which little is assumed to be known about the diffusion network or about the model that determines how information may propagate. In such a highly uncertain environment, one can focus on multi-round diffusion campaigns, with the objective to maximize the number of distinct users that are influenced or activated, starting from a known base of few influential nodes. During a campaign, spread seeds are selected sequentially at consecutive rounds, and feedback is collected in the form of the activated nodes at each round. A round's impact (reward) is then quantified as the number of newly activated nodes. Overall, one must maximize the campaign's total spread, as the sum of rounds' rewards. In this setting, an explore-exploit approach could be used to learn the key underlying diffusion parameters, while running the campaign. We describe and compare two methods of contextual multi-armed bandits, with upper-confidence bounds on the remaining potential of influencers, one using a generalized linear model and the Good-Turing estimator for remaining potential (GLM-GT-UCB), and another one that directly adapts the LinUCB algorithm to our setting (LogNorm-LinUCB). We show that they outperform baseline methods using state-of-the-art ideas, on synthetic and real-world data, while at the same time exhibiting different and complementary behavior, depending on the scenarios in which they are deployed.


title: FARSI: Facebook AR System Investigator for Agile Domain-Specific System-on-Chip Exploration

abstract: Domain-specific SoCs (DSSoCs) are attractive solutions for domains with stringent power/performance/area constraints; however, they suffer from two fundamental complexities. On the one hand, their many specialized hardware blocks result in complex systems and thus high development effort. On the other, their many system knobs expand the complexity of design space, making the search for the optimal design difficult. Thus to reach prevalence, taming such complexities is necessary. This work identifies necessary features of an early-stage design space exploration (DSE) framework that targets the complex design space of DSSoCs and further provides an instance of one called FARSI, (F)acebook (AR) (S)ystem (I)nvestigator. Concretely, FARSI provides an agile system-level simulator with speed up and accuracy of 8,400X and 98.5% comparing to Synopsys Platform Architect. FARSI also provides an efficient exploration heuristic and achieves up to 16X improvementin convergence time comparing to naive simulated annealing (SA). This is done by augmenting SA with architectural reasoning such as locality exploitation and bottleneck relaxation. Furthermore, we embed various co-design capabilities and show that on average, they have a 32% impact on the convergence rate. Finally, we demonstrate that using simple development-cost-aware policies can lower the system complexity, both in terms of the component count and variation by as much as 150% and 118% (e,g., for Network-on-a-Chip subsystem)


title: Density reconstruction from schlieren images through Bayesian nonparametric models

abstract: This study proposes a radically alternate approach for extracting quantitative information from schlieren images. The method uses a scaled, derivative enhanced Gaussian process model to obtain true density estimates from two corresponding schlieren images with the knife-edge at horizontal and vertical orientations. We illustrate our approach on schlieren images taken from a wind tunnel sting model, a supersonic aircraft in flight, and a high-order numerical shock tube simulation.


title: Optimal alphabet for single text compression

abstract: A text written using symbols from a given alphabet can be compressed using the Huffman code, which minimizes the length of the encoded text. It is necessary, however, to employ a text-specific codebook, i.e. the symbol-codeword dictionary, to decode the original text. Thus, the compression performance should be evaluated by the full code length, i.e. the length of the encoded text plus the length of the codebook. We studied several alphabets for compressing texts -- letters, n-grams of letters, syllables, words, and phrases. If only sufficiently short texts are retained, an alphabet of letters or two-grams of letters is optimal. For the majority of Project Gutenberg texts, the best alphabet (the one that minimizes the full code length) is given by syllables or words, depending on the representation of the codebook. Letter 3 and 4-grams, having on average comparable length to syllables/words, perform noticeably worse than syllables or words. Word 2-grams also are never the best alphabet, on the account of having a very large codebook. We also show that the codebook representation is important -- switching from a naive representation to a compact one significantly improves the matters for alphabets with large number of symbols, most notably the words. Thus, meaning-expressing elements of the language (syllables or words) provide the best compression alphabet.


title: Well-posedness of a fully nonlinear evolution inclusion of second order

abstract: The well-posedness of the abstract \textsc{Cauchy} problem for the doubly nonlinear evolution inclusion equation of second order \begin{align*} \begin{cases} u''(t)+\partial \Psi(u'(t))+B(t,u(t))\ni f(t), &\quad t\in (0,T),\, T>0,\\ u(0)=u_0, \quad u'(0)=v_0 \end{cases} \end{align*} in a real separable \textsc{Hilbert} space $\mathscr{H}$, where $u_0\in \mathscr{H}, v_0\in \overline{D(\partial \Psi)}\cap D(\Psi), f\in L^2(0,T;\mathscr{H})$. The functional $\Psi: \mathscr{H} \rightarrow (-\infty,+\infty]$ is supposed to be proper, lower semicontinuous, and convex and the nonlinear operator $B:[0,T]\times \mathscr{H}\rightarrow \mathscr{H}$ is supposed to satisfy a (local) \textsc{Lipschitz} condition. Existence and uniqueness of strong solutions $u\in H^2(0,T^*;\mathscr{H})$ as well as the continuous dependence of solutions from the data re shown by employing the theory of nonlinear semigroups and the Banach fixed-point theorem. If $B$ satisfies a local Lipschitz condition, then the existence of strong local solutions are obtained.


title: A Method for Controlling Extrapolation when Visualizing and Optimizing the Prediction Profiles of Statistical and Machine Learning Models

abstract: We present a novel method for controlling extrapolation in the prediction profiler in the JMP software. The prediction profiler is a graphical tool for exploring high dimensional prediction surfaces for statistical and machine learning models. The profiler contains interactive cross-sectional views, or profile traces, of the prediction surface of a model. Our method helps users avoid exploring predictions that should be considered extrapolation. It also performs optimization over a constrained factor region that avoids extrapolation using a genetic algorithm. In simulations and real world examples, we demonstrate how optimal factor settings without constraint in the profiler are frequently extrapolated, and how extrapolation control helps avoid these solutions with invalid factor settings that may not be useful to the user.


title: Imaging dark charge emitters in diamond via carrier-to-photon conversion

abstract: The application of color centers in wide-bandgap semiconductors to nanoscale sensing and quantum information processing largely rests on our knowledge of the surrounding crystalline lattice, often obscured by the countless classes of point defects the material can host. Here we monitor the fluorescence from a negatively charged nitrogen-vacancy (NV-) center in diamond as we illuminate its vicinity. Cyclic charge state conversion of neighboring point defects sensitive to the excitation beam leads to a position-dependent stream of photo-generated carriers whose capture by the probe NV- leads to a fluorescence change. This "charge-to-photon" conversion scheme allows us to image other individual point defects surrounding the probe NV, including non-fluorescent "single-charge emitters" that would otherwise remain unnoticed. Given the ubiquity of color center photo-chromism, this strategy may likely find extensions to material systems other than diamond.


title: Theoretical description of fission yields: towards a fast and efficient global model

abstract: Background: A quantitative microscopic understanding of the fission-fragment yield distributions represents a major challenge for nuclear theory as it involves the intricate competition between large-amplitude nuclear collective motion and single-particle nucleonic motion. Purpose: A recently proposed approach to global modeling of fission fragment distributions is extended to account for odd-even staggering in charge yields and for neutron evaporation. Method: Fission trajectories are obtained within the density functional theory framework, allowing for a microscopic determination of the most probable fission prefragment configurations. Mass and charge yields distributions are constructed by means of a statistical approach rooted in a microcanonical ensemble. Result: We show that the proposed hybrid model can reproduce experimental mass and charge fragment yields, including the odd-even staggering, for a wide range of fissioning nuclei. Experimental isotopic yields can be described within a simple neutron evaporation scheme. We also explore fission fragment distributions of exotic neutron-rich and superheavy systems, and compare our predictions with other state-of-the art global calculations. Conclusion: Our study suggests that the microscopic rearrangement of nucleons into fission fragments occurs well before the scission, and that the subsequent dynamics is mainly driven by the thermal excitations and bulk features of the nuclear binding. The proposed simple hybrid approach is well suited for large-scale calculations involving hundreds of fissioning nuclei.


title: Disentangling the neighbouring pulsars SXP 15.3 and SXP 305

abstract: SXP 15.3 and SXP 305 are two Be X-ray binaries in the Small Magellanic Cloud that are spatially separated by ~7 arcsec. The small separation between these sources has, in the past, resulted in confusion about the origin of the emission from the combined region. We present long-term optical and X-ray monitoring results of both sources, where we study the historic and recent behaviour. In particular, from data collected as part of the S-CUBED project we see repeating X-ray outbursts from the combined region of the two sources in the recent lightcurve from the Neil Gehrels Swift Observatory, and we investigate the origin of this emission. Using the H-alpha emission line from the Southern African Large Telescope (SALT) and photometric flux from the Optical Gravitational Lensing Experiment (OGLE) to study the changes in the size and structure of the Be disc, we demonstrate that the X-ray emission likely originates from SXP 15.3. Timing analysis reveals unusual behaviour, where the optical outburst profile shows modulation at twice the frequency of the X-ray outbursts. We consider either of these periodicities being the true orbital period in SXP 15.3 and propose models based on the geometric orientations of the Be disc and neutron star to explain the physical origin of the outbursts.


title: Integrated Sensing and Communication with Millimeter Wave Full Duplex Hybrid Beamforming

abstract: Integrated Sensing and Communication (ISAC) has attracted substantial attraction in recent years for spectral efficiency improvement, enabling hardware and spectrum sharing for simultaneous sensing and signaling operations. In-band Full Duplex (FD) is being considered as a key enabling technology for ISAC applications due to its simultaneous transmission and reception capability. In this paper, we present an FD-based ISAC system operating at millimeter Wave (mmWave) frequencies, where a massive Multiple-Input Multiple-Output (MIMO) Base Station (BS) node employing hybrid Analog and Digital (A/D) beamforming is communicating with a DownLink (DL) multi-antenna user and the same waveform is utilized at the BS receiver for sensing the radar targets in its coverage environment. We develop a sensing algorithm that is capable of estimating Direction of Arrival (DoA), range, and relative velocity of the radar targets. A joint optimization framework for designing the A/D transmit and receive beamformers as well as the Self-Interference (SI) cancellation is presented with the objective to maximize the achievable DL rate and the accuracy of the radar target sensing performance. Our simulation results, considering fifth Generation (5G) Orthogonal Frequency Division Multiplexing (OFDM) waveforms, verify our approach's high precision in estimating DoA, range, and velocity of multiple radar targets, while maximizing the DL communication rate.


title: Real-space visualization of short-range antiferromagnetic correlations in a magnetically enhanced thermoelectric

abstract: Short-range magnetic correlations can significantly increase the thermopower of magnetic semiconductors, representing a noteworthy development in the decades-long effort to develop high-performance thermoelectric materials. Here, we reveal the nature of the thermopower-enhancing magnetic correlations in the antiferromagnetic semiconductor MnTe. Using magnetic pair distribution function analysis of neutron scattering data, we obtain a detailed, real-space view of robust, nanometer-scale, antiferromagnetic correlations that persist into the paramagnetic phase above the N\'eel temperature $T_{\mathrm{N}}$ = 307 K. The magnetic correlation length in the paramagnetic state is significantly longer along the crystallographic $c$ axis than within the $ab$ plane, pointing to anisotropic magnetic interactions. Ab initio calculations of the spin-spin correlations using density functional theory in the disordered local moment approach reproduce this result with quantitative accuracy. These findings constitute the first real-space picture of short-range spin correlations in a magnetically enhanced thermoelectric and inform future efforts to optimize thermoelectric performance by magnetic means.


title: Neural Circuit Architectural Priors for Embodied Control

abstract: Artificial neural networks for motor control usually adopt generic architectures like fully connected MLPs. While general, these tabula rasa architectures rely on large amounts of experience to learn, are not easily transferable to new bodies, and have internal dynamics that are difficult to interpret. In nature, animals are born with highly structured connectivity in their nervous systems shaped by evolution; this innate circuitry acts synergistically with learning mechanisms to provide inductive biases that enable most animals to function well soon after birth and learn efficiently. Convolutional networks inspired by visual circuitry have encoded useful biases for vision. However, it is unknown the extent to which ANN architectures inspired by neural circuitry can yield useful biases for other AI domains. In this work, we ask what advantages biologically inspired ANN architecture can provide in the domain of motor control. Specifically, we translate C. elegans locomotion circuits into an ANN model controlling a simulated Swimmer agent. On a locomotion task, our architecture achieves good initial performance and asymptotic performance comparable with MLPs, while dramatically improving data efficiency and requiring orders of magnitude fewer parameters. Our architecture is interpretable and transfers to new body designs. An ablation analysis shows that constrained excitation/inhibition is crucial for learning, while weight initialization contributes to good initial performance. Our work demonstrates several advantages of biologically inspired ANN architecture and encourages future work in more complex embodied control.


title: Unsteady stretching of a glass tube with internal channel pressurisation

abstract: Mathematical modelling is used to examine the unsteady problem of heating and pulling an axisymmetric cylindrical glass tube with an over-pressure applied within the tube to form tapers with a near uniform bore and small wall thickness at the tip. To allow for the dependence of viscosity on temperature, a prescribed axially varying viscosity is assumed. Our motivation is the manufacture of emitter tips for mass spectrometry which provide a continuous fluid flow and do not become blocked. We demonstrate, for the first time, the feasibility of producing such emitters by this process and examine the influence of the process parameters, in particular the pulling force and over-pressure, on the geometry. For small values of these parameters a more uniform bore may be obtained but the geometry is highly sensitive to small fluctuations. Higher values of the parameters result in more variation in the bore size but less sensitivity to fluctuations. The best parameters depend on the accuracy of the puller used to manufacture the tapers and the permissible tolerances on the geometry. The model has wider application to the manufacture of other devices.


title: Beyond chord vocabularies: Exploiting pitch-relationships in a chord estimation metric

abstract: Chord estimation metrics treat chord labels as independent of one another. This fails to represent the pitch relationships between the chords in a meaningful way, resulting in evaluations that must make compromises with complex chord vocabularies and that often require time-consuming qualitative analyses to determine details about how a chord estimation algorithm performs. This paper presents an accuracy metric for chord estimation that compares the pitch content of the estimated chords against the ground truth that captures both the correct notes that are estimated and additional notes that are inserted into the estimate. This is not a stand-alone evaluation protocol but rather a metric that can be integrated as a weighting into existing evaluation approaches.


title: Standardizing Platinum Dainotti-correlated gamma-ray bursts, and using them with standardized Amati-correlated gamma-ray bursts to constrain cosmological model parameters

abstract: We show that the Platinum gamma-ray burst (GRB) data compilation, probing the redshift range $0.553 \leq z \leq 5.0$, obeys a cosmological-model-independent three-parameter fundamental plane (Dainotti) correlation and so is standardizable. While they probe the largely unexplored $z \sim 2.3-5$ part of cosmological redshift space, the GRB cosmological parameter constraints are consistent with, but less precise than, those from a combination of baryon acoustic oscillation (BAO) and Hubble parameter [$H(z)$] data. In order to increase the precision of GRB-only cosmological constraints, we exclude common GRBs from the larger Amati-correlated A118 data set composed of 118 GRBs and jointly analyze the remaining 101 Amati-correlated GRBs with the 50 Platinum GRBs. This joint 151 GRB data set probes the largely unexplored $z \sim 2.3-8.2$ region; the resulting GRB-only cosmological constraints are more restrictive, and consistent with, but less precise than, those from $H(z)$ + BAO data.


title: Asymptotic values of entire functions of infinite order

abstract: We prove that there exists an entire function for which every complex number is an asymptotic value and whose growth is arbitrarily slow subject only to the necessary condition that the function is of infinite order.


title: Multi-agent Motion Planning from Signal Temporal Logic Specifications

abstract: We tackle the challenging problem of multi-agent cooperative motion planning for complex tasks described using signal temporal logic (STL), where robots can have nonlinear and nonholonomic dynamics. Existing methods in multi-agent motion planning, especially those based on discrete abstractions and model predictive control (MPC), suffer from limited scalability with respect to the complexity of the task, the size of the workspace, and the planning horizon. We present a method based on {\em timed waypoints\/} to address this issue. We show that timed waypoints can help abstract nonlinear behaviors of the system as safety envelopes around the reference path defined by those waypoints. Then the search for waypoints satisfying the STL specifications can be inductively encoded as a mixed-integer linear program. The agents following the synthesized timed waypoints have their tasks automatically allocated, and are guaranteed to satisfy the STL specifications while avoiding collisions. We evaluate the algorithm on a wide variety of benchmarks. Results show that it supports multi-agent planning from complex specification over long planning horizons, and significantly outperforms state-of-the-art abstraction-based and MPC-based motion planning methods. The implementation is available at https://github.com/sundw2014/STLPlanning.


title: The Magnetar Connection

abstract: We investigate the combined evolution of the dipolar surface magnetic field (B$_{s}$) and the spin-period (P$_s$) of known magnetars and high magnetic field (B$_s$ $ \gtrsim 10^{13}$~G) radio pulsars. We study the long term behaviour of these objects assuming a simple ohmic dissipation of the magnetic field. Identifying the regions (in the P$_s$-B$_s$ plane) in which these neutron stars would likely move into, before crossing the death-line to enter the pulsar graveyard, we comment upon the possible connection between the magnetars and other classes of neutron stars.


title: An Empirical Study of Least Squares Ratings for USA Ultimate Frisbee

abstract: Ultimate frisbee is one of the fastest-growing sports in the world. In the United States, the governing body USA Ultimate uses a custom power rating system to determine bid allocations for various competitive tournaments. However, this rating system has significant flaws and leaves room for improvement. In this paper, we apply the least squares rating system and demonstrate its improvement over the current system both qualitatively and relative to a number of quantitative metrics.


title: Consistent Approximations in Composite Optimization

abstract: Approximations of optimization problems arise in computational procedures and sensitivity analysis. The resulting effect on solutions can be significant, with even small approximations of components of a problem translating into large errors in the solutions. We specify conditions under which approximations are well behaved in the sense of minimizers, stationary points, and level-sets and this leads to a framework of consistent approximations. The framework is developed for a broad class of composite problems, which are neither convex nor smooth. We demonstrate the framework using examples from stochastic optimization, neural-network based machine learning, distributionally robust optimization, penalty and augmented Lagrangian methods, interior-point methods, homotopy methods, smoothing methods, extended nonlinear programming, difference-of-convex programming, and multi-objective optimization. An enhanced proximal method illustrates the algorithmic possibilities. A quantitative analysis supplements the development by furnishing rates of convergence.


title: Multi-time structure functions and the Lagrangian scaling of turbulence

abstract: We define and characterize multi-time Lagrangian structure functions using data stemming from two swirling flows with mean flow and turbulent fluctuations: A Taylor-Green numerical flow, and a von K\'arm\'an laboratory experiment. Data is obtained from numerical integration of tracers in the former case, and from three-dimensional particle tracking velocimetry measurements in the latter. Multi-time statistics are shown to decrease the contamination of large scales in the inertial range scaling. A time scale at which contamination from the mean flow becomes dominant is identified, with this scale separating two different Lagrangian scaling ranges. The results from the multi-time structure functions also indicate that Lagrangian intermittency is not a result of large-scale flow effects. The multi-time Lagrangian structure functions can be used without prior knowledge of the forcing mechanisms or boundary conditions, allowing their application in different flow geometries.


title: The NP-completeness of Redundant Open-Locating-Dominating Set

abstract: For a graph G, a dominating set D is a subset of vertices in G where each of the vertices in G is in D or adjacent to some vertex in D. An open-locating-dominating (OLD) set models a system with sensors to detect an intruder in a facility or a faulty component in a network of processors. The goal is to detect and pinpoint an intruder's exact location in a system with a minimum number of sensors. In this paper, we focus on a variant of an OLD set called a redundant OLD set and present a proof for the NP-completeness of the problem of a redundant OLD set.


title: Specifying Evacuation Return and Home-switch Stability During Short-term Disaster Recovery Using Location-based Data

abstract: The objectives of this study are: (1) to specify evacuation return and home-switch stability as two critical milestones of short-term recovery during and in the aftermath of disasters; and (2) to understand the disparities among subpopulations in the duration of these critical recovery milestones. Using privacy-preserving fine-resolution location-based data, we examine evacuation return and home move-out rates in Harris County, Texas in the context of the 2017 Hurricane Harvey. For each of the two critical recovery milestones, the results reveal the areas with short- and long-return durations and enable evaluating disparities in evacuation return and home-switch stability patterns. In fact, a shorter duration of critical recovery milestone indicators in flooded areas is not necessarily a positive indication. Shorter evacuation return could be due to barriers to evacuation and shorter home move-out rate return for lower-income residents is associated with living in rental homes. In addition, skewed and non-uniform recovery patterns for both the evacuation return and home-switch stability were observed in all subpopulation groups. All return patterns show a two-phase return progress pattern. The findings could inform disaster managers and public officials to perform recovery monitoring and resource allocation in a more proactive, data-driven, and equitable manner.


title: How Smart Should a Forager Be?

abstract: We introduce an idealized model of an intelligent forager in which higher intelligence corresponds to a larger spatial range over which the forager can detect food. Such a forager diffuses randomly whenever the nearest food is more distant than the forager's detection range, $R$, and moves ballistically towards the nearest food inside its detection range. Concomitantly, the forager's metabolic energy cost per step is an increasing function of its intelligence. A dumb forager wanders randomly and may miss nearby food, thus making it susceptible to starvation. Conversely, a too-smart forager incurs a large metabolic cost per step during its search for food and is again susceptible to starvation. We show that the forager's lifetime is maximized at an optimal, intermediate level of intelligence.


title: Unstable higher Toda brackets III

abstract: We define inductively unstable n-fold Toda brackets for every n>2 in the category of spaces with base points, and then define stable ones.


title: DapStep: Deep Assignee Prediction for Stack Trace Error rePresentation

abstract: The task of finding the best developer to fix a bug is called bug triage. Most of the existing approaches consider the bug triage task as a classification problem, however, classification is not appropriate when the sets of classes change over time (as developers often do in a project). Furthermore, to the best of our knowledge, all the existing models use textual sources of information, i.e., bug descriptions, which are not always available. In this work, we explore the applicability of existing solutions for the bug triage problem when stack traces are used as the main data source of bug reports. Additionally, we reformulate this task as a ranking problem and propose new deep learning models to solve it. The models are based on a bidirectional recurrent neural network with attention and on a convolutional neural network, with the weights of the models optimized using a ranking loss function. To improve the quality of ranking, we propose using additional information from version control system annotations. Two approaches are proposed for extracting features from annotations: manual and using an additional neural network. To evaluate our models, we collected two datasets of real-world stack traces. Our experiments show that the proposed models outperform existing models adapted to handle stack traces. To facilitate further research in this area, we publish the source code of our models and one of the collected datasets.


title: Curvature effects in charge-regulated lipid bilayers

abstract: We formulate a theory of electrostatic interactions in lipid bilayer membranes where both monolayer leaflets contain dissociable moieties that are subject to charge regulation. We specifically investigate the coupling between membrane curvature and charge regulation of a lipid bilayer vesicle using both the linear Debye-H\"uckel (DH) and the non-linear Poisson-Boltzmann (PB) theory. We find that charge regulation of an otherwise symmetric bilayer membrane can induce charge symmetry breaking, non-linear flexoelectricity and anomalous curvature dependence of free energy. The pH effects investigated go beyond the paradigm of electrostatic renormalization of the mechano-elastic properties of membranes.

categories: cond-mat.soft cond-mat.stat-mech physics.bio-ph physics.chem-ph


title: Galaxy Correlation Function and Local Density from Photometric Redshifts Using the Stochastic Order Redshift Technique (SORT)

abstract: The stochastic order redshift technique (SORT) is a simple, efficient, and robust method to improve cosmological redshift measurements. The method relies upon having a small ($\sim$10 per cent) reference sample of high-quality redshifts. Within pencil-beam-like sub-volumes surrounding each galaxy, we use the precise dN/d$z$ distribution of the reference sample to recover new redshifts and assign them one-to-one to galaxies such that the original rank order of redshifts is preserved. Preserving the rank order is motivated by the fact that random variables drawn from Gaussian probability density functions with different means but equal standard deviations satisfy stochastic ordering. The process is repeated for sub-volumes surrounding each galaxy in the survey. This results in every galaxy with an uncertain redshift being assigned multiple "recovered" redshifts from which a new redshift estimate can be determined. An earlier paper applied SORT to a mock Sloan Digital Sky Survey at $z \lesssim$ 0.2 and accurately recovered the two-point correlation function on scales $\gtrsim$4 $h^{-1}$Mpc. In this paper, we test the performance of SORT in surveys spanning the redshift range 0.75$<z<$2.25. We used two mock surveys extracted from the Small MultiDark-Planck and Bolshoi-Planck N-body simulations with dark matter haloes that were populated by the Santa Cruz semi-analytic model. We find that SORT is able to improve redshift estimates and recover distinctive large-scale features of the cosmic web. Further, it provides unbiased estimates of the redshift-space two-point correlation function $\xi(s)$ on scales $\gtrsim$2.5 $h^{-1}$Mpc, as well as local densities in regions of average or higher density. This may allow improved understanding of how galaxy properties relate to their local environments.


title: Quantifying yeast colony morphologies with feature engineering from time-lapse photography

abstract: Baker's yeast (Saccharomyces cerevisiae) is a model organism for studying the morphology that emerges at the scale of multi-cell colonies. To look at how morphology develops, we collect a dataset of time-lapse photographs of the growth of different strains of S. cerevisiae. We discuss the general statistical challenges that arise when using time-lapse photographs to extract time-dependent features. In particular, we show how texture-based feature engineering and representative clustering can be successfully applied to categorize the development of yeast colony morphology using our dataset. The local binary pattern (LBP) from image processing is used to score the surface texture of colonies. This texture score develops along a smooth trajectory during growth. The path taken depends on how the morphology emerges. A hierarchical clustering of the colonies is performed according to their texture development trajectories. The clustering method is designed for practical interpretability; it obtains the best representative colony image for any hierarchical sub-cluster.


title: Generic 1-connectivity of flag domains in Hermitian symmetric spaces

abstract: A flag domain is an open real group orbit in a complex flag manifold. It has been shown that a flag domain is either pseudoconvex or pseudoconcave. Moreover, generically 1-connected flag domains are pseudoconcave. In this study, for flag domains contained in irreducible Hermitian symmetric spaces of type AIII or CI, we determine which pseudoconcave flag domain is generically 1-connected.


title: Adaptive Transfer Learning for Plant Phenotyping

abstract: Plant phenotyping (Guo et al. 2021; Pieruschka et al. 2019) focuses on studying the diverse traits of plants related to the plants' growth. To be more specific, by accurately measuring the plant's anatomical, ontogenetical, physiological and biochemical properties, it allows identifying the crucial factors of plants' growth in different environments. One commonly used approach is to predict the plant's traits using hyperspectral reflectance (Yendrek et al. 2017; Wang et al. 2021). However, the data distributions of the hyperspectral reflectance data in plant phenotyping might vary in different environments for different plants. That is, it would be computationally expansive to learn the machine learning models separately for one plant in different environments. To solve this problem, we focus on studying the knowledge transferability of modern machine learning models in plant phenotyping. More specifically, this work aims to answer the following questions. (1) How is the performance of conventional machine learning models, e.g., partial least squares regression (PLSR), Gaussian process regression (GPR) and multi-layer perceptron (MLP), affected by the number of annotated samples for plant phenotyping? (2) Whether could the neural network based transfer learning models improve the performance of plant phenotyping? (3) Could the neural network based transfer learning be improved by using infinite-width hidden layers for plant phenotyping?


title: Genome-wide networks reveal emergence of epidemic strains of Salmonella Enteritidis

abstract: Objectives: To enhance monitoring of high-burden foodborne pathogens, there is opportunity to combine pangenome data with network analysis. Methods: Salmonella enterica subspecies Enterica serovar Enteritidis isolates were referred to the New South Wales (NSW) Enteric Reference Laboratory between August 2015 and December 2019 (1033 isolates in total), inclusive of a confirmed outbreak. All isolates underwent whole genome sequencing. Distances between genomes were quantified by in silico MLVA as well as core SNPs, which informed construction of undirected networks. Prevalence-centrality spaces were generated from the undirected networks. Components on the undirected SNP network were considered alongside a phylogenetic tree representation. Results: Outbreak isolates were identifiable as distinct components on the MLVA and SNP networks. The MLVA network based centrality/prevalence space did not delineate the outbreak, whereas the outbreak was clearly delineated in the SNP network based centrality/prevalence space. Components on the undirected SNP network showed a high concordance to the SNP clusters based on phylogenetic analysis. Conclusions: Bacterial whole genome data in network based analysis can improve the resolution of population analysis. High concordance of network components and SNP clusters is promising for rapid population analyses of foodborne Salmonella spp. due to the low overhead of network analysis.


title: Supplementing rapid Bayesian parameter estimation schemes with adaptive grids

abstract: In the age of multi-messenger astrophysics, low-latency parameter estimation of gravitational-wave signals is essential for electromagnetic follow-up observations. In this paper, we present a new edition of the Bayesian parameter estimation scheme for compact binaries known as Rapid PE. Rapid PE parallelizes parameter estimation by fixing the intrinsic parameters (such as the masses and spins of the binary) to a grid around a search-identified intrinsic point and by integrating over the extrinsic parameters via Monte Carlo sampling. The original version of Rapid PE calculated an effective Fisher matrix to fix the intrinsic parameters to an ellipsoidal grid. Here we use rectilinear gridding in intrinsic space. The use of rectilinear intrinsic grids provides the option to more easily include spin and allows for adaptive grid refinements to mitigate search biases and map out the posterior more completely without sacrificing sampling efficiency. We demonstrate that this parameter estimation method produces reliable results for simulated gravitational-wave signals from binary neutron star mergers.


title: The Structure of Flavor Mixing and Reconstruction of the Mass Matrix

abstract: The Fermion flavor structure is investigated by bilinear decomposition of the mass matrix after EW symmetry breaking, and the roles of factorized matrices in flavor mixing and mass generation are explored. It is shown that flavor mixing can be addressed as an independent issue. On a new Yukawa basis, the minimal parameterization of flavor mixing is realized containing two relative phases and two free $SO(2)_L$ rotation angles. The validity of the flavor mixing structure is checked in both the lepton and quark sectors. Under the decomposition of flavor mixing, fermion mass matrices are reconstructed under the hierarchy limit. A flat mass matrix with all elements equal to 1 arises naturally from the requirement that homology exists between up-type and down-type fermion mass matrices. Some hints of a flat matrix and flavor breaking are also discussed.


title: Photoinduced enhancement of superconductivity in the plaquette Hubbard model

abstract: Real-time dynamics techniques have proven increasingly useful in understanding strongly correlated systems both theoretically and experimentally. By employing unbiased time-resolved exact diagonalization, we study pump dynamics in the two-dimensional plaquette Hubbard model, where distinct hopping integrals $t_h$ and $t_h^\prime$ are present within and between plaquettes. In the intermediate coupling regime, a significant enhancement of $d$-wave superconductivity is observed and compared with that obtained by simple examination of expectation values with the eigenstates of the Hamiltonian. Our work provides further understanding of superconductivity in the Hubbard model, extends the description of the pairing amplitude to the frequency-anisotropy plane, and offers a promising approach for experimentally engineering emergent out-of-equilibrium states.


title: Model predictive control for robust quantum state preparation

abstract: A critical engineering challenge in quantum technology is the accurate control of quantum dynamics. Model-based methods for optimal control have been shown to be highly effective when theory and experiment closely match. Consequently, realizing high-fidelity quantum processes with model-based control requires careful device characterization. In quantum processors based on cold atoms, the Hamiltonian can be well-characterized. For superconducting qubits operating at milli-Kelvin temperatures, the Hamiltonian is not as well-characterized. Unaccounted for physics (i.e., mode discrepancy), coherent disturbances, and increased noise compromise traditional model-based control. This work introduces model predictive control (MPC) for quantum control applications. MPC is a closed-loop optimization framework that (i) inherits a natural degree of disturbance rejection by incorporating measurement feedback, (ii) utilizes finite-horizon model-based optimizations to control complex multi-input, multi-output dynamical systems under state and input constraints, and (iii) is flexible enough to develop synergistically alongside other modern control strategies. We show how MPC can be used to generate practical optimized control sequences in representative examples of quantum state preparation. Specifically, we demonstrate for a qubit, a weakly-anharmonic qubit, and a system undergoing crosstalk, that MPC can realize successful model-based control even when the model is inadequate. These examples showcase why MPC is an important addition to the quantum engineering control suite.


title: Bi-level Volt/VAR Optimization in Distribution Networks with Smart PV Inverters

abstract: Optimal Volt/VAR control (VVC) in distribution networks relies on an effective coordination between the conventional utility-owned mechanical devices and the smart residential photovoltaic (PV) inverters. Typically, a central controller carries out a periodic optimization and sends setpoints to the local controller of each device. However, instead of tracking centrally dispatched setpoints, smart PV inverters can cooperate on a much faster timescale to reach optimality within a PV inverter group. To accommodate such PV inverter groups in the VVC architecture, this paper proposes a bi-level optimization framework. The upper-level determines the setpoints of the mechanical devices to minimize the network active power losses, while the lower-level represents the coordinated actions that the inverters take for their own objectives. The interactions between these two levels are captured in the bi-level optimization, which is solved using the Karush-Kuhn-Tucker (KKT) conditions. This framework fully exploits the capabilities of the different types of voltage regulation devices and enables them to cooperatively optimize their goals. Case studies on typical distribution networks with field-recorded data demonstrate the effectiveness and advantages of the proposed approach.