The standard Λ Cold Dark Matter (ΛCDM) cosmological model provides a good description of a wide range of astrophysical and cosmological data. However, there are a few big open questions that make the standard model look like an approximation to a more realistic scenario yet to be found. In this paper, we list a few important goals that need to be addressed in the next decade, taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant H₀, the σ₈–S₈ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the 5.0σ tension between the Planck CMB estimate of the Hubble constant H₀ and the SH0ES collaboration measurements. After showing the H0 evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the Planck CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density Ω_m, and the amplitude or rate of the growth of structure (σ₈,fσ₈). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the H₀–S₈ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals. Finally, we give an overview of upgraded experiments and next-generation space missions and facilities on Earth that will be of crucial importance to address all these open questions.

The dynamics of charge qubits coupled to a nanomechanical resonator under the influence of both a phonon bath in contact with the resonator and irreversible decay of the qubits is considered. The focus of our analysis is devoted to multi-particle entanglement and the effects arising from the coupling to the reservoir. Even in the presence of the reservoirs, the inherent entanglement is found to be rather robust. Due to this fact, together with control of system parameters, the system may, therefore, be especially suited for quantum information processing. Our findings also shed light on the evolution of open quantum many-body systems. For instance, due to intrinsic qubit-qubit couplings our model is related to a driven XY spin model.

We present the fourth Fermi Large Area Telescope catalog (4FGL) of gamma-ray sources. Based on the first eight years of science data from the Fermi Gamma-ray Space Telescope mission in the energy range from 50 MeV to 1 TeV, it is the deepest yet in this energy range. Relative to the 3FGL catalog, the 4FGL catalog has twice as much exposure as well as a number of analysis improvements, including an updated model for the Galactic diffuse gamma-ray emission, and two sets of light curves (one-year and two-month intervals). The 4FGL catalog includes 5064 sources above 4 sigma significance, for which we provide localization and spectral properties. Seventy-five sources are modeled explicitly as spatially extended, and overall, 358 sources are considered as identified based on angular extent, periodicity, or correlated variability observed at other wavelengths. For 1336 sources, we have not found plausible counterparts at other wavelengths. More than 3130 of the identified or associated sources are active galaxies of the blazar class, and 239 are pulsars.

Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma-rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four supernova remnants (SNRs), IC 443, W44, W49B, and W51C, with the Fermi Large Area Telescope. This detection provided direct evidence that cosmic-ray protons are (re-)accelerated in SNRs. Here, we present a comprehensive search for low-energy spectral breaks among 311 4FGL catalog sources located within 5° from the Galactic plane. Using 8 yr of data from the Fermi Large Area Telescope between 50 MeV and 1 GeV, we find and present the spectral characteristics of 56 sources with a spectral break confirmed by a thorough study of systematic uncertainty. Our population of sources includes 13 SNRs for which the proton–proton interaction is enhanced by the dense target material; the high-mass gamma-ray binary LS I+61 303; the colliding wind binary η Carinae; and the Cygnus star-forming region. This analysis better constrains the origin of the gamma-ray emission and enlarges our view to potential new cosmic-ray acceleration sites.

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi Large Area Telescope (LAT) catalog of γ-ray sources. Based on the first 12 years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss the degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources.

We show theoretically that excitons can form from spatially separated one-dimensional ground state populations of electrons and holes, and that the resulting excitons can form a quasicondensate. We describe a mean-field Bardeen-Cooper-Schrieffer theory in the low carrier density regime and then focus on the core-shell nanowire giving estimates of the size of the excitonic gap for InAs/GaSb wires and as a function of all the experimentally relevant parameters. We find that optimal conditions for pairing include small overlap of the electron and hole bands, large effective mass of the carriers, and low dielectric constant of the surrounding media. Therefore, one-dimensional systems provide an attractive platform for the experimental detection of excitonic quasicondensation in zero magnetic field.

Previously, graphene nanoribbons set in lateral heterostructures with hexagonal boron nitride were predicted to support topologically protected states at low energy. We investigate how robust the transport properties of these states are against lattice disorder. We find that forms of disorder that do not couple the two valleys of the zigzag graphene nanoribbon do not impact the transport properties at low bias, indicating that these lateral heterostructures are very promising candidates for chip-scale conducting interconnects. Forms of disorder that do couple the two valleys, such as vacancies in the graphene ribbon, or substantial inclusions of armchair edges at the graphene-hexagonal boron nitride interface will negatively affect the transport. However, these forms of disorder are not commonly seen in current experiments.

We highlight the fact that two-dimensional (2D) materials with Dirac-like low energy band structures and spin-orbit coupling (SOC) will produce linearly dispersing topologically protected Jackiw-Rebbi modes at interfaces where the Dirac mass changes sign. These modes may support persistent spin or valley currents parallel to the interface, and the exact arrangement of such topologically protected currents depends crucially on the details of the SOC in the material. As examples, we discuss buckled 2D hexagonal lattices such as silicene or germanene, and transition metal dichalcogenides such as MoS2.

We model optical absorption of monolayer and bilayer graphene on hexagonal boron nitride for the case of closely aligned crystal lattices. We show that perturbations with different spatial symmetry can lead to similar absorption spectra. We suggest that a study of the absorption spectra as a function of the doping for an almost completely full first miniband is necessary to extract meaningful information about the moire characteristics from optical absorption measurements and to distinguish between various theoretical proposals for the physically realistic interaction. Also, for bilayer graphene, the ability to compare spectra for the opposite signs of electric-field-induced interlayer asymmetry might provide additional information about the moire parameters.

When the accretion disc around a weakly magnetised neutron star (NS) meets the stellar surface, it should brake down to match the rotation of the NS, forming a boundary layer. As the mechanisms potentially responsible for this braking are apparently inefficient, it is reasonable to consider this layer as a spreading layer (SL) with negligible radial extent and structure. We perform hydrodynamical 2D spectral simulations of an SL, considering the disc as a source of matter and angular momentum. Interaction of new, rapidly rotating matter with the pre-existing, relatively slow material co-rotating with the star leads to instabilities capable of transferring angular momentum and creating variability on dynamical timescales. For small accretion rates, we find that the SL is unstable for heating instability that disrupts the initial latitudinal symmetry and produces large deviations between the two hemispheres. This instability also results in breaking of the axial symmetry as coherent flow structures are formed and escape from the SL intermittently. At enhanced accretion rates, the SL is prone to shearing instability and acts as a source of oblique waves that propagate towards the poles, leading to patterns that again break the axial symmetry. We compute artificial light curves of an SL viewed at different inclination angles. Most of the simulated light curves show oscillations at frequencies close to 1 kHz. We interpret these oscillations as inertial modes excited by shear instabilities near the boundary of the SL. Their frequencies, dependence on flux, and amplitude variations can explain the high-frequency pair quasi-periodic oscillations observed in many low-mass X-ray binaries.

The GeV-range spectra of blazars are shaped not only by non-thermal emission processes internal to the relativistic jet but also by external pair-production absorption on the thermal emission of the accretion disc and the broad-line region (BLR). For the first time, we compute here the pair-production opacities in the GeV range produced by a realistic BLR accounting for the radial stratification and radiation anisotropy. Using photoionization modelling with the CLOUDY code, we calculate a series of BLR models of different sizes, geometries, cloud densities, column densities and metallicities. The strongest emission features in the model BLR are Ly alpha and He II Ly alpha. Contribution of recombination continua is smaller, especially for hydrogen, because Ly continuum is efficiently trapped inside the large optical depth BLR clouds and converted to Lyman emission lines and higher order recombination continua. The largest effects on the gamma-ray opacity are produced by the BLR geometry and localization of the gamma-ray source. We show that when the gamma-ray source moves further from the central source, all the absorption details move to higher energies and the overall level of absorption drops because of decreasing incidence angles between the gamma-rays and BLR photons. The observed positions of the spectral breaks can be used to measure the geometry and the location of the gamma-ray emitting region relative to the BLR. Strong dependence on geometry means that the soft photons dominating the pair-production opacity may be actually produced by a different population of BLR clouds than the bulk of the observed broad line emission.

Kilohertz-scale quasi-periodic oscillations (kHz QPOs) are a distinct feature of the variability of neutron star low-mass X-ray binaries. Among all the variability modes, they are especially interesting as a probe for the innermost parts of the accretion flow, including the accretion boundary layer (BL) on the surface of the neutron star. All the existing models of kHz QPOs explain only part of their rich phenomenology. Here, we show that some of their properties can be explained by a very simple model of the BL that is spun up by accreting rapidly rotating matter from the disk and spun down by the interaction with the neutron star. In particular, if the characteristic time scales for the mass and the angular momentum transfer from the BL to the star are of the same order of magnitude, our model naturally reproduces the so-called parallel tracks effect, where the QPO frequency is correlated with luminosity at time scales of hours but becomes uncorrelated at time scales of days. The closeness of the two time scales responsible for mass and angular momentum exchange between the BL and the star is an expected outcome of the radial structure of the BL.

We present a search for spatial extension in high-latitude (vertical bar b vertical bar > 5 degrees) sources in recent Fermi point source catalogs. The result is the Fermi High-Latitude Extended Sources Catalog, which provides source extensions (or upper limits thereof) and likelihood profiles for a suite of tested source morphologies. We find 24. extended sources, 19 of which were not previously characterized as extended. These include sources that are potentially associated with supernova remnants and star-forming regions. We also found extended.-ray emission in the vicinity of the Cen. A radio lobes and-at GeV energies for the first time-spatially coincident with the radio emission of the SNR CTA 1, as well as from the Crab Nebula. We also searched for halos around active galactic nuclei, which are predicted from electromagnetic cascades induced by the e(+)e(-) pairs that are deflected in intergalactic magnetic fields. These pairs are produced when gamma-rays interact with background radiation fields. We do not find evidence for extension in individual sources or in stacked source samples. This enables us to place limits on the flux of the extended source components, which are then used to constrain the intergalactic magnetic field to be stronger than 3 x 10(-16) G for a coherence length lambda greater than or similar to 10 kpc, even when conservative assumptions on the source duty cycle are made. This improves previous limits by several orders of magnitude.

The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi-large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission [unresolved gamma-ray background (UGRB)] below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluctuation field. This Letter presents a measurement of the UGRB autocorrelation angular power spectrum based on eight years of Fermi-LAT Pass 8 data products. The analysis is designed to be robust against contamination from resolved sources and noise systematics. The sensitivity to subthreshold sources is greatly enhanced with respect to previous measurements. We find evidence (with similar to 3.7 sigma significance) that the scenario in which two classes of sources contribute to the UGRB signal is favored over a single class. A double power law with exponential cutoff can explain the anisotropy energy spectrum well, with photon indices of the two populations being 2.55 +/- 0.23 and 1.86 +/- 0.15.

Context. Gravitational wave (GW) astronomy has rapidly reached maturity, becoming a fundamental observing window for modern astrophysics. The coalescences of a few tens of black hole (BH) binaries have been detected, while the number of events possibly including a neutron star (NS) is still limited to a few. On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. A preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS.

Aims. In this paper, we present our extensive search campaign aimed at uncovering the potential optical and near infrared electromagnetic counterpart of S190814bv. We found no convincing electromagnetic counterpart in our data. We therefore use our non-detection to place limits on the properties of the putative outflows that could have been produced by the binary during and after the merger.

Methods. Thanks to the three-detector observation of S190814bv, and given the characteristics of the signal, the LIGO and Virgo Collaborations delivered a relatively narrow localisation in low latency - a 50% (90%) credible area of 5 deg(2) (23 deg(2)) - despite the relatively large distance of 26752 Mpc. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical and near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS, and VINROUGE projects also carried out a search on this event. In this paper, we describe the combined observational campaign of these groups.

Results. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN), which was possibly generated by this NS-BH merger, and for the strategy of future searches. The typical depth of our wide-field observations, which cover most of the projected sky localisation probability (up to 99.8%, depending on the night and filter considered), is r similar to 22 (resp. K similar to 21) in the optical (resp. near infrared). We reach deeper limits in a subset of our galaxy-targeted observations, which cover a total similar to 50% of the galaxy-mass-weighted localisation probability. Altogether, our observations allow us to exclude a KN with large ejecta mass M greater than or similar to 0.1 M-circle dot to a high (> 90%) confidence, and we can exclude much smaller masses in a sub-sample of our observations. This disfavours the tidal disruption of the neutron star during the merger.

Conclusions. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv, we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundred megaparsecs will be detected only by large facilities with both a high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.