title: Compositional Servoing by Recombining Demonstrations

abstract: Learning-based manipulation policies from image inputs often show weak task transfer capabilities. In contrast, visual servoing methods allow efficient task transfer in high-precision scenarios while requiring only a few demonstrations. In this work, we present a framework that formulates the visual servoing task as graph traversal. Our method not only extends the robustness of visual servoing, but also enables multitask capability based on a few task-specific demonstrations. We construct demonstration graphs by splitting existing demonstrations and recombining them. In order to traverse the demonstration graph in the inference case, we utilize a similarity function that helps select the best demonstration for a specific task. This enables us to compute the shortest path through the graph. Ultimately, we show that recombining demonstrations leads to higher task-respective success. We present extensive simulation and real-world experimental results that demonstrate the efficacy of our approach.


title: Mechanically-driven Stem Cell Separation in Tissues caused by Proliferating Daughter Cells

abstract: The homeostasis of epithelial tissue relies on a balance between the self-renewal of stem cell populations, cellular differentiation, and loss. Although this balance needs to be tightly regulated to avoid pathologies, such as tumor growth, the regulatory mechanisms, both cell-intrinsic and collective, which ensure tissue steady-state are still poorly understood. Here, we develop a computational model that incorporates basic assumptions of stem cell renewal into distinct populations and mechanical interactions between cells. We find that the model generates unexpected dynamic features: stem cells repel each other in the bulk tissue and are thus found rather isolated, as in a number of in vivo contexts. By mapping the system onto a gas of passive Brownian particles with effective repulsive interactions, that arise from the generated flows of differentiated cells, we show that we can quantitatively describe such stem cell distribution in tissues. The interaction potential between a pair of stem cells decays exponentially with a characteristic length that spans several cell sizes, corresponding to the volume of cells generated per stem cell division. Our findings may help understanding the dynamics of normal and cancerous epithelial tissues.


title: Breakup of particle-laden droplets in airflow

abstract: The atomisation of suspension containing liquid and dispersed particles is prevalent in many applications. Previous studies of droplet breakup mainly focused on homogeneous fluids, and the heterogeneous effect of particles on the breakup progress is unclear. In this study, the breakup of particle-laden droplets in airflow is investigated experimentally. Combining synchronised high-speed images from the side view and the 45$^\circ$ view, we compare the morphology of particle-laden droplets with that of homogeneous fluids in different breakup modes. The results show that the higher effective viscosity of particle-laden droplets affects the initial deformation, and the heterogeneous effect of particles appears in the later breakup stage. To evaluate the heterogeneous effect of particles quantitatively, we eliminate the effect of the higher effective viscosity of particle-laden droplets by comparing cases corresponding to the same inviscid Weber number. The quantitative comparison reveals that the heterogeneous effect of particles accelerates the fragmentation of liquid film and promotes localised rapid piercing. A correlation length that depends on the particle diameter and the volume fraction is proposed to characterise the length scale of the concentration fluctuation under the combined effect of the initial flattening and later stretching during the droplet breakup process. Based on this correlation length, the fragment size distributions are analysed, and the scaling results agree well with the experimental data.


title: Aspect of Clusters Correlation at Light Nuclei Excited State

abstract: The correlation of $\alpha\alpha$ was probed via measuring the transverse momentum $p_{T}$ and width $\delta p_{T}$ of one $\alpha$, for the first time, which represents the spatial and dynamical essentialities of the initial coupling state in $^{8}$Be nucleus. The weighted interaction vertex of 3$\alpha$ reflected by the magnitudes of their relative momentums and relative emission angles proves the isosceles triangle configuration for 3$\alpha$ at the high excited energy analogous Hoyle states.


title: Two-loop five-particle scattering amplitudes

abstract: I discuss the recent advances in the computation of two-loop scattering amplitudes for five-particle processes. The latter are fundamental ingredients to obtain predictions at the next-to-next-to-leading order (NNLO) in QCD for many interesting LHC processes. I discuss the state-of-the-art technology for computing scattering amplitudes analytically, and present new results relevant for the LHC phenomenology.


title: From task structures to world models: What do LLMs know?

abstract: In what sense does a large language model have knowledge? The answer to this question extends beyond the capabilities of a particular AI system, and challenges our assumptions about the nature of knowledge and intelligence. We answer by granting LLMs "instrumental knowledge"; knowledge defined by a certain set of abilities. We then ask how such knowledge is related to the more ordinary, "worldly" knowledge exhibited by human agents, and explore this in terms of the degree to which instrumental knowledge can be said to incorporate the structured world models of cognitive science. We discuss ways LLMs could recover degrees of worldly knowledge, and suggest such recovery will be governed by an implicit, resource-rational tradeoff between world models and task demands.


title: Measurement of the CKM angle $\gamma$ using the $B^{\pm}\rightarrow D^{*} h^{\pm}$ channels

abstract: A measurement of the $CP$-violating observables from $B^{\pm}\rightarrow D^* K^{\pm}$ and $B^{\pm}\rightarrow D^* \pi^{\pm}$ decays is presented, where $D^* (D) $ is an admixture of $D^{*0}$ and $\bar{D}^{*0}$ ($D^0$ and $\bar{D}^0$) states and is reconstructed through the decay chains $ D^* \rightarrow D\pi^0/\gamma$ and $D \to K_S^0 \pi^+\pi^-/K_S^0 K^+K^-$. The measurement is performed by analysing the signal yield variation across the $D$ decay phase space and is independent of any amplitude model. The data sample used was collected by the LHCb experiment in proton-proton collisions and corresponds to a total integrated luminosity of 9 fb$^{-1}$ at centre-of-mass energies of 7, 8 and 13 TeV. The CKM angle $\gamma$ is determined to be $(69^{+13}_{-14})^{\circ}$ using the measured $CP$-violating observables. The hadronic parameters $r^{D^* K^{\pm}}_B, r^{D^* \pi^{\pm}}_B, \delta^{D^* K^{\pm}}_B, \delta^{D^* \pi^{\pm}}_B$, which are the ratios and strong phase differences between favoured and suppressed $B^{\pm}$ decays, are also reported.


title: Spin-vibronic dynamics in open-shell systems beyond the spin Hamiltonian formalism

abstract: Vibronic coupling has a dramatic influence over a large number of molecular processes, ranging from photo-chemistry, to spin relaxation and electronic transport. The simulation of vibronic coupling with multi-reference wavefunction methods has been largely applied to organic compounds, and only early efforts are available for open-shell systems such as transition metal and lanthanide complexes. In this work, we derive a numerical strategy to differentiate the molecular electronic Hamiltonian in the context of multi-reference ab initio methods and inclusive of spin-orbit coupling effects. We then provide a formulation of open quantum system dynamics able to predict the time evolution of the electrons' density matrix under the influence of a Markovian phonon bath up to fourth-order perturbation theory. We apply our method to Co(II) and Dy(III) molecular complexes exhibiting long spin relaxation times and successfully validate our strategy against the use of an effective spin Hamiltonian. Our study shed light on the nature of vibronic coupling, the importance of electronic excited states in spin relaxation, and the need for high-level computational chemistry to quantify it.


title: Zentropy Theory for Quantitative Prediction of Emergent Behaviors through Symmetry-Breaking Configurations

abstract: Density functional theory (DFT) is the de facto approach for predicting self-consistent-field electronic structures of ground-state configurations of complex atoms, molecules, and solids and providing their property data for materials discovery and design. This capability is greatly enabled by the generalized gradient approximation for exchange-correlation interactions with an important set of exchange-correlation functionals developed by John Perdew and his collaborators in last several decades. The scientific community and the present author's group have greatly benefited from this capability. Over the years, the present author's group has integrated the energetics from DFT-based calculations both at zero K and finite temperature into thermodynamic modeling and developed methods to predict tracer diffusivity, elastic coefficients, interfacial energy, and a number of other properties related to the derivatives of free energy. One key outcome is the accurate prediction of free energy of a system through the consideration of both ground-state and stable symmetry-breaking non-ground-state configurations. It is articulated that phonon properties of all individual configurations can be accurately calculated by quasiharmonic approximations in the temperature and volume ranges of interest, and the emergent behaviors and anharmonicity of a system originate primarily from the statistical competition among all the configurations.


title: Multi-Industry Simplex : A Probabilistic Extension of GICS

abstract: Accurate industry classification is a critical tool for many asset management applications. While the current industry gold-standard GICS (Global Industry Classification Standard) has proven to be reliable and robust in many settings, it has limitations that cannot be ignored. Fundamentally, GICS is a single-industry model, in which every firm is assigned to exactly one group - regardless of how diversified that firm may be. This approach breaks down for large conglomerates like Amazon, which have risk exposure spread out across multiple sectors. We attempt to overcome these limitations by developing MIS (Multi-Industry Simplex), a probabilistic model that can flexibly assign a firm to as many industries as can be supported by the data. In particular, we utilize topic modeling, an natural language processing approach that utilizes business descriptions to extract and identify corresponding industries. Each identified industry comes with a relevance probability, allowing for high interpretability and easy auditing, circumventing the black-box nature of alternative machine learning approaches. We describe this model in detail and provide two use-cases that are relevant to asset management - thematic portfolios and nearest neighbor identification. While our approach has limitations of its own, we demonstrate the viability of probabilistic industry classification and hope to inspire future research in this field.


title: The avalanche delay effect in sine-gated single-photon detector based on InGaAs/InP SPADs

abstract: A sine-gated single-photon detector (SPD) intended for use in a quantum key distribution (QKD) system is considered in this paper. An "avalanche delay" effect in the sine-gated SPD is revealed. This effect consists in the appearance of an avalanche triggered at the next gate after the photon arrival gate. It has been determined experimentally that the nature of this effect is not related to the known effects of afterpulsing or charge persistence. This effect negatively affects the overall error rate in the QKD system. The influence of the main detector control parameters, such as temperature, gate amplitude and comparator's threshold voltage, on the avalanche delay effect was experimentally established.


title: Multi-alpha Boson Gas state in Fusion Evaporation Reaction and Three-body Force

abstract: The experimental evidence for the $\alpha$ Boson gas state in the $^{11}$C+$^{12}$C$\rightarrow$$^{23}$Mg$^{\ast}$ fusion evaporation reaction is presented. By measuring the $\alpha$ emission spectrum with multiplicity 2 and 3, we provide insight into the existence of a three-body force among $\alpha$ particles. The observed spectrum exhibited distinct tails corresponding to $\alpha$ particles emitted in pairs and triplets consistent well with the model-calculations of AV18-UX and chiral effective field theory of NV2-3-la*, indicating the formation of $\alpha$ clusters with three-body force in the Boson gas state.


title: On the Error-Propagation of Inexact Hotelling's Deflation for Principal Component Analysis

abstract: Principal Component Analysis (PCA) aims to find subspaces spanned by the so-called principal components that best represent the variance in the dataset. The deflation method is a popular meta-algorithm that sequentially finds individual principal components, starting from the most important ones and working towards the less important ones. However, as deflation proceeds, numerical errors from the imprecise estimation of principal components propagate due to its sequential nature. This paper mathematically characterizes the error propagation of the inexact Hotelling's deflation method. We consider two scenarios: $i)$ when the sub-routine for finding the leading eigenvector is abstract and can represent various algorithms; and $ii)$ when power iteration is used as the sub-routine. In the latter case, the additional directional information from power iteration allows us to obtain a tighter error bound than the sub-routine agnostic case. For both scenarios, we explicitly characterize how the errors progress and affect subsequent principal component estimations.


title: A Finite Element Method Approach for Trajectory Generation via Time-Optimal Control and Model Predictive Control Tracking

abstract: In this paper a framework for solving the time optimal control (TOC) using Galerkin's Weighted Residuals Method (GWRM) and Sequential Convex Programming (SCP) is proposed. The proposed method solves the two-point boundary value problem, avoiding the use of shooting methods that rely heavily on the appropriate initialization of the adjoint state and optimal time. Since TOC yields an open-loop controller, a Model Predictive Control (MPC) scheme is employed to track both the optimal trajectory and controller, allowing the system to reject disturbances. The approach is validated using the Dubins' car dynamics for optimal time trajectory generation.


title: Assessing Robustness via Score-Based Adversarial Image Generation

abstract: Most adversarial attacks and defenses focus on perturbations within small $\ell_p$-norm constraints. However, $\ell_p$ threat models cannot capture all relevant semantic-preserving perturbations, and hence, the scope of robustness evaluations is limited. In this work, we introduce Score-Based Adversarial Generation (ScoreAG), a novel framework that leverages the advancements in score-based generative models to generate adversarial examples beyond $\ell_p$-norm constraints, so-called unrestricted adversarial examples, overcoming their limitations. Unlike traditional methods, ScoreAG maintains the core semantics of images while generating realistic adversarial examples, either by transforming existing images or synthesizing new ones entirely from scratch. We further exploit the generative capability of ScoreAG to purify images, empirically enhancing the robustness of classifiers. Our extensive empirical evaluation demonstrates that ScoreAG matches the performance of state-of-the-art attacks and defenses across multiple benchmarks. This work highlights the importance of investigating adversarial examples bounded by semantics rather than $\ell_p$-norm constraints. ScoreAG represents an important step towards more encompassing robustness assessments.


title: Physics-constrained symbolic model discovery for polyconvex incompressible hyperelastic materials

abstract: We present a machine learning framework capable of consistently inferring mathematical expressions of hyperelastic energy functionals for incompressible materials from sparse experimental data and physical laws. To achieve this goal, we propose a polyconvex neural additive model (PNAM) that enables us to express the hyperelastic model in a learnable feature space while enforcing polyconvexity. An upshot of this feature space obtained via the PNAM is that (1) it is spanned by a set of univariate basis that can be re-parametrized with a more complex mathematical form, and (2) the resultant elasticity model is guaranteed to fulfill the polyconvexity, which ensures that the acoustic tensor remains elliptic for any deformation. To further improve the interpretability, we use genetic programming to convert each univariate basis into a compact mathematical expression. The resultant multi-variable mathematical models obtained from this proposed framework are not only more interpretable but are also proven to fulfill physical laws. By controlling the compactness of the learned symbolic form, the machine learning-generated mathematical model also requires fewer arithmetic operations than its deep neural network counterparts during deployment. This latter attribute is crucial for scaling large-scale simulations where the constitutive responses of every integration point must be updated within each incremental time step. We compare our proposed model discovery framework against other state-of-the-art alternatives to assess the robustness and efficiency of the training algorithms and examine the trade-off between interpretability, accuracy, and precision of the learned symbolic hyperelastic models obtained from different approaches. Our numerical results suggest that our approach extrapolates well outside the training data regime due to the precise incorporation of physics-based knowledge.


title: Spectroscopic variability of massive pre-main-sequence stars in M17

abstract: It is a challenge to study the formation process of massive stars: their formation time is short, they are few, often deeply embedded, and at relatively large distances. Our strategy is to study the outcome of the star formation process and to look for signatures remnant of the formation. We have access to a unique sample of (massive) pre-main-sequence (PMS) stars in the giant HII region M17, showing a photosphere and circumstellar disk. The aim is to determine the variability properties of the hot gaseous disks to understand the physical origin of the emission lines and identify dominant physical processes in these disks. We have obtained multiple-epoch (4-5 epochs) VLT/X-shooter spectra of six young stars in M17 covering about a decade. Using stacked spectra we update the spectral classification and identify circumstellar features. With the temporal variance method (TVS) we determine the extent and amplitude of the spectral line variations. The double-peaked emission lines in the PMS stars with gaseous disks are used to determine peak-to-peak velocities, V/R-ratios and the radial velocity of the systems. We identify many disk features, under which a new detection of CO bandhead and CI emission. In three of the stars we detect spectral variability, mainly in lines originating in the circumstellar disk, in a velocity range up to 320 km/s. In two PMS stars the ratio between the blue and red peaks shows a correlation with the peak-to-peak velocity, possibly explained by a spiral-arm structure in the disk. The PMS stars with variability are at similar positions in the HRD but show significant differences in disk lines and variability. The extent and timescale of the variability differs for each star and per line (sets). We find indications for an accretion flow, slow disk winds and/or disk structures in the hot gaseous inner disk as the cause of the variability in these PMS stars.


title: Searching for Optimal Runtime Assurance via Reachability and Reinforcement Learning

abstract: A runtime assurance system (RTA) for a given plant enables the exercise of an untrusted or experimental controller while assuring safety with a backup (or safety) controller. The relevant computational design problem is to create a logic that assures safety by switching to the safety controller as needed, while maximizing some performance criteria, such as the utilization of the untrusted controller. Existing RTA design strategies are well-known to be overly conservative and, in principle, can lead to safety violations. In this paper, we formulate the optimal RTA design problem and present a new approach for solving it. Our approach relies on reward shaping and reinforcement learning. It can guarantee safety and leverage machine learning technologies for scalability. We have implemented this algorithm and present experimental results comparing our approach with state-of-the-art reachability and simulation-based RTA approaches in a number of scenarios using aircraft models in 3D space with complex safety requirements. Our approach can guarantee safety while increasing utilization of the experimental controller over existing approaches.


title: A Dataset of Anatomical Environments for Medical Robots: Modeling Respiratory Deformation

abstract: Anatomical models of a medical robot's environment can significantly help guide design and development of a new robotic system. These models can be used for benchmarking motion planning algorithms, evaluating controllers, optimizing mechanical design choices, simulating procedures, and even as resources for data generation. Currently, the time-consuming task of generating these environments is repeatedly performed by individual research groups and rarely shared broadly. This not only leads to redundant efforts, but also makes it challenging to compare systems and algorithms accurately. In this work, we present a collection of clinically-relevant anatomical environments for medical robots operating in the lungs. Since anatomical deformation is a fundamental challenge for medical robots operating in the lungs, we describe a way to model respiratory deformation in these environments using patient-derived data. We share the environments and deformation data publicly by adding them to the Medical Robotics Anatomical Dataset (Med-RAD), our public dataset of anatomical environments for medical robots.


title: Model order reduction by convex displacement interpolation

abstract: We present a nonlinear interpolation technique for parametric fields that exploits optimal transportation of coherent structures of the solution to achieve accurate performance. The approach generalizes the nonlinear interpolation procedure introduced in [Iollo, Taddei, J. Comput. Phys., 2022] to multi-dimensional parameter domains and to datasets of several snapshots. Given a library of high-fidelity simulations, we rely on a scalar testing function and on a point set registration method to identify coherent structures of the solution field in the form of sorted point clouds. Given a new parameter value, we exploit a regression method to predict the new point cloud; then, we resort to a boundary-aware registration technique to define bijective mappings that deform the new point cloud into the point clouds of the neighboring elements of the dataset, while preserving the boundary of the domain; finally, we define the estimate as a weighted combination of modes obtained by composing the neighboring snapshots with the previously-built mappings. We present several numerical examples for compressible and incompressible, viscous and inviscid flows to demonstrate the accuracy of the method. Furthermore, we employ the nonlinear interpolation procedure to augment the dataset of simulations for linear-subspace projection-based model reduction: our data augmentation procedure is designed to reduce offline costs -- which are dominated by snapshot generation -- of model reduction techniques for nonlinear advection-dominated problems.


title: Exploring the topological sector optimization on quantum computers

abstract: Optimization problems are the core challenge in many fields of science and engineering, yet general and effective methods are scarce for searching optimal solutions. Quantum computing has been envisioned to help solve such problems, for example, the quantum annealing (QA) method based on adiabatic evolution has been extensively explored and successfully implemented on quantum simulators such as D-wave's annealers and some Rydberg arrays. In this work, we investigate topological sector optimization (TSO) problem, which attracts particular interests in the quantum many-body physics community. We reveal that the topology induced by frustration in the spin model is an intrinsic obstruction for QA and other traditional methods to approach the ground state. We demonstrate that the optimization difficulties of TSO problem are not restricted to the gaplessness, but are also due to the topological nature which are often ignored for the analysis of optimization problems before. To solve TSO problems, we utilize quantum imaginary time evolution (QITE) with a possible realization on quantum computers, which exploits the property of quantum superposition to explore the full Hilbert space and can thus address optimization problems of topological nature. We report the performance of different quantum optimization algorithms on TSO problems and demonstrate that their capability to address optimization problems are distinct even when considering the quantum computational resources required for practical QITE implementations.


title: Towards Foundational Models for Molecular Learning on Large-Scale Multi-Task Datasets

abstract: Recently, pre-trained foundation models have enabled significant advancements in multiple fields. In molecular machine learning, however, where datasets are often hand-curated, and hence typically small, the lack of datasets with labeled features, and codebases to manage those datasets, has hindered the development of foundation models. In this work, we present seven novel datasets categorized by size into three distinct categories: ToyMix, LargeMix and UltraLarge. These datasets push the boundaries in both the scale and the diversity of supervised labels for molecular learning. They cover nearly 100 million molecules and over 3000 sparsely defined tasks, totaling more than 13 billion individual labels of both quantum and biological nature. In comparison, our datasets contain 300 times more data points than the widely used OGB-LSC PCQM4Mv2 dataset, and 13 times more than the quantum-only QM1B dataset. In addition, to support the development of foundational models based on our proposed datasets, we present the Graphium graph machine learning library which simplifies the process of building and training molecular machine learning models for multi-task and multi-level molecular datasets. Finally, we present a range of baseline results as a starting point of multi-task and multi-level training on these datasets. Empirically, we observe that performance on low-resource biological datasets show improvement by also training on large amounts of quantum data. This indicates that there may be potential in multi-task and multi-level training of a foundation model and fine-tuning it to resource-constrained downstream tasks.


title: Measurement of ultrashort bi-photon correlation times with an integrated two-colour broadband SU(1,1)-interferometer

abstract: The bi-photon correlation time, a measure for the conditional uncertainty in the temporal arrival of two photons from a photon pair source, is a key performance identifier for many quantum spectroscopy applications, with shorter correlation times typically yielding better performance. Furthermore, it provides fundamental insight into the effects of dispersion on the bi-photon state. Here, we retrieve ultrashort bi-photon correlation times of around $100\,\mathrm{fs}$ by measuring simultaneously spectral and temporal interferograms at the output of an SU(1,1) interferometer based on an integrated broadband parametric down-conversion source in a $\mathrm{Ti:LiNbO}_3$ waveguide


title: Graph learning in robotics: a survey

abstract: Deep neural networks for graphs have emerged as a powerful tool for learning on complex non-euclidean data, which is becoming increasingly common for a variety of different applications. Yet, although their potential has been widely recognised in the machine learning community, graph learning is largely unexplored for downstream tasks such as robotics applications. To fully unlock their potential, hence, we propose a review of graph neural architectures from a robotics perspective. The paper covers the fundamentals of graph-based models, including their architecture, training procedures, and applications. It also discusses recent advancements and challenges that arise in applied settings, related for example to the integration of perception, decision-making, and control. Finally, the paper provides an extensive review of various robotic applications that benefit from learning on graph structures, such as bodies and contacts modelling, robotic manipulation, action recognition, fleet motion planning, and many more. This survey aims to provide readers with a thorough understanding of the capabilities and limitations of graph neural architectures in robotics, and to highlight potential avenues for future research.


title: Identifying Representations for Intervention Extrapolation

abstract: The premise of identifiable and causal representation learning is to improve the current representation learning paradigm in terms of generalizability or robustness. Despite recent progress in questions of identifiability, more theoretical results demonstrating concrete advantages of these methods for downstream tasks are needed. In this paper, we consider the task of intervention extrapolation: predicting how interventions affect an outcome, even when those interventions are not observed at training time, and show that identifiable representations can provide an effective solution to this task even if the interventions affect the outcome non-linearly. Our setup includes an outcome Y, observed features X, which are generated as a non-linear transformation of latent features Z, and exogenous action variables A, which influence Z. The objective of intervention extrapolation is to predict how interventions on A that lie outside the training support of A affect Y. Here, extrapolation becomes possible if the effect of A on Z is linear and the residual when regressing Z on A has full support. As Z is latent, we combine the task of intervention extrapolation with identifiable representation learning, which we call Rep4Ex: we aim to map the observed features X into a subspace that allows for non-linear extrapolation in A. We show that the hidden representation is identifiable up to an affine transformation in Z-space, which is sufficient for intervention extrapolation. The identifiability is characterized by a novel constraint describing the linearity assumption of A on Z. Based on this insight, we propose a method that enforces the linear invariance constraint and can be combined with any type of autoencoder. We validate our theoretical findings through synthetic experiments and show that our approach succeeds in predicting the effects of unseen interventions.


title: Towards Non-contact 3D Ultrasound for Wrist Imaging

abstract: Objective: The objective of this work is an attempt towards non-contact freehand 3D ultrasound imaging with minimal complexity added to the existing point of care ultrasound (POCUS) systems. Methods: This study proposes a novel approach of using a mechanical track for non-contact ultrasound (US) scanning. The approach thus restricts the probe motion to a linear plane, to simplify the acquisition and 3D reconstruction process. A pipeline for US 3D volume reconstruction employing an US research platform and a GPU-based edge device is developed. Results: The efficacy of the proposed approach is demonstrated through ex-vivo and in-vivo experiments. Conclusion: The proposed approach with the adjustable field of view capability, non-contact design, and low cost of deployment without significantly altering the existing setup would open doors for up gradation of traditional systems to a wide range of 3D US imaging applications. Significance: Ultrasound (US) imaging is a popular clinical imaging modality for the point-of-care bedside imaging, particularly of the wrist/knee in the pediatric population due to its non-invasive and radiation free nature. However, the limited views of tissue structures obtained with 2D US in such scenarios make the diagnosis challenging. To overcome this, 3D US imaging which uses 2D US images and their orientation/position to reconstruct 3D volumes was developed. The accurate position estimation of the US probe at low cost has always stood as a challenging task in 3D reconstruction. Additionally, US imaging involves contact, which causes difficulty to pediatric subjects while monitoring live fractures or open wounds. Towards overcoming these challenges, a novel framework is attempted in this work.


title: A Plug-and-Play Image Registration Network

abstract: Deformable image registration (DIR) is an active research topic in biomedical imaging. There is a growing interest in developing DIR methods based on deep learning (DL). A traditional DL approach to DIR is based on training a convolutional neural network (CNN) to estimate the registration field between two input images. While conceptually simple, this approach comes with a limitation that it exclusively relies on a pre-trained CNN without explicitly enforcing fidelity between the registered image and the reference. We present plug-and-play image registration network (PIRATE) as a new DIR method that addresses this issue by integrating an explicit data-fidelity penalty and a CNN prior. PIRATE pre-trains a CNN denoiser on the registration field and "plugs" it into an iterative method as a regularizer. We additionally present PIRATE+ that fine-tunes the CNN prior in PIRATE using deep equilibrium models (DEQ). PIRATE+ interprets the fixed-point iteration of PIRATE as a network with effectively infinite layers and then trains the resulting network end-to-end, enabling it to learn more task-specific information and boosting its performance. Our numerical results on OASIS and CANDI datasets show that our methods achieve state-of-the-art performance on DIR.


title: Unveiling SU(3) Flux Tubes At Nonzero Temperature: Electric Fields and Magnetic Currents

abstract: We report on the results of measuring the chromoelectric fields in a flux tube created by a static quark-antiquark pair in the finite-temperature SU(3) gauge theory. Below the deconfinement temperature the field behavior is similar to the zero-temperature case. Above the deconfinement temperature the field shape remains the same, but the field values drop when the distance between quark and antiquark increases, thus showing the disappearance of confining potential.


title: Convergent ADMM Plug and Play PET Image Reconstruction

abstract: In this work, we investigate hybrid PET reconstruction algorithms based on coupling a model-based variational reconstruction and the application of a separately learnt Deep Neural Network operator (DNN) in an ADMM Plug and Play framework. Following recent results in optimization, fixed point convergence of the scheme can be achieved by enforcing an additional constraint on network parameters during learning. We propose such an ADMM algorithm and show in a realistic [18F]-FDG synthetic brain exam that the proposed scheme indeed lead experimentally to convergence to a meaningful fixed point. When the proposed constraint is not enforced during learning of the DNN, the proposed ADMM algorithm was observed experimentally not to converge.


title: Quantum kernels for classifying dynamical singularities in a multiqubit system

abstract: Dynamical quantum phase transition is a critical phenomenon involving out-of-equilibrium states and broken symmetries without classical analogy. However, when finite-sized systems are analyzed, dynamical singularities of the rate function can appear, leading to a challenging physical characterization when parameters are changed. Here, we report a quantum support vector machine (QSVM) algorithm that uses quantum Kernels to classify dynamical singularities of the rate function for a multiqubit system. We illustrate our approach using $N$ long-range interacting qubits subjected to an arbitrary magnetic field, which induces a quench dynamics. Inspired by physical arguments, we introduce two different quantum Kernels, one inspired by the ground state manifold and the other based on a single state tomography. Our accuracy and adaptability results show that this quantum dynamical critical problem can be efficiently solved using physically inspiring quantum Kernels. Moreover, we extend our results for the case of time-dependent fields, quantum master equation, and when we increase the number of qubits.


title: Phase-Matched Second-Harmonic Generation from Metasurfaces Inside Multipass Cells

abstract: We demonstrate a simple and scalable approach to increase conversion efficiencies of nonlinear metasurfaces by incorporating them into multipass cells and by letting the pump beam to interact with the metasurfaces multiple times. We experimentally show that by metasurface design, the associated phase-matching criteria can be fulfilled. As a proof of principle, we achieve phase matching of second-harmonic generation (SHG) using a metasurface consisting of aluminium nanoparticles deposited on a glass substrate. The phase-matching condition is verified to be achieved by measuring superlinear dependence of the detected SHG as a function of number of passes. We measure an order of magnitude enhancement in the SHG signal when the incident pump traverses the metasurface up to 9 passes. Results are found to agree well with a simple model developed to estimate the generated SHG signals. We also discuss strategies to further scale-up the nonlinear signal generation. Our approach provides a clear pathway to enhance nonlinear optical responses of metasurface-based devices. The generic nature of our approach holds promise for diverse applications in nonlinear optics and photonics.


title: Living with a Red Dwarf: X-ray, UV, and Ca II Activity-Age Relationships of M Dwarfs

abstract: The vast majority of stars in the nearby stellar neighborhood are M dwarfs. Their low masses and luminosities result in slow rates of nuclear evolution and minimal changes to the star's observable properties, even along astronomical timescales. However, they possess relatively powerful magnetic dynamos and resulting X-ray to UV activity, compared to their bolometric luminosities. This magnetic activity does undergo an observable decline over time, making it an important potential age determinant for M dwarfs. Observing this activity is important for studying the outer atmospheres of these stars, but also for comparing the behaviors of different spectral type subsets of M dwarfs, e.g., those with partially vs. fully convective interiors. Beyond stellar astrophysics, understanding the X-ray to UV activity of M dwarfs over time is also important for studying the atmospheres and habitability of any hosted exoplanets. Earth-sized exoplanets, in particular, are more commonly found orbiting M dwarfs than any other stellar type, and thermal escape (driven by the M dwarf X-ray to UV activity) is believed to be the dominant atmospheric loss mechanism for these planets. Utilizing recently calibrated M dwarf age-rotation relationships, also constructed as part of the $\textit{Living with a Red Dwarf}$ program (Engle & Guinan 2023), we have analyzed the evolution of M dwarf activity over time, in terms of coronal (X-ray), chromospheric (Lyman-$\alpha$, and Ca II), and overall X--UV (5--1700 Angstrom) emissions. The activity-age relationships presented here will be useful for studying exoplanet habitability and atmospheric loss, but also for studying the different dynamo and outer atmospheric heating mechanisms at work in M dwarfs.


title: Kernelization for Counting Problems on Graphs: Preserving the Number of Minimum Solutions

abstract: A kernelization for a parameterized decision problem $\mathcal{Q}$ is a polynomial-time preprocessing algorithm that reduces any parameterized instance $(x,k)$ into an instance $(x',k')$ whose size is bounded by a function of $k$ alone and which has the same yes/no answer for $\mathcal{Q}$. Such preprocessing algorithms cannot exist in the context of counting problems, when the answer to be preserved is the number of solutions, since this number can be arbitrarily large compared to $k$. However, we show that for counting minimum feedback vertex sets of size at most $k$, and for counting minimum dominating sets of size at most $k$ in a planar graph, there is a polynomial-time algorithm that either outputs the answer or reduces to an instance $(G',k')$ of size polynomial in $k$ with the same number of minimum solutions. This shows that a meaningful theory of kernelization for counting problems is possible and opens the door for future developments. Our algorithms exploit that if the number of solutions exceeds $2^{\mathsf{poly}(k)}$, the size of the input is exponential in terms of $k$ so that the running time of a parameterized counting algorithm can be bounded by $\mathsf{poly}(n)$. Otherwise, we can use gadgets that slightly increase $k$ to represent choices among $2^{O(k)}$ options by only $\mathsf{poly}(k)$ vertices.


title: Coding by Design: GPT-4 empowers Agile Model Driven Development

abstract: Generating code from a natural language using Large Language Models (LLMs) such as ChatGPT, seems groundbreaking. Yet, with more extensive use, it's evident that this approach has its own limitations. The inherent ambiguity of natural language presents challenges for complex software designs. Accordingly, our research offers an Agile Model-Driven Development (MDD) approach that enhances code auto-generation using OpenAI's GPT-4. Our work emphasizes "Agility" as a significant contribution to the current MDD method, particularly when the model undergoes changes or needs deployment in a different programming language. Thus, we present a case-study showcasing a multi-agent simulation system of an Unmanned Vehicle Fleet. In the first and second layer of our approach, we constructed a textual representation of the case-study using Unified Model Language (UML) diagrams. In the next layer, we introduced two sets of constraints that minimize model ambiguity. Object Constraints Language (OCL) is applied to fine-tune the code constructions details, while FIPA ontology is used to shape communication semantics and protocols. Ultimately, leveraging GPT-4, our last layer auto-generates code in both Java and Python. The Java code is deployed within the JADE framework, while the Python code is deployed in PADE framework. Concluding our research, we engaged in a comprehensive evaluation of the generated code. From a behavioural standpoint, the auto-generated code aligned perfectly with the expected UML sequence diagram. Structurally, we compared the complexity of code derived from UML diagrams constrained solely by OCL to that influenced by both OCL and FIPA-ontology. Results indicate that ontology-constrained model produce inherently more intricate code, but it remains manageable and low-risk for further testing and maintenance.


title: INVALS: An Efficient Forward Looking Inventory Allocation System

abstract: We design an Inventory Allocation System (INVALS) that, for each item-store combination, plans the quantity to be allocated from a warehouse replenishing multiple stores using trailers, while respecting the typical supply-chain constraints. We formulate a linear objective function which when maximised computes the allocation by considering not only the immediate store needs, but also its future expected demand. Such forward-looking allocation significantly improves the labour and trailer utilisation at the warehouse. To reduce overstocking, we adapt from our objective to prioritise allocating those items in excess which are sold faster at the stores, keeping the days of supply (DOS) to a minimum. For the proposed formulation, which is an instance of Mixed Integer Linear Programming (MILP), we present a scalable algorithm using the concepts of submodularity and optimal transport theory by: (i) sequentially adding trailers to stores based on maximum incremental gain, (ii) transforming the resultant linear program (LP) instance to an instance of capacity constrained optimal transport (COT), solvable using double entropic regularization and incurring the same computational complexity as the Sinkhorn algorithm. When compared against the planning engine that does the allocation only for immediate store needs, INVALS increases on an average the labour utilization by 34.70 and item occupancy in trailers by 37.08. The DOS distribution is also skewed to the left indicating that higher demand items are allocated in excess, reducing the days they are stocked. We empirically observed that for ~ 90% of replenishment cycles, the allocation results from INVALS are identical to the globally optimal MILP solution.


title: Towards A Robust Group-level Emotion Recognition via Uncertainty-Aware Learning

abstract: Group-level emotion recognition (GER) is an inseparable part of human behavior analysis, aiming to recognize an overall emotion in a multi-person scene. However, the existing methods are devoted to combing diverse emotion cues while ignoring the inherent uncertainties under unconstrained environments, such as congestion and occlusion occurring within a group. Additionally, since only group-level labels are available, inconsistent emotion predictions among individuals in one group can confuse the network. In this paper, we propose an uncertainty-aware learning (UAL) method to extract more robust representations for GER. By explicitly modeling the uncertainty of each individual, we utilize stochastic embedding drawn from a Gaussian distribution instead of deterministic point embedding. This representation captures the probabilities of different emotions and generates diverse predictions through this stochasticity during the inference stage. Furthermore, uncertainty-sensitive scores are adaptively assigned as the fusion weights of individuals' face within each group. Moreover, we develop an image enhancement module to enhance the model's robustness against severe noise. The overall three-branch model, encompassing face, object, and scene component, is guided by a proportional-weighted fusion strategy and integrates the proposed uncertainty-aware method to produce the final group-level output. Experimental results demonstrate the effectiveness and generalization ability of our method across three widely used databases.


title: Mean left-right eigenvector self-overlap in the real Ginibre ensemble

abstract: We study analytically the Chalker-Mehlig mean diagonal overlap $\mathcal{O}(z)$ between left and right eigenvectors associated with a complex eigenvalue $z$ of $N\times N$ matrices in the real Ginibre ensemble (GinOE). We first derive a general finite $N$ expression for the mean overlap and then investigate several scaling regimes in the limit $N\rightarrow \infty$. While in the generic spectral bulk and edge of the GinOE the limiting expressions for $\mathcal{O}(z)$ are found to coincide with the known results for the complex Ginibre ensemble (GinUE), in the region of eigenvalue depletion close to the real axis the asymptotic for the GinOE is considerably different. We also study numerically the distribution of diagonal overlaps and conjecture that it is the same in the bulk and at the edge of both the GinOE and GinUE, but essentially different in the depletion region of the GinOE.


title: On the regularity of solutions of some linear parabolic path-dependent PDEs

abstract: We study a class of linear parabolic path-dependent PDEs (PPDEs) defined on the space of c\`adl\`ag paths $x \in D([0,T])$, in which the coefficient functions at time $t$ depend on $x(t)$ and $\int_{0}^{t}x(s)dA_{s}$, for some (deterministic) continuous function $A$ with bounded variations. Under uniform ellipticity and H\"older regularity conditions on the coefficients, together with some technical conditions on $A$, we obtain the existence of a smooth solution to the PPDE by appealing to the notion of Dupire's derivatives. It provides a generalization to the existing literature studying the case where $A_t = t$, and complements our recent work, Bouchard and Tan (2021), on the regularity of approximate viscosity solutions for parabolic PPDEs. As a by-product, we also obtain existence and uniqueness of weak solutions for a class of path-dependent SDEs.


title: Smith-Gysin Sequence

abstract: Starting with a manifold $M$ and a semi-free action of $S^3$ on it, we have the Smith-Gysin sequence: $$ \cdots \to H^{*}( M) \to H^{*-3}(M/S^3, M^{S^3}) \oplus H^{*} (M^{S^3}) \to H^{*+1}(M/S^3, M^{S^3}) \to H^{*+1}(M) \to \cdots $$ In this paper, we construct a Smith-Gysin sequence that does not require the semi-free condition. This sequence includes a new term, referred to as the "exotic term," which depends on the subset $M^{S^1}$: $$ \cdots \to H^{*}(M) \to H^{*-3} (M/S^3, \Sigma/S^3) \oplus H^{*}(M^{S^3}) \oplus \left( H^{*-2}(M^{S^1})\right)^{-\mathbb{Z}_2} \to H^{*+1}(M/S^3,M^{S^3}) \to H^{*+1}(M) \to \cdots $$ Here, $\Sigma \subset M$ is the subset of points in $M$ whose isotropy groups are infinite. The group $\mathbb{Z}_2$ acts on $M^{S^1}$ by $j \in S^3$.


title: Spreading of information on a network: a quantum view

abstract: This paper concerns with the modeling of the spreading of information through a complex, multi-layered network, where the information is transferred from an initial transmitter to a final receiver. The mathematical model is deduced within the framework of operatorial methods, according to the formal mathematical apparatus typical of quantum mechanics. Two different approaches are considered: one based on the ($H,\rho$)-induced dynamics, and one on the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) equation. For each method, numerical results are presented.


title: Distributed Deep Joint Source-Channel Coding with Decoder-Only Side Information

abstract: We consider low-latency image transmission over a noisy wireless channel when correlated side information is present only at the receiver side (the Wyner-Ziv scenario). In particular, we are interested in developing practical schemes using a data-driven joint source-channel coding (JSCC) approach, which has been previously shown to outperform conventional separation-based approaches in the practical finite blocklength regimes, and to provide graceful degradation with channel quality. We propose a novel neural network architecture that incorporates the decoder-only side information at multiple stages at the receiver side. Our results demonstrate that the proposed method succeeds in integrating the side information, yielding improved performance at all channel conditions in terms of the various quality measures considered here, especially at low channel signal-to-noise ratios (SNRs) and small bandwidth ratios (BRs). We have made the source code of the proposed method public to enable further research, and the reproducibility of the results.


title: Bespoke scapegoats: Scientific advisory bodies and blame avoidance in the Covid-19 pandemic and beyond

abstract: Scholars have not asked why so many governments created ad hoc scientific advisory bodies (ahSABs) to address the Covid-19 pandemic instead of relying on existing public health infrastructure. We address this neglected question with an exploratory study of the US, UK, Sweden, Italy, Poland, and Uganda. Drawing on our case studies and the blame-avoidance literature, we find that ahSABs are created to excuse unpopular policies and take the blame should things go wrong. Thus, membership typically represents a narrow range of perspectives. An ahSAB is a good scapegoat because it does little to reduce government discretion and has limited ability to deflect blame back to government. Our explanation of our deviant case of Sweden, that did not create and ahSAB, reinforces our general principles. We draw the policy inference that ahSAB membership should be vetted by the legislature to ensure broad membership.


title: KoMultiText: Large-Scale Korean Text Dataset for Classifying Biased Speech in Real-World Online Services

abstract: With the growth of online services, the need for advanced text classification algorithms, such as sentiment analysis and biased text detection, has become increasingly evident. The anonymous nature of online services often leads to the presence of biased and harmful language, posing challenges to maintaining the health of online communities. This phenomenon is especially relevant in South Korea, where large-scale hate speech detection algorithms have not yet been broadly explored. In this paper, we introduce "KoMultiText", a new comprehensive, large-scale dataset collected from a well-known South Korean SNS platform. Our proposed dataset provides annotations including (1) Preferences, (2) Profanities, and (3) Nine types of Bias for the text samples, enabling multi-task learning for simultaneous classification of user-generated texts. Leveraging state-of-the-art BERT-based language models, our approach surpasses human-level accuracy across diverse classification tasks, as measured by various metrics. Beyond academic contributions, our work can provide practical solutions for real-world hate speech and bias mitigation, contributing directly to the improvement of online community health. Our work provides a robust foundation for future research aiming to improve the quality of online discourse and foster societal well-being. All source codes and datasets are publicly accessible at https://github.com/Dasol-Choi/KoMultiText.


title: Latent Graph Inference with Limited Supervision

abstract: Latent graph inference (LGI) aims to jointly learn the underlying graph structure and node representations from data features. However, existing LGI methods commonly suffer from the issue of supervision starvation, where massive edge weights are learned without semantic supervision and do not contribute to the training loss. Consequently, these supervision-starved weights, which may determine the predictions of testing samples, cannot be semantically optimal, resulting in poor generalization. In this paper, we observe that this issue is actually caused by the graph sparsification operation, which severely destroys the important connections established between pivotal nodes and labeled ones. To address this, we propose to restore the corrupted affinities and replenish the missed supervision for better LGI. The key challenge then lies in identifying the critical nodes and recovering the corrupted affinities. We begin by defining the pivotal nodes as $k$-hop starved nodes, which can be identified based on a given adjacency matrix. Considering the high computational burden, we further present a more efficient alternative inspired by CUR matrix decomposition. Subsequently, we eliminate the starved nodes by reconstructing the destroyed connections. Extensive experiments on representative benchmarks demonstrate that reducing the starved nodes consistently improves the performance of state-of-the-art LGI methods, especially under extremely limited supervision (6.12% improvement on Pubmed with a labeling rate of only 0.3%).


title: Fostering Enterprise Conversations Around Data on Collaboration Platforms

abstract: In enterprise organizations, data-driven decision making processes include the use of business intelligence dashboards and collaborative deliberation on communication platforms such as Slack or Teams. However, apart from those in data analyst roles, there is shallow engagement with dashboard content due to insufficient guidance, context, or access. Through a co-design study with nine enterprise professionals who use dashboard content to communicate with their colleagues, we identified design requirements for sharing selections from dashboards as interactive snapshots on collaboration platforms. We then developed Philo, an interactive demonstration environment centered around the template-based retargeting of dashboard content. Using Philo as a design probe, we interviewed our co-design participants and six additional data professionals, ultimately arriving at a set of design guidelines for fostering conversations around data in enterprise settings.


title: Pairing and Pair Breaking by Gauge Fluctuations in Bilayer Composite Fermion Metals

abstract: We study interlayer pairing of composite fermions in the total $\nu=1/2+1/2$ quantum Hall bilayer as a possible framework for understanding the experimentally observed transition from a compressible state at large layer spacing to a bilayer quantum Hall state at small layer spacing. We consider a model in which the effective interlayer composite fermion pairing interaction mediated by the Chern-Simons gauge fields in the two layers is singular with both attractive (out-of-phase) and repulsive (in-phase) components diverging at low frequency. If only the more singular attractive interaction is included the pairing gap obtained by solving the gap equation is proportional to the inverse of the layer spacing squared. In the so-called local approximation, we find that when the less singular repulsive interactions are also included the pairing gap still falls off as inverse layer spacing squared, consistent with recent analyses, but is strongly suppressed to a degree that may account for the fact that this predicted inverse square dependence is not observed experimentally. The analytically obtained local approximation solutions are then used as a starting point to numerically iterate the full gap equation to assess the validity of the approximation in this limit.


title: A perspective on neuroscience data standardization with Neurodata Without Borders

abstract: Neuroscience research has evolved to generate increasingly large and complex experimental data sets, and advanced data science tools are taking on central roles in neuroscience research. Neurodata Without Borders (NWB), a standard language for neurophysiology data, has recently emerged as a powerful solution for data management, analysis, and sharing. We here discuss our efforts to implement NWB data science pipelines. We describe general principles and specific use cases that illustrate successes, challenges, and non-trivial decisions in software engineering. We hope that our experience can provide guidance for the neuroscience community and help bridge the gap between experimental neuroscience and data science.


title: Model Order Reduction for the 1D Boltzmann-BGK Equation: Identifying Intrinsic Variables Using Neural Networks

abstract: Kinetic equations are crucial for modeling non-equilibrium phenomena, but their computational complexity is a challenge. This paper presents a data-driven approach using reduced order models (ROM) to efficiently model non-equilibrium flows in kinetic equations by comparing two ROM approaches: Proper Orthogonal Decomposition (POD) and autoencoder neural networks (AE). While AE initially demonstrate higher accuracy, POD's precision improves as more modes are considered. Notably, our work recognizes that the classical POD-MOR approach, although capable of accurately representing the non-linear solution manifold of the kinetic equation, may not provide a parsimonious model of the data due to the inherently non-linear nature of the data manifold. We demonstrate how AEs are used in finding the intrinsic dimension of a system and to allow correlating the intrinsic quantities with macroscopic quantities that have a physical interpretation.


title: Towards a quantum field theory description of nonlocal spacetime defects

abstract: We propose an ansatz for encoding the physics of nonlocal spacetime defects in the Green's functions for a scalar field theory defined on a causal set. This allows us to numerically study the effects of nonlocal spacetime defects on the discrete Feynman propagator of the theory defined on the causal set in 1+1 dimensions, and to compare to the defect-free limit. The latter approaches the expected continuum result, on average, when the number of points becomes large. When defects are present, two points with the same invariant spacetime interval can have different propagation amplitudes, depending on whether the propagation is between two ordinary spacetime points, two defects, or a defect and an ordinary point. We show that a coarse-grained description that is only sensitive to the average effect of the defects can be interpreted as a defect-induced mass and wave-function renormalization of the scalar theory.


title: Molecular Structural Dynamics in Water-Ethanol Mixtures: Spectroscopy with Polarized Neutrons Simultaneously Accessing Collective and Self-Diffusion

abstract: Binary mixtures of water with lower alcohols display non-linear phase behaviour upon mixing which are attributed to potential cluster formation at molecular level. Unravelling such elusive structures requires the investigation of hydrogen-bonding sub-nanosecond dynamics. We employ high-resolution neutron time-of-flight spectroscopy with polarization analysis in combination with selective deuteration to study the concentration-dependent structural dynamics, in the water rich part of the phase diagram of water-ethanol mixtures. This method enables the simultaneous access to atomic correlations in space and time, and allows us to separate spatially incoherent scattering probing self-diffusion of the ethanol fraction from the coherent scattering probing collective diffusion of the water network as a whole. Our observations indicate an enhanced rigidity of the hydrogen bond network at mesoscopic lengthscale compared to the intra-molecular scale as the ethanol fraction increases, which is consistent with the hypothesis of clusters.


title: The Value of Ancillary Services for Electrolyzers

abstract: While renewable hydrogen production is envisioned as the main business for electrolyzers, the provision of frequency-supporting ancillary services to the power grid can be seen as their potential by-products. This can generate additional revenue streams for electrolyzers, accelerating their integration by making them more financially attractive. This paper conducts a value analysis for an electrolyzer in Denmark. We develop a mixed-integer linear optimization model, wherein the physical characteristics of the electrolyzer are accurately captured. This model enables the electrolyzer to make informed bidding decisions in the day-ahead electricity as well as the frequency containment reserve (FCR) and the manual frequency restoration reserve (mFRR) markets. For modeling the activation of mFRR services in the bidding stage, we suggest pragmatic approaches. By carrying out an ex-post analysis based on historical activation, we provide lower and upper bounds for the profit of the electrolyzer, such that the true profit could lie in between depending how activation has been modeled in the bidding stage. We find that an electrolyzer in the Danish bidding zone DK1 could have earned 57% more in 2021 by proving ancillary services, while this profit increase in both bidding zones DK1 and DK2 for 2022 could have been drastic.


title: Three-dimensional characterization of the steel-concrete interface by FIB-SEM nanotomography

abstract: While it is widely accepted that the steel-concrete interface (SCI) plays an important role in governing the long-term durability of reinforced concrete structures, understanding about the primary features of the SCI that influence corrosion degradation mechanisms has remained elusive. This lack of knowledge can be attributed, on the one hand, to the complex heterogeneous nature of the SCI, and, on the other hand, the absence of experimental techniques suitable for studying the relevant features of the SCI. Here, we use focused ion beam - scanning electron microscopy (FIB-SEM) nanotomography to obtain high resolution 3D tomograms of the steel-concrete interfacial zone. Five tomograms, spanning volumes ranging from 8,000 to 200,000 cubic micrometer, were acquired for situations representative of both non-corroded and corroded SCIs. The achieved voxel size falls within the range of 30-50 nm, thus providing a resolution clearly surpassing the capabilities of computed X-ray tomography. This resolution enables the 3D characterization of the microstructure at the capillary scale, which is the scale at which relevant corrosion and related mass transport processes occur. Thus, FIB-SEM nanotomography is capable of yielding datasets of the SCI that serve as basis for the generation of digital twins of the interfacial microstructure, thereby enabling future studies about durability and corrosion of reinforced concrete at the pore scale.


title: Adjustable Robust Reinforcement Learning for Online 3D Bin Packing

abstract: Designing effective policies for the online 3D bin packing problem (3D-BPP) has been a long-standing challenge, primarily due to the unpredictable nature of incoming box sequences and stringent physical constraints. While current deep reinforcement learning (DRL) methods for online 3D-BPP have shown promising results in optimizing average performance over an underlying box sequence distribution, they often fail in real-world settings where some worst-case scenarios can materialize. Standard robust DRL algorithms tend to overly prioritize optimizing the worst-case performance at the expense of performance under normal problem instance distribution. To address these issues, we first introduce a permutation-based attacker to investigate the practical robustness of both DRL-based and heuristic methods proposed for solving online 3D-BPP. Then, we propose an adjustable robust reinforcement learning (AR2L) framework that allows efficient adjustment of robustness weights to achieve the desired balance of the policy's performance in average and worst-case environments. Specifically, we formulate the objective function as a weighted sum of expected and worst-case returns, and derive the lower performance bound by relating to the return under a mixture dynamics. To realize this lower bound, we adopt an iterative procedure that searches for the associated mixture dynamics and improves the corresponding policy. We integrate this procedure into two popular robust adversarial algorithms to develop the exact and approximate AR2L algorithms. Experiments demonstrate that AR2L is versatile in the sense that it improves policy robustness while maintaining an acceptable level of performance for the nominal case.


title: Discriminating interacting dark energy models using Statefinder diagnostic

abstract: In the present work, we perform a comparative study of different interacting dark energy (DE) models using the Statefinder diagnostics. In particular, 17 different forms of the energy transfer rate $Q$ between DE and dark matter (DM) were focused on, belonging to the following categories: i) linear models in energy densities of DE and DM, ii) non-linear models, iii) models with a change of direction of energy transfer between DE and DM, iv) models involving derivatives of the energy densities, v) parametrized interactions through a function of the coincidence parameter $\tilde{r}$, and finally we also consider vi) two kinds of models with a self-interaction between DM, without DE. These models have been already studied in the literature and constrained with observational data available at that time. In order to discriminate between them at background level, we use the Statefinder diagnostic, based on the computation and study of the so-called Statefinder parameters $r$, $s$ in addition to the deceleration parameter $q$. We plot the evolution trajectories for the several interacting models on the $r-q$, $r-s$ planes, and we find some distinctive features and departures from $\Lambda$CDM and other DE models, as Quintessence, Chaplygin Gas, running vacuum models (RVM) and Galileon.


title: Constrain relations for superfluid weight and pairings in a chiral flat band superconductor

abstract: Within ten years, flat band (FB) superconductivity has gained a huge interest for its remarkable features and connection to quantum geometry. We investigate the superconductivity in a FB system whose orbitals are inequivalent and in which the gap and the quantum metric are tunable. The key feature of the present theoretical study is to show a unique and simple constrain relation that pairings obey. Furthermore, pairings and superfluid weight in partially filled FB are shown to be controlled by those of the half-filled lattice. We argue that the geometry of the lattice or the complexity of the hopping terms have no impact on the features revealed in this work as far as the system is bipartite.


title: Variability and Predictability of a reduced-order land atmosphere coupled model

abstract: This study delves into the predictability of atmospheric blocking, zonal, and transition patterns utilizing a simplified coupled model. Initially, we comprehensively scrutinize the model's responses to environmental parameters like solar radiation, surface friction, and atmosphere-ground heat exchange. Subsequently, employing Gaussian mixture clustering, we successfully delineate distinct blocking, zonal, and transition flow regimes, unveiling their dependencies on surface friction. To gauge predictability and persistence, we compute the averaged local Lyapunov exponents for each regime. Our investigation uncovers the presence of zonal, blocking, and transition regimes, particularly under conditions of reduced surface friction. As surface friction increases further, the system transitions to a state characterized by two blocking regimes and a transition regime. Intriguingly, periodic behavior emerges under specific surface friction values, returning to patterns observed under low friction coefficients. Model resolution increase impacts the system in a way that only two regimes are then obtained with the clustering: the transition phase disappears and the predictability drops to roughly 2 days for both of the remaining regimes. In accordance with previous research findings, our study underscores that when all three regimes coexist, zonal patterns exhibit a more extended predictability horizon compared to blocking patterns. Remarkably, transition patterns exhibit reduced predictability when coexisting with the other regimes. In addition, within a specified range of surface friction values where two blocking regimes are found, it is observed that blocked atmospheric situations in the west of the applied topography are marked by instabilities and reduced predictability in contrast to the blockings appearing on the eastern side of the topography.


title: Program Synthesis with Best-First Bottom-Up Search

abstract: Cost-guided bottom-up search (BUS) algorithms use a cost function to guide the search to solve program synthesis tasks. In this paper, we show that current state-of-the-art cost-guided BUS algorithms suffer from a common problem: they can lose useful information given by the model and fail to perform the search in a best-first order according to a cost function. We introduce a novel best-first bottom-up search algorithm, which we call Bee Search, that does not suffer information loss and is able to perform cost-guided bottom-up synthesis in a best-first manner. Importantly, Bee Search performs best-first search with respect to the generation of programs, i.e., it does not even create in memory programs that are more expensive than the solution program. It attains best-first ordering with respect to generation by performing a search in an abstract space of program costs. We also introduce a new cost function that better uses the information provided by an existing cost model. Empirical results on string manipulation and bit-vector tasks show that Bee Search can outperform existing cost-guided BUS approaches when employing more complex domain-specific languages (DSLs); Bee Search and previous approaches perform equally well with simpler DSLs. Furthermore, our new cost function with Bee Search outperforms previous cost functions on string manipulation tasks.


title: Robust Losses for Decision-Focused Learning

abstract: Optimization models used to make discrete decisions often contain uncertain parameters that are context-dependent and are estimated through prediction. To account for the quality of the decision made based on the prediction, decision-focused learning (end-to-end predict-then-optimize) aims at training the predictive model to minimize regret, i.e., the loss incurred by making a suboptimal decision. Despite the challenge of this loss function being possibly non-convex and in general non-differentiable, effective gradient-based learning approaches have been proposed to minimize the expected loss, using the empirical loss as a surrogate. However, empirical regret can be an ineffective surrogate because the uncertainty in the optimization model makes the empirical regret unequal to the expected regret in expectation. To illustrate the impact of this inequality, we evaluate the effect of aleatoric and epistemic uncertainty on the accuracy of empirical regret as a surrogate. Next, we propose three robust loss functions that more closely approximate expected regret. Experimental results show that training two state-of-the-art decision-focused learning approaches using robust regret losses improves test-sample empirical regret in general while keeping computational time equivalent relative to the number of training epochs.


title: Pika: Empowering Non-Programmers to Author Executable Governance Policies in Online Communities

abstract: Internet users have formed a wide array of online communities with nuanced and diverse community goals and norms. However, most online platforms only offer a limited set of governance models in their software infrastructure and leave little room for customization. Consequently, technical proficiency becomes a prerequisite for online communities to build governance policies in code, excluding non-programmers from participation in designing community governance. In this paper, we present Pika, a system that empowers non-programmers to author a wide range of executable governance policies. At its core, Pika incorporates a declarative language that decomposes governance policies into modular components, thereby facilitating expressive policy authoring through a user-friendly, form-based web interface. Our user studies with 17 participants show that Pika can empower non-programmers to author governance policies approximately 2.5 times faster than programmers who author in code. We also provide insights about Pika's expressivity in supporting diverse policies that online communities want.


title: Lagrangian modeling of a non-homogeneous turbulent shear flow: Molding homogeneous and isotropic trajectories into a jet

abstract: Turbulence is prevalent in nature and industry, from large-scale wave dynamics to small-scale combustion nozzle sprays. In addition to the multi-scale nonlinear complexity and both randomness and coherent structures in its dynamics, practical turbulence is often non-homogeneous and anisotropic, leading to great modeling challenges. In this letter, an efficient model is proposed to predict turbulent jet statistics with high accuracy. The model leverages detailed knowledge of readily available velocity signals from idealized homogeneous turbulence and transforms them into Lagrangian trajectories of a turbulent jet. The resulting spatio-temporal statistics are compared against experimental jet data showing remarkable agreement at all scales. In particular the intermittency phenomenon is accurately mapped by the model to this inhomogeneous situation, as observed by higher-order moments and velocity increment probability density functions. Crucial to the advancement of turbulence modeling, the transformation is simple to implement, with possible extensions to other inhomogeneous flows such as wind turbine wakes and canopy flows, to name a few.


title: Game, Set, and Conflict: Evaluating Conflict and Game Frames in Indian Election News Coverage

abstract: News frames refer to how journalists organize and present information to convey a particular message or perspective to their readers. When covering elections, these frames shape how the public perceives electoral issues and events. This study examines how news frames, especially conflict and game frames, were employed by news organizations in India to cover the 2014 and 2019 general elections. We analyzed how the frames varied temporally, by region, and by the party being featured in the articles. Key findings include (i) conflict and games frames are employed more often in highly electorally consequential states (higher legislative seats) than in other states (ii) articles featuring challenger parties are more likely to have conflict and game frame articles than those featuring incumbent parties (iii) the national parties (BJP, Bharatiya Janata Party) and (INC, Indian National Congress) disproportionately feature in articles having conflict frames. Overall, our analysis highlights the evolving nature of election campaigns and how conflict and game frames play a major part in them.


title: On the Parameterized Complexity of Multiway Near-Separator

abstract: We study a new graph separation problem called Multiway Near-Separator. Given an undirected graph $G$, integer $k$, and terminal set $T \subseteq V(G)$, it asks whether there is a vertex set $S \subseteq V(G) \setminus T$ of size at most $k$ such that in graph $G-S$, no pair of distinct terminals can be connected by two pairwise internally vertex-disjoint paths. Hence each terminal pair can be separated in $G-S$ by removing at most one vertex. The problem is therefore a generalization of (Node) Multiway Cut, which asks for a vertex set for which each terminal is in a different component of $G-S$. We develop a fixed-parameter tractable algorithm for Multiway Near-Separator running in time $2^{O(k \log k)} * n^{O(1)}$. Our algorithm is based on a new pushing lemma for solutions with respect to important separators, along with two problem-specific ingredients. The first is a polynomial-time subroutine to reduce the number of terminals in the instance to a polynomial in the solution size $k$ plus the size of a given suboptimal solution. The second is a polynomial-time algorithm that, given a graph $G$ and terminal set $T \subseteq V(G)$ along with a single vertex $x \in V(G)$ that forms a multiway near-separator, computes a 14-approximation for the problem of finding a multiway near-separator not containing $x$.


title: Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations

abstract: While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed-phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute-solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.


title: Saliency-Guided Hidden Associative Replay for Continual Learning

abstract: Continual Learning is a burgeoning domain in next-generation AI, focusing on training neural networks over a sequence of tasks akin to human learning. While CL provides an edge over traditional supervised learning, its central challenge remains to counteract catastrophic forgetting and ensure the retention of prior tasks during subsequent learning. Amongst various strategies to tackle this, replay based methods have emerged as preeminent, echoing biological memory mechanisms. However, these methods are memory intensive, often preserving entire data samples, an approach inconsistent with humans selective memory retention of salient experiences. While some recent works have explored the storage of only significant portions of data in episodic memory, the inherent nature of partial data necessitates innovative retrieval mechanisms. Current solutions, like inpainting, approximate full data reconstruction from partial cues, a method that diverges from genuine human memory processes. Addressing these nuances, this paper presents the Saliency Guided Hidden Associative Replay for Continual Learning. This novel framework synergizes associative memory with replay-based strategies. SHARC primarily archives salient data segments via sparse memory encoding. Importantly, by harnessing associative memory paradigms, it introduces a content focused memory retrieval mechanism, promising swift and near-perfect recall, bringing CL a step closer to authentic human memory processes. Extensive experimental results demonstrate the effectiveness of our proposed method for various continual learning tasks.


title: A detailed study of the polarisation convention of the Giant Metrewave Radio Telescope

abstract: Our work aims to investigate the polarisation convention of the Giant Metrewave Radio Telescope (GMRT) radio telescope and understand whether the telescope follows the standard IAU/IEEE convention. The GMRT antennas are prime focus antennas, i.e. the radiation falling on the antenna feed reverses its circular polarisation. If this reflection is not taken into account, it will result in a reversal of circular polarisation. We carried out several tests to understand the GMRT polarisation convention. The observations were carried out on several strong and highly polarised pulsars with known polarisation properties at GMRT wavelengths, mainly covering frequency bands 2, 3 and 4. In addition, we tracked the signal from the feed to the fibre optic system, and fibre optic system to the user end. We also tracked satellites of known polarisations as well as utilised right and left circularly polarised helical antennas to study the polarisation convention. Our study shows that the GMRT channels 1 and 2 are true R and L, however, GMRT being a prime focus instrument, the reflection due to the dish reverses the sense of polarisation and converts the right circular polarisation (RCP) into left circular polarisation (LCP) and vice versa. This has not been taken into account and thus the GMRT Stokes V and U signs need to be reversed to make them consistent with the IAU/IEEE convention. This objective can be achieved by reassigning channel 1 to L and channel 2 to R for all circular feeds, i.e. bands 2, 3, and 4. The study remains inconclusive for the GMRT band 5, which has linear feeds.


title: Applying Reinforcement Learning to Option Pricing and Hedging

abstract: This thesis provides an overview of the recent advances in reinforcement learning in pricing and hedging financial instruments, with a primary focus on a detailed explanation of the Q-Learning Black Scholes approach, introduced by Halperin (2017). This reinforcement learning approach bridges the traditional Black and Scholes (1973) model with novel artificial intelligence algorithms, enabling option pricing and hedging in a completely model-free and data-driven way. This paper also explores the algorithm's performance under different state variables and scenarios for a European put option. The results reveal that the model is an accurate estimator under different levels of volatility and hedging frequency. Moreover, this method exhibits robust performance across various levels of option's moneyness. Lastly, the algorithm incorporates proportional transaction costs, indicating diverse impacts on profit and loss, affected by different statistical properties of the state variables.


title: In Vitro and In Silico Characterization of the Aggregation of Thrombi on Ventricular Assist Device Cannula

abstract: The unacceptably high stroke rate of HeartMate III VAD without signs of adherent pump thrombosis is hypothesized to be the result of the thrombi originating on the inflow cannula, ingesting and ejecting emboli from the VAD. Therefore, inflow cannula thrombosis has been an emerging focus. The inflow cannula of contemporary VADs, which incorporate both polished and rough regions serve as useful benchmarks to study the effects of roughness and shear on thrombogenesis. An in vitro study was conducted to emulate the micro-hemodynamic condition on a sintered inflow cannula, and to observe the deposition and detachment patterns. Together with a computational fluid dynamic tool, this study aimed to provide insight into the optimization of inflow cannula and potentially reducing adverse neurological events due to upstream thrombus.


title: Near optimal bounds for weak and strong spatial mixing for the anti-ferromagnetic Potts model on trees

abstract: We show that the anti-ferromagnetic Potts model on trees exhibits strong spatial mixing for a near-optimal range of parameters. Our work complements recent results of Chen, Liu, Mani, and Moitra [arXiv.2304.01954] who showed this to be true in the infinite temperature setting, corresponding to uniform proper colorings. We furthermore prove weak spatial mixing results complementing results in [arXiv.2304.01954].


title: Kagome KMn$_3$Sb$_5$ metal: Magnetism, lattice dynamics, and anomalous Hall conductivity

abstract: Kagome metals are reported to exhibit remarkable properties, including superconductivity, charge density wave order, and a large anomalous Hall conductivity, which facilitate the implementation of spintronic devices. In this work, we study a novel kagome metal based on Mn magnetic sites in a KMn$_3$Sb$_5$ stoichiometry. By means of first-principles density functional theory calculations, we demonstrate that the studied compound is dynamically stable, locking the ferromagnetic order as the ground state configuration, thus preventing the charge-density-wave state as reported in its vanadium-based counterpart KV$_3$Sb$_5$. Our calculations predict that KMn$_3$Sb$_5$ exhibits an out-of-plane (001) ferromagnetic response as the ground state, allowing for the emergence of topologically protected Weyl nodes near the Fermi level and nonzero anomalous Hall conductivity ($\sigma_{ij}$) in this centrosymmetric system. We obtain a tangible $\sigma_{xy} = 314$ S$\cdot$cm$^{-1}$ component, which is comparable to that of other kagome metals. Finally, we explore the effect of the on-site Coulomb repulsion ($+U$) on the structural and electronic properties and find that, although the lattice parameters and $\sigma_{xy}$ moderately vary with increasing $+U$, KMn$_3$Sb$_5$ stands as an ideal stable ferromagnetic kagome metal with a large anomalous Hall conductivity response.


title: On Tractable Convex Relaxations of Standard Quadratic Optimization Problems under Sparsity Constraints

abstract: Standard quadratic optimization problems (StQPs) provide a versatile modelling tool in various applications. In this paper, we consider StQPs with a hard sparsity constraint, referred to as sparse StQPs. We focus on various tractable convex relaxations of sparse StQPs arising from a mixed-binary quadratic formulation, namely, the linear optimization relaxation given by the reformulation-linearization technique, the Shor relaxation, and the relaxation resulting from their combination. We establish several structural properties of these relaxations in relation to the corresponding relaxations of StQPs without any sparsity constraints, and pay particular attention to the rank-one feasible solutions retained by these relaxations. We then utilize these relations to establish several results about the quality of the lower bounds arising from different relaxations. We also present several conditions that ensure the exactness of each relaxation.


title: Holomorphic bundles framed along a real hypersurface and the Riemann-Hilbert problem

abstract: Let $X$ be a connected, compact complex manifold and $S\subset X$ a separating real hypersurface, so that $X$ decomposes as a union of compact complex manifolds with boundary $\bar X^\pm$. Let $\mathcal{M}$ be the moduli space of $S$-framed holomorphic bundles, i.e. of pairs $(E,\theta)$ of fixed topological type consisting of a holomorphic bundle $E$ on $X$ and a trivialization $\theta$ - belonging to a fixed H\"older regularity class $\mathcal{C}^{\kappa+1}$ - of its restriction to $S$. The restrictions to $\bar X^\pm$ of an $S$-framed holomorphic bundle $(E,\theta)$ are boundary framed formally holomorphic bundles $(E^\pm,\theta^\pm)$ which induce, via $\theta^\pm$, the same tangential Cauchy-Riemann operators on the trivial bundle on $S$, so one obtains a natural map from $\mathcal{M}$ into the fiber product $\mathcal{M}^-\times_\mathcal{C}\mathcal{M}^+$ over the space $\mathcal{C}$ of Cauchy-Riemann operators on the trivial bundle on $S$. Our main result states: this map is a homeomorphism for $\kappa\in (0,\infty]\setminus\mathbb{N}$. The proof is based on a gluing principle for formally holomorphic bundles along a real hypersurface. This principle can also be used to give a complex geometric interpretation of the space of solutions of a large class of Riemann-Hilbert type problems. The results generalize in two directions: first one can replace the decomposition associated with a separating hypersurface by the the manifold with boundary $\widehat X_S$ obtained by cutting $X$ along any oriented hypersurface $S$. Second one can consider principal $G$ bundles for an arbitrary complex Lie group $G$. We give explicit examples of moduli spaces of (boundary) framed holomorphic bundles and explicit formulae for the homeomorphisms provided by the general results.


title: Minerva: Decentralized Collaborative Query Processing over InterPlanetary File System

abstract: Data silos create barriers in accessing and utilizing data dispersed over networks. Directly sharing data easily suffers from the long downloading time, the single point failure and the untraceable data usage. In this paper, we present Minerva, a peer-to-peer cross-cluster data query system based on InterPlanetary File System (IPFS). Minerva makes use of the distributed Hash table (DHT) lookup to pinpoint the locations that store content chunks. We theoretically model the DHT query delay and introduce the fat Merkle tree structure as well as the DHT caching to reduce it. We design the query plan for read and write operations on top of Apache Drill that enables the collaborative query with decentralized workers. We conduct comprehensive experiments on Minerva, and the results show that Minerva achieves up to $2.08 \times$ query performance acceleration compared to the original IPFS data query, and could complete data analysis queries on the Internet-like environments within an average latency of $0.615$ second. With collaborative query, Minerva could perform up to $1.39 \times$ performance acceleration than centralized query with raw data shipment.


title: Functional Geometry Guided Protein Sequence and Backbone Structure Co-Design

abstract: Proteins are macromolecules responsible for essential functions in almost all living organisms. Designing reasonable proteins with desired functions is crucial. A protein's sequence and structure are strongly correlated and they together determine its function. In this paper, we propose NAEPro, a model to jointly design Protein sequence and structure based on automatically detected functional sites. NAEPro is powered by an interleaving network of attention and equivariant layers, which can capture global correlation in a whole sequence and local influence from nearest amino acids in three dimensional (3D) space. Such an architecture facilitates effective yet economic message passing at two levels. We evaluate our model and several strong baselines on two protein datasets, $\beta$-lactamase and myoglobin. Experimental results show that our model consistently achieves the highest amino acid recovery rate, TM-score, and the lowest RMSD among all competitors. These findings prove the capability of our model to design protein sequences and structures that closely resemble their natural counterparts. Furthermore, in-depth analysis further confirms our model's ability to generate highly effective proteins capable of binding to their target metallocofactors. We provide code, data and models in Github.


title: Exploring Transition from Stability to Chaos through Random Matrices

abstract: This study explores the application of random matrices to track chaotic dynamics within the Chirikov standard map. Our findings highlight the potential of matrices exhibiting Wishart-like characteristics, combined with statistical insights from their eigenvalue density, as a promising avenue for chaos monitoring. Inspired by a technique originally designed for detecting phase transitions in spin systems, we successfully adapt and apply it to identify analogous transformative patterns in the context of the Chirikov standard map. Leveraging the precision previously demonstrated in localizing critical points within magnetic systems in our prior research, our method accurately pinpoints the Chirikov resonance-overlap criterion for the chaos boundary at $K\approx 2.43$, reinforcing its effectiveness.


title: Neur2RO: Neural Two-Stage Robust Optimization

abstract: Robust optimization provides a mathematical framework for modeling and solving decision-making problems under worst-case uncertainty. This work addresses two-stage robust optimization (2RO) problems (also called adjustable robust optimization), wherein first-stage and second-stage decisions are made before and after uncertainty is realized, respectively. This results in a nested min-max-min optimization problem which is extremely challenging computationally, especially when the decisions are discrete. We propose Neur2RO, an efficient machine learning-driven instantiation of column-and-constraint generation (CCG), a classical iterative algorithm for 2RO. Specifically, we learn to estimate the value function of the second-stage problem via a novel neural network architecture that is easy to optimize over by design. Embedding our neural network into CCG yields high-quality solutions quickly as evidenced by experiments on two 2RO benchmarks, knapsack and capital budgeting. For knapsack, Neur2RO finds solutions that are within roughly $2\%$ of the best-known values in a few seconds compared to the three hours of the state-of-the-art exact branch-and-price algorithm; for larger and more complex instances, Neur2RO finds even better solutions. For capital budgeting, Neur2RO outperforms three variants of the $k$-adaptability algorithm, particularly on the largest instances, with a 10 to 100-fold reduction in solution time. Our code and data are available at https://github.com/khalil-research/Neur2RO.


title: A Serverless Architecture for Efficient and Scalable Monte Carlo Markov Chain Computation

abstract: Computer power is a constantly increasing demand in scientific data analyses, in particular when Markov Chain Monte Carlo (MCMC) methods are involved, for example for estimating integral functions or Bayesian posterior probabilities. In this paper, we describe the benefits of a parallel computation of MCMC using a cloud-based, serverless architecture: first, the computation time can be spread over thousands of processes, hence greatly reducing the time the user should wait to have its computation completed. Second, the overhead time required for running in parallel several processes is minor and grows logarithmically with respect to the number of processes. Third, the serverless approach does not require time-consuming efforts for maintaining and updating the computing infrastructure when/if the number of walkers increases or for adapting the code to optimally use the infrastructure. The benefits are illustrated with the computation of the posterior probability distribution of a real astronomical analysis.