The Vera C. Rubin Observatory will soon survey the southern sky, delivering a depth and sky coverage that is unprecedented in time domain astronomy. As part of commissioning, Data Preview 1 (DP1) has been released. It comprises a LSSTComCam observing campaign between November and December 2024 with multi-band imaging of seven fields, covering roughly 0.4 square degrees each, providing a first glimpse into the data products that will become available once the Legacy Survey of Space and Time begins. In this work, we search three fields for extragalactic transients. We identify eight new likely supernovae, and three known ones from a sample of 369,644 difference image analysis objects. Photometric classification using Superphot+ assigns sub-classes with >95% confidence to only one SN Ia and one SN II in this sample. Our findings are in agreement with supernova detection rate predictions of $15\pm4$ supernovae from simulations using simsurvey. The supernova detection rate in the data is possibly affected by the lack of suitable templates. Nevertheless, this work demonstrates the quality of the data products delivered in DP1 and indicates that the Rubin Observatory's Legacy Survey of Space and Time (LSST) is well placed to fulfill its discovery potential in time domain astronomy.

We present a Sugawara-type construction for boundary charges in 4d BF theory and in a general family of related TQFTs. Starting from the underlying current Lie algebra of boundary symmetries, this gives rise to well-defined quadratic charges forming an algebra of vector fields. In the case of 3d BF theory (i.e. 3d gravity), it was shown in [PRD 106 (2022), arXiv:2012.05263 [hep-th]] that this construction leads to a two-dimensional family of diffeomorphism charges which satisfy a certain modular duality. Here we show that adapting this construction to 4d BF theory first requires to split the underlying gauge algebra. Surprisingly, the space of well-defined quadratic generators can then be shown to be once again two-dimensional. In the case of tangential vector fields, this canonically endows 4d BF theory with a $\mathrm{diff}(S^2)\times\mathrm{diff}(S^2)$ or $\mathrm{diff}(S^2)\ltimes\mathrm{vect}(S^2)_\mathrm{ab}$ algebra of boundary symmetries depending on the gauge algebra. The prospect is to then understand how this can be reduced to a gravitational symmetry algebra by imposing Plebański simplicity constraints.

We introduce SpectraPyle, a versatile spectral stacking pipeline developed for the Euclid mission's NISP spectroscopic surveys, aimed at extracting faint emission lines and spectral features from large galaxy samples in the Wide and Deep Surveys. Designed for computational efficiency and flexible configuration, SpectraPyle supports the processing of extensive datasets critical to Euclid's non-cosmological science goals. We validate the pipeline using simulated spectra processed to match Euclid's expected final data quality. Stacking enables robust recovery of key emission lines, including Halpha, Hbeta, [O III], and [N II], below individual detection limits. However, the measurement of galaxy properties such as star formation rate, dust attenuation, and gas-phase metallicity are biased at stellar mass below log10(M*/Msol) ~ 9 due to the flux-limited nature of Euclid spectroscopic samples, which cannot be overcome by stacking. The SFR-stellar mass relation of the parent sample is recovered reliably only in the Deep survey for log10(M*/Msol) > 10, whereas the metallicity-mass relation is recovered more accurately over a wider mass range. These limitations are caused by the increased fraction of redshift measurement errors at lower masses and fluxes. We examine the impact of residual redshift contaminants that arises from misidentified emission lines and noise spikes, on stacked spectra. Even after stringent quality selections, low-level contamination (< 6%) has minimal impact on line fluxes due to the systematically weaker emission of contaminants. Percentile-based analysis of stacked spectra provides a sensitive diagnostic for detecting contamination via coherent spurious features at characteristic wavelengths. While our simulations include most instrumental effects, real Euclid data will require further refinement of contamination mitigation strategies.

Primordial black holes are under intense scrutiny since the detection of gravitational waves from mergers of solar-mass black holes in 2015. More recently, the development of numerical tools and the precision observational data have rekindled the effort to constrain the black hole abundance in the lower mass range, that is M &lt; 10^{23}g. In particular, primordial black holes of asteroid mass $M \sim 10^{17}-10^{23}\,$g may represent 100\% of dark matter. While the microlensing and stellar disruption constraints on their abundance have been relieved, Hawking radiation of these black holes seems to be the only detection (and constraining) mean. Hawking radiation constraints on primordial black holes date back to the first papers by Hawking. Black holes evaporating in the early universe may have generated the baryon asymmetry, modified big bang nucleosynthesis, distorted the cosmic microwave background, or produced cosmological backgrounds of stable particles such as photons and neutrinos. At the end of their lifetime, exploding primordial black holes would produce high energy cosmic rays that would provide invaluable access to the physics at energies up to the Planck scale. In this review, we describe the main principles of Hawking radiation, which lie at the border of general relativity, quantum mechanics and statistical physics. We then present an up-to-date status of the different constraints on primordial black holes that rely on the evaporation phenomenon, and give, where relevant, prospects for future work. In particular, non-standard black holes and emission of beyond the Standard Model degrees of freedom is currently a hot subject.

Markarian 231 (Mrk 231) is one of the brightest ultraluminous infrared galaxies (ULIRGs) known to date. It displays a unique optical-UV spectrum, characterized by a strong and perplexing attenuation in the near-UV, associated with a sudden polarization peak. Building on previous spectro-photometric modeling, we investigated the hypothesis that the core of Mrk 231 may host a binary SMBH system. In this scenario, the accretion disk of the primary, more massive SMBH is responsible for the optical-UV spectrum. The disk of the secondary, less massive SMBH, would be expected to essentially emit in the far UV. We applied this model to archival photometric and polarimetric data of Mrk 231 and tried to obtain the best fit possible. To support our findings, we performed radiative transfer calculations to determine the spatial disposition of each main component constituting Mrk 231. We find that a binary SMBH model can reproduce both the observed flux and polarization of Mrk 231 remarkably well. We infer that the core potentially hosts a binary SMBH system, with a primary SMBH of about 1.6x10^8 solar masses and a secondary of about 1.1x10^7 solar masses , separated by a semimajor axis of 146 AU.The secondary SMBH drives a degree of polarization of 3 % between 0.1 and 0.2 {\mu}m, with a corresponding polarization position angle of about 134{\deg} , which is consistent with scattering within an accretion disk. The primary SMBH and the structure around it are responsible for a degree of polarization of 23 % between 0.3 and 0.4 {\mu}m with a corresponding polarization position angle of about 96{\deg} , that is possibly attributed to scattering within the quasar's wind. Finally, our model predicts the existence of a second peak in polarized flux in the far-ultraviolet, a telltale signature that could definitively prove the presence of a binary SMBH.

Photometric galaxy surveys probe the late-time Universe where the density field is highly non-Gaussian. A consequence is the emergence of the super-sample covariance (SSC), a non-Gaussian covariance term that is sensitive to fluctuations on scales larger than the survey window. In this work, we study the impact of the survey geometry on the SSC and, subsequently, on cosmological parameter inference. We devise a fast SSC approximation that accounts for the survey geometry and compare its performance to the common approximation of rescaling the results by the fraction of the sky covered by the survey, $f_\mathrm{SKY}$, dubbed 'full-sky approximation'. To gauge the impact of our new SSC recipe, dubbed 'partial-sky', we perform Fisher forecasts on the parameters of the $(w_0,w_a)$-CDM model in a 3x2 points analysis, varying the survey area, the geometry of the mask and the galaxy distribution inside our redshift bins. The differences in the marginalised forecast errors, with the full-sky approximation performing poorly for small survey areas but excellently for stage-IV-like areas, are found to be absorbed by the marginalisation on galaxy bias nuisance parameters. For large survey areas, the unmarginalised errors are underestimated by about 10% for all probes considered. This is a hint that, even for stage-IV-like surveys, the partial-sky method introduced in this work will be necessary if tight priors are applied on these nuisance parameters.

Nanofluidics shows great promise for energy conversion and desalination applications. The performance of nanofluidic devices is controlled by liquid-solid friction, quantified by the Navier friction coefficient (FC). Despite decades of research, there is no well-established generic framework to determine the frequency dependent Navier FC from atomistic simulations. Here, we have derived analytical expressions to connect the Navier FC to the random force autocorrelation on the confining wall, from the observation that the random force autocorrelation can be related to the hydrodynamic boundary condition, where the Navier FC appears. The analytical framework is generic in the sense that it explicitly includes the system size dependence and also the frequency dependence of the FC, which enabled us to address (i) the long-standing plateau issue in the evaluation of the FC and (ii) the non-Markovian behavior of liquid-solid friction of a Lennard-Jones liquid and of water on various walls and at various temperatures, including the supercooled regime. This new framework opens the way to explore the frequency dependent FC for a wide range of complex liquids.

Convection is a ubiquitous process driving geophysical/astrophysical fluid flows, which are typically strongly constrained by planetary rotation on large scales. A celebrated model of such flows, rapidly rotating Rayleigh-Bénard convection, has been extensively studied in direct numerical simulations (DNS) and laboratory experiments, but the parameter values attainable by state-of-the-art methods are limited to moderately rapid rotation (Ekman numbers $Ek\gtrsim10^{-8}$), while realistic geophysical/astrophysical $Ek$ are significantly smaller. Asymptotically reduced equations of motion, the nonhydrostatic quasi-geostrophic equations (NHQGE), describing the flow evolution in the limit $Ek\to 0$, do not apply at finite rotation rates. The geophysical/astrophysical regime of small but finite $Ek$ therefore remains currently inaccessible. Here, we introduce a new, numerically advantageous formulation of the Navier-Stokes-Boussinesq equations informed by the scalings valid for $Ek\to0$, the \textit{Rescaled Rapidly Rotating incompressible Navier-Stokes Equations} (RRRiNSE). We solve the RRRiNSE using a spectral quasi-inverse method resulting in a sparse, fast algorithm to perform efficient DNS in this previously unattainable parameter regime. We validate our results against the literature across a range of $Ek$ and demonstrate that the algorithmic approaches taken remain accurate and numerically stable at $Ek$ as low as $10^{-15}$. Like the NHQGE, the RRRiNSE derive their efficiency from adequate conditioning, eliminating spurious growing modes that otherwise induce numerical instabilities at small $Ek$. We show that the time derivative of the mean temperature is inconsequential for accurately determining the Nusselt number in the stationary state, significantly reducing the required simulation time, and demonstrate that full DNS using RRRiNSE agree with the NHQGE at very small $Ek$.

Type Ia Supernovae have yet again the opportunity to revolutionize the field of cosmology as the new generation of surveys are acquiring thousands of nearby SNeIa opening a new era in cosmology: the direct measurement of the growth of structure parametrized by $fD$. This method is based on the SNeIa peculiar velocities derived from the residual to the Hubble law as direct tracers of the full gravitational potential caused by large scale structure. With this technique, we could probe not only the properties of dark energy, but also the laws of gravity. In this paper we present the analytical framework and forecasts. We show that ZTF and LSST will be able to reach 5\% precision on $fD$ by 2027. Our analysis is not significantly sensitive to photo-typing, but known selection functions and spectroscopic redshifts are mandatory. We finally introduce an idea of a dedicated spectrograph that would get all the required information in addition to boost the efficiency to each SNeIa so that we could reach the 5\% precision within the first two years of LSST operation and the few percent level by the end of the survey.

The Landauer principle states that at least $k_B T \ln 2$ of energy is required to erase a 1-bit memory, with $k_B T$ the thermal energy of the system. We study the effects of inertia on this bound using as one-bit memory an underdamped micro-mechanical oscillator confined in a double-well potential created by a feedback loop. The potential barrier is precisely tunable in the few $k_B T$ range. We measure, within the stochastic thermodynamic framework, the work and the heat of the erasure protocol. We demonstrate experimentally and theoretically that, in this underdamped system, the Landauer bound is reached with a 1 % uncertainty, with protocols as short as 100 ms.

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a $1.4M_\odot$ neutron star to be $R=11.98^{+0.35}_{-0.40}$km.

Type Ia Supernovae have yet again the opportunity to revolutionize the field of cosmology as the new generation of surveys are acquiring thousands of nearby SNeIa opening a new era in cosmology: the direct measurement of the growth of structure parametrized by $fD$. This method is based on the SNeIa peculiar velocities derived from the residual to the Hubble law as direct tracers of the full gravitational potential caused by large scale structure. With this technique, we could probe not only the properties of dark energy, but also the laws of gravity. In this paper we present the analytical framework and forecasts. We show that ZTF and LSST will be able to reach 5\% precision on $fD$ by 2027. Our analysis is not significantly sensitive to photo-typing, but known selection functions and spectroscopic redshifts are mandatory. We finally introduce an idea of a dedicated spectrograph that would get all the required information in addition to boost the efficiency to each SNeIa so that we could reach the 5\% precision within the first two years of LSST operation and the few percent level by the end of the survey.

The control of particle trajectories in structured microfluidic environments has significantly advanced sorting technologies, most notably through deterministic lateral displacement (DLD). While previous work has largely targeted rigid, near-spherical particles, the sorting of flexible, anisotropic objects such as fibers remains largely unexplored. Here, we combine experiments and simulations to demonstrate how tilted pillar arrays enable efficient, length-based separation of flexible fibers. We discover that these arrays act as band-pass filters, selectively inducing lateral migration in fibers whose lengths are close to the array period. Fibers significantly shorter or longer exhibit minimal lateral deviation. This migration arises from the interplay of fluid-structure interactions between fibers and the complex flow and steric interactions with the pillars. Depending on their length, fibers exhibit distinct transport regimes: short fibers zigzag in between pillars following the flow, intermediate length fibers exhibit wrapping and jumping from one pillar to another, leading to lateral displacement, and long fibers deform extensively, following mixed zigzag-jump trajectories with minimal lateral migration. We identify the mechanical tension that develops in the fiber when wrapped around the pillars as the driving mechanism of cross-streamline transport. Leveraging this band-pass effect, we designed a highly efficient separation device to collect monodisperse fiber suspensions. Our findings not only expand the functional scope of DLD-like systems but also open new avenues for understanding transport of anisotropic objects in porous media.

We address the question of the role of low-energy nuclear physics data in constraining neutron star global properties, e.g., masses, radii, angular momentum, and tidal deformability, in the absence of a phase transition in dense matter. To do so, we assess the capacity of 415 relativistic mean field and non-relativistic Skyrme-type interactions to reproduce the ground state binding energies, the charge radii and the giant monopole resonances of a set of spherical nuclei. The interactions are classified according to their ability to describe these characteristics and we show that a tight correlation between the symmetry energy and its slope is obtained providing $N=Z$ and $N\ne Z$ nuclei are described with the same accuracy (mainly driven by the charge radius data). By additionally imposing the constraints from isobaric analog states and neutron skin radius in $^{208}$Pb, we obtain the following estimates: $E_{sym,2}=31.8\pm 0.7$ MeV and $L_{sym,2}=58.1\pm 9.0$ MeV. We then analyze predictions of neutron star properties and we find that the 1.4$M_\odot$ neutron star (NS) radius lies between 12 and 14 km for the "better" nuclear interactions. We show that i) the better reproduction of low-energy nuclear physics data by the nuclear models only weakly impacts the global properties of canonical mass neutron stars and ii) the experimental constraint on the symmetry energy is the most effective one for reducing the uncertainties in NS matter. However, since the density region where constraints are required are well above densities in finite nuclei, the largest uncertainty originates from the density dependence of the EDF, which remains largely unknown.

In the next decade, many optical surveys will aim to tackle the question of dark energy nature, measuring its equation of state parameter at the permil level. This requires trusting the photometric calibration of the survey with a precision never reached so far, controlling many sources of systematic uncertainties. The measurement of the on-site atmospheric transmission for each exposure, or on average for each season or for the full survey, can help reach the permil precision for magnitudes. This work aims at proving the ability to use slitless spectroscopy for standard star spectrophotometry and its use to monitor on-site atmospheric transmission as needed, for example, by the Vera C. Rubin Observatory Legacy Survey of Space and Time supernova cosmology program. We fully deal with the case of a disperser in the filter wheel, which is the configuration chosen in the Rubin Auxiliary Telescope. The theoretical basis of slitless spectrophotometry is at the heart of our forward model approach to extract spectroscopic information from slitless data. We developed a publicly available software called Spectractor (this https URL) that implements each ingredient of the model and finally performs a fit of a spectrogram model directly on image data to get the spectrum. We show on simulations that our model allows us to understand the structure of spectrophotometric exposures. We also demonstrate its use on real data, solving specific issues and illustrating how our procedure allows the improvement of the model describing the data. Finally, we discuss how this approach can be used to directly extract atmospheric transmission parameters from data and thus provide the base for on-site atmosphere monitoring. We show the efficiency of the procedure on simulations and test it on the limited data set available.

We present the Olsson$.$wl Mathematica package which aims to find linear transformations for some classes of multivariable hypergeometric functions. It is based on a well-known method developed by P. O. M. Olsson in J. Math. Phys. 5, 420 (1964) in order to derive the analytic continuations of the Appell $F_1$ double hypergeometric series from the linear transformations of the Gauss $_2F_1$ hypergeometric function. We provide a brief description of Olsson's method and demonstrate the commands of the package, along with examples. We also provide a companion package, called ROC2$.$wl and dedicated to the derivation of the regions of convergence of double hypergeometric series. This package can be used independently of Olsson$.$wl.

We present $griz$ photometric light curves for the full 5 years of the Dark Energy Survey Supernova program (DES-SN), obtained with both forced Point Spread Function (PSF) photometry on Difference Images (DIFFIMG) performed during survey operations, and Scene Modelling Photometry (SMP) on search images processed after the survey. This release contains $31,636$ DIFFIMG and $19,706$ high-quality SMP light curves, the latter of which contains $1635$ photometrically-classified supernovae that pass cosmology quality cuts. This sample spans the largest redshift ($z$) range ever covered by a single SN survey ($0.1

In the tetrad formulation of gravity, the so-called simplicity constraints play a central role. They appear in the Hamiltonian analysis of the theory, and in the Lagrangian path integral when constructing the gravity partition function from topological BF theory. We develop here a systematic analysis of the corner symplectic structure encoding the symmetry algebra of gravity, and perform a thorough analysis of the simplicity constraints. Starting from a precursor phase space with Poincaré and Heisenberg symmetry, we obtain the corner phase space of BF theory by imposing kinematical constraints. This amounts to fixing the Heisenberg frame with a choice of position and spin operators. The simplicity constraints then further reduce the Poincaré symmetry of the BF phase space to a Lorentz subalgebra. This picture provides a particle-like description of (quantum) geometry: The internal normal plays the role of the four-momentum, the Barbero-Immirzi parameter that of the mass, the flux that of a relativistic position, and the frame that of a spin harmonic oscillator. Moreover, we show that the corner area element corresponds to the Poincaré spin Casimir. We achieve this central result by properly splitting, in the continuum, the corner simplicity constraints into first and second class parts. We construct the complete set of Dirac observables, which includes the generators of the local $\mathfrak{sl}(2,\mathbb{C})$ subalgebra of Poincaré, and the components of the tangential corner metric satisfying an $\mathfrak{sl}(2,\mathbb{R})$ algebra. We then present a preliminary analysis of the covariant and continuous irreducible representations of the infinite-dimensional corner algebra. Moreover, as an alternative path to quantization, we also introduce a regularization of the corner algebra and interpret this discrete setting in terms of an extended notion of twisted geometries.

We summarize two geometrical approaches to analytically evaluate higher-fold Mellin-Barnes (MB) integrals in terms of hypergeometric functions. The first method is based on intersections of conic hulls, while the second one, which is more recent, relies on triangulations of a set of points. We demonstrate that, once automatized, the triangulation approach is computationally more efficient than the conic hull approach. As an application of this triangulation approach, we describe how one can derive simpler hypergeometric solutions of the conformal off-shell massless two-loop double box and one-loop hexagon Feynman integrals than those previously obtained from the conic hull approach. Lastly, by applying the above techniques on the MB representation of multiple polylogarithms, we show how to obtain new convergent series representations for these functions. These new analytic expressions were numerically cross-checked with GINAC.