Search

Stochastic differential equations (SDEs) are a fundamental tool for modelling dynamic processes, including gene regulatory networks (GRNs), contaminant transport, financial markets, and image generation. However, learning the underlying SDE from observational data is a challenging task, especially when individual trajectories are not observable. Motivated by burgeoning research in single-cell datasets, we present the first comprehensive approach for jointly estimating the drift and diffusion of an SDE from its temporal marginals. Assuming linear drift and additive diffusion, we prove that these parameters are identifiable from marginals if and only if the initial distribution is not invariant to a class of generalized rotations, a condition that is satisfied by most distributions. We further prove that the causal graph of any SDE with additive diffusion can be recovered from the SDE parameters. To complement this theory, we adapt entropy-regularized optimal transport to handle anisotropic diffusion, and introduce APPEX (Alternating Projection Parameter Estimation from $X_0$), an iterative algorithm designed to estimate the drift, diffusion, and causal graph of an additive noise SDE, solely from temporal marginals. We show that each of these steps are asymptotically optimal with respect to the Kullback-Leibler divergence, and demonstrate APPEX's effectiveness on simulated data from linear additive noise SDEs.

We introduce CODES, a benchmark for comprehensive evaluation of surrogate architectures for coupled ODE systems. Besides standard metrics like mean squared error (MSE) and inference time, CODES provides insights into surrogate behaviour across multiple dimensions like interpolation, extrapolation, sparse data, uncertainty quantification and gradient correlation. The benchmark emphasizes usability through features such as integrated parallel training, a web-based configuration generator, and pre-implemented baseline models and datasets. Extensive documentation ensures sustainability and provides the foundation for collaborative improvement. By offering a fair and multi-faceted comparison, CODES helps researchers select the most suitable surrogate for their specific dataset and application while deepening our understanding of surrogate learning behaviour., Comment: 12 pages, 10 figures, accepted for the Machine Learning and the Physical Sciences workshop at NeurIPS 2024, source code available on GitHub at https://github.com/robin-janssen/CODES-Benchmark

The Workflows Community Summit gathered 111 participants from 18 countries to discuss emerging trends and challenges in scientific workflows, focusing on six key areas: time-sensitive workflows, AI-HPC convergence, multi-facility workflows, heterogeneous HPC environments, user experience, and FAIR computational workflows. The integration of AI and exascale computing has revolutionized scientific workflows, enabling higher-fidelity models and complex, time-sensitive processes, while introducing challenges in managing heterogeneous environments and multi-facility data dependencies. The rise of large language models is driving computational demands to zettaflop scales, necessitating modular, adaptable systems and cloud-service models to optimize resource utilization and ensure reproducibility. Multi-facility workflows present challenges in data movement, curation, and overcoming institutional silos, while diverse hardware architectures require integrating workflow considerations into early system design and developing standardized resource management tools. The summit emphasized improving user experience in workflow systems and ensuring FAIR workflows to enhance collaboration and accelerate scientific discovery. Key recommendations include developing standardized metrics for time-sensitive workflows, creating frameworks for cloud-HPC integration, implementing distributed-by-design workflow modeling, establishing multi-facility authentication protocols, and accelerating AI integration in HPC workflow management. The summit also called for comprehensive workflow benchmarks, workflow-specific UX principles, and a FAIR workflow maturity model, highlighting the need for continued collaboration in addressing the complex challenges posed by the convergence of AI, HPC, and multi-facility research environments.

We simulate the charging of a single electrolyte-filled pore using the modified Poisson-Nernst-Planck and Navier-Stokes equations. We find that electroconvection, previously ignored in this context, can substantially speed up the charging dynamics. We derive an analytical model that describes the induced fluid velocity and the electric current arising due to convection. Our findings suggest that convection becomes significant beyond a certain threshold voltage that is an inherent electrolyte property.

The first very long baseline interferometry (VLBI) detections at 870$\mu$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescopes in Chile, Hawaii, and Spain, obtained during observations in October 2018. The longest-baseline detections approach 11$\,$G$\lambda$ corresponding to an angular resolution, or fringe spacing, of 19$\mu$as. The Allan deviation of the visibility phase at 870$\mu$m is comparable to that at 1.3$\,$mm on the relevant integration time scales between 2 and 100$\,$s. The detections confirm that the sensitivity and signal chain stability of stations in the Event Horizon Telescope (EHT) array are suitable for VLBI observations at 870$\mu$m. Operation at this short wavelength, combined with anticipated enhancements of the EHT, will lead to a unique high angular resolution instrument for black hole studies, capable of resolving the event horizons of supermassive black holes in both space and time., Comment: Corresponding author: S. Doeleman

We present observations of CO(1--0) and CO(2--1) lines from the Institut de radioastronomie millim\'etrique (IRAM) 30m telescope toward 20 nearby, optically luminous type 2 quasars (QSO2s) and observations of [C II] 158$\mu$m line from the Stratospheric Observatory For Infrared Astronomy (SOFIA) for 5 QSO2s in the CO sample and 5 type 1 quasars (QSO1s). In the traditional evolutionary scenario explaining different types of QSOs, obscured QSO2s emerge from gas-rich mergers observed as luminous infrared galaxies (LIRGs) and then turn into unobscured QSO1s as the black holes clear out the obscuring material in a blow-out phase. We test the validity of this theoretical prediction by comparing the gas fractions and star formation efficiencies among LIRGs and QSOs. We find that CO luminosity, CO-derived gas masses and gas fractions in QSO1s are consistent with those estimated for QSO2s, while LIRGs exhibit a closer resemblance to QSO2s in terms of CO-derived gas masses and gas fractions, and [C II] luminosity. However, comparisons between [C II] luminosity and star formation tracers such as the CO and infrared luminosity imply additional sources of [C II] emission in QSO1s likely tracing neutral atomic or ionized gas. All three types of galaxies have statistically indistinguishable distributions of star formation efficiency. Our results are consistent with part of the evolutionary scenario where nearby QSO2s could emerge from LIRGs, but they are unlikely to be the precursors of nearby QSO1s., Comment: 32 pages, 10 figures, 7 tables; the complete set of Figure 1 is appended to the end of document. Under review (revisions submitted) by ApJ

Bayesian optimization (BO) is a powerful and data-efficient method for iterative materials discovery and design, particularly valuable when prior knowledge is limited, underlying functional relationships are complex or unknown, and the cost of querying the materials space is significant. Traditional BO methodologies typically utilize conventional Gaussian Processes (cGPs) to model the relationships between material inputs and properties, as well as correlations within the input space. However, cGP-BO approaches often fall short in multi-objective optimization scenarios, where they are unable to fully exploit correlations between distinct material properties. Leveraging these correlations can significantly enhance the discovery process, as information about one property can inform and improve predictions about others. This study addresses this limitation by employing advanced kernel structures to capture and model multi-dimensional property correlations through multi-task (MTGPs) or deep Gaussian Processes (DGPs), thus accelerating the discovery process. We demonstrate the effectiveness of MTGP-BO and DGP-BO in rapidly and robustly solving complex materials design challenges that occur within the context of complex multi-objective optimization -- carried out by leveraging the pyiron workflow manager over FCC FeCrNiCoCu high entropy alloy (HEA) spaces, where traditional cGP-BO approaches fail. Furthermore, we highlight how the differential costs associated with querying various material properties can be strategically leveraged to make the materials discovery process more cost-efficient.

Dark Matter (DM) remains poorly probed on critical, sub-galactic scales, where predictions from different models diverge in terms of abundance and density profiles of halos. Gravitational lens systems on milli-arcsecond scales (milli-lenses) are expected for a population of dense DM halos, or free-floating supermassive black holes (SMBHs), that might be comprised of primordial black holes (PBHs), in the mass range of $10^6$ to $10^9 M_\odot$. In this paper, we aim to look for milli-lens systems via a systematic search in a large sample of radio-loud AGN observed with very-long-baseline interferometry (VLBI). We present the observational strategy to discriminate milli-lenses from contaminant objects mimicking a milli-lens morphology. In a pilot project, we have investigated VLBI images from 13,828 sources from the Astrogeo VLBI image database and reduced the number of candidates to 40 in a first step. We present here the images and analysis of sensitive follow-up observations with the EVN at 5 and 22 GHz, and streamline our analysis to reject milli-lens candidates. Using constraints such as the surface brightness ratio, conservation of spectral shape, stability of flux ratios over time, and changes in morphology, we can confidently discriminate between milli-lenses and their mimickers. Using the above constraints, we rule out 32 out of our initial 40 candidates as milli-lenses, demonstrating the power of our approach. Also, we find many new candidates for compact symmetric objects, that are thought to be short-lived, jetted radio sources. This serves as a pathfinder for the final sample used for the SMILE (Search for MIlli-LEnses) project, which will allow us to constrain DM models by comparing the results to theoretical predictions. This SMILE sample will consist of $\sim$5,000 sources based on the VLA CLASS survey, including many observations obtained for this project specifically., Comment: 35 pages, 36 figures, submitted for publication in Astronomy & Astrophysics

Current instruction-tuned language models are exclusively trained with textual preference data and thus are often not aligned with the unique requirements of other modalities, such as speech. To better align language models with the speech domain, we explore (i) prompting strategies grounded in radio-industry best practices and (ii) preference learning using a novel speech-based preference data of 20K samples, generated with a wide spectrum of prompts that induce varying dimensions of speech-suitability and labeled by annotators who listen to response pairs. Both human and automatic evaluation show that both prompting and preference learning increase the speech-suitability of popular instruction-tuned LLMs. Interestingly, we find that prompting and preference learning can be additive; combining them achieves the best win rates in head-to-head comparison, resulting in responses that are preferred or tied to the base model in 76.2% of comparisons on average. Lastly, we share lexical, syntactical, and qualitative analyses to showcase how each method contributes to improving the speech-suitability of generated responses., Comment: EMNLP2024

The candidate Kitaev magnet Na$_2$Co$_2$TeO$_6$ exhibits a characteristic ferrimagnetic response at low temperatures, with a finite residual magnetization that changes sign at a compensation point located at around half the ordering temperature. We argue that the behavior can be naturally understood to arise in this material as a consequence of a noncollinear triple-$\mathbf{q}$ magnetic ground state. Using large-scale classical Monte Carlo simulations, we study the finite-temperature response of the pertinent honeycomb Heisenberg-Kitaev-$\Gamma$-$\Gamma'$ model in weak training fields. Our model features all symmetry-allowed nearest-neighbor exchange interactions, as well as sublattice-dependent next-nearest-neighbor interactions, consistent with the reported crystal structure of the material. We also consider a six-spin ring exchange perturbation, which allows us to tune between the two different magnetic long-range orders that have been suggested for this material in the literature, namely, a collinear single-$\mathbf{q}$ zigzag state and a noncollinear triple-$\mathbf{q}$ state. We demonstrate that the experimentally-observed ferrimagnetic response of Na$_2$Co$_2$TeO$_6$ can be well described within our modeling if the magnetic ground state features noncollinear triple-$\mathbf{q}$ order. By contrast, a classical Heisenberg-Kitaev-$\Gamma$-$\Gamma'$-type model with collinear zigzag ground state is insufficient even in principle to describe the observed behavior. Our results illustrate the unconventional physics of noncollinear magnetic long-range orders hosted by frustrated magnets with bond-dependent interactions., Comment: 12 pages, 6 figures

Efficient attention deployment in visual search is limited by human visual memory, yet this limitation can be offset by exploiting the environment's structure. This paper introduces a computational cognitive model that simulates how the human visual system uses visual hierarchies to prevent refixations in sequential attention deployment. The model adopts computational rationality, positing behaviors as adaptations to cognitive constraints and environmental structures. In contrast to earlier models that predict search performance for hierarchical information, our model does not include predefined assumptions about particular search strategies. Instead, our model's search strategy emerges as a result of adapting to the environment through reinforcement learning algorithms. In an experiment with human participants we test the model's prediction that structured environments reduce visual search times compared to random tasks. Our model's predictions correspond well with human search performance across various set sizes for both structured and unstructured visual layouts. Our work improves understanding of the adaptive nature of visual search in hierarchically structured environments and informs the design of optimized search spaces.

Pyrochlore magnets have proven to provide an excellent arena for the realization of a variety of many-body phenomena such as classical and quantum order-by-disorder, as well as spin liquid phases described by emergent gauge field theories. These phenomena arise from the competition between different symmetry-breaking magnetic orders. In this work, we consider a subspace of the most general bilinear nearest-neighbor Hamiltonian on the pyrochlore lattice, parameterized by the local interaction parameter $J_{z\pm}$, where three symmetry-breaking phases converge. We demonstrate that for small values of $|J_{z\pm}|$, a conventional $\mathbf q=0$ ordered phase is selected by a thermal order-by-disorder mechanism. For $|J_{z\pm}|$ above a certain finite threshold, a novel spin-nematic phase is stabilized at low temperatures. Instead of the usual Bragg peaks, the spin-nematic phase features lines of high intensity in the spin structure factor. At intermediate temperatures above the low-temperature orders, a rank-2 U(1) classical spin liquid is realized for all $J_{z\pm} \neq 0$. We fully characterize all phases using classical Monte-Carlo simulations and a self-consistent Gaussian approximation., Comment: 17 pages, 17 figures

MICADO is a first light instrument for the Extremely Large Telescope (ELT), set to start operating later this decade. It will provide diffraction limited imaging, astrometry, high contrast imaging, and long slit spectroscopy at near-infrared wavelengths. During the initial phase operations, adaptive optics (AO) correction will be provided by its own natural guide star wavefront sensor. In its final configuration, that AO system will be retained and complemented by the laser guide star multi-conjugate adaptive optics module MORFEO (formerly known as MAORY). Among many other things, MICADO will study exoplanets, distant galaxies and stars, and investigate black holes, such as Sagittarius A* at the centre of the Milky Way. After their final design phase, most components of MICADO have moved on to the manufacturing and assembly phase. Here we summarize the final design of the instrument and provide an overview about its current manufacturing status and the timeline. Some lessons learned from the final design review process will be presented in order to help future instrumentation projects to cope with the challenges arising from the substantial differences between projects for 8-10m class telescopes (e.g. ESO-VLT) and the next generation Extremely Large Telescopes (e.g. ESO-ELT). Finally, the expected performance will be discussed in the context of the current landscape of astronomical observatories and instruments. For instance, MICADO will have similar sensitivity as the James Webb Space Telescope (JWST), but with six times the spatial resolution., Comment: Proceedings of the SPIE, Volume 13096, id. 1309611 11 pp. (2024)

Cell motility in dense cell collectives is pivotal in various diseases like cancer metastasis and asthma. A central aspect in these phenomena is the heterogeneity in cell motility, but identifying the motility of individual cells is challenging. Previous work has established the importance of the average cell shape in predicting cell dynamics. Here, we aim to identify the importance of individual cell shape features, rather than collective features, to distinguish between high-motility (active) and low-motility (passive) cells in heterogeneous cell layers. Employing the Cellular Potts Model, we generate simulation snapshots and extract static features as inputs for a simple machine-learning model. Our results show that when the passive cells are non-motile, this machine-learning model can accurately predict whether a cell is passive or active using only single-cell shape features. Furthermore, we explore scenarios where passive cells also exhibit some degree of motility, albeit less than active cells. In such cases, our findings indicate that a neural network trained on shape features can accurately classify cell motility, particularly when the number of active cells is low, and the motility of active cells is significantly higher compared to passive cells. This work offers potential for physics-inspired predictions of single-cell properties with implications for inferring cell dynamics from static histological images.

The mode-coupling theory (MCT) of the glass transition ranks among the most successful first-principles kinetic theories to describe glassy dynamics. However, MCT does not fully account for crucial aspects of the dynamics near the glass transition. To facilitate improving the theory, we critically test the approximations inherent in MCT for a supercooled mixture using Brownian dynamics simulations. Although each MCT approximation significantly impacts the predicted dynamics, our findings show that long-time cancellations of errors occur due to static and dynamic approximations of four-point correlation functions, with the validity of these approximations remaining relatively constant across different temperatures. Notably, the MCT form of the memory functional maintains remarkably high accuracy even in the supercooled regime when evaluated with the intermediate scattering function from simulations. The primary discrepancies between theoretical predictions and experimental results arise from the self-consistent nature of the MCT equations, which amplify minor errors in the memory kernel. This suggests that minor corrections to the memory functional could substantially enhance the theory's predictive accuracy.

Principal component analysis (PCA) is one of the most popular dimension reduction techniques in statistics and is especially powerful when a multivariate distribution is concentrated near a lower-dimensional subspace. Multivariate extreme value distributions have turned out to provide challenges for the application of PCA since their constraint support impedes the detection of lower-dimensional structures and heavy-tails can imply that second moments do not exist, thereby preventing the application of classical variance-based techniques for PCA. We adapt PCA to max-stable distributions using a regression setting and employ max-linear maps to project the random vector to a lower-dimensional space while preserving max-stability. We also provide a characterization of those distributions which allow for a perfect reconstruction from the lower-dimensional representation. Finally, we demonstrate how an optimal projection matrix can be consistently estimated and show viability in practice with a simulation study and application to a benchmark dataset., Comment: 34 pages, 11 Figures

Critical single-particle fluctuations associated with particle displacements are inherent to simple glass-forming liquids in the limit of large dimensions and leave a pseudo-critical trace across all finite dimensions. This characteristic could serve as a crucial test for distinguishing between theories of glass formation. We here examine these critical fluctuations, as captured by the well-established non-Gaussian parameter, within both mode-coupling theory (MCT) and dynamical mean-field theory (DMFT) across dimensions for hard sphere liquids and for the minimally structured Mari--Kurchan model. We establish general scaling laws relevant to any liquid dynamics theory in large dimensions and show that the dimensional scalings predicted by MCT are inconsistent with those from DMFT. Simulation results for hard spheres in moderately high dimensions align with the DMFT scenario, reinforcing the relevance of mean-field theory for capturing glass physics in finite dimensions. We identify potential adjustments to MCT to account for certain mean-field physics. Our findings also highlight that local structure and spatial dimensionality can affect single-particle critical fluctuations in non-trivial ways., Comment: 42 pages, 9 figures

In this study, we focus on different strategies drivers use in terms of interleaving between driving and non-driving related tasks (NDRT) while taking back control from automated driving. We conducted two driving simulator experiments to examine how different cognitive demands of texting, priorities, and takeover time budgets affect drivers' takeover strategies. We also evaluated how different takeover strategies affect takeover performance. We found that the choice of takeover strategy was influenced by the priority and takeover time budget but not by the cognitive demand of the NDRT. The takeover strategy did not have any effect on takeover quality or NDRT engagement but influenced takeover timing.

$\delta$-Covering, for some covering range $\delta>0$, is a continuous facility location problem on undirected graphs where all edges have unit length. The facilities may be positioned on the vertices as well as on the interior of the edges. The goal is to position as few facilities as possible such that every point on every edge has distance at most $\delta$ to one of these facilities. For large $\delta$, the problem is similar to dominating set, which is hard to approximate, while for small $\delta$, say close to $1$, the problem is similar to vertex cover. In fact, as shown by Hartmann et al. [Math. Program. 22], $\delta$-Covering for all unit-fractions $\delta$ is polynomial time solvable, while for all other values of $\delta$ the problem is NP-hard. We study the approximability of $\delta$-Covering for every covering range $\delta>0$. For $\delta \geq 3/2$, the problem is log-APX-hard, and allows an $\mathcal O(\log n)$ approximation. For every $\delta < 3/2$, there is a constant factor approximation of a minimum $\delta$-cover (and the problem is APX-hard when $\delta$ is not a unit-fraction). We further study the dependency of the approximation ratio on the covering range $\delta < 3/2$. By providing several polynomial time approximation algorithms and lower bounds under the Unique Games Conjecture, we narrow the possible approximation ratio, especially for $\delta$ close to the polynomial time solvable cases., Comment: 22 pages, 6 figures, extended abstract accepted at WAOA24

The Probe far-Infrared Mission for Astrophysics (PRIMA) is a proposed space observatory which will use arrays of thousands of kinetic inductance detectors (KIDs) to perform low- and moderate-resolution spectroscopy throughout the far-infrared. The detectors must have noise equivalent powers (NEPs) at or below 0.1 aW/sqrt(Hz) to be subdominant to noise from sky backgrounds and thermal noise from PRIMA's cryogenically cooled primary mirror. Using a Radio Frequency System on a Chip for multitone readout, we measure the NEPs of detectors on a flight-like array designed to observe at a wavelength of 210 microns. We find that 92% of the KIDs measured have an NEP below 0.1 aW/sqrt(Hz) at a noise frequency of 10 Hz., Comment: 7 pages, 7 figures

Ultra diffuse galaxies, characterized by their low surface brightness and large physical size, constitute a subclass of dwarf galaxies that challenge our current understanding of galaxy formation and evolution. In this paper, we probe the properties of 74 UDGs, identified in the MATLAS survey, based on a comprehensive study of their globular cluster (GC) populations. We obtained high resolution HST imaging of these galaxies using the ACS F606W and F814W filters, allowing us to select GCs based on color and concentration index. After background subtraction and completeness correction, we calculate an overall total of 387 GCs. The number of GCs per galaxy ranges from 0 to 38, with the majority (64%) having low counts (0-2 GCs). On average, the more massive UDGs host a larger number of GCs. We find that our UDGs have specific frequencies (S_N) ranging from 0 to 91, with a small population (9%) with S_N > 30. The median S_N of our sample is similar to the one for the Perseus cluster UDGs, despite the fact that our UDGs are found in lower density environments. The S_N measurements for individual galaxies can extend beyond those found in Perseus, but remain below the values found for UDGs in the Virgo and Coma cluster. Based on a trending analysis of the S_N values with the host galaxy properties, we find trends with host galaxy size, roundness, color, and local density. For the UDGs with sufficiently high statistics, we study 2D density maps of the GC distributions, which show a variety of appearances: symmetric, asymmetric, off-center, and elongated. The UDGs with disturbed density maps also show disturbed stellar light morphologies. We further quantify the distribution by modeling it with a Sersic profile, finding R_{e,GC}/R_{e,gal} ~ 1.0, which indicates that the GCs follow the stellar light of the host galaxy., Comment: 20 pages, 17 figures, 2 tables, accepted for publication in A&A

We show that the Kontsevich-Segal-Witten (KSW) criterion applied to the no-boundary state constrains anisotropic deformations of de Sitter space. We consider squashed $S^3$ and $S^1 \times S^2$ boundaries and find that in both models, the KSW criterion excludes a significant range of homogeneous but anisotropic configurations. For squashed $S^3$ boundaries, the excluded range includes all surface geometries with negative scalar curvature, in line with dS/CFT reasoning. For $S^1 \times S^2$ boundaries, we find that KSW selects the low-temperature regime of configuration space where the $S^1$ is sufficiently large compared to the $S^2$. In both models, the KSW criterion renders the semiclassical wave function normalizable, up to one-loop effects., Comment: 33 pages

This work introduces Dareplane, a modular and broad technology agnostic open source software platform for brain-computer interface research with an application focus on adaptive deep brain stimulation (aDBS). While the search for suitable biomarkers to inform aDBS has provided rich results over the last two decades, development of control strategies is not progressing at the same pace. One difficulty for investigating control approaches resides with the complex setups required for aDBS experiments. The Dareplane platform supports aDBS setups, and more generally brain computer interfaces, by providing a modular, technology-agnostic, and easy-to-implement software platform to make experimental setups more resilient and replicable. The key features of the platform are presented and the composition of modules into a full experimental setup is discussed in the context of a Python-based orchestration module. The performance of a typical experimental setup on Dareplane for aDBS is evaluated in three benchtop experiments, covering (a) an easy-to-replicate setup using an Arduino microcontroller, (b) a setup with hardware of an implantable pulse generator, and (c) a setup using an established and CE certified external neurostimulator. Benchmark results are presented for individual processing steps and full closed-loop processing. The results show that the microcontroller setup in (a) provides timing comparable to the realistic setups in (b) and (c). The Dareplane platform was successfully used in a total of 19 open-loop DBS sessions with externalized DBS and electrocorticography (ECoG) leads. In addition, the full technical feasibility of the platform in the aDBS context is demonstrated in a first closed-loop session with externalized leads on a patient with Parkinson's disease receiving DBS treatment.

By 2050, the global population will exceed 9 billion, demanding a 70% increase in food production. Animal proteins alone may not suffice and contribute to global warming. Alternative proteins such as legumes, algae, and insects are being explored, but their health impacts are largely unknown. For this, three-week-old rats were fed diets containing 20% protein from various sources for six weeks. A casein-based control diet was compared to soy isolate, spirulina powder, chickpea isolate, chickpea flour, and fly larvae powder. Except for spirulina, alternative protein groups showed comparable growth patterns to the casein group. Morphological and mechanical tests of femur bones matched growth patterns. Caecal 16S analysis highlighted the impact on gut microbiota diversity. Chickpea flour showed significantly lower $\alpha$-diversity compared with casein and chickpea isolate groups while chickpea flour, had the greatest distinction in $\beta$-diversity. Alternative protein sources supported optimal growth, but quality and health implications require further exploration.

The Ornstein-Zernike equation is a powerful tool in liquid state theory for predicting structural and thermodynamic properties of fluids. Combined with a suitable closure, it has been shown to reproduce e.g. the static structure factor, pressure, and compressibility of liquids to a great degree of accuracy. However, out of the multitude of closures that exist for the Ornstein-Zernike equation, it is hard to predict a priori which closure will give the most accurate predictions for the system at hand. To alleviate this problem, we compare the predictive power of many closures on a curated set of representative benchmark systems, including those with hard-sphere, inverse power-law, Gaussian core, and Lennard-Jones particles, in three and two dimensions. For example, we find that the well-known and highly used Percus-Yevick closure gives significantly worse predictions than lesser-known closures of equal complexity in all cases studied. We anticipate that the trends observed in our results will aid in making more informed decisions regarding closure choices. To facilitate the adoption of more modern closure theories, we also have packaged, documented, and distributed the code necessary to numerically solve the equations for a given closure and pair interaction potential.

We report on the effects of cosmic ray interactions with the Kinetic Inductance Detector (KID) based focal plane array for the Terahertz Intensity Mapper (TIM). TIM is a NASA-funded balloon-borne experiment designed to probe the peak of the star formation in the Universe. It employs two spectroscopic bands, each equipped with a focal plane of four $\sim\,$900-pixel, KID-based array chips. Measurements of an 864-pixel TIM array shows 791 resonators in a 0.5$\,$GHz bandwidth. We discuss challenges with resonator calibration caused by this high multiplexing density. We robustly identify the physical positions of 788 (99.6$\,$%) detectors using a custom LED-based identification scheme. Using this information we show that cosmic ray events occur at a rate of 2.1$\,\mathrm{events/min/cm^2}$ in our array. 66$\,$% of the events affect a single pixel, and another 33$\,$% affect $<\,$5 KIDs per event spread over a 0.66$\,\mathrm{cm^2}$ region (2 pixel pitches in radius). We observe a total cosmic ray dead fraction of 0.0011$\,$%, and predict that the maximum possible in-flight dead fraction is $\sim\,$0.165$\,$%, which demonstrates our design will be robust against these high-energy events., Comment: 14 pages, 5 figures. Accepted for the publication in Journal of Low Temperature Physics (2024)

Self-calibration methods with the CLEAN algorithm have been widely employed in Very Long Baseline Interferometry (VLBI) data processing in order to correct antenna-based amplitude and phase corruptions present in the data. However, human interaction during the conventional CLEAN self-calibration process can impose a strong effective prior, which in turn may produce artifacts within the final image and hinder the reproducibility of final results. In this work, we aim to demonstrate a combined self-calibration and imaging method for VLBI data in a Bayesian inference framework. The method corrects for amplitude and phase gains for each antenna and polarization mode by inferring the temporal correlation of the gain solutions. We use Stokes I data of M87 taken with the Very Long Baseline Array (VLBA) at 43GHz, pre-calibrated using the rPICARD CASA-based pipeline. For antenna-based gain calibration and imaging, we use the Bayesian imaging software resolve. To estimate gain and image uncertainties, we use a Variational Inference method. We obtain a high-resolution M87 Stokes I image at 43GHz in conjunction with antenna-based gain solutions using our Bayesian self-calibration and imaging method. The core with counter-jet structure is better resolved, and extended jet emission is better described compared to the CLEAN reconstruction. Furthermore, uncertainty estimation of the image and antenna-based gains allows us to quantify the reliability of the result. Our Bayesian self-calibration and imaging method is able to reconstruct robust and reproducible Stokes I images and gain solutions with uncertainty estimation by taking into account the uncertainty information in the data., Comment: accepted for publication in A&A

The wavelength dependent refraction of light in the atmosphere causes the chromatic dispersion of a target on the focal plane of an instrument. This is known as atmospheric dispersion, with one of the consequences being wavelength dependent flux losses which are difficult to minimise, requiring analysis in both instrument design and operations. We present Atmosphyre, a novel python package developed to characterise the impact of atmospheric dispersion on a spectrograph, with a focus on fibre multi-object spectrographs (MOS) which will be at the forefront of ground-based astronomy for the next few decades. We show example simulations and provide recommendations for minimising fibre MOS flux losses. We conclude that the guiding wavelength should typically be bluer than the observing band mid-wavelength, around 25-45% of the way through the band. The aperture should be centred on this wavelength's location on the focal plane. This wavelength/position remains constant for all reasonable declinations and target hour angles. We also present an application of the package to MOSAIC, the ELT's multi-object spectrograph. We find that differential losses greater than 10% are unavoidable for 1h observations that are a) after a local hour angle of 2.5h, or b) at declinations below -60 degrees and above 10 degrees. We identify that the introduction of an atmospheric dispersion corrector (ADC) would result in the significant reduction of spectral distortions, a gain in survey speed for many observations, and enable the implementation of wider visible observing bands; as a result, there has been a proposal to adopt ADCs at a positioner level for MOSAIC. Future work includes adding field differential refraction to Atmosphyre, important for future wide-field multi-object spectrograph projects such as the proposed WST., Comment: 10 pages, 6 figures, submitted to SPIE Astronomical Telescopes & Instrumentation, Ground-based and Airborne Instrumentation for Astronomy X (Conference 13096, Paper 219)

We characterize, by means of large-scale fermion quantum Monte Carlo simulations, metallic and deconfined quantum phase transitions in a bilayer honeycomb model in terms of their quantum critical and finite-temperature properties.The model features three different phases at zero temperature as function of interaction strength. At weak interaction, a fully symmetric Dirac semimetal state is realized. At intermediate and strong interaction, respectively, two long-range-ordered phases that break different symmetries are stabilized. The ordered phases feature partial and full, respectively, gap openings in the fermion spectrum. The first transition between the disordered and long-range-ordered semimetallic phases has previously been argued to be described by the $(2+1)$-dimensional Gross-Neveu-SO(3) field theory. By performing simulations with an improved symmetric Trotter decomposition, we further substantiate this claim by computing the critical exponents $1/\nu$, $\eta_\phi$, and $\eta_\psi$, which turn out to be consistent with the field-theoretical expectation within numerical and analytical uncertainties. The second transition between the two long-range-ordered phases has previously been proposed as a possible instance of a metallic deconfined quantum critical point. We further develop this scenario by analyzing the spectral functions in the single-particle, particle-hole, and particle-particle channels. Our results indicate gapless excitations with a unique velocity, supporting the emergence of Lorentz symmetry at criticality. We also compute the finite-temperature phase boundaries of the ordered states above the fully gapped state at large interaction. The phase boundaries vanishes smoothly in the vicinity of the putative metallic deconfined quantum critical point, in agreement with the expectation for a continuous or weakly-first-order transition., Comment: 13 pages, 9 figures

Circulating tumor cell clusters play an important role in the metastatic cascade. These clusters can acquire a migratory and more invasive phenotype, and coordinate their motion to migrate as a collective. Before such clusters can form by collectively detaching from a primary tumor, however, the cluster must first aggregate in the tumor interior. The mechanism of this cluster formation process is still poorly understood. One of the possible ways for cells to cluster is by aligning their direction of motion with their neighboring cells. This work aims to investigate the role of this cell-cell alignment interaction on the formation of motile cell clusters inside the bulk of a tumor using computer simulations. We employ a Cellular Potts model in which we model a two-dimensional heterogeneous confluent layer containing both motile and non-motile cells. Our results indicate that the degree of clustering is governed by two distinct processes: the formation of clusters due to the presence of cell-cell alignment interactions among motile cells, and the suppression of clustering due to the presence of the dynamic cellular environment (comprised of the non-motile cells). We find that the largest motile clusters are formed for intermediate alignment strengths, contrary to what is observed for motile cells in free space (that is, unimpeded by a dense cellular environment), in which case stronger cell-cell alignment always leads to larger clustering. Our findings suggest that the presence of a densely-packed cellular environment and strong cell-cell alignment inhibits the formation of large migratory clusters within the primary tumor, providing physical insight into potential factors at play during the early stages of metastasis., Comment: 12 pages, 6 figures

We extend the analytic description of complex no-boundary solutions in the context of inflation to large field models. We discuss the Kontsevich-Segal-Witten (KSW) criterion and find it is satisfied in small field models, while in large field models it depends on an integral involving $V'(\phi)$ over the range of inflation. It follows that the criterion does not truly constrain inflationary phenomenology since one can complete any inflaton potential beyond observable scales so as to satisfy KSW., Comment: 15 pages

We consider a recently proposed approach to graph signal processing based on graphons. We show how the graphon-based approach to GSP applies to graphs sampled from a stochastic block model. We obtain a basis for the graphon Fourier transform on such samples directly from the link probability matrix and the block sizes of the model. This formulation allows us to bound the sensitivity of the Fourier transform to small changes in block sizes. We then focus on the case where the probability matrix corresponds to a (weighted) Cayley graph. If block sizes are equal, a nice Fourier basis can be derived from the underlying group. We explore how, in the case where block sizes are not equal, some or all nice properties of the group basis can be maintained. We complement the theoretical results with simulations., Comment: 23 pages