Papers for Thursday, Mar 10 2022

Recent observations and simulations indicate substantial evolution in the properties of galaxies with time, wherein rotationally-supported and steady thin discs (like those frequently observed in the local universe) emerge from galaxies that are clumpy, irregular, and have bursty star formation rates (SFRs). To better understand the progenitors of local disc galaxies we carry out an analysis of three FIRE-2 simulated galaxies with a mass similar to the Milky Way at redshift z=0. We show that all three galaxies transition from bursty to steady SFRs at a redshift between z=0.5 and z=0.8, and that this transition coincides with a rapid (< ~1 Gyr) emergence of a rotationally-supported interstellar medium (ISM).In the late phase with steady SFR, the rotational energy comprises > ~90% of the total kinetic + thermal energy in the ISM, and is roughly half the gravitational energy. By contrast, during the early phase with bursty star formation, the ISM has a quasi-spheroidal morphology and its energy budget is dominated by quasi-isotropic flows including turbulence and coherent inflows/outflows. This result, that rotational support is subdominant at early times, challenges the common application of equilibrium disc models to the high-redshift progenitors of Milky Way-like galaxies. We further find that the formation of a rotation-supported ISM coincides with the formation of a thermal energy-supported inner circumgalactic medium (CGM). Before this transition, the inner CGM is also supported by turbulence and coherent flows, indicating that at early times there is no clear boundary between the ISM and inner CGM.

We characterize the average X-ray and radio properties of quiescent galaxies (QGs) with $\log{(M_\star/M_\odot)}>10$ at $0<z<5$. QGs are photometrically selected from the latest COSMOS2020 catalog. We conduct the stacking analysis of X-ray images of the Chandra COSMOS Legacy Survey for individually undetected QGs. Thanks to the large sample and deep images, the stacked X-ray signal is significantly detected up to $z\sim5$. The average X-ray luminosity can not be explained by the X-ray luminosity of X-ray binaries, suggesting that the low-luminosity active galactic nuclei (AGNs) ubiquitously exist in QGs. Moreover, the X-ray AGN luminosity of QGs at $z>1.5$ is higher than that of star-forming galaxies (SFGs), derived in the same manner as QGs. The stacking analysis of the VLA-COSMOS images is conducted for the identical sample, and the radio signal for QGs is also detected up to $z\sim5$. We find that the radio AGN luminosity of QGs at $z>1.5$ is also higher than SFGs, which is in good agreement with the X-ray analysis. The enhanced AGN activity in QGs suggested by the individual analysis in the X-ray and radio wavelength supports its important role for quenching at high redshift. Their enhanced AGN activity is less obvious at $z<1.5$, which can be interpreted as an increasing role of others at lower redshifts, such as environmental quenching.

In order to explore how the early internal rotational properties of star clusters are affected by the external potential of their host galaxies, we have run a suite of $N$-body simulations following the early dynamical evolution and violent relaxation of rotating star clusters embedded in a tidal field. Our study focuses on models for which the cluster's rotation axis has a generic orientation relative to the torque of the tidal field. The interaction between the violent relaxation process, angular momentum of the cluster, and the external torque creates a complex kinematic structure within the cluster, most prominently a radial variation in the position of the rotation axis along both the polar and azimuthal directions. We also examine the cluster's velocity dispersion anisotropy and show that the projected anisotropy may be affected by the variation of the rotation axis directions within the cluster; the combination of projection effects and the complex kinematical features may result in the measurement of tangential anisotropy in the cluster's inner regions. We also characterize the structural properties of our clusters as a function of their initial rotation and virial ratio and find that clusters may develop a triaxial morphology and a radial variation of the minor axis not necessarily aligned with the rotation axis. Finally, we examine the long-term evolution of these complex kinematic features.

Progenitor models for the "luminous" subclass of Fast Blue Optical Transients (LFBOTs; prototype: AT2018cow) are challenged to simultaneously explain all of their observed properties: fast optical rise times < days; peak luminosities >1e44 erg/s; low yields <0.1 Msun of 56Ni; aspherical ejecta with a wide velocity range (<3000 km/s to > 0.1-0.5 c with increasing polar latitude); presence of hydrogen-depleted-but-not-free dense circumstellar material (CSM) on radial scales from ~1e14 cm to ~3e16 cm; embedded variable source of non-thermal X-ray/gamma-rays, suggestive of a compact object. We show that all of these properties are consistent with the tidal disruption and hyper-accretion of a Wolf-Rayet (WR) star by a black hole (BH) or neutron star (NS) binary companion. In contrast with related previous models, the merger occurs with a long delay following the common envelope (CE) event responsible for birthing the binary, as a result of gradual angular momentum loss to a relic circumbinary disk. Disk-wind outflows from the merger-generated accretion flow generate the 56Ni-poor aspherical ejecta with the requisite velocity range. The optical light curve is powered primarily by reprocessing X-rays from the inner accretion flow/jet, though CSM shock interaction also contributes. Primary CSM sources include mass-loss from WR star (e.g., from the L2 point) during the earliest stages of the merger (~<1e14 cm) and the relic CE disk and its photoevaporation-driven wind (>~ 1e16 cm). Longer delayed mergers may instead give rise to supernovae Type Ibn/Icn (depending on the WR evolutionary state), potentially connecting these transient classes with LFBOTs.

We observed the high-mass protostellar core G335.579-0.272 ALMA1 at ${\sim}200$ au (0.05") resolution with the Atacama Large Millimeter/submillimeter Array (ALMA) at 226 GHz (with a mass sensitivity of $5\sigma=0.2$ M$_\odot$ at 10 K). We discovered that at least a binary system is forming inside this region, with an additional nearby bow-like structure (${\lesssim}1000$ au) that could add an additional member to the stellar system. These three sources are located at the center of the gravitational potential well of the ALMA1 region and the larger MM1 cluster. The emission from CH$_3$OH (and many other tracers) is extended ($>1000$ au), revealing a common envelope toward the binary system. We use CH$_2$CHCN line emission to estimate an inclination angle of the rotation axis of $26^\circ$ with respect to the line of sight based on geometric assumptions and derive a kinematic mass of the primary source (protostar+disk) of 3.0 M$_\odot$ within a radius of 230 au. Using SiO emission, we find that the primary source drives the large scale outflow revealed by previous observations. Precession of the binary system likely produces a change in orientation between the outflow at small scales observed here and large scales observed in previous works. The bow structure may have originated by entrainment of matter into the envelope due to widening or precession of the outflow, or, alternatively, an accretion streamer dominated by the gravity of the central sources. An additional third source, forming due to instabilities in the streamer, cannot be ruled out as a temperature gradient is needed to produce the observed absorption spectra.

We analyse observational signatures of magnetic fields for simulations of a Milky-Way like disc with supernova-driven interstellar turbulence and self-consistent chemical processes. In particular, we post-process two simulations data sets of the SILCC Project for two initial amplitudes of the magnetic field, $B_0 = $ 3 and 6 $\mu$G, to study the evolution of Faraday rotation measures (RM) and synchrotron luminosity. For calculating the RM, three different models of the electron density $n_e$ are considered. A constant electron density, and two estimations based on the density of ionized species and the fraction of the total gas, respectively. Our results show that the RM profiles are extremely sensitive to the $n_e$ models, which assesses the importance of accurate electron distribution observations/estimations for the magnetic fields to be probed using Faraday RMs. As a second observable of the magnetic field, we estimate the synchrotron luminosity in the simulations using a semi-analytical cosmic ray model. We find that the synchrotron luminosity decreases over time, which is connected to the decay of magnetic energy in the simulations. The ratios between the magnetic, the cosmic ray, and the thermal energy density indicate that the assumption of equipartition does not hold for most regions of the ISM. In particular, for the ratio of the cosmic ray to the magnetic energy the assumption of equipatition could lead to a wrong interpretation of the observed synchrotron emission.

Magnetic reconnection in relativistic plasmas is well established as a fast and efficient particle accelerator, capable of explaining the most dramatic astrophysical flares. With particle-in-cell simulations, we demonstrate the importance of non-ideal fields for the early stages ("injection") of particle acceleration. Most of the particles ending up with high energies (near or above the mean magnetic energy per particle) must have passed through non-ideal regions where the assumptions of ideal magnetohydrodynamics are broken (i.e., regions with $E>B$ or nonzero $E_\parallel={\bf E}\cdot {\bf B}/B$), whereas particles that do not experience non-ideal fields end up with Lorentz factors of order unity. Thus, injection by non-ideal fields is a necessary prerequisite for further acceleration. Our results have important implications for the origin of nonthermal particles in high-energy astrophysical sources.

Neutrino-cooled accretion flow around a black hole, produced by a compact binary merger, is a promising scenario for jet formation and magnetic-driven winds to explain short duration gamma ray bursts (GRBs) central engine and kilonovae based on GW170817 gravitational wave observation. Magnetorotational instability (MRI) turbulence and Blandford-Znajek (BZ) mechanism are expected to play key roles in the thermal equilibrium of the disk (balancing neutrino cooling) and in driving accretion and creating jets. Using the open-source GRMHD HARM-COOL code, we study the magnetically-driven evolution of an accretion disk with realistic equation of state in the fixed curved space-time background. We identify the effects of the neutrino cooling and the magnetic field, paying particular attention to the dynamical, thermal and composition evolution of the disk and outflows.

The chromospheric H-alpha spectral line is a strong line in the spectrum of the Sun and other stars. In the stellar regime, this spectral line is already used as a powerful tracer of stellar activity. For the Sun, other tracers (i.e, CaII K) are typically used to monitor solar activity. We used observations of full-disk H-alpha filtergrams of the Chromospheric Telescope (ChroTel) to extract the imaging H-alpha excess and deficit, which are related to bright features (plage regions) and dark absorption features (filaments and sunspots), respectively. The aim of this study is to introduce the imaging H-alpha excess and deficit as tracers of solar activity and compare them to other established indicators: the relative sunspot number, the F10.7cm radio flux, and the MgII index. The H-alpha excess and deficit follow the behavior of the solar activity over the course of the cycle, whereby the peak of the H-alpha deficit is shortly after the solar maximum. The H-alpha excess is closely correlated to the chromospheric MgII index. The highest correlation of the H-alpha deficit is found with the F10.7cm radio flux. The H-alpha deficit reflects the cyclic behavior of polar crown filaments and their disappearance shortly before the solar maximum. We investigated the mean intensity distribution for H-alpha excess regions for solar minimum and maximum, whereby the shape of the distributions is very similar, but with different amplitudes. Furthermore, we investigate whether the area coverage fraction or the changing H-alpha excess in the active regions dominates temporal variability in solar H-alpha observations. The area coverage fraction and the H-alpha excess are strongly correlated, whereas the weak correlation between the area coverage fraction and mean intensity leaves us pessimistic that the degeneracy between these two quantities can be broken for the modeling of unresolved stellar surfaces.

The radio emission from very inclined air showers (with zenith angles $\theta \gtrsim 65^\circ$) illuminates large footprints from a few up to hundreds of square kilometers. Inclined air showers can be detected with kilometer-spaced radio-antenna arrays, suitable to measure cosmic rays up to the highest energies ($10^{18.4}\,\text{eV} \lesssim E \sim 10^{20}\,\text{eV}$). Radio antennas are sensitive to the electromagnetic particle cascades of the showers, which are absorbed in the atmosphere before reaching the ground. Hence, combining radio antennas with particle detectors, which, in turn, measure the muons reaching the ground, offers the unique potential to discriminate between showers initiated by light and heavy nuclei and study hadronic interactions. A precise understanding and description of the radio-emission footprint of inclined air showers is indispensable for any analyses of radio data. Such a model has to master the description of complex asymmetry patterns imprinted in the radio-emission. We propose a model that describes the radio-emission from inclined air showers by explicitly describing the dominant, rotationally symmetric geomagnetic emission in combination with additional effects which disturb the symmetric emission pattern. We use a comprehensive set of CoREAS simulations to exploit correlations between the radio emission and air shower properties. We find a description which is based on two observables only: the geometrical distance $d_\mathrm{max}$ between the shower core at the ground and the shower maximum and the geomagnetic radiation energy $E_\mathrm{geo}$. Hence, this model can be used in a robust and efficient reconstruction algorithms. We demonstrate that with this model, the electromagnetic shower energy can be reconstructed by sparse antenna arrays with an intrinsic resolution of only 5\% and a negligible bias.

Filaments are an ubiquitous feature of molecular clouds, and appear to play a critical role in assembling the material to form stars. The dominant filaments are observed to have a rather narrow range of widths around $\sim 0.1$ pc, and to be preferentially aligned perpendicularly to the direction of the local magnetic field. We have previously argued that the observed filament widths can be explained if filaments are formed by converging, mildly supersonic flows, resulting from large-scale turbulent motions in the parent molecular cloud. Here we demonstrate that the introduction of a magnetic field perpendicular to the filament long axis does not greatly alter this conclusion, as long as the mass-to-flux ratio is supercritical. The distribution of widths for supercritical magnetised filaments formed via this mechanism is peaked at slightly higher values, and is slightly broader, than for non-magnetised filaments, but still reproduces the basic properties of the width distributions derived from far-infrared observations of molecular clouds. In contrast, subcritical filaments have width distributions with a fundamentally different shape, and typically have much larger widths than those observed. Both subcritical and supercritical filaments are consistent with the observed lack of correlation between filament widths and filament surface densities.

In coronal loop modeling, it is commonly assumed that the loops are semi-circular with a uniform cross-sectional area. However, observed loops are rarely semi-circular, and extrapolations of the magnetic field show that the field strength decreases with height, implying that the cross-sectional area should expand with height. In this work, we examine these two assumptions directly to understand how they affect the hydrodynamic and radiative response to strong, impulsive heating events. Both the magnitude and rate of area expansion impact the dynamics directly, and we show that an expanding cross-section significantly lengthens the time for a loop to cool and drain, increases upflow durations, and suppresses sound waves. An increase in the eccentricity of loops, on the other hand, only increases the draining timescale, and is a minor effect in general. Spectral line intensities are also strongly impacted by the variation in the cross-sectional area since they depend on both the volume of the emitting region as well as the density and ionization state. With a larger expansion, the density is reduced, so the lines at all heights are relatively reduced in intensity and, because of the increase of cooling times, the hottest lines remain bright for significantly longer. Future modeling work needs to include area expansion for an accurate picture of the hydrodynamics, and future observations are needed to provide tighter constraints on the magnitude, rate, and location of the expansion or lack thereof.

Using photometry and proper motions from Gaia Early Data Release 3, we detect a 45 degree-long trailing stellar debris stream associated with the old, metal-poor globular cluster NGC 7089. With a width on the order of 100 pc, the extended stream appears to be as dynamically cold as the coldest known streams found to date. There is some evidence for an extended leading tail extending between 28 and 37 degrees from the cluster, though the greater distance of this tail, combined with proper motions that are virtually indistinguishable from those of foreground stars, make the detection much less certain. The proper motion profile and the path on the sky of the trailing tail are not well matched using a simple Galactic potential comprised purely of a disk, bulge, and spherical halo. However, the addition of a moving, massive (M = 1.88 x 10^(11) solar masses) Large Magellanic Cloud brings the model predictions into much better agreement with the observables. We provide tables of the most highly ranked candidate stream stars for follow-up by ongoing and future spectroscopic surveys.

Recent advances in sub-millimeter observations of young circumstellar nebulae have opened an unprecedented window into the structure of protoplanetary disks, which has revealed the surprising ubiquity of broken and misaligned disks. In this work, we demonstrate that such disks are capable of torquing the spin axis of their host star, representing a hitherto unexplored pathway by which stellar obliquities may be generated. The basis of this mechanism is a crossing of the stellar spin precession and inner disk regression frequencies, resulting in adiabatic excitation of the stellar obliquity. We derive analytical expressions for the characteristic frequencies of the inner disk and star as a function of the disk gap boundaries, and place an approximate limit on the disk architectures for which frequency crossing and resulting obliquity excitation are expected, thereby illustrating the efficacy of this model. Cumulatively, our results support the emerging concensus that significant spin-orbit misalignments are an expected outcome of planet formation.

On the one hand, several types of global-scale inertial modes of oscillation have been observed on the Sun. They include the equatorial Rossby modes, critical-latitude modes, and high-latitude modes. On the other hand, the columnar convective modes (predicted by simulations; also known as banana cells or thermal Rossby waves) remain elusive. We aim to investigate the influence of turbulent diffusivities, non-adiabatic stratification, differential rotation, and a latitudinal entropy gradient on the linear global modes of the rotating solar convection zone. We solve numerically for the eigenmodes of a rotating compressible fluid inside a spherical shell. We identify modes in the inertial frequency range including the columnar convective modes, as well as modes of mixed character. The corresponding mode dispersion relations and eigenfunctions are computed for azimuthal orders $m \leq 16$. The three main results are as follows. Firstly, we find that, for $m \gtrsim 5$, the radial dependence of the equatorial Rossby modes with no radial node ($n=0$) is radically changed from the traditional expectation ($r^m$) for turbulent diffusivities $\gtrsim 10^{12}$ cm$^2$ s$^{-1}$. Secondly, we find mixed modes, i.e. modes that share properties of the equatorial Rossby modes with one radial node ($n=1$) and the columnar convective modes. Thirdly, we show that the $m=1$ high-latitude mode in the model is consistent with the solar observations when the latitudinal entropy gradient corresponding to a thermal wind balance is included (baroclinally unstable mode). To our knowledge, this work is the first realistic eigenvalue calculation of the global modes of the rotating solar convection zone. This calculation reveals a rich spectrum of modes in the inertial frequency range, which can be directly compared to the observations. In turn, the observed modes can inform us about the solar convection zone.

The hot Neptune desert, a distinct lack of highly irradiated planets in the size range of Neptune, remains one of the most intriguing results of exoplanet population studies. A deeper understanding of the atmosphere of exoplanets sitting at the edge or even within the Neptune desert will allow us to better understand if planetary formation or evolution processes are at the origin of the desert. A detection of sodium in WASP-166b was presented previously with tentative line broadening at the 3.4 sigma with the HARPS spectrograph. We update this result with two transits observed with the ESPRESSO spectrograph, confirming the detection in each night and the broadened character of the line. This result marks the first confirmed resolved sodium detection within the Neptune desert. In this work, we additionally highlight the importance of treating low-SNR spectral regions, particularly where absorption lines of stellar sodium and planetary sodium overlap at mid-transit - an important caveat for future observations of the system.

Abundances of fluorine ($^{19}$F), as well as isotopic ratios of $^{16}$O/$^{17}$O, are derived in a sample of luminous young ($\sim$10$^{7}$--10$^{8}$ yrs) red giants in the Galactic center (with galactocentric distances ranging from 0.6--30 pc), using high-resolution infrared spectra and vibration-rotation lines of H$^{19}$F near $\lambda$2.3$\mu$m. Five of the six red giants are members of the Nuclear star cluster that orbits the central supermassive black hole. Previous investigations of the chemical evolution of $^{19}$F in Galactic thin and thick disk stars have revealed that the nucleosynthetic origins of $^{19}$F may be rather complex, resulting from two, or more, astrophysical sites; fluorine abundances behave as a primary element with respect to Fe abundances for thick disk stars and as a secondary element in thin disk stars. The Galactic center red giants analyzed fall within the thin disk relation of F with Fe, having near-solar, to slightly larger, abundances of Fe ($<$[Fe/H]$>$=+0.08$\pm$0.04), with a slight enhancement of the F/Fe abundance ratio ($<$[F/Fe]$>$=+0.28$\pm$0.17). In terms of their F and Fe abundances, the Galactic center stars follow the thin disk population, which requires an efficient source of $^{19}$F that could be the winds from core-He burning Wolf Rayet stars, or thermally-pulsing AGB stars, or a combination of both. The observed increase of [F/Fe] with increasing [Fe/H] found in thin disk and Galactic center stars is not predicted by any published chemical evolution models that are discussed, thus a quantitative understanding of yields from the various possible sources of $^{19}$F remains unknown.

Contemporary reverberation mapping campaigns are employing wide-area photometric data and high-multiplex spectroscopy to efficiently monitor hundreds of active galactic nuclei (AGN). However, the interaction of the window function(s) imposed by the observation cadence with the reverberation lag and AGN variability time scales (intrinsic to each source over a range of luminosities) impact our ability to recover these fundamental physical properties. Time dilation effects due to the sample source redshift distribution introduces added complexity. We present comprehensive analysis of the implications of observational cadence, seasonal gaps and campaign baseline duration (i.e., the survey window function) for reverberation lag recovery. We find the presence of a significant seasonal gap dominates the efficacy of any given campaign strategy for lag recovery across the parameter space, particularly for those sources with lags above $\sim$100 days in the observed-frame. Using the OzDES survey as a baseline, we consider the implications of this analysis for the 4MOST/TiDES campaign providing concurrent follow-up of the LSST deep-drilling fields, as well as upcoming programs. We conclude the success of such surveys will be critically limited by the seasonal visibility of some potential field choices, but show significant improvement from extending the baseline. Optimising the sample selection to fit the window function will improve survey efficacy.

The standard alpha-disc model predicts an anti-correlation between the density of the inner accretion disc and the black hole mass times square of the accretion rate, as seen in higher mass ($M_{\rm BH}>10^{6} M_{\odot}$) active galactic nuclei (AGNs). In this work, we test the predictions of the alpha-disc model and study the properties of the inner accretion flow for the low-mass end ($M_{\rm BH}\approx 10^{5-6}M_{\odot}$) of AGNs. We utilize a new high-density disc reflection model where the density parameter varies from $n_{\rm e}=10^{15}$ to $10^{20}$ cm$^{-3}$ and apply it to the broadband X-ray (0.3-10 keV) spectra of the low-mass AGN sample. The sources span a wide range of Eddington fractions and are consistent with being sub-Eddington or near-Eddington. The X-ray spectra reveal a soft X-ray excess below $\sim 1.5$ keV which is well modeled by high-density reflection from an ionized accretion disc of density $n_{\rm e}\sim 10^{18}$ cm$^{-3}$ on average. The results suggest a radiation pressure-dominated disc with an average of 70% fraction of the disc power transferred to the corona, consistent with that observed in higher mass AGNs. We show that the disc density higher than $10^{15}$ cm$^{-3}$ can result from the radiation pressure compression when the disc surface does not hold a strong magnetic pressure gradient. We find tentative evidence for a drop in black hole spin at low-mass regimes.

A well-established constraint on the size of non-porous olivine grains or the porosity of aggregates consisting of small olivine grains from prominent narrow peaks in thermal infrared spectra characteristic of crystalline silicates is reexamined. To thoroughly investigate thermal infrared peaks, we make theoretical argument for the absorption and scattering of light by non-porous, non-spherical olivine particles, which is followed by numerical verification. Our study provides perfectly rational explanations of the physics behind the small-particle effect of emission peaks in the framework of classical electrodynamics and convincing evidence of small-particle's emission peaks in the literature. While resonant absorption excited by surface roughness on the order of submicrometer scales can be identified even for non-porous olivine particles with a radius of $10~{\rm \mu m}$, it makes only a negligible contribution to thermal infrared spectra of the particles. In contrast, the porosity of non-spherical particles has a significant impact on the strength and wavelength of the peaks, while the resonant absorption excited by an ensemble of small grains takes place at a wavelength different than one expects for surface roughness. We finally reaffirm that twin peaks of olivine in thermal infrared spectra of dust particles in astronomical environments are the intrinsic diagnostic characters of submicrometer-sized small grains and their aggregate particles in fluffy and porous configurations.

In astronomy, if we denote the dimension of data as $d$ and the number of samples as $n$, we often meet a case with $n \ll d$. Traditionally, such a situation is regarded as ill-posed, and there was no choice but to throw away most of the information in data dimension to let $d < n$. The data with $n \ll d$ is referred to as high-dimensional low sample size (HDLSS). {}To deal with HDLSS problems, a method called high-dimensional statistics has been developed rapidly in the last decade. In this work, we first introduce the high-dimensional statistical analysis to the astronomical community. We apply two representative methods in the high-dimensional statistical analysis methods, the noise-reduction principal component analysis (NRPCA) and regularized principal component analysis (RPCA), to a spectroscopic map of a nearby archetype starburst galaxy NGC 253 taken by the Atacama Large Millimeter/Submillimeter Array (ALMA). The ALMA map is a typical HDLSS dataset. First we analyzed the original data including the Doppler shift due to the systemic rotation. The high-dimensional PCA could describe the spatial structure of the rotation precisely. We then applied to the Doppler-shift corrected data to analyze more subtle spectral features. The NRPCA and RPCA could quantify the very complicated characteristics of the ALMA spectra. Particularly, we could extract the information of the global outflow from the center of NGC 253. This method can also be applied not only to spectroscopic survey data, but also any type of data with small sample size and large dimension.

Naphthalene (C10H8) is the simplest polycyclic aromatic hydrocarbon (PAH) and an important component in a series of astrochemical reactions involving hydrocarbons. Its molecular charge state affects the stability of its isomeric structures, which is specially relevant in ionised astrophysical environments. We thus perform an extensive computational search for low-energy molecular structures of neutral, singly, and multiply charged naphthalene and its isomers with charge states +q = 0-4 and investigate their geometric properties and bonding situations. We find that isomerisation reactions should be frequent for higher charged states and that open chains dominate their low-energy structures. We compute both the scaled-harmonic and anharmonic infrared spectra of selected low-energy species and provide the calculated scaling factors for the naphthalene neutral, cation, and dication global minima. All simulated spectra reproduce satisfactorily the experimental data and, thus, are adequate for aiding observations. Moreover, the potential presence of these species in the emission spectra of the circumnuclear regions of active galactic nuclei (AGNs), with high energetic X-ray photon fluxes, is explored using the experimental value of the naphthalene photodissociation cross-section, \sigma_{ph-d}, to determine its half-life, t_{1/2}, at a photon energy of 2.5 keV in a set of relevant sources. Finally, we show that the computed IR bands of the triply and quadruply charged species are able to reproduce some features of the selected AGN sources.

AstroSat is India's first Ultra-violet (UV) and X-ray astronomy observatory in space. The satellite was launched by the Indian Space Research Organisation on a Polar Satellite Launch Vehicle on 28 September 2015 from Sriharikota Range north of Chennai on the eastern coast of India. AstroSat carries five scientific instruments and one auxiliary instrument. Four of these consist of co-aligned telescopes and detectors mounted on a common deck of the satellite to observe stars and galaxies simultaneously in the near- and far-UV wavelengths and a broad range of X-ray energies (0.3 to 80 keV). The fifth instrument consists of three X-ray detectors and is mounted on a rotating platform on a side that is oriented 90 degrees with respect to the other instruments to scan the sky for X-ray transients. An auxiliary instrument monitors the charged particle environment in the path of the satellite.

We present the complete Bayesian statistical analysis of the HArps-n red Dwarf Exoplanet Survey (HADES), which monitored the radial velocities of a large sample of M dwarfs with HARPS-N at TNG, over the last 6 years. The targets were selected in a narrow range of spectral types from M0 to M3, $0.3$ M$_\odot < M_\star < 0.71$ M$_\odot$, in order to study the planetary population around a well-defined class of host stars. We take advantage of Bayesian statistics to derive an accurate estimate of the detectability function of the survey. Our analysis also includes the application of Gaussian Process approach to take into account stellar activity induced radial velocity variations, and improve the detection limits, around the most-observed and most-active targets. The Markov chain Monte Carlo and Gaussian process technique we apply in this analysis has proven very effective in the study of M-dwarf planetary systems, helping the detection of most of the HADES planets. From the detectability function we can calculate the occurrence rate of small mass planets around early-M dwarfs, either taking into account only the 11 already published HADES planets or adding also the 5 new planetary candidates discovered in this analysis, and compare them with the previous estimates of planet occurrence around M-dwarf or Solar-type stars: considering only the confirmed planets, we find the highest frequency for low-mass planets ($1$ M$_\oplus < m_p \sin i < 10$ M$_\oplus$) with periods $10$ d$ < P < 100$ d, $f_\text{occ} = 85^{+5}_{-19}\%$, while for short-period planets ($1$ d$ < P < 10$ d) we find a frequency of $f_\text{occ} = 10.3^{+8.4}_{-3.3}\%$, significantly lower than for later-M dwarfs. These results, and their comparison with other surveys focused on different stellar types, confirms the central role that stellar mass plays in the formation and evolution of planetary systems.

Context. Emission lines formed in the transition region and corona show dominantly redshifts and blueshifts, respectively. Aims. We investigate the Doppler shifts in a 3D radiation MHD model of the quiet Sun and compare these to observed properties. We concentrate on Si IV 1394 A originating in the transition region and examine the Doppler shifts of several other spectral lines at different formation temperatures. Methods. We construct a radiation MHD model extending from the upper convection zone to the lower corona using the MURaM code. In this quiet Sun model the magnetic field is self-consistently maintained by the action of a small-scale dynamo. We synthesize the profiles of several optically thin emission lines, formed at temperatures from the transition region into the corona. We investigate the spatial structure and coverage of red- and blueshifts and how this changes with line-formation temperature. Results. The model successfully reproduces the observed change of average net Doppler shifts from red- to blueshifted from the transition region into the corona. In particular, the model shows a clear imbalance of area coverage of red- vs. blueshifts in the transition region of ca. 80% to 20%. We determine that (at least) four processes generate the systematic Doppler shifts in our model, including pressure enhancement in the transition region, transition region brightenings unrelated to coronal emission, boundaries between cold and hot plasma, and siphon-type flows. Conclusions. We show that there is not a single process that is responsible for the observed net Doppler shifts in the transition region and corona. Because current 3D MHD models do not yet fully capture the evolution of spicules, one of the key ingredients of the chromosphere, most probably these have still to be added to the list of processes responsible for the persistent Doppler shifts.

We present an overview of PION, an open-source software project for solving radiation-magnetohydrodynamics equations on a nested grid, aimed at modelling asymmetric nebulae around massive stars. A new implementation of hybrid OpenMP/MPI parallel algorithms is briefly introduced, and improved scaling is demonstrated compared with the current release version. Three-dimensional simulations of an expanding nebula around a Wolf-Rayet star are then presented and analysed, similar to previous 2D simulations in the literature. The evolution of the emission measure of the gas and the X-ray surface brightness are calculated as a function of time, and some qualitative comparison with observations is made.

The standard method to generate dynamical models with a finite extent is to apply a truncation in binding energy to the distribution function. This approach has the disadvantages that one cannot choose the density to start with, that the important dynamical quantities cannot be calculated analytically, and that a fraction of the possible bound orbits are excluded a priori. We explore another route and start from a truncation in radius rather than a truncation in binding energy. We focus on the simplest truncated density profile, the uniform density sphere. We explore the most common inversion techniques to generate distribution functions for the uniform density sphere, corresponding to a large range of possible anisotropy profiles. We find that the uniform density sphere cannot be supported by the standard isotropic, constant anisotropy or Osipkov-Merritt models, as all these models are characterised by negative distribution functions. We generalise the Cuddeford inversion method to models with a tangential anisotropy and present a one-parameter family of dynamical models for the uniform density sphere. Each member of this family is characterised by an anisotropy profile that smoothly decreases from an arbitrary value $\beta_0\leqslant0$ at the centre to completely tangential at the outer radius. All models have a positive distribution function over the entire phase space, and a nonzero occupancy of all possible bound orbits. This shows that one can generate nontrivial self-consistent dynamical models based on preset density profile with a finite extent.

The paper presents the results of the analysis of the geometric characteristics of sunspots for the period of 19-24 cycles of activity. The shape of sunspots was studied on the basis of the method of normalization of images of sunspots to study the average profile of the spot. The deviation of the shape of sunspots from the axisymmetric configuration is investigated. It was found that the spots, as a rule, have an ellipsoid shape, and the major axis of the ellipse has a predominant inclination to the equator, opposite in the Northern and Southern hemispheres. The angle of inclination of the sunspot axis corresponds to the angle of inclination of the bipoles in the activity cycles. The relationship between the shape of sunspots in the current cycle and the amplitude of the next cycle of activity is found. The greater the elongation along the longitude of the current cycle of spots, the higher the next cycle of activity will be.

We present a simple but general argument that strongly limits the abundance of Primordial Black Holes (PBHs) (or other unknown population of compact objects) with masses similar to those determined by LIGO/Virgo from BH binary mergers. We show that quasar microlensing can be very sensitive to the mass of the lenses, and that it is able to distinguish between stars and BHs of high mass, when the finite size of the source is taken into account. A significant presence of massive BHs would produce frequent high flux magnifications (except for unrealistically large sources) which have been very rarely observed. On the contrary, a typical stellar population would induce flux magnifications consistent with the observations. This result excludes PBHs (or any type of compact object) in the mass range determined by LIGO/Virgo as the main dark matter constituents in the lens galaxies.

Assuming a population of Black Holes (BHs) with masses in the range inferred by LIGO/Virgo from BH mergers, we use quasar microlensing observations to estimate their abundances. We consider a mixed population of stars and BHs and the presence of a smooth dark matter component. We adopt reverberation mapping estimates of the quasar size. According to a Bayesian analysis of the measured microlensing magnifications, a population of BHs with masses $\sim$ 30$M_{\odot}$ constitutes less than 0.4 % of the total matter at 68 % confidence level (less than 0.9 % at 90 % confidence). We have explored the whole mass range of LIGO/Virgo BHs finding that this upper limit ranges from 0.5 % to 0.4 % at 68 % C.L. (from 1.1 % to 0.9 % at 90 % C.L.) when the BHs mass change from 10 to 60$M_{\odot}$. We estimate a 16 % contribution from the stars, in agreement with previous studies based on a single mass population that do not consider explicitly the presence of BHs. These results are consistent with the estimates of BH abundances from the statistics of LIGO/Virgo mergers and rule out that PBHs (or any type of compact objects), in this mass range constitute a significant fraction of the dark matter.

We present a new method to obtain the interstellar radiation field (ISRF) of regions that contain stars, dust and gas. Starting from stellar spectra, we first run the radiative transfer code Skirt to compute the mean intensity field due to the stellar component alone. Then, we import that result to the photoionisation code Cloudy to compute the emissivity and opacity of the given mixture of gas and dust as a function of wavelength. Finally, we ask Skirt again to compute the mean intensity field, adding the total contribution of stars, gas and dust. This process is repeated iteratively, calling both codes sequentially in order to obtain increasingly accurate estimates. We have designed a first test, reminiscent of an HII region, that consists of a B star, approximated as a black-body, surrounded by a spherical shell of gas and dust with uniform density. We find that the results of our three-dimensional radiative transfer method are in excellent agreement with a spherically symmetric Cloudy simulation. As a realistic scientific application, we calculate the interstellar radiation field of a Milky Way-like galaxy based on two different chemical evolution models. Both of them give results broadly consistent with previous ones reported in the literature for the interstellar radiation field of our Galaxy, albeit they systematically underestimate the mid-infrared emission, with significant differences in this range, as well as in the ultraviolet, stemming from the input stellar and ISM properties. These results show the feasibility of our method to incorporate radiative transfer to chemical evolution models, increasing their predictive power and using this interstellar radiation field to further constrain their parameters. Python source code to implement our method is publicly available at https://github.com/MarioRomeroC/Mixclask.

Following the previous research on epicyclic oscillations of accretion disks around black holes (BHs) and neutron stars (NSs), a new model of high-frequency quasi-periodic oscillations (QPOs) has been proposed (CT model), which deals with oscillations of fluid in marginally overflowing accretion tori (i.e., tori terminated by cusps). According to preliminary investigations, the model provides better fits of the NS QPO data compared to the relativistic precession (RP) model. It also implies a significantly higher upper limit on the Galactic microquasars BH spin. A short analytic formula has been noticed to well reproduce the model's predictions on the QPO frequencies in Schwarzschild spacetimes. Here we derive an extended version of this formula that applies to rotating compact objects. We start with the consideration of Kerr spacetimes and derive a formula that is not restricted to a particular specific angular momentum distribution of the inner accretion flow, such as Keplerian or constant. Finally, we consider Hartle-Thorne spacetimes and include corrections implied by the NS oblateness. For a particular choice of a single parameter, our relation provides frequencies predicted by the CT model. For another value, it provides frequencies predicted by the RP model. We conclude that the formula is well applicable for rotating oblateness NSs and both models. We briefly illustrate application of our simple formula on several NS sources and confirm the expectation that the CT model is compatible with realistic values of the NS mass and provides better fits of data than the RP model.

It is shown that gas-dust perturbations in a disc with dust settling to the disc midplane exhibit the non-linear three-wave resonant interactions between streaming dust wave (SDW) and two inertial waves (IW). In the particular case considered in this paper, SDW at the wavenumber $k^\bullet = 2\kappa / (g_z t_s)$, where $\kappa$, $g_z$ and $t_s$ are, respectively, epicyclic frequency, vertical gravitational acceleration and particle's stopping time, interacts with two IW at the lower wavenumbers $k^\prime$ and $k^{\prime\prime}$ such that $k^\prime < k_{\rm DSI} < k^{\prime\prime} < k^\bullet$, where $k_{\rm DSI} = \kappa / (g_z t_s)$ is the wavenumber of the linear resonance between SDW and IW associated with the previously discovered linear dust settling instability. The problem is solved analytically in the limit of the small dust fraction. As soon as the dynamical dust back reaction on gas is taken into account, ${\bf k}^\bullet$, ${\bf k}^\prime$ and ${\bf k}^{\prime\prime}$ become slightly non-collinear and the emerging interaction of waves leads to simultaneous explosive growth of their amplitudes. This growth is explained by the conservative exchange with energy between the waves. The amplitudes of all three waves grow because the negative energy SDW transfers its energy to the positive energy IW. The product of the dimensionless amplitude of initially dominant wave and the time of explosion can be less than Keplerian time in a disc. It is shown that, generally, the three-wave resonance of an explosive type exists in a wide range of wavenumbers $0 < k^\bullet \leq 2\kappa / (g_z t_s)$. An explosive instability of gas-dust mixture may facilitate the dust clumping and the subsequent formation of planetesimals in young protoplanetary discs.

Residual foreground contamination by thermal Sunyaev-Zeldovich (tSZ) effect from galaxy clusters in cosmic microwave background (CMB) maps propagates into the reconstructed CMB lensing field, and thus biases the intrinsic cross-correlation between CMB lensing and large-scale structure (LSS). Through stacking analysis, we show that residual tSZ contamination causes an increment of lensing convergence in the central part of the clusters and a decrement of lensing convergence in the cluster outskirts. We quantify the impact of residual tSZ contamination on cross-correlations between the Planck 2018 CMB lensing convergence maps and the SDSS-IV galaxy density data through cross-power spectrum computation. In contrast with the Planck 2018 tSZ-deprojected SMICA lensing map, our analysis using the tSZ-contaminated SMICA lensing map measures a $\sim2.5\%$ negative bias at multipoles $\ell\lesssim 500$ and transits to a $\sim9\%$ positive bias at $\ell\gtrsim1500$, which validates earlier theoretical predictions of the overall shape of such tSZ-induced spurious cross-correlation. The tSZ-induced lensing convergence field in Planck CMB data is detected with more than $1\sigma$ significance at $\ell\lesssim 500$ and more than $14\sigma$ significance at $\ell\gtrsim1500$, yielding an overall $14.8\sigma$ detection. We also show that masking galaxy clusters in CMB data is not sufficient to eliminate the spurious lensing signal, still detecting a non-negligible bias with $5.5\sigma$ significance on cross-correlations with galaxy density fields. Our results emphasize how essential it is to deproject the tSZ effect from CMB maps at the component separation stage and adopt tSZ-free CMB lensing maps for cross-correlations with LSS data.

Decaying cold dark matter (CDM) has been considered as a mechanism to tackle the tensions in the Hubble expansion rate and the clustering of matter. However, polarization measurements of the cosmic microwave background (CMB) severely constrain the fraction of dark matter decaying before recombination, and lensing of the CMB anisotropies by large-scale structure set strong constraints on dark matter decaying after recombination. Together, these constraints make an explanation of the Hubble tension in terms of decaying dark matter unlikely. In response to this situation, we investigate whether a dark matter ensemble with CDM particles decaying into free streaming dark radiation in different epochs can alleviate the problem. We find that it does not.

We present the median-stacked Lyman-$\alpha$ surface brightness profile of 968 spectroscopically selected Lyman-$\alpha$ emitting galaxies (LAEs) at redshifts $1.9<z<3.5$ in the early data of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The selected LAEs are high-confidence Lyman-$\alpha$ detections with large signal-to-noise ratios observed with good seeing conditions (point-spread-function full-width-at-half-maximum $<1.4"$), excluding active galactic nuclei (AGN). The Lyman-$\alpha$ luminosities of the LAEs are $10^{42.4}-10^{43}\, \mathrm{erg}\, \mathrm{s}^{-1}$. We detect faint emission in the median-stacked radial profiles at the level of $(3.6\pm 1.3)\times 10^{-20}\,\mathrm{erg}\,\mathrm{s}^{-1}\,\mathrm{cm}^{-2}\,\mathrm{arcsec}^{-2}$ from the surrounding Lyman-$\alpha$ halos out to $r\simeq 160$ kpc (physical). The shape of the median-stacked radial profile is consistent at $r<80\,\mathrm{kpc}$ with that of much fainter LAEs at $3<z<4$ observed with the Multi Unit Spectroscopic Explorer (MUSE), indicating that the median-stacked Lyman-$\alpha$ profiles have similar shapes at redshifts $2<z<4$ and across a factor of $10$ in Lyman-$\alpha$ luminosity. While we agree with the results from the MUSE sample at $r<80\,\mathrm{kpc}$, we extend the profile over a factor of two in radius. At $r>80\,\mathrm{kpc}$, our profile is flatter than the MUSE model. The measured profile agrees at most radii with that of galaxies in the Byrohl et al. (2021) cosmological radiative transfer simulation at $z=3$. This suggests that the surface brightness of a Lyman-$\alpha$ halo at $r\lesssim 100$ kpc is dominated by resonant scattering of Lyman-$\alpha$ photons from star-forming regions in the central galaxy, whereas at $r > 100$ kpc it is dominated by photons from galaxies in surrounding dark matter halos.

We search for the signature of parity-violating physics in the Cosmic Microwave Background using Planck polarization data from the Public Release 4 (PR4 or $\mathtt{NPIPE}$). For nearly full-sky data, we initially find a birefringence angle $\beta=0.30^\circ\pm0.11^\circ$ ($68\%$~C.L.). We also find that the values of $\beta$ decrease as we enlarge the Galactic mask, which can be interpreted as the effect of polarized foreground emission. We use two independent approaches to model this effect and mitigate its impact on $\beta$. Although results are promising, and the good agreement between both models is encouraging, we do not assign cosmological significance to the measured value of $\beta$ until we improve our knowledge of the foreground polarization. Acknowledging that the miscalibration of polarization angles is not the only instrumental systematic that can create spurious TB and EB correlations, we also perform a detailed study of $\mathtt{NPIPE}$ end-to-end simulations to prove that our measurements of $\beta$ are not significantly affected by any of the known systematics.

We investigate the Hubble tension in the non-flat Super-$\Lambda {\rm CDM}$ model. The non-flat Super-$\Lambda {\rm CDM}$ model extends the Super-$\Lambda {\rm CDM}$ model by including the spatial curvature as a free parameter. The Super-$\Lambda {\rm CDM}$ model extends the standard $\Lambda {\rm CDM}$ model of cosmology through additional parameters accounting for the possible effect of a trispectrum in the primordial fluctuations. In the cosmic microwave background data, this effect can be parameterized using parameters that change the observed angular power spectrum from the theoretical power spectrum due to a trispectrum that couples long and short wavelength modes. In this work, we perform Markov Chain Monte Carlo (MCMC) data analysis on the recent Planck 2018 temperature and polarization fluctuations data and the local Hubble constant measurements from supernovae data assuming a non-flat Super-$\Lambda {\rm CDM}$ model. We find that there is a preference for non-zero values of the spatial curvature parameter $\Omega_k$ and the Super-$\Lambda {\rm CDM}$ parameter $A_0$ at a level of $\Delta \chi^2$ improvement of approximately 23.

Polarized component maps in the Northern Sky are derived from the QUIJOTE-MFI wide survey data at 11 and 13 GHz, the WMAP K and Ka bands and all Planck polarized channels using the parametric component separation method B-SeCRET. The addition of QUIJOTE-MFI data significantly improves the uncertainty in the parameter estimation of the low frequency dominant foreground, in particular the estimation of the synchrotron spectral index. We find statistically significant spatial variability across the sky. A power law model of the synchrotron emission provides a good fit of the data outside the galactic plane but fails to track the complexity of this region. Moreover, when we assume a synchrotron model with uniform curvature we find, in the 95% confidence region, a non-zero $c_s$ value. However, there is not sufficient statistical significance to determine which model is favoured.

Fast radio bursts (FRBs) represent one of the most exciting astrophysical discoveries of the recent past. The study of their low-frequency emission, which was only effectively picked up about ten years after their discovery, has helped shape the field thanks to some of the most important detections to date. Observations between 400 and 800 MHz, carried out by the CHIME/FRB telescope, in particular, have led to the detection of ~500 FRBs in little more than 1 year and, among them, ~20 repeating sources. Detections at low frequencies have uncovered a nearby population that we can study in detail via continuous monitoring and targeted campaigns. The latest, most important discoveries include: periodicity, both at the days level in repeaters and at the millisecond level in apparently non-repeating sources; the detection of an FRB-like burst from a galactic magnetar; and the localisation of an FRB inside a globular cluster in a nearby galaxy. The systematic study of the population at low frequencies is important for the characterisation of the environment surrounding the FRBs and, at a global level, to understand the environment of the local universe. This review is intended to give an overview of the efforts leading to the current rich variety of low-frequency studies and to put into a common context the results achieved in order to trace a possible roadmap for future progress in the field.

In astronomy, upcoming space telescopes with wide-field optical instruments have a spatially varying point spread function (PSF). Certain scientific goals require a high-fidelity estimation of the PSF at target positions where no direct measurement of the PSF is provided. Even though observations of the PSF are available at some positions of the field of view (FOV), they are undersampled, noisy, and integrated in wavelength in the instrument's passband. PSF modeling requires building a model from these observations that can infer a super-resolved PSF at any wavelength and any position in the FOV. Current data-driven PSF models can tackle spatial variations and super-resolution, but are not capable of capturing chromatic variations. Our model, coined WaveDiff, proposes a paradigm shift in the data-driven modeling of the point spread function field of telescopes. By adding a differentiable optical forward model into the modeling framework, we change the data-driven modeling space from the pixels to the wavefront. The proposed model relies on efficient automatic differentiation technology as well as modern stochastic first-order optimization techniques recently developed by the thriving machine-learning community. Our framework paves the way to building powerful models that are physically motivated and do not require special calibration data. This paper demonstrates the WaveDiff model on a simplified setting of a space telescope. The proposed framework represents a performance breakthrough with respect to existing data-driven approaches. The pixel reconstruction errors decrease 6-fold at observation resolution and 44-fold for a 3x super-resolution. The ellipticity errors are reduced by a factor of at least 20 and the size error by a factor of more than 250. By only using noisy broad-band in-focus observations, we successfully capture the PSF chromatic variations due to diffraction.

In the local universe, galaxies in clusters show different properties compared to more isolated systems. Understanding how this difference originates and whether it is already in place at high redshift is still a matter of debate. Thanks to uniquely deep optical spectra from the VANDELS survey, we investigate environmental effects on the stellar mass-metallicity relation (MZR) for a sample of ~1000 star-forming galaxies at 2<z<4. We complement our dataset with MOSFIRE follow-up of 21 galaxies to study the environmental dependence of the gas-phase MZR. Robust stellar and gas metallicities are derived, respectively, from well-calibrated photospheric absorptions features at 1501 and 1719 \AA in the stacked spectra, and from optical emission lines ([OII]3726-3729, [OIII]5007, and Hbeta) in individual systems. We characterize the environment through multiple criteria by using the local galaxy density maps previously derived in VANDELS. We find that environmental effects are weak at these redshifts, and more important around the densest overdensity structures, where galaxies have a lower stellar Z (by 0.2 dex) and a lower gas-phase Z (by 0.1 dex) compared to the field, with a significance of 1 and 2 sigma, respectively. Crucially, this offset cannot be explained by a selection effect due to a higher SFR, a fainter UV continuum, or different dust attenuations and stellar ages. Despite the still low S/N of our results, we propose a combination of increased mergers and high-speed encounters, more efficient AGN feedback in dense cores, and cold gas inflows as viable mechanisms diluting the metal content of overdense galaxies or expelling their metals to the IGM. Finally, some tensions remain with semi-analytic models and hydrodynamical simulations, which predict no significant offset as a function of host halo mass, suggesting that an explicit implementation of environmental processes is needed.

The current cosmological structure formation models predict that galaxies evolve through frequent mergers. During these events, the supermassive black holes (SMBHs) residing in the centres of the galaxies shrink to the central region, while losing energy via dynamical friction. Detection of SMBHs in these galaxy mergers is straightforward if they are actively accreting matter from their surroundings as active galactic nuclei (AGN). Currently, only a few dual AGN are known. One way to identify dual AGN candidates is by detecting double-peaked emission lines in their spectra. If these are broad spectral lines, it may indicate the existence of two broad line regions associated with two AGN. 2MASS J165939.7+183436 is a merging system, where the detected double-peaked broad emission lines can be explained by a dual AGN with an estimated separation of $0.085"$. Radio emission from this object was detected in the Faint Images of the Radio Sky at Twenty-centimeters survey with a flux density of 2.6 mJy. We used the European Very Long Baseline Interferometer Network and the enhanced Multi-Element Remotely Linked Interferometer Network to image 2MASS J165939.7+183436 at 1.7 GHz. We did not detect compact radio emission from the source.

We use $\sim 1,000$ X-ray sources in the COSMOS-Legacy survey and study the position of the AGN relative to the star forming main sequence (MS). We also construct a galaxy (non-AGN) reference sample that includes $\sim 90,000$ sources. We apply the same photometric selection criteria on both datasets and construct their spectral energy distributions (SEDs) using optical to far-infrared photometry compiled by the HELP project. We perform SED fitting, using the X-CIGALE algorithm and the same parametric grid for both datasets, to measure the star formation rate (SFR) and stellar mass of the sources. The mass completeness of the data is calculated at different redshift intervals and is applied on both samples. We define our own main sequence, based on the distributions of the specific SFR at different redshift ranges and exclude quiescent galaxies from our analysis. These allow us to compare the SFR of the two populations in a uniform manner, minimizing systematics and selection effects. Our results show that at low to moderate X-ray luminosities, AGN tend to have lower or, at most, equal star formation rates compared to non-AGN systems with similar stellar mass and redshift. At higher ($\rm L_{X,2-10keV} > 2-3\times 10^{44}\,erg\,s^{-1}$), we observe an increase of the SFR of AGN, for systems that have $\rm 10.5 < log\,[M_*(M_\odot)] < 11.5$.

The detection of exoplanets is hindered by the presence of complex astrophysical and instrumental noises. Given the difficulty of the task, it is important to ensure that the data are exploited to their fullest potential. In the present work, we search for an optimal exoplanet detection criterion. We adopt a general Bayesian multiple hypothesis testing framework, where the hypotheses are indexed by continuous variables. This framework is adaptable to the different observational methods used to detect exoplanets as well as other data analysis problems. We describe the data as a combination of several parametrized patterns and nuisance signals. We wish to determine which patterns are present, and for a detection to be valid, the parameters of the claimed pattern have to correspond to a true one with a certain accuracy. We search for a detection criterion minimizing false and missed detections, either as a function of their relative cost, or when the expected number of false detections is bounded. We find that if the patterns can be separated in a technical sense, the two approaches lead to the same optimal procedure. We apply it to the retrieval of periodic signals in unevenly sampled time series, emulating the search for exoplanets in radial velocity data. We show on a simulation that, for a given tolerance to false detections, the new criterion leads to 15 to 30\% more true detections than other criteria, including the Bayes factor.

An updated, binary-coded message has been developed for transmission to extraterrestrial intelligences in the Milky Way galaxy. The proposed message includes basic mathematical and physical concepts to establish a universal means of communication followed by information on the biochemical composition of life on Earth, the Solar System's time-stamped position in the Milky Way relative to known globular clusters, as well as digitized depictions of the Solar System, and Earth's surface. The message concludes with digitized images of the human form, along with an invitation for any receiving intelligences to respond. Calculation of the optimal timing during a given calendar year is specified for potential future transmission from both the Five-hundred-meter Aperture Spherical radio Telescope in China and the SETI Institute's Allen Telescope Array in northern California to a selected region of the Milky Way which has been proposed as the most likely for life to have developed. These powerful new beacons, the successors to the Arecibo radio telescope which transmitted the 1974 message upon which this expanded communication is in part based, can carry forward Arecibo's legacy into the 21st century with this equally well-constructed communication from Earth's technological civilization.

Polarization measurements of thermal radiation from magnetic white dwarf (MWD) stars have been proposed as a probe of axion-photon mixing. The radiation leaving the surface of the MWD is unpolarized, but if low-mass axions exist then photons polarized parallel to the direction of the MWD's magnetic field may convert into axions, which induces a linear polarization dependent on the strength of the axion-photon coupling $g_{a\gamma\gamma}$. We model this process by using the formalism of axion-photon mixing in the presence of strong-field vacuum birefringence to show that of all stellar types MWDs are the most promising targets for axion-induced polarization searches. We then consider linear polarization data from multiple MWDs, including SDSS J135141 and Grw+70$^\circ$8247, to show that after rigorously accounting for astrophysical uncertainties the axion-photon coupling is constrained to $|g_{a\gamma\gamma}| \lesssim 5.4 \times 10^{-12}$ GeV$^{-1}$ at 95% confidence for axion masses $m_a \lesssim 3 \times 10^{-7}$ eV. This upper limit puts in tension the previously-suggested explanation of the anomalous transparency of the Universe to TeV gamma-rays in terms of axions. We identify MWD targets for which future data and modeling efforts could further improve the sensitivity to axions.

Axion-like particles (ALPs) are a broad class of pseudo-scalar bosons that generically arise from broken symmetries in extensions of the standard model. In many scenarios, ALPs can mix with photons in regions with high magnetic fields. Photons from distant sources can mix with ALPs, which then travel unattenuated through the Universe, before they mix back to photons in the Milky Way galactic magnetic field. Thus, photons can traverse regions where their signals would normally be blocked or attenuated. In this paper, we study TeV $\gamma$-ray observations from distant blazars, utilizing the significant $\gamma$-ray attenuation expected from such signals to look for excess photon fluxes that may be due to ALP-photon mixing. We find no such excesses among a stacked population of seven blazars and constrain the ALP-photon coupling constant to fall below $\sim$3$\times$10$^{-11}$ GeV$^{-1}$ for ALP masses below 300 neV. These results are competitive with, or better than, leading terrestrial and astrophysical constraints in this mass range.

Using the exact WKB analysis of the higher-order differential equation, we analyze the mechanism of asymmetric preheating of a complex scalar field. The solution requires the Stokes phenomena of the fourth-order differential equation. We have identified the Stokes phenomena, which are crucial for the matter-antimatter asymmetry. The new Stokes lines and the virtual turning points, which are developed to analyze the Stokes phenomena of the higher-order differential equations, are crucial for global consistency but do not introduce a new effect in this model.

Gravitational-wave observations provide a unique opportunity to test General Relativity (GR) in the strong-field and highly-dynamical regime of the theory. Parameterized tests of GR are one well-known approach for quantifying violations of GR. This approach constrains deviations in the coefficients of the post-Newtonian (PN) phasing formula, which describes the gravitational-wave phase evolution of a compact binary as it inspirals. Current bounds from this test using LIGO/Virgo observations assume that binaries are circularized by the time they enter the detector frequency band. Here, we investigate the impact of residual binary eccentricity on the parameterized tests. We study the systematic biases in the parameter bounds when a phasing based on the circular orbit assumption is employed for a system that has some small residual eccentricity. We find that a systematic bias (for example, on the leading Newtonian deformation parameter) becomes comparable to the statistical errors for even moderate eccentricities of $\sim 0.04$ at $10$ Hz in LIGO/Virgo band. This happens at an even lower value of orbital eccentricity ($\sim 0.005$ at $10$ Hz) in the frequency band of third-generation (3G) detectors like Cosmic Explorer. These results demonstrate that incorporating physical effects like eccentricity in waveform models is important for accurately extracting science results from future detectors.

Minisuperpace Quantum Cosmology is an approach by which it is possible to infer initial conditions for dynamical systems which can suitably represent observable and non-observable universes. Here we discuss theories of gravity which, from various points of view, extend Einstein's General Relativity. Specifically, the Hamiltonian formalism for $f(R)$, $f(T)$ and $f(\mathcal{G})$ gravity, with $R$, $T$, and $\mathcal{G}$ being the curvature, torsion and Gauss--Bonnet scalars, respectively, is developed starting from the Arnowitt-Deser-Misner approach. The Minisuperspace Quantum Cosmology is derived for all these models and cosmological solutions are obtained thanks to the existence of Noether symmetries. The Hartle criterion allows the interpretation of solutions in view of observable universes.

This essay is a nontechnical primer for a broader audience, in which I paint a broad-brush picture of modern cosmology. I begin by reviewing the evidence for the big bang, including the expansion of our Universe, the cosmic microwave background, and the primordial abundances of the light elements. Next, I discuss how these and other cosmological observations can be well explained by means of the concordance model of cosmology, putting a particular emphasis on the composition of the cosmic energy budget in terms of visible matter, dark matter, and dark energy. This sets the stage for a short overview of the history of the Universe from the earliest moments of its existence all the way to the accelerated expansion at late times and beyond. Finally, I summarize the current status of the field, including the challenges it is currently facing such as the Hubble tension, and conclude with an outlook onto the bright future that awaits us in the coming years and decades. The text is complemented by an extensive bibliography serving as a guide for readers who wish to delve deeper.

The stability and equation of state of quark matter are studied within both two-flavor and (2+1)-flavor Nambu-Jona-Lasinio (NJL) models including the vector interactions. With a free parameter $\alpha$, the Lagrangian is constructed by two parts, the original NJL Lagrangian and the Fierz transformation of it, as $L=(1-\alpha) L_{\rm{NJL}}+\alpha L_{\rm{Fierz}}$. We find that there is a possibility for both $ud$ nonstrange and $uds$ strange matter being absolute stable, depending on the interplay of the confinement with quark vector interaction and the exchange interaction channels. The calculated quark star properties can reconcile with the recently measured masses and radii of PSR J0030+0451 and PSR J0740+6620, as well as the tidal deformability of GW170817. Furthermore, the more strongly-interacting quark matter in the nonstrange stars allows a stiffer equation of state and consequently a higher maximum mass ($\sim2.7\, M_{\odot}$) than the strange ones ($\sim2.1\, M_{\odot}$). The sound velocities in strange and nonstrange quark star matter are briefly discussed compared to those of neutron star matter.

For sufficiently wide ranges of applied control signals (control voltages), MEMS and piezoelectric Deformable Mirrors (DMs), exhibit nonlinear behavior. The nonlinear behavior manifests itself in nonlinear actuator couplings, nonlinear actuator deformation characteristics, and in the case of piezoelectric DMs, hysteresis. Furthermore, in a number of situations, DM behavior can change over time, and this requires a procedure for updating the DM models on the basis of the observed data. If not properly modeled and if not taken into account when designing control algorithms, nonlinearities, and time-varying DM behavior, can significantly degrade the achievable closed-loop performance of Adaptive Optics (AO) systems. Widely used approaches for DM control are based on pre-estimated linear time-invariant DM models in the form of influence matrices. Often, these models are not being updated during system operation. Consequently, when the nonlinear DM behavior is excited by control signals with wide operating ranges, or when the DM behavior changes over time, the state-of-the-art DM control approaches relying upon linear control methods, might not be able to produce a satisfactory closed-loop performance of an AO system. Motivated by these key facts, we present a novel method for data-driven DM control. Our approach combines a simple open-loop control method with a recursive least squares method for dynamically updating the DM model. The DM model is constantly being updated on the basis of the dynamically changing DM operating points. That is, the proposed method updates both the control actions and the DM model during the system operation. We experimentally verify this approach on a Boston Micromachines MEMS DM with 140 actuators. Preliminary experimental results reported in this manuscript demonstrate good potential for using the developed method for DM control.

The Big Bang initial singularity problem can be solved by means of bouncing solutions. In the context of extended theories of gravity, we will look for covariant effective actions whose field equations contain up to fourth-order derivatives of the metric tensor. In finding such bouncing solutions, we will make use of an order reduction technique based on a perturbative approach. Reducing the order of the field equations to second-order, we are able to find solutions which are perturbatively close to General Relativity. We will build the covariant effective actions of the resulting order reduced theories.