Papers for Monday, May 22 2023

Machine learning techniques have been previously used to model and predict column densities in the TMC-1 dark molecular cloud. In interstellar sources further along the path of star formation, such as those where a protostar itself has been formed, the chemistry is known to be drastically different from that of largely quiescent dark clouds. To that end, we have tested the ability of various machine learning models to fit the column densities of the molecules detected in source B of the Class 0 protostellar binary IRAS 16293-2422. By including a simple encoding of isotopic composition in our molecular feature vectors, we also examine for the first time how well these models can replicate the isotopic ratios. Finally, we report the predicted column densities of the chemically relevant molecules that may be excellent targets for radioastronomical detection in IRAS 16293-2422B.

The exploration and study of exoplanets remain at the frontier of astronomical research, challenging scientists to continuously innovate and refine methodologies to navigate the vast, complex data these celestial bodies produce. This literature the review aims to illuminate the emerging trends and advancements within this sphere, specifically focusing on the interplay between exoplanet detection, classification, and visualization, and the the increasingly pivotal role of machine learning and computational models. Our journey through this realm of exploration commences with a comprehensive analysis of fifteen meticulously selected, seminal papers in the field. These papers, each representing a distinct facet of exoplanet research, collectively offer a multi-dimensional perspective on the current state of the field. They provide valuable insights into the innovative application of machine learning techniques to overcome the challenges posed by the analysis and interpretation of astronomical data. From the application of Support Vector Machines (SVM) to Deep Learning models, the review encapsulates the broad spectrum of machine learning approaches employed in exoplanet research. The review also seeks to unravel the story woven by the data within these papers, detailing the triumphs and tribulations of the field. It highlights the increasing reliance on diverse datasets, such as Kepler and TESS, and the push for improved accuracy in exoplanet detection and classification models. The narrative concludes with key takeaways and insights, drawing together the threads of research to present a cohesive picture of the direction in which the field is moving. This literature review, therefore, serves not just as an academic exploration, but also as a narrative of scientific discovery and innovation in the quest to understand our cosmic neighborhood.

We present the first simulations of core-collapse supernovae (CCSNe) in axial symmetry (2D) with feedback from fast neutrino flavor conversion (FFC). Our schematic treatment of FFCs assumes instantaneous flavor equilibration under the constraint of lepton-number conservation. Systematically varying the spatial domain where FFCs are assumed to occur, we find that they facilitate SN explosions in low-mass (9-12 solar masses) progenitors that otherwise explode with longer time delays, whereas FFCs weaken the tendency to explode of higher-mass (around 20 solar masses) progenitors.

We use the TNG50 cosmological magnetohydrodynamical simulation of the IllustrisTNG project to show that magnetic fields in the circumgalactic medium (CGM) have significant angular structure. This azimuthal anisotropy at fixed distance is driven by galactic feedback processes that launch strong outflows into the halo, preferentially along the minor axes of galaxies. These feedback-driven outflows entrain strong magnetic fields from the interstellar medium, dragging fields originally amplified by small-scale dynamos into the CGM. At the virial radius, $z=0$ galaxies with M$_\star \sim 10^{10}\,\rm{M_\odot}$ show the strongest anisotropy ($\sim 0.35$ dex). This signal weakens with decreasing impact parameter, and is also present but weaker for lower mass as well as higher mass galaxies. Creating mock Faraday rotation measure (RM) sightlines through the simulated volume, we find that the angular RM trend is qualitatively consistent with recent observational measurements. We show that rich structure is present in the circumgalactic magnetic fields of galaxies. However, TNG50 predicts small RM amplitudes in the CGM that make detection difficult as a result of other contributions along the line of sight.

We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback and magnetic fields, to investigate the effectiveness of different fitting methods for volume and column density probability distribution functions (PDFs). The first method fits a piecewise lognormal and power-law (PL) function to recover PDF parameters such as the PL slope and transition density. The second method fits a polynomial spline function and examines the first and second derivatives of the spline to determine the PL slope and the functional transition density. We demonstrate that fitting a spline allows us to directly determine if the data has multiple PL slopes. The first PL (set by the transition between lognormal and PL function) can also be visualized in the derivatives directly. In general, the two methods produce fits that agree reasonably well for volume density but vary for column density, likely due to the increased statistical noise in column density maps as compared to volume density. We test a well-known conversion for estimating volume density PL slopes from column density slopes and find that the spline method produces a better match (\c{hi}2 of 2.38 vs \c{hi}2 of 5.92), albeit with a significant scatter. Ultimately, we recommend the use of both fitting methods on column density data to mitigate the effects of noise.

The flow of gas into and out of galaxies leaves traces in the circumgalactic medium which can then be studied using absorption lines towards background quasars. We analyse 27 log(N_HI) > 18.0 HI absorbers at z = 0.2 to 1.4 from the MUSE-ALMA Halos survey with at least one galaxy counterpart within a line of sight velocity of +/-500 km s^{-1}. We perform 3D kinematic forward modelling of these associated galaxies to examine the flow of dense, neutral gas in the circumgalactic medium. From the VLT/MUSE, HST broadband imaging and VLT/UVES and Keck/HIRES high-resolution UV quasar spectroscopy observations, we compare the impact parameters, star-formation rates and stellar masses of the associated galaxies with the absorber properties. We find marginal evidence for a bimodal distribution in azimuthal angles for strong HI absorbers, similar to previous studies of the MgII and OVI absorption lines. There is no clear metallicity dependence on azimuthal angle and we suggest a larger sample of absorbers are required to fully test the relationship predicted by cosmological hydrodynamical simulations. A case-by-case study of the absorbers reveals that ten per cent of absorbers are consistent with gas accretion, up to 30 per cent trace outflows while the remainder trace gas in the galaxy disk, the intragroup medium and low-mass galaxies below the MUSE detection limit. Our results highlight that the baryon cycle directly affects the dense neutral gas required for star-formation and plays a critical role in galaxy evolution.

$clustertools$ is a Python package for analyzing star cluster simulations. The package is built around the $StarCluster$ class, which stores all data read in from the snapshot of a given model star cluster. The package contains functions for loading data from commonly used $N$-body codes, generic snapshots, and software for generating initial conditions. All operations and functions within $clustertools$ are then designed to act on a $StarCluster$. $clustertools$ can be used for unit and coordinate transformations, the calculation of key structural and kinematic parameters, analysis of the cluster's orbit and tidal tails, and measuring common cluster properties like its mass function, density profile, and velocity dispersion profile (among others). While originally designed with star clusters in mind, $clustertools$ can be used to study other types of $N$-body systems, including stellar streams and dark matter sub-halos.

Extremely Red Quasars (ERQs) are thought to represent a brief episode of young quasar and galactic evolution characterized by rapid outflows and obscured growth due to dusty environments. We use new redshift measurements from CO and Ly$\alpha$ emission-lines to better constrain outflow velocities from previous line measurements. We present sample of 82 ERQs, and the analysis confirms that ERQs have a higher incidence of large CIV blueshifts, accompanied by large Rest Equivalent Widths (REWs) and smaller line widths than blue quasars. We find that strong blueshifts (>2000 km s$^{-1}$) are present in 12/54 (22.22 per cent) of ERQs with the most robust redshift indicators. At least 4 out of 15 ERQs in the sample also have blueshifts in their H$\beta$ and low-ionization UV lines ranging from $-$500 to $-$1500 km s$^{-1}$. ERQs with strong CIV blueshifts are substantially offset in CIV REW and Full-Width at Half-Maximum (FWHM) from typical blue quasars in the same velocity range. ERQs have average values of REW = 124 A and FWHM = 5274 km s$^{-1}$, while blue quasars have REW = 24 A and FWHM = 6973 km s$^{-1}$. The extreme nature of the outflows in ERQs might explain some of their other spectral properties, such as the large CIV REWs and peculiar wingless profiles owing to more extended broad-line regions participating in outflows. The physical reasons for the extreme outflow properties of ERQs are unclear; however, larger Eddington ratios and/or softer ionizing spectra incident on the outflow gas cannot be ruled out.

We use medium- and high-resolution spectroscopy of close pairs of quasars to analyze the circumgalactic medium (CGM) surrounding 32 damped Ly$\alpha$ absorption systems (DLAs). The primary quasar sightline in each pair probes an intervening DLA in the redshift range $1.6<z_\text{abs}<3.5$, such that the secondary sightline probes absorption from Ly$\alpha$ and a large suite of metal-line transitions (including $~\rm OI$, $~\rm CII$, $~\rm CIV$, $~\rm SiII$, and $~\rm SiIV$) in the DLA host galaxy's CGM at transverse distances $24\ \text{kpc}\le R_\bot\le284~\rm kpc$. Analysis of Ly$\alpha$ in the CGM sightlines shows an anti-correlation between $R_\bot$ and $~\rm HI$ column density ($N_\text{HI}$) with 99.8$\%$ confidence, similar to that observed around luminous galaxies. The incidences of $~\rm CII$ and $~\rm SiII$ with $N>10^{13}~\rm cm^{-2}$ within 100 kpc of DLAs are larger by $2\sigma$ than those measured in the CGM of Lyman break galaxies (C$_f(N_\text{CII})>0.89$ and C$_f(N_\text{SiII})=0.75_{-0.17}^{+0.12}$). Metallicity constraints derived from ionic ratios for nine CGM systems with negligible ionization corrections and $N_\text{HI}>10^{18.5}~\rm cm^{-2}$ show a significant degree of scatter (with metallicities/limits across the range $-2.06\lesssim\log Z/Z_{\odot}\lesssim-0.75$), suggesting inhomogeneity in the metal distribution in these environments. Velocity widths of $\text{CIV}\lambda1548$ and low-ionization metal species in the DLA vs. CGM sightlines are strongly ($>2\sigma$) correlated, suggesting they trace the potential well of the host halo over $R_\bot\lesssim300$ kpc scales. At the same time, velocity centroids for $\text{CIV}\lambda1548$ differ in DLA vs. CGM sightlines by $>100~\rm km\ s^{-1}$ for $\sim50\%$ of velocity components, but few components have velocities that would exceed the escape velocity assuming dark matter host halos of $\ge10^{12}M_\odot$.

Prolate rotation is characterized by a significant stellar rotation around a galaxy's major axis, which contrasts with the more common oblate rotation. Prolate rotation is thought to be due to major mergers and thus studies of prolate-rotating systems can help us better understand the hierarchical process of galaxy evolution. Dynamical studies of such galaxies are important to find their gravitational potential profile, total mass, and dark matter fraction. Recently, it has been shown in a cosmological simulation that it is possible to form a prolate-rotating dwarf galaxy following a dwarf-dwarf merger event. The simulation also shows that the unusual prolate rotation can be time enduring. In this particular example, the galaxy continued to rotate around its major axis for at least $7.4$\,Gyr (from the merger event until the end of the simulation). In this project, we use mock observations of the hydro-dynamically simulated prolate-rotating dwarf galaxy to fit various stages of its evolution with Jeans dynamical models. The Jeans models successfully fit the early oblate state before the major merger event, and also the late prolate stages of the simulated galaxy, recovering its mass distribution, velocity dispersion, and rotation profile. We also ran a prolate-rotating N-body simulation with similar properties to the cosmologically simulated galaxy, which gradually loses its angular momentum on a short time scale $\sim100$\,Myr. More tests are needed to understand why prolate rotation is time enduring in the cosmological simulation, but not in a simple N-body simulation.

Identifying the astrophysical sources responsible for the high-energy cosmic neutrinos has been a longstanding challenge. In a previous work, we report evidence for a spatial correlation between blazars from the 5th Roma-BZCat catalog and neutrino data of the highest detectable energies, i.e. >0.1 PeV, collected by the IceCube Observatory in the southern celestial hemisphere. The statistical significance is found at the level of 2 x 10^{-6} post-trial. In this work we test whether a similar correlation exists in the northern hemisphere, were IceCube is mostly sensitive to <0.1 PeV energies. We find a consistent correlation between blazars and northern neutrino data at the pre-trial p-value of 5.12 x 10^{-4}, and a post-trial chance probability of 6.79 x 10^{-3}. Combining the post-trial probabilities observed for the southern and northern experiments yields a global post-trial chance probability of 2.59 x 10^{-7} for the genuineness of such correlation. This implies that the spatial correlation is highly unlikely to arise by chance. Our studies push forward an all-sky subset of 52 objects as highly likely PeVatron extragalactic accelerators.

Kilometre-sized planetesimals form from pebbles of a range of sizes. We present the first simulations of the streaming instability that begin with a realistic, peaked size distribution, as expected from grain growth predictions. Our 3D numerical simulations directly form planetesimals via the gravitational collapse of pebble clouds. Models with multiple grain sizes show spatially distinct dust populations. The smallest grains in the size distribution do not participate in the formation of filaments or the planetesimals that are formed by the remaining ~80% of the dust mass. This implies a size cutoff for pebbles incorporated into asteroids and comets. Observations cannot resolve this dust clumping. However, we show that clumping, combined with optical depth effects, can cause significant underestimates of the dust mass, with 20%-80% more dust being present even at moderate optical depths if the streaming instability is active.

Fast radio bursts (FRBs) are bright radio transients of micro-to-millisecond duration and unknown extragalactic origin. Central to the mystery of FRBs are their extremely high characteristic energies, which surpass the typical energies of other radio transients of similar duration, like Galactic pulsar and magnetar bursts, by orders of magnitude. Calibration of FRB-detecting telescopes for burst flux and fluence determination is crucial for FRB science, as these measurements enable studies of the FRB energy and brightness distribution in comparison to progenitor theories. The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a radio interferometer of cylindrical design. This design leads to a high FRB detection rate but also leads to challenges for CHIME/FRB flux calibration. This paper presents a comprehensive review of these challenges, as well as the automated flux calibration software pipeline that was developed to calibrate bursts detected in the first CHIME/FRB catalog, consisting of 536 events detected between July 25th, 2018 and July 1st, 2019. We emphasize that, due to limitations in the localization of CHIME/FRB bursts, flux and fluence measurements produced by this pipeline are best interpreted as lower limits, with uncertainties on the limiting value.

We present observations of polarized dust emission at 850 $\mu$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $\mu$G in the main starless core and up to $\sim$90$\pm$40 $\mu$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfv\'enic. We also present a new method of data reduction for these denser but fainter objects like starless cores.

The Galactic dwarf spheroidal galaxies (dSphs) provide valuable insight into dark matter (DM) properties and its role in galaxy formation. Their close proximity enables the measurement of line-of-sight velocities for resolved stars, which allows us to study DM halo structure. However, uncertainties in DM mass profile determination persist due to the degeneracy between DM mass density and velocity dispersion tensor anisotropy. Overcoming this requires large kinematic samples and identification of foreground contamination. With 1.25 deg$^2$ and 2394 fibers, PFS plus pre-imaging with Hyper Suprime Cam will make significant progress in this undertaking.

Fits the infrared spectra from the nuclear regions of AGN can place constraints on the dust properties, distribution, and geometry by comparison with models. However, none of the currently available models fully describe the observations of AGN currently available. Among the aspects least explored, here we focus on the role of dust grain size. We offer the community a new spectral energy distribution (SED) library, hereinafter [GoMar23] model, which is based on the two-phase torus model developed before with the inclusion of the grain size as a model parameter, parameterized by the maximum grain size Psize or equivalently the mass-weighted average grain size < P >. We created 691,200 SEDs using the SKIRT code, where the maximum grain size can vary within the range Psize = 0.01 - 10.0um ( < P >= 0.007 - 3.41um). We fit this new and several existing libraries to a sample of 68 nearby and luminous AGNs with Spitzer/IRS spectra dominated by AGN-heated dust. We find that the [GoMar23] model can adequately reproduce up to 85-88% of the spectra. The dust grain size parameter significantly improves the final fit in up to 90% of these spectra. Statistical tests indicate that the grain size is the third most important parameter in the fitting procedure (after the size and half opening angle of the torus). The requirement of a foreground extinction by our model is lower compared to purely clumpy models. We find that 41% of our sample requires that the maximum dust grain size is as large as Psize =10um (< P >= 3.41um). Nonetheless, we also remark that disk+wind and clumpy torus models are still required to reproduce the spectra of a non-negligible fraction of objects, suggesting the need for several dust geometries to explain the infrared continuum of AGN. This work provides tentative evidence for dust grain growth in the proximity of the AGN.

Circumstellar discs likely have a short window when they are self-gravitating and prone to the effects of disc instability, but during this time the seeds of planet formation can be sown. It has long been argued that disc fragmentation can form large gas giant planets at wide orbital separations, but its place in the planet formation paradigm is hindered by a tendency to form especially large gas giants or brown dwarfs. We instead suggest that planet formation can occur early in massive discs, through the gravitational collapse of dust which can form the seeds of giant planets. This is different from the usual picture of self-gravitating discs, in which planet formation is considered through the gravitational collapse of the gas disc into a gas giant precursor. It is familiar in the sense that the core is formed first, and gas is accreted thereafter, as is the case in the core accretion scenario. However, by forming a $\sim 1 M_{\oplus}$ seed from the gravitational collapse of dust within a self-gravitating disc there exists the potential to overcome traditional growth barriers and form a planet within a few times $10^5$ years. The accretion of pebbles is most efficient with centimetre-sized dust, but the accretion of millimetre sizes can also result in formation within a Myr. Thus, if dust can grow to these sizes, planetary seeds formed within very young, massive discs could drastically reduce the timescale of planet formation and potentially explain the observed ring and gap structures in young discs.

While various indirect methods are used to detect exoplanets, one of the most effective and accurate methods is the transit method, which measures the brightness of a given star for periodic dips when an exoplanet is passing in front of the parent star. For systems with multiple transiting planets, the gravitational perturbations between planets affect their transit times. The difference in transit times allows a measurement of the planet masses and orbital eccentricities. These parameters help speculating on the formation, evolution and stability of the system. Using Transit Timing Variations (TTVs), we measure the masses and eccentricities of two planets orbiting K2-21, a relatively bright K7 dwarf star. These two planets exhibit measurable TTVs, have orbital periods of about 9.32 days and 15.50 days, respectively, and a period ratio of about 1.66, which is relatively near to the 5:3 mean motion resonance. We report that the inner and outer planets in the K2-21 system have properties consistent with the presence of a hydrogen and helium dominated atmospheres, as we estimate their masses to be 1.59^{+0.52}_{-0.44} M_E and 3.88^{+1.22}_{-1.07} M_E and densities of 0.22^{+0.05}_{-0.04} rho_E and 0.34^{+0.08}_{-0.06} rho_E, respectively (M_E and rho_E are the mass and density of Earth, respectively). Our results show that the inner planet is less dense than the outer planet; one more counter-intuitive exoplanetary system such as Kepler-105, LTT 1445, TOI-175 and Kepler-279 systems.

The relation between the event rate of long Gamma-Ray Bursts at low redshift and the star formation rate is still controversial, especially in the low-redshift end. Dong et al. confirmed that the Gamma-Ray Burst rate always exceeds the star formation rate at low-redshift of z < 1 in despite of the sample completeness. However, the reason of low-redshift excess is still unclear. Considering low-luminosity bursts with smaller redshift generally, we choose three Swift long burst samples and classify them into low- and high-luminosity bursts in order to check whether the low-redshift excess is existent and if the excess is biased by the sample size and completeness. To degenerate the redshift evolution from luminosity, we adopt the non-parametric method to study the event rate of the two types of long bursts in each sample. It is found that the high-luminosity burst rates are consistent with the star formation rate within the whole redshift range while the event rates of low-luminosity bursts exceed the star formation rate at low redshift of z < 1. Consequently, we conclude that the low-redshift excess is contributed by the low-luminosity bursts with possibly new origins unconnected with the star formation, which is also independent of the sample size and the sample completeness.

The main aim of this paper is to present the multi scalar field components as candidates to be the dark energy of the universe and their observational constraints. We start with the canonical Quintessence and Phantom fields with quadratic potentials and show that a more complex model should bear in mind to satisfy current cosmological observations. Then we present some implications for a combination of two fields, named as Quintom models. We consider two types of models, one as the sum of the quintessence and phantom potentials and other including an interacting term between fields. We find that adding one degree of freedom, by an interacting term, the dynamics enriches considerably and could lead to an improvement in the fit of $-2\ln\Delta \Like_{\rm max}= 5.13$, compared to $\Lambda$CDM. The resultant effective equation of state is now able to cross the phantom divide line, and in several cases present an oscillatory or discontinuous behavior, depending on the interaction value. The parameter constraints of the scalar field models (quintessence, phantom, quintom and interacting quintom) were performed using Cosmic Chronometers, Supernovae Ia and Baryon Acoustic Oscillations data; and the Log-Bayes factors were computed to compare the performance of the models. We show that single scalar fields may face serious troubles and hence the necessity of a more complex models, i.e. multiple fields.

Using quantum chemical calculations, we model the pathways for synthesizing two purine nucleobases, adenine and guanine, in the gas-phase interstellar environment, surrounded by neutral atomic hydrogen (HI). HI is found active in facilitating a series of fundamental proton transfer processes of organic synthesis, including bond formation, cyclization, dehydrogenation, and H migration. The reactive potential barriers were significantly reduced in the alternative pathways created by HI, leading to a remarkable increase in the reaction rate. The presence of HI also lowered the reactive activation temperature from 757.8 K to 131.5-147.0 K, indicating the thermodynamic feasibility of these pathways in star-forming regions where some of the reactants have been astronomically detected. Our findings suggest that HI may serve as an effective catalyst for interstellar organic synthesis.

We calculate the upper bounds of the population of theoretically stable Centaur orbits between Uranus and Neptune. These small bodies are on low-eccentricity, low-inclination orbits in two specific bands of semi-major axis, centred at $\sim$24.6 au and $\sim$25.6 au. They exhibit unusually long Gyr-stable lifetimes in previously published numerical integrations, orders of magnitude longer than that of a typical Centaur. Despite the increased breadth and depth of recent solar system surveys, no such objects have been found. Using the Outer Solar System Origins Survey (OSSOS) survey simulator to calculate the detection efficiency for these objects in an ensemble of fully characterised surveys, we determine that a population of 72 stable Centaurs with absolute magnitude $H_{r}\leq10$ ($95\%$ confidence upper limit) could remain undetected. The upcoming Legacy Survey of Space and Time (LSST) will be able to detect this entire intrinsic population due to its complete coverage of the ecliptic plane. If detected, these objects will be interesting dynamically-accessible mission targets -- especially as comparison of the stable Centaur orbital phase space to the outcomes of several modern planetary migration simulations suggests that these objects could be close to primordial in nature.

We perform high-resolution hydrodynamical simulations using the framework of {\it MACER} to investigate supermassive black hole (SMBH) feeding and feedback in a massive compact galaxy, which has a small effective radius but a large stellar mass, with a simulation duration of 10 Gyr. We compare the results with a reference galaxy with a similar stellar mass but a less concentrated stellar density distribution, as typically found in local elliptical galaxies. We find that about 10% of the time, the compact galaxy develops multi-phase gas within a few kpc, but the accretion flow through the inner boundary below the Bondi radius is always a single phase. The inflow rate in the compact galaxy is several times larger than in the reference galaxy, mainly due to the higher gas density caused by the more compact stellar distribution. Such a higher inflow rate results in stronger SMBH feeding and feedback and a larger fountain-like inflow-outflow structure. Compared to the reference galaxy, the star formation rate in the compact galaxy is roughly two orders of magnitude higher but is still low enough to be considered quiescent. Over the whole evolution period, the black hole mass grows by $\sim$50% in the compact galaxy, much larger than the value of $\sim$ 3% in the reference galaxy.

Stellar mass black holes in X-ray binaries (XRBs) are known to display different states characterized by different spectral and timing properties, understood in the framework of a hot corona coexisting with a thin accretion disk whose inner edge is truncated. There are several open questions related to the nature and properties of the corona, the thin disk, and dynamics behind the hard state. This motivated us to perform two-dimensional hydrodynamical simulations of accretion flows onto a 10 solar masses black hole. We consider a two-temperature plasma, incorporate radiative cooling with bremmstrahlung, synchrotron and comptonization losses and approximate the Schwarzschild spacetime via a pseudo-Newtonian potential. We varied the mass accretion rate in the range 0.02 <= Mdot/Mdot_Edd <= 0.35. Our simulations show the natural emergence of a colder truncated thin disk embedded in a hot corona, as required to explain the hard state of XRBs. We found that as Mdot increases, the corona contracts and the inner edge of the thin disk gets closer to the event horizon. At a critical accretion rate 0.02 <= Mdot_crit\Mdot_Edd <= 0.06, the thin disk disappears entirely. We discuss how our simulations compare with XRB observations in the hard state.

We develop a deep learning technique to reconstruct the dark matter density field from the redshift-space distribution of dark matter halos. We implement a UNet-architecture neural network and successfully trained it using the COLA fast simulation, which is an approximation of the N-body simulation with $512^3$ particles in a box size of $500 h^{-1}{\rm {Mpc}}$. We evaluate the resulting UNet model not only using the training-like test samples, but also using the typical N-body simulations, including the Jiutian simulation which has $6144^3$ particles in a box size of $1000 h^{-1}{\rm {Mpc}}$, and the ELUCID simulation which has a different cosmology. The real-space dark matter density fields in the three simulations can all be recovered consistently with only a small reduction of the cross-correlation power spectrum at 1\% and 10\% levels at $k=0.1$ and $0.3~h\mathrm{Mpc^{-1}}$, respectively. It is evident that the reconstruction helps to correct for the redshift-space distortions and is unaffected by the different cosmologies between the training sample ({\bf Planck2018}) and the test sample ({\bf WMAP5}). In addition, we tested the application of the UNet-reconstructed density field to recover the velocity \& tidal field and found it outperforms the traditional approach based on the linear bias model, showing a 12.2 percent improvement in the correlation slope and a 21.1 percent reduction in the scatter between the predicted and the true velocities. As a result, our method is highly efficient and has an outstanding level of extrapolation reliability beyond the training set. This offers an optimal solution that determines the three-dimensional underlying density field from the abundant galaxy survey data.

The 21cm global signal is an important probe to reveal the properties of the first astrophysical objects and the processes of the structure formation from which one can constrain astrophysical and cosmological parameters. To extract the information of such parameters, one needs to efficiently evaluate the 21cm global signal for statistical analysis. First we developed an artificial neural network-based emulator to predict the 21cm global signal, which works with significantly less computational cost and high precision. Then we apply our emulator to demonstrate the parameter estimation based on the Bayesian analysis by using the publicly available EDGES low-band data. We find that the result is sensitive to the foreground model, the assumption of noise, and the frequency range used in the analysis. The Bayesian evidence suggests the models with higher order polynomial function and enhanced noise are preferred. We also compare models suggested from the EDGES low-band data and the ones from recent JWST measurements of the galaxy luminosity function at $z=16$. We find that the model which produces the 21cm absorption line at $z\approx15$ is well consistent with the central value of the observed luminosity function at $z=16$.

Temporal broadening is a commonly observed property of fast radio bursts (FRBs), associated with turbulent media which cause radiowave scattering. Similarly to dispersion, scattering is an important probe of the media along the line of sight to an FRB source, such as the circum-burst or circum-galactic mediums (CGM). Measurements of characteristic scattering times alone are insufficient to constrain the position of the dominant scattering media along the line of sight. However, where more than one scattering screen exists, Galactic scintillation can be leveraged to form strong constraints. We quantify the scattering and scintillation in 10 FRBs with 1) known host galaxies and redshifts and 2) captured voltage data enabling high time resolution analysis, obtained from the Commensal Real-time ASKAP (Australian Square Kilometre Array Pathfinder) Fast Transient survey science project (CRAFT). We find strong evidence for two screens in three cases. For FRBs 20190608B and 20210320C, we find evidence for scattering screens less than approximately 16.7 and 3000 kpc respectively, from their sources. For FRB 20201124A we find evidence for a scattering screen at $\approx$26 kpc. Each of these measures is consistent with the scattering occurring in the host ISM (inter-stellar medium) or CGM. If pulse broadening is assumed to be contributed by the host galaxy ISM or circum-burst environment, the definitive lack of observed scintillation in four FRBs in our sample suggests that existing models may be over-estimating scattering times associated with the Milky Way's ISM, similar to the anomalously low scattering observed for FRB 20201124A.

Radiometer experiments to detect 21-cm Hydrogen line emission from the Cosmic Dawn and Epoch of Reionization rely upon precise absolute calibration. During calibration, noise generated by amplifiers within the radiometer receiver must be accounted for; however, it is difficult to measure as the noise power varies with source impedance. In this letter, we introduce a convenient method to measure the noise parameters of a receiver system, which is practical for low-frequency receivers used in global 21-cm experiments.

HEALPix by G\'orski et. al. (2005) is de-facto standard for Cosmic Microwave Background (CMB) data storage and analysis, and is widely used in current and upcoming CMB experiments. Almost all the datasets in Legacy Archive for Microwave Background Data Analysis (LAMBDA) use HEALPix as a format of choice. Visualizing the data plays important role in research, and several toolsets were developed to do that with HEALPix maps, most notably original Fortran facilities and Python integration with healpy. With the current state of GPU performance, it is now possible to visualize extremely large maps in real time on a laptop or a tablet. HEALPix Viewer described here is developed for macOS, and takes full advantage of GPU acceleration to handle extremely large datasets in real time. It compiles natively on Intel and Arm64 architectures, and uses Metal framework for high-performance GPU computations. The aim of this project is to reduce the effort required for interactive data exploration, as well as time overhead for producing publication-quality maps. Drag and drop integration with Keynote and Powerpoint makes creating presentations easy. The main codebase is written in Swift, a modern and efficient compiled language, with high-performance computing parts delegated entirely to GPU, and a few inserts in C interfacing to cfitsio library for I/O. Graphical user interface is written in SwiftUI, a new declarative UI framework based on Swift. Most common spherical projections and colormaps are supported out of the box, and the available source code makes it easy to customize the application and to add new features if desired.

The Sun produces the most powerful explosions in the solar system, solar flares, that can also be accompanied by large eruptions of magnetised plasma, coronal mass ejections (CMEs). These processes can accelerate electron beams up to relativistic energies through magnetic reconnection processes during solar flares and CME-driven shocks. Energetic electron beams can in turn generate radio bursts through the plasma emission mechanism. CME shocks, in particular, are usually associated with type II solar radio bursts. However, on a few occasions, type II bursts have been reported to occur either in the absence of CMEs or shown to be more likely related with the flaring process. It is currently an open question how a shock generating type II bursts forms without the occurrence of a CME eruption. Here, we aim to determine the physical mechanism responsible for a type II burst which occurs in the absence a CME. By using radio imaging from the Nan{\c c}ay Radioheliograph, combined with observations from the Solar Dynamics Observatory and the Solar Terrestrial Relations Observatory spacecraft, we investigate the origin of a type II radio burst that appears to have no temporal association with a white-light CME. We identify a typical type II radio burst with band-split structure that is associated with a C-class solar flare. The type II burst source is located above the flaring active region and ahead of disturbed coronal loops observed in extreme ultraviolet images. The type II is also preceded by type III radio bursts, some of which are in fact J-bursts indicating that accelerated electron beams do not all escape along open field lines. The type II sources show single-frequency movement towards the flaring active region. The type II is located ahead of a faint extreme-ultraviolet (EUV) front propagating through the corona.

Dense star clusters are spectacular self-gravitating stellar systems in our Galaxy and across the Universe - in many respects. They populate disks and spheroids of galaxies as well as almost every galactic center. In massive elliptical galaxies nuclear clusters harbor supermassive black holes, which might influence the evolution of their host galaxies as a whole. The evolution of dense star clusters is not only governed by the aging of their stellar populations and simple Newtonian dynamics. For increasing particle number, unique gravitational effects of collisional many-body systems begin to dominate the early cluster evolution. As a result, stellar densities become so high that stars can interact and collide, stellar evolution and binary stars change the dynamical evolution, black holes can accumulate in their centers and merge with relativistic effects becoming important. Recent high-resolution imaging has revealed even more complex structural properties with respect to stellar populations, binary fractions and compact objects as well as - the still controversial - existence of intermediate mass black holes in clusters of intermediate mass. Dense star clusters therefore are the ideal laboratory for the concomitant study of stellar evolution and Newtonian as well as relativistic dynamics. Not only the formation and disruption of dense star clusters has to be considered but also their galactic environments in terms of initial conditions as well as their impact on galactic evolution. This review deals with the specific computational challenges for modelling dense, gravothermal star clusters.

Context: The study of high-mass stars found to be isolated in the field of the Milky Way may help to probe the feasibility of the core-accretion mechanism in the case of massive star formation. The existence of truly isolated stars may efficiently probe the possibility that individual massive stars can be born in isolation. Aims: We observed WR67a (hereafter Sapaki), an O3If* star that appears to be isolated close to the center of a well-developed giant cavity that is aptly traced by 8.0 $\mu$m hot dust emission. Methods: We acquired medium-resolution ($R=4100$) and moderate signal-to-noise ($S/N = 95$ at 4500 \r{A}) spectra for Sapaki in the range of 3800-10500 \r{A} with the Magellan Echellette (MagE) at Las Campanas Observatory. We computed the line-of-sight total extinctions. Additionally, we restricted its heliocentric distance by using a range of different estimators. Moreover, we measured its radial velocity from several lines in its spectrum. Finally, we analyzed its proper motions from Gaia to examine its possible runaway status. Results: The star has been classified as having the spectral type O3If* given its resemblance to standard examples of the class. In addition, we found that Sapaki is highly obscured, reaching a line-of-sight extinction value of $A_{V} = 7.87$. We estimated the heliocentric distance to be in the range of $d = 4-7$ kpc. We also estimated its radial velocity to be $V_{r} = -34.2 \pm 15.6$ km/s. We may also discard its runaway status solely based on its 2D kinematics. Furthermore, by analyzing proper motions and parallaxes provided by Gaia, we found only one other star with compatible measurements. Conclusions: Given its apparent non-runaway status and the absence of clustering, Sapaki appears to be a solid candidate for isolated high-mass star formation in the Milky Way.

Data-driven analysis methods can help to infer physical properties of red giant stars where "gold-standard" asteroseismic data are not available. The study of optical and infrared spectra of red giant stars with data-driven analyses has revealed that differences in oscillation frequencies and their separations are imprinted in said spectra. This makes it possible to confidently differentiate core-helium burning red clump stars (RC) from those that are still on their first ascent of the red giant branch (RGB). We extend these studies to a tenfold larger wavelength range of 0.33 to 2.5 microns with the moderate-resolution VLT/X-shooter spectrograph. Our analysis of 49 stars with asteroseismic data from the K2 mission confirms that CN, CO and CH features are indeed the primary carriers of spectroscopic information on the evolutionary stages of red giant stars. We report 215 informative features for differentiating the RC from the RGB within the range of 0.33 to 2.5 microns. This makes it possible for existing and future spectroscopic surveys to optimize their wavelength regions to deliver both a large variety of elemental abundances and reliable age estimates of luminous red giant stars.

Methylamine has been the only simple alkylamine detected in the interstellar medium for a long time. With the recent secure and tentative detections of vinylamine and ethylamine, respectively, dimethylamine has become a promising target for searches in space. Its rotational spectrum, however, has been known only up to 45 GHz until now. Here we investigate the rotation-tunneling spectrum of dimethylamine in selected regions between 76 and 1091 GHz using three different spectrometers in order to facilitate its detection in space. The quantum number range is extended to $J = 61$ and $K_a = 21$, yielding an extensive set of accurate spectroscopic parameters. To search for dimethylamine, we refer to the spectral line survey ReMoCA carried out with the Atacama Large Millimeter/submillimeter Array toward the high-mass star-forming region Sagittarius B2(N) and a spectral line survey of the molecular cloud G+0.693$-$0.027 employing the IRAM 30 m and Yebes 40 m radio telescopes. We report nondetections of dimethylamine toward the hot molecular cores Sgr B2(N1S) and Sgr B2(N2b) as well as G+0.693$-$0.027 which imply that dimethylamine is at least 14, 4.5 and 39 times less abundant than methylamine toward these sources, respectively. The observational results are compared to computational results from a gas-grain astrochemical model. The modeled methylamine to dimethylamine ratios are compatible with the observational ratios. However, the model produces too much ethylamine compared with methylamine which could mean that the already fairly low levels of dimethylamine in the models may also be too high.

We present a novel Hall magnetohydrodynamics (HMHD) numerical simulation of a three-dimensional (3D) magnetic flux rope (MFR) -- generated by magnetic reconnections from an initial 3D bipolar sheared field. Magnetic reconnections during the HMHD evolution are compared with the MHD. In both simulations, the MFRs generate as a consequence of the magnetic reconnection at null points which has not been realized in contemporary simulations. Interestingly, the evolution is faster and more intricate in the HMHD simulation. Repetitive development of the twisted magnetic field lines (MFL) in the vicinity of 3D nulls (reconnection site) is unique to the HMHD evolution of the MFR. The dynamical evolution of magnetic field lines around the reconnection site being affected by the Hall forcing, correspondingly affects the large-scale structures.

Mass loss from massive stars plays a determining role in their evolution through the upper Hertzsprung-Russell diagram. The hydrodynamic theory that describes their steady-state winds is the line-driven wind theory (m-CAK). From this theory, the mass loss rate and the velocity profile of the wind can be derived, and estimating these properly will have a profound impact on quantitative spectroscopy analyses from the spectra of these objects. Currently, the so-called beta-law, which is an approximation for the fast solution, is widely used instead of m-CAK hydrodynamics, and when the derived value is beta greater than 1.2, there is no hydrodynamic justification for these values. This review focuses on (1) a detailed topological analysis of the equation of motion (EoM), (2) solving the EoM numerically for all three different (fast and two slow) wind solutions, (3) deriving analytical approximations for the velocity profile via the LambertW function and (4) presenting a discussion of the applicability of the slow solutions.

We have studied the effect of changing the density and magnetic field strength in the coherent pulses that are emitted as energetic showers develop in the atmosphere. For this purpose we have developed an extension of ZHS, a program to calculate coherent radio pulses from electromagnetic showers in homogeneous media, to account for the Lorentz force due to a magnetic field. This makes it possible to perform quite realistic simulations of radio pulses from air showers in a medium similar to the atmosphere but without variations of density with altitude. The effects of independently changing the density, the refractive index and the magnetic field strength are studied in the frequency domain for observers in the Cherenkov direction at far distances from the shower. This approach is particularly enlightening providing an explanation of the spectral behavior of the induced electric field in terms of shower development parameters. More importantly, it clearly displays the complex scaling properties of the pulses as density and magnetic field intensity are varied. The usually assumed linear behavior of electric field amplitude with magnetic field intensity is shown to hold up to a given magnetic field strength at which the extra time delays due to the deflection in the magnetic field break it. Scaling properties of the pulses are obtained as the density of air decreases relative to sea level. A remarkably accurate scaling law is obtained that relates the spectra of pulses obtained when reducing the density and increasing the magnetic field.

The evolutionary relationships and mechanisms governing the behavior of the wide variety of luminous stars populating the upper H-R diagram are not well established. Luminous blue variables (LBVs) are particularly rare, with only a few dozen identified in the Milky Way and nearby galaxies. Since 2012, the Barber Observatory Luminous Stars Survey has monitored more than 100 luminous targets in M33, including M33C-4119 which has recently undergone photometric and spectroscopic changes consistent with an S Doradus eruption of an LBV.

Coronal rain is the most dramatic cooling phenomenon of the solar corona and an essential diagnostic tool for the coronal heating properties. A puzzling feature of the solar corona, besides the heating, is its EUV filamentary structure and variability. We aim to identify observable features of the TNE-TI scenario underlying coronal rain at small and large spatial scales, to understand the role it plays in the solar corona. We use EUV datasets at unprecedented spatial resolution of ~240 km from EUI/HRIEUV and SPICE of Solar Orbiter from the spring 2022 perihelion. EUV absorption features produced by coronal rain are detected at scales as small as 260 km. As the rain falls, heating and compression is produced immediately downstream, leading to a small EUV brightening accompanying the fall and producing a "fireball" phenomenon. Just prior to impact, a flash-like EUV brightening downstream of the rain, lasting a few minutes is observed for the fastest events. For the first time, we detect the atmospheric response to the rain's impact on the chromosphere and consists of upward propagating rebound shocks and flows partly reheating the loop. The observed widths of the rain clumps are 500 +- 200 km. They exhibit a broad velocity distribution of 10 - 150 km s^-1, peaking below 50 km s^-1. Coronal strands of similar widths are observed along the same loops co-spatial with cool filamentary structure, which we interpret as the CCTR. Matching with the expected cooling, prior to the rain appearance sequential loop brightenings are detected in gradually cooler lines from corona to chromospheric temperatures. Despite the large rain showers, most cannot be detected in AIA 171 in quadrature, indicating that LOS effects play a major role in coronal rain visibility. Still, AIA 304 and SPICE observations reveal that only a small fraction of the rain can be captured by HRIEUV.

The term `$\alpha$-meteoroid' was introduced to describe a group of micrometeoroids with certain dynamical properties, which -- alongside the group of the $\beta$-meteoroids -- had been identified by the first generation of reliable in-situ dust detectors in interplanetary space. In recent years, use of the term $\alpha$-meteoroid has become more frequent again, under a subtly but crucially altered definition. This work shall bring attention to the discrepancy between the term's original and newly established meaning, and spotlight the now-overlooked group of particles that the term used to describe. We review past and present pertinent literature around the term $\alpha$-meteoroid, and assess the dynamics of the originally referred-to particles with respect to possible sources, showing that their formation is the expected consequence of collisional grinding of the zodiacal cloud at short heliocentric distances. The abundance of the original $\alpha$-meteoroids, which are essentially `bound $\beta$-meteoroids', makes them relevant to all in-situ dust experiments in the inner solar system. Due to the change of the term's meaning, however, they are not considered by contemporary studies. The characterization of this particle population could elucidate the processing of the innermost zodiacal cloud, and should thus be objective of upcoming in-situ dust experiments. The attained ambiguity of the term $\alpha$-meteoroid is not easily resolved, warranting great care and clarity going forward.

We combine Hubble Space Telescope (HST) Fine Guidance Sensor, Hipparcos, and Gaia DR3 astrometric observations of the K0 V star 14 Her with the results of an analysis of extensive ground-based radial velocity data to determine perturbation orbits and masses for two previously known companions, 14 Her b and c. Radial velocities obtained with the Hobby-Eberly Telescope and from the literature now span over twenty five years. With these data we obtain improved RV orbital elements for both the inner companion, 14 Her b and the long-period outer companion, 14 Her c. We also find evidence of an additional RV signal with P $/sim$ 3789d. We then model astrometry from Hipparcos, HST, and Gaia with RV results to obtain system parallax and proper motion, perturbation periods, inclinations, and sizes due to 14 Her b and c. We find P_b = 1767.6 +/- 0.2 d, perturbation semi-major axis {\alpha}_b = 1.3 +/- 0.1 mas, and inclination i_b = 36 +/- 3 degrees, P_c = 52160 +/- 1028 d, perturbation semi-major axis {\alpha}_c = 10.3 +/- 0.7 mas, and inclination i_c = 82 +/- 14 degrees. In agreement with a past investigation, the 14 Her b, c orbits exhibit significant mutual inclination. Assuming a primary mass M = 0.98 +/- 0.04Msun, we obtain companion masses M_b = 8.5 +/- 1.0Mjup and M_c = 7.1 +/- 1.0Mjup.

Models of planet-disk interaction are mainly based on 2D and 3D viscous hydrodynamical simulations. Accretion is classically prescribed by an alpha parameter which characterizes the turbulent radial transport of angular momentum (AM) in the disk. This accretion scenario has been questioned for a few years and an alternative paradigm has been proposed that involves the vertical transport of AM by MHD winds. We revisit planet-disk interaction in such context, with a focus on the planet's ability to open a gap and produce meridional flows. Accretion, magnetic field and wind torque in the gap are also explored, as well as the gravitational torque exerted by the disk onto the planet. We carry out high-resolution 3D global non-ideal MHD simulations of a gaseous disk threaded by a large-scale vertical magnetic field harboring a planet in a fixed circular orbit using the GPU-accelerated code Idefix. We consider various planet masses and disk magnetizations. We find that gap-opening always occurs for sufficiently massive planets, with deeper gaps when the planet mass increases and when the initial magnetization decreases. We propose a gap opening criterion when accretion is dominated by MHD winds. We show that accretion is unsteady and comes from surface layers in the outer disk, bringing material directly towards the planet poles. A planet gap is a privileged region for magnetic field accumulation, leading to nearly sonic accretion stream through the gap. For massive planets, the wind torque induces an asymmetric gap, both in depth and in width, that gradually erodes the outer gap edge, reducing the outer Lindblad torque and potentially reversing the migration direction of Jovian planets in magnetized disks after a few hundreds of orbits. For low-mass planets, we find strongly fluctuating gravitational torques that are mostly positive on average, indicating a stochastic outward migration.

While supermassive black holes are ubiquitous features of galactic nuclei, only a small minority are observed during episodes of luminous accretion. The physical mechanism(s) driving the onset of fueling and ignition in these active galactic nuclei (AGN) are still largely unknown for many galaxies and AGN-selection criteria. Attention has focused on AGN triggering by means of major galaxy mergers gravitationally funneling gas towards the galactic center, with evidence both for and against this scenario. However, several recent studies have found that radio-loud AGN overwhelmingly reside in ongoing or recent major galaxy mergers. In this study, we test the hypothesis that major galaxy mergers are important triggers for radio-loud AGN activity in powerful quasars during cosmic noon (1 < z < 2). To this end, we compare Hubble Space Telescope WFC3/IR observations of the z > 1 3CR radio-loud broad-lined quasars to three matched radio-quiet quasar control samples. We find strong evidence for major-merger activity in nearly all radio-loud AGN, in contrast to the much lower merger fraction in the radio-quiet AGN. These results suggest major galaxy mergers are key ingredients to launching powerful radio jets. Given many of our radio-loud quasars are blue, our results present a possible challenge to the "blow-out" paradigm of galaxy evolution models in which blue quasars are the quiescent end result following a period of red quasar feedback initiated by a galaxy merger. Finally, we find a tight correlation between black hole mass and host galaxy luminosity for these different high-redshift AGN samples inconsistent with those observed for local elliptical galaxies.

We outline and discuss a model for the enhanced abundances of trans-Fe elements in impulsive Solar Energetic Particle (SEP) events, where large mass dependent abundance enhancements are frequently seen. It comes about as a variation of the ponderomotive force model for the First Ionization Potential (FIP) Effect, i.e. the increase in coronal abundance of elements like Fe, Mg, and Si that are ionized in the solar chromosphere relative to those that are neutral. In this way, the fractionation region is placed in the chromosphere, and is connected to the solar envelope allowing the huge abundance variations to occur, that might otherwise be problematic with a coronal fractionation site. The principal mechanism behind the mass-independent FIP fractionation becoming the mass dependent impulsive SEP fractionation is the suppression of acoustic waves in the chromosphere. The ponderomotive force causing the fractionation must be due to torsional Alfven waves, which couple much less effectively to slow modes than do shear waves, and upward propagating acoustic waves deriving from photospheric convection must be effectively mode converted to fast modes at the chromospheric layer where Alfven and sound speeds are equal, and subsequently totally internally reflected. We further discuss observations of the environments thought to be the source of impulsive SEPs, and the extent to which the real Sun might meet these conditions.

We develop statistical methods within a Bayesian framework to infer the star formation history from photometric surveys of pre-main sequence populations. Our procedures include correcting for biases due to extinction in magnitude-limited surveys, and using distributions from subsets of stars with individual extinction measurements. We also make modest corrections for unresolved binaries. We apply our methods to samples of populations with Gaia photometry in the Orion A molecular cloud. Using two well-established sets of evolutionary tracks, we find that, although our sample is incomplete at youngest ages due to extinction, star formation has proceeded in Orion A at a relatively constant rate between ages of about 0.3 and 5 Myr, in contrast to other studies suggesting multiple epochs of star formation. Similar results are obtained for a set of tracks that attempt to take the effects of strong magnetic fields into account. We also find no evidence for a well-constrained "birthline" that would result from low-mass stars appearing first along the deuterium-burning main sequence, especially using the magnetic evolutionary tracks. While our methods have been developed to deal with Gaia data, they may be useful for analyzing other photometric surveys of star-forming regions.

The novel Gaia Multi Peak (GMP) technique has proven to be able to successfully select dual and lensed AGN candidates at sub-arcsec separations. Both populations are important because dual AGNs represent one of the central, still largely untested, predictions of lamdaCDM cosmology, and compact lensed quasars allow to probe the central regions of the lensing galaxies. In this work, we present high spatial resolution spectroscopy of twelve GMP-selected systems. We use the the adaptive-optics assisted integral-field spectrograph MUSE at VLT to resolve each system and study the nature of each component. All the targets reveal the presence of two components confirming the GMP selection. We classify five targets as dual AGNs, two as lensed systems, and five as a chance alignment of a star and and AGN. Having separations between 0.30" and 0.86", these dual and lensed systems are, to date, among the most compact ever discovered at z >0.3. This is the largest sample of distant dual AGNs with sub-arcsec separations ever presented in a single paper.

Upcoming photometric surveys will discover tens of thousands of Type Ia supernovae (SNe Ia), vastly outpacing the capacity of our spectroscopic resources. In order to maximize the science return of these observations in the absence of spectroscopic information, we must accurately extract key parameters, such as SN redshifts, with photometric information alone. We present Photo-zSNthesis, a convolutional neural network-based method for predicting full redshift probability distributions from multi-band supernova lightcurves, tested on both simulated Sloan Digital Sky Survey (SDSS) and Vera C. Rubin Legacy Survey of Space and Time (LSST) data as well as observed SDSS SNe. We show major improvements over predictions from existing methods on both simulations and real observations as well as minimal redshift-dependent bias, which is a challenge due to selection effects, e.g. Malmquist bias. The PDFs produced by this method are well-constrained and will maximize the cosmological constraining power of photometric SNe Ia samples.

We consider optical fibers as detectors for scalar ultralight dark matter (UDM) and propose using a fiber-based interferometer to search for scalar UDM with particle mass in the range $10^{-17} - 10^{-13}$ eV/$c^2$ $\left(10^{-3}- 10 \text{ Hz}\right)$. Composed of a solid core and a hollow core fiber, the proposed detector would be sensitive to relative oscillations in the fibers' refractive indices due to scalar UDM-induced modulations in the fine-structure constant $\alpha$. We predict that, implementing detector arrays or cryogenic cooling, the proposed optical fiber-based scalar UDM search has the potential to reach new regions of the parameter space. Such a search would be particularly well-suited to probe for a Solar halo of dark matter with a sensitivity exceeding that of previous DM searches over the particle mass range $7\times 10^{-17} - 2\times 10^{-14}$ eV/$c^2$.

In asymptotically flat spacetimes, bearing the null geodesics reaching the future null infinity in mind, we propose new concepts, the ``dark horizons'' as generalizations of the photon sphere. They are defined in terms of the structure of escape/capture cones of photons with respect to a unit timelike vector field. More specifically, considering a two-sphere that represents a set of emission directions of photons, the dark horizons are located at positions where a hemisphere is marginally included in the capture and escape cones, respectively. We show that both of them are absent in the Minkowski spacetime, while they exist in spacetimes with black hole(s) under a certain condition. We derive the general properties of the dark horizons in spherically symmetric spacetimes and explicitly calculate the locations of the dark horizons in the Vaidya spacetime and the Kerr spacetime.

Cosmological correlators from inflation are often generated at tree level and hence loop contributions are bounded to be small corrections by perturbativity. Here we discuss a scenario where this is not the case. Recently, it has been shown that for any number of scalar fields of any mass, the parity-odd trispectrum of a massless scalar must vanish in the limit of exact scale invariance due to unitarity and the choice of initial state. By carefully handling UV-divergences, we show that the one-loop contribution is non-vanishing and hence leading. Surprisingly, the one-loop parity-odd trispectrum is simply a rational function of kinematics, which we compute explicitly in a series of models, including single-clock inflation. Although the loop contribution is the leading term in the parity-odd sector, its signal-to-noise ratio is typically bounded from above by that of a corresponding tree-level parity-even trispectrum, unless instrumental noise and systematics for the two observables differ. Furthermore, we identify a series of loop contributions to the wavefunction that cancel exactly when computing correlators, suggesting a more general phenomenon.

Evidence Networks can enable Bayesian model comparison when state-of-the-art methods (e.g. nested sampling) fail and even when likelihoods or priors are intractable or unknown. Bayesian model comparison, i.e. the computation of Bayes factors or evidence ratios, can be cast as an optimization problem. Though the Bayesian interpretation of optimal classification is well-known, here we change perspective and present classes of loss functions that result in fast, amortized neural estimators that directly estimate convenient functions of the Bayes factor. This mitigates numerical inaccuracies associated with estimating individual model probabilities. We introduce the leaky parity-odd power (l-POP) transform, leading to the novel ``l-POP-Exponential'' loss function. We explore neural density estimation for data probability in different models, showing it to be less accurate and scalable than Evidence Networks. Multiple real-world and synthetic examples illustrate that Evidence Networks are explicitly independent of dimensionality of the parameter space and scale mildly with the complexity of the posterior probability density function. This simple yet powerful approach has broad implications for model inference tasks. As an application of Evidence Networks to real-world data we compute the Bayes factor for two models with gravitational lensing data of the Dark Energy Survey. We briefly discuss applications of our methods to other, related problems of model comparison and evaluation in implicit inference settings.

An oscillating inflaton field induces small amplitude oscillations of the Hubble parameter at the end of inflation. These Hubble parameter induced oscillations, in turn, trigger parametric particle production of all light fields, even if they are not directly coupled to the inflaton. We here study the induced particle production for a light scalar field (e.g. the Standard Model Higgs field) after inflation as a consequence of this effect. Our analysis yields a model-independent lower bound on the efficiency of energy transfer from the inflaton condensate to particle excitations.

Rapid progress in electromagnetic black hole observation presents a theoretical challenge: how can the universal signatures of extreme gravitational lensing be distilled from stochastic astrophysical signals? With this motivation, the two-point correlation function of specific intensity fluctuations across image positions, times, and frequencies is here considered. The contribution of strongly deflected light rays, those which make up the photon ring, is analytically computed for a Kerr black hole illuminated by a simple geometric-statistical emission model. We subsequently integrate over the image to yield a spectro-temporal correlation function which is relevant for unresolved sources. Finally, some observational aspects are discussed and a preliminary assessment of detectability with current and upcoming missions is provided.

Experiments aimed at searching for variations in the fine-structure constant $\alpha$ are based on spectroscopy of transitions in microscopic bound systems, such as atoms and ions, or resonances in optical cavities. The sensitivities of these systems to variations in $\alpha$ are typically on the order of unity and are fixed for a given system. For heavy atoms, highly charged ions and nuclear transitions, the sensitivity can be increased by benefiting from the relativistic effects and favorable arrangement of quantum states. This article proposes a new method for controlling the sensitivity factor of macroscopic physical systems. Specific concepts of optical cavities with tunable sensitivity to $\alpha$ are described. These systems show qualitatively different properties from those of previous studies of the sensitivity of macroscopic systems to variations in $\alpha$, in which the sensitivity was found to be fixed and fundamentally limited to an order of unity. Although possible experimental constraints attainable with the specific optical cavity arrangements proposed in this article do not yet exceed the present best constraints on $\alpha$ variations, this work paves the way for developing new approaches to searching for variations in the fundamental constants of physics.

Single field models of inflation capable to produce primordial black holes usually require a significant departure from the standard, perturbative slow-roll regime. In fact, in many of these scenarios, the size of the slow-roll parameter $|\eta|$ becomes larger than one during a short phase of inflationary evolution. In order to develop an analytical control on these systems, we explore the limit of $|\eta|$ large, and promote $1/|\eta|$ to a small quantity to be used for perturbative expansions. Formulas simplify, and we obtain analytic expressions for the two and three point functions of curvature fluctuations, which share some of the features found in realistic inflationary models generating primordial black holes. We study one-loop corrections in this framework: we discuss criteria for adsorbing ultraviolet divergences into the available parameters, leaving log-enhanced infrared contributions of controllable size.

Cosmic string cusps are sources of short-lived, linearly polarised gravitational wave bursts which can be searched for in gravitational wave detectors. We assess the capability of LISA to detect these bursts using the latest LISA configuration and operational assumptions. For such short bursts, we verify that LISA can be considered as ``frozen", namely that one can neglect LISA's orbital motion. We consider two models for the network of cosmic string loops, and estimate that LISA should be able to detect 1-3 bursts per year assuming a string tension $G\mu \approx 10^{-11} - 10^{-10.5}$ and detection threshold $\rm{SNR} \ge 20$. Non-detection of these bursts would constrain the string tension to $G\mu\lesssim 10^{-11}$ for both models.

We explore the gravitational wave spectrum generated by string-wall structures in an $SO(10)$ ($Spin(10)$) based scenario of pseudo-Goldstone boson dark matter (pGDM) particle. This dark matter candidate is a linear combination of the Standard Model (SM) singlets present in the 126 and 16 dimensional Higgs fields. The Higgs $126$-plet vacuum expectation value (VEV) $\left<126_H\right>$ leaves unbroken the $\mathbb{Z}_2$ subgroup of $\mathbb{Z}_4$, the center of $SO(10)$. Among other things, this yields topologically stable cosmic strings with a string tension $\mu \sim \left<126_H\right>^2$. The subsequent (spontaneous) breaking of $\mathbb{Z}_2$ at a significantly lower scale by the $16$-plet VEV $\left<16_H\right>$ leads to the appearance of domain walls bounded by the strings produced earlier. We display the gravitational wave spectrum for $G \mu$ values varying between $10^{-15}$ and $10^{-9}$ ($\left<126_H\right>\sim 10^{11}$ - $10^{14}$ GeV), and $\left<16_H\right>\sim 0.1$ - $10^3$ TeV range ($G$ denotes Newton's constant.) These predictions can be tested, as we show, by a variety of (proposed) experiments including LISA, ET, CE and others.