Search

In [1] the asymptotic charges of p-form gauge theories in any dimension are studied. Here we prove an existence and uniqueness theorem for the duality map linking asymptotic electric-like charges of the dual descriptions and we give it an algebraic topology interpretation. As a result the duality map has a topological nature and ensures the charge of a description has information of the dual description. The result of the theorem could be generalized to more generic gauge theories where the gauge field is a mixed symmetry tensor leading to a deeper understanding of gauge theories, of the non-trivial charges associated to them and of the duality of their observable.

Deep Learning-based Reduced Order Models (DL-ROMs) provide nowadays a well-established class of accurate surrogate models for complex physical systems described by parametrized PDEs, by nonlinearly compressing the solution manifold into a handful of latent coordinates. Until now, design and application of DL-ROMs mainly focused on physically parameterized problems. Within this work, we provide a novel extension of these architectures to problems featuring geometrical variability and parametrized domains, namely, we propose Continuous Geometry-Aware DL-ROMs (CGA-DL-ROMs). In particular, the space-continuous nature of the proposed architecture matches the need to deal with multi-resolution datasets, which are quite common in the case of geometrically parametrized problems. Moreover, CGA-DL-ROMs are endowed with a strong inductive bias that makes them aware of geometrical parametrizations, thus enhancing both the compression capability and the overall performance of the architecture. Within this work, we justify our findings through a thorough theoretical analysis, and we practically validate our claims by means of a series of numerical tests encompassing physically-and-geometrically parametrized PDEs, ranging from the unsteady Navier-Stokes equations for fluid dynamics to advection-diffusion-reaction equations for mathematical biology., Comment: 30 pages, 15 figures

We compute the surface charges associated to $p-$form gauge fields in arbitrary spacetime dimension for large values of the radial coordinate. In the critical dimension where radiation and Coulomb falloff coincide we find asymptotic charges involving asymptotic parameters, i.e. parameters with a component of order zero in the radial coordinate. However, in different dimensions we still find nontrivial asymptotic charges now involving parameters that are not asymptotic times the radiation-order fields. For $p$=1 and $D>4$, our charges thus differ from those presented in the literature. We then show that under Hodge duality electric charges for $p-$forms are mapped to magnetic charges for the dual $q-$forms, with $q = D-p-2$. For charges involving fields with radiation falloffs the duality relates charges that are finite and nonvanishing. For the case of Coulomb falloffs, above or below the critical dimension, Hodge duality exchanges overleading charges in one theory with subleading ones in its dual counterpart., Comment: 30 pages

We propose the use of the Extended Kalman Filter (EKF) for online data assimilation and update of a dynamic model, preliminary identified through the Sparse Identification of Nonlinear Dynamics (SINDy). This data-driven technique may avoid biases due to incorrect modelling assumptions and exploits SINDy to approximate the system dynamics leveraging a predefined library of functions, where active terms are selected and weighted by a sparse set of coefficients. This results in a physically-sound and interpretable dynamic model allowing to reduce epistemic uncertainty often affecting machine learning approaches. Treating the SINDy model coefficients as random variables, we propose to update them while acquiring (possibly noisy) system measurements, thus enabling the online identification of time-varying systems. These changes can stem from, e.g., varying operational conditions or unforeseen events. The EKF performs model adaptation through joint state-parameters estimation, with the Jacobian matrices required to computed the model sensitivity inexpensively evaluated from the SINDy model formulation. The effectiveness of this approach is demonstrated through three case studies: (i) a Lokta-Volterra model in which all parameters simultaneously evolve during the observation period; (ii) a Selkov model where the system undergoes a bifurcation not seen during the SINDy training; (iii) a MEMS arch exhibiting a 1:2 internal resonance. The ability of EKF of recovering inactivated functional terms from the SINDy library, or discarding unnecessary contribution, is also highlighted. Based on the presented applications, this method shows strong promise for handling time-varying nonlinear dynamic systems possibly experiencing bifurcating behaviours., Comment: 20 pages, 8 figures

Effectively controlling systems governed by Partial Differential Equations (PDEs) is crucial in several fields of Applied Sciences and Engineering. These systems usually yield significant challenges to conventional control schemes due to their nonlinear dynamics, partial observability, high-dimensionality once discretized, distributed nature, and the requirement for low-latency feedback control. Reinforcement Learning (RL), particularly Deep RL (DRL), has recently emerged as a promising control paradigm for such systems, demonstrating exceptional capabilities in managing high-dimensional, nonlinear dynamics. However, DRL faces challenges including sample inefficiency, robustness issues, and an overall lack of interpretability. To address these issues, we propose a data-efficient, interpretable, and scalable Dyna-style Model-Based RL framework for PDE control, combining the Sparse Identification of Nonlinear Dynamics with Control (SINDy-C) algorithm and an autoencoder (AE) framework for the sake of dimensionality reduction of PDE states and actions. This novel approach enables fast rollouts, reducing the need for extensive environment interactions, and provides an interpretable latent space representation of the PDE forward dynamics. We validate our method on two PDE problems describing fluid flows - namely, the 1D Burgers equation and 2D Navier-Stokes equations - comparing it against a model-free baseline, and carrying out an extensive analysis of the learned dynamics.

In this paper tackle the problem of computing the ranks of certain eulerian magnitude homology groups of a graph G. First, we analyze the computational cost of our problem and prove that it is #W[1]-complete. Then we develop the first diagonal algorithm, a breadth-first-search-based algorithm parameterized by the diameter of the graph to calculate the ranks of the homology groups of interest. To do this, we leverage the close relationship between the combinatorics of the homology boundary map and the substructures appearing in the graph. We then discuss the feasibility of the presented algorithm and consider future perspectives.

The superradiance instability of rotating black holes leads to the formation of an ultralight boson cloud with distinctive observational signatures, making black holes an effective probe of ultralight dark matter. However, around black holes in a binary system, the superradiance effect of such clouds can be terminated by tidal perturbations from the companion, leading to cloud depletion. In this study, we perform the first analysis of the impact of this termination effect on superradiant black hole binaries which are realistically modelled after their statistics in our Galaxy. Working with a dataset of approximately $10^7$ black hole binaries simulated using the Stellar EVolution for N-body (SEVN) population synthesis code, we identify the superradiant candidates and those that manage to survive the termination effect. We then calculate the cloud survival rate for various boson masses and black hole spin models. Our findings reveal that the $l=m=1$ cloud modes are generally stable against termination, whereas the $l=m=2$ modes can be significantly affected, with survival rates dropping below $10\%$ for boson masses below approximately $0.5\times 10^{-12}$ eV. In addition, our analysis indicates that clouds that overcome termination typically exhibit a higher superradiant growth rate and therefore a higher detectability., Comment: 13 pages, 10 figures

Selecting distinct quantum pathways in two-dimensional electronic spectroscopy (2DES) can give particularly deep insights into coherent and incoherent interactions and quantum dynamics in various materials. This includes isolating rephasing and non-rephasing pathways for conventional single-quantum 2DES, but also the ability to record double- and zero-quantum spectra. Such advanced 2DES schemes usually require phase-cycling when performed in a partially or fully collinear geometry. A particularly simple and effective implementation of 2DES utilizes an in-line birefringent interferometer, the Translating-Wedge-based Identical pulses eNcoding System (TWINS), for the generation of an inherently phase-stable collinear excitation pulse pair. Here, we demonstrate how the TWINS can be adapted to allow for phase-cycling and experimental access to isolated quantum pathways. These new capabilities are demonstrated by recording rephasing, non-rephasing, zero-quantum and double-quantum 2DES on a molecular J-aggregate. This easy-to-implement extension opens up new experimental possibilities for TWINS-based 2DES in multidimensional all-optical and photoemission spectroscopy and microscopy., Comment: 13 pages, 4 figures

ANNZ is a fast and simple algorithm which utilises artificial neural networks (ANNs), it was known as one of the pioneers of machine learning approaches to photometric redshift estimation decades ago. We enhanced the algorithm by introducing new activation functions like tanh, softplus, SiLU, Mish and ReLU variants; its new performance is then vigorously tested on legacy samples like the Luminous Red Galaxy (LRG) and Stripe-82 samples from SDSS, as well as modern galaxy samples like the Physics of the Accelerating Universe Survey (PAUS). This work focuses on testing the robustness of activation functions with respect to the choice of ANN architectures, particularly on its depth and width, in the context of galaxy photometric redshift estimation. Our upgraded algorithm, which we named ANNZ+, shows that the tanh and Leaky ReLU activation functions provide more consistent and stable results across deeper and wider architectures with > 1 per cent improvement in root-mean-square error ($\sigma_{\textrm{RMS}}$) and 68th percentile error ($\sigma_{68}$) when tested on SDSS data sets. While assessing its capabilities in handling high dimensional inputs, we achieved an improvement of 11 per cent in $\sigma_{\textrm{RMS}}$ and 6 per cent in $\sigma_{68}$ with the tanh activation function when tested on the 40-narrowband PAUS dataset; it even outperformed ANNZ2, its supposed successor, by 44 per cent in $\sigma_{\textrm{RMS}}$. This justifies the effort to upgrade the 20-year-old ANNZ, allowing it to remain viable and competitive within the photo-z community today. The updated algorithm ANNZ+ is publicly available at https://github.com/imdadmpt/ANNzPlus., Comment: 37 pages, 9 figures, submitted to JCAP

Steering a system towards a desired target in a very short amount of time is challenging from a computational standpoint. Indeed, the intrinsically iterative nature of optimal control problems requires multiple simulations of the physical system to be controlled. Moreover, the control action needs to be updated whenever the underlying scenario undergoes variations. Full-order models based on, e.g., the Finite Element Method, do not meet these requirements due to the computational burden they usually entail. On the other hand, conventional reduced order modeling techniques such as the Reduced Basis method, are intrusive, rely on a linear superimposition of modes, and lack of efficiency when addressing nonlinear time-dependent dynamics. In this work, we propose a non-intrusive Deep Learning-based Reduced Order Modeling (DL-ROM) technique for the rapid control of systems described in terms of parametrized PDEs in multiple scenarios. In particular, optimal full-order snapshots are generated and properly reduced by either Proper Orthogonal Decomposition or deep autoencoders (or a combination thereof) while feedforward neural networks are exploited to learn the map from scenario parameters to reduced optimal solutions. Nonlinear dimensionality reduction therefore allows us to consider state variables and control actions that are both low-dimensional and distributed. After (i) data generation, (ii) dimensionality reduction, and (iii) neural networks training in the offline phase, optimal control strategies can be rapidly retrieved in an online phase for any scenario of interest. The computational speedup and the high accuracy obtained with the proposed approach are assessed on different PDE-constrained optimization problems, ranging from the minimization of energy dissipation in incompressible flows modelled through Navier-Stokes equations to the thermal active cooling in heat transfer.

In this work, we present the novel mathematical framework of latent dynamics models (LDMs) for reduced order modeling of parameterized nonlinear time-dependent PDEs. Our framework casts this latter task as a nonlinear dimensionality reduction problem, while constraining the latent state to evolve accordingly to an (unknown) dynamical system. A time-continuous setting is employed to derive error and stability estimates for the LDM approximation of the full order model (FOM) solution. We analyze the impact of using an explicit Runge-Kutta scheme in the time-discrete setting, resulting in the $\Delta\text{LDM}$ formulation, and further explore the learnable setting, $\Delta\text{LDM}_\theta$, where deep neural networks approximate the discrete LDM components, while providing a bounded approximation error with respect to the FOM. Moreover, we extend the concept of parameterized Neural ODE - recently proposed as a possible way to build data-driven dynamical systems with varying input parameters - to be a convolutional architecture, where the input parameters information is injected by means of an affine modulation mechanism, while designing a convolutional autoencoder neural network able to retain spatial-coherence, thus enhancing interpretability at the latent level. Numerical experiments, including the Burgers' and the advection-reaction-diffusion equations, demonstrate the framework's ability to obtain, in a multi-query context, a time-continuous approximation of the FOM solution, thus being able to query the LDM approximation at any given time instance while retaining a prescribed level of accuracy. Our findings highlight the remarkable potential of the proposed LDMs, representing a mathematically rigorous framework to enhance the accuracy and approximation capabilities of reduced order modeling for time-dependent parameterized PDEs., Comment: 43 pages

In this paper we describe a numerical framework for achieving passive thermal cloaking of arbitrary shapes in both static and transient regimes. The design strategy is cast as the solution of an optimal control problem (OCP) for the heat equation where the coefficients of the thermal diffusivity matrix take the role of control functions and the distance between the uncloaked and the cloaked field is minimized in a suitable observation domain. The control actions enter bilinearly in the heat equation, thus making the resulting OCP nonlinear, and its analysis nontrivial. We show that optimal diffusivity coefficients exist both for the static and the transient case; we derive a system of first-order necessary optimality conditions; finally, we carry out their numerical approximation using the Finite Element Method. A series of numerical test cases assess the capability of our strategy to tackle passive thermal cloaking of arbitrarily complex two-dimensional objects., Comment: Accepted for publication in a future issue of the Proceedings of the Royal Society A

The simulation of many complex phenomena in engineering and science requires solving expensive, high-dimensional systems of partial differential equations (PDEs). To circumvent this, reduced-order models (ROMs) have been developed to speed up computations. However, when governing equations are unknown or partially known, typically ROMs lack interpretability and reliability of the predicted solutions. In this work we present a data-driven, non-intrusive framework for building ROMs where the latent variables and dynamics are identified in an interpretable manner and uncertainty is quantified. Starting from a limited amount of high-dimensional, noisy data the proposed framework constructs an efficient ROM by leveraging variational autoencoders for dimensionality reduction along with a newly introduced, variational version of sparse identification of nonlinear dynamics (SINDy), which we refer to as Variational Identification of Nonlinear Dynamics (VINDy). In detail, the method consists of Variational Encoding of Noisy Inputs (VENI) to identify the distribution of reduced coordinates. Simultaneously, we learn the distribution of the coefficients of a pre-determined set of candidate functions by VINDy. Once trained offline, the identified model can be queried for new parameter instances and new initial conditions to compute the corresponding full-time solutions. The probabilistic setup enables uncertainty quantification as the online testing consists of Variational Inference naturally providing Certainty Intervals (VICI). In this work we showcase the effectiveness of the newly proposed VINDy method in identifying interpretable and accurate dynamical system for the R\"ossler system with different noise intensities and sources. Then the performance of the overall method - named VENI, VINDy, VICI - is tested on PDE benchmarks including structural mechanics and fluid dynamics.

Narrowband galaxy surveys have recently gained interest as a promising method to achieve the necessary accuracy on the photometric redshift estimate of individual galaxies for stage-IV cosmological surveys. One key advantage is the ability to provide higher spectral resolution information about galaxies that should allow a more accurate and precise estimation of galaxy stellar population properties. However, the impact of adding narrow-band photometry on the stellar population properties estimate is largely unexplored. The scope of this work is two-fold: on one side, leveraging the predictive power of broad-band and narrow-band data to infer galaxy physical properties such as stellar masses, ages, star formation rates and metallicities. On the other hand, evaluating the improvement of performance in estimating galaxy properties when we use narrow-band data instead of broad-band. In this work we measure the stellar population properties of a sample of galaxies in the COSMOS field for which both narrowband and broadband data are available. In particular, we employ narrowband data from PAUS and broad-band data from CFHTLS. We use two different spectral energy distribution fitting codes to measure galaxy properties, namely CIGALE and Prospector. We find that the increased spectral resolution of narrow-band photometry does not yield a substantial improvement on constraining galaxy properties using spectral energy distribution fitting. Still we find that we obtain a more diverse distribution of metallicities and dust optical depths with cigale when employing the narrowband data. The effect is not as prominent as expected, which we relate this to the low narrowband SNR of a majority of the galaxies, the respective drawbacks of both codes as well as the coverage only in the optical regime. The measured properties are afterwards compared to the COSMOS2020 catalogue, showing good agreement., Comment: Published in A&A. 18 pages, 16 figures, 2 tables

Digital twins require computationally-efficient reduced-order models (ROMs) that can accurately describe complex dynamics of physical assets. However, constructing ROMs from noisy high-dimensional data is challenging. In this work, we propose a data-driven, non-intrusive method that utilizes stochastic variational deep kernel learning (SVDKL) to discover low-dimensional latent spaces from data and a recurrent version of SVDKL for representing and predicting the evolution of latent dynamics. The proposed method is demonstrated with two challenging examples -- a double pendulum and a reaction-diffusion system. Results show that our framework is capable of (i) denoising and reconstructing measurements, (ii) learning compact representations of system states, (iii) predicting system evolution in low-dimensional latent spaces, and (iv) quantifying modeling uncertainties.

This paper proposes a novel waveform design approach named COSMIC (Connectivity-Oriented Sensing Method for Imaging and Communication). The generation of radio images of the environment is enabled by imposing an extended orthogonality condition among waveforms. In contrast to traditional time, frequency, or space multiplexing systems, COSMIC achieves orthogonality through algebraic precoding of the transmitted waveforms from all antennas. Furthermore, by leveraging the degrees of freedom inherent in the assumption that the imaging field of view is significantly smaller than the length of the transmitted signals, the resulting waveforms are designed to convey communication symbols without impeding the sensing operations. We analytically demonstrate that the efficiency of imaging and communication sub-systems can be tuned according to the specific use case requirements. Simulation results indicate that COSMIC waveforms facilitate accurate environmental imaging while preserving near-optimal communication performance in terms of error rate and capacity.

In the context of structural health monitoring (SHM), the selection and extraction of damage-sensitive features from raw sensor recordings represent a critical step towards solving the inverse problem underlying the structural health identification. This work introduces a new way to enhance stochastic approaches to SHM through the use of deep neural networks. A learnable feature extractor and a feature-oriented surrogate model are synergistically exploited to evaluate a likelihood function within a Markov chain Monte Carlo sampling algorithm. The feature extractor undergoes a supervised pairwise training to map sensor recordings onto a low-dimensional metric space, which encapsulates the sensitivity to structural health parameters. The surrogate model maps the structural health parameters onto their feature description. The procedure enables the updating of beliefs about structural health parameters, effectively replacing the need for a computationally expensive numerical (finite element) model. A preliminary offline phase involves the generation of a labeled dataset to train both the feature extractor and the surrogate model. Within a simulation-based SHM framework, training vibration responses are cost-effectively generated by means of a multi-fidelity surrogate modeling strategy to approximate sensor recordings under varying damage and operational conditions. The multi-fidelity surrogate exploits model order reduction and artificial neural networks to speed up the data generation phase while ensuring the damage-sensitivity of the approximated signals. The proposed strategy is assessed through three synthetic case studies, demonstrating remarkable results in terms of accuracy of the estimated quantities and computational efficiency.

The coupling of Proper Orthogonal Decomposition (POD) and deep learning-based ROMs (DL-ROMs) has proved to be a successful strategy to construct non-intrusive, highly accurate, surrogates for the real time solution of parametric nonlinear time-dependent PDEs. Inexpensive to evaluate, POD-DL-ROMs are also relatively fast to train, thanks to their limited complexity. However, POD-DL-ROMs account for the physical laws governing the problem at hand only through the training data, that are usually obtained through a full order model (FOM) relying on a high-fidelity discretization of the underlying equations. Moreover, the accuracy of POD-DL-ROMs strongly depends on the amount of available data. In this paper, we consider a major extension of POD-DL-ROMs by enforcing the fulfillment of the governing physical laws in the training process -- that is, by making them physics-informed -- to compensate for possible scarce and/or unavailable data and improve the overall reliability. To do that, we first complement POD-DL-ROMs with a trunk net architecture, endowing them with the ability to compute the problem's solution at every point in the spatial domain, and ultimately enabling a seamless computation of the physics-based loss by means of the strong continuous formulation. Then, we introduce an efficient training strategy that limits the notorious computational burden entailed by a physics-informed training phase. In particular, we take advantage of the few available data to develop a low-cost pre-training procedure; then, we fine-tune the architecture in order to further improve the prediction reliability. Accuracy and efficiency of the resulting pre-trained physics-informed DL-ROMs (PTPI-DL-ROMs) are then assessed on a set of test cases ranging from non-affinely parametrized advection-diffusion-reaction equations, to nonlinear problems like the Navier-Stokes equations for fluid flows., Comment: 38 pages

We develop a novel deep learning technique, termed Deep Orthogonal Decomposition (DOD), for dimensionality reduction and reduced order modeling of parameter dependent partial differential equations. The approach consists in the construction of a deep neural network model that approximates the solution manifold through a continuously adaptive local basis. In contrast to global methods, such as Principal Orthogonal Decomposition (POD), the adaptivity allows the DOD to overcome the Kolmogorov barrier, making the approach applicable to a wide spectrum of parametric problems. Furthermore, due to its hybrid linear-nonlinear nature, the DOD can accommodate both intrusive and nonintrusive techniques, providing highly interpretable latent representations and tighter control on error propagation. For this reason, the proposed approach stands out as a valuable alternative to other nonlinear techniques, such as deep autoencoders. The methodology is discussed both theoretically and practically, evaluating its performances on problems featuring nonlinear PDEs, singularities, and parametrized geometries.

Measured data from a dynamical system can be assimilated into a predictive model by means of Kalman filters. Nonlinear extensions of the Kalman filter, such as the Extended Kalman Filter (EKF), are required to enable the joint estimation of (possibly nonlinear) system dynamics and of input parameters. To construct the evolution model used in the prediction phase of the EKF, we propose to rely on the Sparse Identification of Nonlinear Dynamics (SINDy). SINDy enables to identify the evolution model directly from preliminary acquired data, thus avoiding possible bias due to wrong assumptions and incorrect modelling of the system dynamics. Moreover, the numerical integration of a SINDy model leads to great computational savings compared to alternate strategies based on, e.g., finite elements. Last, SINDy allows an immediate definition of the Jacobian matrices required by the EKF to identify system dynamics and properties, a derivation that is usually extremely involved with physical models. As a result, combining the EKF with SINDy provides a data-driven computationally efficient, easy-to-apply approach for the identification of nonlinear systems, capable of robust operation even outside the range of training of SINDy. To demonstrate the potential of the approach, we address the identification of a linear non-autonomous system consisting of a shear building model excited by real seismograms, and the identification of a partially observed nonlinear system. The challenge arising from the use of SINDy when the system state is not entirely accessible has been relieved by means of time-delay embedding. The great accuracy and the small uncertainty associated with the state identification, where the state has been augmented to include system properties, underscores the great potential of the proposed strategy, paving the way for the setting of predictive digital twins in different fields., Comment: 26 pages, 11 figures

Gravitational solitons (gravisolitons) are particular exact solutions of Einstein field equation in vacuum build on a given background solution. Their interpretation is not yet fully clear but they contain many of the physically relevant solutions low $N$-solitons solutions. However, a systematic study and characterization of gravisolitons solution for every $N$ is lacking and their relevance in a theory of quantum gravity is not fully understood. This work aims to investigate and characterize some properties of $N$-axialsoliton solutions such as their asymptotically behaviour and asymptotic symmetries given minimal assumptions on the background metric. We develop an explicit systematic asymptotically expansion for the $N$-axialsoliton solution and we compute the leading order of the asymptotic killing vectors. Moreover, in the perspective to better understand the role of gravisolitons in quantum gravity we make a link, and a one of the first explicit test, to the corner symmetry proposal deriving which subalgebra of the universal corner symmetry algebra is generated by the asymptotic Killing vectors of $N$-axialsoliton solution. In the spirit of the corner proposal, the axialgravisoliton corner symmetry algebra (\textfrak{agcsa}) can be useful for the quantization of the non-asymptotically flat sector of gravity while, in the spirit of IR triangle, new soft theorems and memory effects could emerge.

Tailored light-matter interactions in the strong coupling regime enable the manipulation and control of quantum systems with up to unit efficiency, with applications ranging from quantum information to photochemistry. While strong light-matter interactions are readily induced at the valence electron level using long-wavelength radiation, comparable phenomena have been only recently observed with short wavelengths, accessing highly-excited multi-electron and inner-shell electron states. However, the quantum control of strong-field processes at short wavelengths has not been possible, so far, due to the lack of pulse shaping technologies in the extreme ultraviolet (XUV) and X-ray domain. Here, exploiting pulse shaping of the seeded free-electron laser (FEL) FERMI, we demonstrate the strong-field quantum control of ultrafast Rabi dynamics in helium atoms with high fidelity. Our approach unravels a strong dressing of the ionization continuum, otherwise elusive to experimental observables. The latter is exploited to achieve control of the total ionization rate, with prospective applications in many XUV and soft X-ray experiments. Leveraging recent advances in intense few-femtosecond to attosecond XUV to soft X-ray light sources, our results open an avenue to the efficient manipulation and selective control of core electron processes and electron correlation phenomena in real time.

On September 26, 2022, NASA's Double Asteroid Redirection Test (DART) mission successfully impacted Dimorphos, the natural satellite of the binary near-Earth asteroid (65803) Didymos. Numerical simulations of the impact provide a means to explore target surface material properties and structures, consistent with the observed momentum deflection efficiency, ejecta cone geometry, and ejected mass. Our simulation, which best matches observations, indicates that Dimorphos is weak, with a cohesive strength of less than a few pascals (Pa), similar to asteroids (162173) Ryugu and (101955) Bennu. We find that a bulk density of Dimorphos, rhoB, lower than 2400 kg/m3, and a low volume fraction of boulders (<40 vol%) on the surface and in the shallow subsurface, are consistent with measured data from the DART experiment. These findings suggest Dimorphos is a rubble pile that might have formed through rotational mass shedding and re-accumulation from Didymos. Our simulations indicate that the DART impact caused global deformation and resurfacing of Dimorphos. ESA's upcoming Hera mission may find a re-shaped asteroid, rather than a well-defined crater.

In this work we analyze the effectiveness of the Sparse Identification of Nonlinear Dynamics (SINDy) technique on three benchmark datasets for nonlinear identification, to provide a better understanding of its suitability when tackling real dynamical systems. While SINDy can be an appealing strategy for pursuing physics-based learning, our analysis highlights difficulties in dealing with unobserved states and non-smooth dynamics. Due to the ubiquity of these features in real systems in general, and control applications in particular, we complement our analysis with hands-on approaches to tackle these issues in order to exploit SINDy also in these challenging contexts., Comment: Submitted to IFAC SYSID 2024

A growing number of contributions in the literature have identified a puzzle in the European carbon allowance (EUA) market. Specifically, a persistent cost-of-carry spread (C-spread) over the risk-free rate has been observed. We are the first to explain the anomalous C-spread with the credit spread of the corporates involved in the emission trading scheme. We obtain statistical evidence that the C-spread is cointegrated with both this credit spread and the risk-free interest rate. This finding has a relevant policy implication: the most effective solution to solve the market anomaly is including the EUA in the list of European Central Bank eligible collateral for refinancing operations. This change in the ECB monetary policy operations would greatly benefit the carbon market and the EU green transition.

L\'evy-driven Ornstein-Uhlenbeck (OU) processes represent an intriguing class of stochastic processes that have garnered interest in the energy sector for their ability to capture typical features of market dynamics. However, in the current state of play, Monte Carlo simulations of these processes are not straightforward for two main reasons: i) algorithms are available only for some specific processes within this class; ii) they are often computationally expensive. In this paper, we introduce a new simulation technique designed to address both challenges. It relies on the numerical inversion of the characteristic function, offering a general methodology applicable to all L\'evy-driven OU processes. Moreover, leveraging FFT, the proposed methodology ensures fast and accurate simulations, providing a solid basis for the widespread adoption of these processes in the energy sector. Lastly, the algorithm allows an optimal control of the numerical error. We apply the technique to the pricing of energy derivatives, comparing the results with the existing benchmarks. Our findings indicate that the proposed methodology is at least one order of magnitude faster than the existing algorithms, while maintaining an equivalent level of accuracy.

This paper illustrates the potential of an Integrated Sensing and Communication (ISAC) system, operating in the sub-6 GHz frequency range, for Synthetic Aperture Radar (SAR) imaging via an Unmanned Aerial Vehicle (UAV) employed as an aerial base station. The primary aim is to validate the system's ability to generate SAR imagery within the confines of modern communication standards, including considerations like power limits, carrier frequency, bandwidth, and other relevant parameters. The paper presents two methods for processing the signal reflected by the scene. Additionally, we analyze two key performance indicators for their respective fields, the Noise Equivalent Sigma Zero (NESZ) and the Bit Error Rate (BER), using the QUAsi Deterministic RadIo channel GenerAtor (QuaDRiGa), demonstrating the system's capability to image buried targets in challenging scenarios. The paper shows simulated Impulse Response Functions (IRF) as possible pulse compression techniques under different assumptions. An experimental campaign is conducted to validate the proposed setup by producing a SAR image of the environment captured using a UAV flying with a Software-Defined Radio (SDR) as a payload.

A Particle Swarm Optimizer for the search of balanced Boolean functions with good cryptographic properties is proposed in this paper. The algorithm is a modified version of the permutation PSO by Hu, Eberhart and Shi which preserves the Hamming weight of the particles positions, coupled with the Hill Climbing method devised by Millan, Clark and Dawson to improve the nonlinearity and deviation from correlation immunity of Boolean functions. The parameters for the PSO velocity equation are tuned by means of two meta-optimization techniques, namely Local Unimodal Sampling (LUS) and Continuous Genetic Algorithms (CGA), finding that CGA produces better results. Using the CGA-evolved parameters, the PSO algorithm is then run on the spaces of Boolean functions from $n=7$ to $n=12$ variables. The results of the experiments are reported, observing that this new PSO algorithm generates Boolean functions featuring similar or better combinations of nonlinearity, correlation immunity and propagation criterion with respect to the ones obtained by other optimization methods., Comment: Extended version of the poster paper "Heuristic Search by Particle Swarm Optimization of Boolean Functions for Cryptographic Applications" published in GECCO 2015

We present photometric redshifts (photo-$z$) for the deep wide fields of the Physics of the Accelerating Universe Survey (PAUS), covering an area of $\sim$50 deg$^{2}$, for $\sim$1.8 million objects up to $i_{\textrm{AB}}<23$. The PAUS deep wide fields overlap with the W1 and W3 fields from CFHTLenS and the G09 field from KiDS/GAMA. Photo-$z$ are estimated using the 40 narrow bands (NB) of PAUS and the broad bands (BB) of CFHTLenS and KiDS. We compute the redshifts with the SED template-fitting code BCNZ, with a modification in the calibration technique of the zero-point between the observed and the modelled fluxes, that removes any dependence on spectroscopic redshift samples. We enhance the redshift accuracy by introducing an additional photo-$z$ estimate ($z_{\textrm{b}}$), obtained through the combination of the BCNZ and the BB-only photo-$z$. Comparing with spectroscopic redshifts estimates ($z_{\textrm{s}}$), we obtain a $\sigma_{68} \simeq 0.019$ for all galaxies with $i_{\textrm{AB}}<23$ and a typical bias $|z_{\textrm{b}}-z_{\textrm{s}}|$ smaller than 0.01. For $z_{\textrm{b}} \sim (0.10-0.75)$ we find $\sigma_{68} \simeq (0.003-0.02)$, this is a factor of $10-2$ higher accuracy than the corresponding BB-only results. We obtain similar performance when we split the samples into red (passive) and blue (active) galaxies. We validate the redshift probability $p(z)$ obtained by BCNZ and compare its performance with that of $z_{\textrm{b}}$. These photo-$z$ catalogues will facilitate important science cases, such as the study of galaxy clustering and intrinsic alignment at high redshifts ($z \lesssim 1$) and faint magnitudes., Comment: 24 pages, 26 figures, submitted to MNRAS

With the recent success of artificial intelligence in neuroscience, a number of deep learning (DL) models were proposed for classification, anomaly detection, and pattern recognition tasks in electroencephalography (EEG). EEG is a multi-channel time-series that provides information about the individual brain activity for diagnostics, neuro-rehabilitation, and other applications (including emotions recognition). Two main issues challenge the existing DL-based modeling methods for EEG: the high variability between subjects and the low signal-to-noise ratio making it difficult to ensure a good quality in the EEG data. In this paper, we propose two variational autoencoder models, namely vEEGNet-ver3 and hvEEGNet, to target the problem of high-fidelity EEG reconstruction. We properly designed their architectures using the blocks of the well-known EEGNet as the encoder, and proposed a loss function based on dynamic time warping. We tested the models on the public Dataset 2a - BCI Competition IV, where EEG was collected from 9 subjects and 22 channels. hvEEGNet was found to reconstruct the EEG data with very high-fidelity, outperforming most previous solutions (including our vEEGNet-ver3 ). Furthermore, this was consistent across all subjects. Interestingly, hvEEGNet made it possible to discover that this popular dataset includes a number of corrupted EEG recordings that might have influenced previous literature results. We also investigated the training behaviour of our models and related it with the quality and the size of the input EEG dataset, aiming at opening a new research debate on this relationship. In the future, hvEEGNet could be used as anomaly (e.g., artefact) detector in large EEG datasets to support the domain experts, but also the latent representations it provides could be used in other classification problems and EEG data generation.

We use the GALFORM semi-analytical galaxy formation model implemented in the Planck Millennium N-body simulation to build a mock galaxy catalogue on an observer's past lightcone. The mass resolution of this N-body simulation is almost an order of magnitude better than in previous simulations used for this purpose, allowing us to probe fainter galaxies and hence build a more complete mock catalogue at low redshifts. The high time cadence of the simulation outputs allows us to make improved calculations of galaxy properties and positions in the mock. We test the predictions of the mock against the Physics of the Accelerating Universe Survey, a narrow band imaging survey with highly accurate and precise photometric redshifts, which probes the galaxy population over a lookback time of 8 billion years. We compare the model against the observed number counts, redshift distribution and evolution of the observed colours and find good agreement; these statistics avoid the need for model-dependent processing of the observations. The model produces red and blue populations that have similar median colours to the observations. However, the bimodality of galaxy colours in the model is stronger than in the observations. This bimodality is reduced on including a simple model for errors in the GALFORM photometry. We examine how the model predictions for the observed galaxy colours change when perturbing key model parameters. This exercise shows that the median colours and relative abundance of red and blue galaxies provide constraints on the strength of the feedback driven by supernovae used in the model.

Electroencephalografic (EEG) data are complex multi-dimensional time-series that are very useful in many applications, from diagnostics to driving brain-computer interface systems. Their classification is still a challenging task, due to the inherent within- and between-subject variability and their low signal-to-noise ratio. On the other hand, the reconstruction of raw EEG data is even more difficult because of the high temporal resolution of these signals. Recent literature has proposed numerous machine and deep learning models that could classify, e.g., different types of movements, with an accuracy in the range 70% to 80% (with 4 classes). On the other hand, a limited number of works targeted the reconstruction problem, with very limited results. In this work, we propose vEEGNet, a DL architecture with two modules, i.e., an unsupervised module based on variational autoencoders to extract a latent representation of the data, and a supervised module based on a feed-forward neural network to classify different movements. To build the encoder and the decoder of VAE we exploited the well-known EEGNet network. We implemented two slightly different architectures of vEEGNet, thus showing state-of-the-art classification performance, and the ability to reconstruct both low-frequency and middle-range components of the raw EEG. Although preliminary, this work is promising as we found out that the low-frequency reconstructed signals are consistent with the so-called motor-related cortical potentials, well-known motor-related EEG patterns and we could improve over previous literature by reconstructing faster EEG components, too. Further investigations are needed to explore the potentialities of vEEGNet in reconstructing the full EEG data, generating new samples, and studying the relationship between classification and reconstruction performance.

Coherent multistatic radio imaging represents a pivotal opportunity for forthcoming wireless networks, which involves distributed nodes cooperating to achieve accurate sensing resolution and robustness. This paper delves into cooperative coherent imaging for vehicular radar networks. Herein, multiple radar-equipped vehicles cooperate to improve collective sensing capabilities and address the fundamental issue of distinguishing weak targets in close proximity to strong ones, a critical challenge for vulnerable road users protection. We prove the significant benefits of cooperative coherent imaging in the considered automotive scenario in terms of both probability of correct detection, evaluated considering several system parameters, as well as resolution capabilities, showcased by a dedicated experimental campaign wherein the collaboration between two vehicles enables the detection of the legs of a pedestrian close to a parked car. Moreover, as \textit{coherent} processing of several sensors' data requires very tight accuracy on clock synchronization and sensor's positioning -- referred to as \textit{phase synchronization} -- (such that to predict sensor-target distances up to a fraction of the carrier wavelength), we present a general three-step cooperative multistatic phase synchronization procedure, detailing the required information exchange among vehicles in the specific automotive radar context and assessing its feasibility and performance by hybrid Cram\'er-Rao bound., Comment: 13 pages

The near infrared background (NIRB) is the collective light from unresolved sources observed in the band 1-10 $\mu$m. The measured NIRB angular power spectrum on angular scales $\theta \gtrsim 1$ arcmin exceeds by roughly two order of magnitudes predictions from known galaxy populations. The nature of the sources producing these fluctuations is still unknown. Here we test primordial black holes (PBHs) as sources of the NIRB excess. Considering PBHs as a cold dark matter (DM) component, we model the emission of gas accreting onto PBHs in a cosmological framework. We account for both accretion in the intergalactic medium (IGM) and in DM haloes. We self consistently derive the IGM temperature evolution, considering ionization and heating due to X-ray emission from PBHs. Besides $\Lambda$CDM, we consider a model that accounts for the modification of the linear matter power spectrum due to the presence of PBHs; we also explore two PBH mass distributions, i.e. a $\delta$-function and a lognormal distribution. For each model, we compute the mean intensity and the angular power spectrum of the NIRB produced by PBHs with mass 1-$10^3~\mathrm{M}_{\odot}$. In the limiting case in which the entirety of DM is made of PBHs, the PBH emission contributes $<1$ per cent to the observed NIRB fluctuations. This value decreases to $<0.1$ per cent if current constraints on the abundance of PBHs are taken into account. We conclude that PBHs are ruled out as substantial contributors to the NIRB., Comment: Accepted for publication in MNRAS

Deep Learning is having a remarkable impact on the design of Reduced Order Models (ROMs) for Partial Differential Equations (PDEs), where it is exploited as a powerful tool for tackling complex problems for which classical methods might fail. In this respect, deep autoencoders play a fundamental role, as they provide an extremely flexible tool for reducing the dimensionality of a given problem by leveraging on the nonlinear capabilities of neural networks. Indeed, starting from this paradigm, several successful approaches have already been developed, which are here referred to as Deep Learning-based ROMs (DL-ROMs). Nevertheless, when it comes to stochastic problems parameterized by random fields, the current understanding of DL-ROMs is mostly based on empirical evidence: in fact, their theoretical analysis is currently limited to the case of PDEs depending on a finite number of (deterministic) parameters. The purpose of this work is to extend the existing literature by providing some theoretical insights about the use of DL-ROMs in the presence of stochasticity generated by random fields. In particular, we derive explicit error bounds that can guide domain practitioners when choosing the latent dimension of deep autoencoders. We evaluate the practical usefulness of our theory by means of numerical experiments, showing how our analysis can significantly impact the performance of DL-ROMs.

Many applications in computational physics involve approximating problems with microstructure, characterized by multiple spatial scales in their data. However, these numerical solutions are often computationally expensive due to the need to capture fine details at small scales. As a result, simulating such phenomena becomes unaffordable for many-query applications, such as parametrized systems with multiple scale-dependent features. Traditional projection-based reduced order models (ROMs) fail to resolve these issues, even for second-order elliptic PDEs commonly found in engineering applications. To address this, we propose an alternative nonintrusive strategy to build a ROM, that combines classical proper orthogonal decomposition (POD) with a suitable neural network (NN) model to account for the small scales. Specifically, we employ sparse mesh-informed neural networks (MINNs), which handle both spatial dependencies in the solutions and model parameters simultaneously. We evaluate the performance of this strategy on benchmark problems and then apply it to approximate a real-life problem involving the impact of microcirculation in transport phenomena through the tissue microenvironment.

High-fidelity numerical simulations of partial differential equations (PDEs) given a restricted computational budget can significantly limit the number of parameter configurations considered and/or time window evaluated for modeling a given system. Multi-fidelity surrogate modeling aims to leverage less accurate, lower-fidelity models that are computationally inexpensive in order to enhance predictive accuracy when high-fidelity data are limited or scarce. However, low-fidelity models, while often displaying important qualitative spatio-temporal features, fail to accurately capture the onset of instability and critical transients observed in the high-fidelity models, making them impractical as surrogate models. To address this shortcoming, we present a new data-driven strategy that combines dimensionality reduction with multi-fidelity neural network surrogates. The key idea is to generate a spatial basis by applying the classical proper orthogonal decomposition (POD) to high-fidelity solution snapshots, and approximate the dynamics of the reduced states - time-parameter-dependent expansion coefficients of the POD basis - using a multi-fidelity long-short term memory (LSTM) network. By mapping low-fidelity reduced states to their high-fidelity counterpart, the proposed reduced-order surrogate model enables the efficient recovery of full solution fields over time and parameter variations in a non-intrusive manner. The generality and robustness of this method is demonstrated by a collection of parametrized, time-dependent PDE problems where the low-fidelity model can be defined by coarser meshes and/or time stepping, as well as by misspecified physical features. Importantly, the onset of instabilities and transients are well captured by this surrogate modeling technique.

Computer-based simulations of non-invasive cardiac electrical outputs, such as electrocardiograms and body surface potential maps, usually entail severe computational costs due to the need of capturing fine-scale processes and to the complexity of the heart-torso morphology. In this work, we model cardiac electrical outputs by employing a coupled model consisting of a reaction-diffusion model - either the bidomain model or the most efficient pseudo-bidomain model - on the heart, and an elliptic model in the torso. We then solve the coupled problem with a segregated and staggered in-time numerical scheme, that allows for independent and infrequent solution in the torso region. To further reduce the computational load, main novelty of this work is in introduction of an interpolation method at the interface between the heart and torso domains, enabling the use of non-conforming meshes, and the numerical framework application to realistic cardiac and torso geometries. The reliability and efficiency of the proposed scheme is tested against the corresponding state-of-the-art bidomain-torso model. Furthermore, we explore the impact of torso spatial discretization and geometrical non-conformity on the model solution and the corresponding clinical outputs. The investigation of the interface interpolation method provides insights into the influence of torso spatial discretization and of the geometrical non-conformity on the simulation results and their clinical relevance., Comment: 26 pages,11 figures, 3 tables

Mesh-based simulations play a key role when modeling complex physical systems that, in many disciplines across science and engineering, require the solution of parametrized time-dependent nonlinear partial differential equations (PDEs). In this context, full order models (FOMs), such as those relying on the finite element method, can reach high levels of accuracy, however often yielding intensive simulations to run. For this reason, surrogate models are developed to replace computationally expensive solvers with more efficient ones, which can strike favorable trade-offs between accuracy and efficiency. This work explores the potential usage of graph neural networks (GNNs) for the simulation of time-dependent PDEs in the presence of geometrical variability. In particular, we propose a systematic strategy to build surrogate models based on a data-driven time-stepping scheme where a GNN architecture is used to efficiently evolve the system. With respect to the majority of surrogate models, the proposed approach stands out for its ability of tackling problems with parameter dependent spatial domains, while simultaneously generalizing to different geometries and mesh resolutions. We assess the effectiveness of the proposed approach through a series of numerical experiments, involving both two- and three-dimensional problems, showing that GNNs can provide a valid alternative to traditional surrogate models in terms of computational efficiency and generalization to new scenarios. We also assess, from a numerical standpoint, the importance of using GNNs, rather than classical dense deep neural networks, for the proposed framework.

The digital twin concept represents an appealing opportunity to advance condition-based and predictive maintenance paradigms for civil engineering systems, thus allowing reduced lifecycle costs, increased system safety, and increased system availability. This work proposes a predictive digital twin approach to the health monitoring, maintenance, and management planning of civil engineering structures. The asset-twin coupled dynamical system is encoded employing a probabilistic graphical model, which allows all relevant sources of uncertainty to be taken into account. In particular, the time-repeating observations-to-decisions flow is modeled using a dynamic Bayesian network. Real-time structural health diagnostics are provided by assimilating sensed data with deep learning models. The digital twin state is continually updated in a sequential Bayesian inference fashion. This is then exploited to inform the optimal planning of maintenance and management actions within a dynamic decision-making framework. A preliminary offline phase involves the population of training datasets through a reduced-order numerical model and the computation of a health-dependent control policy. The strategy is assessed on two synthetic case studies, involving a cantilever beam and a railway bridge, demonstrating the dynamic decision-making capabilities of health-aware digital twins.

Reaction systems are discrete dynamical systems that model biochemical processes in living cells using finite sets of reactants, inhibitors, and products. We investigate the computational complexity of a comprehensive set of problems related to the existence of fixed points and attractors in two constrained classes of reaction systems, in which either reactants or inhibitors are disallowed. These problems have biological relevance and have been extensively studied in the unconstrained case; however, they remain unexplored in the context of reactantless or inhibitorless systems. Interestingly, we demonstrate that although the absence of reactants or inhibitors simplifies the system's dynamics, it does not always lead to a reduction in the complexity of the considered problems., Comment: 32 pages. The previous version contained major mistakes. The wrong results are the former Corollaries 9, 10, 15, 17, 23, 25, 25, 32 and 36, Proposition 31 and Remark 30 (all following from a single mistake). We have now corrected all these results

We present a new simulation for Higgs boson production in association with bottom quarks ($b\bar{b}H$) at next-to-leading order (NLO) accuracy matched to parton showers in hadronic collisions. Both contributions, the standard one proportional to the bottom-quark Yukawa coupling and the loop-induced one proportional to the top-quark Yukawa coupling from the gluon-fusion process, are taken into account in a scheme with massive bottom quarks. Therefore, we provide the full simulation of the $b\bar{b}H$ final state in the Standard Model, which constitutes also a crucial background to measurements for Higgs-boson pair ($HH$) production at the Large Hadron Collider when at least one of the Higgs bosons decays to bottom quarks. So far, the modeling of the $b\bar{b}H$ final state induced one of the dominant theoretical uncertainties to $HH$ measurements, as the gluon-fusion component was described only at the leading order (LO) with uncertainties of $\mathcal{O}(100\%)$. Including NLO corrections in its simulation allows us to reduce the scale dependence to $\mathcal{O}(50\%)$ so that it becomes subdominant with respect to other systematic uncertainties. As a case study, we provide an in-depth analysis of the $b\bar{b}H$ background to $HH$ measurements with realistic selection cuts in the $2b2\gamma$ channel. We also compare our novel simulation with the currently-employed ones, discussing possible issues and shortcomings of a scheme with massless bottom quarks. Finally, we propagate the effect of the new $b\bar{b}H$ simulation to $HH$ searches in the $2b2\gamma$ and $2b2\tau$ final states, and we find an improvement of up to 10% (20%) on the current (HL-LHC) limits on the $HH$ cross section., Comment: 25 pages, 5 figures, 3 tables

Geometric Semantic Geometric Programming (GSGP) is one of the most prominent Genetic Programming (GP) variants, thanks to its solid theoretical background, the excellent performance achieved, and the execution time significantly smaller than standard syntax-based GP. In recent years, a new mutation operator, Geometric Semantic Mutation with Local Search (GSM-LS), has been proposed to include a local search step in the mutation process based on the idea that performing a linear regression during the mutation can allow for a faster convergence to good-quality solutions. While GSM-LS helps the convergence of the evolutionary search, it is prone to overfitting. Thus, it was suggested to use GSM-LS only for a limited number of generations and, subsequently, to switch back to standard geometric semantic mutation. A more recently defined variant of GSGP (called GSGP-reg) also includes a local search step but shares similar strengths and weaknesses with GSM-LS. Here we explore multiple possibilities to limit the overfitting of GSM-LS and GSGP-reg, ranging from adaptive methods to estimate the risk of overfitting at each mutation to a simple regularized regression. The results show that the method used to limit overfitting is not that important: providing that a technique to control overfitting is used, it is possible to consistently outperform standard GSGP on both training and unseen data. The obtained results allow practitioners to better understand the role of local search in GSGP and demonstrate that simple regularization strategies are effective in controlling overfitting., Comment: 22 pages, 4 figures, 1 table

Networked sensing refers to the capability of properly orchestrating multiple sensing terminals to enhance specific figures of merit, e.g., positioning accuracy or imaging resolution. Regarding radio-based sensing, it is essential to understand \textit{when} and \textit{how} sensing terminals should be orchestrated, namely the best cooperation that trades between performance and cost (e.g., energy consumption, communication overhead, and complexity). This paper addresses networked sensing from a physics-driven perspective, aiming to provide a general theoretical benchmark to evaluate its \textit{imaging} performance bounds and to guide the sensing orchestration accordingly. Diffraction tomography theory (DTT) is the method to quantify the imaging resolution of any radio sensing experiment from inspection of its spectral (or wavenumber) content. In networked sensing, the image formation is based on the back-projection integral, valid for any network topology and physical configuration of the terminals. The \textit{wavefield networked sensing} is a framework in which multiple sensing terminals are orchestrated during the acquisition process to maximize the imaging quality (resolution and grating lobes suppression) by pursuing the deceptively simple \textit{wavenumber tessellation principle}. We discuss all the cooperation possibilities between sensing terminals and possible killer applications. Remarkably, we show that the proposed method allows obtaining high-quality images of the environment in limited bandwidth conditions, leveraging the coherent combination of multiple multi-static low-resolution images., Comment: 13 pages

POD-DL-ROMs have been recently proposed as an extremely versatile strategy to build accurate and reliable reduced order models (ROMs) for nonlinear parametrized partial differential equations, combining (i) a preliminary dimensionality reduction obtained through proper orthogonal decomposition (POD) for the sake of efficiency, (ii) an autoencoder architecture that further reduces the dimensionality of the POD space to a handful of latent coordinates, and (iii) a dense neural network to learn the map that describes the dynamics of the latent coordinates as a function of the input parameters and the time variable. Within this work, we aim at justifying the outstanding approximation capabilities of POD-DL-ROMs by means of a thorough error analysis, showing how the sampling required to generate training data, the dimension of the POD space, and the complexity of the underlying neural networks, impact on the solution accuracy. This decomposition, combined with the constructive nature of the proofs, allows us to formulate practical criteria to control the relative error in the approximation of the solution field of interest, and derive general error estimates. Furthermore, we show that, from a theoretical point of view, POD-DL-ROMs outperform several deep learning-based techniques in terms of model complexity. Finally, we validate our findings by means of suitable numerical experiments, ranging from parameter-dependent operators analytically defined to several parametrized PDEs.

Nonlinear spectroscopy with quantum entangled photons is an emerging field of research that holds the promise to achieve a superior signal-to-noise ratio and effectively isolate many-body interactions. Photon sources used for this purpose however lack the frequency tunability and spectral bandwidth demanded by contemporary molecular materials. Here, we present design strategies for efficient spontaneous parametric downconversion to generate biphoton states with adequate spectral bandwidth and at visible wavelengths. Importantly, we demonstrate, by suitable design of the nonlinear optical interaction, the scope to engineer the degree of spectral correlations between the photons of the pair. We also present an experimental methodology to effectively characterize such spectral correlations. Importantly, we believe that such a characterization tool can be effectively adapted as a spectroscopy platform to optically probe system-bath interactions in materials.

The High Granularity Timing Detector (HGTD) will be installed in the ATLAS experiment to mitigate pile-up effects during the High Luminosity (HL) phase of the Large Hadron Collider (LHC) at CERN. Low Gain Avalanche Detectors (LGADs) will provide high-precision measurements of the time of arrival of particles at the HGTD, improving the particle-vertex assignment. To cope with the high-radiation environment, LGADs have been optimized by adding carbon in the gain layer, thus reducing the acceptor removal rate after irradiation. Performances of several carbon-enriched LGAD sensors from different vendors, and irradiated with high fluences of 1.5 and 2.5 x 10^15 neq/cm2, have been measured in beam test campaigns during the years 2021 and 2022 at CERN SPS and DESY. This paper presents the results obtained with data recorded by an oscilloscope synchronized with a beam telescope which provides particle position information within a resolution of a few um. Collected charge, time resolution and hit efficiency measurements are presented. In addition, the efficiency uniformity is also studied as a function of the position of the incident particle inside the sensor pad.

Most of the approaches published in the literature to construct S-boxes via Cellular Automata (CA) work by either iterating a finite CA for several time steps, or by a one-shot application of the global rule. The main characteristic that brings together these works is that they employ a single CA rule to define the vectorial Boolean function of the S-box. In this work, we explore a different direction for the design of S-boxes that leverages on Orthogonal CA (OCA), i.e. pairs of CA rules giving rise to orthogonal Latin squares. The motivation stands on the facts that an OCA pair already defines a bijective transformation, and moreover the orthogonality property of the resulting Latin squares ensures a minimum amount of diffusion. We exhaustively enumerate all S-boxes generated by OCA pairs of diameter $4 \le d \le 6$, and measure their nonlinearity. Interestingly, we observe that for $d=4$ and $d=5$ all S-boxes are linear, despite the underlying CA local rules being nonlinear. The smallest nonlinear S-boxes emerges for $d=6$, but their nonlinearity is still too low to be used in practice. Nonetheless, we unearth an interesting structure of linear OCA S-boxes, proving that their Linear Components Space (LCS) is itself the image of a linear CA, or equivalently a polynomial code. We finally classify all linear OCA S-boxes in terms of their generator polynomials., Comment: 22 pages. Extended version of "On the Linear Components Space of S-boxes Generated by Orthogonal Cellular Automata" arXiv:2203.14365v, presented at ACRI 2022. Currently under submission at Natural Computing

Uncertainty quantification (UQ) tasks, such as sensitivity analysis and parameter estimation, entail a huge computational complexity when dealing with input-output maps involving the solution of nonlinear differential problems, because of the need to query expensive numerical solvers repeatedly. Projection-based reduced order models (ROMs), such as the Galerkin-reduced basis (RB) method, have been extensively developed in the last decades to overcome the computational complexity of high fidelity full order models (FOMs), providing remarkable speedups when addressing UQ tasks related with parameterized differential problems. Nonetheless, constructing a projection-based ROM that can be efficiently queried usually requires extensive modifications to the original code, a task which is non-trivial for nonlinear problems, or even not possible at all when proprietary software is used. Non-intrusive ROMs - which rely on the FOM as a black box - have been recently developed to overcome this issue. In this work, we consider ROMs exploiting proper orthogonal decomposition to construct a reduced basis from a set of FOM snapshots, and Gaussian process regression (GPR) to approximate the RB projection coefficients. Two different approaches, namely a global GPR and a tensor-decomposition-based GPR, are explored on a set of 3D time-dependent solid mechanics examples. Finally, the non-intrusive ROM is exploited to perform global sensitivity analysis (relying on both screening and variance-based methods) and parameter estimation (through Markov chain Monte Carlo methods), showing remarkable computational speedups and very good accuracy compared to high-fidelity FOMs.

The design of binary error-correcting codes is a challenging optimization problem with several applications in telecommunications and storage, which has also been addressed with metaheuristic techniques and evolutionary algorithms. Still, all these efforts focused on optimizing the minimum distance of unrestricted binary codes, i.e., with no constraints on their linearity, which is a desirable property for efficient implementations. In this paper, we present an Evolutionary Strategy (ES) algorithm that explores only the subset of linear codes of a fixed length and dimension. To that end, we represent the candidate solutions as binary matrices and devise variation operators that preserve their ranks. Our experiments show that up to length $n=14$, our ES always converges to an optimal solution with a full success rate, and the evolved codes are all inequivalent to the Best-Known Linear Code (BKLC) given by MAGMA. On the other hand, for larger lengths, both the success rate of the ES as well as the diversity of the evolved codes start to drop, with the extreme case of $(16,8,5)$ codes which all turn out to be equivalent to MAGMA's BKLC., Comment: 15 pages, 3 figures, 3 tables