Papers for Wednesday, May 31 2023

The aperiodic variability ubiquitously observed from accreting black hole X-ray binary systems is generally analysed within the framework of the so-called ``theory of propagating fluctuations''. In this paper we derive the Fourier transforms of the Green's function solutions of the thin disc equations. These solutions suffice to describe all possible solutions through standard convolution techniques. Solutions are found for both Newtonian discs and general relativistic solutions with a vanishing ISCO stress. We use this new relativistic theory to highlight the Kerr black hole spin dependence of a number of observable variability properties of black hole discs. The phase lags, coherence, and power density spectra of Kerr discs are shown to be strong functions of black hole spin. Observations of the aperiodic variability of black hole accretion sources may now, at least in principle, offer a new avenue to directly constrain black hole spins.

An overview of the methods used to estimate the mass of the Galaxy and the results obtained by various authors recently according to modern data is given. In particular, the estimates obtained based on the analysis of the galactic rotation curve, on the kinematics of the Galactic dwarf satellites and globular clusters, on the streams of such dwarf galaxies, on escape speed, as well as on halo stars are considered. Estimates of the Galaxy mass in the form $M (<r)$, $M_{\rm 200}$ and $M_{\rm vir}$ are considered. According to 20 individual estimates, the average value was found $\overline M_{\rm 200}=0.88\times 10^{12}~M_\odot$ with a dispersion of $0.24\times 10^{12}~M_\odot$ and a weighted average error of $0.06\times 10^{12}~M_\odot$. According to 25 individual estimates, $\overline M_{\rm vir}=1.02\times10^{12}~M_\odot$ was obtained with a dispersion of $0.41\times 10^{12}~M_\odot$ and a weighted average error of $0.09\times10^{12}~M_\odot$.

The Kepler space telescope was responsible for the discovery of over 2,700 confirmed exoplanets, more than half of the total number of exoplanets known today. These discoveries took place during both Kepler's primary mission, when it spent 4 years staring at the same part of the sky, and its extended K2 mission, when a mechanical failure forced it to observe different parts of the sky along the ecliptic. At the very end of the mission, when Kepler was exhausting the last of its fuel reserves, it collected a short set of observations known as K2 Campaign 19. So far, no planets have been discovered in this dataset because it only yielded about a week of high-quality data. Here, we report some of the last planet discoveries made by Kepler in the Campaign 19 dataset. We conducted a visual search of the week of high-quality Campaign 19 data and identified three possible planet transits. Each planet candidate was originally identified with only one recorded transit, from which we were able to estimate the planets' radii and estimate the semimajor axes and orbital periods. Analysis of lower-quality data collected after low fuel pressure caused the telescope's pointing precision to suffer revealed additional transits for two of these candidates, allowing us to statistically validate them as genuine exoplanets. We also tentatively confirm the transits of one planet with TESS. These discoveries demonstrate Kepler's exoplanet detection power, even when it was literally running on fumes.

Large-scale structure formation is studied in a kinetic theory approach, extending the standard perfect pressureless fluid description for dark matter by including the velocity dispersion tensor as a dynamical degree of freedom. The evolution of power spectra for density, velocity and velocity dispersion degrees of freedom is investigated in a non-perturbative approximation scheme based on the Dyson$\unicode{x2013}$Schwinger equation. In particular, the generation of vorticity and velocity dispersion is studied and predictions for the corresponding power spectra are made, which qualitatively agree well with results obtained from $N$-body simulations. It is found that velocity dispersion grows strongly due to non-linear effects and at late times its mean value seems to be largely independent of the initial conditions. By taking this into account, a rather realistic picture of non-linear large-scale structure formation can be obtained, albeit the numerical treatment remains challenging, especially for very cold dark matter models.

We present a spatially resolved study of stellar populations in 6 galaxies with stellar masses $M_*\sim10^{10}M_\odot$ at $z\sim3.7$ using 14-filter JWST/NIRCam imaging from the JADES and JEMS surveys. The 6 galaxies are visually selected to have clumpy substructures with distinct colors over rest-frame $3600-4100$\r{A}, including a bright dominant stellar core that is close to their stellar-light centroids. With 23-filter photometry from HST to JWST, we measure the stellar-population properties of individual structural components via SED fitting using Prospector. We find that the central stellar cores are $\gtrsim2$ times more massive than the Toomre mass, indicating they may not form via in-situ fragmentation. The stellar cores have stellar ages of $0.4-0.7$ Gyr that are similar to the timescale of clump inward migration due to dynamical friction, suggesting that they likely instead formed through the coalescence of giant stellar clumps. While they have not yet quenched, the 6 galaxies are below the star-forming main sequence by $0.2-0.7$ dex. Within each galaxy, we find that the specific star formation rate is lower in the central stellar core, and the stellar-mass surface density of the core is already similar to quenched galaxies of the same masses and redshifts. Meanwhile, the stellar ages of the cores are either comparable to or younger than the extended, smooth parts of the galaxies. Our findings are consistent with model predictions of the gas-rich compaction scenario for the buildup of galaxies' central regions at high redshifts. We are likely witnessing the coeval formation of dense central cores, along with the onset of galaxy-wide quenching at $z>3$.

The spatial resolution of astronomical images is limited by atmospheric turbulence and diffraction in the telescope optics, resulting in blurred images. This makes it difficult to accurately measure the brightness of blended objects because the contributions from adjacent objects are mixed in a time-variable manner due to changes in the atmospheric conditions. However, this effect can be corrected by characterizing the Point Spread Function (PSF), which describes how a point source is blurred on a detector. This function can be estimated from the stars in the field of view, which provides a natural sampling of the PSF across the entire field of view. Once the PSF is estimated, it can be removed from the data through the so-called deconvolution process, leading to images of improved spatial resolution. The deconvolution operation is an ill-posed inverse problem due to noise and pixelization of the data. To solve this problem, regularization is necessary to guarantee the robustness of the solution. Regularization can take the form of a sparse prior, meaning that the recovered solution can be represented with only a few basis eigenvectors. STARRED is a Python package developed in the context of the COSMOGRAIL collaboration and applies to a vast variety of astronomical problems. It proposes to use an isotropic wavelet basis, called Starlets, to regularize the solution of the deconvolution problem. This family of wavelets has been shown to be well-suited to represent astronomical objects. STARRED provides two modules to first reconstruct the PSF, and then perform the deconvolution. It is based on two key concepts: i) the image is reconstructed in two separate channels, one for the point sources and one for the extended sources, and ii) the code relies on the deliberate choice of not completely removing the effect of the PSF, but rather bringing the image to a higher resolution.

Radio searches for extraterrestrial intelligence have mainly targeted the discovery of narrowband continuous-wave beacons and artificially dispersed broadband bursts. Periodic pulse trains, in comparison to the above technosignature morphologies, offer an energetically efficient means of interstellar transmission. A rotating beacon at the Galactic Center (GC), in particular, would be highly advantageous for galaxy-wide communications. Here, we present blipss, a CPU-based open-source software that uses a fast folding algorithm (FFA) to uncover channel-wide periodic signals in radio dynamic spectra. Running blipss on 4.5 hours of 4-8 GHz data gathered with the Robert C. Byrd Green Bank Telescope, we searched the central 6' of our Galaxy for kHz-wide signals with periods between 11-100 s and duty cycles ($\delta$) between 10-50%. Our searches, to our knowledge, constitute the first FFA exploration for periodic alien technosignatures. We report a non-detection of channel-wide periodic signals in our data. Thus, we constrain the abundance of 4-8 GHz extraterrestrial transmitters of kHz-wide periodic pulsed signals to fewer than one in about 600,000 stars at the GC above a 7$\sigma$ equivalent isotropic radiated power of $\approx 2 \times 10^{18}$ W at $\delta \simeq 10\%$. From an astrophysics standpoint, blipss, with its utilization of a per-channel FFA, can enable the discovery of signals with exotic radio frequency sweeps departing from the standard cold plasma dispersion law.

Swift has now completed 18 years of mission, during which it discovered thousands of gamma-ray bursts (GRBs) as well as new classes of high-energy transient phenomena. Its first breakthrough result was the localization of short duration GRBs, which enabled for redshift measurements and kilonova searches. Swift, in synergy with the Hubble Space Telescope and a wide array of ground-based telescopes, provided the first tantalizing evidence of a kilonova in the aftermath of a short GRB. In 2017, Swift observations of the gravitational wave event GW170817 captured the early UV photons from the kilonova AT2017gfo, opening a new window into the physics of kilonovae. Since then, Swift has continued to expand the sample of known kilonovae, leading to the surprising discovery of a kilonova in a long duration GRB. This article will discuss recent advances in the study of kilonovae driven by the fundamental contribution of Swift.

While supermassive binary black holes (SMBBHs) inspiral toward merger they may also accrete significant amounts of matter. To study the dynamics of such a system requires simultaneously describing the evolving spacetime and the dynamics of magnetized plasma. Here we present the first relativistic calculation simulating two equal-mass, non-spinning black holes as they inspiral from an initial separation of $20M$ ($G=c=1$) almost to merger, $\simeq 9M$, while accreting gas from a surrounding disk, where $M$ is the total binary mass. We find that the accretion rate $\dot M$ onto the black holes first decreases during this period and then reaches a plateau, dropping by only a factor of $\sim 3$ despite its rapid inspiral. An estimated bolometric light curve follows the same profile. The minidisks through which the accretion reaches the black holes are very non-standard. Reynolds, not Maxwell, stresses dominate, and they oscillate between two distinct structural states. In one part of the cycle, ``sloshing" streams transfer mass from one minidisk to the other through the L1 point at a rate $\sim 0.1\times$ the accretion rate, carrying kinetic energy at a rate that can be as large as the peak minidisk bolometric luminosity. We also discover that the minidisks have time-varying tilts with respect to the orbital plane similar in magnitude to the circumbinary disk's aspect ratio. The unsigned poloidal flux on the black hole event horizon is roughly constant at a dimensionless level $\phi\sim 2-3$, but doubles just before merger; if the black holes had significant spin, this flux could support jets whose power could approach the radiated luminosity. This simulation is the first to employ our multipatch infrastructure \pwmhd, decreasing computational expense per physical time to $\sim 3\%$ of similar runs using conventional single-grid methods.

Gravitational-wave observations have revealed sources whose unusual properties challenge our understanding of compact-binary formation. Inferring the formation processes that are best able to reproduce such events may therefore yield key astrophysical insights. A common approach is to simulate a population of mergers and count the fraction of these synthetic events that lie within a chosen region in the measured parameters of a real event. Though appealing owing to its simplicity, this approach is flawed because it neglects the full posterior information, depends on a heuristically constructed confidence region, and fails for high signal-to-noise detections. We point out that the statistically consistent solution is to compute the posterior odds between two simulated populations, which crucially is a relative measure, and show how to include the effect of observational biases by conditioning on source detectability. Applying the approach to several gravitational-wave events and astrophysical populations, we assess the degree to which we can conclude model preference not just between distinct formation pathways but also between subpopulations within a given pathway.

We investigate the imprints of accretion and dynamical friction on the gravitational-wave signals emitted by binary black holes embedded in a scalar dark matter cloud. As a key feature in this work, we focus on scalar fields with a repulsive self-interaction that balances against the self-gravity of the cloud. To a first approximation, the phase of the gravitational-wave signal receives extra correction terms at $-4$PN and $-5.5$PN orders, relative to the prediction of vacuum general relativity, due to accretion and dynamical friction, respectively. Future observations by LISA and B-DECIGO have the potential to detect these effects for a large range of scalar masses~$m_\mathrm{DM}$ and self-interaction couplings~$\lambda_4$; observations by ET and Advanced~LIGO could also detect these effects, albeit in a more limited region of parameter space. Crucially, we find that even if a dark matter cloud has a bulk density~$\rho_0$ that is too dilute to be detected via the effects of dynamical friction, the imprints of accretion could still be observable because it is controlled by the independent scale $\rho_a = 4 m_{\rm DM}^4 c^3/(3 \lambda_4 \hbar^3)$. In the models we consider, the infalling dark matter increases in density up to this characteristic scale $\rho_a$ near the Schwarzschild radius, which sets the accretion rate and its associated impact on the gravitational~waveform.

TOI 1416 (BD+42 2504, HIP 70705) is a V=10 late G or early K-type dwarf star with transits detected by TESS. Radial velocities verify the presence of the transiting planet TOI-1416 b, with a period of 1.07d, a mass of $3.48 M_{Earth}$ and a radius of $1.62 R_{Earth}$, implying a slightly sub-Earth density of $4.50$ g cm$^{-3}$. The RV data also further indicate a tentative planet c with a period of 27.4 or 29.5 days, whose nature cannot be verified due to strong suspicions about contamination by a signal related to the Moon's synodic period of 29.53 days. The near-USP (Ultra Short Period) planet TOI-1416 b is a typical representative of a short-period and hot ($T_{eq} \approx$ 1570 K) super-Earth like planet. A planet model of an interior of molten magma containing a significant fraction of dissolved water provides a plausible explanation for its composition, and its atmosphere could be suitable for transmission spectroscopy with JWST. The position of TOI-1416 b within the radius-period distribution corroborates that USPs with periods of less than one day do not form any special group of planets. Rather, this implies that USPs belong to a continuous distribution of super-Earth like planets with periods ranging from the shortest known ones up to ~ 30 days, whose period-radius distribution is delimitated against larger radii by the Neptune desert and by the period-radius valley that separates super-Earths from sub-Neptune planets. In the abundance of small-short periodic planets against period, a plateau between periods of 0.6 to 1.4 days has however become notable that is compatible with the low-eccentricity formation channel. For the Neptune desert, its lower limits required a revision due to the increasing population of short period planets and new limits are provided. These limits are also given in terms of the planets' insolation and effective temperatures.

We investigate the dark energy phenomenology in an extended parameter space where we allow the curvature density of our universe as a free-to-vary parameter. The inclusion of the curvature density parameter is motivated from the recently released observational evidences indicating the closed universe model at many standard deviations. Here we assume that the dark energy equation-of-state follows the PADE approximation, a generalized parametrization that may recover a variety of existing dark energy models. Considering three distinct PADE parametrizations, labeled as PADE-I, SPADE-I and PADE-II, we first constrain the cosmological scenarios driven by them using the joint analyses of a series of recently available cosmological probes, namely, Pantheon sample of Supernovae Type Ia, baryon acoustic oscillations, big bang nucleosynthesis, Hubble parameter measurements from cosmic chronometers, cosmic microwave background distance priors from Planck 2018 and then we include the future Gravitational Waves standard sirens (GWSS) data from the Einstein telescope with the combined analyses of these current cosmological probes. We find that the current cosmological probes indicate a very strong evidence of a dynamical dark energy at more than 99\% CL in both PADE-I, and PADE-II, but no significant evidence for the non-flat universe is found in any of these parametrizations. Interestingly, when the future GWSS data from the Einstein telescope are included with the standard cosmological probes an evidence of a non-flat universe is found in all three parametrizations together with a very strong preference of a dynamical dark energy at more than 99\% CL in both PADE-I, and PADE-II. Although from the information criteria analysis, namely, AIC, BIC, DIC, the non-flat $\Lambda$-Cold Dark Matter model remains the best choice, but from DIC, PADE parametrizations are appealing.

The Galileo Project's acoustic monitoring, omni-directional system (AMOS) aids in the detection and characterization of aerial phenomena. It uses a multi-band microphone suite spanning infrasonic to ultrasonic frequencies, providing an independent signal modality for validation and characterization of detected objects. The system utilizes infrasonic, audible, and ultrasonic systems to cover a wide range of sounds produced by both natural and man-made aerial phenomena. Sound signals from aerial objects can be captured given certain conditions, such as when the sound level is above ambient noise and isn't excessively distorted by its transmission path. Findings suggest that audible sources can be detected up to 1 km away, infrasonic sources can be detected over much longer distances, and ultrasonic at shorter ones. Initial data collected from aircraft recordings with spectral analysis will help develop algorithms and software for quick identification of known aircraft. Future work will involve multi-sensor arrays for sound localization, larger data sets analysis, and incorporation of machine learning and AI for detection and identification of more types of phenomena in all frequency bands.

To date, there are little reliable data on the position, velocity and acceleration characteristics of Unidentified Aerial Phenomena (UAP). The dual hardware and software system described in this document provides a means to address this gap. We describe a weatherized multi-camera system which can capture images in the visible, infrared and near infrared wavelengths. We then describe the software we will use to calibrate the cameras and to robustly localize objects-of-interest in three dimensions. We show how object localizations captured over time will be used to compute the velocity and acceleration of airborne objects.

(Abridged) The Galileo Project aims to improve the detection and classification of aerial objects by using SkyWatch, a passive multistatic radar system. This system utilizes commercial FM radio transmitters, providing critical measurements of range, location, and kinematics, which can help distinguish known aerial objects from those with anomalous movements. SkyWatch can estimate the 3D position and velocity of objects at altitudes up to 80km, distances up to 150km, and speeds up to {\textpm}2{\textpm}2km/s ({\textpm}6{\textpm}6Mach). The radar system's effectiveness is currently being tested with several receivers positioned around the Harvard-Smithsonian Center for Astrophysics. Results will inform the design of a compact receiver for a large-scale network, enhancing the detection of Unidentified Aerial Phenomena (UAP).

(Abridged) The Galileo Project aims to investigate Unidentified Aerial Phenomena (UAP) using an integrated software and instrumentation system for multimodal tracking of aerial phenomena and anomaly recognition. This study aims to highlight outlier events within the high-dimensional parameter space of collected measurements. The project employs various instruments including wide-field cameras, narrow-field instruments, passive multistatic radar arrays, radio spectrum analyzers, microphones, and environmental sensors for a comprehensive collection of data. These tools enable us to corroborate and verify true detections while recognizing artifacts. The team is developing data processing pipelines to fuse multi-sensor data, track hypotheses, classify semi-supervised data, and detect outliers.

The gamma-ray binary LS I+61 303 has been widely monitored at different wavelengths since its discovery more than sixty years ago. However, the nature of the compact object and the peculiar behavior of the system are still largely debated. Aimed at investigating the rapid X-ray variability of LS I+61 303, we have analysed all the archival RXTE/PCA data of the source, taken between 1996 and 2011. The timing analysis yields a periodicity of $P\sim 26.6\pm 0.3$ days, which is statistically compatible with several periodicities reported in the literature for LS I+61 303. Using this period, we performed a data phase-resolved analysis to produce a set of phase-bin-averaged energy spectra and power density spectra. These power density spectra are dominated by weak red noise below 0.1 Hz, and show no signal above this frequency. The amplitude of the red noise varies mildly with the phase, and shows a maximum that coincides with a dip of the X-ray flux and a softer photon index. Aside from low-frequency noise, this analysis does not provide any statistically significant periodic or quasi-periodic timing feature in the RXTE/PCA data of LS I+61 303.

Problems with the concordance cosmology $\Lambda$CDM as the cosmological constant problem, coincidence problems and Hubble tension has led to many proposed alternatives, as the $\Lambda(t)$CDM, where the now called $\Lambda$ cosmological term is allowed to vary due to an interaction with pressureless matter. Here, we analyze one class of these proposals, namely, $\Lambda=\alpha'a^{-2}+\beta H^2+\lambda_*$, based on dimensional arguments. Using SNe Ia, cosmic chronometers data plus constraints on $H_0$ from SH0ES and Planck satellite, we constrain the free parameters of this class of models. By using the Planck prior over $H_0$, we conclude that the $\lambda_*$ term can not be discarded by this analysis, thereby disfavouring models only with the time-variable terms. The SH0ES prior over $H_0$ has an weak evidence in this direction. The subclasses of models with $\alpha'=0$ and with $\beta=0$ can not be discarded by this analysis.

We present 1-5um spectroscopy of the young planetary mass companion TWA 27B (2M1207B) performed with NIRSpec on board the James Webb Space Telescope. In these data, the fundamental band of CH_4 is absent and the fundamental band of CO is weak. The nondetection of CH_4 reinforces a previously observed trend of weaker CH_4 with younger ages among L dwarfs, which has been attributed to enhanced non-equilibrium chemistry among young objects. The weakness of CO may reflect an additional atmospheric property that varies with age, such as the temperature gradient or cloud thickness. We are able to reproduce the broad shape of the spectrum with an ATMO cloudless model that has T=1300 K, non-equilibrium chemistry, and a temperature gradient reduction caused by fingering convection. However, the fundamental bands of CH_4 and CO are somewhat stronger in the model. In addition, the model temperature of 1300 K is higher than expected from evolutionary models given the luminosity and age of TWA 27B (T=1200 K). Previous models of young L-type objects suggest that the inclusion of clouds could potentially resolve these issues; it remains to be seen whether cloudy models can provide a good fit to the 1-5um data from NIRSpec. TWA 27B exhibits emission in Paschen transitions and the He I triplet at 1.083um, which are signatures of accretion that provide the first evidence of a circumstellar disk. We have used the NIRSpec data to estimate the bolometric luminosity of TWA 27B (log L/L_sun=-4.466+/-0.014), which implies a mass of 5-6 MJup according to evolutionary models.

We present evidence of hitherto undiscovered global-scale oscillations in the near-surface shear layer of the Sun. These oscillations are seen as large scale variations of radial shear in both the zonal and meridional flows relative to their mean values. The variations cover all or most of a visible hemisphere, and reverse with a timescale on the order of a solar rotation. A large annual variation in the meridional shear anomaly is understandable in terms of the tilt of the rotation axis, but the rapid oscillations of the shear anomalies in both zonal and the meridional directions appear to be modulated in a more complex, not-quite annual way, although the latter are also strongly modulated by the projected rotational axis angle. Small-scale anomalies in the neighborhood of active regions lend support to their solar origin and physical interpretation. These results were obtained by analyzing ring-diagram fits of low-order modes in high-resolution Doppler data from the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory.

Aims. Fast Radio Bursts are bright radio transients whose origin has not yet explained. The search for a multi-wavelength counterpart of those events can put a tight constrain on the emission mechanism and the progenitor source. Methods. We conducted a multi-wavelength observational campaign on FRB 20180916B between October 2020 and August 2021 during eight activity cycles of the source. Observations were led in the radio band by the SRT both at 336 MHz and 1547 MHz and the uGMRT at 400 MHz. Simultaneous observations have been conducted by the optical telescopes Asiago (Galileo and Copernico), CMO SAI MSU, CAHA 2.2m, RTT-150 and TNG, and X/Gamma-ray detectors on board the AGILE, Insight-HXMT, INTEGRAL and Swift satellites. Results. We present the detection of 14 new bursts detected with the SRT at 336 MHz and seven new bursts with the uGMRT from this source. We provide the deepest prompt upper limits in the optical band fro FRB 20180916B to date. In fact, the TNG/SiFAP2 observation simultaneous to a burst detection by uGMRT gives an upper limit E_optical / E_radio < 1.3 x 10^2. Another burst detected by the SRT at 336 MHz was also co-observed by Insight-HMXT. The non-detection in the X-rays yields an upper limit (1-30 keV band) of E_X-ray / E_radio in the range of (0.9-1.3) x 10^7, depending on which model is considered for the X-ray emission.

The presence of supermassive black holes at redshift z > 6 raises some questions about their formation and growth in the early universe. Due to the construction of new telescopes like the ELT to observe and detect SMBHs, it will be useful to derive theoretical estimates for the population and to compare observations and model predictions in the future. In consequence our main goal is to estimate the population of SMBHs using a semi-analytic code known as Galacticus which is a code for the formation and evolution of galaxies where we are about to include different scenarios for SMBHs formation indicating the initial mass of the black hole seed, its formation conditions and recipes for the evolution of the components of the galaxies. We found that the principal mechanism of growing SMBHs is is via galaxy mergers and accretion of matter. For the comparison of our results with observations, we calculate the radius of influence of the black hole to estimate which part of the population could be detected, leading to relations similar to the observed ones.

Three major planets, Venus, Earth, and Mercury formed out of the solar nebula. A fourth planetesimal, Theia, also formed near Earth where it collided in a giant impact, rebounding as the planet Mars. During this impact Earth lost ${\approx}4$\% of its crust and mantle that is now is found on Mars and the Moon. At the antipode of the giant impact, $\approx$60\% of Earth's crust, atmosphere, and a large amount of mantle were ejected into space forming the Moon. The lost crust never reformed and became the Earth's ocean basins. The Theia impact site corresponds to Indian Ocean gravitational anomaly on Earth and the Hellas basin on Mars. The dynamics of the giant impact are consistent with the rotational rates and axial tilts of both Earth and Mars. The giant impact removed sufficient CO$_2$ from Earth's atmosphere to avoid a runaway greenhouse effect, initiated plate tectonics, and gave life time to form near geothermal vents at the continental margins. Mercury formed near Venus where on a close approach it was slingshot into the Sun's convective zone losing 94\% of its mass, much of which remains there today. Black carbon, from CO$_2$ decomposed by the intense heat, is still found on the surface of Mercury. Arriving at 616 km/s, Mercury dramatically altered the Sun's rotational energy, explaining both its anomalously slow rotation rate and axial tilt. These results are quantitatively supported by mass balances, the current locations of the terrestrial planets, and the orientations of their major orbital axes.

In the present article we put up for discussion the idea of there existing several versions, phases, of the vacuum, in the spirit in which we have long worked on this idea, namely the Multiple Point Criticality Principle, which also says that these different vacuum phases have the same energy density. We mention that we indeed predicted the Higgs mass to be 135 plus minus 10 GeV, which when measured turned out to be 125 GeV, using the assumption of this Multiple Point Criticality Principle. We consider the possibility that there is one type of vacuum in the galaxy clusters (the usual vacuum) and another type of vacuum in the voids. The hope that there could indeed be such a low tension S of the domain wall between these two phases, that it would not totally upset cosmology is based on our dark matter model. In this model dark matter consists of pearls or bubbles of a new vacuum phase, with ordinary matter inside it under very high pressure. The order of magnitude of cubic root S of order MeV or 100 MeV could make such domain walls astronomically viable. We successfully estimate the order of magnitude of the variations in the fine structure constant in different places astronomically, but the similar variations in proton mass over electron mass should have been much bigger than seen experimentally in our model. The Universe s surprisingly early galaxies seen by JWST, James Webb telescope, may agree well with our model. Replacing the usual cosmological constant by domain walls in the standard cosmological model would lead to a cubic root of the tension being cubic root of S of order 30 MeV.

Enceladus is a prime target in the search for life in our solar system, having an active plume likely connected to a large liquid water subsurface ocean. Using the sensitive NIRSpec instrument onboard JWST, we searched for organic compounds and characterized the plume's composition and structure. The observations directly sample the fluorescence emissions of H2O and reveal an extraordinarily extensive plume (up to 10,000 km or 40 Enceladus radii) at cryogenic temperatures (25 K) embedded in a large bath of emission originating from Enceladus' torus. Intriguingly, the observed outgassing rate (300 kg/s) is similar to that derived from close-up observations with Cassini 15 years ago, and the torus density is consistent with previous spatially unresolved measurements with Herschel 13 years ago, suggesting that the vigor of gas eruption from Enceladus has been relatively stable over decadal timescales. This level of activity is sufficient to maintain a derived column density of 4.5x1017 m-2 for the embedding equatorial torus, and establishes Enceladus as the prime source of water across the Saturnian system. We performed searches for several non-water gases (CO2, CO, CH4, C2H6, CH3OH), but none were identified in the spectra. On the surface of the trailing hemisphere, we observe strong H2O ice features, including its crystalline form, yet we do not recover CO2, CO nor NH3 ice signatures from these observations. As we prepare to send new spacecraft into the outer solar system, these observations demonstrate the unique ability of JWST in providing critical support to the exploration of distant icy bodies and cryovolcanic plumes.

A comprehensive understanding of molecular clumps is essential for investigating star formation. We present an algorithm for molecular clump detection, called FacetClumps. This algorithm uses a morphological approach to extract signal regions from the original data. The Gaussian Facet model is employed to fit the signal regions, which enhances the resistance to noise and the stability of the algorithm in diverse overlapping areas. The introduction of the extremum determination theorem of multivariate functions offers theoretical guidance for automatically locating clump centers. To guarantee that each clump is continuous, the signal regions are segmented into local regions based on gradient, and then the local regions are clustered into the clump centers based on connectivity and minimum distance to identify the regional information of each clump. Experiments conducted with both simulated and synthetic data demonstrate that FacetClumps exhibits great recall and precision rates, small location error and flux loss, a high consistency between the region of detected clump and that of simulated clump, and is generally stable in various environments. Notably, the recall rate of FacetClumps in the synthetic data, which comprises $^{13}CO$ ($J = 1-0$) emission line of the MWISP within $11.7^{\circ} \leq l \leq 13.4^{\circ}$, $0.22^{\circ} \leq b \leq 1.05^{\circ}$ and 5 km s$^{-1}$ $\leq v \leq$ 35 km s$^{-1}$ and simulated clumps, reaches 90.2\%. Additionally, FacetClumps demonstrates satisfactory performance when applied to observational data.

Interstellar amides have attracted significant attentions as they are potential precursors for a wide variety of organics essential to life. However, our current understanding of their formation in space is heavily based on observations in star-forming regions and hence the chemical networks lack the constraints on their early origin. In this work, unbiased sensitive spectral surveys with IRAM 30m and Yebes 40m telescopes are used to systematically study a number of amides towards a quiescent Galactic Centre molecular cloud, G+0.693-0.027. We report the first detection of acetamide (CH3C(O)NH2) and trans-N-methylformamide (CH3NHCHO) towards this cloud. In addition, with the wider frequency coverage of the survey, we revisited the detection of formamide (NH2CHO) and urea (carbamide; NH2C(O)NH2), which had been reported previously towards G+0.693-0.027. Our results are compared with those present in the literature including recent laboratory experiments and chemical models. We find constant abundance ratios independently of the evolutionary stages, suggesting that amides related chemistry is triggered in early evolutionary stages of molecular cloud and remain unaffected by the warm-up phase during the star formation process. Although a correlation between more complex amides and NH2CHO have been suggested, alternative formation routes involving other precursors such as acetaldehyde (CH3CHO), methyl isocyanate (CH3NCO) and methylamine (CH3NH2) may also contribute to the production of amides. Observations of amides together with these species towards a larger sample of sources can help to constrain the amide chemistry in the interstellar medium.

In the paper, a comparative primary cosmic rays (PCR) comparative analysis by $E_0$ and the spectra of variable stars by periods is carried out in order to establish the causes of irregularities in the spectrum of PCR by $E_0$. The study was performed using the public database of the KASCADE-Grande experiment and GCVS and ZTF variable star catalogues. It has been suggested that the acceleration of PCR to high and super-high energies occurs not only on the shock waves of supernovae, but also in bursts of giants and super-giants. The relationship between the periods of variable stars and the maximum energy $E_0$ of the nuclei of PCRs generated by these types of stars is shown. Irregularities in the PCR spectrum by $E_0$ are associated with the transition from one dominant stars type to another as $E_0$ increases. The knee in the PCR spectrum at $E_0~=~3-5~PeV$ is associated with a decrease in the contribution of SRB variability stars and a further increase in the contribution of Mira variable stars to the PCR flux. The bump in the PCR spectrum with a maximum at $E_0~=~80~PeV$, established in the KASCADE-Grande experiment, is formed by giant stars and super-giants of the Mira and SRC variability.

Binary molecules such as CO, OH, CH, CN, and C$_2$ are often used as abundance indicators in stars. These species are usually assumed to be formed in chemical equilibrium. The time-dependent effects of hydrodynamics can affect the formation and dissociation of these species and may lead to deviations from chemical equilibrium. We aim to model departures from chemical equilibrium in dwarf stellar atmospheres by considering time-dependent chemical kinetics alongside hydrodynamics and radiation transfer. We examine the effects of a decreasing metallicity and an altered C/O ratio on the chemistry when compared to the equilibrium state. We used the radiation-(magneto)hydrodynamics code CO5BOLD, and its own chemical solver to solve for the chemistry of 15 species and 83 reactions. The species were treated as passive tracers and were advected by the velocity field. The steady-state chemistry was also computed to isolate the effects of hydrodynamics. In most of the photospheres in the models we present, the mean deviations are smaller than $0.2$ dex, and they generally appear above $\log{\tau} = -2$. The deviations increase with height because the chemical timescales become longer with decreasing density and temperature. A reduced metallicity similarly results in longer chemical timescales and in a reduction in yield that is proportional to the drop in metallicity; a decrease by a factor $100$ in metallicity loosely corresponds to an increase by factor $100$ in chemical timescales. As both CH and OH are formed along reaction pathways to CO, the C/O ratio means that the more abundant element gives faster timescales to the constituent molecular species. Overall, the carbon enhancement phenomenon seen in very metal-poor stars is not a result of an improper treatment of molecular chemistry for stars up to a metallicity as low as [Fe/H] = $-3.0$.

Recently, the cosmological tensions, $H_0$ and $S_8$ in particular, have inspired modification of both pre and post recombination physics simultaneously. Early dark energy is a promising pre-recombination solution of the $H_0$ tension known to be compatible with CMB. However, the compatibility of early dark energy, as well as general early resolutions, with CMB is no longer obvious if the late Universe is also modified. Aside from cosmological parameters, the main channel through which late Universe physics affects CMB observable is gravitational lensing. Using a new Gaussian Process function sampling method, we obtained the early Universe (CMB) only constraints on the full shape of the lensing potential, without relying on observation data of the late Universe. It is found that CMB prefers a lensing potential shape that is $\Lambda$CDM-like in $80<L<400$ but with enhanced amplitude beyond this range. The obtained shape constraints can serve as a CMB-compatibility guideline for both late and early Universe model building that modifies the lensing potential.

Following from our recent work, we present here a detailed analysis of dust and star-formation scaling relations, done on a representative sample of nearby galaxies. H$\alpha$ images are analysed in order to derive the integrated flux / luminosity for each galaxy, used as a more instantenous and accurate star-formation rate (SFR) tracer, and the relevant photometric and structural parameters. Dust and inclination corrected H$\alpha$ luminosities and SFRs are subsequently determined using a method that circumvents the assumption of a dust attenuation curve and the use of the Balmer decrements or other hydrogen recombination lines in order to estimate the dust attenuation, which have been shown to be affected by various biases or being inconsistent between different types of galaxies. We investigate the extent to which dust and inclination effects bias the specific parameters of these relations, the scatter and degree of correlation between the parameters, and which relations are fundamental or are just a consequence of others. Our results are consistent within errors with other similar studies. By comparing the scalelengths of the B band optical and H$\alpha$ (star-forming) discs, we found on average, the distribution of star-formation to be more extended than the stellar continuum emission one (the ratio being 1.10), this difference increasing with stellar mass. Similarly, more massive galaxies have a more compact stellar emission surface density than the star-formation one (average ratio of 0.77) for our sample. The method proposed can be applied in larger scale studies of star-formation and ISM evolution, at low to intermediate redshifts.

Magnetic reconnection is a fundamental mechanism in astrophysics. A common challenge in mimicking this process numerically in particular for the Sun is that the solar electrical resistivity is small compared to the diffusive effects caused by the discrete nature of codes. We aim to study different anomalous resistivity models and their respective effects on simulations related to magnetic reconnection in the Sun. We used the Bifrost code to perform a 2D numerical reconnection experiment in the corona that is driven by converging opposite polarities at the solar surface. This experiment was run with three different commonly used resistivity models: 1) the hyper-diffusion model originally implemented in Bifrost, 2) a resistivity proportional to the current density, and 3) a resistivity proportional to the square of the electron drift velocity. The study was complemented with a 1D experiment of a Harris current sheet with the same resistivity models. The 2D experiment shows that the three resistivity models are capable of producing results in satisfactory agreement with each other in terms of the current sheet length, inflow velocity, and Poynting influx. Even though Petschek-like reconnection occurred with the current density-proportional resistivity while the other two cases mainly followed plasmoid-mediated reconnection, the large-scale evolution of thermodynamical quantities such as temperature and density are quite similar between the three cases. For the 1D experiment, some recalibration of the diffusion parameters is needed to obtain comparable results. Specifically the hyper-diffusion and the drift velocity-dependent resistivity model needed only minor adjustments, while the current density-proportional model needed a rescaling of several orders of magnitude.

The non-linear dynamics of scalar fields coupled to matter and gravity can lead to remarkable density-dependent screening effects. In this short review we present the main classes of screening mechanisms, and discuss their tests in laboratory and astrophysical systems. We particularly focus on reviewing numerical and technical aspects involved in modeling the non-linear dynamics of screening. In this review, we focus on tests using laboratory experiments and astrophysical systems, such as stars, galaxies and dark matter halos.

We solve the full quasilinear kinetic equation governing nonresonant interactions of Alfv\'en waves with relativistic plasmas. This work was motivated by the need to determine the energy available for the synchrotron maser in the context of Fast Radio Bursts (FRBs). This interaction can result in plasma heating and the formation of population inversions necessary for the maser. We find that population inversions containing $\sim 1-10\%$ of the distribution's energy form in the relativistic regime, providing an explanation for the formation of the inversion in the environment expected near FRBs.

Observing the spatial distribution and excitation processes of atomic and molecular gas in the inner regions (< 20 au) of young (< 10 Myr) protoplanetary disks helps us to understand the conditions for the formation and evolution of planetary systems. In the framework of the PENELLOPE and ULLYSES projects, we aim to characterize the atomic and molecular component of protoplanetary disks in a sample of 11 Classical T Tauri Stars (CTTs) of the Orion OB1 and $\sigma$-Orionis associations. We analyzed the flux-calibrated optical-forbidden lines and the fluorescent ultraviolet $\rm H_2$ progressions using spectra acquired with ESPRESSO at VLT, UVES at VLT and HST-COS. Line morphologies were characterized through Gaussian decomposition. We then focused on the properties of the narrow low-velocity (FWHM < 40 $km$ $s^{-1}$ and |$v_p$| < 30 $km$ $s^{-1}$) component (NLVC) of the [OI] 630 nm line, compared with the properties of the UV-$\rm H_2$ lines. We found that the [OI]630 NLVC and the UV-$\rm H_2$ lines are strongly correlated in terms of peak velocities, full width at half maximum, and luminosity. The luminosities of the [OI]630 NLVC and UV-$\rm H_2$ correlate with the accretion luminosity with a similar slope, as well as with the luminosity of the CIV 154.8, 155 nm doublet. We discuss such correlations in the framework of the currently suggested excitation processes for the [OI]630 NLVC. Our results can be interpreted in a scenario in which the [OI]630 NLVC and UV-$\rm H_2$ have a common disk origin with a partially overlapped radial extension. We also suggest that the excitation of the [OI] NLVC is mainly induced by stellar FUV continuum photons more than being of thermal origin. This study demonstrates the potential of contemporaneous wide-band high-resolution spectroscopy in linking different tracers of protoplanetary disks.

Metre wavelength type II solar radio bursts are believed to be the signatures of shock-accelerated electrons in the corona. Studying these bursts can give information about the initial kinematics, dynamics and energetics of CMEs in the absence of white-light observations. In this study, we investigate the occurrence of type II bursts in solar cycles 23 and 24 and their association with coronal mass ejections (CMEs). We also explore the possibility of occurrence of type II bursts in the absence of a CME. We performed statistical analysis of type II bursts that occurred between 200 - 25 MHz in solar cycle 23 and 24 and found the temporal association of these radio bursts with CMEs. We categorised the CMEs based on their linear speed and angular width, and studied the distribution of type II bursts with `fast' ($speed ~\geq 500 km/s$), `slow' ($speed ~< 500 km/s$), `wide' ($width ~\geq 60^o$) and `narrow' ($width ~< 60^o$) CMEs. We explored the type II bursts occurrence dependency with solar cycle phases. Our results suggest that type II bursts dominate at heights $\approx 1.7 - 2.3 \pm 0.3 ~R_{\odot}$ with a clear majority having an onset height around 1.7 $\pm 0.3~R_{\odot}$ assuming the four-fold Newkirk model. The results indicate that most of the type II bursts had a white-light CME counterpart, however there were a few type II which did not have a clear CME association. There were more CMEs in cycle 24 than cycle 24. However, the number of type II radio bursts were less in cycle 24 compared to cycle 23. The onset heights of type IIs and their association with wide CMEs reported in this study indicate that the early CME lateral expansion may play a key role in the generation of these radio bursts.

On long enough timescales, chaotic diffusion has the potential to significantly alter the appearance of a dynamical system. The solar system is no exception: diffusive processes take part in the transportation of small bodies and provide dynamical pathways even for the distant trans-Neptunian objects to reach the inner solar system. In this Letter, we carry out a thorough investigation of the nature of chaotic diffusion. We analyze the temporal evolution of the mean squared displacement of ten thousand ensembles of test particles and quantify in each case the diffusion exponent (enabling the classification between normal, sub-, and super-diffusion), the generalized diffusion coefficient, and a characteristic diffusion timescale, too. This latter quantity is compared with an entropy-based timescale, and the two approaches are studied in light of direct computations as well. Our results are given in the context of two-dimensional maps, thereby facilitating the understanding of the relationship between the typical phase space structures and the properties of chaotic diffusion.

Aims. We aim to characterize the properties of the stellar populations in the central few hundred parsecs of nearby galactic nuclei; specifically their age, mass, and 3D geometry. Methods. We use spatially resolved spectroscopic observations of NGC 1433, NGC 1566, and NGC 1808 obtained with SINFONI to constrain a 3D model composed of a spherically symmetric nuclear star cluster (NSC) and an extended thick stellar disk. We computed UV to mid-infrared single stellar population (UMISSP) spectra to determine the age of the stellar populations and construct synthetic observations for our model. To overcome degeneracies between key parameters, we simultaneously fit the spatially resolved line-of-sight velocity, line-of-sight-velocity-dispersion, low-spectral-resolution NIR continuum, and high-spectral-resolution CO absorption features for each pixel. Results. For the three objects, we derive the age and mass of the young and old stellar populations in the NSC and surrounding disk, as well as their 3D geometry: radius for the NSC; thickness, inclination, and position angle for the disk. These results are consistent with published independent measurements when available. Conclusions. The proposed method allows us to derive a consistent 3D model of the stellar populations in nearby galactic centers solely based on a near-infrared IFU observation.

The Lyman-$\alpha$ offers a unique avenue for studying the distribution of matter in the high redshift universe and extracting precise constraints on the nature of dark matter, neutrino masses, and other $\Lambda$CDM extensions. However, interpreting this observable requires accurate modelling of the thermal and ionisation state of the intergalactic medium, and therefore resorting to computationally expensive hydrodynamical simulations. In this work, we build a neural network that serves as a surrogate model for rapid predictions of the one-dimensional \lya flux power spectrum ($P_{\rm 1D}$), thereby making Bayesian inference feasible for this observable. Our emulation technique is based on modelling $P_{\rm 1D}$ as a function of the slope and amplitude of the linear matter power spectrum rather than as a function of cosmological parameters. We show that our emulator achieves sub-percent precision across the full range of scales ($k_{\parallel }=0.1$ to 4Mpc$^{-1}$) and redshifts ($z=2$ to 4.5) considered, and also for three $\Lambda$CDM extensions not included in the training set: massive neutrinos, running of the spectral index, and curvature. Furthermore, we show that it performs at the 1% level for ionisation and thermal histories not present in the training set and performs at the percent level when emulating down to $k_{\parallel}$=8Mpc$^{-1}$. These results affirm the efficacy of our emulation strategy in providing accurate predictions even for cosmologies and reionisation histories that were not explicitly incorporated during the training phase, and we expect it to play a critical role in the cosmological analysis of the DESI survey.

This paper examines options for orbit configurations for a space interferometer. In contrast to previously presented concepts for space very long baseline interferometry, we propose a combination of regular and retrograde near-Earth circular orbits in order to achieve a faster filling of $(u,v)$ coverage. With the rapid relative motion of the telescopes, it will be possible to quickly obtain high quality images of supermassive black holes. As a result of such an approach, it will be possible for the first time to conduct high quality studies of the supermassive black hole close surroundings in dynamics.

Large-scale magnetic fields in the nuclear regions of protogalaxies can trigger the formation and early growth of supermassive black holes (SMBHs) by direct collapse and boosted accretion. Turbulence associated with gravitational infall and star formation can drive the rms field strength toward equipartition with the mean gas kinetic energy; this field has a generic tendency to self-organize into large, coherent structures. If the poloidal component of the field (relative to the rotational axis of a star-forming disc) becomes organized on scales $\lesssim r$ and attains an energy of order a few percent of the turbulent energy in the disc, then dynamo effects are expected to generate magnetic torques capable of boosting the inflow speed and thickening the disk. The accretion flow can transport matter toward the center of mass at a rate adequate to create and grow a massive direct-collapse black hole (DCBH) seed and fuel the subsequent AGN at a high rate, without becoming gravitationally unstable. Fragmentation and star formation are thus suppressed and do not necessarily deplete the mass supply for the accretion flow, in contrast to prevailing models for growing and fueling SMBHs through disc accretion.

It is plausible that most of the Stars in the Milky Way (MW) Galaxy, like the Sun, consist of planetary systems, instead of a single planet. Out of the estimately discovered 3,950 planet-hosting stars, about 860 of them are known to be multiplanetary systems (as of March, 2023). Gravitational microlensing, which is the magnification in the light of a source star, due to a single or several lenses, has proven to be one of the most useful Astrophysical phenomena with many applications. Until now, many extrasolar planets (exoplanets) have been discovered through binary microlensing, where the lens system consists of a star with one planet. In this paper, we discuss and explore the detection of multi-planetary systems that host two exoplanets via microlensing. This is done through the analysis and modeling of possible triple lens configurations (one star and two planets) of a microlensing event. Furthermore, we examine different magnifications and caustic areas of the second planet, by comparing the magnification maps of triple and binary models in different settings. We also discuss the possibility of detecting the corresponding light curves of such planetary systems with the future implementation of the Nancy Grace Roman (Roman) Space Telescope and its Galactic Time Domain survey.

Recent observations from B-star Exoplanet Abundance Study (BEAST) have illustrated the existence of sub-stellar companions around very massive stars. In this paper, we present the detection of two lower mass companions to a relatively nearby ($148.7^{+1.5}_{-1.3}$ pc), young ($17^{+3}_{-4}$ Myr), bright (V=$6.632\pm0.006$ mag), $2.58\pm0.06~ M_{\odot}$ B9V star HIP 81208 residing in the Sco-Cen association, using the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument at the Very Large Telescope (VLT) in Chile. Analysis of the photometry obtained gives mass estimates of $67^{+6}_{-7}~M_J$ for the inner companion and $0.135^{+0.010}_{-0.013}~M_{\odot}$ for the outer companion, indicating the former to be most likely a brown dwarf and the latter to be a low-mass star. The system is compact but unusual, as the orbital planes of the two companions are likely close to orthogonal. The preliminary orbital solutions we derived for the system indicate that the star and the two companions are likely in a Kozai resonance, rendering the system dynamically very interesting for future studies.

We revisit the precessing black hole binary model, a candidate to explain the bizarre quasi-periodic optical flares in OJ-287's light curve, from first principles. We deviate from existing work in three significant ways: 1) Including crucial aspects of relativistic dynamics related to the accretion disc's gravitational moments. 2) Adopting a model-agnostic prescription for the disc's density and scale height. 3) Using monte-carlo Markhov-chain methods to recover reliable system parameters and uncertainties. We showcase our model's predictive power by timing the 2019 great Eddington flare within 40 hr of the observed epoch, exclusively using data available prior to it. Additionally, we obtain a novel direct measurement of OJ-287's disc mass and quadrupole moment exclusively from the optical flare timings. Our improved methodology can uncover previously unstated correlations in the parameter posteriors and patterns in the flare timing uncertainties. In contrast to the established literature, we predict the 26th optical flare to occur on the 21st of August 2023 $\pm$ 32 days, shifted by almost a year with respect to the alleged "missing" flare of October 2022.

Massive elliptical galaxies, that serve as lenses in gravitational lensing time delay measurements of the Hubble parameter $H_0$, often reside in a host group. We consider degeneracies in the modeling of the group halo. When the group effect on imaging can be summarized by its flexion (the next order term beyond shear in the tidal expansion), the posterior likelihood map can develop disjoint local minima, associated with an approximate discrete symmetry of a dominant flexion term. Monte-Carlo Markov Chain (MCMC) algorithms that are not designed to explore a rich posterior landscape can miss some of the minima, introducing systematic bias. We study mock data and demonstrate that the bias in $H_0$ can exceed $10\%$, and pulls the inference value of $H_0$ above its truth value, for a reason that can be traced to the structure of a mismodeled flexion term. MCMC algorithms that are designed to cope with a rich posterior landscape can uncover the structure. If the group is X-ray bright enough, X-ray data may also help to resolve the degeneracy, by pinpointing the group's center of mass. Finally, we show that some implementations in the literature used an inaccurate kinematical prior, mis-modeling the group velocity dispersion by as much as $20\%$

We present FORECAST, a new flexible and adaptable software package that performs forward modeling of the output of any cosmological hydrodynamical simulations to create a wide range of realistic synthetic astronomical images. With customizable options for filters, field of view size and survey parameters, it allows users to tailor the synthetic images to their specific requirements. FORECAST constructs light-cone exploiting the output snapshots of a simulation and computes the observed flux of each simulated stellar element, modeled as a Single Stellar Population, in any chosen set of pass-band filters, including k-correction, IGM absorption and dust attenuation. As a first application, we emulated the GOODS-South field as observed for the CANDELS survey exploiting the IllustrisTNG simulation. We produce images of 200 sq. arcmin., in 13 bands (eight Hubble Space Telescope optical and near-infrared bands from ACS B435 to WFC3 H160, the VLT HAWK-I Ks band, and the four IRAC filters from Spitzer), with depths consistent with the real data. We analysed the images with the same processing pipeline adopted for real data in CANDELS and ASTRODEEP publications, and we compared the results against both the input data used to create the images, and real data, generally finding good agreement with both, with some interesting exceptions which we discuss. As part of this work, we release the FORECAST code and two datasets: the CANDELS dataset analyzed in this study, and 10 JWST CEERS survey-like images (8 NIRCam and 2 MIRI) in a field of view of 200 sq. arcmin. between z=0-20. FORECAST is a flexible tool: it creates images that can then be processed and analysed using standard photometric algorithms, allowing for a consistent comparison among observations and models, and for a direct estimation of the biases introduced by such techniques.

Observations of low-ionization state (LIS) metal lines provide crucial insights into the interstellar medium of galaxies, yet, disentangling the physical processes responsible for the emerging line profiles is difficult. This work investigates how mock spectra generated using a single galaxy in a radiation-hydrodynamical simulation can help us interpret observations of a real galaxy. We create 22,500 C II and Si II spectra from the virtual galaxy at different times and through multiple lines of sight and compare them with the 45 observations of low-redshift star-forming galaxies from the COS Legacy Spectroscopic SurveY (CLASSY). We find that the mock profiles provide accurate replicates to the observations of 38 galaxies with a broad range of stellar masses ($10^6$ to $10^9$ $M_\odot$) and metallicities (0.02 to 0.55 $Z_\odot$). Additionally, we highlight that aperture losses explain the weakness of the fluorescent emission in several CLASSY spectra and must be accounted for when comparing simulations to observations. Overall, we show that the evolution of a single simulated galaxy can produce a large diversity of spectra whose properties are representative of galaxies of comparable or smaller masses. Building upon these results, we explore the origin of the continuum, residual flux, and fluorescent emission in the simulation. We find that these different spectral features all emerge from distinct regions in the galaxy's ISM, and their characteristics can vary as a function of the viewing angle. While these outcomes challenge simplified interpretations of down-the-barrel spectra, our results indicate that high-resolution simulations provide an optimal framework to interpret these observations.

We investigated the giant molecular clouds (GMCs) in M74 (NGC 628) obtained by the PHANGS project. We applied the GMC Types according to the activity of star formation: Type I without star formation, Type II with H$\alpha$ luminosity ($L_{\mathrm{H\alpha}}$) smaller than $10^{37.5} \mathrm{erg~s^{-1}}$, and Type III with $L_{\mathrm{H\alpha}}$ greater than $10^{37.5} \mathrm{erg~s^{-1}}$. In total, 432 GMCs were identified, where the individual GMC Types are 65, 203, and 164, for Type I, Type II, and Type III, respectively. The size and mass of the GMCs range from 23 - 237 pc and $10^{4.9}$ - $10^{7.1}$ M$_{\odot}$, showing a trend that mass and radius increase from Type I to II to III. Clusters younger than 4 Myr and HII regions are found to be concentrated within 150 pc of a GMC, indicating a tight association of these young objects with the GMCs. The virial ratio tends to decrease from Type I to III, indicating that Type III GMCs are most relaxed gravitationally among the three. We interpret that GMCs evolve from Type I to III, as previously found in the LMC. The evolutionary timescales of the three Types are estimated to be 2 Myr, 6 Myr, and 4 Myr, respectively, on a steady state assumption, where we assume the timescale of Type III is equal to the age of the associated clusters, indicating a GMC lifetime of 12 Myr or longer. Chevance et al. (2020) investigated GMCs using the same PHANGS dataset of M74, while these authors did not define a GMC, reaching an evolutionary picture with a 20 Myr duration of the non-star forming phase, five times longer than 4 Myr. We compare the present results with those by Chevance et al. (2020) and argue that defining individual GMCs is essential to understanding GMC evolution.

The tail of the distribution of primordial fluctuations (corresponding to the likelihood of realization of large fluctuations) is of interest, from both theoretical and observational perspectives. In particular, it is relevant for the accurate evaluation of the primordial black hole (PBH) abundance. In this paper, we first analyze the non-perturbative $\delta N$ formalism as a method to non-perturbatively estimate the probability distribution function (PDF) of primordial fluctuations, discuss its underlying assumptions and deal with several subtleties that may arise as a result of considering large fluctuations. Next, we employ the method to study several non-attractor single-field inflationary models as the simplest examples that may lead to the abundant production of PBHs. We conclude that the Gaussian extrapolation from linear perturbation theory may fail drastically to predict the likelihood of large fluctuations. Specifically, we show that a truncation of the tail, a power-law tail, a double-exponential tail, and a doubly peaked distribution can all be realized for the curvature perturbation in the single-field non-attractor models of inflation. We thus show that there is a diverse zoo of possible tails from inflation so that a model-dependent, non-perturbative study of the distribution of the primordial fluctuations seems inevitable concerning PBH abundance.

We show that one-loop corrections to the large-scale power spectrum from small-scale modes in non-slow-roll dynamics are always negligible, namely they are volume suppressed by the ratio of the short to long distance scales. One-loop contributions proportional to the long wavelength tree-level power spectrum, and not sharing this suppression, appear only when considering a subset of vertexes, but they cancel exactly when all relevant interactions are taken into account. We prove the previous statement in two different ways, i.e. by using two equivalent forms of the interaction Hamiltonian. Contributions from boundary terms to equal time correlators are included when necessary.

We analyze the classical linear gravitational effect of idealized pion-like dynamical systems, consisting of light quarks connected by attractive gluonic material with a stress-energy $p=-\rho c^2$ in one or more dimensions. In one orbit of a system of total mass $M$, quarks of mass $m<<M$ expand apart initially with $v/c\sim 1$, slow due to the gluonic attraction, reach a maximum size $R_0 \sim \hbar/ Mc$, then recollapse. We solve the linearized Einstein equations and derive the effect on freely falling bodies for two systems: a gluonic bubble model where uniform gluonic stress-energy fills a spherical volume bounded by a 2D surface comprising the quarks' rest mass, and a gluonic string model where a thin string connects two pointlike quarks. The bubble model is shown to produce a secular mean outward residual velocity of test particles that lie within its orbit. It is shown that the mean gravitational repulsion of bubble-like virtual-pion vacuum fluctuations agrees with the measured value of the cosmological constant, for a bubble with a radius equal to about twice the pion de Broglie length. These results support the view that the gravity of standard QCD vacuum fluctuations is the main source of cosmic acceleration.

We present a new avenue to black hole evaporation using a heat-kernel approach analogous as for the Schwinger effect. Applying this method to an uncharged massless scalar field in a Schwarzschild spacetime, we show that spacetime curvature takes a similar role as the electric field strength in the Schwinger effect. We interpret our results as local pair production in a gravitational field and derive a radial production profile. The resulting emission peaks near the unstable photon orbit. Comparing the particle number and energy flux to the Hawking case, we find both effects to be of similar order. However, our pair production mechanism itself does not explicitly make use of the presence of a black hole event horizon.

X-ray astronomy provides information regarding the electromagnetic emission of active galactic nuclei and X-ray binaries. These events provide details regarding the astrophysical environment of black holes and stars, and help us understand gamma-ray bursts. They produce estimates for the maximum mass of neutron stars and eventually will contribute to the discovery of their equation of state. Thus, it is crucial to study these configurations to increase the yield of X-ray astronomy when combined with multimessenger gravitational-wave astrophysics and black hole shadows. Unfortunately, an exact solution of the field equations does not exist for neutron stars. Nevertheless, there exist a variety of approximate compact objects that may characterize massive or neutron stars. The most studied approximation is the Hartle-Thorne metric that represents slowly-rotating compact objects, like massive stars, white dwarfs and neutron stars. Recent investigations of photon orbits and shadows of such metric revealed that it exhibits chaos close to resonances. Here, we thoroughly investigate particle orbits around the Hartle-Thorne spacetime. We perform an exhaustive analysis of bound motion, by varying all parameters involved in the system. We demonstrate that chaotic regions, known as Birkhoff islands, form around resonances, where the ratio of the radial and polar frequency of geodesics, known as the rotation number, is shared throughout the island. This leads to the formation of plateaus in rotation curves during the most prominent $2/3$ resonance, which designate nonintegrability. We measure their width and show how each parameter affects it. The nonintegrability of Hartle-Thorne metric may affect quasiperiodic oscillations of low-mass X-ray binaries, when chaos is taken into account, and improve estimates of mass, angular momentum and multipole moments of astrophysical compact objects.

The radial direction of the Peccei--Quinn field can drive cosmic inflation, given a non-minimal coupling to gravity. This scenario has been considered to simultaneously explain inflation, the strong $CP$ problem, and dark matter. We argue that Peccei--Quinn inflation is extremely sensitive to higher-dimensional operators. Further combining with the discussion on the axion quality required for solving the strong $CP$ problem, we examine the validity of this scenario. We also show that after Peccei--Quinn inflation, resonant amplifications of the field fluctuation is inevitably triggered.

In an axiverse with numerous axions, the cosmological moduli problem poses a significant challenge because the abundance of axions can easily exceed that of dark matter. The well-established stochastic axion scenario offers a simple solution, relying on relatively low-scale inflation. However, axions are typically subject to mixing due to mass and kinetic terms, which can influence the solution using stochastic dynamics. Focusing on the fact that the QCD axion has a temperature-dependent mass, unlike other axions, we investigate the dynamics of the QCD axion and another axion with mixing. We find that the QCD axion abundance is significantly enhanced and becomes larger than that of the other axion for a certain range of parameters. This enhancement widens the parameter regions accounting for dark matter. In addition, we also find a parameter region in which both axions have enhanced abundances of the same order, which result in multi-component dark matter.

The luminosity distance is a key observable of gravitational-wave observations. We demonstrate how one can correctly retrieve the luminosity distance of compact binary coalescences if the gravitational-wave signal is strongly lensed. We perform a proof-of-concept parameter estimation for the luminosity distance supposing (i) strong lensing produces two lensed gravitational-wave signals, (ii) the advanced LIGO-Virgo network detects both lensed signals as independent events, and (iii) the two events are identified as strongly lensed signals originated from a single compact binary coalescence. Focusing on the maximum magnification allowed in the given lensing scenario, we find that the strong lensing can improve the precision of the distance estimation by up to a factor of two compared to that can be expected for the signal experiencing no lensing. Our results imply that strong lensing of gravitational waves can be helpful for better constraining the distance to the source, and furthermore, the Hubble constant.

In this chapter we review the recent developments of realizing $R^2$-like inflation in the framework of a most general UV nonlocal extension of Einstein's general theory of relativity (GR). It is a well-motivated robust approach towards quantum gravity. In the past decades, nonlocal gravitational theories which are quadratic in curvature have been understood to be ghost-free and super-renormalizable around maximally symmetric spacetimes. However, in the context of early Universe cosmology we show that one must go beyond the quadratic curvature nonlocal gravity in order to achieve a consistent ghost-free framework of Universe evolution from quasi de Sitter to Minkowski spacetime. In this regard, we discuss a construction of a most general nonlocal gravity action that leads to $R^2$-like inflation and discuss the corresponding observational predictions for the scalar and tensor spectral tilts, tensor-to-scalar ratio, and the primordial non-Gaussianities. We present an analysis of how the nonlocal inflationary cosmology goes beyond the established notions of effective field theories of inflation. Finally, we comment on some open questions and prospects of higher curvature nonlocal gravity on its way of achieving the UV completion.

We investigate the inflation and reheating phenomenology in scalar-Einstein-Gauss-Bonnet theory of gravity where a scalar field non-minimally couples with the Gauss-Bonnet (GB) curvature term. Regarding the inflationary phenomenology, we find -- (1) the inflation starts with a quasi de-Sitter phase and has an exit at a finite e-fold, (2) the scalar and tensor perturbations prove to be ghost free and do not suffer from gradient instability, (3) the curvature perturbation amplitude as well as its tilt and the tensor-to-scalar ratio turn out to be simultaneously compatible with the recent Planck data for suitable values of the parameters. After the inflation ends, the scalar field starts to decay to radiation with a constant decay width. For our considered scalar potential and the GB coupling function, the model results to an analytic power law solution of the Hubble parameter and a logarithmic solution of the scalar field during the reheating era, where the exponent of the Hubble parameter determines the effective EoS parameter ($w_\mathrm{eff}$) during the same. The stability of such reheating dynamics is examined by dynamical analysis which ensures that $w_\mathrm{eff}$ can go beyond unity and reach up-to the maximum value of $\mathrm{max}(w_\mathrm{eff}) = 1.56$. The scenario with $w_\mathrm{eff} > 1$ proves to be purely due to the presence of the GB coupling function, which in turn may have important consequences on enhancing the primordial gravitational waves' amplitude observed today. The inflationary e-fold number gets further constrained by the input of the reheating stage. We finally construct the complete forms of scalar potential ($V(\phi)$) and the GB coupling ($\xi(\phi)$) function that smoothly transits from inflation to reheating, and numerically solve the Hubble parameter and the scalar field for such complete forms of $V(\phi)$ and $\xi(\phi)$.

Bayesian inference remains one of the most important tool-kits for any scientist, but increasingly expensive likelihood functions are required for ever-more complex experiments, raising the cost of generating a Monte Carlo sample of the posterior. Recent attention has been directed towards the use of emulators of the posterior based on Gaussian Process (GP) regression combined with active sampling to achieve comparable precision with far fewer costly likelihood evaluations. Key to this approach is the batched acquisition of proposals, so that the true posterior can be evaluated in parallel. This is usually achieved via sequential maximization of the highly multimodal acquisition function. Unfortunately, this approach parallelizes poorly and is prone to getting stuck in local maxima. Our approach addresses this issue by generating nearly-optimal batches of candidates using an almost-embarrassingly parallel Nested Sampler on the mean prediction of the GP. The resulting nearly-sorted Monte Carlo sample is used to generate a batch of candidates ranked according to their sequentially conditioned acquisition function values at little cost. The final sample can also be used for inferring marginal quantities. Our proposed implementation (NORA) demonstrates comparable accuracy to sequential conditioned acquisition optimization and efficient parallelization in various synthetic and cosmological inference problems.