Search

We analyze JWST MIRI/MRS observations of the infrared PAH bands in the nuclear and circumnuclear regions of local AGN from the GATOS Survey. In this work, we examine the PAH properties in the circumnuclear regions of AGN and AGN-outflows, and compare them to those in star-forming regions and the innermost regions of AGN. This study employs 4.9-28.1 micron sub-arcsecond angular resolution data to investigate the properties of PAH in three nearby sources (DL~30-40 Mpc). Our findings align with previous JWST studies, showing that the central regions of AGN show a larger fraction of neutral PAH molecules (i.e. elevated 11.3/6.2 and 11.3/7.7 PAH ratios) compared to star-forming galaxies. We find that the AGN might affect not only the PAH population in the innermost region but also in the extended regions up to ~kpc scales. By comparing our observations to PAH diagnostic diagrams, we find that, in general, regions located in the projected direction of the AGN-outflow occupy similar positions on the PAH diagnostic diagrams as those of the innermost regions of AGN. Star-forming regions that are not affected by the AGN in these galaxies share the same part of the diagram as Star-forming galaxies. We examine the potential of the PAH-H2 diagram to disentangle AGN versus star-forming activity. Our results suggest that in Sy-like AGN, illumination and feedback from the AGN might affect the PAH population at nuclear and kpc scales, in particular, the ionization state of the PAH grains. However, PAH sizes are rather similar. The carriers of the ionized PAH bands (6.2 and 7.7 micron) are less resilience than those of neutral PAH bands (11.3 micron), which might be particularly important for strongly AGN-host coupled systems. Therefore, caution must be applied when using PAH bands as star-formation rate indicators in these systems even at kpc scales, with the ionized ones being more affected by the AGN., Comment: Accepted for publication in A&A. 21 pages, 13 Figures

We present observations of the type-2 Seyfert NGC7172 obtained with the medium-resolution spectrometer (MRS) of the Mid-Infrared Instrument (MIRI) on board of the JWST. This galaxy hosts one of the lowest ionised gas mass outflow rates (Mout~0.005 M/yr) in a sample of six AGN with similar bolometric luminosities (log Lbol~44erg/s) within the Galactic Activity, Torus and Outflow Survey (GATOS). We aim to understand the properties of the ionised gas outflow, mainly using the emission lines from the neon transitions, that cover a broad range of ionisation potentials (IP) from ~20 eV to ~130 eV. We applied parametric and non-parametric methods to characterise the line emission and kinematics. The low excitation lines (IP<25eV, e.g.[NeII]) trace the rotating disc emission. The high excitation lines (IP>90eV, e.g.[NeV]), which are likely photoionised exclusively by the AGN, are expanding in the direction nearly perpendicular to the disc of the galaxy, with maximum projected velocities of ~350-500 km/s. In particular, [NeV] and [NeVI] lines reveal a biconical ionised gas outflow emerging N-S from the nuclear region, extending at least ~2.5"N and 3.8"S (projected distance of ~450 and 680 pc). Most of the emission arising in the northern part of the cone was not previously detected due to obscuration. Given the almost face-on orientation of the outflow and the almost edge-on orientation of the galaxy, NGC7172 may be a case of weak coupling. Nevertheless, we found evidence for positive feedback in two distinct outflowing clumps at projected distances of 3.1" and 4.3" (i.e. ~560 and 780 pc) SW from the AGN. We estimated a star formation rate in these regions using the [NeII] and [NeIII] luminosities of 0.08 M/yr, that is ~10% of that found in the circumnuclear ring. The star formation activity might have been triggered by the interaction between the ionised gas outflow and the ISM of the galaxy., Comment: 20 pages, 6 figures. Accepted for publication in A&A. Abstract adapted for the arxiv version

The ionisation cones of NGC5728 have a deficit of molecular gas based on millimetre observations of CO(2-1) emission. Although photoionisation from the active nucleus may lead to suppression of this transition, warm molecular gas can still be present. We report the detection of eight mid-infrared rotational H$_2$ lines throughout the central kiloparsec, including the ionisation cones, using integral field spectroscopic observations with JWST/MIRI MRS. The H$_2$ line ratios, characteristic of a power-law temperature distribution, indicate that the gas is warmest where it enters the ionisation cone through disk rotation, suggestive of shock excitation. In the nucleus, where the data can be combined with an additional seven ro-vibrational H$_2$ transitions, we find that moderate velocity (30 km s$^{-1}$) shocks in dense ($10^5$ cm$^{-3}$) gas, irradiated by an external UV field ($G_0 = 10^3$), do provide a good match to the full set. The warm molecular gas in the ionisation cone that is traced by the H$_2$ rotational lines has been heated to temperatures $>200$ K. Outside of the ionisation cone the molecular gas kinematics are undisturbed. However, within the ionisation cone, the kinematics are substantially perturbed, indicative of a radial flow, but one that is quantitatively different from the ionised lines. We argue that this outflow is in the plane of the disk, implying a short 50 pc acceleration zone up to speeds of about 400 km s$^{-1}$ followed by an extended deceleration over $\sim$700 pc where it terminates. The deceleration is due to both the radially increasing galaxy mass, and mass-loading as ambient gas in the disk is swept up., Comment: A&A accepted; 16 pages

We use a sample of 64 nearby (D=7-45 Mpc) disk galaxies including 45 AGN and 19 non-AGN, that have high spatial resolution multiline CO observations obtained with the ALMA and/or PdBI arrays to study the distribution of cold molecular gas in their circumunuclear disks (CND). We analyze whether the concentration of cold molecular gas changes as a function of the X-ray luminosity in the 2-10 keV range ($L_{\rm X}$). We also study the concentration of the hot molecular gas using NIR data obtained for the H2 1-0S(1) line. We find a turnover in the distribution of the cold molecular gas concentration as a function of $L_{\rm X}$ with a breakpoint which divides the sample into two branches: the AGN build-up branch ($L_{\rm X}\leq10^{41.5\pm0.3}$erg/s) and the AGN feedback branch ($L_{\rm X}\geq10^{41.5\pm0.3}$erg/s) . Lower luminosity AGN and non-AGN of the AGN build-up branch show high cold molecular gas concentrations and centrally peaked radial profiles on nuclear ($r\leq50$~pc) scales. Higher luminosity AGN of the AGN feedback branch, show a sharp decrease in the concentration of molecular gas and flat or inverted radial profiles. The cold molecular gas concentration index ($CCI$), defined as the ratio of surface densities at $r\leq50$~pc and $r\leq200$~pc , namely $CCI \equiv$~log$_{\rm 10}(\Sigma^{\rm gas}_{\rm 50}/\Sigma^{\rm gas}_{\rm 200}$), spans a factor ~4-5 between the galaxies lying at the high end of the AGN build-up branch and the galaxies of the AGN feedback branch. The concentration and radial distributions of the hot molecular gas in our sample follow less extreme trends as a function of the X-ray luminosity. These observations confirm, on a three times larger sample, previous evidence found by the GATOS survey that the imprint of AGN feedback on the CND-scale distribution of molecular gas is more extreme in higher luminosity Seyfert galaxies of the local universe., Comment: Accepted by Astronomy & Astrophysics (A&A) (14/06/2024), 26 pages, 16 figures

Adapters provide an efficient and lightweight mechanism for adapting trained transformer models to a variety of different tasks. However, they have often been found to be outperformed by other adaptation mechanisms, including low-rank adaptation. In this paper, we provide an in-depth study of adapters, their internal structure, as well as various implementation choices. We uncover pitfalls for using adapters and suggest a concrete, improved adapter architecture, called Adapter+, that not only outperforms previous adapter implementations but surpasses a number of other, more complex adaptation mechanisms in several challenging settings. Despite this, our suggested adapter is highly robust and, unlike previous work, requires little to no manual intervention when addressing a novel scenario. Adapter+ reaches state-of-the-art average accuracy on the VTAB benchmark, even without a per-task hyperparameter optimization., Comment: To appear at CVPR 2024. Code: https://github.com/visinf/adapter_plus

The superb image quality, stability and sensitivity of the JWST permit deconvolution techniques to be pursued with a fidelity unavailable to ground-based observations. We present an assessment of several deconvolution approaches to improve image quality and mitigate effects of the complex JWST point spread function (PSF). The optimal deconvolution method is determined by using WebbPSF to simulate JWST's complex PSF and MIRISim to simulate multi-band JWST/Mid-Infrared Imager Module (MIRIM) observations of a toy model of an active galactic nucleus (AGN). Five different deconvolution algorithms are tested: (1) Kraken deconvolution, (2) Richardson-Lucy, (3) Adaptive Imaging Deconvolution Algorithm, (4) Sparse regularization with the Condat-V\~u algorithm, and (5) Iterative Wiener Filtering and Thresholding. We find that Kraken affords the greatest FWHM reduction of the nuclear source of our MIRISim observations for the toy AGN model while retaining good photometric integrity across all simulated wavebands. Applying Kraken to Galactic Activity, Torus, and Outflow Survey (GATOS) multi-band JWST/MIRIM observations of the Seyfert 2 galaxy NGC 5728, we find that the algorithm reduces the FWHM of the nuclear source by a factor of 1.6-2.2 across all five filters. Kraken images facilitate detection of a SE to NW $\thicksim$2".5 ($\thicksim$470 pc, PA $\simeq$115\deg) extended nuclear emission, especially in the longest wavelengths. We demonstrate that Kraken is a powerful tool to enhance faint features otherwise hidden in the complex JWST PSF., Comment: 32 pages, 23 figures, published in AJ 2024 February 7

The usual approach to classify the ionizing source using optical spectroscopy is based on the use of diagnostic diagrams that compares the relative strength of pairs of collisitional metallic lines (e.g., [O iii] and [N ii]) with respect to recombination hydrogen lines (e.g., H{\beta} and H{\alpha}). Despite of being accepted as the standard procedure, it present known problems, including confusion regimes and/or limitations related to the required signal-to-noise of the involved emission lines. These problems affect not only our intrinsic understanding of inter-stellar medium and its poroperties, but also fundamental galaxy properties, such as the star-formation rate and the oxygen abundance, and key questions just as the fraction of active galactic nuclei, among several others. We explore the existing alternatives in the literature to minimize the confusion among different ionizing sources and proposed a new simple diagram that uses the equivalent width and the velocity dispersion from one single emission line, H{\alpha}, to classify the ionizing sources. We use aperture limited and spatial resolved spectroscopic data in the nearby Universe (z{\sim}0.01) to demonstrate that the new diagram, that we called WHaD, segregates the different ionizing sources in a more efficient way that previously adopted procedures. A new set of regions are defined in this diagram to select betweeen different ionizing sources. The new proposed diagram is well placed to determine the ionizing source when only H{\alpha} is available, or when the signal-to-noise of the emission lines involved in the classical diagnostic diagrams (e.g., H{\beta})., Comment: 10 pages, 5 figures, accepted for publishing in A&A

We present MEGARA (Multi-Espectr\'ografo en GTC de Alta Resoluci\'on para Astronom\'ia) Integral Field Unit (IFU) observations of 5 local type-2 quasars (QSO2s, z $\sim 0.1$) from the Quasar Feedback (QSOFEED) sample. These active galactic nuclei (AGN) have bolometric luminosities of 10$^{45.5-46}$ erg/s and stellar masses of $\sim$10$^{11}$ M$_{\odot}$. We explore the kinematics of the ionized gas through the [O~III]$\lambda$5007 $\r{A}$ emission line. The nuclear spectra of the 5 QSO2s, extracted in a circular aperture of $\sim$ 1.2" ($\sim$ 2.2 kpc) in diameter, show signatures of high velocity winds in the form of broad (full width at half maximum; 1300$\leq$FWHM$\leq$2240 km/s and blueshifted components. We find that 4 out of the 5 QSO2s present outflows that we can resolve with our seeing-limited data, and they have radii ranging from 3.1 to 12.6 kpc. In the case of the two QSO2s with extended radio emission, we find that it is well-aligned with the outflows, suggesting that low-power jets might be compressing and accelerating the ionized gas in these radio-quiet QSO2s. In the four QSO2s with spatially resolved outflows, we measure ionized mass outflow rates of 3.3-6.5 Msun/yr when we use [S~II]-based densities, and of 0.7-1.6 Msun/yr when trans-auroral line-based densities are considered instead. We compare them with the corresponding molecular mass outflow rates (8 - 16 Msun/yr), derived from CO(2-1) ALMA observations at 0.2" resolution. Both phases show lower outflow mass rates than those expected from observational scaling relations where uniform assumptions on the outflow properties were adopted. This might be indicating that the AGN luminosity is not the only driver of massive outflows and/or that these relations need to be re-scaled using accurate outflow properties. We do not find a significant impact of the outflows on the global star formation rates., Comment: 24 pages, 13 figures and 4 tables. Accepted for publication in A&A; A&A 665, A55 (2023); doi: 10.1051/0004-6361/202347715

We use JWST/MIRI MRS spectroscopy of a sample of six local obscured type 1.9/2 active galactic nuclei (AGN) to compare their nuclear mid-IR absorption bands with the level of nuclear obscuration traced by X-rays. This study is the first to use sub-arcsecond angular resolution data of local obscured AGN to investigate the nuclear mid-IR absorption bands with a wide wavelength coverage (4.9-28.1 $\mu$m). All the nuclei show the 9.7 $\mu$m silicate band in absorption. We compare the strength of the 9.7 and 18 $\mu$m silicate features with torus model predictions. The observed silicate features are generally well explained by clumpy and smooth torus models. We report the detection of the 6 $\mu$m dirty water ice band (i.e., a mix of water and other molecules such as CO and CO$_2$) at sub-arcsecond scales ($\sim$0.26 arcsec at 6 $\mu$m; inner $\sim$50 pc) in a sample of local AGN with different levels of nuclear obscuration in the range log N$_{\rm H}^{\rm X-Ray}$(cm$^{-2}$)$\sim22-25$. We find a good correlation between the 6 $\mu$m water ice optical depths and N$_{\rm H}^{\rm X-Ray}$. This result indicates that the water ice absorption might be a reliable tracer of the nuclear intrinsic obscuration in AGN. The weak water ice absorption in less obscured AGN (log N$_H^{X-ray}$ (cm$^{-2}$)$\lesssim$23.0 cm$^{-2}$) might be related to the hotter dust temperature ($>$T$_{sub}^{H_2O}\sim$110 K) expected to be reached in the outer layers of the torus due to their more inhomogeneous medium. Our results suggest it might be necessary to include the molecular content, such as, H$_2$O, aliphatic hydrocarbons (CH-) and more complex PAH molecules in torus models to better constrain key parameters such as the torus covering factor (i.e. nuclear obscuration)., Comment: Accepted for publication in A&A. 13 pages, 12 Figures

Introduction. There is currently no guidance on how to assess the calibration of multistate models used for risk prediction. We introduce several techniques that can be used to produce calibration plots for the transition probabilities of a multistate model, before assessing their performance in the presence of non-informative and informative censoring through a simulation. Methods. We studied pseudo-values based on the Aalen-Johansen estimator, binary logistic regression with inverse probability of censoring weights (BLR-IPCW), and multinomial logistic regression with inverse probability of censoring weights (MLR-IPCW). The MLR-IPCW approach results in a calibration scatter plot, providing extra insight about the calibration. We simulated data with varying levels of censoring and evaluated the ability of each method to estimate the calibration curve for a set of predicted transition probabilities. We also developed evaluated the calibration of a model predicting the incidence of cardiovascular disease, type 2 diabetes and chronic kidney disease among a cohort of patients derived from linked primary and secondary healthcare records. Results. The pseudo-value, BLR-IPCW and MLR-IPCW approaches give unbiased estimates of the calibration curves under non-informative censoring. These methods remained unbiased in the presence of informative censoring, unless the mechanism was strongly informative, with bias concentrated in the areas of predicted transition probabilities of low density. Conclusions. We recommend implementing either the pseudo-value or BLR-IPCW approaches to produce a calibration curve, combined with the MLR-IPCW approach to produce a calibration scatter plot, which provides additional information over either of the other methods., Comment: Pre-print for article currently under review

We present a polarization variability analysis of a sample of 26 $\gamma$-ray blazars monitored by the Steward Observatory between 2008 and 2018 in the optical band. We investigate the properties and long-term variability of their optical polarization, searching for differences between blazar types. We observe that BL Lac objects are typically less polarized and less variable than flat spectrum radio quasars (FSRQs). Moreover, BL Lacs display a distribution of their polarization angle typically oriented in a preferential direction, contrary to the rather random distribution of FSRQs. For the latter blazar type, as well as those sources showing a bright stellar emission, we take into account the depolarizing effect introduced by the broad line region and the host galaxy on the measured polarization degree. In this sample we also observe that BL Lacs present an uncorrelated evolution of the flux and the polarization. Contrary, FSRQs show a correlation before the depolarization correction, that is lost however after considering this effect. In addition, we study the behaviour of the polarization angle, searching for angle rotations in its long-term evolution. We derive that the FSRQs studied here show rotations more frequently than BL Lac objects by a factor $\sim$1.5. During these periods we also observe a systematic decrease of the polarization fraction, as well as a marginal flux increase, not significant however to connect rotations with optical flares. We interpret these results within the extended shock-in-jet scenario, able to explain the overall features observed here for the polarization of the blazar sample., Comment: Accepted for publication in MNRAS

We present the analysis of archival XMM-Newton and Chandra observations of CH Cyg, one of the most studied symbiotic stars (SySts). The combination of the high-resolution XMM-Newton RGS and Chandra HETG X-ray spectra allowed us to obtain reliable estimates of the chemical abundances and to corroborate the presence of multi-temperature X-ray-emitting gas. Spectral fitting of the medium-resolution XMM-Newton MOS (MOS1+MOS2) spectrum required the use of an additional component not seen in previous studies in order to fit the 2.0-4.0 keV energy range. Detailed spectral modelling of the XMM-Newton MOS data suggests the presence of a reflection component, very similar to that found in active galactic nuclei. The reflection component is very likely produced by an ionised disk (the accretion disk around the white dwarf) and naturally explains the presence of the fluorescent Fe emission line at 6.4 keV while also contributing to the soft and medium energy ranges. The variability of the global X-ray properties of CH Cyg are discussed as well as the variation of the three Fe lines around the 6-7 keV energy range. We conclude that reflection components are needed to model the hard X-ray emission and may be present in most $\beta/\delta$-type SySt., Comment: 12 pages, 7 figures, 5 tables; To appear in MNRAS

Engineered dissipative reservoirs have the potential to steer many-body quantum systems toward correlated steady states useful for quantum simulation of high-temperature superconductivity or quantum magnetism. Using up to 49 superconducting qubits, we prepared low-energy states of the transverse-field Ising model through coupling to dissipative auxiliary qubits. In one dimension, we observed long-range quantum correlations and a ground-state fidelity of 0.86 for 18 qubits at the critical point. In two dimensions, we found mutual information that extends beyond nearest neighbors. Lastly, by coupling the system to auxiliaries emulating reservoirs with different chemical potentials, we explored transport in the quantum Heisenberg model. Our results establish engineered dissipation as a scalable alternative to unitary evolution for preparing entangled many-body states on noisy quantum processors.

Undesired coupling to the surrounding environment destroys long-range correlations on quantum processors and hinders the coherent evolution in the nominally available computational space. This incoherent noise is an outstanding challenge to fully leverage the computation power of near-term quantum processors. It has been shown that benchmarking Random Circuit Sampling (RCS) with Cross-Entropy Benchmarking (XEB) can provide a reliable estimate of the effective size of the Hilbert space coherently available. The extent to which the presence of noise can trivialize the outputs of a given quantum algorithm, i.e. making it spoofable by a classical computation, is an unanswered question. Here, by implementing an RCS algorithm we demonstrate experimentally that there are two phase transitions observable with XEB, which we explain theoretically with a statistical model. The first is a dynamical transition as a function of the number of cycles and is the continuation of the anti-concentration point in the noiseless case. The second is a quantum phase transition controlled by the error per cycle; to identify it analytically and experimentally, we create a weak link model which allows varying the strength of noise versus coherent evolution. Furthermore, by presenting an RCS experiment with 67 qubits at 32 cycles, we demonstrate that the computational cost of our experiment is beyond the capabilities of existing classical supercomputers, even when accounting for the inevitable presence of noise. Our experimental and theoretical work establishes the existence of transitions to a stable computationally complex phase that is reachable with current quantum processors.

Access to knowledge of the point spread function (PSF) of adaptive optics(AO)-assisted observations is still a major limitation when processing AO data. This limitation is particularly important when image analysis requires the use of deconvolution methods. As the PSF is a complex and time-varying function, reference PSFs acquired on calibration stars before or after the scientific observation can be too different from the actual PSF of the observation to be used for deconvolution, and lead to artefacts in the final image. We improved the existing PSF-estimation method based on the so-called marginal approach by enhancing the object prior in order to make it more robust and suitable for observations of resolved extended objects. Our process is based on a two-step blind deconvolution approach from the literature. The first step consists of PSF estimation from the science image. For this, we made use of an analytical PSF model, whose parameters are estimated based on a marginal algorithm. This PSF was then used for deconvolution. In this study, we first investigated the requirements in terms of PSF parameter knowledge to obtain an accurate and yet resilient deconvolution process using simulations. We show that current marginal algorithms do not provide the required level of accuracy, especially in the presence of small objects. Therefore, we modified the marginal algorithm by providing a new model for object description, leading to an improved estimation of the required PSF parameters. Our method fulfills the deconvolution requirement with realistic system configurations and different classes of Solar System objects in simulations. Finally, we validate our method by performing blind deconvolution with SPHERE/ZIMPOL observations of the Kleopatra asteroid., Comment: accepted for publication in Astronomy & Astrophysics (A&A)

Mrk 477 is the closest type-2 quasar (QSO2), at a distance of 163 Mpc. This makes it an ideal laboratory for studying the interplay between nuclear activity and star formation with a great level of detail and signal-to-noise. In this Letter we present new mid-infrared (mid-IR) imaging and spectroscopic data with an angular resolution of 0.4 arcsec (~300 pc) obtained with the Gran Telescopio Canarias (GTC) instrument CanariCam. The N-band (8-13 micron) spectrum of the central ~400 pc of the galaxy reveals [S IV]10.51 micron emission, but no 8.6 or 11.3 micron PAH features, which are commonly used as tracers of recent star formation. This is in stark contrast with the presence of a nuclear starburst of ~300 pc in size, an age of 6 Myr, and a mass of 1.1x10^8 Msun, as constrained from ultraviolet Hubble Space Telescope observations. Considering this, we argue that even the more resilient, neutral molecules that mainly produce the 11.3 micron PAH band are most likely being destroyed in the vicinity of the active nucleus despite the relatively large X-ray column density, of log N_H=23.5 cm^-2, and modest X-ray luminosity, of 1.5x10^43 erg/s. This highlights the importance of being cautious when using PAH features as star formation tracers in the central region of galaxies to evaluate the impact of feedback from active galactic nuclei., Comment: 7 pages, 3 figures. A&A Letters accepted version including proof corrections

An interesting feature in active galactic nuclei (AGN) accreting at low rate is the weakness of the reflection features in their X-ray spectra, which can result from the gradual disappearance of the torus with decreasing accretion rates. It has been suggested that low luminosity AGN (LLAGN) would have a different reflector configuration compared with high luminosity AGN, either covering a smaller fraction of the sky or simply having less material. Additionally, we note that the determination of the spectral index ($\Gamma$) and the cut-off energy of the primary power-law emission is affected by the inclusion of reflection models, showing the importance of using them to study the accretion mechanism, especially in the case of the LLAGN that have previously shown a high dispersion on the relation between $\Gamma$ and the accretion rate. Our purpose is to constrain the geometry and column density of the reflector in a sample of LLAGN covering a broad X-ray range of energy combining data from XMM-Newton + NuSTAR + Swift of a hard X-ray-flux limited sample of 17 LLAGN from BASS/DR2 with accretion rates $\lambda_{Edd}$=L$_{\rm Bol}$/L$_{\rm Edd}$<10$^{-3}$. We fit all spectra using the reflection model for torus (borus02) and accretion disk (Xillver) reflectors. We found a tentative correlation between the torus column density and the accretion rate, LLAGN shows a lower column density compared with the high-luminosity objects. We also confirm the relation between $\Gamma$ and $\lambda_{Edd}$, with a smaller scatter than previously reported, thanks to the inclusion of high-energy data and the reflection models. Our results are consistent with a break at $\lambda_{Edd}\sim10^{-3}$, suggestive of a different accretion mechanism compared with higher accretion AGN., Comment: Accepted for publication in Astronomy & Astrophysics

Characterizing the PSF of adaptive optics instruments is of paramount importance both for instrument design and observation planning/optimization. Simulation software, such as PASSATA, have been successfully utilized for PSF characterization in instrument design, which make use of standardized atmospheric turbulence profiles to produce PSFs that represent the typical instrument performance. In this contribution we study the feasibility of using such tool for nowcast application (present-time forecast), such as the characterization of an on-sky measured PSF in real observations. Specifically we will analyze the performance of the simulation software in characterizing the real-time PSF of two different state-of-the-art SCAO adaptive optics instruments: SOUL at the LBT, and SAXO at the VLT. The study will make use of on-sky measurements of the atmospheric turbulence and compare the results of the simulations to the measured PSF figures of merit (namely the FHWM and the Strehl Ratio) retrieved from the instrument telemetry in real observations. Our main goal in this phase is to quantify the level of uncertainly of the AO simulations in reproducing real on-sky observed PSFs with an end-to-end code (PASSATA). In a successive phase we intend to use a faster analytical code (TIPTOP). This work is part of a wider study which aims to use simulation tools joint to atmospheric turbulence forecasts performed nightly to forecast in advance the PSF and support science operations of ground-based telescopes facilities. The 'PSF forecast' option might therefore be added to ALTA Center or the operational forecast system that will be implemented soon at ESO., Comment: 10 pages, 6 figures

An important measure of the development of quantum computing platforms has been the simulation of increasingly complex physical systems. Prior to fault-tolerant quantum computing, robust error mitigation strategies are necessary to continue this growth. Here, we study physical simulation within the seniority-zero electron pairing subspace, which affords both a computational stepping stone to a fully correlated model, and an opportunity to validate recently introduced ``purification-based'' error-mitigation strategies. We compare the performance of error mitigation based on doubling quantum resources in time (echo verification) or in space (virtual distillation), on up to $20$ qubits of a superconducting qubit quantum processor. We observe a reduction of error by one to two orders of magnitude below less sophisticated techniques (e.g. post-selection); the gain from error mitigation is seen to increase with the system size. Employing these error mitigation strategies enables the implementation of the largest variational algorithm for a correlated chemistry system to-date. Extrapolating performance from these results allows us to estimate minimum requirements for a beyond-classical simulation of electronic structure. We find that, despite the impressive gains from purification-based error mitigation, significant hardware improvements will be required for classically intractable variational chemistry simulations., Comment: 10 pages, 13 page supplementary material, 12 figures. Experimental data available at https://doi.org/10.5281/zenodo.7225821

We present a statistical analysis of the properties of the obscuring material around active galactic nuclei (AGN). This study represents the first of its kind for an ultra-hard X-ray (14-195keV; Swift/BAT) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample) using high angular resolution infrared data and various torus models: smooth, clumpy and two-phase torus models and clumpy disc+wind models. We find that the smooth, clumpy and two-phase torus models (i.e. without including the polar dusty wind component) and disc+wind models provide best fits for a comparable number of galaxies, 8/24 (33.3%) and 9/24 (37.5%), respectively. We find that the best-fitted models depend on the hydrogen column density (NH), which is related to the X-ray (unobscured/obscured) and/or optical (Sy1/Sy2) classification. In particular, smooth, clumpy and two-phase torus models best reproduce the infrared (IR) emission of AGN with relatively high hydrogen column density (median value of log (NH)=23.5+-0.8; i.e. Sy2). However, clumpy disc+wind models provide the best fits to the nuclear IR spectral energy distributions (SEDs) of Sy1/1.8/1.9 (median value of log (NH)=21.0+-1.0), specifically in the near-IR (NIR) range. The success of the disc+wind models in fitting the NIR emission of Sy1 galaxies is due to the combination of adding large graphite grains to the dust composition and self-obscuration effects caused by the wind at intermediate inclinations. In general, we find that the Seyfert galaxies having unfavourable (favourable) conditions, i.e. nuclear hydrogen column density and Eddington ratio, for launching IR dusty polar outflows are best-fitted with smooth, clumpy and two-phase torus (disk+wind) models confirming the predictions from simulations. Therefore, our results indicate that the nature of the inner dusty structure in AGN depend on the intrinsic AGN properties., Comment: Accepted for publication in A&A, 43 pages, 43 figures, 16 tables

The Multi Unit Spectroscopic Explorer (MUSE) is an integral field spectrograph on the Very Large Telescope Unit Telescope 4, capable of laser guide star assisted and tomographic adaptive optics using the GALACSI module. Its observing capabilities include a wide field (1 square arcmin), ground layer AO mode (WFM-AO) and a narrow field (7.5"x7.5"), laser tomography AO mode (NFM-AO). The latter has had several upgrades in the 4 years since commissioning, including an optimisation of the control matrices for the AO system and a new sub-electron noise detector for its infra-red low order wavefront sensor. We set out to quantify the NFM-AO system performance by analysing $\sim$230 spectrophotometric standard star observations taken over the last 3 years. To this end we expand upon previous work, designed to facilitate analysis of the WFM-AO system performance. We briefly describe the framework that will provide a user friendly, semi-automated way for system performance monitoring during science operations. We provide the results of our performance analysis, chiefly through the measured Strehl ratio and full width at half maximum (FWHM) of the core of the point spread function (PSF) using two PSF models, and correlations with atmospheric conditions. These results will feed into a range of applications, including providing a more accurate prediction of the system performance as implemented in the exposure time calculator, and the associated optimization of the scientific output for a given set of limiting atmospheric conditions., Comment: SPIE proceedings (2022), Observatory Operations: Strategies, Processes, and Systems IX

We have studied the nuclear region of the previously detected dual AGN system in the galaxy pair IRAS 05589+2828 and 2MASX J06021107+2828382 through new optical spectroscopy observations, along with radio and X-ray archival data. Our multiwavelength data strongly suggest that the Sy1 \iras\, (z=0.0330$\pm$0.0002) conforms to a dual AGN system with the Sy2 \twomas\, (z=0.0334$\pm$0.0001) with a projected separation obtained from the radio data of 20.08\arcsec\, ($\sim$13.3\,kpc). Analysis of the optical spectra reveals a faint narrow extended emission from H$\alpha$ and [OIII] amidst the two AGN, supporting evidence for an ongoing merger. \iras\, is a double component narrow emission line AGN, with complex broad Balmer emission line profiles that clearly show a strong red-peaklet with a velocity shift of $\sim$3500\,km\,s$^{-1}$. The black hole mass estimates of \iras\, and \twomas\, are log\,M$\rm_{BH}$\,=\,8.59\,$\pm$\,0.14 (M$_\odot$) and log\,M$\rm_{BH}$\,=\,8.21$\pm$0.2 (M$_\odot$), respectively. In the X-ray bands, \iras\, is compatible with a Type 1 object, showing both spectral and flux variability. \chandra\, data of 2MASX\,J06021107+2828382 allowed us to measure a high hardness ratio in this source, providing evidence for a Type 2 AGN. The 22 GHz image obtained with the {\it Karl G. Jansky Very Large Array} has revealed that both AGN are compact radio objects with spectral indices -0.26$\pm$0.03 and -0.70$\pm$0.11, confirming for the first time its dual AGN nature in the radio bands., Comment: Accepted for publication in MNRAS. 20 pages, 15 figures

Observations across the heliosphere typically rely on in situ spacecraft observations producing time-series data. While often this data is analysed visually, it lends itself more naturally to our sense of sound. The simplest method of converting oscillatory data into audible sound is audification -- a one-to-one mapping of data samples to audio samples -- which has the benefit that no information is lost, thus is a true representation of the original data. However, audification can make some magnetospheric ULF waves observations pass by too quickly for someone to realistically be able to listen to effectively. For this reason, we detail various existing audio time scale modification techniques developed for music, applying these to ULF wave observations by spacecraft and exploring how they affect the properties of the resulting audio. Through a public dialogue we arrive at recommendations for ULF wave researchers on rendering these waves audible and discuss the scientific and educational possibilities of these new methods.

Educational research highlights that improved careers education is needed to increase participation in science, technology, engineering, and mathematics (STEM). Current UK careers resources concerning the space sector, however, are found to perhaps not best reflect the diversity of roles present and may in fact perpetuate misconceptions about the usefulness of science. We, therefore, compile a more diverse set of space-related jobs, which will be used in the development of a new space careers resource.

To maximize the economic benefits of geothermal energy production, it is essential to optimize geothermal reservoir management strategies, in which geologic uncertainty should be considered. In this work, we propose a closed-loop optimization framework, based on deep learning surrogates, for the well control optimization of geothermal reservoirs. In this framework, we construct a hybrid convolution-recurrent neural network surrogate, which combines the convolution neural network (CNN) and long short-term memory (LSTM) recurrent network. The convolution structure can extract spatial information of geologic parameter fields and the recurrent structure can approximate sequence-to-sequence mapping. The trained model can predict time-varying production responses (rate, temperature, etc.) for cases with different permeability fields and well control sequences. In the closed-loop optimization framework, production optimization based on the differential evolution (DE) algorithm, and data assimilation based on the iterative ensemble smoother (IES), are performed alternately to achieve real-time well control optimization and geologic parameter estimation as the production proceeds. In addition, the averaged objective function over the ensemble of geologic parameter estimations is adopted to consider geologic uncertainty in the optimization process. Several geothermal reservoir development cases are designed to test the performance of the proposed production optimization framework. The results show that the proposed framework can achieve efficient and effective real-time optimization and data assimilation in the geothermal reservoir production process., Comment: 37 pages, 24 figures

Blazars optical emission is generally dominated by relativistic jets, although the host galaxy, accretion disk and broad line region (BLR) may also contribute significantly. Disentangling their contributions has been challenging for years due to the dominance of the jet. To quantify the contributions to the spectral variability, we use the statistical technique for dimensionality reduction Non-Negative Matrix Factorization on a spectroscopic data set of 26 $\gamma$-ray blazars. This technique allows to model large numbers of spectra in terms of a reduced number of components.We use a priori knowledge to obtain components associated to meaningful physical processes. The sources are classified according to their optical spectrum as host-galaxy dominated BL Lac objects (BL Lacs), BL Lacs, or Flat Spectrum Radio Quasars (FSRQs). Host-galaxy sources show less variability, as expected, and bluer-when-brighter trends, as the other BL Lacs. For FSRQs, more complicated colour-flux behaviours are observed: redder-when-brighter for low states saturating above a certain level and, in some cases, turning to bluer-when-brighter. We are able to reproduce the variability observed during 10 years using only 2 to 4 components, depending on the type. The simplest scenario corresponds to host-galaxy blazars, whose spectra are reconstructed using the stellar population and a power law for the jet. BL Lac spectra are reproduced using from 2 to 4 power laws. Different components can be associated to acceleration/cooling processes taking place in the jet. The reconstruction of FSRQs also incorporates a QSO-like component to account for the BLR, plus a very steep power law, associated to the accretion disk., Comment: Accepted for publication in MNRAS

Recent advances in multimodal vision and language modeling have predominantly focused on the English language, mostly due to the lack of multilingual multimodal datasets to steer modeling efforts. In this work, we address this gap and provide xGQA, a new multilingual evaluation benchmark for the visual question answering task. We extend the established English GQA dataset to 7 typologically diverse languages, enabling us to detect and explore crucial challenges in cross-lingual visual question answering. We further propose new adapter-based approaches to adapt multimodal transformer-based models to become multilingual, and -- vice versa -- multilingual models to become multimodal. Our proposed methods outperform current state-of-the-art multilingual multimodal models (e.g., M3P) in zero-shot cross-lingual settings, but the accuracy remains low across the board; a performance drop of around 38 accuracy points in target languages showcases the difficulty of zero-shot cross-lingual transfer for this task. Our results suggest that simple cross-lingual transfer of multimodal models yields latent multilingual multimodal misalignment, calling for more sophisticated methods for vision and multilingual language modeling., Comment: Findings of ACL 2022

Reasoning over multiple modalities, e.g. in Visual Question Answering (VQA), requires an alignment of semantic concepts across domains. Despite the widespread success of end-to-end learning, today's multimodal pipelines by and large leverage pre-extracted, fixed features from object detectors, typically Faster R-CNN, as representations of the visual world. The obvious downside is that the visual representation is not specifically tuned to the multimodal task at hand. At the same time, while transformer-based object detectors have gained popularity, they have not been employed in today's multimodal pipelines. We address both shortcomings with TxT, a transformer-based crossmodal pipeline that enables fine-tuning both language and visual components on the downstream task in a fully end-to-end manner. We overcome existing limitations of transformer-based detectors for multimodal reasoning regarding the integration of global context and their scalability. Our transformer-based multimodal model achieves considerable gains from end-to-end learning for multimodal question answering., Comment: To appear at the 43rd DAGM German Conference on Pattern Recognition (GCPR) 2021

To study the nuclear ($\lesssim1\,$kpc) dust of nearby ($z<0.1$) type 1 Quasi Stellar Objects (QSOs) we obtained new near-infrared (NIR) high angular resolution ($\sim0.3$ arcsec) photometry in the H and Ks bands, for 13 QSOs with available mid-infrared (MIR) high angular resolution spectroscopy ($\sim7.5-13.5\,\mu$m). We find that in most QSOs the NIR emission is unresolved. We subtract the contribution from the accretion disk, which decreases from NIR ($\sim35\%$) to MIR ($\sim2.4\%$). We also estimate these percentages assuming a bluer accretion disk and find that the contibution in the MIR is nearly seven time larger. We find that the majority of objects ($64\%$, 9/13) are better fitted by the Disk+Wind H17 model \citep[][]{Hoenig17}, while others can be fitted by the Smooth F06 \citep[$14\%$, 2/13,][]{Fritz06}, Clumpy N08 \citep[$7\%$, 1/13,][]{Nenkova08a,Nenkova08b}, Clumpy H10 \citep[$7\%$, 1/13,][]{Hoenig10b}, and Two-Phase media S16 \citep[$7\%$, 1/13,][]{Stalev16} models. However, if we assume the bluer accretion disk, the models fit only 2/13 objects. We measured two NIR to MIR spectral indexes, $\alpha_{NIR-MIR(1.6,8.7\,\mu\text{m})}$ and $\alpha_{NIR-MIR(2.2,8.7\,\mu\text{m})}$, and two MIR spectral indexes, $\alpha_{MIR(7.8, 9.8\,\mu\text{m})}$ and $\alpha_{MIR(9.8, 11.7\,\mu\text{m})}$, from models and observations. From observations, we find that the NIR to MIR spectral indexes are $\sim-1.1$ and the MIR spectral indexes are $\sim-0.3$. Comparing the synthetic and observed values, we find that none of the models simultaneously match the measured NIR to MIR and $7.8-9.8\,\mu$m slopes. However, we note that measuring the $\alpha_{MIR(7.8, 9.8\,\mu\text{m})}$ on the starburst-subtracted {\it Spitzer}/IRS spectrum, gives values of the slopes ($\sim-2$) that are similar to the synthetic values obtained from the models., Comment: Accepted in ApJ

We study the problem of deconvolution for light-sheet microscopy, where the data is corrupted by spatially varying blur and a combination of Poisson and Gaussian noise. The spatial variation of the point spread function (PSF) of a light-sheet microscope is determined by the interaction between the excitation sheet and the detection objective PSF. First, we introduce a model of the image formation process that incorporates this interaction, therefore capturing the main characteristics of this imaging modality. Then, we formulate a variational model that accounts for the combination of Poisson and Gaussian noise through a data fidelity term consisting of the infimal convolution of the single noise fidelities, first introduced in L. Calatroni et al. "Infimal convolution of data discrepancies for mixed noise removal", SIAM Journal on Imaging Sciences 10.3 (2017), 1196-1233. We establish convergence rates in a Bregman distance under a source condition for the infimal convolution fidelity and a discrepancy principle for choosing the value of the regularisation parameter. The inverse problem is solved by applying the primal-dual hybrid gradient (PDHG) algorithm in a novel way. Finally, numerical experiments performed on both simulated and real data show superior reconstruction results in comparison with other methods., Comment: 34 pages, 13 figures

The magnetopause marks the outer edge of the Earth's magnetosphere and a distinct boundary between solar wind and magnetospheric plasma populations. In this letter, we use global magnetohydrodynamic simulations to examine the response of the terrestrial magnetopause to fast-forward interplanetary shocks of various strengths and compare to theoretical predictions. The theory and simulations indicate the magnetopause response can be characterised by three distinct phases; an initial acceleration as inertial forces are overcome, a rapid compressive phase comprising the majority of the distance travelled, and large-scale damped oscillations with amplitudes of the order of an Earth radius. The two approaches agree in predicting subsolar magnetopause oscillations with frequencies 2-13 mHz but the simulations notably predict larger amplitudes and weaker damping rates. This phenomenon is of high relevance to space weather forecasting and provides a possible explanation for magnetopause oscillations observed following the large interplanetary shocks of August 1972 and March 1991., Comment: 9 pages, 3 figures, 1 table. Accepted as a Geophysical Research Letter on 09 July 2021

We compare mid-IR and ALMA far-IR images of 12 nearby Seyferts selected from GATOS. The mid-IR unresolved emission contributes more than 60% of the nuclear emission in most galaxies. By contrast, the ALMA 870micron continuum emission is mostly resolved and typically along the torus equatorial direction (Paper I, Garcia-Burillo et al. 2021). The Eddington ratios and nuclear hydrogen column densities NH of half the sample are favorable to launching polar and/or equatorial dusty winds, according to simulations. Six show mid-IR extended emission in the polar direction as traced by the NLR and perpendicular to the ALMA emission. In a few, the nuclear NH might be too high to uplift large quantities of dusty material along the polar direction. Five galaxies have low NH and/or Eddington ratios and thus polar dusty winds are not likely. We generate new CAT3D-WIND disk-wind model images. At low wind-to-disk cloud ratios the far-IR model images have disk- and ring-like morphologies. The X-shape associated with dusty winds is seen better in the far-IR at intermediate-high inclinations for the extended-wind configurations. In most models, the mid-IR emission comes from the inner part of the disk/cone. Extended bi-conical and one-sided polar mid-IR emission is seen in extended-wind configurations and high wind-to-disk cloud ratios. When convolved to our resolution, the model images reproduce qualitative aspects of the observed morphologies. Low-intermediate wind-to-disk ratios are required to account for the large fractions of unresolved mid-IR emission. This work and Paper I provide observational support for the torus+wind scenario. The wind component is more relevant at high Eddington ratios and/or AGN luminosities, and polar dust emission is predicted at NH of up to $10^{24}$cm$^{-2}$. The torus/disk component, on the other hand, prevails at low luminosities and/or Eddington ratios. (Abridged), Comment: 27 pages, 27 figures. Accepted for publication in A&A

This paper presents a deep learning-based de-homogenization method for structural compliance minimization. By using a convolutional neural network to parameterize the mapping from a set of lamination parameters on a coarse mesh to a one-scale design on a fine mesh, we avoid solving the least square problems associated with traditional de-homogenization approaches and save time correspondingly. To train the neural network, a two-step custom loss function has been developed which ensures a periodic output field that follows the local lamination orientations. A key feature of the proposed method is that the training is carried out without any use of or reference to the underlying structural optimization problem, which renders the proposed method robust and insensitive wrt. domain size, boundary conditions, and loading. A post-processing procedure utilizing a distance transform on the output field skeleton is used to project the desired lamination widths onto the output field while ensuring a predefined minimum length-scale and volume fraction. To demonstrate that the deep learning approach has excellent generalization properties, numerical examples are shown for several different load and boundary conditions. For an appropriate choice of parameters, the de-homogenized designs perform within $7-25\%$ of the homogenization-based solution at a fraction of the computational cost. With several options for further improvements, the scheme may provide the basis for future interactive high-resolution topology optimization., Comment: 28 pages, 16 figures

We present the first results of the Galaxy Activity, Torus and Outflow Survey (GATOS), a project aimed at understanding the properties of the dusty molecular tori and their connection to the host galaxy in nearby Seyfert galaxies. Our project expands the range of AGN luminosities and Eddington ratios covered by previous surveys of Seyferts conducted by ALMA and allows us to study the gas feeding and feedback cycle in a combined sample of 19 Seyferts. We used ALMA to obtain new images of the emission of molecular gas and dust using the CO(3-2) and HCO+(4-3) lines as well as their underlying continuum emission at 870 microns with high spatial resolutions (0.1'' ~ 7 - 13 pc) in the CND of 10 nearby (D < 28 Mpc) Seyfert galaxies. Our new ALMA observations detect 870 micron continuum and CO line emission from spatially resolved disks located around the AGN in all the sources. The bulk of the continuum flux can be accounted for by thermal emission from dust in the majority of the targets. For most of the sources the disks show a preponderant orientation perpendicular to the AGN wind axes, as expected for dusty molecular tori. The median diameters and molecular gas masses of the tori are ~ 42 pc, and ~ 6 x 10**5 Msun, respectively. We find a positive correlation between the line-of-sight gas column densities responsible for the absorption of X-rays and the molecular gas column densities derived from CO towards the AGN in our sources. The radial distributions of molecular gas in the CND of our combined sample show signs of nuclear-scale molecular gas deficits. We also detect molecular outflows in the sources that show the most extreme nuclear-scale gas deficits in our sample. These observations find for the first time supporting evidence that the imprint of AGN feedback is more extreme in higher luminosity and/or higher Eddington ratio Seyfert galaxies., Comment: 41 pages, new version revised and accepted by A&A on the 22nd of June 2021

Societal biases are a major issue in school students' access to and interaction with science. Schools engagement programmes in science from universities, like independent research projects, which could try and tackle these problems are, however, often inequitable. We evaluate these concerns applied to one such programme, `Physics Research in School Environments' (PRiSE), which features projects in space science, astronomy, and particle physics. Comparing the schools involved with PRiSE to those of other similar schemes and UK national statistics, we find that PRiSE has engaged a much more diverse set of schools with significantly more disadvantaged groups than is typical. While drop-off occurs within the protracted programme, we find no evidence of systematic biases present. The majority of schools that complete projects return for multiple years of the programme, with this repeated buy-in from schools again being unpatterned by typical societal inequalities. Therefore, schools' ability to succeed at independent research projects appears independent of background within the PRiSE framework. Qualitative feedback from teachers show that the diversity and equity of the programme, which they attribute to the level of support offered through PRiSE's framework, is valued and they have highlighted further ways of making the projects potentially even more accessible. Researcher-involvement, uncommon in many other programmes, along with teacher engagement and communication are found to be key elements to success in independent research projects overall.

Physics in schools is distinctly different from, and struggles to capture the excitement of, university research-level work. Initiatives where students engage in independent research linked to cutting-edge physics within their school over several months might help mitigate this, potentially facilitating the uptake of science in higher education. However, how such initiatives are best supported remains unclear and understudied. This paper evaluates a provision framework, `Physics Research in School Environments' (PRiSE), using survey data from participating 14-18 year-old students and their teachers to understand their experience of the programme. The results show that PRiSE appears to provide much more positive experiences than typical university outreach initiatives due to the nature of the opportunities afforded over several months, which schools would not be able to provide without external input. The intensive support offered is deemed necessary, with all elements appearing equally important. Based on additional feedback from independent researchers and engagement professionals, we also suggest the framework could be adopted at other institutions and applied to their own areas of scientific research, something which has already started to occur.

Using 6 years of evaluation data we assess the medium- and long-term impacts upon a diverse range of students, teachers, and schools from participating in a programme of protracted university-mentored projects based in cutting-edge space science, astronomy, and particle physics research. After having completed their 6-month-long projects, the 14-18 year-old school students report having substantially increased in confidence relating to relevant scientific topics and methods as well as having developed numerous skills, outcomes which are corroborated by teachers. There is evidence that the projects helped increase students' aspirations towards physics, whereas science aspirations (generally high to begin with) were typically maintained or confirmed through their involvement. Longitudinal evaluation 3 years later has revealed that these projects have been lasting experiences for students which they have benefited and drawn upon in their subsequent university education. Data on students' destinations suggests that their involvement in research projects has made them more likely to undertake physics and STEM degrees than would otherwise be expected. Cases of co-created novel physics research resulting from PRiSE also has seemed to have a powerful effect, not only on the student co-authors but participating students from other schools also. Teachers have also been positively affected through participating, with the programme having influenced their own knowledge, skills, and pedagogy, as well as having advantageous effects felt across their wider schools. These impacts suggest that similar `research in schools' initiatives may have a role to play in aiding the increased uptake and diversity of physics and/or STEM in higher education as well as meaningfully enhancing the STEM environment within schools.

The discovery of topological order has revolutionized the understanding of quantum matter in modern physics and provided the theoretical foundation for many quantum error correcting codes. Realizing topologically ordered states has proven to be extremely challenging in both condensed matter and synthetic quantum systems. Here, we prepare the ground state of the toric code Hamiltonian using an efficient quantum circuit on a superconducting quantum processor. We measure a topological entanglement entropy near the expected value of $\ln2$, and simulate anyon interferometry to extract the braiding statistics of the emergent excitations. Furthermore, we investigate key aspects of the surface code, including logical state injection and the decay of the non-local order parameter. Our results demonstrate the potential for quantum processors to provide key insights into topological quantum matter and quantum error correction., Comment: 6 pages 4 figures, plus supplementary materials

In this article we provide formulations of energy flux and radiation stress consistent with the scaling regime of the Korteweg-de Vries (KdV) equation. These quantities can be used to describe the shoaling of cnoidal waves approaching a gently sloping beach. The transformation of these waves along the slope can be described using the shoaling equations, a set of three nonlinear equations in three unknowns: the wave height H, the set-down and the elliptic parameter m. We define a numerical algorithm for the efficient solution of the shoaling equations, and we verify our shoaling formulation by comparing with experimental data from two sets of experiments as well as shoaling curves obtained in previous works., Comment: 19 pages, 7 figures

Background SAMtools and BCFtools are widely used programs for processing and analysing high-throughput sequencing data. Findings The first version appeared online twelve years ago and has been maintained and further developed ever since, with many new features and improvements added over the years. The SAMtools and BCFtools packages represent a unique collection of tools that have been used in numerous other software projects and countless genomic pipelines. Conclusion Both SAMtools and BCFtools are freely available on GitHub under the permissive MIT licence, free for both non-commercial and commercial use. Both packages have been installed over a million times via Bioconda. The source code and documentation are available from http://www.htslib.org.

The optimization and scheduling of scientific observations done with instrumentation supported by adaptive optics could greatly benefit from the forecast of PSF figures of merit (FWHM, Strehl Ratio, Encircle Energy and contrast), that depend on the AO instrument, the scientific target and turbulence conditions during the observing night. In this contribution we explore the the possibility to forecast a few among the most useful PSF figures of merit (SR and FWHM). To achieve this goal, we use the optical turbulence forecasted by the mesoscale atmospheric model Astro-Meso-NH on a short timescale as an input for PSF simulation software developed and tailored for specific AO instruments. A preliminary validation will be performed by comparing the results with on-sky measured PSF figures of merit obtained on specific targets using the SCAO systems SOUL (FLAO upgrade) feeding the camera LUCI at LBT and SAXO, the extreme SCAO system feeding the high resolution SPHERE instrument at VLT. This study will pave the way to the implementation of an operational forecasts of such a figure of merits on the base of existing operational forecast system of the atmosphere (turbulence and atmospheric parameters). In this contribution we focus our attention on the forecast of the PSF on-axis., Comment: 9 pages, 3 figures, 2 tables

A promising approach to study condensed-matter systems is to simulate them on an engineered quantum platform. However, achieving the accuracy needed to outperform classical methods has been an outstanding challenge. Here, using eighteen superconducting qubits, we provide an experimental blueprint for an accurate condensed-matter simulator and demonstrate how to probe fundamental electronic properties. We benchmark the underlying method by reconstructing the single-particle band-structure of a one-dimensional wire. We demonstrate nearly complete mitigation of decoherence and readout errors and arrive at an accuracy in measuring energy eigenvalues of this wire with an error of ~0.01 rad, whereas typical energy scales are of order 1 rad. Insight into this unprecedented algorithm fidelity is gained by highlighting robust properties of a Fourier transform, including the ability to resolve eigenenergies with a statistical uncertainty of 1e-4 rad. Furthermore, we synthesize magnetic flux and disordered local potentials, two key tenets of a condensed-matter system. When sweeping the magnetic flux, we observe avoided level crossings in the spectrum, a detailed fingerprint of the spatial distribution of local disorder. Combining these methods, we reconstruct electronic properties of the eigenstates where we observe persistent currents and a strong suppression of conductance with added disorder. Our work describes an accurate method for quantum simulation and paves the way to study novel quantum materials with superconducting qubits.

We study the multi-phase feedback processes in the central ~3 kpc of the barred Sy 2 galaxy NGC 5643. We use observations of the cold molecular gas (ALMA CO(2-1)) and ionized gas (MUSE IFU). We study different regions along the outflow zone which extends out to ~2.3 kpc in the same direction (east-west) as the radio jet, as well as nuclear/circumnuclear regions in the host galaxy disk. The deprojected outflowing velocities of the cold molecular gas (median Vcentral~189 km s^-1) are generally lower than those of the outflowing ionized gas, which reach deprojected velocities of up to 750 km s^-1 close to the AGN, and their spatial profiles follow those of the ionized phase. This suggests that the outflowing molecular gas in the galaxy disk is being entrained by the AGN wind. We derive molecular and ionized outflow masses of ~5.2x10^7 Msun and 8.5x10^4 Msun and molecular and ionized outflow mass rates of ~51 Msun yr^-1 and 0.14 Msun yr^-1. Therefore, the molecular phase dominates the outflow mass and outflow mass rate, while the outflow kinetic power and momentum are similar in both phases. However, the wind momentum load for the molecular and ionized outflow phases are ~27-5 and <1, which suggests that the molecular phase is not momentum conserving while the ionized one most certainly is. The molecular gas content (~1.5x10^7 Msun) of the eastern spiral arm is approximately 50-70% of the content of the western one. We interpret this as destruction/clearing of the molecular gas produced by the AGN wind impacting in the host galaxy. The increase of the molecular phase momentum implies that part of the kinetic energy from the AGN wind is transmitted to the molecular outflow. This suggest that in Sy-like AGN such as NGC 5643, the radiative/quasar and the kinetic/radio AGN feedback modes coexist and may shape the host galaxies even at kpc-scales via both positive and (mild) negative feedback., Comment: 26 pages, 21 figures, 5 tables. Astronomy and Astrophysics, Accepted 2020 September 11, in press

This work investigates the performance of the on-demand machine learning (ODML) algorithm introduced in Leal et al. (2020) when applied to different reactive transport problems in heterogeneous porous media. ODML was devised to accelerate the computationally expensive geochemical reaction calculations in reactive transport simulations. We demonstrate that the ODML algorithm speeds up these calculations by one to three orders of magnitude. Such acceleration, in turn, significantly accelerates the entire reactive transport simulation. The numerical experiments are performed by implementing the coupling of two open-source software packages: Reaktoro (Leal, 2015) and Firedrake (Rathgeber et al., 2016)., Comment: 37 pages, 21 figures, 6 tables

One distinctive feature of low-luminosity active galactic nuclei (LLAGN) is the relatively weak reflection features they may display in the X-ray spectrum, which can result from the disappearance of the torus with decreasing accretion rates. Some material, however, must surround the active nucleus, i.e., the accretion flow itself and, possibly, a flattened-out or thinned torus. In this work, we study whether reflection is indeed absent or undetectable due to its intrinsically weak features together with the low statistics inherent to LLAGN. Here we focus on NGC 3718 ($L/L_{\rm Edd}\sim10^{-5}$) combining observations from XMM--Newton and the deepest to date NuSTAR (0.5--79 keV) spectrum of a LLAGN, to constrain potential reflectors, and analyze how the fitted coronal parameters depend on the reflection model. We test models representing both an accretion disc (Relxill) and a torus-like (MYTorus and Borus) neutral reflector. From a statistical point of view, reflection is not required, but its inclusion allows to place strong constraints on the geometry and physical features of the surroundings: both neutral reflectors (torus) tested should be Compton thin ($N_H<10^{23.2}$cm$^{-2}$) and preferentially cover a large fraction of the sky. If the reflected light instead arises from an ionized reflector, a highly ionized case is preferred. These models produce an intrinsic power-law spectral index in the range [1.81--1.87], where the torus models result in steeper slopes. The cut-off energy of the power-law emission also changes with the inclusion of reflection models, resulting in constrained values for the disc reflectors and unconstrained values for torus reflectors., Comment: 15 pages and 16 Figures

We investigate the relation between the detection of the $11.3\,\mu$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO(2-1) transition. Galaxies with a nuclear detection of the $11.3\,\mu$m PAH feature contain more cold molecular gas (median $1.6\times 10^7\,M_\odot$) and have higher column densities ($N({\rm H}_2) = 2 \times 10^{23}\,{\rm cm}^{-2}$) over the regions sampled by the mid-IR slits than those without a detection. This suggests that molecular gas plays a role in shielding the PAH molecules in the harsh environments of Seyfert nuclei. Choosing the PAH molecule naphthalene as an illustration, we compute its half-life in the nuclear regions of our sample when exposed to 2.5keV hard X-ray photons. We estimate shorter half-lives for naphthalene in nuclei without a $11.3\,\mu$m PAH detection than in those with a detection. The Spitzer/IRS PAH ratios on circumnuclear scales ($\sim$ 4" $\sim$ 0.25-1.3kpc) are in between model predictions for neutral and partly ionized PAHs. However, Seyfert galaxies in our sample with the highest nuclear H$_2$ column densities are not generally closer to the neutral PAH tracks. This is because in the majority of our sample galaxies, the CO(2-1) emission in the inner $\sim$ 4" is not centrally peaked and in some galaxies traces circumnuclear sites of strong star formation activity. Spatially resolved observations with the MIRI medium-resolution spectrograph (MRS) on the James Webb Space Telescope will be able to distinguish the effects of an active galactic nucleus (AGN) and star formation on the PAH emission in nearby AGN., Comment: Accepted for publication in A&A

The ultralow frequency analogues of sound waves in Earth's magnetosphere play a crucial role in space weather, however, the public is largely unaware of this risk to our everyday lives and technology. As a way of potentially reaching new audiences, SSFX made 8 years of satellite wave recordings audible to the human ear with the aim of using it to create art. Partnering with film industry professionals, the standard processes of international film festivals were adopted by the project in order to challenge independent filmmakers to incorporate these sounds into short films in creative ways. Seven films covering a wide array of topics/genres (despite coming from the same sounds) were selected for screening at a special film festival out of 22 submissions. The works have subsequently been shown at numerous established film festivals and screenings internationally. These events have attracted diverse non-science audiences resulting in several unanticipated impacts upon them, thereby demonstrating how working with the art world can open up dialogues with both artists and audiences who would not ordinarily engage with science., Comment: Geoscience Communication

Depleted oil fields in the Gulf of Suez (Egypt) can serve as geothermal reservoirs for power generation using a CO2-Plume Geothermal (CPG) system, while geologically sequestering CO2. This entails the injection of a substantial amount of CO2 into the highly permeable brine-saturated Nubian Sandstone. Numerical models of two-phase flow processes are indispensable for predicting the CO2-plume migration at a representative geological scale. Such models require reliable constitutive relationships, including relative permeability and capillary pressure curves. In this study, quasi-static pore-network modelling has been used to simulate the equilibrium positions of fluid-fluid interfaces, and thus determine the capillary pressure and relative permeability curves. Three-dimensional images with a voxel size of 0.65 micro m3 of a Nubian Sandstone rock sample have been obtained using Synchrotron Radiation X-ray Tomographic Microscopy. From the images, topological properties of pores/throats were constructed. Using a pore-network model, we performed a sequential primary drainage-main imbibition cycle of quasi-static invasion in order to quantify (1) the CO2 and brine relative permeability curves, (2) the effect of initial wetting-phase saturation (i.e. the saturation at the point of reversal from drainage to imbibition) on the residual-trapping potential, and (3) study the relative permeability-saturation hysteresis. The results improve our understanding of the potential magnitude of capillary trapping in Nubian Sandstone, essential for future field-scale simulations. Further, an initial basin-scale assessment of CO2 storage capacity, which incorporates capillary trapping, yields a range of 14-49 GtCO2 in Nubian Sandstone, Gulf of Suez Basin.

We present the highlights from our recent study of 22 local (z$<$0.025) type-1 LINERs from the Palomar Survey, on the basis of optical long-slit spectroscopic observations taken with TWIN/CAHA, ALFOSC/NOT and HST/STIS (Cazzoli et al. 2018). Our goals were threefold: (a) explore the AGN- nature of these LINERs by studying the broad (BLR-originated) H$\alpha$ component; (b) derive a reliable interpretation for the multiple narrow components of emission lines by studying their kinematics and ionisation mechanism (via standard BPTs); (c) probe the neutral gas in the nuclei of these LINERs for the first time. Hence, kinematics and fluxes of a set of emission lines, from H$\beta$ to [SII], and the NaD doublet in absorption have been modelled and measured, after the subtraction of the underlying light from the stellar component., Comment: Proceedings of the IAU Symposium 356, 'Nuclear activity in galaxies across cosmic time'. Based on Cazzoli et al. 2018, MNRAS, 480, 1106 available at https://academic.oup.com/mnras/article/480/1/1106/5050381

Precise photometric and astrometric measurements on astronomical images require an accurate knowledge of the Point Spread Function (PSF). When the PSF cannot be modelled directly from the image, PSF-reconstruction techniques become the only viable solution. So far, however, their performance on real observations has rarely been quantified. Aims. In this Letter, we test the performance of a novel hybrid technique, called PRIME, on Adaptive Optics-assisted SPHERE/ZIMPOL observations of the Galactic globular cluster NGC6121. Methods. PRIME couples PSF-reconstruction techniques, based on control-loop data and direct image fitting performed on the only bright point-like source available in the field of view of the ZIMPOL exposures, with the aim of building the PSF model. Results. By exploiting this model, the magnitudes and positions of the stars in the field can be measured with an unprecedented precision, which surpasses that obtained by more standard methods by at least a factor of four for on-axis stars and by up to a factor of two on fainter, off-axis stars. Conclusions. Our results demonstrate the power of PRIME in recovering precise magnitudes and positions when the information directly coming from astronomical images is limited to only a few point-like sources and, thus, paving the way for a proper analysis of future Extremely Large Telescope observations of sparse stellar fields or individual extragalactic objects., Comment: 5 pages, 4 figures, accepted for publication by A&A Letters