Search

An automated framework is presented for the numerical solution of optimal control problems with PDEs as constraints, in both the stationary and instationary settings. The associated code can solve both linear and non-linear problems, and examples for incompressible flow equations are considered. The software, which is based on a Python interface to the Firedrake system, allows for a compact definition of the problem considered by providing a few lines of code in a high-level language. The software is provided with efficient iterative linear solvers for optimal control problems with PDEs as constraints. The use of advanced preconditioning techniques results in a significant speed-up of the solution process for large-scale problems. We present numerical examples of the applicability of the software on classical control problems with PDEs as constraints.

Valley photonic crystals provide efficient designs for the routing of light through channels in extremely compact geometries. The topological origin of the robust transport and the specific geometries under which it can take place have been questioned in recent works. In this article, we introduce a design for valley photonic crystals with richer arrangement possibilities than the standard valley photonic crystals based on two holes of different sizes in the unit cell. Our approach is based on the permutation of three sets of rhombi in an hexagonal lattice to investigate the interplay between Berry curvature, valley Chern number and chirality of interfaces to achieve robust edge-modes propagation along domain walls. We study three types of interfaces with different symmetries: the non-chiral interface with glide-mirror symmetry commonly used in honeycomb-type valley crystals, and two chiral interfaces with or without inversion symmetry of the adjacent bulk lattices. In the latter case, no valley topology is expected. We show that for the three families, edges preserving the shape of the interface through 120{\deg} sharp corners can sustain edge-modes with comparable robustness. Moreover, interfaces with glide-mirror symmetry offer promising performances in circuits with more exotic configurations, like 60{\deg} and 90{\deg} corners or arbitrary curves in which valley preservation is not guaranteed. Our work raises questions about the topological origin of the robustness of transport in valley photonic crystals, discusses the role of the chirality of the interfaces in the propagation around sharp corners, and provides a lattice scheme with broad design possibilities., Comment: 2 pages, 6 figures, 1 supplemental materials file, 2 supplemental materials movies

We address the solution of the distributed control problem for the steady, incompressible Navier--Stokes equations. We propose an inexact Newton linearization of the optimality conditions. Upon discretization by a finite element scheme, we obtain a sequence of large symmetric linear systems of saddle-point type. We use an augmented Lagrangian-based block triangular preconditioner in combination with the flexible GMRES method at each Newton step. The preconditioner is applied inexactly via a suitable multigrid solver. Numerical experiments indicate that the resulting method appears to be fairly robust with respect to viscosity, mesh size, and the choice of regularization parameter when applied to 2D problems.

AI assistants are being increasingly used by students enrolled in higher education institutions. While these tools provide opportunities for improved teaching and education, they also pose significant challenges for assessment and learning outcomes. We conceptualize these challenges through the lens of vulnerability, the potential for university assessments and learning outcomes to be impacted by student use of generative AI. We investigate the potential scale of this vulnerability by measuring the degree to which AI assistants can complete assessment questions in standard university-level STEM courses. Specifically, we compile a novel dataset of textual assessment questions from 50 courses at EPFL and evaluate whether two AI assistants, GPT-3.5 and GPT-4 can adequately answer these questions. We use eight prompting strategies to produce responses and find that GPT-4 answers an average of 65.8% of questions correctly, and can even produce the correct answer across at least one prompting strategy for 85.1% of questions. When grouping courses in our dataset by degree program, these systems already pass non-project assessments of large numbers of core courses in various degree programs, posing risks to higher education accreditation that will be amplified as these models improve. Our results call for revising program-level assessment design in higher education in light of advances in generative AI., Comment: 20 pages, 8 figures

We show that the complete photoemission dynamics in situations of electron-ion entanglement can be retrieved from photoelectron spectral measurements without information on the ion. To this end, we develop an energy-time analysis of the photoelectron's reduced density matrix based on first principles. We test and assess our approach with numerical simulations on a low dimensional model molecule in interaction with broadband composite pulses occulting the vibrational resolution. Our method is directly applicable to recent experimental schemes measuring the photoelectron reduced density matrices in atomic and molecular photoemission. Therefore, it opens a new window on the dynamics of decoherence and entanglement at the attosecond timescale., Comment: 7 pages, 2 figures

Although the sequential tsunami scenario detection framework was validated in our previous work, several tasks remain to be resolved from a practical point of view. This study aims to evaluate the performance of the previous tsunami scenario detection framework using a diverse database consisting of complex fault rupture patterns with heterogeneous slip distributions. Specifically, we compare the effectiveness of scenario superposition to that of the previous most likely scenario detection method. Additionally, how the length of the observation time window influences the accuracy of both methods is analyzed. We utilize an existing database comprising 1771 tsunami scenarios targeting the city Westport (WA, U.S.), which includes synthetic wave height records and inundation distributions as the result of fault rupture in the Cascadia subduction zone. The heterogeneous patterns of slips used in the database increase the diversity of the scenarios and thus make it a proper database for evaluating the performance of scenario superposition. To assess the performance, we consider various observation time windows shorter than 15 minutes and divide the database into five testing and learning sets. The evaluation accuracy of the maximum offshore wave, inundation depth, and its distribution is analyzed to examine the advantages of the scenario superposition method over the previous method. We introduce the dynamic time warping (DTW) method as an additional benchmark and compare its results to that of the Bayesian scenario detection method., Comment: 15 pages, 12 figures

We consider the two-photon ionization of Hydrogen-like atoms. We find an approximate expression of the long-range phase based on an asymptotic expansion of the continuum eigenfunctions within the Wentzel-Kramers-Brillouin approximation. Combined with commonly used perturbative approaches, the resulting analytic formalism can treat, at the same time, the two-photon propensity rules, the anisotropy in the continuum-continuum photoionization time delay and the soft-photon regime.

We investigate ultrafast vibronic dynamics triggered by intense femtosecond infrared pulses in small molecules. Our study is based on numerical simulations performed with 2D model molecules, and analyzed in the perspective of the renown Lochfrass and Bond-Softening models. We give a new interpretation of the observed nuclear wave packet dynamics, with a focus on the phase of the bond oscillations. Our simulations also reveal intricate features in the field-induced nuclear motion that are not accounted for by existing models. Our analyses assign these features to strong dynamic correlations between the active electron and the nuclei, which significantly depend on the carrier envelope phase of the pulse, even for relatively ``long'' pulses, which should make them experimentally observable., Comment: 34 pages, 14 figures

We investigate how vibration affects molecular photoemission dynamics, through simulations on two-dimension asymmetric model molecules including the electronic and nuclear motions in a fully correlated way. We show that a slight anisotropy in the electron-ion momentum sharing is sufficient to prevent one from unambigously characterizing the vibrationnaly averaged photoemission dynamics in terms of stereo Wigner delays. We further show that vibrational resolution can be retrieved in fixed-nuclei simulations, using effective molecular conformations that are specific to each vibrational channel. The optimal internuclear distances found empirically in 1-photon processes can be identified a priori using simple physical arguments. They also turn out to be efficient to simulate vibrationnally-resolved \rabbit measurements and to account for interchannel coherences in broadband 1-photon ionization., Comment: 17 pages, 13 figures

We present the results of an optical link to a corner cube on board a tethered balloon at 300 m altitude including a Tip/Tilt compensation for the balloon tracking. Our experiment measures the carrier phase of a 1542 nm laser, which is the useful signal for frequency comparison of distant clocks. An active phase noise compensation of the carrier is implemented, demonstrating a fractional frequency stability of 8x10-19 after 16 s averaging, which slightly (factor ~3) improves on best previous links via an airborne platform. This state-of-the-art result is obtained with a transportable set-up that enables a fast field deployment., Comment: accepted in Optics Express

We investigate signatures of anisotropy on the dynamics of time-resolved near-threshold molecular photoemission, through simulations on a one-dimension asymmetric model molecule. More precisely, we study the relationship between the fundamental Wigner delays that fully characterizes the dynamics of one-photon ionization, and the delays inferred from two-color interferometric rabbit measurements. Our results highlights two different properties pertaining to each of these delays. The first one is related to the inherent necessity to fix an arbitrary electron position origin to define and compute de Wigner-delay. The second one is the dependency of the rabbit delay on the frequency of the probe laser. Our results show that the angular variations of both delays converge for a specific choice of the position origin and in the limit of a vanishing probe frequency., Comment: 13 pages, 9 figures

We present cross sections for the formation of positronium chloride (PsCl) in its ground state from the charge exchange between positronium (Ps) and chloride (Cl$^-$) in the range of 10 meV - 100 eV Ps energy. We have used theoretical models based on the first Born approximation in its three-body formulation. We simulated the collisions between Ps and Cl$^-$ using ab-initio methods at both mean-field and correlated levels extrapolated to the complete basis set limit. We have investigated Ps excited states up to ${n=4}$. The results suggest that the channel Ps(${n=2}$) is of particular interest for the production of PsCl in the ground state, and shows that an accurate treatment of the electronic correlation leads to a significant change in the global shape of the PsCl production cross section with respect to the mean-field level., Comment: 13 Pages, 7 Figures, 3 Tables

We present an algorithm to solve the dispersive depth-averaged Serre-Green-Naghdi (SGN) equations using patch-based adaptive mesh refinement. These equations require adding additional higher derivative terms to the nonlinear shallow water equations. This has been implemented as a new component of the open source GeoClaw software that is widely used for modeling tsunamis, storm surge, and related hazards, improving its accuracy on shorter wavelength phenomena. We use a formulation that requires solving an elliptic system of equations at each time step, making the method implicit. The adaptive algorithm allows different time steps on different refinement levels, and solves the implicit equations level by level. Computational examples are presented to illustrate the stability and accuracy on a radially symmetric test case and two realistic tsunami modeling problems, including a hypothetical asteroid impact creating a short wavelength tsunami for which dispersive terms are necessary., Comment: 25 pages, 16 figures, submitted for publication (revised)

We demonstrate that strong-field ionization of atoms driven by circularly polarized light becomes an adiabatic process when described in the frame rotating with the laser field. As a direct consequence, a conservation law emerges: in the rotating frame the energy of the tunneling electron is conserved for rotationally invariant potentials. This conservation law, arising from a classical picture, is retrieved through a proper classical-quantum correspondence when considering the full quantum system, beyond the Strong Field Approximation.

Important applications of single-particle tracking (SPT) aim at deciphering the diffusion properties of single fluorescent nanoparticles immersed in heterogeneous environments, such as multi-cellular biological tissues. To maximize the particle localization precision in such complex environments, high numerical aperture objectives are often required, which intrinsically restrict depth-of-focus (DOF) to less than a micrometer and impedes recording long trajectories when particles escape the plane of focus. In this work, we show that a simple binary phase mask can work with the spherical aberration inevitably induced by thick sample inhomogeneities, to extend the DOF of a single-molecule fluorescence microscope over more than 4 {\mu}m. The effect of point-spread-function (PSF) engineering over spherical aberration regularizes inhomogeneities of the PSF along the optical axis by restricting it to a narrow distribution. This allows the use of a single fitting function (i.e. Gaussian function) to localize single emitters over the whole extended DOF. Application of this simple approach on diffusing nanoparticles demonstrate that SPT trajectories can be recorded on significantly longer times., Comment: 18 pages, 4 figures

The Hall-Vinen-Bekharevich-Khalatnikov (HVBK) model is widely used to numerically study quantum turbulence in superfluid helium. Based on the two-fluid model of Tisza and Landau, the HVBK model describes the normal (viscous) and superfluid (inviscid) components of the flow using two Navier-Stokes type equations, coupled through a mutual friction force term. This feature makes the HVBK model very appealing in applying statistical tools used in classical turbulence to study properties of quantum turbulence. A large body of literature used low-order statistics (spectra, or second-order structure functions in real space) to unravel exchanges between the two fluids at several levels. The novelty in this study is to use a theoretical approach based on first principles to derive transport equations for the third-order moments for each component of velocity., Comment: 27 pages, 10 figures, accepted for publication in Journal of Fluid Mechanics

In this article, we derive fast and robust parallel-in-time preconditioned iterative methods for the all-at-once linear systems arising upon discretization of time-dependent PDEs. The discretization we employ is based on a Runge--Kutta method in time, for which the development of parallel solvers is an emerging research area in the literature of numerical methods for time-dependent PDEs. By making use of classical theory of block matrices, one is able to derive a preconditioner for the systems considered. The block structure of the preconditioner allows for parallelism in the time variable, as long as one is able to provide an optimal solver for the system of the stages of the method. We thus propose a preconditioner for the latter system based on a singular value decomposition (SVD) of the (real) Runge--Kutta matrix $A_{\mathrm{RK}} = U \Sigma V^\top$. Supposing $A_{\mathrm{RK}}$ is invertible, we prove that the spectrum of the system for the stages preconditioned by our SVD-based preconditioner is contained within the right-half of the unit circle, under suitable assumptions on the matrix $U^\top V$ (the assumptions are well posed due to the polar decomposition of $A_{\mathrm{RK}}$). We show the numerical efficiency of our SVD-based preconditioner by solving the system of the stages arising from the discretization of the heat equation and the Stokes equations, with sequential time-stepping. Finally, we provide numerical results of the all-at-once approach for both problems, showing the speed-up achieved on a parallel architecture.

The linear scaling of the mixed third-order moment of the magnetohydrodynamic fluctuations is used to estimate the energy transfer rate of the turbulent cascade in the expanding solar wind. In 1976 the Helios 2 spacecraft measured three samples of fast solar wind originating from the same coronal hole, at different distance from the sun. Along with the adjacent slow solar wind streams, these represent a unique database for studying the radial evolution of turbulence in samples of undisturbed solar wind. A set of direct numerical simulations of the MHD equations performed with the Lattice-Boltzmann code FLAME is also used for interpretation. We show that the turbulence energy transfer rate decays approximately as a power law of the distance, and that both the amplitude and decay law correspond to the observed radial temperature profile in the fast wind case. Results from magnetohydrodynamic numerical simulations of decaying magnetohydrodynamic turbulence show a similar trend for the total dissipation, suggesting an interpretation of the observed dynamics in terms of decaying turbulence, and that multi-spacecraft studies of the solar wind radial evolution may help clarifying the nature of the evolution of the turbulent fluctuations in the ecliptic solar wind., Comment: In press on Astron. Astrophys

In this work, we use valley-topological triangular resonators coupled to an input waveguide to evaluate the quality of the topological protection. To that purpose, we first analyze via numerical simulations the existence of backward scattering at cavity corners or transmission with pseudo-spin conversion at the splitter between the input waveguide and the cavity. We evidence that a breakdown of topological protection takes place, in particular at sharp corners, which results in transmission minima and split-resonances, otherwise non-existent. In order to evaluate the small coupling coefficients associated to this breakdown, a phenomenological model based on an exact parameterization of scattering matrices at splitters and corners of the resonators is then introduced. By comparison with the numerical simulations, we are able to quantify the loss of topological protection at sharp bends and splitters. Finally, we use the obtained set of phenomenological parameters to compare the predictions of the phenomenological model with full numerical simulations for fractal-inspired cavities based on the Sierpi\'nski triangle construction. We show that the agreement is overall good, but shows more differences for the cavity composed of the smallest triangles. Our results suggest that even in a system exempt of geometrical and structural defects, topological protection is not complete at corners, sharp bends and splitters. However, simpler but predictive calculations can be realized with a phenomenological approach, allowing simulations of very large devices beyond the reach of standard simulation methods, which is crucial to design photonic devices which gather compactness and low losses through topological conduction of electromagnetic waves., Comment: 24 pages, 7 figures, one supplementary informations file

Simulating plasmas in the Hall-MagnetoHydroDynamics (Hall-MHD) regime represents a valuable {approach for the investigation of} complex non-linear dynamics developing in astrophysical {frameworks} and {fusion machines}. Taking into account the Hall electric field is {computationally very challenging as} it involves {the integration of} an additional term, proportional to $\bNabla \times ((\bNabla\times\mathbf{B})\times \mathbf{B})$ in the Faraday's induction {law}. {The latter feeds back on} the magnetic field $\mathbf{B}$ at small scales (between the ion and electron inertial scales), {requiring} very high resolution{s} in both space and time {in order to properly describe its dynamics.} The computational {advantage provided by the} kinetic Lattice Boltzmann (LB) approach is {exploited here to develop a new} code, the \textbf{\textsc{F}}ast \textbf{\textsc{L}}attice-Boltzmann \textbf{\textsc{A}}lgorithm for \textbf{\textsc{M}}hd \textbf{\textsc{E}}xperiments (\textsc{flame}). The \textsc{flame} code integrates the plasma dynamics in lattice units coupling two kinetic schemes, one for the fluid protons (including the Lorentz force), the other to solve the induction equation describing the evolution of the magnetic field. Here, the newly developed algorithm is tested against an analytical wave-solution of the dissipative Hall-MHD equations, pointing out its stability and second-order convergence, over a wide range of the control parameters. Spectral properties of the simulated plasma are finally compared with those obtained from numerical solutions from the well-established pseudo-spectral code \textsc{ghost}. Furthermore, the LB simulations we present, varying the Hall parameter, highlightthe transition from the MHD to the Hall-MHD regime, in excellent agreement with the magnetic field spectra measured in the solar wind.

Quality assessment is a key element for the evaluation of hardware and software involved in image and video acquisition, processing, and visualization. In the medical field, user-based quality assessment is still considered more reliable than objective methods, which allow the implementation of automated and more efficient solutions. Regardless of increasing research in this topic in the last decade, defining quality standards for medical content remains a non-trivial task, as the focus should be on the diagnostic value assessed from expert viewers rather than the perceived quality from na\"{i}ve viewers, and objective quality metrics should aim at estimating the first rather than the latter. In this paper, we present a survey of methodologies used for the objective quality assessment of medical images and videos, dividing them into visual quality-based and task-based approaches. Visual quality based methods compute a quality index directly from visual attributes, while task-based methods, being increasingly explored, measure the impact of quality impairments on the performance of a specific task. A discussion on the limitations of state-of-the-art research on this topic is also provided, along with future challenges to be addressed., Comment: Submitted for peer review at Multimedia Tools and Applications

We study the effects of rotations on a cold atom accelerometer onboard a Nadir pointing satellite. A simulation of the satellite attitude combined with a calculation of the phase of the cold atom interferometer allow us to evaluate the noise and bias induced by rotations. In particular, we evaluate the effects associated to the active compensation of the rotation due to Nadir pointing. This study was realized in the context of the preliminary study phase of the CARIOQA Quantum Pathfinder Mission.

A strong potential gain for space applications is expected from the anticipated performances of inertial sensors based on cold atom interferometry (CAI) that measure the acceleration of freely falling independent atoms by manipulating them with laser light. In this context, CNES and its partners initiated a phase 0 study, called CARIOQA, in order to develop a Quantum Pathfinder Mission unlocking key features of atom interferometry for space and paving the way for future ambitious space missions utilizing this technology. As a cornerstone for the implementation of quantum sensors in space, the CARIOQA phase 0 aimed at defining the Quantum Pathfinder Mission's scenario and associated performance objectives. To comply with these objectives, the payload architecture has been designed to achieve long interrogation time and active rotation compensation on a BEC-based atom interferometer. A study of the satellite architecture, including all the subsystems, has been conducted. Several technical solutions for propulsion and attitude control have been investigated in order to guarantee optimal operating conditions (limitation of micro-vibrations, maximization of measurement time). A preliminary design of the satellite platform was performed., Comment: Proceedings of International Conference on Space Optics (ICSO) 2022; 3-7 October 2022; Dubrovnik; Croatia

In this work, we investigate the black hole (BH) population of globular clusters (GCs) in Milky Way- (MW) and Andromeda- (M31) like galaxies. We combine the population synthesis code MASinGa and the MOCCA-Survey Database I to infer the properties of GCs harbouring a BH subsystem (BHS), an IMBH, or neither of those. We find that the typical number of GCs with a BHS, an IMBH, or none become comparable in the galactic outskirts, whilst the inner galactic regions are dominated by GCs without a significant dark component. Our models suggest that GCs harbouring a BHS are slightly heavier and with larger half-mass radii compared to the overall population. We retrieve the properties of binary BHs (BBHs) that have either merged in the last 3 Gyr or survived in their parent cluster until present-day. We find that around 80\% of the merging BBHs form due to dynamical interactions while the remaining originate from evolution of primordial binaries. We infer a merger rate for BBHs in the local Universe of $1.0-23\,\,\rm{yr^{-1}\,Gpc^{-3}}$, depending on the adopted assumptions. We find around 100-240 BBHs survive until present-day and are mostly concentrated in the inner few kpc of the galaxy. We estimate also the number of BHs transported into the galactic nucleus by infalling star clusters, finding around 1,000-3,000 BHs and 100-200 BBHs are transported over a time span of 12 Gyr. This enables us to constrains the total amount of BHs and BBHs binaries lurking in nuclear star cluster, i.e. $N_{BHs}=(1.4-2.2)\times10^4$ and $N_{BBHs}=700-1,100$., Comment: 18 pages,15 figures

A comprehensive study of the co-evolution of globular cluster systems (GCS) in galaxies requires the ability to model both the large scale dynamics (0.01 - 10 kpc) regulating their orbital evolution, and the small scale dynamics (sub-pc - AU) regulating the internal dynamics of each globular cluster (GC). In this work we present a novel method that combine semi-analytic models of GCS with fully self-consistent Monte Carlo models to simultaneously evolve large GCSs. We use the population synthesis code MASinGa and the MOCCA-Survey Database I to create synthetic GC populations aimed at representing the observed features of GCs in the Milky Way (MW) and Andromeda (M31). Our procedure enables us to recover the spatial and mass distribution of GCs in such galaxies, and to constrain the amount of mass that GCs left either in the halo as dispersed debris, or in the galactic centre, where they can contribute to the formation of a nuclear star cluster (NSC) and can bring stellar and possibly intermediate mass black holes there. The final masses reported by our simulations are of a few order of magnitudes smaller than the observed values. These differences show that mass build-up of a NSC and central BHs in galaxies like MW and M31 cannot be solely explained by the infalling GC scenario. This build-up is likely to depend on the interplay between interactions and mergers of infalling GCs and gas. The latter can contribute to both in-situ star formation in the NSC and growth of the central BH., Comment: 17 pages, 10 figures

We introduce a new prescription for the evolution of globular clusters (GCs) during the initial embedded gas phase into a Monte Carlo method. With a simplified version of the Monte Carlo MOCCA code embedded in the AMUSE framework, we study the survival of GCs after the removal of primordial gas. We first test our code and show that our results for the evolution of mass and Lagrangian radii are in good agreement with those obtained with N-body simulations. The Monte Carlo code enables a more rapid exploration of the evolution of systems with a larger number of stars than N-body simulations. We have carried out a new survey of simulations to explore the evolution of globular clusters with up to $N = 500000$ stars for a range of different star formation efficiencies and half-mass radii. Our study shows the range of initial conditions leading to the clusters' dissolution and those for which the clusters can survive this early evolutionary phase., Comment: 13 pages, 6 figures

For a fixed graph H, the H-Recoloring problem asks whether for two given homomorphisms from a graph G to H, we can transform one into the other by changing the image of a single vertex of G in each step and maintaining a homomorphism from G to H throughout. We extend an algorithm of Wrochna for H-Recoloring where H is a square-free loopless undirected graph to the more general setting of directed graphs. We obtain a polynomial-time algorithm for H-Recoloring in this setting whenever H is a loopless digraph that does not contain a 4-cycle of algebraic girth zero and whenever H is a reflexive digraph that contains neither a 3-cycle of algebraic girth 1 nor a 4-cycle of algebraic girth zero.

We present an upgraded version of the \MOCCA code for the study of dynamical evolution of globular clusters (GCs) and its first application to the study of evolution of multiple stellar populations. We explore initial conditions spanning different structural parameters for the first (FG) and second generation of stars (SG) and we analyze their effect on the binary dynamics and survival. Here, we focus on the number ratio of FG and SG binaries, its spatial variation, and the way their abundances are affected by various cluster initial properties. We find that present-day SG stars are more abundant in clusters that were initially tidally filling. Conversely, FG stars stay more abundant in clusters that were initially tidally underfilling. We find that the ratio between binary fractions is not affected by the way we calculate these fractions (e.g. only main-sequence binaries (MS) or observational binaries, i.e. MS stars $> 0.4 M_{\odot}$ mass ratios $> 0.5$). This implies that the MS stars themselves are a very good proxy for probing entire populations of FG and SG. We also discuss how it relates to the observations of Milky Way GCs. We show that \MOCCA models are able to reproduce the observed range of SG fractions for Milky Way GCs for which we know these fractions. We show how the SG fractions depend on the initial conditions and provide some constraints for the initial conditions to have more numerous FG or SG stars at the Hubble time., Comment: Acceped to MNRAS

We present a description of the GBAR positron (e+) trapping apparatus, which consists of a three stage Buffer Gas Trap (BGT) followed by a High Field Penning Trap (HFT), and discuss its performance. The overall goal of the GBAR experiment is to measure the acceleration of the neutral antihydrogen (H) atom in the terrestrial gravitational field by neutralising a positive antihydrogen ion (H+), which has been cooled to a low temperature, and observing the subsequent H annihilation following free fall. To produce one H+ ion, about 10^10 positrons, efficiently converted into positronium (Ps), together with about 10^7 antiprotons (p), are required. The positrons, produced from an electron linac-based system, are accumulated first in the BGT whereafter they are stacked in the ultra-high vacuum HFT, where we have been able to trap 1.4(2) x 10^9 positrons in 1100 seconds.

Structured illumination in Single Molecule Localization Microscopy provides new information on the position of molecules and thus improves the localization precision compared to standard localization methods. Here, we used a time-shifted sinusoidal excitation pattern to modulate the fluorescence signal of the molecules whose position information is carried by the phase and recovered by synchronous demodulation. We designed two flexible fast demodulation systems located upstream of the camera, allowing us to overcome the limiting camera acquisition frequency and thus to maximize the collection of photons in the demodulation process. The temporally modulated fluorescence signal was then sampled synchronously on the same image, repeatedly during acquisition. This microscopy, called ModLoc, allows to experimentally improve the localization precision by a factor of 2.4 in one direction, compared to classical Gaussian fitting methods. A temporal study and an experimental demonstration both show that the short lifetimes of the molecules in blinking regimes impose a modulation frequency in the kilohertz range, which is beyond the reach of current cameras. A demodulation system operating at these frequencies would thus be necessary to take full advantage of this new localization approach.

We summarize the discussions at a virtual Community Workshop on Cold Atoms in Space concerning the status of cold atom technologies, the prospective scientific and societal opportunities offered by their deployment in space, and the developments needed before cold atoms could be operated in space. The cold atom technologies discussed include atomic clocks, quantum gravimeters and accelerometers, and atom interferometers. Prospective applications include metrology, geodesy and measurement of terrestrial mass change due to, e.g., climate change, and fundamental science experiments such as tests of the equivalence principle, searches for dark matter, measurements of gravitational waves and tests of quantum mechanics. We review the current status of cold atom technologies and outline the requirements for their space qualification, including the development paths and the corresponding technical milestones, and identifying possible pathfinder missions to pave the way for missions to exploit the full potential of cold atoms in space. Finally, we present a first draft of a possible road-map for achieving these goals, that we propose for discussion by the interested cold atom, Earth Observation, fundamental physics and other prospective scientific user communities, together with ESA and national space and research funding agencies., Comment: Summary of the Community Workshop on Cold Atoms in Space and corresponding Road-map: https://indico.cern.ch/event/1064855/

The long-term goal of this work is the development of high-fidelity simulation tools for dispersive tsunami propagation. A dispersive model is especially important for short wavelength phenomena such as an asteroid impact into the ocean, and is also important in modeling other events where the simpler shallow water equations are insufficient. Adaptive simulations are crucial to bridge the scales from deep ocean to inundation, but have difficulties with the implicit system of equations that results from dispersive models. We propose a fractional step scheme that advances the solution on separate patches with different spatial resolutions and time steps. We show a simulation with 7 levels of adaptive meshes and onshore inundation resulting from a simulated asteroid impact off the coast of Washington. Finally, we discuss a number of open research questions that need to be resolved for high quality simulations., Comment: 16 pages, 5 figures, submitted to Proc. International Congress of Mathematicians, 2022

A dominating set in a directed graph is a set of vertices $S$ such that all the vertices that do not belong to $S$ have an in-neighbour in $S$. A locating set $S$ is a set of vertices such that all the vertices that do not belong to $S$ are characterized uniquely by the in-neighbours they have in $S$, i.e. for every two vertices $u$ and $v$ that are not in $S$, there exists a vertex $s\in S$ that dominates exactly one of them. The size of a smallest set of a directed graph $D$ which is both locating and dominating is denoted by $\gamma^{LD}(D)$. Foucaud, Heydarshahi and Parreau proved that any twin-free digraph $D$ satisfies $\gamma^{LD}(D)\leq \frac{4n} 5 +1$ but conjectured that this bound can be lowered to $\frac{2n} 3$. The conjecture is still open. They also proved that if $D$ is a tournament, i.e. a directed graph where there is one arc between every pair of vertices, then $\gamma^{LD}(D)\leq \lceil \frac{n}{2}\rceil$. The main result of this paper is the generalization of this bound to connected local tournaments, i.e. connected digraphs where the in- and out-neighbourhoods of every vertex induce a tournament. We also prove $\gamma^{LD}(D)\leq \frac{2n} 3$ for all quasi-twin-free digraphs $D$ that admit a supervising vertex (a vertex from which any vertex is reachable). This class of digraphs generalizes twin-free acyclic graphs, the most general class for which this bound was known.

We derive novel, fast, and parameter-robust preconditioned iterative methods for steady and time-dependent Navier--Stokes control problems. Our approach may be applied to time-dependent problems which are discretized using backward Euler or Crank--Nicolson, and is also a valuable candidate for Stokes control problems discretized using Crank--Nicolson. The key ingredients of the solver are a saddle-point type approximation for the linear systems, an inner iteration for the $(1,1)$-block accelerated by a preconditioner for convection--diffusion control, and an approximation to the Schur complement based on a potent commutator argument applied to an appropriate block matrix. A range of numerical experiments validate the effectiveness of our new approach.

We study a one-dimensional array of bosons with infinite-range interactions mediated by a laser-driven dissipative optical cavity. The cavity-mediated infinite-range interactions open up a new pathway to fermionization, hitherto only known for dipolar bosons due to their long-range interactions. In parameter ranges attainable in state-of-the-art experiments, we systematically compare observables for bosons and fermions with infinite-range interactions. At large enough laser pump power, many observables, including density distributions in real and momentum space, correlation functions, eigenvalues of the one-body density matrix, and superradiance order parameter, become identical for bosons and fermions. We map out the emergence of this cavity-induced fermionization as a function of pump power and contact interactions. We discover that cavity-mediated interactions can compensate a reduction by several orders of magnitude in the strength of the contact interactions needed to trigger fermionization., Comment: main text: 7 pages, 3 figures - supplement: 7 pages, 4 figures

Non-Gaussian statistics of large-scale fields are routinely observed in data from atmospheric and oceanic campaigns and global models. Recent direct numerical simulations (DNSs) showed that large-scale intermittency in stably stratified flows is due to the emergence of sporadic, extreme events in the form of bursts in the vertical velocity and the temperature. This phenomenon results from the interplay between waves and turbulent motions, affecting mixing. We provide evidence of the enhancement of the classical small-scale (or internal) intermittency due to the emergence of large-scale drafts, connecting large- and small-scale bursts. To this aim we analyze a large set of DNSs of the stably stratified Boussinesq equations over a wide range of values of the Froude number ($Fr\approx 0.01-1$). The variation of the buoyancy field kurtosis with $Fr$ is similar to (though with smaller values than) the kurtosis of the vertical velocity, both showing a non-monotonic trend. We present a mechanism for the generation of extreme vertical drafts and vorticity enhancements which follows from the exact equations for field gradients., Comment: 7 pages, 4 figures

We present the implementation of updated stellar evolution recipes in the codes \texttt{Nbody6++GPU, MOCCA} and \texttt{McLuster}. We test them through numerical simulations of star clusters containing $1.1\times 10^5$ stars (with $2.0\times 10^4$ in primordial hard binaries) performing high-resolution direct $N$-body (\texttt{Nbody6++GPU}) and Monte-Carlo (\texttt{MOCCA}) simulations to an age of 10~Gyr. We compare models implementing either delayed or core-collapse supernovae mechanisms, a different mass ratio distribution for binaries, and white dwarf natal kicks enabled/disabled. Compared to \texttt{Nbody6++GPU}, the \texttt{MOCCA} models appear to be denser, with a larger scatter in the remnant masses, and a lower binary fraction on average. The \texttt{MOCCA} models produce more black holes (BHs) and helium white dwarfs (WDs), whilst \texttt{Nbody6++GPU} models are characterised by a much larger amount of WD-WD binaries. The remnant kick velocity and escape speed distributions are similar for the BHs and neutron stars (NSs), and some NSs formed via electron-capture supernovae, accretion-induced collapse or merger-induced collapse escape the cluster in all simulations. The escape speed distributions for the WDs, on the other hand, are very dissimilar. We categorise the stellar evolution recipes available in \texttt{Nbody6++GPU}, \texttt{MOCCA} and \texttt{Mcluster} into four levels: the one implemented in previous \texttt{Nbody6++GPU} and \texttt{MOCCA} versions (\texttt{level A}), state-of-the-art prescriptions (\texttt{level B}), some in a testing phase (\texttt{level C}), and those that will be added in future versions of our codes., Comment: 31 pages, 15 figures, 5 tables, Accepted by MNRAS

A proof labelling scheme for a graph class $\mathcal{C}$ is an assignment of certificates to the vertices of any graph in the class $\mathcal{C}$, such that upon reading its certificate and the certificates of its neighbors, every vertex from a graph $G\in \mathcal{C}$ accepts the instance, while if $G\not\in \mathcal{C}$, for every possible assignment of certificates, at least one vertex rejects the instance. It was proved recently that for any fixed surface $\Sigma$, the class of graphs embeddable in $\Sigma$ has a proof labelling scheme in which each vertex of an $n$-vertex graph receives a certificate of at most $O(\log n)$ bits. The proof is quite long and intricate and heavily relies on an earlier result for planar graphs. Here we give a very short proof for any surface. The main idea is to encode a rotation system locally, together with a spanning tree supporting the local computation of the genus via Euler's formula., Comment: 10 pages, 5 figures - v4: final version

We consider the dynamics of two-dimensional interacting ultracold bosons triggered by suddenly switching on an artificial gauge field. The system is initialized in the ground state of a harmonic trapping potential. As a function of the strength of the applied artificial gauge, we analyze the emergent dynamics by monitoring the angular momentum, the fragmentation as well the entropy and variance of the entropy of absorption or single-shot images. We solve the underlying time-dependent many-boson Schr\"odinger equation using the multiconfigurational time-dependent Hartree method for indistinguishable particles (MCTDH-X). We find that the artificial gauge field implants angular momentum in the system. Fragmentation -- multiple macroscopic eigenvalues of the reduced one-body density matrix -- emerges in sync with the dynamics of angular momentum: the bosons in the many-body state develop non-trivial correlations. Fragmentation and angular momentum are experimentally difficult to assess; here, we demonstrate that they can be probed by statistically analyzing the variance of the image entropy of single-shot images that are the standard projective measurement of the state of ultracold atomic systems., Comment: 14 pages, 6 figures, MCTDH-X software hosted at http://ultracold.org

Quasi all water soluble composites use graphene oxide (GO) or reduced graphene oxide (rGO) as graphene based additives despite the long and harsh conditions required for their preparation. Herein, polyvinyl alcohol (PVA) films containing few layer graphene (FLG) are prepared by the co-mixing of aqueous colloids and casting, where the FLG colloid is first obtained via an effcient, rapid, simple, and bio-compatible exfoliation method providing access to relatively large FLG flakes. The enhanced mechanical, electrical conductivity, and O2 barrier properties of the films are investigated and discussed together with the structure of the films. In four different series of the composites, the best Young modulus is measured for the films containing around 1% of FLG. The most significant enhancement is obtained for the series with the largest FLG sheets contrary to the elongation at break which is well improved for the series with the lowest FLG sheets. Relatively high one-side electrical conductivity and low percolation threshold are achieved when compared to GO(rGO) composites (almost 10-3 S.cm-1 for 3% of FLG and transport at 0.5% FLG), while the conductivity is affected by the formation of a macroscopic branched FLG network. The composites demonstrate a reduction of O2 transmission rate up to 60%.

The emergence of quantum technologies, including cold atom based accelerometers, offers an opportunity to improve the performances of space geodesy missions. In this context, CNES initiated an assessment study called GRICE (GRadiom\'etrie \`a Interf\'erom\`etres quantiques Corr\'el\'es pour l'Espace) in order to evaluate the contribution of cold atom technologies to space geodesy and to the end users of geodetic data. In this paper, we present mission scenario for gravity field mapping based on a long baseline gradiometer. The mission is based on a constellation of two satellites, flying at an altitude of 373 km, each equipped with a cold atom accelerometer with a sensitivity of $6 \times 10^{-10}$~m.s$^{-2}$.$\mathrm{\tau}^{-1/2}$. A laser link measures the distance between the two satellites and couples these two instruments in order to produce a correlated differential acceleration measurement. The main parameters, determining the performances of the payload, have been investigated. We carried out a general study of satellite architecture and simulations of the mission performances in terms of restitution of the gravity field. The simulations show that this concept would give its best performance in terms of monthly gravity fields recovery under 1000~km resolution. In the resolution band between 1000 and 222~km, the improvement of the GRICE gradient approach over the traditional range-rate approach is globally in the order of 10 to 25\%.

Single-shot images are the standard readout of experiments with ultracold atoms -- the tarnished looking glass into their many-body physics. The efficient extraction of observables from single-shot images is thus crucial. Here, we demonstrate how artificial neural networks can optimize this extraction. In contrast to standard averaging approaches, machine learning allows both one- and two-particle densities to be accurately obtained from a drastically reduced number of single-shot images. Quantum fluctuations and correlations are directly harnessed to obtain physical observables for bosons in a tilted double-well potential at an unprecedented accuracy. Strikingly, machine learning also enables a reliable extraction of momentum-space observables from real-space single-shot images and vice versa. This obviates the need for a reconfiguration of the experimental setup between in-situ and time-of-flight imaging, thus potentially granting an outstanding reduction in resources., Comment: 7+8 pages, 3+8 figures, software available at http://ultracold.org

Research on image quality assessment (IQA) remains limited mainly due to our incomplete knowledge about human visual perception. Existing IQA algorithms have been designed or trained with insufficient subjective data with a small degree of stimulus variability. This has led to challenges for those algorithms to handle complexity and diversity of real-world digital content. Perceptual evidence from human subjects serves as a grounding for the development of advanced IQA algorithms. It is thus critical to acquire reliable subjective data with controlled perception experiments that faithfully reflect human behavioural responses to distortions in visual signals. In this paper, we present a new study of image quality perception where subjective ratings were collected in a controlled lab environment. We investigate how quality perception is affected by a combination of different categories of images and different types and levels of distortions. The database will be made publicly available to facilitate calibration and validation of IQA algorithms.