Papers for Thursday, Jan 19 2023

The origin of stellar-mass black hole mergers discovered through gravitational waves is being widely debated. Mergers in the disks of active galactic nuclei (AGN) represent a promising source of origin, with possible observational clues in the gravitational wave data. Beyond gravitational waves, a unique signature of AGN-assisted mergers is electromagnetic emission from the accreting black holes. Here we show that jets launched by accreting black holes merging in an AGN disk can be detected as peculiar transients by infrared, optical, and X-ray observatories We further show that this emission mechanism can explain the possible associations between gravitational wave events and the optical transient ZTF19abanrhr and the proposed gamma-ray counterparts GW150914-GBM and LVT151012-GBM. We demonstrate how these associations, if genuine, can be used to reconstruct the properties of these events' environments. Searching for infrared and X-ray counterparts to similar electromagnetic transients in the future, once host galaxies are localized by optical observations, could provide a smoking gun signature of the mergers' AGN origin.

The flattened, possibly co-rotating plane of satellite galaxies around Centaurus A, if more than a fortuitous alignment, adds to the pre-existing tension between the well-studied Milky Way and M31 planes and the $\Lambda$CDM model of structure formation. It was recently reported that the Centaurus A satellite plane (CASP) may be rotationally supported, but a further understanding of the system's kinematics is elusive in the absence of full three-dimensional velocities. We constrain the transverse velocities of 27 satellites that would rotationally stabilise the Centaurus A plane, and classify the satellites by whether their possible orbits are consistent with the CASP. Five satellites are identified to be unlikely to participate in the plane, two of which are clearly non-members. Despite their previously reported line-of-sight velocity trend suggestive of a common co-rotating motion, 17 out of 22 potential CASP members are consistent with either orbital direction within both the full range of possible kinematics as well as when limiting orbits to those within the plane. On the other hand, disregarding the 5 off-plane satellites found to be inconsistent with CASP membership enhances the significance of the CASP's line-of-sight velocity trend fivefold. Our results are robust with different mass estimates of the Centaurus A halo, and the adoption of either spherical or triaxial NFW potentials.

Ultralight dark matter (ULDM) is usually taken to be a single scalar field. Here we explore the possibility that ULDM consists of $N$ light scalar fields with only gravitational interactions. This configuration is more consistent with the underlying particle physics motivations for these scenarios than a single ultralight field. ULDM halos have a characteristic granular structure that increases stellar velocity dispersion and can be used as observational constraints on ULDM models. In multifield simulations, we find that inside a halo the amplitude of the total density fluctuations decreases as $1/\sqrt{N}$ and that the fields do not become significantly correlated over cosmological timescales. Smoother halos heat stellar orbits less efficiently, reducing the velocity dispersion relative to the single field case and thus weakening the observational constraints on the field mass. Analytically, we show that for $N$ equal-mass fields with mass $m$ the ULDM contribution to the stellar velocity dispersion scales as $1/(N m^3)$. Lighter fields heat the most efficiently and if the smallest mass $m_L$ is significantly below the other field masses the dispersion scales as $1/(N^2 m_L^3)$.

General relativity and its Newtonian weak field limit are not sufficient to explain the observed phenomenology in the Universe, from the formation of large-scale structures to the dynamics of galaxies, with the only presence of baryonic matter. The most investigated cosmological model, the $\Lambda$CDM, accounts for the majority of observations by introducing two dark components, dark energy and dark matter, which represent $\sim$95% of the mass-energy budget of the Universe. Nevertheless, the $\Lambda$CDM model faces important challenges on the scale of galaxies. For example, some very tight relations between the properties of dark and baryonic matters in disk galaxies, such as the baryonic Tully-Fisher relation (BTFR), the mass discrepancy-acceleration relation (MDAR), and the radial acceleration relation (RAR), which see the emergence of the acceleration scale $a_0 \simeq 1.2 \times 10^{-10}$ m s$^{-2}$, cannot be intuitively explained by the CDM paradigm, where cosmic structures form through a stochastic merging process. An even more outstanding coincidence is due to the fact that the acceleration scale $a_0$, emerging from galaxy dynamics, also seems to be related to the cosmological constant $\Lambda$. Another challenge is provided by dwarf galaxies, which are darker than what is expected in their innermost regions. These pieces of evidence can be more naturally explained, or sometimes even predicted, by modified theories of gravity, that do not introduce any dark fluid. I illustrate possible solutions to these problems with the modified theory of gravity MOND, which departs from Newtonian gravity for accelerations smaller than $a_0$, and with Refracted Gravity, a novel classical theory of gravity introduced in 2016, where the modification of the law of gravity is instead regulated by a density scale.

Minerals are solid state nuclear track detectors - nuclear recoils in a mineral leave latent damage to the crystal structure. Depending on the mineral and its temperature, the damage features are retained in the material from minutes (in low-melting point materials such as salts at a few hundred degrees C) to timescales much larger than the 4.5 Gyr-age of the Solar System (in refractory materials at room temperature). The damage features from the $O(50)$ MeV fission fragments left by spontaneous fission of $^{238}$U and other heavy unstable isotopes have long been used for fission track dating of geological samples. Laboratory studies have demonstrated the readout of defects caused by nuclear recoils with energies as small as $O(1)$ keV. This whitepaper discusses a wide range of possible applications of minerals as detectors for $E_R \gtrsim O(1)$ keV nuclear recoils: Using natural minerals, one could use the damage features accumulated over $O(10)$ Myr$-O(1)$ Gyr to measure astrophysical neutrino fluxes (from the Sun, supernovae, or cosmic rays interacting with the atmosphere) as well as search for Dark Matter. Using signals accumulated over months to few-years timescales in laboratory-manufactured minerals, one could measure reactor neutrinos or use them as Dark Matter detectors, potentially with directional sensitivity. Research groups in Europe, Asia, and America have started developing microscopy techniques to read out the $O(1) - O(100)$ nm damage features in crystals left by $O(0.1) - O(100)$ keV nuclear recoils. We report on the status and plans of these programs. The research program towards the realization of such detectors is highly interdisciplinary, combining geoscience, material science, applied and fundamental physics with techniques from quantum information and Artificial Intelligence.

One of the primary goals when studying galaxy formation is to understand how the luminous component of the Universe, galaxies, relates to the growth of structure which is dominated by the gravitational collapse of dark matter haloes. The stellar-to-halo mass relation probes how galaxies occupy dark matter haloes and what that entails for their star formation history. We deliver the first self-consistent empirical model that can place constraints on the stellar-to-halo mass relation down to log stellar mass $\log_{10}(m^*/\mathrm{M}_{\odot}) \leq 5.0$ by fitting our model directly to Local Group dwarf data. This is accomplished by penalising galaxy growth in late-forming, low-mass haloes by mimicking the effects of reionization. This process serves to regulate the number density of galaxies by altering the scatter in halo peak mass $M^{\mathrm{peak}}_{h}$ at fixed stellar mass, creating a tighter scatter than would otherwise exist without a high-$z$ quenching mechanism. Our results indicate that the previously established double-power law stellar-to-halo mass relation can be extended to include galaxies with $\log_{10}(M^{\mathrm{peak}}_{\mathrm{h}}/\mathrm{M}_{\odot})\gtrsim 10.0$. Furthermore, we show that haloes with $\log_{10}(M^{\mathrm{peak}}_{\mathrm{h}}/\mathrm{M}_{\odot})\lesssim 9.3$ by $z=4$ are unlikely to host a galaxy with $\log_{10}(m^*/\mathrm{M}_{\odot}) > 5.0$.

The abundance of carbon relative to oxygen (C/O) is a promising probe of star formation history in the early universe, as these elements are produced on different timescales. We present a measurement of $\log{\mathrm{(C/O)}} = -1.01\pm0.12$ (stat) $\pm0.10$ (sys) in a $z=6.23$ galaxy observed as part of the GLASS-JWST Early Release Science Program. Notably, we achieve good precision thanks to the detection of the rest-frame ultraviolet O III], C III], and C IV emission lines delivered by JWST/NIRSpec. The C/O abundance is $\sim$0.8 dex lower than the solar value and is consistent with the expected yield from core-collapse supernovae, indicating negligible carbon enrichment from intermediate-mass stars. This in turn implies rapid buildup of a young stellar population with age $\lesssim$100 Myr in a galaxy seen $\sim$900 million years after the Big Bang. Our chemical abundance analysis is consistent with spectral energy distribution modeling of JWST/NIRCam photometric data, which indicates a current stellar mass $\log\,\mathrm{M}_* / \mathrm{Msun} = 8.4^{+0.4}_{-0.2}$ and specific star formation rate sSFR $\simeq 20$ Gyr$^{-1}$. These results showcase the value of chemical abundances and C/O in particular to study the earliest stages of galaxy assembly.

The evolution of the luminosity function (LF) of Active Galactic Nuclei (AGNs) at $z \gtrsim 5$ represents a key constraint to understand their contribution to the ionizing photon budget necessary to trigger the last phase transition in the Universe, i.e. the epoch of Reionization. Recent searches for bright high-z AGNs suggest that the space densities of this population at $z>4$ has to be revised upwards, and sparks new questions about their evolutionary paths. Gas accretion is the key physical mechanism to understand both the distribution of luminous sources and the growth of central Super-Massive Black Holes (SMBHs). In this work, we model the high-z AGN-LF assuming that high-z luminous AGN shine at their Eddington limit: we derive the expected evolution as a function of the ``duty-cycle'' ($f_{\rm dc}$), i.e. the fraction of life-time that a given SMBH spends accreting at the Eddington rate. Our results show that intermediate values ($f_{\rm dc} \simeq 0.1$) predict the best agreement with the ionizing background and photoionization rate, but do not provide enough ionizing photons to account for the observed evolution of the hydrogen neutral fraction. Smaller values ($f_{\rm dc} \lesssim 0.05$) are required for AGNs to be the dominant population responsible for Hydrogen reionization in the Early Universe. We then show that this low-$f_{\rm dc}$ evolution can be reconciled with the current constraints on Helium reionization, although it implies a relatively large number of inactive SMBHs at $z\gtrsim5$, in tension with SMBH growth models based on heavy seeding.

Massive stars are thought to be progenitors of Long Gamma Ray Bursts (GRBs), most likely with a bias favouring low metallicity progenitors. Because galaxies do not have a constant metallicity throughout, the combination of line-of-sight absorption metallicity inferred from GRB afterglow spectroscopy and of host galaxy global metallicity derived from emission lines diagnostics represents a powerful way to probe both the bias function for GRB progenitors, and the chemical inhomogeneities across star forming regions. In this study, we predict the relationship between Zabs and Zem using three different hydrodynamical cosmological simulations: Illustris, EAGLE, and IllustrisTNG. We find that while the qualitative shape of the curve relating emission versus absorption metallicity remains the same, the predicted relationship between these two observables is significantly different between the simulations. Using data for the host galaxy of GRB121024A for which both Zabs and Zem have been measured, we find marginal support for the Illustris simulation as producing the most-realistic internal metallicity distributions within star-forming galaxies at cosmic noon. Overall, all simulations predict similar properties for the bulk of the GRB host galaxy population, but each has distinct features in the tail of the Zabs-Zem distribution that in principle allow to discriminate between models if a sufficiently large sample of observations are available (i.e. N>11 on average). Substantial progress is expected in the near future, with upcoming JWST/NIRspec observations of 10 GRB host galaxies for which absorption metallicity from the afterglow spectra exists.

Black widows are extreme millisecond pulsar binaries where the pulsar wind ablates their low-mass companion stars. Their optical light curves vary periodically due to the high irradiation and tidal distortion of the companion, which allows us to infer the binary parameters. We present simultaneous multi-band observations obtained with the HIPERCAM instrument at the 10.4-m GTC telescope for six of these systems. The combination of this five-band fast photometer with the world's largest optical telescope enables us to inspect the light curve range near minima. We present the first light curve for PSR J1641+8049, as well as attain a significant increase in signal-to-noise and cadence compared with previous publications for the remaining 5 targets: PSR J0023+0923, PSR J0251+2606, PSR J0636+5129, PSR J0952-0607 and PSR J1544+4937. We report on the results of the light curve modelling with the Icarus code for all six systems, which reveals some of the hottest and densest companion stars known. We compare the parameters derived with the limited but steadily growing black widow population for which optical modelling is available. We find some expected correlations, such as that between the companion star mean density and the orbital period of the system, but also a puzzling positive correlation between the orbital inclination and the irradiation temperature of the companion. We propose such a correlation would arise if pulsars with magnetic axis orthogonal to their spin axis are capable of irradiating their companions to a higher degree.

The hard X-ray power law, prominent in the hard state in black hole X-ray binaries, is generally due to thermal Comptonization in the corona. Optically thin synchrotron emission from compact jets is commonly seen at infrared wavelengths in the hard state. The extent of this spectrum to higher energies remains uncertain. Here, a multi-wavelength study of GX 339-4 is presented. The IR to X-ray spectral index is measured and compared to the X-ray spectral index fitted separately. On some dates in which the jet dominates the IR emission, the X-ray power law and the IR to X-ray power law spectral indices are both in the range alpha = -0.7 +/- 0.2 (where F_nu ~ nu^alpha), i.e. photon index, Gamma = 1.7 +/- 0.2. This suggests they could be the same power law with the same origin, or that this is a coincidence. On other dates in the hard state, alpha_{IR-X} < alpha_{X}, ruling out a common origin. It is likely that Comptonization dominates on most dates, as expected. However, the X-ray power law never appears to be fainter than the jet power law extrapolated from IR to X-ray, implying that the jet contribution imposes a lower limit to the X-ray flux. If confirmed, this would imply the cooling break in the synchrotron spectrum probably resides at X-ray or higher energies. It is suggested that X-ray spectral fitting should include an extra power law with a break (ideally fit to IR too).

Pulsar timing arrays (PTAs) are Galactic-scale gravitational wave (GW) detectors consisting of precisely-timed pulsars distributed across the sky. Within the decade, PTAs are expected to detect the nanohertz GWs emitted by close-separation supermassive black hole binaries (SMBHBs), thereby opening up the low frequency end of the GW spectrum for science. Individual SMBHBs which power active galactic nuclei are also promising multi-messenger sources; they may be identified via theoretically predicted electromagnetic (EM) signatures and be followed up by PTAs for GW observations. In this work, we study the detection and parameter estimation prospects of a PTA which targets EM-selected SMBHBs. Adopting a simulated Galactic millisecond pulsar population, we envisage three different pulsar timing campaigns which observe three mock sources at different sky locations. We find that an all-sky PTA which times the best pulsars is an optimal and feasible approach to observe EM-selected SMBHBs and measure their source parameters to high precision (i.e., comparable to or better than conventional EM measurements). We discuss the implications of our findings in the context of the future PTA experiment with the planned Deep Synoptic Array-2000 and the multi-messenger studies of SMBHBs such as the well-known binary candidate OJ 287.

We present the analysis of multiepoch observations of a set of 12 variable, Compton-thin, local (z<0.1) active galactic nuclei (AGN) selected from the 100-month BAT catalog. We analyze all available X-ray data from \chandra, \xmm, and \nustar, adding up to a total of 53 individual observations. This corresponds to between 3 and 7 observations per source, probing variability timescales between a few days and $\sim 20$~yr. All sources have at least one \nustar observation, ensuring high-energy coverage, which allows us to disentangle the line-of-sight and reflection components in the X-ray spectra. For each source, we model all available spectra simultaneously, using the physical torus models \myt, \bor, and \uxc. The simultaneous fitting, along with the high-energy coverage, allows us to place tight constraints on torus parameters such as the torus covering factor, inclination angle, and average column density. We also estimate the line-of-sight column density ($N_{\rm H}$) for each individual observation. Within the 12 sources, we detect clear line-of-sight $N_{\rm H}$ variability in 5, non-variability in 5, and for 2 of them it is not possible to fully disentangle intrinsic-luminosity and $N_{\rm H}$ variability. We observe large differences between the average values of line-of-sight $N_{\rm H}$ (or $N_{\rm H}$ of the obscurer) and the average $N_{\rm H}$ of the torus (or $N_{\rm H}$ of the reflector), for each source, by a factor between $\sim2$ to $>100$. This behavior, which suggests a physical disconnect between the absorber and the reflector, is more extreme in sources that present $N_{\rm H}$ variability. $N_{\rm H}$-variable AGN also tend to present larger obscuration and broader cloud distributions than their non-variable counterparts. We observe that large changes in obscuration only occur at long timescales, and use this to place tentative lower limits on torus cloud sizes.

The unified model of active galactic nuclei (AGN) includes a toroidal obscuring structure to explain the differences between Type I and Type II AGN as an effect of inclination angle. This toroidal structure is thought to be 'clumpy' as the line-of-sight column density, $N_{H}$, has been observed to vary with time in many sources. We present a new method which uses a variation in hardness ratio to predict whether an AGN will have experienced $N_H$ variability across different observations. We define two sets of hard and soft bands that are chosen to be sensitive to the energies most affected by changes in $N_H$. We calculate these ratios for Chandra and XMM-Newton observations on a sample of 12 sources with multiple observations, and compare the predictions of this method with the $N_H$ values obtained from spectral fitting. We find that the method proposed in this work is effective in preselecting sources for variability studies.

The growing number of exoplanet discoveries and advances in machine learning techniques allow us to find, explore, and understand characteristics of these new worlds beyond our Solar System. We analyze the dataset of 762 confirmed exoplanets and eight Solar System planets using efficient machine-learning approaches to characterize their fundamental quantities. By adopting different unsupervised clustering algorithms, the data are divided into two main classes: planets with $\log R_{p}\leq0.91R_{\oplus}$ and $\log M_{p}\leq1.72M_{\oplus}$ as class 1 and those with $\log R_{p}>0.91R_{\oplus}$ and $\log M_{p}>1.72M_{\oplus}$ as class 2. Various regression models are used to reveal correlations between physical parameters and evaluate their performance. We find that planetary mass, orbital period, and stellar mass play preponderant roles in predicting exoplanet radius. The validation metrics (RMSE, MAE, and $R^{2}$) suggest that the Support Vector Regression has, by and large, better performance than other models and is a promising model for obtaining planetary radius. Not only do we improve the prediction accuracy in logarithmic space, but also we derive parametric equations using the M5P and Markov Chain Monte Carlo methods. Planets of class 1 are shown to be consistent with a positive linear mass-radius relation, while for planets of class 2, the planetary radius represents a strong correlation with their host stars' masses.

External accretion events such as a galaxy merger or the accretion of gas from the immediate environment of a galaxy, can create a large misalignment between the gas and the stellar kinematics. Numerical simulations have suggested that misaligned structures may promote the inflow of gas to the nucleus of the galaxy and the accretion of gas by the central supermassive black hole. We show for the first time that galaxies with a strong misalignment between the ionised gas and stellar kinematic angles have a higher observed fraction of active black holes than galaxies with aligned rotation of gas and stars. The increase in black hole activity suggests that the process of formation and/or presence of misaligned structures is connected with the fuelling of active supermassive black holes.

Attempts to understand the dynamics of warped astrophysical discs have garnered significant attention, largely motivated by the growing catalogue of observed distorted systems. Previous studies have shown that the evolution of the warp is crucially regulated by the internal flow fields established by the undulating geometry. These are typically modelled as laminar horizontal, shearing flows which oscillate back and forth at approximately the orbital frequency. However this shearing motion is known to be susceptible to a hydrodynamic, parametric instability of inertial waves which might modify the warped dynamics. Whilst the linear growth phase is well understood, the subsequent nonlinear saturation combined with the self-consistent feedback onto the warp has not been studied. In this work, we implement a novel numerical setup using the recent ring model framework of Fairbairn and Ogilvie, within the Lagrangian code GIZMO. We formally identify several locally growing modes in the simulation, as predicted by a three-mode coupling analysis of the instability, and find decent agreement with the theoretical growth rates. We understand the saturation mechanism as a wave breaking process which suppresses the growth of shorter wavelength parametric couplings first, whilst allowing the longest mode to dominate the final quasi-steady, wavelike turbulence. The Reynolds stresses, transporting energy from the warp to the small scales, can be effectively modelled using a time-dependent, anisotropic viscous alpha model which closely captures the amplitude and phase evolution of the warp. Consequently, this model might help inform future global studies which are commonplace but typically don't resolve the parametric instability.

The ultra-luminous infrared galaxy Arp2 20 is a late-stage merger with several tidal structures in the outskirts and two very compact, dusty nuclei that show evidence for extreme star formation and host at least one AGN. New and archival high-resolution images taken by the Hubble Space Telescope provide a state-of-the-art view of the structures, dust, and stellar clusters in Arp 220. We find that ~90% of the Halpha emission arises from a shock-ionized bubble emanating from the AGN in the western nucleus, while the nuclear disks dominate the Pbeta emission. Four very young (~3-6 Myr) but lower mass (< 10^4 Msun) clusters are detected in Halpha within a few arcsec of the nuclei, but produce less than 1% of the line emission. We see little evidence for a population of massive clusters younger than 100Myr anywhere in Arp 220. From the masses and ages of the detected clusters, we find that star formation took place more-or-less continuously starting ~few Gyr ago with a rate between ~3-12 Msun/yr. Approximately 100Myr ago, star formation shut off suddenly everywhere, except in the nuclear disks. A very recent flicker of weak star formation produced the four young, low-mass clusters, while the rest of the galaxy appears to have remained in a post-starburst state. Cluster ages indicate that the tidal structures on the west side of the galaxy are older than those on the east side, but all appear to pre-date the shutoff of star formation. Arp 220 has many of the characteristics expected of a 'Shocked Post-Starburst Galaxy' or SPOG, since most of the system has been in a post-starburst state for the past ~100Myr and the detected Halpha emission arises from shocked rather than photo-ionized gas.

The fifth iteration of the Sloan Digital Sky Survey (SDSS-V) is set to obtain optical and near-infrared spectra of $\sim$5 million stars of all ages and masses throughout the Milky Way. As a part of these efforts, APOGEE & BOSS Young Star Survey (ABYSS) will observe $\sim10^5$ stars with ages $<$30 Myr that have been selected using a set of homogeneous selection functions that make use of different tracers of youth. The ABYSS targeting strategy we describe in this paper is aimed to provide the largest spectroscopic census of young stars to-date. It consists of 8 different types of selection criteria that take the position on the HR diagram, infrared excess, variability, as well as the position in phase space in consideration. The resulting catalog of $\sim$200,000 sources (of which a half are expected to be observed) provides representative coverage of the young Galaxy, including both nearby diffuse associations as well as more distant massive complexes, reaching towards the inner Galaxy and the Galactic center.

We study point-like explosive events (EE), characterized by emission in the far wings of spectral lines, in a quiet region near the South Pole, using Interface Region Imaging Spectrograph (IRIS) spectra at two slit positions, slit-jaw (SJ) observations and Atmospheric Imaging Assembly (AIA) images. The events were best visible in SiIV spectra; they were weak in SJs, occasionally visible in 1600 A and 304 A AIA images, and invisible in higher temperature AIA images. We identified EEs from position--time images in the far wings of the SiIV lines and measured their distance from the limb. A Gaussian model of the height distribution showed that EEs occur in a narrow (0.9") height range, centered at 3.2" above the continuum limb at 2832.0 A. On the disk, we found that they occur in network boundaries. Further, we studied the line profiles of two bright EEs above the limb and one on the disk. We found that what appears as broad-band emission, is actually a superposition of 2--3 narrow-band Gaussian components with well-separated line profiles, indicating that material is expelled towards and/or away from the observer in discrete episodes in time and in space. The expelled plasma accelerates quickly, reaching line-of-sight (LOS) velocities up to 90 km/s. Overall, the motion was practically along the LOS, as the velocity on the plane of sky was small. In some cases tilted spectra were observed that could be interpreted in terms of rotating motions of up to 30 km/s. We did not find any strong absorption features in the wing of the SiIV lines, although in one case a very weak absorption feature was detected. No motions, indicative of jets, were detected in SJ or AIA images. Reconnection in an asymmetric magnetic-field geometry, in the middle or near the top of small loops is a plausible explanation of their observational characteristics.

3C 223 is a radio loud, Type 2 quasar at $z=0.1365$ with an intriguing XMM-Newton spectrum that implicated it as a rare, Compton-thick ($N_{\rm H} \gtrsim 1.25 \times 10^{24}$ cm$^{-2}$) active galactic nucleus (AGN). We obtained contemporaneous XMM-Newton and NuSTAR spectra to fit the broad-band X-ray spectrum with the physically-motivated MYTorus and borus02 models. We confirm earlier results that the obscuring gas is patchy with both high (though not Compton-thick) levels of obscuration ($N_{\rm H} > 10^{23}$ cm$^{-2}$) and gas clouds with column densities up to an order of magnitude lower. The spectral fitting results indicate additional physical processes beyond those modeled in the spectral grids of MYTorus and borus02 impact the emergent spectrum: the Compton-scattering region may be extended beyond the putative torus; a ring of heavy Compton-thick material blocks most X-ray emission along the line of sight; or the radio jet is beamed, boosting the production of Fe K$\alpha$ line photons in the global medium compared with what is observed along the line of sight. We revisit a recent claim that no radio loud Compton-thick AGN have yet been conclusively shown to exist, finding three reported cases of radio loud AGN with global average (but not line-of-sight) column densities that are Compton-thick. Now that it is possible to separately determine line-of-sight and global column densities, inhomogeneity in the obscuring medium has consequences for how we interpet the spectrum and classify an AGN as "Compton-thick."

We present a set of recommended best practices for JWST data collection for members of the community focussed on the direct imaging and spectroscopy of exoplanetary systems. These findings and recommendations are based on the early analysis of the JWST Early Release Science Program 1386, "High-Contrast Imaging of Exoplanets and Exoplanetary Systems with JWST." Our goal is for this information to be useful for observers in preparation of JWST proposals for Cycle 2 and beyond. In addition to compiling a set of best practices from our ERS program, in a few cases we also draw on the expertise gained within the instrument commissioning programs, as well as include a handful of data processing best practices. We anticipate that this document will be regularly updated and resubmitted to arXiv.org to ensure that we have distributed our knowledge of best-practices for data collection as widely and efficiently as possible.

We study the statistical properties of the eigenvalues of the primordial tidal and deformation tensor for random Gaussian cosmic density fields. With the tidal and deformation tensors, Hessians of the gravitational and velocity potential, being Gaussian, the corresponding eigenvalue fields are distinctly non-Gaussian. Following the extension of the Doroshkevich formula for the joined distribution of eigenvalues to two-dimensional fields, we evaluate the two- and three-point correlation functions of the eigenvalue fields. In addition, we assess the number densities of singular points of the eigenvalue fields and find their corresponding two- and three-point correlation functions. The role of tidal forces and the resulting mass element deformation in shaping the prominent anisotropic wall-like and filamentary components of the cosmic web has since long been recognized based on the Zel'dovich approximation. Less well-known is that the weblike spatial pattern is already recognizable in the primordial tidal and deformation eigenvalue field, even while the corresponding Gaussian density and the potential field appear merely as a spatially incoherent and unstructured random field. Furthermore, against the background of a full phase-space assessment of structure formation in the Universe, the caustic skeleton theory entails a fully analytical framework for the nonlinear evolution of the cosmic web. It describes the folding of the dark matter sheet and the emerging caustic singularities, fully specified by the deformation eigenvalues and eigenvectors. Finally, tidal tensor eigenvalues are of central importance, and understanding their distribution is critical in predicting the resulting rotation and orientation. The current study applies to two-dimensional Gaussian random fields and will be generalized to a three-dimensional analysis in an upcoming study.

Soon after the Gaia data release (DR) 3 in June 2022, some candidates (and one confirmed) of detached black hole (BH) - luminous companion (LC) binaries have been reported. Existing and future detections of astrometric BH-LC binaries will shed light on the spatial distribution of these systems, which can deepen our understanding of the natal kicks and the underlying formation mechanism of BHs. By tracking Galactic orbits of BH-LC binaries obtained from BSE, we find that distributions of BH mass and the height from the Galactic plane |z| would help us give a constraint on supernova model. We also indicate that the correlations of (i) orbital periods and eccentricities, and (ii) BH mass and |z| could be clues for the strength of natal kick. We also discuss the possibility of forming BH-LC binaries like the BH binary candidates reported in Gaia DR3 and Gaia BH 1, finding that if the candidates as well as the confirmed binary originate from isolated binaries, they favor models which produce low-mass BHs and have high common envelope efficiencies exceeding unity.

We present criteria for separately classifying stars and unresolved background galaxies in photometric catalogs generated with the point spread function (PSF) fitting photometry software DOLPHOT from images taken of Draco II, WLM, and M92 with the Near Infrared Camera (NIRCam) on JWST. Photometric quality metrics from DOLPHOT in one or two filters can recover a pure sample of stars. Conversely, colors formed between short-wavelength (SW) and long-wavelength (LW) filters can be used to effectively identify pure samples of galaxies. Our results highlight that the existing DOLPHOT output parameters can be used to reliably classify stars in our NIRCam data without the need to resort to external tools or more complex heuristics.

The Sun reveals itself in the 385.8-2.439-nHz band of polar ({\phi}Sun>|70{\deg}|) fast (>700 km s^-1) solar wind's decade-scale dynamics as a globally completely vibrating, revolving-field magnetoalternator rather than a proverbial engine. Thus North-South separation of 1994-2008 Ulysses <10 nT wind polar samplings spanning ~1.6 10^7-2.5 10^9-erg base energies reveals Gauss-Vanicek spectral signatures of an entirely >99%-significant Sun-borne global sharp Alfven resonance (AR), Pi=PS/i, imprinted into the winds to the order n=100+ and co-triggered by the PS=~11-yr Schwabe global mode northside, its ~10-yr degeneration equatorially, and ~9-yr degeneration southside. The Sun is a typical ~3-dB-attenuated ring-system of differentially rotating and contrarily (out-of-phase) vibrating conveyor belts and layers with a continuous spectrum and resolution (<81.3 nHz (S), <55.6 nHz (N)) in lowermost frequencies (<2 {\mu}Hz in most modes). AR is accompanied by an also sharp symmetrical antiresonance P(-) whose both N/S tailing harmonics P(-17) are the well-known PR=~154-day Rieger period dominating planetary dynamics and space weather. Unlike a resonating motor restrained from separating its casing, the freely resonating Sun exhausts the wind in an axial shake-off beyond L1 at highly coherent discrete wave modes generated in the Sun, so to understand solar-type stars, only global decadal scales matter. The result verified against remote data and the experiment, so it instantly replaces dynamo with magnetoalternator and advances Standard Stellar Models, improving fundamental understanding of billions of trillions of solar-type stars. Gauss-Vanicek spectral analysis revolutionizes planetary & space sciences by rigorously simulating multiple spacecraft or fleet formations from a single spacecraft and physics by directly computing nonlinear global dynamics (rendering spherical approximation obsolete).

We present Spitzer 10-34 {\mu}m spectroscopic observations of the diffuse gas in the inner region of the star-forming region IC 348 of the Perseus Molecular Cloud. We find evidence for the strongest mid-IR bands of common molecules as H\textsubscript{2}, OH, H\textsubscript{2}O,CO\textsubscript{2} and NH\textsubscript{3} and of several carbonaceous molecules which may play an important role in the production of more complex hydrocarbons: HCN, C\textsubscript{2}H\textsubscript{2}, C\textsubscript{4}H\textsubscript{2}, HC\textsubscript{3}N, HC\textsubscript{5}N, C\textsubscript{2}H\textsubscript{6}, C\textsubscript{6}H\textsubscript{2}, C\textsubscript{6}H\textsubscript{6}. The excitation diagram of H\textsubscript{2} reveals the presence of warm gas (270 +- 30 K) at the observed locations. Assuming this temperature, the derived abundances of CO\textsubscript{2} and NH\textsubscript{3} relative to H\textsubscript{2} are 10\textsuperscript{-8} and 10\textsuperscript{-7}, respectively. From the water lines we obtain an abundance of order 10\textsuperscript{-6} and higher gas temperatures. The abundances derived for HCN and C\textsubscript{2}H\textsubscript{2}, key molecules in the development of prebiotic building blocks, are of order 10\textsuperscript{-7} and 10\textsuperscript{-9}, respectively. More complex molecules such as PAHs and the fullerenes C\textsubscript{60} and C\textsubscript{70} are also present. IC 348 appears to be very rich and diverse in molecular content. The JWST spectroscopic capabilities may provide details on the spatial distribution of all these molecules and extend the present search to more complex hydrocarbons.

We investigate spherical domain walls~(DWs) nucleated via quantum tunneling in multifield inflationary models and curvature perturbations induced by the inhomogeneous distribution of those DWs. We consider the case that the Euclidean action $S_{E}$ of DWs changes with time during inflation so that most of DWs nucleate when $S_{E}$ reaches the minimum value and the radii of DWs are almost the same. When the Hubble horizon scale exceeds the DW radius after inflation, DWs begin to annihilate and release their energy into background radiation. Because of the random nature of the nucleation process, the statistics of DWs is of the Poisson type and the power spectrum of curvature perturbations has a characteristic slope ${\cal P}_{\cal R}(k)\propto k^{3}$. The amplitude of ${\cal P}_{\cal R}(k)$ depends on the tension and abundance of DWs at the annihilation time while the peak mode depends on the mean separation of DWs. We also numerically obtain the energy spectra of scalar-induced gravitational waves from predicted curvature perturbations which are expected to be observed in multiband gravitational-wave detectors.

The origin of the stellar initial mass function (IMF) and its relation with the core mass function (CMF) are actively debated issues with important implications in astrophysics. Recent observations in the W43 molecular complex of top-heavy CMFs, with an excess of high-mass cores compared to the canonical mass distribution, raise questions about our understanding of the star formation processes and their evolution in space and time. We aim to compare populations of protostellar and prestellar cores in three regions imaged in the ALMA-IMF Large Program. We created an homogeneous core catalogue in W43, combining a new core extraction in W43-MM1 with the catalogue of W43-MM2&MM3 presented in a previous work. Our detailed search for protostellar outflows enabled us to identify between 23 and 30 protostellar cores out of 127 cores in W43-MM1 and between 42 and 51 protostellar cores out of 205 cores in W43-MM2&MM3. Cores with neither outflows nor hot core emission are classified as prestellar candidates. We found a similar fraction of cores which are protostellar in the two regions, about 35%. This fraction strongly varies in mass, from 15-20% at low mass, between 0.8 and 3$M_{\odot} $ up to about 80% above 16$M_{\odot}$. Protostellar cores are found to be, on average, more massive and smaller in size than prestellar cores. Our analysis also revealed that the high-mass slope of the prestellar CMF in W43, $\alpha=-1.46_{-0.19}^{+0.12}$, is consistent with the Salpeter slope, and thus the top-heavy form measured for the global CMF, $\alpha=-0.96$, is due to the protostellar core population. Our results could be explained by clump-fed models in which cores grow in mass, especially during the protostellar phase, through inflow from their environment. The difference between the slopes of the prestellar and protostellar CMFs moreover implies that high-mass cores grow more in mass than low-mass cores.

The next generation of galaxy surveys will provide highly accurate measurements of the large-scale structure of the Universe, allowing for more stringent tests of gravity on cosmological scales. Higher order statistics are a valuable tool to study the non-Gaussianities in the matter field and to break degeneracies between modified gravity and other physical or nuisance parameters. However, understanding from first principles the behaviour of these correlations is essential to characterise deviations from General Relativity (GR), and the purpose of this work. This work uses contemporary ideas of Standard Perturbation Theory on biased tracers to characterize the three point correlation function (3PCF) at tree level for Modified Gravity models with a scale-dependent gravitational strength, and applies the theory to two specific models ($f(R)$ and DGP) that are representative for Chameleon and Vainshtein screening mechanisms. Additionally, we use a multipole decomposition, which apart from speeding up the algorithm to extract the signal from data, also helps to visualize and characterize GR deviations.

We use medium resolution JWST/NIRSpec observations from the Cosmic Evolution Early Release Science (CEERS) Survey to place the first constraints on dust attenuation and star formation based on the Paschen lines for a sizable sample of 63 galaxies at redshifts z=1.0-3.1. Our analysis indicates strong correlations between the Balmer decrement, Ha/Hb, and line ratios that include the Paschen lines (i.e., Paa/Hb, Pab/Hb, and the Paschen decrement, Paa/Pab), suggesting that the former is sensitive to the overall dust obscuration towards HII regions in high-redshift galaxies. The line ratios are used to derive the nebular reddening, E(B-V)neb, and star-formation rates (SFRs). There is marginal evidence that SFRs deduced from the Paschen lines may exceed by ~25% those derived from the Balmer lines alone, suggesting the presence of star formation that is optically thick in the Balmer lines, though deeper observations are needed to confirm this result. Using the Paschen-line constraints on bolometric SFRs, we reevaluate the relationship between dust obscuration and UV spectral slope, and find a reddening of the UV continuum that, on average, follows the SMC extinction curve. This analysis highlights the need for deeper spectroscopy of more representative samples to evaluate nebular dust attenuation and bolometric SFRs in high-redshift galaxies, and their relationship to the reddening of the UV continuum.

Amino acids are building-blocks of proteins, basic constituents of all organisms and essential to life on Earth. They are present in carbonaceous chondrite meteorites and comets, but their origin is still unknown. We present Spitzer spectroscopic observations in the star-forming region IC 348 of the Perseus Molecular Cloud showing the possible detections of mid-IR emission lines consistent with the most intense laboratory bands of the three aromatic amino acids, tyrosine, phenylalanine and tryptophan and the aliphatic amino acids isoleucine and glycine. Based on these tentative identifications, preliminary estimates of column densities give values 10-100 times higher for isoleucine and glycine than for the aromatic amino acids as in some meteorites. Potential counterparts of the strongest laboratory bands of each amino acid are also found in the combined spectrum of 32 interstellar locations obtained in diverse unrelated star-forming regions.

One of the frontiers for advancing what is known about dark matter lies in using strong gravitational lenses to characterize the population of the smallest dark matter halos. There is a large volume of information in strong gravitational lens images -- the question we seek to answer is to what extent we can refine this information. To this end, we forecast the detectability of a mixed warm and cold dark matter scenario using the anomalous flux ratio method from strong gravitational lensed images. The halo mass function of the mixed dark matter scenario is suppressed relative to cold dark matter but still predicts numerous low-mass dark matter halos relative to warm dark matter. Since the strong lens signal is a convolution over a range of dark matter halo masses and since the signal is sensitive to the specific configuration of dark matter halos, not just the halo mass function, degeneracies between different forms of suppression in the halo mass function, relative to cold dark matter, can arise. We find that, with a set of lenses with different configurations of the main deflector and hence different sensitivities to different mass ranges of the halo mass function, the different forms of suppression of the halo mass function between the warm dark matter model and the mixed dark matter model can be distinguished with $40$ lenses with Bayesian odds of 29.4:1.

The radio galaxy NGC 1052 casts absorption features of sulfur-bearing molecules, H$_2$S, SO, SO$_2$, and CS toward the radio continuum emission from the core and jets. Using ALMA, we have measured the equivalent widths of SO absorption features in multiple transitions and determined the temperatures of $344 \pm 43$ K and $26 \pm 4$ K in sub-millimeter and millimeter wavelengths, respectively. Since sub-mm and mm continuum represents the core and jets, the high and low temperatures of the absorbers imply warm environment in the molecular torus and cooler downstream flows. The high temperature in the torus is consistent with the presence of 22-GHz H$_2$O maser emission, vibrationally excited HCN and HCO$^+$ absorption lines, and sulfur-bearing molecules in gas phase released from dust. The origin of the sulfur-bearing gas is ascribed to evaporation of icy dust component through jet-torus interaction. Shock heating is the sole plausible mechanism to maintain such high temperature of gas and dust in the torus. Implication of jet-torus interaction also supports collimation of the sub-relativistic jets by gas pressure of the torus.

Kilonovae are generally believed to originate from the ejecta of binary neutron stars (NSs) or black hole and NS mergers. Free neutrons might be retained in the outermost layer of the ejecta to produce a precursor via $\beta$-decay. During the propagation of kilonovae to observers, a small percentage of them might be gravitationally lensed by foreground objects. In this paper, three lens models, i.e., the point-mass model, the singular isothermal sphere (SIS) model, and the Chang-Refsdal model, were taken into consideration to explore the light curves and polarizations of gravitationally lensed kilonovae. We found that if the time delay between two images exceeds the ejecta heating timescale for the lens mass $\sim 10^{10}~M_\odot$ in the SIS model, a tiny bump-like signal will be generated in the light curve, and the total luminosity will be magnified in all cases. The polarization of lensed kilonovae is significantly enhanced in most cases. Future detections of lensed kilonovae will impose constraints on the morphology of the ejecta and aid in the determination of the nature of compact object mergers and the search for strong gravitational lenses.

Cross-matching operation, which is to find corresponding data for the same celestial object or region from multiple catalogues,is indispensable to astronomical data analysis and research. Due to the large amount of astronomical catalogues generated by the ongoing and next-generation large-scale sky surveys, the time complexity of the cross-matching is increasing dramatically. Heterogeneous computing environments provide a theoretical possibility to accelerate the cross-matching, but the performance advantages of heterogeneous computing resources have not been fully utilized. To meet the challenge of cross-matching for substantial increasing amount of astronomical observation data, this paper proposes Heterogeneous-computing-enabled Large Catalogue Cross-matcher (HLC2), a high-performance cross-matching framework based on spherical position deviation on CPU-GPU heterogeneous computing platforms. It supports scalable and flexible cross-matching and can be directly applied to the fusion of large astronomical cataloguesfrom survey missions and astronomical data centres. A performance estimation model is proposed to locate the performance bottlenecks and guide the optimizations. A two-level partitioning strategy is designed to generate an optimized data placement according to the positions of celestial objects to increase throughput. To make HLC2 a more adaptive solution, the architecture-aware task splitting, thread parallelization, and concurrent scheduling strategies are designed and integrated. Moreover, a novel quad-direction strategy is proposed for the boundary problem to effectively balance performance and completeness. We have experimentally evaluated HLC2 using public released catalogue data. Experiments demonstrate that HLC2 scales well on different sizes of catalogues and the cross-matching speed is significantly improved compared to the state-of-the-art cross-matchers.

On April 27, 2020, the soft gamma ray repeater SGR J1935+2154 entered its intense outburst episode again. Insight-HXMT carried out about one month observation of the source. A total number of 75 bursts were detected during this activity episode by Insight-HXMT, and persistent emission data were also accumulated. We report on the spin period search result and the phase distribution of burst start times and burst photon arrival times of the Insight-HXMT high energy detectors and Fermi Gamma-ray Burst Monitor (GBM). We find that the distribution of burst start times is uniform within its spin phase for both Insight-HXMT and Fermi-GBM observations, whereas the phase distribution of burst photons is related to the type of a burst's energy spectrum. The bursts with the same spectrum have different distribution characteristics in the initial and decay episodes for the activity of magnetar SGR J1935+2154.

A combined uniform and long-time series of Ca-K images from the Kodaikanal Observatory (KO), Mount Wilson Observatory (MWO), and Mauna Loa Solar Observatory (MLSO) were used to identify and study the Ca-K small-scale features and their solar cycle variations over a century. The small scale features are classified into three distinct categories: enhanced network (EN), active network (AN), and quiet network (QN). All these features show that their areas vary according to the 11-year solar cycle. The relative amplitude of the Ca-K network variations agree with that of the sunspot cycle. The total area of these small-scale features varies from about 5% during the minimum phase of the solar cycle to about 20% during its maximum phase. Considering the average intensity and the amplitude of their area variations, we find that the total contribution of EN, AN and QN to the irradiance variation of the Sun is about 3%.

The computation of the Love numbers for a spherically symmetric self-gravitating viscoelastic Earth is a classical problem in global geodynamics. Here we revisit the problem of the numerical evaluation of loading and tidal Love numbers in the static limit for an incompressible planetary body, adopting a Laplace inversion scheme based upon the Post-Widder formula as an alternative to the {traditional viscoelastic normal modes method. We also consider, whithin the same framework, complex-valued, frequency-dependent Love numbers that describe the response to a periodic forcing, which are paramount in the study of the tidal deformation of planets. Furthermore, we numerically obtain the time-derivatives of Love numbers, suitable for modeling geodetic signals in response to surface loads variations. A number of examples are shown, in which time and frequency-dependent Love numbers are evaluated for the Earth and planets adopting realistic rheological profiles. The numerical solution scheme is implemented in ALMA${}^3$ (the plAnetary Love nuMbers cAlculator, version 3), an upgraded open-source Fortran 90 program that computes the Love numbers for radially layered planetary bodies with a wide range of rheologies, including transient laws like Andrade or Burgers.

The powerful high-redshift quasar J2102+6015 (at z=4.575) may provide useful information for studying supermassive black hole growth, galaxy evolution and feedback in the early Universe. The source has so far been imaged with very long baseline interferometry (VLBI) at 2/8 GHz (S/X) bands only, showing complex compact structure. Its total radio spectrum peaks at ~6 GHz in the rest frame. There is no sign of Doppler-boosted jet emission, and the separation of the two major features in its east-west oriented structure spanning ~10 milliarcsec does not change significantly on a timescale longer than a decade. However, VLBI astrometric monitoring observations suggest quasi-periodic (~3 yr) variation in its absolute position. J2102+6015 is presumably a young radio source with jets misaligned with respect to the line of sight. Here we briefly report on our new high-resolution imaging observations made with the European VLBI Network (EVN) at 5 and 22 GHz frequencies in 2021 June, and give an overview of what is currently known about this peculiar distant jetted active galactic nucleus.

Long-term observations of synchrotron emission from supernovae (SNe), covering more than a year after the explosion, provide a unique opportunity to study the poorly-understood evolution of massive stars in the final millennium of their lives via changes in the mass-loss rate. Here, we present a result of our long-term monitoring of a peculiar type IIL SN 2018ivc, using the Atacama Large Millimeter/submillimeter Array (ALMA). Following the initial decay, it showed unprecedented rebrightening starting at ~ a year after the explosion. This is one of the rare examples showing such rebrightening in the synchrotron emission, and the first case at millimeter wavelengths. We find it to be in the optically-thin regime unlike the optically-thick centimeter emission. As such, we can robustly reconstruct the distribution of the circumstellar matter (CSM) and thus the mass-loss history in the final ~1,000 years. We find that the progenitor of SN 2018ivc had experienced a very high mass-loss rate >~10^{-3} Msun/yr ~1,500 years before the explosion, which was followed by a moderately high mass-loss rate (>~10^{-4} Msun/yr) up until the explosion. From this behavior, we suggest SN 2018ivc represents an extreme version of a binary evolution toward SNe IIb, which bridges the hydrogen-poor SNe (toward SNe Ib/c, without a hydrogen envelope) and hydrogen-rich SNe (SNe IIP, with a massive envelope).

In this paper, we present a Ly$\alpha$ halo extended to $\sim200$ kpc identified by stacking $\sim 3300$ Ly$\alpha$ emitters at $z=2.2-2.3$. We carry out imaging observations and data reduction with Subaru/Hyper Suprime-Cam (HSC). Our total survey area is $\sim12$ deg$^2$ and imaging depths are $25.5-27.0$ mag. Using the imaging data, we select 1,240 and 2,101 LAE candidates at $z=2.2$ and 2.3, respectively. We carry out spectroscopic observations of our LAE candidates and data reduction with Magellan/IMACS to estimate the contamination rate of our LAE candidates. We find that the contamination rate of our sample is low (8%). We stack our LAE candidates with a median stacking method to identify the Ly$\alpha$ halo at $z=2$. We show that the Ly$\alpha$ halo is extended to $\sim200$ kpc at a surface brightness level of $10^{-20}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. Comparing to previous studies, our Ly$\alpha$ halo is more extended at radii of $\sim25-100$ kpc, which is not likely caused by the contamination in our sample but by different redshifts and fields instead. To investigate how central galaxies affect surrounding LAHs, we divide our LAEs into subsamples based on the Ly$\alpha$ luminosity ($L_{\rm Ly\alpha}$), rest-frame Ly$\alpha$ equivalent width (EW$_0$), and UV magnitude (M$_{\rm uv}$). We stack the subsamples and find that higher $L_{\rm Ly\alpha}$, lower EW$_0$, and brighter M$_{\rm uv}$ cause more extended halos. Our results suggest that more massive LAEs generally have more extended Ly$\alpha$ halos.

We report the large scale structure and clustering analysis of Ly$\alpha$ emitters (LAEs) and Ly$\alpha$ blobs (LABs) at $z=2.2-2.3$. Using 3,341 LAEs, 117 LABs, and 58 bright (Ly$\alpha$ luminosity $L_{\rm Ly\alpha}>10^{43.4}$ erg s$^{-1}$) LABs at $z=2.2-2.3$ selected with Subaru/Hyper Suprime-Cam (HSC), we calculate the LAE overdensity to investigate the large scale structure at $z=2$. We show that 74% LABs and 78% bright LABs locate in overdense regions, which is consistent with the trend found by previous studies that LABs generally locate in overdense regions. We find that one of our 8 fields dubbed J1349 contains $39/117\approx33\%$ of our LABs and $22/58\approx38\%$ of our bright LABs. A unique and overdense $24'\times12'$ ($\approx 40\times20$ comoving Mpc$^2$) region in J1349 has 12 LABs (8 bright LABs). By comparing to SSA22 that is one of the most overdense LAB regions found by previous studies, we show that the J1349 overdense region contains $\geq 2$ times more bright LABs than the SSA22 overdense region. We calculate the angular correlation functions (ACFs) of LAEs and LABs in the unique J1349 field and fit the ACFs with a power-law function to measure the slopes. The slopes of LAEs and LABs are similar, while the bright LABs show a $\approx 2$ times larger slope suggesting that bright LABs are more clustered than faint LABs and LAEs. We show that the amplitudes of ACFs of LABs are higher than LAEs, which suggests that LABs have a $\approx 10$ times larger galaxy bias and field-to-field variance than LAEs. The strong field-to-field variance is consistent with the large differences of LAB numbers in our 8 fields.

Cosmic variance introduces significant uncertainties into galaxy number density properties when surveying the high-z Universe with a small volume, such uncertainties produce the field-to-field variance of galaxy number $\sigma_{g}$ in observational astronomy. This uncertainty significantly affects the Luminosity Functions (LF) measurement of Lya Emitters (LAEs). For most previous Lya LF studies, $\sigma_{g}$ is often adopted from predictions by cosmological simulations, but barely confirmed by observations. Measuring cosmic variance requires a huge sample over a large volume, exceeding the capabilities of most astronomical instruments. In this study, we demonstrate an observational approach for measuring the cosmic variance contribution for $z\approx2.2$ Lya LFs. The LAE candidates are observed using narrowband and broadband of the Subaru/Hyper Suprime-Cam (HSC), with 8 independent fields, making the total survey area $\simeq11.62$deg$^2$ and a comoving volume of $\simeq8.71\times10^6$Mpc$^3$. These eight fields are selected using the project of MAMMOTH. We report a best-fit Schechter function with parameters $\alpha=-1.75$ (fixed), $L_{Ly\alpha}^{*}=5.18_{-0.40}^{+0.43} \times 10^{42}$erg s$^{-1}$ and $\phi_{Lya}^{*}=4.87_{-0.55}^{+0.54}\times10^{-4}$Mpc$^{-3}$ for the overall Lya LFs. After clipping out the regions that can bias the cosmic variance measurements, we calculate $\sigma_{g}$, by sampling LAEs within multiple pointings assigned on the field image. We investigate the relation between $\sigma_{g}$ and survey volume $V$, and fit a simple power law: $\sigma_g=k\times(\frac{V_{\rm eff}}{10^5 {\rm Mpc}^3})^{\beta}$. We find best-fit values of $-1.209_{-0.106}^{+0.106}$ for $\beta$ and $0.986_{-0.100}^{+0.108}$ for k. We compare our measurements with predictions from simulations and find that the cosmic variance of LAEs might be larger than that of general star-forming galaxies.

The quantitative characterization of amplitudes and periods of perturbations in orbital elements can help to our understanding the minor bodies dynamics, especially in case motion in vicinity of resonances. The main goal of this study is to present the new way of the orbital elements approximation by the results of numerical integration. We develop the approximation of orbital elements of the small body by the perturbations with combinational frequencies. We use a criterion of minimum of the standard error in our approximation. In result we give the approximation of orbital elements for the selected members of the Hilda group in the 3:2 mean motion resonance with Jupiter.

HESS J1809$-$193 is one of the unidentified very-high-energy gamma-ray sources in the H.E.S.S. Galactic Plane Survey (HGPS). It is located in a rich environment, with an energetic pulsar and associated X-ray pulsar wind nebula, several supernova remnants, and molecular clouds in the vicinity. Furthermore, HESS J1809$-$193 was recently detected at energies above 56 TeV with HAWC, which makes it a PeVatron candidate, that is, a source capable of accelerating cosmic rays up to PeV energies. We present a new analysis of the TeV gamma-ray emission of HESS J1809$-$193 with H.E.S.S., based on improved analysis techniques. We find that the emission is best described by two components with distinct morphologies and energy spectra. We complement this study with an analysis of Fermi-LAT data in the same region. Finally, taking into account further multi-wavelength data, we interpret our results both in a hadronic and leptonic framework.

We present average flux density measurements of 151 radio pulsars at 1.4 GHz with the Parkes 'Murriyang' radio telescope. We recommend our results be included in the next version of the ATNF pulsar catalogue. The large sample of pulsars together with their wide dispersion measure (DM) range make this data set useful for studying variability of flux density, pulsar spectra, and interstellar medium (ISM). We derive the modulation indices and structure-function from the flux density time series for 95 and 54 pulsars, respectively. We suggest the modulation index also be included in the next version of the pulsar catalogue to manifest the variability of pulsar flux density. The modulation index of flow density and DM are negatively correlated. The refractive scintillation (RISS) timescales or its lower bound for a set of 15 pulsars are derived. They are very different from theoretical expectations, implying the complicated properties of the ISM along different lines of sight. The structure-function for other pulsars is flat. The RISS parameters for some of these pulsars possibly could be derived with different observing strategies in the future.

A model-independent or non-parametric approach for modeling a database has been widely used in cosmology. In these scenarios, the data has been used directly to reconstruct an underlying function. In this work, we introduce a novel semi-model-independent method to do the task. The new approach not only removes some drawbacks of previous methods but also has some remarkable advantages. We combine the well-known Gaussian linear model with a neural network and introduce a procedure for the reconstruction of an arbitrary function. In the scenario, the neural network produces some arbitrary base functions which subsequently are fed to the Gaussian linear model. Given a prior distribution on the free parameters, the Gaussian linear model provides a close form for the posterior distribution as well as the Bayesian evidence. In addition, contrary to other methods, it is straightforward to compute the uncertainty.

Binary black holes emit gravitational waves as they inspiral towards coalescence. Searches for electromagnetic counterparts to these gravitational waves rely on looking for common sources producing both signals. In this paper, we take a different approach: we investigate the impact of radiation zone effects, including retardation effects and gravitational wave propagation onto the circumbinary disk around stellar-mass, spinning black holes, using general relativistic hydrodynamical simulations. Then we used a general relativistic ray-tracing code to extract its X-ray spectrum and lightcurve. This allowed us to show that radiation zone effects leave an imprint onto the disk, leading to quasi-periodic patterns in the X-ray lightcurve. The amplitude of the modulation is weak (<1%) but increases with time and is strongly dependent on the inclination angle.

Coronal mass ejections (CMEs) are energetic expulsions of organized magnetic features from the Sun. The study of CME quasi-periodicity helps establish a possible relationship between CMEs, solar flares, and geomagnetic disturbances. We used the angular width of CMEs as a criterion for classifying the CMEs in the study. Based on 25 years of observational data, we systematically analyzed the quasi-periodic variations corresponding to the CME occurrence rate of different angular widths in the northern and southern hemispheres, using frequency and time-frequency analysis methods. There are various periods for CMEs of different angular widths: 9 months, 1.7 years, and 3.3-4.3 years. Compared with previous studies based on the occurrence rate of CMEs, we obtained the same periods of 1.2(+-0.01) months, 3.1(+-0.04) months, ~6.1(+-0.4) months, 1.2(+-0.1) years, and 2.4(+-0.4) years. We also found additional periods of all CMEs that appear only in one hemisphere or during a specific solar cycle. For example, 7.1(+-0.2) months and 4.1(+-0.2) years in the northern hemisphere, 1(+-0.004) months, 5.9(+-0.2) months, 1(+-0.1) years, 1.4(+-0.1) years, and 2.4(+-0.4) years in the southern hemisphere, 6.1(+-0.4) months in solar cycle 23 (SC23) and 6.1(+-0.4) months, 1.2(+-0.1) years, and 3.7(+-0.2) years in solar cycle 24 (SC24). The analysis shows that quasi-periodic variations of the CMEs are a link among oscillations in coronal magnetic activity, solar flare eruptions, and interplanetary space.

Measurements of the Elongation Rate of the depth of shower maximum above 1 EeV are reviewed. There is evidence, from four independent estimates of this rate, made in the two hemispheres using three different techniques, for a decrease in the Elongation Rate above ~3 EeV, as first discovered by the Pierre Auger Collaboration over 15 year ago. Unless there is a dramatic change in the hadronic physics above this energy, the mean mass of the primary cosmic rays must increase as a function of energy, well into the decade beyond 10 EeV. To estimate the mass, the use of hadronic models is required, the accuracy of which remains uncertain. However, the possibility of a dramatic change in the hadronic physics appears unlikely, and would be inconsistent with data from the Auger Collaboration on the mass composition in the range 3 to 10 EeV, and on the anisotropy of arrival directions above 8 EeV. Both of these conclusions are insensitive to uncertainties in the shower models. Some remarks are made about the belief, still firmly held by some, that the highest-energy cosmic rays are dominantly protons.

Supernova remnants are commonly considered to produce most of the Galactic cosmic rays via diffusive shock acceleration. However, many questions about the physical conditions at shock fronts, such as the magnetic-field morphology close to the particle acceleration sites, remain open. Here we report the detection of a localized polarization signal from some synchrotron X-ray emitting regions of Tycho's supernova remnant made by the Imaging X-ray Polarimetry Explorer. The derived polarization degree of the X-ray synchrotron emission is 9+/-2% averaged over the whole remnant, and 12+/-2% at the rim, higher than the 7-8% polarization value observed in the radio band. In the west region the polarization degree is 23+/-4%. The X-ray polarization degree in Tycho is higher than for Cassiopeia A, suggesting a more ordered magnetic-field or a larger maximum turbulence scale. The measured tangential polarization direction corresponds to a radial magnetic field, and is consistent with that observed in the radio band. These results are compatible with the expectation of turbulence produced by an anisotropic cascade of a radial magnetic-field near the shock, where we derive a magnetic-field amplification factor of 3.4+/-0.3. The fact that this value is significantly smaller than those expected from acceleration models is indicative of highly anisotropic magnetic-field turbulence, or that the emitting electrons either favor regions of lower turbulence, or accumulate close to where the magnetic-field orientation is preferentially radially oriented due to hydrodynamical instabilities.

Characterizing prestellar cores in star-forming regions is an important step towards the validation of theoretical models of star formation. Thanks to their sub-arcsecond resolution, ALMA observations can potentially provide samples of prestellar cores up to distances of a few kpc, where regions of massive star formation can be targeted. However, the extraction of real cores from dust-continuum observations of turbulent star-forming clouds is affected by complex projection effects. In this work, we study the problem of core extraction both in the idealized case of column-density maps and in the more realistic case of synthetic 1.3 mm ALMA observations. The analysis is carried out on 12 regions of high column density from our 250 pc simulation. We find that derived core masses are highly unreliable, with only a weak correlation between the masses of cores selected in the synthetic ALMA maps and those of the corresponding three-dimensional cores. The fraction of real three-dimensional cores detected in the synthetic maps increases monotonically with mass and remains always below 50%. Above 1 M$_{\odot}$, the core mass function derived from the column-density maps is very steep, while the core mass function from the synthetic ALMA maps has a shallower slope closer to that of the real three-dimensional cores. Because of the very large mass uncertainties, proper guidance from realistic simulations is essential if ALMA observations of protoclusters at kpc distances are to be used to test star-formation models.

Radio Frequency Interference (RFI) of impulsive nature is created by sources like sparking on high-power transmission lines due to gap or corona discharge and automobile sparking, and it affects the entire observing frequency bands of low-frequency radio telescopes. Such RFI is a significant problem at the Upgraded Giant Metrewave Radio Telescope (uGMRT). A real-time RFI filtering scheme has been developed and implemented to mitigate the effect on astronomical observations. The scheme works in real-time on pre-correlation data from each antenna and allows the detection of RFI based on median absolute deviation statistics. The samples are identified as RFI based on user-defined thresholds and are replaced by digital noise, a constant or zeros. We review the testing and implementation of this system at the uGMRT. We illustrate the effectiveness of the filtering for continuum, spectral line and time-domain data. The real-time filter is released for regular observations in the bands falling in 250 - 1450 MHz, and recent observing cycles show growing usage. Further, we explain the relevance of the released system to the Square Kilometer Array (SKA) receiver chain and possible ways of implementation to meet the computational requirements.

When analysing data from air-shower arrays, it has become common practice to use the signal at a considerable distance from the shower axis ($r_\text{opt}$) as a surrogate for the size of the shower. This signal, $S(r_\text{opt}$), can then be related to the primary energy in a variety of ways. After a brief review of the reasons behind the introduction of $r_\text{opt}$ laid out in a seminal paper by Hillas in 1969, it will be shown that $r_\text{opt}$, is a more effective tool when detectors are laid out on a triangular grid than when detectors are deployed on a square grid. This result may have implications for explaining the differences between the flux observed by the Auger and Telescope collaborations above 10\,EeV and should be kept in mind when designing new shower arrays.

In the paper we develop multi-class classification of Fermi-LAT gamma-ray sources using machine learning with hierarchical determination of classes. One of the main challenges in the multi-class classification of the Fermi-LAT sources is that the size of some of the classes is relatively small, for example with less than 10 associated sources belonging to a class. In the paper we propose a hierarchical structure for the determination of the classes. This enables us to have control over the size of classes and to compare the performance of the classification for different numbers of classes. In particular, the class probabilities in the two-class case can be computed either directly by the two-class classification or by summing probabilities of children classes in multi-class classification. We find that the classifications with few large classes have comparable performance with classifications with many smaller classes. Thus, on the one hand, the few-class classification can be recovered by summing probabilities of classification with more classes while, on the other hand, the classification with many classes gives a more detailed information about the physical nature of the sources. As a result of this work, we construct three probabilistic catalogs. Two catalogs are based on Gaussian mixture model for the determination of the groups of classes: one with the requirement for the minimal number of sources in a group to be larger than 100, which results in six final groups of physical classes, while the other one has the requirement that the minimal number of sources in a group is larger than 15, which results in nine groups of physical classes. The third catalog is based on a random forest determination of the groups of classes with the requirement that the minimal number of sources in a group is larger than 100, which results in six groups of classes. (abridged)

Aims. Gaia DR3 GSP-Phot and GSP-Spec parameters of known K- and M-type stars with luminosity class I are examined and compared with parameters collected from the literature, to assess their accuracy and their potential for stellar classification of unknown intrinsically bright late-types. Gaia DR3 GSP-Phot and GSP-Spec parameters were generated by the Astrophysical Parameters Inference Software (Apsis). Methods. In the Gaia DR3 catalog, there are about 40,000 entries with Apsis parameters similar to those of known red supergiants, RSGs, good parallaxes, and infrared 2MASS and WISE data. By using parallactic distances, infrared photometry, and variability information, only 203 new entries are found with luminosities and temperatures similar to that of known RSGs and G-band amplitudes smaller than 0.5 mag. Their low-resolution BP/RP spectra are compared with an empirically made spectral library of BP/RP spectra of known bright late-type stars (C-rich, S-type, O-rich asymptotic giant branch stars (AGBs), and RSGs) to obtain their spectral types. Results. Among them, 15 S-type stars are identified by peculiar absorption features due to ZrO and LaO visible in their BP/RP spectra, one S/C star, and nine C-rich stars by their strong CN absorption bands. K- and M-types can be reproduced with an accuracy of two subtypes. 20 new RSGs are confirmed, of which six have bolometric magnitudes brighter than those of the AGB limit: 2MASS J21015501+4517205, 2MASS J16291280-4956384, 2MASS J10192621-5818105, 2MASS J20230860+3651450, 2MASS J17084131-4026595, and 2MASS J16490055-4217328. The flag for C-rich stars of the Gaia DR3 LPV pipeline is erroneously true for some RSGs and a visual inspection of the BP/RP spectra is mandatory.

HD$\,$172555 is a young ($\sim$20$\,$Myr) A7V star surrounded by a 10$\,$au wide debris disk suspected to be replenished partly by collisions between large planetesimals. Small evaporating transiting bodies, exocomets, have also been detected in this system by spectroscopy. After $\beta\,$Pictoris, this is another example of a system possibly witnessing a phase of heavy bombardment of planetesimals. In such system, small bodies trace dynamical evolution processes. We aim at constraining their dust content by using transit photometry. We performed a 2-day-long photometric monitoring of HD$\,$172555 with the CHEOPS space telescope in order to detect shallow transits of exocomets with a typical expected duration of a few hours. The large oscillations in the lightcurve indicate that HD$\,$172555 is a $\delta\,$Scuti pulsating star. Once removing those dominating oscillations, we find a hint for a transient absorption. If fitted with an exocomet transit model, it corresponds to an evaporating body passing near the star at a distance of $6.8\pm1.4\,$R$_\star$ (or $0.05\pm 0.01\,$au) with a radius of 2.5 km. These properties are comparable to those of the exocomets already found in this system using spectroscopy, as well as those found in the $\beta\,$Pic system. The nuclei of solar system's Jupiter family comets, with radii of 2-6$\,$km, are also comparable in size. This is the first evidence for an exocomet photometric transit detection in the young system of HD$\,$172555.

We perform a comprehensive analysis of a population of 19 X-ray bright tidal disruption events (TDEs), fitting their X-ray spectra with a new, physically self consistent, relativistic accretion disc model. Not all of the TDEs inhabit regions of parameter space where the model is valid, or have sufficient data for a detailed analysis, and physically interpretable parameters for a sub-sample of 11 TDEs are determined. These sources have thermal (power-law free) X-ray spectra. The radial sizes measured from these spectra lie at values consistent with the inner-most stable circular orbit of black holes with masses given by the $M_{\rm BH}-\sigma$ relationship, and can be used as an independent measurement of $M_{\rm BH}$. The bolometric disc luminosity can also be inferred from X-ray data. All of the TDEs have luminosities which are sub-Eddington ($L_{\rm bol, disc} \lesssim L_{\rm edd}$), and larger than the typical hard-state transitional luminosity of X-ray binary discs ($L_{\rm bol, disc} \gtrsim 0.01 L_{\rm edd}$). The {\it peak} bolometric luminosity is found to be linearly correlated with the $M_{\rm BH}-\sigma$ mass. The TDE X-ray-to-bolometric correction can reach values up to $\sim 100$, and grows exponentially at late times, resolving the missing energy problem. We show that the peak disc luminosities of some TDEs are smaller than their observed optical luminosities, implying that not all of the early time optical emission can be sourced from reprocessed disc emission. Our results are supportive of the hypothesis that thermal X-ray bright TDEs are in accretion states analogous to the ``soft'' accretion state of X-ray binaries, and that black hole accretion processes are scale (mass) invariant.

We present measurements of [NeIII]{\lambda}3869 emission in z~1 low-mass galaxies taken from the Keck/DEIMOS spectroscopic surveys HALO7D and DEEPWinds. We identify 167 individual galaxies with significant [NeIII] emission lines, including 112 "dwarf" galaxies with log(M_{\star}/M_{\odot}) < 9.5, with 0.3 < z < 1.4. We also measure [NeIII] emission from composite spectra derived from all [OII]{\lambda}{\lambda}3727,3729 line emitters in this range. This provides a unique sample of [NeIII]-emitters in the gap between well-studied emitters at z = 0 and 2 < z < 3. To study evolution in ionization conditions in the ISM over this time, we analyze the log([NeIII]{\lambda}3869/[OII]{\lambda}{\lambda}3727,3729) ratio (Ne3O2) as a function of the stellar mass and of the log([OIII]{\lambda}{\lambda}4959,5007/[OII]{\lambda}{\lambda}3727,3729) ratio (O32). We find that the typical star-forming dwarf galaxy at this redshift, as measured from the composite spectra, shares the Ne3O2-M_{\star} relation with local galaxies, but have higher O32 at given Ne3O2. This finding implies that the ionization and metallicity characteristics of the z~1 dwarf population do not evolve substantially from z~1 to z=0, suggesting that the known evolution in those parameter from z~2 has largely taken place by z~1. Individual [NeIII]-detected galaxies have emission characteristics situated between local and z~2 galaxies, with elevated Ne3O2 and O32 emission potentially explained by variations in stellar and nebular metallicity. We also compare our dwarf sample to similarly low-mass z > 7 galaxies identified in JWST Early Release Observations, finding four HALO7D dwarfs with similar size, metallicity, and star formation properties.

Rossby waves arise due to the conservation of total vorticity in rotating fluids and may govern the large-scale dynamics of stellar interiors. Recent space missions collected huge information about the light curves and activity of many stars, which triggered observations of Rossby waves in stellar surface and interiors. We aim to study the theoretical properties of Rossby waves in stratified interiors of uniformly rotating radiative stars with sub-adiabatic vertical temperature gradient. We use the equatorial beta-plane approximation and linear vertical gradient of temperature to study the linear dynamics of equatorially trapped Rossby and inertia-gravity waves in interiors of radiative stars. The governing equation is solved by the method of separation of variables in the vertical and latitudinal directions. Vertical and latitudinal solutions of the waves are found to be governed by Bessel functions and Hermite polynomials, respectively. Appropriate boundary conditions at stellar surface and poles define analytical dispersion relations for Rossby, Rossby-gravity and inertia-gravity waves. The waves are confined in surface layers of 30-50 $H_0$, where $H_0$ is surface density scale height, and trapped between the latitudes of $\pm 60^0$. Observable frequencies (normalized by the angular frequency of the stellar rotation) of Rossby waves with $m=1$ ($m=2$), where $m$ is the toroidal wavenumber, are in the interval of 0.65-1 (1.4-2), depending on stellar rotation, radius and surface temperature. Rossby-type waves can be systematically observed using light curves of Kepler and TESS stars. Observations and theory then can be used for the sounding of stellar interiors.

The abundance ratios between isomers of a COM observed in the ISM provides valuable information about the chemistry and physics of the gas and eventually, the history of molecular clouds. In this context, the origin of an abundance of c-HCOOH acid of only 6% the isomer c-HCOOH abundance in cold cores, remains unknown. Herein, we explain the presence of c-HCOOH in dark molecular clouds through the destruction and back formation of c-HCOOH and t-HCOOH in a cyclic process that involves HCOOH and highly abundant molecules such as HCO+ and NH3. We use high-level ab initio methods to compute the potential energy profiles for the cyclic destruction/formation routes of c-HCOOH and t-HCOOH. Global rate constants and branching ratios were calculated based on the transition state theory and the master equation formalism under the typical conditions of the ISM. The destruction of HCOOH by reaction with HCO+ in the gas phase leads to three isomers of the cation HC(OH)2+. The most abundant cation can react in a second step with other abundant molecules of the ISM like NH3 to form back c-HCOOH and t-HCOOH. This mechanism explains the formation of c-HCOOH in dark molecular clouds. Considering this mechanism, the fraction of c-HCOOH with respect t-HCOOH is 25.7%. To explain the 6% reported by the observations we propose that further destruction mechanisms of the cations of HCOOH should be taken into account. The sequential acid-base (SAB) mechanism proposed in this work involves fast processes with very abundant molecules in the ISM. Thus, HCOOH very likely suffers our proposed transformations in the conditions of dark molecular clouds. This is a new approach in the framework of the isomerism of organic molecules in the ISM which has the potential to try to explain the ratio between isomers of organic molecules detected in the ISM.

The length of the solar activity cycle fluctuates considerably. The temporal evolution of the corresponding cycle phase, that is, the deviation of the epochs of activity minima or maxima from strict periodicity, provides relevant information concerning the physical mechanism underlying the cyclic magnetic activity. An underlying strictly periodic process (akin to a perfect "clock"), with the observer seeing a superposition of the perfect clock and a small random phase perturbation, leads to long-term phase stability in the observations. Such behavior would be expected if cycles were synchronized by tides caused by orbiting planets or by a hypothetical torsional oscillation in the solar radiative interior. Alternatively, in the absence of such synchronization, phase fluctuations accumulate and a random walk of the phase ensues, which is a typical property of randomly perturbed dynamo models. Based on the sunspot record and the reconstruction of solar cycles from cosmogenic C14, we carried out rigorous statistical tests in order to decipher whether there exists phase synchronization or random walk. Synchronization is rejected at significance levels of between 95% (28 cycles from sunspot data) and beyond 99% (84 cycles reconstructed from C14, while the existence of random walk in the phases is consistent with all data sets. This result strongly supports randomly perturbed dynamo models with little inter-cycle memory.

Solar magnetic activity is expressed via variations of sunspots and active regions varying on different timescales. The most accepted is an 11-year period supposedly induced by the electromagnetic solar dynamo mechanism. There are also some shorter or longer timescales detected: the biennial cycle (2-2.7 years), Gleisberg cycle (80-100 years), and Hallstatt's cycle (2100-2300 years). Recently, using Principal Component Analysis (PCA) of the observed solar background magnetic field (SBMF), another period of 330-380 years, or Grand Solar Cycle (GSC), was derived from the summary curve of two eigenvectors of SBMF. In this paper, a spectral analysis of the averaged sunspot numbers, solar irradiance, and the summary curve of eigenvectors of SBMF was carried out using Morlet wavelet and Fourier transforms. We detect a 10.7-year cycle from the sunspots and modulus summary curve of eigenvectors as well a 22 years cycle and the grand solar cycle of 342-350-years from the summary curve of eigenvectors. The Gleissberg centennial cycle is only detected on the full set of averaged sunspot numbers for 400 years or by adding a quadruple component to the summary curve of eigenvectors. Another period of 2200-2300 years is detected in the Holocene data of solar irradiance measured from the abundance of $^{14}$C isotope. This period was also confirmed with the period of 2100 years derived from a baseline of the summary curve, supposedly, caused by the solar inertial motion (SIM) induced by the gravitation of large planets. The implication of these findings for different deposition of solar radiation into the northern and southern hemispheres of the Earth caused by the combined effects of the solar activity and solar inertial motion on the terrestrial atmosphere are also discussed.

We consider resolved imaging of faint sources with the solar gravitational lens (SGL) while treating the Sun as an extended gravitating body. We use our new diffraction integral that describes how a spherical electromagnetic wave is modified by the static gravitational field of an extended body, represented by series of multipole moments characterizing its interior mass distribution. Dominated by the solar quadrupole moment, these deviations from spherical symmetry significantly perturb the image that is projected by the Sun into its focal region, especially at solar equatorial latitudes. To study the optical properties of the quadrupole SGL, we develop an approximate solution for the point spread function of such an extended lens. We also derive semi-analytical expressions to estimate signal levels from extended targets. With these tools, we study the impact of solar oblateness on imaging with the SGL. Given the small value of the solar quadrupole moment, the majority of the signal photons arriving from an extended target still appear within the image area projected by the monopole lens. However, these photons are scrambled, thus reducing the achievable signal-to-noise ratio during image recovery process (i.e., deconvolution). We also evaluate the spectral sensitivity for high-resolution remote sensing of exoplanets with the extended SGL. We assess the impact on image quality and demonstrate that despite the adverse effects of the quadrupole moment, the SGL remains uniquely capable of delivering high-resolution imaging and spectroscopy of faint, small and distant targets, notably terrestrial exoplanets within ~30--100 parsec from us.

Magma oceans are episodes of large-scale melting of the mantle of terrestrial planets. The energy delivered by the Moon-forming impact induced a deep magma ocean on the young Earth, corresponding to the last episode of core-mantle equilibration. The crystallization of this magma ocean led to the outgassing of volatiles initially present in the Earth's mantle, resulting in the formation of a secondary atmosphere. During outgassing, the magma ocean acts as a chemical buffer for the atmosphere via the oxygen fugacity, set by the equilibrium between ferrous- and ferric-iron oxides in the silicate melts. By tracking the evolution of the oxygen fugacity during magma ocean solidification, we model the evolving composition of a C-O-H atmosphere. We use the atmosphere composition to calculate its thermal structure and radiative flux. This allows us to calculate the lifetime of the terrestrial magma ocean. We find that, upon crystallizing, the magma ocean evolves from a mildly reducing to a highly oxidized redox state, thereby transiting from a CO- and H2-dominated atmosphere to a CO2- and H2O-dominated one. We find the overall duration of the magma ocean crystallization to depend mostly on the bulk H content of the mantle, and to remain below 1.5 millions years for up to 9 Earth's water oceans' worth of H. Our model also suggests that reduced atmospheres emit lower infrared radiation than oxidized ones, despite of the lower greenhouse effect of reduced species, resulting in a longer magma ocean lifetime in the former case. Although developed for a deep magma ocean on Earth, the framework applies to all terrestrial planet and exoplanet magma oceans, depending on their volatile budgets.

Recently, Rozelot & Eren pointed out that the first solar gravitational moment (J2) might exhibit a temporal variation. The suggested explanation is through the temporal variation of the solar rotation with latitude. This issue is deeper developed due to an accurate knowledge of the long-term variations in solar differential rotation regarding solar activity. Here we analyze solar cycles 12-24, investigating the long-term temporal variations in solar differential rotation. It is shown that J2 exhibits a net modulation over the 13 studied cycles of approximately (89.6 +- 0.1) yr, with a peak-to-peak amplitude of approximately 0.1 x 10-7 for a reference value of 2.07 x 10-7). Moreover, J2 exhibits a positive linear trend in the period of minima solar activity (sunspot number up to around 40) and a marked declining trend in the period of maxima (sunspot number above 50). In absolute magnitude, the mean value of J2 is more significant during periods of minimum than in periods of maximum. These findings are based on observational results that are not free of errors and can be refined further by considering torsional oscillations for example. They are comforted by identifying a periodic variation of the J2 term evidenced through the analysis of the perihelion precession of planetary orbits either deduced from ephemerides or computed in the solar equatorial coordinate system instead of the ecliptic coordinate one usually used.

We show that the energy-weighted angular (zenith, azimuth) distribution of extensive air showers (EAS), produced by Ultra High Energy (UHE) cosmic rays at the Pierre Auger Observatory (PAO), has a thrust axis almost transverse to the interplanetary magnetic field (IMF), with a thrust value $Tp \geq 0.64$ ( where 1.0 means a perfect alignment and 0.5 isotropy). This behavior strongly suggests an effect of the IMF on the charged shower particles, producing additional lateral scattering. We discuss the weakening of the Earth's magnetic field during geomagnetic storms (30\% of observational time) when the IMF becomes preponderant, strengthening the alignment.

The LIGO-India project to build and operate an advanced LIGO (aLIGO) gravitational wave (GW) detector in India in collaboration with LIGO-USA was considered and initiated as an Indian national megascience project in 2011. Procedural formalities and site selection efforts progressed since then and the provisional approval for the Indian national project was obtained in 2016, immediately following the first direct detection of gravitational waves with the aLIGO detectors. With KAGRA GW detector in Japan being tuned to be part of the GW detector network, it is now the occasion to assess the progress of LIGO-India project, and evaluate its relevance and scope for gravitational wave science and astronomy. Various key factors like human-power, management, funding, schedule etc., in the implementation of the project are reassessed in the backdrop of the evolution of the global GW detector sensitivity. In what I consider as a realistic estimate, it will take more than a decade, beyond 2032, to commission the detector even with a fraction of the projected design sensitivity. I estimate that the budget for implementation will be more than doubled, to about Rs. 35 billion (> $430 million). The detrimental consequences for the project are discussed, from my personal point of view. However, a revamped action plan with urgency and the right leadership can make LIGO-India a late but significant success for multi-messenger astronomy for several years after 2032, because of its design similitude to the operational aLIGO detectors. For achieving this, it is imperative that the LIGO-India detector is replanned and launched in the post-O5 upgraded A# version, similar to the projected LIGO-USA detectors.

The habitability of planets is closely connected with the stellar activity, mainly the frequency of flares and the distribution of flare energy. Kepler and TESS find many flaring stars are detected via precise time-domain photometric data, and the frequency and energy distribution of stellar flares on different types of stars are studied statistically. However, the completeness and observational bias of detected flare events from different missions (e.g. Kepler and TESS) vary a lot. We use a unified data processing and detection method for flares events based on the light curve from Kepler and TESS. Then we perform injection and recovery tests in the original light curve of each star for each flare event to correct the completeness and energy of flares. Three samples of flaring stars are selected from Kepler and TESS, with rotating periods from 1 to $\sim$ 5 days. Adopting a hot-blackbody assumption, our results show that the cumulative flare frequency distributions (FFDs) of the same stars in Kepler and TESS bands tend to be consistent after correction, revealing a more natural flaring frequency and energy distribution. Our results also extend the low-energy limit in cumulative FFD fitting to $10^{31.5-33}$ erg on different types of stars. For solar-type stars, the average power-law index of cumulative FFD ($\alpha_{\rm cum}$) is $-0.84$, which indicates that low-energy flares contribute less to the total flare energy. With a piecewise correlation between $\alpha_{\rm cum}$ and $T_{\rm eff}$, $\alpha_{\rm cum}$ first rises with $T_{\rm eff}$ from M2 to K1 stars, then slightly decreases for stars hotter than K1.

Microwave Kinetic Inductance Detectors (MKIDs) are a class of superconducting cryogenic detectors that simultaneously exhibit energy resolution, time resolution and spatial resolution. The pixel yield of MKID arrays is usually a critical figure of merit in the characterisation of an MKIDs array. Currently, for MKIDs intended for the detection of optical and near-infrared photons, only the best arrays exhibit a pixel yield as high as 75-80%. The uniformity of the superconducting film used for the fabrication of MKIDs arrays is often regarded as the main limiting factor to the pixel yield of an array. In this paper we will present data on the uniformity of the TiN/Ti/TiN multilayers deposited at the Tyndall National Institute and compare these results with a statistical model that evaluates how inhomogeneities affect the pixel yield of an array.

We present oxygen abundances relative to hydrogen (O/H) in the narrow line regions (NLRs) gas phases of Seyferts 1 (Sy 1s) and Seyferts 2 (Sy 2s) Active Galactic Nuclei (AGNs). We used fluxes of the optical narrow emission line intensities [$3\,500<\lambda($\AA$)<7\,000$] of 561 Seyfert nuclei in the local universe ($z\lesssim0.31$) from the second catalog and data release (DR2) of the BAT AGN Spectroscopic Survey, which focuses on the \textit{Swift}-BAT hard X-ray ($\gtrsim10$ keV) detected AGNs. We derived O/H from relative intensities of the emission lines via the strong-line methods. We find that the AGN O/H abundances are related to their hosts stellar masses and that they follow a downward redshift evolution. The derived O/H together with the hard X-ray luminosity ($L_{\rm X}$) were used to study the X-ray luminosity-metallicity ($L_{\rm X}$-$Z_{\rm NLR}$) relation for the first time in Seyfert galaxies. In contrast to the broad-line focused ($L_{\rm X}$-$Z_{\rm BLR}$) studies, we find that the $L_{\rm X}$-$Z_{\rm NLR}$ exhibit significant anti-correlations with the Eddington ratio ($\lambda_{\rm Edd}$) and these correlations vary with redshifts. This result indicates that the low-luminous AGNs are more actively undergoing Interstellar Medium (ISM) enrichment through star formation in comparison with the more luminous X-ray sources. Our results suggest that the AGN is somehow driving the galaxy chemical enrichment, as a result of the inflow of pristine gas that is diluting the metal rich gas, together with a recent cessation on the circumnuclear star-formation.

Summed many years of work at Pulkovo, the orbits of 67 wide pairs of visual double and multiple stars (included in 64 systems) which were obtained by the Apparent Motion Parameters (AMP) method are presented. This short arc orbit determination method is based on the most reliable astrometric and astrophysical data corresponding to one instant of time. The rest of the observations accumulated in the world serve to control the quality of the orbit and refine some parameters. All early determined AMP-orbits were compared with new observations, some of them were recalculated, new ones were added. For the stars of Pulkovo program of observations with a 26-inch refractor, the Gaia DR2 data were analised. Based on these data, the orbits of 16 stars were calculated. In 20 cases from 67, the quasi-instant motion according to the Gaia DR2 data at the instant 2015.5 contradicts the motion according to all-world observations. A possible reason is the presence of inner subsystems. The orientation of the obtained orbits in the galactic coordinate system is also given.

We report what is believed to be the first example of fully continuous, 24-hour vertical monitoring of atmospheric optical turbulence. This is achieved using a novel instrument, the 24-hour Shack-Hartmann Image Motion Monitor (24hSHIMM). Optical turbulence is a fundamental limitation for applications such as free-space optical communications, where it limits the achievable bandwidth, and ground-based optical astronomy, restricting the observational precision. Knowledge of the turbulence enables us to select the best sites, design optical instrumentation and optimise the operation of ground-based optical systems. The 24hSHIMM estimates the vertical optical turbulence coherence length, time, angle and Rytov variance from the measurement of a four-layer vertical turbulence profile and a wind speed profile retrieved from meteorological forecasts. To illustrate our advance we show the values of these parameters recorded during a 35-hour, continuous demonstration of the instrument. Due to its portability and ability to work in stronger turbulence, the 24hSHIMM can also operate in urban locations, providing the field with a truly continuous, versatile turbulence monitor for all but the most demanding of applications.

We investigate whether there is a variation in the orbital period of the short-period brown dwarf-mass KELT-1\,b, which is one of the best candidates to observe orbital decay. We obtain 19 high-precision transit light curves of the target using six different telescopes. We add all precise and complete transit light curves from open databases and the literature, as well as the available TESS observations from sectors 17 and 57, to form a transit timing variation (TTV) diagram spanning more than 10 years of observations. The analysis of the TTV diagram, however, is inconclusive in terms of a secular or periodic variation, hinting that the system might have synchronized. We update the transit ephemeris and determine an informative lower limit for the reduced tidal quality parameter of its host star of Q$_{\star}^{\prime} > (8.5 \pm 3.9) \times 10^{6}$ assuming that the stellar rotation is not yet synchronised. Using our new photometric observations, published light curves, the TESS data, archival radial velocities and broadband magnitudes, we also update the measured parameters of the system. Our results are in good agreement with those found in previous analyses.

We present a method of conducting data-driven simulations of solar active regions and flux emergence with the MURaM radiative magnetohydrodynamics (MHD) code. The horizontal electric field derived from the full velocity and magnetic vectors, is implemented at the photospheric (bottom) boundary to drive the induction equation. The energy equation accounts for thermal conduction along magnetic fields, optically-thin radiative loss, and heating of coronal plasma by viscous and resistive dissipation, which allows for a realistic presentation of the thermodynamic properties of coronal plasma that are key to predicting the observational features of solar active regions and eruptions. To validate the method, the photospheric data from a comprehensive radiative MHD simulation of solar eruption (the ground truth) are used to drive a series of numerical experiments. The data-driven simulation reproduces the accumulation of free magnetic energy over the course of flux emergence in the ground truth with an error of 3\%. The onset time is approximately 8\,min delayed compared to the ground truth. However, a precursor-like signature can be identified at the correct onset time. The data-driven simulation captures key eruption-related emission features and plasma dynamics of the ground truth flare over a wide temperature span from $\log_{10}T{=}4.5$ to $\log_{10}T{>}8$. The evolution of the flare and coronal mass ejection as seen in synthetic extreme ultraviolet images is also reproduced with high fidelity. The method helps to understand the evolution of magnetic field in a more realistic coronal environment and to link the magnetic structures to observable diagnostics.

We present a survey of the Perseus molecular cloud in the J $=$ 1$\rightarrow$0 transition of HCN, a widely used tracer of dense molecular gas. The survey was conducted with the CfA 1.2 m telescope, which at 89 GHz has a beam width of 11' and a spectral resolution of 0.85 km s$^{-1}$. A total of 8.1 deg$^2$ was surveyed on a uniform 10' grid to a sensitivity of 14 mK per channel. The survey was compared with similar surveys of CO and dust in order to study and calibrate the HCN line as a dense gas tracer. We find the HCN emission to extend over a considerable fraction of the cloud. We show that the HCN intensity remains linear with H$_2$ column density well into the regime where the CO line saturates. We use radiative transfer modeling to show that this likely results from subthermal excitation of HCN in a cloud where the column and volume densities of H$_2$ are positively correlated. To match our HCN observations the model requires an exponential decrease in HCN abundance with increasing extinction, consistent with HCN depletion onto grains. The modeling also reveals that the mean volume density of H$_2$ in the HCN emitting regions is $\sim$ 10$^4$ cm$^{-3}$, well below the HCN critical density. For the first time, we obtain a direct measurement of the ratio of dense gas mass to HCN luminosity for an entire nearby molecular cloud: $\alpha$(HCN) $=$ 92 M$_\odot$/(K km s$^{-1}$ pc$^2$).

Solaris is a transformative Solar Polar Discovery-class mission concept to address crucial outstanding questions that can only be answered from a polar vantage. Solaris will image the Sun's poles from ~75 degree latitude, providing new insight into the workings of the solar dynamo and the solar cycle, which are at the foundation of our understanding of space weather and space climate. Solaris will also provide enabling observations for improved space weather research, modeling and prediction, revealing a unique, new view of the corona, coronal dynamics and CME eruptions from above.

The eighteenth data release of the Sloan Digital Sky Surveys (SDSS) is the first one for SDSS-V, the fifth generation of the survey. SDSS-V comprises three primary scientific programs, or "Mappers": Milky Way Mapper (MWM), Black Hole Mapper (BHM), and Local Volume Mapper (LVM). This data release contains extensive targeting information for the two multi-object spectroscopy programs (MWM and BHM), including input catalogs and selection functions for their numerous scientific objectives. We describe the production of the targeting databases and their calibration- and scientifically-focused components. DR18 also includes ~25,000 new SDSS spectra and supplemental information for X-ray sources identified by eROSITA in its eFEDS field. We present updates to some of the SDSS software pipelines and preview changes anticipated for DR19. We also describe three value-added catalogs (VACs) based on SDSS-IV data that have been published since DR17, and one VAC based on the SDSS-V data in the eFEDS field.

Milky Way-mass galaxies in the FIRE-2 simulations demonstrate two main modes of star formation. At high redshifts star formation occurs in a series of short and intense bursts, while at low redshifts star formation proceeds at a steady rate with a transition from one mode to another at times ranging from 3 to 7 Gyr ago for different galaxies. We analyse how the mode of star formation affects iron and alpha-element abundance. We find that the early bursty regime imprints a measurable pattern in stellar elemental abundances in the form of a "sideways chevron" shape on the [Fe/H] - [O/Fe] plane and the scatter in [O/Fe] at a given stellar age is higher than when a galaxy is in the steady regime. That suggests that the evolution of [O/Fe] scatter with age provides an estimate of the end of the bursty phase. We investigate the feasibility of observing of this effect by adding mock observational errors to a simulated stellar survey and find that the transition between the bursty and steady phase should be detectable in the Milky Way, although larger observational uncertainties make the transition shallower. We apply our method to observations of the Milky Way from the Second APOKASC Catalog and estimate that the transition to steady star formation in the Milky Way happened 7-8 Gyrs ago, earlier than transition times measured in the simulations.

A Faraday rotation measure (RM) catalogue, or RM Grid, is a valuable resource for the study of cosmic magnetism. Using the second data release (DR2) from the LOFAR Two-metre Sky Survey (LoTSS), we have produced a catalogue of 2461 extragalactic high-precision RM values across 5720 deg$^{2}$ of sky (corresponding to a polarized source areal number density of $\sim$0.43 deg$^{-2}$). The linear polarization and RM properties were derived using RM synthesis from the Stokes $Q$ and $U$ channel images at an angular resolution of 20'' across a frequency range of 120 to 168 MHz with a channel bandwidth of 97.6 kHz. The fraction of total intensity sources ($>1$ mJy beam$^{-1}$) found to be polarized was $\sim$0.2%. The median detection threshold was 0.6 mJy beam$^{-1}$ ($8\sigma_{QU}$), with a median RM uncertainty of 0.06 rad m$^{-2}$ (although a systematic uncertainty of up to 0.3 rad m$^{-2}$ is possible, after the ionosphere RM correction). The median degree of polarization of the detected sources is 1.8%, with a range of 0.05% to 31%. Comparisons with cm-wavelength RMs indicate minimal amounts of Faraday complexity in the LoTSS detections, making them ideal sources for RM Grid studies. Host galaxy identifications were obtained for 88% of the sources, along with redshifts for 79% (both photometric and spectroscopic), with the median redshift being 0.6. The focus of the current catalogue was on reliability rather than completeness, and we expect future versions of the LoTSS RM Grid to have a higher areal number density. In addition, 25 pulsars were identified, mainly through their high degrees of linear polarization.

We study the generation of helical magnetic fields during inflation by considering a model which does not suffer from strong coupling or large back-reaction. Electromagnetic conformal invariance is broken only during inflation by coupling the gauge-invariants $F_{\mu\nu}F^{\mu\nu}$ and $F_{\mu\nu}{\tilde{F}}^{\mu\nu}$ to a time-dependent function $I$ with a sharp transition during inflation. The magnetic power spectrum is scale-invariant up to the transition and very blue-shifted after that. The subsequent evolution of the helical magnetic field is subjected to magneto-hydrodynamical processes, resulting in far larger coherence lengths than those occurring after adiabatic decay. Scale-invariant quadratic gravity is a suitable framework to test the model, providing a natural physical interpretation. We show that fully helical magnetic fields are generated with values in agreement with the lower bounds on fields in the Intergalactic Medium derived from blazar observations. This model holds even at large/intermediate energy scales of inflation, contrary to what has been found in previous works.

By studying the normalized cumulative radial distribution (nCRD) of the stars in the central region of a Monte Carlo-simulated globular cluster, we recently defined three parameters able to pinpoint the stage of internal dynamical evolution reached by the system: $A_5$ (i.e., the area subtended by the nCRD within 5$\%$ the half-mass radius, $r_h$), $P_5$ (the value of the nCRD at 5$\%$ $r_h$), and $S_{2.5}$ (the slope of the nCRD at 2.5$\%$ $r_h$). Here we extend the analysis and explore the effects that different fractions (0$\%$, 10$\%$, and 20$\%$) of primordial binaries and stellar-mass black holes (BHs) induce on the dynamical history of the system. As expected, the gradual contraction of the cluster becomes milder and core collapse shallower for increasing binary fraction. Nevertheless, the cluster dynamical evolution is still properly traced by the three parameters. For models with a larger initial retention of stellar mass BHs the evolution depends on the timescale of their subsequent dynamical ejection. An early dynamical ejection of BHs results in a long-term evolution of the three parameters similar to that found in systems with no initial BH retention. Conversely, in the model that retains a large number of BHs for extended time (slow dynamical ejection of BHs), the system is characterized by a less concentrated structure and by the lack of significant temporal evolution of the three parameters. The smaller values of the three parameters found in this case might be used to indirectly infer the possible presence of BHs in the cluster.

We investigate the backreaction of nonlinear perturbations on the global evolution of the Universe within the cosmic screening approach. To this end, we have considered the second-order scalar perturbations. An analytical study of these perturbations followed by a numerical evaluation shows that, first, the corresponding average values have a negligible backreaction effect on the Friedmann equations and, second, the second-order correction to the gravitational potential is much less than the first-order quantity. Consequently, the expansion of perturbations into orders of smallness in the cosmic screening approach is correct.

For about the last 60 years the search for extraterrestrial intelligence has been monitoring the sky for evidence of remotely detectable technological life beyond Earth, with no positive results to date. While the lack of detection can be attributed to the highly incomplete sampling of the search space, technological emissions may be actually rare enough that we are living in a time when none cross the Earth. This possibility has been considered in the past, but not to quantitatively assess its consequences on the galactic population of technoemissions. Here we derive the likelihood of the Earth not being crossed by signals for at least 60 years to infer upper bounds on their rate of emission. We found less than about one to five emissions per century generated from the Milky Way (95 % credible level), implying optimistic waiting times until the next crossing event of no less than 60 to 1,800 years with a 50 % probability. A significant fraction of highly directional signals increases the emission rates upper bounds, but without systematically changing the waiting time. Our results provide a benchmark for assessing the lack of detection and may serve as a basis to form optimal strategies for the search for extraterrestrial intelligence.

On December 31 2022, the Vlaamse Radio- en Televisieomroeporganisatie (VRT), the national public-service broadcaster for the Flemish Community of Belgium, recovered a video recording of a 1964 interview of Georges Lema\^itre. Up until now, that footage was thought to have been lost. This footage represents a unique insight into the views of the physicist often coined as the "father of the Big Bang". The interview was conducted in French and is available online with Flemish subtitles. In an effort to make this treasure broadly available, we provide in this paper some brief context, an English translation of the interview as well as the French transcript for reference.

During binary-black-hole (BBH) mergers, energy and momenta are carried away from the binary system as gravitational radiation. Access to the radiated energy and momenta allows us to accurately predict the properties of the remnant black hole. We develop a python package gw_remnant to efficiently extract the remnant mass, remnant spin, peak luminosity and the final kick imparted on the remnant black hole directly from the gravitational radiation. We then compute the remnant properties of the final black hole in case of non-spinning BBH mergers with mass ratios ranging from q=2.5 to q=1000 using waveforms generated from BHPTNRSur1dq1e4, a recently developed numerical relativity informed surrogate model based on black-hole perturbation theory framework. We validate our results against the remnant properties estimated from numerical relativity (NR) surrogate models in the comparable mass ratio regime and against recently available high-mass ratio RIT NR simulations at q=[15,32,64]. We find that our remnant property estimates match very closely to the estimates obtained from NR surrogate model and the NR data respectively in both the regimes. We then present BHPTNR_Remnant, a surrogate model for the properties of the remnant black hole in BBH mergers with q=2.5 to q=1000, using Gaussian process regression fitting methods. Finally, we comment on the possible implication of remnant information in gravitational waveform modelling. We make both the gw_remnant and BHPTNR_Remnant packages publicly available.

We examine the behavior of charged pions in neutron-rich matter using heavy-baryon chiral perturbation theory. This study is motivated by the prospect that pions, or pion-like, excitations, may be relevant in neutron-rich matter encountered in core-collapse supernovae and neutron star mergers. We find, as previously expected, that the $\pi^-$ mass increases with density and precludes s-wave condensation, and the mass of the $\pi^+$ mode decreases with increasing density. The energy difference between neutrons and protons in neutron-rich matter related to the nuclear symmetry energy alters the power counting. It enhances higher-order contributions to the pion self-energy. Previously unimportant but attractive diagrams are now enhanced. The net effect of this is the appearance of a new collective mode with the quantum numbers of the $\pi^+$, and a pronounced reduction of the $\pi^+$ mass.

Microlensed gravitational waves (GWs) are likely observable by recognizing the signature of interference caused by $\sim\!\mathcal{O}(10\textrm{--}100)~\textrm{ms}$ time delays between multiple lensed signals. However, the shape of the anticipated microlensed GW signals might be confused with the modulation appearing in the waveform of GWs from precessing compact binary mergers. Their morphological similarity may be an obstacle to template-based searches to correctly identifying the origin of observed GWs and it seamlessly raises a fundamental question, \emph{can we discern microlensed GW signals from the signal of precessing compact binary mergers?} We discuss the feasibility of distinguishing those GWs via examining simulated GW signals with and without the presence of noise. We find that it is certainly possible if we compare signal-to-noise ratios (SNRs) computed with templates of different hypotheses for a given target signal. We show that proper parameter estimation for lensed GWs enables us to identify the targets of interest by focusing on a half number of assumptions for the target signal than the SNR-based test.

We consider stability of non-rotating viscous gaseous stars modeled by the Navier-Stokes-Poisson system. Under general assumptions on the equations of states, we proved that the number of unstable modes for the linearized Navier-Stokes-Poisson system equals that of the linearized Euler-Poisson system modeling inviscid gaseous stars. In particular, the turning point principle holds true for non-rotating stars with or without viscosity. That is, the transition of stability only occurs at the extrema of the total mass and the number of unstable modes is determined by the mass-radius curve. For the proof, we establish an infinite dimensional Kelvin-Tait-Chetaev theorem for a class of linear second order PDEs with dissipation. Moreover, we prove that linear stability implies nonlinear asymptotic stability and linear instability implies nonlinear instability for Navier-Stokes-Poisson system.

It is usually thought that the efolds number of inflation must be bounded by its de Sitter entropy, otherwise we will have an information paradox. However, in light of the island rule for computing the entanglement entropy, we show that such a bound might be nonexistent, while the information flux of primordial perturbation modes the observer after inflation is able to detect follows a Page curve. In corresponding eternally inflating spacetime, it seems that our inflating patch must be accompanied with a neighbourly collapsed patch (eventually developing into a black hole) so that its Hawking radiation might be just our primordial perturbations. Accordingly, the perturbation spectrum we observed will present a ``Page-like" suppression at large scale.

We study the evolution of abelian electromagnetic as well as non-abelian gauge fields, in the presence of space-time oscillations. In the non-abelian case, we consider linear approximation, to analyse the time evolution of the field modes. In both abelian and non-abelian, the mode equations, show the presence of the same parametric resonant spatial modes. The large growth of resonant modes induces large fluctuations in physical observables including those that break the $CP-$symmetry. We also evolve small random fluctuations of fields, using numerical simulations in $2+1$ dimensions. These simulations help study non-linear effects $vs$ the gauge coupling, in the non-abelian case. Our results show that there is an increase in energy density with the coupling, at late times. These results suggest that gravitational waves may excite non-abelian gauge fields more efficiently than electromagnetic fields. Also, gravitational waves in the early Universe and from the merger of neutron stars, black holes etc. may enhance $CP-$violation and generate an imbalance in chiral charge distributions, magnetic fields etc.

We propose a method for detecting light dark matter particles via their creation of chiral phonons in standard model matter. We suggest metal-organic frameworks (MOFs) as candidate materials for such detectors, as their structural flexibility yields low-energy chiral phonons with magnetic moments that are large enough to detect using sensitive magnetometers, and their anisotropy leads to directional sensitivity, which mitigates background contamination. To demonstrate our proposal, we calculate the phononic structure of the MOF InF$_3$($4,4'$-bipyridine), and show that it has highly chiral acoustic phonons. Detection of such chiral phonons via their magnetic moments would dramatically lower the excitation energy threshold to the energy of a single phonon. We show that single phonon detection in a MOF would extend detector reach ten or more orders of magnitude below current limits, enabling exploration of a multitude of as-yet-unprobed dark matter candidates.

The dynamical stability of strange-dwarf hybrid stars and strange planets, constituted by strange-quark-matter cores and dilute-nuclear-matter crusts, is revisited by analyzing the fundamental mode eigenfrequencies of the radial oscillation equations with non-trivial boundary conditions for slow and fast conversions characterizing distinct microphysical scales originating at the density-discontinuous interface. Our calculations indicate that in the case of rapid conversions the so-called {\it reaction mode} plays the fundamental role in these non-compact objects and allow their existence in nature. Interestingly, slow conversions display the same stability window as the seminal work of Glendenning-Kettner-Weber. The robustness of our findings is demonstrated for different transition densities and also using an equation of state from perturbative QCD for the ultra-dense core.

A likely source of a gravitational-wave background (GWB) in the frequency band of the Advanced LIGO, Virgo and KAGRA detectors is the superposition of signals from the population of unresolvable stellar-mass binary-black-hole (BBH) mergers throughout the Universe. Since the duration of a BBH merger in band ($\sim\!1~{\rm s}$) is much shorter than the expected separation between neighboring mergers ($\sim\!10^3~{\rm s}$), the observed signal will be "popcorn-like" or intermittent with duty cycles of order $10^{-3}$. However, the standard cross-correlation search for stochastic GWBs currently performed by the LIGO-Virgo-KAGRA collaboration is based on a continuous-Gaussian signal model, which does not take into account the intermittent nature of the background. The latter is better described by a Gaussian mixture-model, which includes a duty cycle parameter that quantifies the degree of intermittence. Building on an earlier paper by Drasco and Flanagan, we propose a stochastic-signal-based search for intermittent GWBs. For such signals, this search performs better than the standard continuous cross-correlation search. We present results of our stochastic-signal-based approach for intermittent GWBs applied to simulated data for some simple models, and compare its performance to the other search methods, both in terms of detection and signal characterization. Additional testing on more realistic simulated data sets, e.g., consisting of astrophysically-motivated BBH merger signals injected into colored detector noise containing noise transients, will be needed before this method can be applied with confidence on real gravitational-wave data.

Strong-field QED (SFQED) processes are central in determining the dynamics of particles and plasmas in extreme electromagnetic fields such as those present in the vicinity of compact astrophysical objects or generated with ultraintense lasers. SFQEDtoolkit is an open source library designed to allow users for a straightforward implementation of SFQED processes in existing particle-in-cell (PIC) and Monte Carlo codes. Through advanced function approximation techniques, high-energy photon emission and electron-positron pair creation probability rates and energy distributions are calculated within the locally-constant-field approximation (LCFA) as well as with more advanced models [Phys. Rev. A 99, 022125 (2019)]. SFQEDtoolkit is designed to provide users with high-performance and high-accuracy, and neat examples showing its usage are provided. In the near future, SFQEDtoolkit will be enriched to model the angular distribution of the generated particles, i.e., beyond the commonly employed collinear emission approximation, as well as to model spin and polarization dependent SFQED processes. Notably, the generality and flexibility of the presented function approximation approach makes it suitable to be employed in other areas of physics, chemistry and computer science.