Papers for Wednesday, Oct 12 2022

Quasar absorption line analysis is critical for studying gas and dust components and their physical and chemical properties as well as the evolution and formation of galaxies in the early universe. Ca II absorbers, which are one of the dustiest absorbers and are located at lower redshifts than most other absorbers, are especially valuable when studying physical processes and conditions in recent galaxies. However, the number of known quasar Ca II absorbers is relatively low due to the difficulty of detecting them with traditional methods. In this work, we developed an accurate and quick approach to search for Ca II absorption lines using deep learning. In our deep learning model, a convolutional neural network, tuned using simulated data, is used for the classification task. The simulated training data are generated by inserting artificial Ca II absorption lines into original quasar spectra from the Sloan Digital Sky Survey (SDSS) whilst an existing Ca II catalog is adopted as the test set. The resulting model achieves an accuracy of 96% on the real data in the test set. Our solution runs thousands of times faster than traditional methods, taking a fraction of a second to analyze thousands of quasars while traditional methods may take days to weeks. The trained neural network is applied to quasar spectra from SDSS's DR7 and DR12 and discovered 399 new quasar Ca II absorbers. In addition, we confirmed 409 known quasar Ca II absorbers identified previously by other research groups through traditional methods.

Both the origin of, and the population demographics of, massive black holes (MBHs) remains an open question in modern day astrophysics. Here we introduce the BlackDemon suite of cosmological simulations using the Enzo code. The suite consists of three, high resolution, distinct regions, each with a side length of 1 $\rm{h}^{-1}$ Mpc, evolving within a larger overdense region. The simulation suite has spatial and mass resolution capable of resolving the formation of the first galaxies and MBHs within each region. We report here, as the first in a series of papers, the evolution of the simulation suite up to the point where the first halo in each distinct region is undergoing gravitational collapse. The first halo to collapse has a total virial mass of M$_{\rm{halo}} \sim 1 \times 10^7 \ \rm{M}_{\odot}$ and collapses at $z \sim 23$. The other two haloes (each in a separate region) have masses of M$_{\rm{halo}} \sim 2 \times 10^6 \ \rm{M}_{\odot}$ and M$_{\rm{halo}} \sim 7 \times 10^6 \ \rm{M}_{\odot}$ respectively. We find that a significant number of haloes across each region are undergoing rapid assembly and experiencing extremely high rapid mass inflow which in turn is leading to longer cooling times. Given the high mass inflow rates we estimate that stars with masses between $10^3$ and $10^5$ $\rm{M}_{\odot}$ are likely to be a generic outcome in these high accretion rate haloes which form in rare overdense regions. Our probing of the high-mass end tail of the population III initial mass function means we are coming full circle with the initial predicted masses of the first stars by Abel, Bryan & Norman (2002) and Bromm, Coppi & Larson (2002) again being realised. These massive population III stars will go on to form MBHs of comparable masses which in turn means that the first black holes were also likely massive with masses between $10^3$ and $10^5$ M$_{\odot}$.

The Nancy Grace Roman Space Telescope, NASA's next flagship observatory, will redefine deep-field galaxy survey with a field of view two orders of magnitude larger than Hubble and an angular resolution of matching quality. These future deep-wide galaxy surveys necessitate new simulations to forecast their scientific output and to optimise survey strategies. In this work, we present five realizations of 2-deg^2 lightcones, containing a total of >25 million simulated galaxies with -16 < MUV < -25 spanning z ~ 0 to 10. This dataset enables a new set of experiments with the impacts of survey size on the derived galaxy formation and cosmological constraints. The intrinsic and observable galaxy properties are predicted using a well-established, physics-based semi-analytic modelling approach. We provide forecasts for number density, cosmic SFR, field-to-field variance, and angular two-point correlation functions, and demonstrate how the future wide-field surveys will be able to improve these measurements relative to current generation surveys. We also present a comparison between these lightcones and others that have been constructed with empirical models. The mock lightcones are designed to facilitate the exploration of multi-instrument synergies and connecting with current generation instruments and legacy surveys. In addition to Roman, we also provide photometry for a number of other instruments on upcoming facilities, including Euclid and Rubin, as well as the instruments that are part of many legacy surveys. Full object catalogues and data tables for the results presented in this work are made available through a web-based, interactive portal https://www.simonsfoundation.org/semi-analytic-forecasts.

Antimatter cosmic-rays are used to probe new phenomena in physics, including dark matter annihilation. We use the cosmic-ray positron fraction spectrum by the Alpha Magnetic Spectrometer, to search for such an annihilation signal in the Galaxy. We focus on dark matter with mass between 5 and 120 GeV, producing high-energy electrons and positrons. In these cosmic-ray energies the interplay of multiple astrophysical sources and phenomena, makes this search highly sensitive to the underlying astrophysical background assumptions. We use a vast public library of astrophysical models for the cosmic-ray positron fraction background, to derive robust upper limits on the dark matter's annihilation cross section for a number of annihilation channels. This library accounts for different types of cosmic-ray sources and uncertainties on their distribution in space and time. Also, it accounts for uncertainties on those sources' output, their injected into the interstellar medium cosmic-ray spectra and for uncertainties on cosmic-ray propagation. For any given dark matter particle mass and annihilation channel, upper limits on the annihilation cross section are given by bands that stretch a full order of magnitude in its value. Our work provides weaker limits compared to earlier results, that are however robust to all the relevant astrophysical uncertainties. Between 5 and 15 GeV, we find indications for a possible excess flux of cosmic-ray electrons and positrons. That excess is found for most, but not all of our astrophysical background parameter space, and its significance can vary appreciably. Further scrutiny is necessary to improve the understanding of these lower energy cosmic rays. Finally, we note that even if an excess signal is found in these energies, the current background uncertainties do not allow us to accurately deduce its underlying particle properties.

Primordial black holes (PBHs) form from large-amplitude initial density fluctuations and may comprise some or all of the dark matter. If PBHs have an extended mass spectrum, or in mixed PBH-particle dark matter scenarios, the extreme density fluctuations necessary to produce PBHs also lead to the formation of a much greater abundance of dark matter minihaloes with internal densities potentially of order $10^{12}$ M$_\odot$ pc$^{-3}$. We develop an analytic description of the formation of these ultradense haloes and use it to quantitatively compare PBH and halo distributions. PBHs that contribute only a per cent level fraction of the dark matter are accompanied by ultradense haloes that nevertheless comprise an order-unity fraction. This finding significantly alters the predictions of many PBH scenarios, enabling a variety of new observational tests.

The recent detection of a repeating fast radio burst (FRB) in an old globular cluster in M81 challenges traditional FRB formation mechanisms based on magnetic activity in young neutron stars formed recently in core-collapse supernovae. Furthermore, the detection of this repeater in such a nearby galaxy implies a high local universe rate of similar events in globular clusters. Building off the properties inferred from the M81 FRB, we predict the number of FRB sources in nearby ($d \lesssim 20\,$Mpc) galaxies with large globular cluster systems known. Incorporating the uncertain burst energy distribution, we estimate the rate of bursts detectable in these galaxies by radio instruments such as FAST and MeerKat. Of all local galaxies, we find M87 is the best candidate for FRB detections. We predict M87's globular cluster system contains $\mathcal{O}(10)$ FRB sources at present and that a dedicated radio survey (by FAST) of $\mathcal{O}(10)\,$hr has a $90\%$ probability of detecting a globular cluster FRB in M87. The detection of even a handful of additional globular cluster FRBs would provide invaluable constraints on FRB mechanisms and population properties. Previous studies have demonstrated young neutron stars formed following collapse of dynamically-formed massive white dwarf binary mergers may provide the most natural mechanism for these bursts. We explore the white dwarf merger scenario using a suite of $N$-body cluster models, focusing in particular on such mergers in M87 clusters. We describe a number of outstanding features of this scenario that in principle may be testable with an ensemble of observed FRBs in nearby globular clusters.

HE 0107$-$5240 is a hyper metal-poor star with $\rm [Fe/H]=-5.39$. We performed high-res observations with the ESPRESSO spectrograph at the VLT to constrain the kinematical properties of the binary system HE 0107$-$5240 and to probe the binarity of the sample of 8 most metal-poor stars with $\rm [Fe/H]<-4.5$. Radial velocities are obtained by using cross-correlation in the interval 4200$-$4315A, which contains the strong CH band, against a template in an iterative process. A Bayesian method is applied to calculate the orbit by using the ESPRESSO measurements and others from the literature. A chemical analysis has also been performed in HE0107$-$5240 by means of spectral synthesis. Observations of HE 0107$-$5240 spanning more than 3 years show a monotonic decreasing trend in radial velocity at a rate of approximately by 0.5 ms$^{-1}$d$^{-1}$. The period is constrained at $P_{\rm orb} = 13009_{-1370}^{+1496}$d. New more stringent upper-limits have been found for several elements: a)[Sr/Fe] and [Ba/Fe] are lower than $-0.76$ and $+0.2$ respectively, confirming the star is a CEMP-no; b)$A(Li)< 0.5$ is well below the plateau at $A(Li)=1.1$ found in the Lower Red Giant Branch stars, suggesting Li was originally depleted; and c)the isotopic ratio $^{12}$C/$^{13}$C is 87$\pm6$ showing very low $^{13}$C in contrast to what expected from a spinstar progenitor. We confirm that HE 0107$-$5240 is a binary star with a long period of about 13000d ($\sim36$ years).The carbon isotopic ratio excludes the possibility that the companion has gone through the AGB phase and transferred mass to the currently observed star. The binarity of HE 0107$-$5240 implies some of the first generations of low-mass stars form in multiple systems and indicates that the low metallicity does not preclude the formation of binaries. Finally, a solid indication of $v_{ rad}$ variation has been found also in SMSS 1605$-$1443.

We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $\Delta^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that $\Delta^2 (k = 0.36$ $h$ Mpc$^{-1}) \leq 3,496$ mK$^2$ at $z = 10.4$, an improvement by a factor of 2.1 and 2.6 respectively. These limits are mostly consistent with thermal noise over a wide range of $k$ after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration (2022b), we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as $z = 10.4$, ruling out a broad set of so-called "cold reionization" scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result's 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars.

Magnetized radio filaments are found in abundance in the inner few hundred pc of our Galaxy. Progress in understanding this population of filaments has been slow, in part due to a lack of detection elsewhere in the Galaxy or in external galaxies. Recent highly sensitive radio continuum observations of radio galaxies in galaxy clusters have revealed remarkable isolated filamentary structures in the ICM that are linked to radio jets, tails and lobes. The origin of this class of filaments is not understood either. Here, we argue that the underlying physical mechanisms responsible for the creation of the two populations are the same because of their similarity in morphology, spacing between the filaments, aspect ratio, magnetic energy densities to thermal pressure of the medium, and that both populations have undergone synchrotron aging. These similarities provide an opportunity to investigate the physical processes in the ISM and ICM for the first time. We consider that the origin of the filaments in both the GC and ICM is a result of the interaction of a large-scale wind with clouds, or the filaments arise through the stretching and collection of field lines by turbulence in weakly magnetized medium. We examine these ideas toward four radio galaxy filaments associated with four radio galaxies IC 40B, IC 4496, J1333-3141, ESO137-006 and argue that much can be understood in the future by comparing these two populations of filaments.

The formation history of globular clusters (GCs) at redshift $z > 4$ remains an unsolved problem. In this work, we use the cosmological, $N$-body hydrodynamical ``zoom-in'' simulation GigaEris to study the properties and formation of proto-GC candidates in the region surrounding the progenitor of a Milky Way-sized galaxy. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion and allows to resolve systems at the scale of star clusters. We define proto-GC candidate systems as gravitationally bound stellar systems with baryonic mass fraction $F_{\rm b} \geq 0.75$ and stellar velocity dispersion $\sigma_{\star} < 20$ km s$^{-1}$. At $z=4.4$ we identify 9 systems which satisfy our criteria, all of which form between 10 kpc to 30 kpc from the centre of the main host. Their baryonic masses are in the range $10^5$- $10^7$ M$_{\odot}$. By the end of the simulation, they still have a relatively low stellar mass ($M_{\star} \sim 10^4$--$10^5$ M$_{\odot}$) and a metallicity ($-1.8 \lesssim {\rm [Fe/H]} \lesssim -0.8$) similar to the blue Galactic GCs. All of the identified systems except one appear to be associated with gas filaments accreting onto the main galaxy in the circum-galactic region, and formed at $z=5-4$. The exception is the oldest object, which appears to be a stripped compact dwarf galaxy that has interacted with the main halo between $z = 5.8$ and $z=5.2$ and has lost its entire dark matter content due to tidal mass loss.

We continue to explore the relationship between globular cluster total number, $N_{\rm GC}$, and central black hole mass, $M_\bullet$, in spiral galaxies. We present here results for the Sab galaxies NGC 3368, NGC 4736 (M 94) and NGC 4826 (M 64), and the Sm galaxy NGC 4395. The globular cluster (GC) candidate selection is based on the ($u^*$ - $i^\prime$) versus ($i^\prime$ - $K_s$) color-color diagram, and $i^\prime$-band shape parameters. We determine the $M_\bullet$ versus $N_{\rm GC}$ correlation for these spirals, plus NGC 4258, NGC 253, M 104, M 81, M 31, and the Milky Way. We also redetermine the correlation for the elliptical sample in Harris, Poole, & Harris (2014), with updated galaxy types from Sahu et al. 2019b. Additionally, we derive total stellar galaxy mass, $M_\ast$, from its two-slope correlation with $N_{\rm GC}$ (Hudson, Harris, & Harris 2014), and fit $M_\bullet$ versus $M_\ast$ for both spirals and ellipticals. We obtain log $M_\bullet \propto$ (1.01 $\pm$ 0.13) log $N_{\rm GC}$ for ellipticals, and log $M_\bullet \propto$ (1.64 $\pm$ 0.24) log $N_{\rm GC}$ for late type galaxies (LTG). The linear $M_\bullet$ versus $N_{\rm GC}$ correlation in ellipticals could be due to statistical convergence through mergers, but not the much steeper correlation for LTG. However, in the $M_\bullet$ versus total stellar mass ($M_\ast$) parameter space, with $M_\ast$ derived from its correlation with $N_{\rm GC}$, $M_\bullet \propto$ (1.48 $\pm$ 0.18) log $M_\ast$ for ellipticals, and $M_\bullet \propto$ (1.21 $\pm$ 0.16) log $M_\ast$ for LTG. The observed agreement between ellipticals and LTG in this parameter space may imply that black holes and galaxies co-evolve through "calm" accretion, AGN feedback, and other secular processes.

Many studies have looked at the impact of magnetic fields on star formation in molecular clouds and Milky Way like galaxies, concluding that the field suppresses star formation. However, most of these studies are based on fully developed fields that have reached the saturation level, with little work on investigating how the growth phase of a primordial field affects star formation in low metallicity environments. In this paper, we investigate the impact of the growth phase of a primordial field on low metallicity dwarf galaxies. We perform high-resolution AREPO simulations of 5 isolated dwarf galaxies. Two models are hydrodynamical, two start with a primordial B-field of 10^-6 micro G, and one with a saturated B-field of 10^-2 micro G. All models include a non-equilibrium, time-dependent chemical network that includes the effects of gas shielding from the ambient UV field. Sink particles form directly from the gravitational collapse of gas and are treated as star-forming clumps that can accrete gas. We vary the metallicity, UV-field, and cosmic ray ionization rate between 0.01 and 0.10 of solar values. We find that the magnetic field has little impact on the star formation rate, which is in tension with previously published results. We show that an increase in the mass fractions of both molecular hydrogen and cold gas, along with changes in the perpendicular gas velocity dispersion's and the B-field acting in the weak-field model overcomes the expected suppression in star formation.

Continuum reverberation mapping of AGN can provide new insight into the nature and geometry of the accretion flow. Some of the X-rays from the central corona irradiating the disc are absorbed, increasing the local disc temperature. This gives an additional re-processed contribution to the spectral energy distribution (SED) which is lagged and smeared relative to the driving X-ray lightcurve. We directly calculate this reverberation from the accretion disc, creating fully time dependent SEDs for a given X-ray light-curve. We apply this to the intensive monitoring data on Faraill 9 (Hern\'andez Santisteban et al. 2020), and find that it is not possible to produce the observed UV variability by X-ray reprocessing of the observed lightcurve from the disc. Instead, we find that the majority of the variability must be intrinsic to the UV emission process, adding to evidence from changing-look AGN that this region has a structure which is quite unlike a Shakura-Sunyaev disc. We filter out this long timescale variability and find that reprocessing alone is still insufficient to explain even the fast variability in our assumed geometry of a central source illuminating a flat disc. The amplitude of reprocessing can be increased by any vertical structure such as the BLR and/or an inner disc wind, giving a better match. Fundamentally though the model is missing the major contributor to the variabilty, which is from the UV/EUV emission region itself rather than reprocessing.

This White Paper summarises potential key science topics to be achieved with Thai National Radio Telescope (TNRT). The commissioning phase has started in mid 2022. The key science topics consist of "Pulsars and Fast Radio Bursts (FRBs)", "Star Forming Regions (SFRs)", "Galaxy and Active Galactic Nuclei (AGNs)", "Evolved Stars", "Radio Emission of Chemically Peculiar (CP) Stars", and "Geodesy", covering a wide range of observing frequencies in L/C/X/Ku/K/Q/W-bands (1-115 GHz). As a single-dish instrument, TNRT is a perfect tool to explore time domain astronomy with its agile observing systems and flexible operation. Due to its ideal geographical location, TNRT will significantly enhance Very Long Baseline Interferometry (VLBI) arrays, such as East Asian VLBI Network (EAVN), Australia Long Baseline Array (LBA), European VLBI Network (EVN), in particular via providing a unique coverage of the sky resulting in a better complete "uv" coverage, improving synthesized-beam and imaging quality with reducing side-lobes. This document highlights key science topics achievable with TNRT in single-dish mode and in collaboration with VLBI arrays.

Beginning in 2016, the IceCube Neutrino Observatory has sent out alerts in real time containing the information of high-energy ($E \gtrsim 100$~TeV) neutrino candidate events with moderate-to-high ($\gtrsim 30$\%) probability of astrophysical origin. In this work, we use a recent catalog of such alert events, which, in addition to events announced in real-time, includes events that were identified retroactively, and covers the time period of 2011-2020. We also search for additional, lower-energy, neutrinos from the arrival directions of these IceCube alerts. We show how performing such an analysis can constrain the contribution of rare populations of cosmic neutrino sources to the diffuse astrophysical neutrino flux. After searching for neutrino emission coincident with these alert events on various timescales, we find no significant evidence of either minute-scale or day-scale transient neutrino emission or of steady neutrino emission in the direction of these alert events. This study also shows how numerous a population of neutrino sources has to be to account for the complete astrophysical neutrino flux. Assuming sources have the same luminosity, an $E^{-2.5}$ neutrino spectrum and number densities that follow star-formation rates, the population of sources has to be more numerous than $7\times 10^{-9}~\textrm{Mpc}^{-3}$. This number changes to $3\times 10^{-7}~\textrm{Mpc}^{-3}$ if number densities instead have no cosmic evolution.

A collision-induced magnetic reconnection (CMR) mechanism was recently proposed to explain the formation of a filament in the Orion A molecular cloud. In this mechanism, a collision between two clouds with antiparallel magnetic fields produces a dense filament due to the magnetic tension of the reconnected fields. The filament contains fiber-like sub-structures and is confined by a helical magnetic field. To show whether the dense filament is capable of forming stars, we use the \textsc{Arepo} code with sink particles to model star formation following the formation of the CMR-filament. First, the CMR-filament formation is confirmed with \textsc{Arepo}. Second, the filament is able to form a star cluster after it collapses along its main axis. Compared to the control model without magnetic fields, the CMR model shows two distinctive features. First, the CMR-cluster is confined to a factor of $\sim4$ smaller volume. The confinement is due to the combination of the helical field and gravity. Second, the CMR model has a factor of $\sim2$ lower star formation rate. The slower star formation is again due to the surface helical field that hinders gas inflow from larger scales. Mass is only supplied to the accreting cluster through streamers.

Phosphorus (P) is considered to be one of the key elements for life, making it an important element to look for in the abundance analysis of spectra of stellar systems. Yet, there exists only a handful of spectroscopic studies to estimate the P abundances and investigate its trend across a range of metallicities. We have observed full HK band spectra at a spectral resolving power of R=45,000 with IGRINS instrument. Abundances are determined using SME in combination with 1D MARCS stellar atmosphere models. The investigated sample of stars have reliable stellar parameters estimated using optical FIES spectra (GILD; J\"onsson et al. in prep.). In order to determine the P abundances from the 16482.92 Angstrom P line, we take special care of the CO($\nu=7-4$) blend. We determine the C, N, O abundances from atomic carbon and a range of non-blended molecular lines (CO, CN, OH) which are aplenty in the H band region of K giant stars, assuring an appropriate modelling of the blending CO($\nu=7-4$) line. We present [P/Fe] vs [Fe/H] trend for 38 K giant stars in the metallicity range of -1.2 dex $<$ [Fe/H] $<$ 0.4 dex. We find that our trend matches well with the compiled literature sample of prominently dwarf stars and limited number of giant stars. Our trend is found to be higher by $\sim$ 0.05 - 0.1 dex compared to the theoretical chemical evolution trend in Cescutti et al. 2012 resulting from core collapse supernova (type II) of massive stars with the P yields from Kobayashi et al. (2006) arbitrarily increased by a factor of 2.75. Thus the enhancement factor might need to be $\sim$ 0.05 - 0.1 dex higher to match our trend. We also find an empirically determined primary behaviour for phosphorus. Furthermore, the phosphorus abundance is found to be elevated by $\sim$ 0.6 - 0.9 dex in two metal poor s-enriched stars compared to the theoretical chemical evolution trend.

During the Noachian, Mars' crust may have provided a favorable environment for microbial life. The porous brine-saturated regolith would have created a physical space sheltered from UV and cosmic radiations and provided a solvent, while the below-ground temperature and diffusion of a dense reduced atmosphere may have supported simple microbial organisms that consume H2 and CO2 as energy and carbon sources and produce methane as a waste. On Earth, hydrogenotrophic methanogenesis was among the earliest metabolisms but its viability on early Mars has never been quantitatively evaluated. Here we present a probabilistic assessment of Mars' Noachian habitability to H2-based methanogens, and quantify their biological feedback on Mars' atmosphere and climate. We find that subsurface habitability was very likely, and limited mainly by the extent of surface ice coverage. Biomass productivity could have been as high as in early Earth's ocean. However, the predicted atmospheric composition shift caused by methanogenesis would have triggered a global cooling event, ending potential early warm conditions, compromising surface habitability and forcing the biosphere deep into the Martian crust. Spatial projections of our predictions point to lowland sites at low-to-medium latitudes as good candidates to uncover traces of this early life at or near the surface.

It has long been known that the icy shell of Jupiter's moon Europa may drift non-synchronously due to tidal torques. Here we argue that torques applied by the underlying ocean are also important, and can result in retrograde non-synchronous rotation (NSR). We drive an ice shell rotation model with the results of a high-resolution state-of-the-art ocean general circulation model and take into account the viscoelastic deformation of the ice shell. We use the icy shell model results together with observed limits on the icy shell drift speed to constrain icy shell parameters such as its effective viscosity, which is currently unknown by at least four orders of magnitude. Our results suggest at best sluggish ice convection. Depending on the relaxation time scale of the ice shell and on the ocean currents, the icy shell may exhibit negligible drift, constant drift, or oscillatory drift superimposed on random fluctuations. Future spacecraft observations will test these predictions and yield insight into the properties of the ice shell and underlying ocean.

We aim to compare variations in the full-UV dust extinction curve (912-3000 Angstrom), with the HI/H$_2$/total H content along diffuse Milky Way sightlines, to investigate possible connections between ISM conditions and dust properties. We combine an existing sample of 75 UV extinction curves based on IUE and FUSE data, with atomic and molecular column densities measured through UV absorption. The H$_2$ column density data are based on existing Lyman-Werner absorption band models from earlier work on the extinction curves. Literature values for the HI column density were compiled, and improved for 23 stars by fitting a Ly$\alpha$ profile to archived spectra. We discover a strong correlation between the H$_2$ column and the far-UV extinction, and the underlying cause is a linear relationship between H$_2$ and the strength of the far-UV rise feature. This extinction does not scale with HI, and the total H column scales best with $A(V)$ instead. The carrier of the far-UV rise therefore coincides with molecular gas, and further connections are shown by comparing the UV extinction features to the molecular fraction. Variations in the gas-to-extinction ratio $N(\rm{H})/A(V)$ correlate with the UV-to-optical extinction ratio, and we speculate this could be due to coagulation or shattering effects. Based on the H$_2$ temperature, the strongest far-UV rise strengths are found to appear in colder and denser sightlines.

Westerlund 1 (Wd 1) is one of the most relevant star clusters in the Milky Way to study massive star formation, although it is still poorly known. Here, we used photometric and spectroscopic data to model the eclipsing binary W36, showing that its spectral type is O6.5 III + O9.5 IV, hotter and more luminous than thought before. Its distance $d_{\rm W36}$ $=$ 4.03$\pm$0.25 kpc agrees, within the errors, with three recent Gaia-EDR3-based distances reported in Paper I, Beasor & Davies, and by Negueruela's group. However, they follow different approaches to fix the zero-points for red sources such as those in Wd 1 and to select the best approach, we used an accurate modelling of W36. The weighted mean distance of our parallax (Paper I) and binary distances results in $d_{\rm wd1}$ = 4.05 $\pm$0.20 kpc, with an unprecedented accuracy of 5%. We adopted isochrones based on the Geneva code with supersolar abundances to infer the age of W36B as 6.4 $\pm$ 0.7 Myr. This object seems to be part of the prolific star formation burst represented by OB giants and supergiants that occurred at 7.1 $\pm$ 0.5 Myr ago, which coincides with the recently published PMS isochrone with age 7.2 Myr. Other BA-type luminous evolved stars and Yellow Hypergiants spread in the age range of 8--11 Myr. The four Red Supergiants discussed in paper I represent the oldest population of the cluster with an age of 10.7 $\pm$ 1 Myr. The multiple episodes of star formation in Wd 1 are reminiscent of that reported for the R136/30 Dor LMC cluster.

We present here a simple hydrodynamic model based on a sequence of steady states of the inner sub-Keplerian accretion disc to model its spectral states. Correlations between different hydrodynamic steady states are studied with a goal to understand the origin of, e.g., the aperiodic variabilities. The plausible source of corona/outflow close to the central compact object is shown to be a consequence of steady state transition in the underlying accretion flow. We envisage that this phenomenological model can give insight on the influence of environment on the inner sub-Keplerian accretion disc.

We describe the discovery of a solar neighborhood ($d=474$~pc) binary system consisting of a main-sequence sunlike star and a massive non-interacting black hole candidate. We selected this system from the \textit{Gaia} DR3 binary catalog based on its high mass ratio and location close to the main sequence. The spectral energy distribution (SED) of the visible star is well described by a single stellar model, indicating that no contribution from another luminous source is needed to fit the observed photometry. We derive stellar parameters from a high S/N Magellan/MIKE spectrum, classifying the star as a main-sequence star with $T_{\rm eff} = 5972~\rm K$, $\log{g} = 4.54$, and $M = 0.91$~\msun. The spectrum also shows no indication of a second luminous component. We have measured radial velocities of this system with the Automated Planet Finder, Magellan, and Keck over the past three months, which we use to determine the spectroscopic orbit of the binary. We show that the velocity data are consistent with the \textit{Gaia} astrometric orbit and provide independent evidence for a massive dark companion. From a combined fit of the astrometric and spectroscopic data, we derive a companion mass of $11.9^{+2.0}_{-1.6}$\msun. We conclude that this binary system harbors a massive black hole on an eccentric $(e =0.45 \pm 0.02)$, long-period ($185.4 \pm 0.1$ d) orbit. The main-sequence star that orbits this black hole is moderately metal-poor ($\mbox{[Fe/H]} = -0.30$), on a Galactic orbit similar to thin disk stars. Our conclusions are independent of \cite{ElBadry2022Disc}, who recently reported the discovery of the same system, and find a marginally lower companion mass than we do here.

The gravitational collapse of rapidly rotating massive stars can lead to the onset of the low $T/\|W\|$ instability within the central proto-neutron star (PNS), which leaves strong signatures in both the gravitational wave (GW) and neutrino emission. Strong large-scale magnetic fields are usually invoked to explain outstanding stellar explosions of rapidly rotating progenitors, but their impact on the growth of such instability has not yet been cleared. We analyze a series of three-dimensional magnetohydrodynamic models to characterize the effects of different magnetic configurations on the development of the low $T/\|W\|$ and the related multi-messenger features. In the absence of magnetic fields, we observe the growth on dynamical time scales of the low $T/\|W\|$, associated with a strong burst of GW and a correlated modulation of the neutrino emission. However, models with a strong magnetic field show a quenching of the low $T/\|W\|$, due to a flattening of the rotation profile in the first $\sim100$ ms after shock formation caused by the magnetic transport of angular momentum. The associated GW emission is weakened by an order of magnitude, exhibits a broader spectral shape, and has no dominant feature associated with the PNS large-scale oscillation modes. Neutrino luminosities are damped along the equatorial plane due to a more oblate PNS, and the only clear modulation in the signal is due to SASI activity. Finally, magnetized models produce lower luminosities for $\nu_e$ than for $\bar{\nu}_e$, which is connected to a higher concentration of neutron-rich material in the PNS surroundings.

In this paper, we solve the complete coordination problem of robotic fiber positioners using supervisory control theory. In particular, we model positioners and their behavioral specifications as discrete-event systems by the discretization of their motion spaces. We synthesize a coordination supervisor associated with a specific set of positioners. In particular, the coordination supervisor includes the solutions to the complete coordination problem of its corresponding positioners. Then, we use the backtracking forcibility technique of supervisory control theory to present an algorithm based on a completeness condition to solve the coordination problem similar to a reconfiguration problem. We illustrate the functionality of our method using an example.

Gravitational waves from merging binary black holes can be used to shed light on poorly understood aspects of massive binary stellar evolution, such as the evolution of massive stars (including their mass-loss rates), the common envelope phase, and the rate at which massive stars form throughout the cosmic history of the Universe. In this paper we explore the \emph{correlated} impact of these phases on predictions for the merger rate and chirp mass distribution of merging binary black holes, aiming to identify possible degeneracies between model parameters. In many of our models, a large fraction (more than 70% of detectable binary black holes) arise from the chemically homogeneous evolution scenario; these models tend to over-predict the binary black hole merger rate and produce systems which are on average too massive. Our preferred models favour enhanced mass-loss rates for helium rich Wolf--Rayet stars, in tension with recent theoretical and observational developments. We identify correlations between the impact of the mass-loss rates of Wolf--Rayet stars and the metallicity evolution of the Universe on the rates and properties of merging binary black holes. Based on the observed mass distribution, we argue that the $\sim 10\%$ of binary black holes with chirp masses greater than $40$ M$_\odot$ (the maximum predicted by our models) are unlikely to have formed through isolated binary evolution, implying a significant contribution (> 10%) from other formation channels such as dense star clusters or active galactic nuclei. Our models will enable inference on the uncertain parameters governing binary evolution in the near future.

We evaluate the single event tolerance of the X-ray silicon-on-insulator (SOI) pixel sensor named XRPIX, developed for the future X-ray astronomical satellite FORCE. In this work, we measure the cross-section of single event upset (SEU) of the shift register on XRPIX by irradiating heavy ion beams with linear energy transfer (LET) ranging from 0.022 MeV/(mg/cm2) to 68 MeV/(mg/cm2). From the SEU cross-section curve, the saturation cross-section and threshold LET are successfully obtained to be $3.4^{+2.9}_{-0.9}\times 10^{-10}~{\rm cm^2/bit}$ and $7.3^{+1.9}_{-3.5}~{\rm MeV/(mg/cm^2)}$, respectively. Using these values, the SEU rate in orbit is estimated to be $\lesssim$ 0.1 event/year primarily due to the secondary particles induced by cosmic-ray protons. This SEU rate of the shift register on XRPIX is negligible in the FORCE orbit.

The three-dimensional correlation function offers an effective way to summarize the correlation of the large-scale structure even for imaging galaxy surveys. We have applied the projected three-dimensional correlation function, $\xi_{\rm p}$ to measure the Baryonic Acoustic Oscillations (BAO) scale on the first-three years Dark Energy Survey data. The sample consists of about 7 million galaxies in the redshift range $ 0.6 < z_{\rm p } < 1.1 $ over a footprint of $4108 \, \mathrm{deg}^2 $. Our theory modeling includes the impact of realistic true redshift distributions beyond Gaussian photo-$z$ approximation. To increase the signal-to-noise of the measurements, a Gaussian stacking window function is adopted in place of the commonly used top-hat. Using the full sample, $ D_{\rm M}(z_{\rm eff} ) / r_{\rm s} $, the ratio between the comoving angular diameter distance and the sound horizon, is constrained to be $ 19.00 \pm 0.67 $ (top-hat) and $ 19.15 \pm 0.58 $ (Gaussian) at $z_{\rm eff} = 0.835$. The constraint is weaker than the angular correlation $w$ constraint ($18.84 \pm 0.50$) because the BAO signals are heterogeneous across redshift. When a homogeneous BAO-signal sub-sample in the range $ 0.7 < z_{\rm p } < 1.0 $ ($z_{\rm eff} = 0.845$) is considered, $\xi_{\rm p} $ yields $ 19.80 \pm 0.67 $ (top-hat) and $ 19.84 \pm 0.53 $ (Gaussian). The latter is mildly stronger than the $w$ constraint ($19.86 \pm 0.55 $). We find that the $\xi_{\rm p} $ results are more sensitive to photo-$z$ errors than $w$ because $\xi_{\rm p}$ keeps the three-dimensional clustering information causing it to be more prone to photo-$z$ noise. The Gaussian window gives more robust results than the top-hat as the former is designed to suppress the low signal modes. $\xi_{\rm p}$ and the angular statistics such as $w$ have their own pros and cons, and they serve an important crosscheck with each other.

We study the kinematics of emission lines that arise from many physical processes in NGC 6153 based upon deep, spatially-resolved, high resolution spectra acquired with the UVES spectrograph at the ESO VLT. Our most basic finding is that the plasma in NGC 6153 is complex, especially its temperature structure. The kinematics of most emission lines defines a classic expansion law, with the outer part expanding fastest (normal nebular plasma). However, the permitted lines of \ion{O}{1}, \ion{C}{2}, \ion{N}{2}, \ion{O}{2}, and \ion{Ne}{2} present a constant expansion velocity that defines a second kinematic component (additional plasma component). The physical conditions imply two plasma components, with the additional plasma component having lower temperature and higher density. The [\ion{O}{2}] density and the [\ion{N}{2}] temperature are anomalous, but may be understood considering the contribution of recombination to these forbidden lines. The two plasma components have very different temperatures. The normal nebular plasma appears to be have temperature fluctuations in part of its volume (main shell), but only small fluctuations elsewhere. The additional plasma component contains about half of the mass of the N$^{2+}$ and O$^{2+}$ ions, but only $3-5$\% of the mass of H$^+$ ions, so the two plasma components have very different chemical abundances. We estimate abundances of $12+\log(\mathrm O^{2+}/\mathrm H^+)\sim 9.2$\,dex and $\mathrm{He}/\mathrm H\sim 0.13$. Although they are all complications, multiple plasma components, temperature fluctuations, and the contributions of multiple physical processes to a given emission line are all part of the reality in NGC 6153, and should generally be taken into account.

We present a new Australian Square Kilometre Array Pathfinder (ASKAP) sample of 14 radio Supernova Remnant (SNR) candidates in the Large Magellanic Cloud (LMC). This new sample is a significant increase to the known number of older, larger and low surface brightness LMC SNRs. We employ a multi-frequency search for each object and found possible traces of optical and occasionally X-ray emission in several of these 14 SNR candidates. One of these 14 SNR candidates (MCSNR J0522-6543) has multi-frequency properties that strongly indicate a bona fide SNR. We also investigate a sample of 20 previously suggested LMC SNR candidates and confirm the SNR nature of MCSNR J0506-6815. We detect lower surface brightness SNR candidates which were likely formed by a combination of shock waves and strong stellar winds from massive progenitors (and possibly surrounding OB stars). Some of our new SNR candidates are also found in a lower density environments in which SNe type Ia explode inside a previously excavated interstellar medium (ISM).

Planetary climates are strongly affected by planetary orbital parameters such as obliquity, eccentricity, and precession. In exoplanetary systems, exo-terrestrial planets should have various obliquities. High-obliquity planets would have extreme seasonal cycles due to the seasonal change of the distribution of the insolation. Here, we introduce the Non-hydrostatic ICosahedral Atmospheric Model(NICAM), a global cloud-resolving model, to investigate the climate of high-obliquity planets. This model can explicitly simulate a three-dimensional cloud distribution and vertical transports of water vapor. We simulated exo-terrestrial climates with high resolution using the supercomputer FUGAKU. We assumed aqua-planet configurations with 1 bar of air as a background atmosphere, with four different obliquities ($0^{\circ}$, $23.5^{\circ}$, $45^{\circ}$, and $60^{\circ}$). We ran two sets of simulations: 1) low-resolution (~ 220 km-mesh as the standard resolution of a general circulation model for exoplanetary science) with parametrization for cloud formation, and 2) high-resolution (~ 14 km-mesh) with an explicit cloud microphysics scheme. Results suggest that high-resolution simulations with an explicit treatment of cloud microphysics reveal warmer climates due to less low cloud fraction and a large amount of water vapor in the atmosphere. It implies that treatments of cloud-related processes lead to a difference between different resolutions in climatic regimes in cases with high obliquities.

We used the times of occultations and eclipses between the components of the 65803 Didymos binary system observed in its lightcurves from 2003-2021 to estimate the orbital parameters of Dimorphos relative to Didymos. We employed a weighted least-squares approach and a modified Keplerian orbit model in order to accommodate the effects from non-gravitational forces such as Binary YORP that could cause a linear change in mean motion over time. We estimate that the period of the mutual orbit at the epoch 2022 September 26.0 TDB, the day of the DART impact, is $11.9214869 \pm 0.000028$~h ($1\sigma$) and that the mean motion of the orbit is changing at a rate of $(5.0\pm 1.0)\times 10^{-18}$~rad s$^{-2}$ $(1\sigma$). The formal $3\sigma$ uncertainty in orbital phase of Dimorphos during the planned Double Asteroid Redirection Test (DART) mission is $5.4^\circ$. Observations from July to September 2022, a few months to days prior to the DART impact, should provide modest improvements to the orbital phase uncertainty and reduce it to about $4.2^\circ$. These results, generated using a relatively simple model, are consistent with those generated using the more sophisticated model of \citet{scheirich22}, which demonstrates the reliability of our method and adds confidence to these mission-critical results.

We present pre- and post-explosion observations of the Type II-P supernova (SN II-P) 2019mhm located in NGC 6753. Based on optical spectroscopy and photometry, we show that SN 2019mhm exhibits broad lines of hydrogen with a velocity of $-8500\pm200$ km s$^{-1}$ and a $111\pm2$ day extended plateau in its luminosity, typical of the Type II-P subclass. We also fit its late-time bolometric light curve and infer that it initially produced a ${}^{56}$Ni mass of $5.7\substack{+7.0\\-3.0} \times 10^{-3}$ M$_{\odot}$. Using imaging from the Wide Field Planetary Camera 2 on the Hubble Space Telescope obtained 19 years before explosion, we aligned to a post-explosion Wide Field Camera 3 image and demonstrate that there is no detected counterpart to the SN to a limit of $>$24.53~mag in F814W, corresponding to an absolute magnitude limit of $M_{\rm F814W} < -7.7$~mag. Comparing to massive-star evolutionary tracks, we determine that the progenitor star had a maximum zero-age main sequence mass $<$17.5 M$_{\odot}$, consistent with other SN~II-P progenitor stars. SN 2019mhm can be added to the growing population of SNe II-P with both direct constraints on the brightness of their progenitor stars and well-observed SN properties.

In astronomy and cosmology, significant effort is devoted to characterizing and understanding spatial cross-correlations between points - e.g. galaxy positions, high energy neutrino arrival directions, X-ray and AGN sources, and continuous field - e.g. weak lensing and Cosmic Microwave Background (CMB) maps. Recently, we introduced the $k$-nearest neighbor formalism to better characterize the clustering of discrete (point) datasets. Here we extend it to the point-field cross-correlation analysis. It combines $k$NN measurements of the point data set with measurements of the field smoothed on many scales. The resulting statistics are sensitive to all orders in the joint clustering of the points and the field. We demonstrate that this approach, unlike the 2-pt cross-correlation, can measure the statistical dependence of two datasets even when there are no linear (Gaussian) correlations. We further demonstrate that this framework is far more effective than the two-point function in detecting cross-correlations when the continuous field is contaminated by high levels of noise. For a particularly high level of noise, the cross-correlations between halos and the underlying matter field in a cosmological simulation, between $10h^{-1}{\rm Mpc}$ and $30h^{-1}{\rm Mpc}$, is detected at $>5\sigma$ significance using the technique presented here, when the two-point cross-correlation significance is $\sim 1\sigma$. Finally, we show that the $k$NN cross-correlations of halos and the matter field can be well-modeled on quasilinear scales by the Hybrid Effective Field Theory (HEFT) framework, with the same set of bias parameters as are used for the two-point cross-correlations. The substantial improvement in the statistical power of detecting cross-correlations with this method makes it a promising tool for various cosmological applications.

The advent of large-scale stellar spectroscopic surveys naturally leads to the implementation of machine learning techniques to isolate, for example, small sub-samples of potentially interesting stars from the full data set. A recent example is the application of the t-SNE statistical method to $\sim$600,000 stellar spectra from the GALAH survey in order to identify a sample of candidate extremely metal-poor (EMP, [Fe/H] $\leq$ -3) stars. We report the outcome of low-resolution spectroscopic follow-up of 83 GALAH EMP candidates that lack any previous metallicity estimates. Overall, the statistical selection is found to be efficient ($\sim$one-third of the candidates have [Fe/H] $\leq$ -2.75) with low contamination ($<$10% have [Fe/H] $>$ -2), and with a metallicity distribution function that is consistent with previous work. Five stars are found to have [Fe/H] $\leq$ -3.0, one of which is a main sequence turnoff star. Two other stars are revealed as likely CEMP-$s$ stars and a known carbon star is re-discovered. The results indicate that the statistical selection approach employed is valid, and therefore it can be applied to forthcoming even larger stellar spectroscopic surveys with the expectation of similar positive outcomes.

In this paper, we report on the hard X-ray observation of the X-ray pulsar 1E 1145.1-6141 performed with the Nuclear Spectroscopic Telescope Array mission (NuSTAR). The coherent pulsation of the source with a period of $\sim296.653\;\pm\;0.021\;s $ is detected. The source may be in the equilibrium phase, according to the most recent measurements of its pulse period. The pulse profile reveals a mild energy dependence and generally hints at a pencil-beam pattern. The pulse profile have evolved with time. The Pulse fraction is found to depend on energy with a fall in the value at $\sim 32\; keV$. The NuSTAR spectra can be approximated by a composite model with two continuum components, a blackbody emission, cut-off powerlaw, and a discrete component in the form of gaussian to account for the emission line of iron. The estimated absorbed flux of the source is $\sim6\times10^{-10}\;erg\;cm^{-2}\;s^{-1}$ which corresponds to a luminosity of $\sim5\times 10^{36}\;erg\;s^{-1}$. Pulse phase-resolved spectroscopy were performed to understand the evolution of spectral parameters with pulse phase. The estimated blackbody radius is found to be consistent with the size of the theoretical prediction.

This study aims to characterise, for the first time, intranight optical variability (INOV) of low-mass active galactic nuclei (LMAGN) which host a black hole (BH) of mass $M_{BH} \sim 10^6 M_{\odot}$, i.e., even less massive than the Galactic centre black hole Sgr A* and 2-3 orders of magnitude below the supermassive black holes (SMBH, $M_{BH}$ $\sim$ $10^8 - 10^9 M_{\odot}$) which are believed to power quasars. Thus, LMAGN are a crucial subclass of AGN filling the wide gap between SMBH and stellar-mass BHs of Galactic X-ray binaries. We have carried out a 36-session campaign of intranight optical monitoring of a well-defined, representative sample of 12 LMAGNs already detected in X-ray and radio bands. This set of LMAGN is found to exhibit INOV at a level statistically comparable to that observed for blazars (M$_{BH} \gtrsim$ 10$^{8-9}$ M$_{\odot}$) and for the $\gamma$-ray detected Narrow-line Seyfert1 galaxies (M$_{BH}\sim 10^7$ M$_{\odot}$) which, too, are believed to have relativistic jets. This indicates that the blazar-level activity can even be sustained by central engines with black holes near the upper limit for Intermediate Mass Black Holes ($M_{BH}$ $\sim$ $10^3 - 10^6 M_{\odot}$).

Comet C/2014 UN$_{271}$, alternative designation \emph{BB} after its discoverers \emph{Bernardinelli/Bernstein}, and commonly referred to as UN$_{271}$, is an extreme case on two fronts, firstly its solar distance on discovery ($>$ 29 $\si{au}$) and secondly the size of its nucleus (137$\pm$ 15 $\si{km}$). With an aphelion distance of $\sim$33,000 $\si{au}$ (w.r.t. the solar system barycentre) and an orbital period $\sim$2 million $\si{years}$, it is definitely an object from the solar system's \emph{Oort cloud}, and also by a good measure the largest Oort cloud object ever observed. \emph{In situ} observation of UN$_{271}$ would be of considerable scientific importance. Unlike most Oort cloud comets which have been discovered for the first time only as they near the inner solar system, UN$_{271}$ was discovered early enough to provide adequate advanced warning to plan for such a mission. In this paper we describe the various methods for reaching UN$_{271}$ during the period around its perihelion and ecliptic plane passage, with both flyby and rendezvous options; exploiting direct transfers, Jupiter powered gravitational assists (GA) or alternatively a series of GAs of the inner planets. Viable flyby and rendezvous trajectories are found, especially using the NASA Space Launch System (SLS) as the launch vehicle.

We establish a new cosmological-model-independent method to determine the Hubble constant $H_0$ from the localized FRBs and the Hubble parameter measurements and obtain a first such determination $H_0=70.60\pm2.11~\mathrm{km/s/Mpc}$ of about 3.00\% uncertainty with data from the eighteen localized FRBs and nineteen Hubble parameter measurements in the redshift range $0<z\leq0.66$. This value, which is independent of the cosmological model, lies in the middle of the results from the Planck CMB observations and the nearby type Ia supernovae (SN Ia) data. Simulations show that the uncertainty of $H_0$ can be decreased to the level of that from the nearby SN Ia when mock data from 500 localized FRBs with 25 Hubble parameter in the redshift range of $0<z\leq1$ are used. Since localized FRBs are expected to be detected in large quantities, our method will be able to give a reliable and more precise determination of $H_0$ in the very near future, which will help us to figure out the possible origin of the Hubble constant tension.

We study the radio variability of galaxies with and without sources of hydroxyl (OH) megamaser radiation based on the continuum radio measurements conducted in 2019-2022 with the radio telescope RATAN-600 at frequencies of 2.3, 4.7, 8.2, and 11.2 GHz. Presumably, radio continuum emission significantly affects the megamaser radiation brightness, therefore, such a characteristic as the variability of radio emission is important for determining the OHM galaxies parameters. With additional data from the literature, the parameters of radio variability on a time scale up to 30 years were estimated. The median values of the variability index for 48 OHM galaxies are in the range $V_{S}=0.08$-$0.17$, and for 30 galaxies without OH emission they are $V_{S}=0.08$-$0.28$. For some individual galaxies in both samples, flux density variations reach 30-50%. These sources either are commonly associated with AGNs or reveal active star formation. Generally, the variability of luminous infrared galaxies with and without OH megamaser emission is moderate and of the same order of magnitude on long time scales. From estimating the spectral energy distribution parameters in a broad frequency range (from MHz to THz), we determined the spectral index below 50 GHz and the color temperatures of dust components for megamaser and control sample galaxies. At a level of $\rho<0.05$, there are no statistically significant differences in the distribution of these parameters for the two samples, as well there are no statistically significant correlations between the dust color temperatures and the variability index or luminosity in the OH line.

Based on galaxies in the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and dark matter haloes in the dark matter only, cosmological and constrained ELUCID simulation, we investigate the relation between the observed radii of central galaxies with stellar mass $\gtrsim 10^{8} h^{-2}{\rm M}_\odot$ and the virial radii of their host dark matter haloes with virial mass $\gtrsim 10^{10.5} h^{-1}{\rm M}_\odot$, and the dependence of galaxy-halo size relation on the halo spin and concentration. Galaxies in observation are matched to dark matter (sub-)haloes in the ELUCID simulation using a novel neighborhood subhalo abundance matching method. For galaxy 2D half-light radii $R_{50}$, we find that early- and late-type galaxies have the same power-law index 0.55 with $R_{50} \propto R_{\rm vir}^{0.55}$, although early-type galaxies have smaller 2D half-light radii than late-type galaxies at fixed halo virial radii. When converting the 2D half-light radii $R_{50}$ to 3D half-mass radii $r_{1/2}$, both early- and late-type galaxies display similar galaxy-halo size relations with $\log r_{1/2} = 0.55 \log (R_{\rm vir}/210 h^{-1}{\rm kpc}) + 0.39$. We find that the galaxy-halo size ratio $r_{1/2}/ R_{\rm vir}$ decreases with increasing halo mass. At fixed halo mass, there is no significant dependence of galaxy-halo size ratio on the halo spin or concentration.

Context. A universally accepted definition of what a vortex is has not yet been reached. Therefore, we lack an unambiguous and rigorous method for the identification of vortices in fluid flows. Such a method would be necessary to conduct robust statistical studies on vortices in highly dynamical and turbulent systems, such as the solar atmosphere. Aims. We aim to develop an innovative and robust automated methodology for the identification of vortices based on local and global characteristics of the flow. Moreover, the use of a threshold that could potentially prevent the detection of weak vortices in the identification process should be avoided. Methods. We present a new method that combines the rigor of mathematical criteria with the global perspective of morphological techniques. The core of the method consists in the estimation of the center of rotation for every point of the flow that presents some degree of curvature in its neighborhood. For that, we employ the Rortex criterion and combine it with morphological considerations of the velocity field. We then identify coherent vortical structures by clusters of estimated centers of rotation. Results. We demonstrate that the Rortex is a more reliable criterion than are the swirling strength and the vorticity for the extraction of physical information from vortical flows, because it measures the rigid-body rotational part of the flow alone and is not biased by the presence of pure or intrinsic shears. We show that the method performs well on a simplistic test case composed of two Lamb-Oseen vortices. We combine the proposed method with a state of the art clustering algorithm to build an automated vortex identification algorithm. (Abridged)

Black hole X-ray binaries (BHXBs) show rich phenomenology in the spectral and timing properties. We collected the spectral data of 20 BHXBs from the literature across different spectral states. The spectral properties are studied in the forms of the inner disc temperature ($T_{\rm in}$), photon index ($\Gamma$), hot electron temperature ($kT_{\rm e}$), X-ray flux ($F_{\rm X}$) and luminosity ($L_{\rm X}$). We studied various correlations among different spectral parameters to understand the accretion process on a global scale. In the thermal soft states (TSS), we find most of the sources followed $F_{\rm disc} \propto T_{\rm in}^4$ relation. A `V'-shaped correlation is found between $\Gamma$ and total luminosity ($L_{\rm tot}$) in the hard Comptonized state (HCS). The Comptonized luminosity is observed to be correlated with the disc luminosity in the HCS and TSS. No notable correlation is observed in the intermediate state (IMS). The evolution of the inner disc radius ($R_{\rm in}$) is unclear in the HCS and IMS. We also discuss how the hot electron temperature changes with other spectral parameters. We observe that the iron line flux correlates with disc and Comptonized fluxes. The strength of the reprocessed emission is found to vary across spectral states.

Cluster galaxies exhibit substantially lower star formation rates than field galaxies today, but it is conceivable that clusters were sites of more active star formation in the early universe. Herein, we present an interpretation of the star formation history (SFH) of group/cluster galaxies based on the large-scale cosmological hydrodynamic simulation, Horizon-AGN. We find that massive galaxies in general have small values of e-folding timescales of star formation decay (i.e., ``mass quenching'') regardless of their environment, whilst low-mass galaxies exhibit prominent environmental dependence. In massive host halos (i.e., clusters), the e-folding timescales of low-mass galaxies are further decreased if they reside in such halos for a longer period of time. This ``environmental quenching'' trend is consistent with the theoretical expectation from ram pressure stripping. Furthermore, we define a ``transition epoch'' as where cluster galaxies become less star-forming than field galaxies. The transition epoch of group/cluster galaxies varies according to their stellar and host cluster halo masses. Low-mass galaxies in massive clusters show the earliest transition epoch of $\sim 7.6$ Gyr ago in lookback time. However, it decreases to $\sim 5.2$ Gyr for massive galaxies in low-mass clusters. Based on our findings, we can describe cluster galaxy's SFH with regard to the cluster halo-to-stellar mass ratio.

KIC 7955301 is a hierarchical triple system with eclipse timing and depth variations discovered by the Kepler mission. It is composed of a non-eclipsing primary star at the bottom of the red giant branch on a 209-day orbit with a K/G-type main-sequence inner eclipsing binary, orbiting in 15.3 days. This system was noted for the large amplitude of its eclipse timing variations (4 hours), and the clear solar-like oscillations of the red-giant component, including p-modes of degree up to l=3 and mixed l=1 modes. The system is a single-lined spectroscopic triple. We perform a dynamical model by combining the Kepler photometric data, eclipse timing variations, and radial-velocity data obtained at Apache Point (ARCES) and Haute Provence (SOPHIE) observatories. The dynamical mass of the red-giant is determined with a 2% precision at 1.30 (+0.03,-0.02) solar mass. We perform asteroseismic modeling based on the global seismic parameters and on the individual frequencies. Both methods lead to a mass of the red giant that matches the dynamical mass within the uncertainties. Asteroseismology also reveals the rotation rate of the core (15 days), the envelope (150 days), and the inclination (75 deg) of the red giant. Three different approaches lead to an age between 3.3 and 5.8 Gyr, which highlights the difficulty of determining stellar ages despite the exceptional wealth of available information. On short timescales, the inner binary exhibits eclipses with varying depths during a 7.3-year long interval, and no eclipses during the consecutive 11.9 years. This is why Kepler could detect its eclipses, TESS will not, and the future ESA PLATO mission should. Over the long term, the system owes its evolution to the evolution of its individual components. It could end its current smooth evolution by merging by the end of the red giant or the asymptotic giant branch of the primary star.

Aims. We analyze whether the effects of secular perturbations, originating from a substellar companion, on the dust dynamics in a debris disk can be investigated with spatially resolved observations. Methods. We numerically simulated the collisional evolution of narrow and eccentric cold planetesimal belts around a star of spectral type A3V that are secularly perturbed by a companion that orbits either closer to or farther from the star than the belt. Based on the resulting spatial dust distributions, we simulated spatially resolved maps of their surface brightness in the $K$, $N$, and $Q$ bands and at wavelengths of 70$\mu$m and 1300$\mu$m. Results. Assuming a nearby debris disk seen face-on, we find that the brightness distribution varies significantly with observing wavelength, for example between the $N$ and $Q$ band. This can be explained by the varying relative contribution of the emission of the smallest grains near the blowout limit. The orbits of both the small grains that form the halo and the large grains close to the parent belt precess due to the secular perturbations induced by a companion orbiting inward of the belt. The halo, being composed of older grains, trails the belt. The magnitude of the trailing decreases with increasing perturber mass and hence with increasing strength of the perturbation, a trend we recovered in synthetic maps of surface brightness by fitting ellipses to lines of constant brightness. Systems with an outer perturber do not show a uniform halo precession since the orbits of small grains are strongly altered. We identified features of the brightness distributions suitable for distinguishing between systems with a potentially detectable inner or outer perturber, especially with a combined observation with JWST/MIRI in the $Q$ band tracing small grain emission and with ALMA at mm wavelengths tracing the position of the parent planetesimal belt.

This study aims at timing the eclipses of the binary star TU UMi. The times of minima are taken from the literature, from our observations in April 2004 and from TESS observations between 2019 and 2022. The orbital period analysis of the system indicates that there is a cyclic oscillation with an amplitude of 0.0081d and a period of 9.03 yr, accompanied by a continuous decrease at a rate of $dP/dt =-1.12 \times 10^{-7}\,\mathrm{d}\, \mathrm{yr}^{-1}$. We study the secular evolution of the orbital period of the system and the possibility of the existence of a third companion or the magnetic activity cycle of the primary component in the system.

We explore the bounds that can be placed on interactions between cold dark matter and vacuum energy, with equation of state $w=-1$, using state-of-the-art cosmological observations. We consider linear perturbations about a simple background model where the energy transfer per Hubble time, $Q/H$, is a general linear function of the dark matter density, $\rho_c$, and vacuum energy, $V$. We explain the parameter degeneracies found when fitting cosmic microwave background (CMB) anisotropies alone, and show how these are broken by the addition of supernovae data, baryon acoustic oscillations (BAO) and redshift-space distortions (RSD). In particular, care must be taken when relating redshift-space distortions to the growth of structure in the presence of non-zero energy transfer. Interactions in the dark sector can alleviate the tensions between low-redshift measurements of the Hubble parameter, $H_0$, or weak-lensing, $S_8$, and the values inferred from CMB data. However these tensions return when we include constraints from supernova and BAO-RSD datasets. In the general linear interaction model we show that, while it is possible to relax both the Hubble and weak-lensing tensions simultaneously, the reduction in these tensions is modest (reduced to less slightly than $4\sigma$ and $2\sigma$ respectively).

We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We stack ALMA $^{12}$CO(J=2-1) spectra coherently using \hi\ velocities from the VIVA survey to detect faint CO emission out to galactocentric radii $r_{\rm gal} \sim 1.2\,r_{25}$. We determine the scale-lengths for the molecular and stellar components, finding a $\sim$3:5 relation between them, and indicating that CO emission tends to be more centrally concentrated than stellar mass in VERTICO galaxies when compared to field galaxies. While the spatially-resolved $R_{\rm mol}$ shows a decreasing trend with radius similar to field samples, the mean molecular-to-atomic gas ratio within the stellar effective radius $R_{\rm e}$, $R_{\rm mol}(r<R_{\rm e})$, shows a systematic increase with the level of \hi\, truncation and/or asymmetry (perturbation). Analysis of the molecular- and the atomic-to-stellar mass ratios within $R_{\rm e}$, $R^{\rm mol}_{\star}(r<R_{\rm e})$ and $R^{\rm atom}_{\star}(r<R_{\rm e})$, respectively, shows that while VERTICO galaxies do not exhibit significant variations in $R^{\rm mol}_{\star}(r<R_{\rm e})$ with \hi\, perturbation and when compared to galaxies from the field, they do show low $R^{\rm atom}_{\star}(r<R_{\rm e})$ (i.e., perturbation is enough to lower their $R^{\rm atom}_{\star}(r<R_{\rm e})$). We measure a systematic decrease of the SFE$_{\rm mol}$ within $R_{\rm e}$, SFE$_{\rm mol}(r<R_{\rm e})$, with increasingly perturbed \hi. Compared to field galaxies, the increasingly perturbed atomic gas in VERTICO galaxies increases their $R_{\rm mol}$ and decreases the efficiency with which their molecular gas forms stars.

We present a catalogue of over 7000 sources from the GLEAM survey which have significant structure on sub-arcsecond scales at 162MHz. The compact nature of these sources was detected and quantified via their Interplanetary Scintillation (IPS) signature, measured in interferometric images from the Murchison Widefield Array. The advantage of this approach is that all sufficiently compact sources across the survey area are included down to a well-defined flux density limit. The survey is based on $\sim$250$\times$ 10-minute observations, and the area covered is somewhat irregular, but the area within 1hr<RA<11hr; $-10^\circ<$Decl.$<+20^\circ$ is covered entirely, and over 85% of this area has a detection limit for compact structure below 0.2Jy. 7839 sources clearly showing IPS were detected ($>5\sigma$ confidence), with a further 5550 tentative ($>2\sigma$ confidence) detections. Normalised Scintillation Indices (NSI; a measure of the fraction of flux density coming from a compact component) are reported for these sources. Robust and informative upper limits on the NSI are reported for a further 31081 sources. This represents the largest survey of compact sources at radio frequencies ever undertaken.

This seventh paper of the MUse RAdio Loud Emission lines Snapshot (MURALES) project presents the results of the observations obtained with the VLT/MUSE integral field spectrograph of 3C radio sources and discusses the optical spectral properties of the nuclei of 26 objects with 0.3<z<0.82 (median redshift 0.51). At these redshifts the H$\alpha$ and [NII] emission lines are not covered by optical spectra and alternative diagnostic diagrams are needed to separate the different spectroscopic sub-classes. We derive a robust spectroscopic classification into high and low excitation galaxies (HEGs and LEGs) by only using ratios of emission lines in the rest frame UV and blue portion of the spectra. A key result is that FRII/LEGs are found also at the highest level of radio power (up to L$_{178} \sim 2\times 10^{35}$ erg/s/Hz), among the most luminous radio sources in the Universe. Furthermore, their fraction within the FRII RGs population does not strongly depend on radio luminosity. This suggests that the jet properties in powerful FRII radio sources do not depend on the accretion mode and on the structure of the accretion disk as expected if the jet launching process is due to the extraction of the rotational energy of the supermassive black hole. The alternative possibility of recurrent transitions between a LEG and a HEG phase is disfavored based on the variation timescales of the various AGN components.

Asteroseismic scaling relations can provide high-precision measurements of mass and radius for red giant (RG) stars displaying solar-like oscillations. Their accuracy can be validated and potentially improved using independent and accurate observations of mass, radius, effective temperature and metallicity. We seek to achieve this using long period SB2 eclipsing binaries hosting oscillating RGs. We explore KIC 8430105, for which a previous study found significant asteroseismic overestimation of mass and radius when compared with eclipsing binary measurements. We measured dynamical masses and radii for both components to be significantly lower than previously established, increasing the discrepancy between asteroseismic and dynamical measurements. Our dynamical measurements of the RG component were compared to corresponding measurements of mass and radius using asteroseismic scaling relations. Uncorrected scaling relations overestimated the mass of the RG by 26%, the radius by 11%, and the average density by 7%, in agreement with studies for other systems. However, using a theoretical correction to $\Delta \nu$, we managed to obtain an asteroseismic average density that is $1\sigma$ consistent with our dynamical result. We obtained several measurements of $\nu_{max}$ that are not fully consistent. With $\nu_{max} = 76.78 \pm 0.81\mu $Hz, the $\Delta \nu$ correction provided $2 \sigma$ consistent mass and radius for the giant. The age of the system was estimated to be $3.7 \pm 0.4$ Gyr.

We assess the dependence of Earth's disk-integrated mid-infrared thermal emission spectrum on observation geometries and investigate which and how spectral features are impacted by seasonality on Earth. We compiled an exclusive dataset containing 2690 disk-integrated thermal emission spectra for four different full-disk observing geometries (North & South Pole centered and Africa & Pacific centred equatorial views) over four consecutive years. The spectra were derived from 2378 spectral channels in the wavelength range from 3.75 to 15.4 micron (nominal resolution $\approx$ 1200) and were recorded by the Atmospheric Infrared Sounder aboard the Aqua satellite. We learned that there is significant seasonal variability in Earth's thermal emission spectrum, and the strength of spectral features of bio-indicators, such as N2O, CH4, O3 and CO2 depends strongly on both season and viewing geometry. In addition, we found a strong spectral degeneracy with respect to the latter two indicating that multi-epoch measurements and time-dependent signals may be required in order to fully characterize planetary environments. Even for Earth and especially for equatorial views, the variations in flux and strength of absorption features in the disk-integrated data are small and typically $\leq$ 10%. Disentangling these variations from the noise in future exoplanet observations will be a challenge. However, irrespectively of when the planet will be measured (i.e., day or night or season) the results from mid-infrared observations will remain the same to the zeroth order which is an advantage over reflected light observations.

We develop a general framework to compute photometric variations induced by the oblique rotation of a star with an axisymmetric inhomogeneous surface. We apply the framework to compute lightcurves of white dwarfs adopting two simple models of their surface inhomogeneity. Depending on the surface model and the location of the observer, the resulting lightcurve exhibits a departure from a purely sinusoidal curve that are observed for a fraction of white dwarfs. As a specific example, we fit our model to the observed phase-folded lightcurve of a fast-spinning white dwarf ZTF J190132.9+145808.7 (with the rotation period of 419s). We find that both the size and obliquity angle of the spot responsible for the photometric variation are very large, implying an interesting constraint on the surface distribution of the magnetic field on white dwarfs.

Lithosphere is the outer rigid part of terrestrial body, usually consisting of the crust and part of the mantle. Characterizing the physical properties of lithosphere is critical in the investigation to its evolution history. Through the modeling of mass-related loads within lithosphere, physical parameters such as elastic thickness of lithosphere can be inferred by gravity and topography data. At impact basin region, however, the low correlation between topography and gravity makes this model inapplicable. In this work, we proposed a loading model incorporated with the mantle uplift structure commonly formed at impact basin. The resulting deflection caused by this mantle uplift structure is also modeled in the governing equation of the thin elastic shell. Gravity anomaly of the deflected lithosphere is calculated at the surface and the crustal-mantle boundary, then the theoretical gravity admittance and correlation can be compared with observed data. The application of the mantle loading model at four large impact basins on Mars show better fit to the observed admittance and correlation compared with loading model that without this initial mantle plug. Our work suggests that proper modeling of impact-induced load and its resulting deflection is the key to the understanding of physical properties of planetary lithosphere at basin region.

I describe a method to estimate response matrices of Cosmic Microwave Background (CMB) lensing power spectra estimators to the true sky power under realistic conditions. Applicable to all lensing reconstruction pipelines based on quadratic estimators (QE), it uses a small number of Gaussian CMB Monte-Carlos and specially designed QE's in order to obtain sufficiently accurate matrices with little computational effort. This method may be used to improve the modelling of CMB lensing band-powers by incorporating at least some of the non-idealities encountered in CMB lensing reconstruction. These non-idealities always include masking, and often inhomogeneous filtering, either in the harmonic domain or pixel space. I obtain these matrices for Planck latest lensing reconstructions, and then show that the residual couplings induced by masking explain very well the residual multiplicative bias seen on the Planck simulations, removing the need for an empirical correction.

The 2-m aperture Chinese Space Station Telescope (CSST), which observes at wavelengths ranging from 255 to 1000 nm, is expected to start science operations in 2024. An ultra-deep field observation program covering approximately 10 square degrees is proposed with supernovae (SNe) and other transients as one of its primary science drivers. This paper presents the simulated detection results of type Ia supernovae (SNe Ia) and explores the impact of new datasets on the determinations of cosmological parameters. The simulated observations are conducted with an exposure time of 150 s and cadences of 10, 20, and 30 days. The survey mode covering a total of 80 observations but with a random cadence in the range of 4 to 14 days is also explored. Our simulation results indicate that the CSST can detect up to $\sim 1800$ SNe Ia at z $<$ 1.3. The simulated SNe Ia are then used to constrain the cosmological parameters. The constraint on $\Omega_m$ can be improved by 37.5% using the 10-day cadence sample in comparison with the Pantheon sample. A deeper measurement simulation with a 300 s exposure time together with the Pantheon sample improves the current constraints on $\Omega_m$ by 58.3% and $\omega$ by 47.7%. Taking future ground-based SNe Ia surveys into consideration, the constraints on $\omega$ can be improved by 59.1%. The CSST ultra-deep field observation program is expected to discover large amounts of SNe Ia over a broad redshift span and enhance our understanding of the nature of dark energy.

Solar flares and coronal mass ejections (CMEs) cause immediate and adverse effects on the interplanetary space and geospace. The deeper understanding of the mechanisms that produce them and the construction of efficient prediction schemes are of paramount importance. The source regions of flares and CMEs exhibit some common morphological characteristics associated with strongly sheared magnetic polarity inversion lines, indicative of the complex magnetic configurations that store huge amounts of free magnetic energy and helicity. This knowledge is transformed into parameters that can help us distinguish efficiently between quiet, flare-, and CME-productive active regions. Nonetheless, flare and CME prediction still faces a number of challenges. The magnetic field information is constrained at the photosphere and accessed only from one vantage point of observation; the dynamic behavior of active regions is still not fully incorporated into predictions; the stochasticity of flares and CMEs renders their prediction probabilistic. To meet these challenges, new properties have been put forward to describe different aspects of magnetic energy storage mechanisms in active regions and offer the opportunity of parametric studies for over an entire solar cycle. This inventory of predictors now includes information from flow fields, transition region/coronal spectroscopy, data-driven modeling of the coronal magnetic field, as well as parameterizations of dynamic effects from time series. Further work towards these directions may help alleviate the current limitations in observing the magnetic field of higher atmospheric layers. This paper reviews these efforts as well as the importance of transforming new knowledge into more efficient predictors and including new types of data.

The {\sl TESS} space mission has recently demonstrated its great potential to discover new pulsating white dwarf and pre-white dwarf stars, and to detect periodicities with high precision in already known white-dwarf pulsators. We report the discovery of two new pulsating He-rich atmosphere white dwarfs (DBVs) and present a detailed asteroseismological analysis of three already known DBV stars employing observations collected by the {\sl TESS} mission along with ground-based data. We extracted frequencies from the {\sl TESS} light curves of these DBV stars using a standard pre-whitening procedure to derive the potential pulsation frequencies. All the oscillation frequencies that we found are associated with $g$-mode pulsations with periods spanning from $\sim 190$ s to $\sim 936$ s. We find hints of rotation from frequency triplets in some of the targets, including the two new DBVs. For three targets, we find constant period spacings, which allowed us to infer their stellar masses and constrain the harmonic degree $\ell$ of the modes. We also performed period-to-period fit analyses and found an asteroseismological model for three targets, with stellar masses generally compatible with the spectroscopic masses. Obtaining seismological models allowed us to estimate the seismological distances and compare them with the precise astrometric distances measured with {\it Gaia}. We find a good agreement between the seismic and the astrometric distances for three stars (PG~1351+489, EC~20058$-$5234, and EC~04207$-$4748), although for the other two stars (WD~J152738.4$-$50207 and WD~1708$-$871), the discrepancies are substantial. The high-quality data from the {\sl TESS} mission continue to provide important clues to determine the internal structure of pulsating pre-white dwarf and white dwarf stars through the tools of asteroseismology.

In this work, we assess the combined use of Gaia photometry and astrometry with infrared data from CatWISE in improving the identification of extragalactic sources compared to the classification obtained using Gaia data. We evaluate different input feature configurations and prior functions, with the aim of presenting a classification methodology integrating prior knowledge stemming from realistic class distributions in the universe. In our work, we compare different classifiers, namely Gaussian Mixture Models (GMMs), XGBoost and CatBoost, and classify sources into three classes - star, quasar, and galaxy, with the target quasar and galaxy class labels obtained from SDSS16 and the star label from Gaia EDR3. In our approach, we adjust the posterior probabilities to reflect the intrinsic distribution of extragalactic sources in the universe via a prior function. We introduce two priors, a global prior reflecting the overall rarity of quasars and galaxies, and a mixed prior that incorporates in addition the distribution of the these sources as a function of Galactic latitude and magnitude. Our best classification performances, in terms of completeness and purity of the galaxy and quasar classes, are achieved using the mixed prior for sources at high latitudes and in the magnitude range G = 18.5 to 19.5. We apply our identified best-performing classifier to three application datasets from Gaia DR3, and find that the global prior is more conservative in what it considers to be a quasar or a galaxy compared to the mixed prior. In particular, when applied to the pure quasar and galaxy candidates samples, we attain a purity of 97% for quasars and 99.9% for galaxies using the global prior, and purities of 96% and 99% respectively using the mixed prior. We conclude our work by discussing the importance of applying adjusted priors portraying realistic class distributions in the universe.

Protoplanetary disks surrounding young stars are the birth place of planets. Of particular interest are the transition disks with large inner dust cavities of tens of au, hinting at the presence of massive companions. These cavities were first recognized by a deficit in their Spectral Energy Distribution (SED), later confirmed by millimeter interferometry observations. The Atacama Large Millimeter/submillimeter Array (ALMA) has truly revolutionized the field of spatially resolved imaging of protoplanetary disks in both dust and gas, providing important hints for the origin of gaps and cavities. At the same time, new types of substructures have been revealed. Also infrared observations show a large range of substructures both in resolved imaging, interferometry and spectroscopy. Since the last review paper of transition disks (Protostars and Planets VI), a huge amount of data has been taken, which led to numerous new insights in the origin of transition disks. In this review I will summarize the observational efforts from the past decade, compare their insights with the predictions from SED modeling, analyze the properties of the transition disk population and discuss their role in general disk evolution.

Lorentz invariance violation~(LIV) can change the threshold behavior predicted by special relativity and cause threshold anomalies which affect the propagation of cosmic photons. In this work, we focus on the threshold anomaly effect on cosmic photon attenuations by extragalactic background light~(EBL) and discuss how to identify LIV from observations of very high energy~(VHE) photons propagated from long distance in the universe. We point out that the Large High Altitude Air Shower Observatory~(LHAASO), one of the most sensitive gamma-ray detector arrays currently operating at TeV and PeV energies, is an ideal facility for performing such LIV searching.

We present the full data release of 12CO (3-2) High-Resolution Survey (COHRS), which has mapped the inner Galactic plane over the range of 9.5$^{\circ}$ $\le$ l $\le$ 62.3$^{\circ}$ and $|b| \le 0.5^{\circ}$. The COHRS has been carried out using the Heterodyne Array Receiver Program (HARP) on the 15 m James Clerk Maxwell Telescope (JCMT) in Hawaii. The released data are smoothed to have a spatial resolution of 16.6 arcsec and a velocity resolution of 0.635 km/s, achieving a mean root-mean-square of $\sim 0.6$ K on $T_\mathrm{A}^*$. The COHRS data are useful for investigating detailed three-dimensional structures of individual molecular clouds and large-scale structures such as spiral arms in the Galactic plane. Furthermore, data from other available public surveys of different CO isotopologues and transitions with similar angular resolutions to this survey, such as FUGIN, SEDIGISM, and CHIMPS/CHIMPS2, allow studying the physical properties of molecular clouds and comparing their states with each other. In this paper, we report further observations on R2 and improved data reduction since the original COHRS release. We discuss the characteristics of the COHRS data and present integrated-emission images and a position-velocity (PV) map of the region covered. The PV map shows a good match with the spiral-arm traces from the existing CO and HI surveys. We also obtain and compare integrated one-dimensional distributions of 12CO (1-0) and (3-2) and those of star-forming populations to each other.

The size frequency distribution of exoplanet radii between 1 and 4$R_{\oplus}$ is bimodal with peaks at $\sim$1.4 $R_{\oplus}$ and $\sim$2.4 $R_{\oplus}$, and a valley at $\sim$1.8$R_{\oplus}$. This radius valley separates two classes of planets -- usually referred to as "super-Earths" and "mini-Neptunes" -- and its origin remains debated. One model proposes that super-Earths are the outcome of photo-evaporation or core-powered mass-loss stripping the primordial atmospheres of the mini-Neptunes. A contrasting model interprets the radius valley as a dichotomy in the bulk compositions, where super-Earths are rocky planets and mini-Neptunes are water-ice rich worlds. In this work, we test whether the migration model is consistent with the radius valley and how it distinguishes these views. In the migration model, planets migrate towards the disk inner edge forming a chain of planets locked in resonant configurations. After the gas disk dispersal, orbital instabilities "break the chains" and promote late collisions. This model broadly matches the period-ratio and planet-multiplicity distributions of Kepler planets, and accounts for resonant chains such as TRAPPIST-1, Kepler-223, and TOI-178. Here, by combining the outcome of planet formation simulations with compositional mass-radius relationships, and assuming complete loss of primordial H-rich atmospheres in late giant-impacts, we show that the migration model accounts for the exoplanet radius valley and the intra-system uniformity ("peas-in-a-pod") of Kepler planets. Our results suggest that planets with sizes of $\sim$1.4 $R_{\oplus}$ are mostly rocky, whereas those with sizes of $\sim$2.4 $R_{\oplus}$ are mostly water-ice rich worlds. Our results do not support an exclusively rocky composition for the cores of mini-Neptunes.

We report spectral analysis of the persistent black hole X-ray binary, 4U 1957+115, using AstroSat, Swift and NuSTAR observations carried out between 2016-2019. Modelling with a disk emission, thermal comptonization and blurred reflection components revealed that the source was in the high soft state with the disk flux $\sim 87$ % of the total and high energy photon index $\sim2.6$. There is an evidence that either the inner disk radius varied by $\sim25$ % or the colour hardening factor changed by $\sim12$ %. The values of the inner disk radius imply that for a non-spinning black hole, the black hole mass is $<7$ M$_\odot$ and the source is located $>30 $ kpc away. On the other hand, a rapidly spinning black hole would be consistent with the more plausible black hole mass of $< 10$ M$_\odot$ and a source distance of $\sim10$ kpc. Fixing the distance to $10$ kpc and using a relativistic accretion disk model, constrained the black hole mass to 6 M$_\odot$ and inclination angle to 72$^{\circ}$. A positive correlation is detected between the accretion rate and inner radii or equivalently between the accretion rate and colour factor.

Direct collapse black holes (DCBHs) are the leading candidates for the origin of the first supermassive black holes. However, the role of magnetic fields during their formation is still unclear as none of the previous studies has been evolved long enough to assess their impact during the accretion phase. Here, we report the results from a suite of 3D cosmological magneto-hydrodynamic (MHD) simulations which are evolved for 1.6 Myrs comparable to the expected lifetime of supermassive stars (SMSs). Our findings suggest that magnetic fields are rapidly amplified by strong accretion shocks irrespective of the initial magnetic field strength and reach the saturation state. They stabilize the accretion disks and significantly reduce fragmentation by enhancing the Jeans mass in comparison with pure hydrodynamical runs. Although the initial clump masses are larger in MHD runs, the rapid coalescence of clumps in non-MHD cases due to the higher degree of fragmentation results in similar masses. Overall, the central clumps have masses of $\rm 10^5~M_{\odot}$ and the mean mass accretion rates of $\rm \sim 0.1 ~M_{\odot}/yr$ are similar in both MHD and non-MHD cases. The multiplicity of SMSs is significantly reduced in MHD simulations. Such strongly amplified magnetic fields are expected to launch Jets and outflows which may be detected with upcoming radio telescopes.

Conformity denotes the correlation of properties between pairs of galaxies as a function of separation. Correlations between properties such as star formation rate (SFR), stellar mass, and specific star formation rate (sSFR) have implications for the impact of environment upon galaxy formation and evolution. Conformity between primary galaxies and satellites within the same dark matter halo has been well documented in simulations and observations. However, the existence of conformity at greater distances - known as two-halo conformity - remains uncertain. We investigate whether galaxies in the Local Volume to a distance of 4 Mpc show conformity by examining SFR, sSFR, stellar mass, and quenched fraction as a function of physical separation. Making use of the star formation histories of these galaxies, we then extend this analysis back in time to offer the first probe of conformity inside our past light cone. At the present day, we find that the stellar mass or sSFR of a galaxy correlates with the median SFR of neighboring galaxies at a separation of 2 to 3 Mpc. At a lookback time of 1 Gyr, we find a correlation with the quenched fraction of neighboring galaxies, again at 2 to 3 Mpc separation. These signals of conformity likely arise from the differences between the recent star formation histories of Local Group dwarf galaxies and those outside the Local Group. As current and future missions including JWST, Rubin, and Roman expand the sample of Local Volume galaxies, tests of conformity using star formation histories will provide an important tool for exploring spatio-temporal correlations between galaxies.

We show that the very recent detection of the gamma ray burst GRB221009A by the LHAASO experiment -- with its very-high-energy tail extending up to 18 TeV -- provides a hint at the existence of axion-like-particles with mass ${\cal O} (10^{-10}) \, {\rm eV}$ and two-photon coupling ${\cal O} (10^{-11}) \, {\rm GeV}^{- 1}$.

We report on recent progress in the search for dark matter particles with masses from 1 MeV to 1 GeV. Several dark matter candidates in this mass range are expected to generate measurable electronic-recoil signals in direct-detection experiments. We focus on dark matter particles scattering with electrons in semiconductor detectors since they have fundamentally the highest sensitivity due to their low ionization threshold. Charge-coupled device (CCD) silicon detectors are the leading technology, with significant progress expected in the coming years. We present the status of the CCD program and briefly report on other efforts.

Recent gravitational wave observations of neutron star-neutron star and neutron star-black hole binaries appear to indicate that massive neutron stars may not be too uncommon in merging systems. These discoveries have led to an increased interest in the simulation of merging compact binaries involving massive stars. In this manuscript, we present a first set of evolution of massive neutron star binaries using Monte-Carlo radiation transport for the evolution of neutrinos. We study a range of systems, from nearly symmetric binaries that collapse to a black hole before forming a disk or ejecting material, to more asymmetric binaries in which tidal disruption of the lower mass star leads to the production of more interesting post-merger remnants. For the latter type of systems, we additionally study the impact of viscosity on the properties of the outflows, and compare our results to two recent simulations of identical binaries performed with the WhiskyTHC code. We find excellent agreement on the black hole properties, disk mass, and mass and velocity of the outflows, and some minor but noticeable differences in the evolution of the electron fraction when using a subgrid viscosity model. The method used to account for r-process heating in the determination of the outflow properties appears to have a larger impact on our result than those differences between numerical codes. We also take advantage of the use of a Monte-Carlo code to study in more detail the neutrino energy spectrum, and use the simulation with the most ejected material to verify that our newly implemented Lagrangian tracers provide a reasonable sampling of the matter outflows as they leave the computational grid.

The Migdal effect has received much attention from the dark matter direct detection community, in particular due to its power in setting limits on sub-GeV particle dark matter. Currently, there is no experimental confirmation of the Migdal effect through nuclear scattering using Standard Model probes. In this work, we extend existing calculations of the Migdal effect to the case of neutron-nucleus scattering, with a particular focus on neutron scattering angle distributions in silicon. We identify kinematic regimes wherein the assumptions present in current calculations of the Migdal effect hold for neutron scattering, and demonstrate that these include many viable neutron calibration schemes. We then apply this framework to propose an experimental strategy to measure the Migdal effect in cryogenic silicon detectors using the NEXUS facility at Fermilab.

The spatial covariant gravities provide a natural way to including odd-order spatial derivative terms into the gravitational action, which breaks the parity symmetry at gravitational sector. A lot of parity-violating scalar-tensor theories can be mapped to the spatial covariant framework by imposing the unitary gauge. This provides us a general framework for exploring the parity-violating effects in the primordial gravitational waves (PGWs). The main purpose of this paper is to investigate the polarization of PGWs in the spatial covariant gravities and their possible observational effects. To this end, we first construct the approximate analytical solution to the mode function of the PGWs during the slow-roll inflation by using the uniform asymptotic approximation. With the approximate solution, we calculate explicitly the power spectrum and the corresponding circular polarization of the PGWs analytically. It is shown that the new contributions to power spectrum from spatial covariant gravities contain two parts, one from the parity-preserving terms and the other from the parity-violating terms. While the parity-preserving terms can only affect the overall amplitudes of PGWs, the parity-violating terms induce nonzero circular polarization of PGWs, i.e., the left-hand and right-hand polarization modes of GWs have different amplitudes. The observational implications of this nonzero circular polarization is also briefly discussed.

The analysis of the properties of the X-ray radiation emitted from geometrically thin accretion disks around black holes can be a powerful tool to test General Relativity in the strong field regime. This chapter reviews the state-of-the-art of gravity tests with black hole X-ray data. So far, most efforts have been devoted to test the Kerr hypothesis - namely that the spacetime around astrophysical black holes is described by the Kerr solution - and X-ray data can currently provide among the most stringent constraints on possible deviations from the Kerr geometry. As of now, all X-ray analyses are consistent with the predictions of General Relativity.

Gravitational wave (GW) astronomy offers the potential to probe the wave-optics regime of gravitational lensing. Wave optics (WO) effects are relevant at low frequencies, when the wavelength is comparable to the characteristic lensing time delay multiplied by the speed of light, and are thus often negligible for electromagnetic signals. Accurate predictions require computing the conditionally convergent diffraction integral, which involves highly oscillatory integrands and is numerically difficult. We develop and implement several methods to compute lensing predictions in the WO regime valid for general gravitational lenses. First, we derive approximations for high and low frequencies, obtaining explicit expressions for several analytic lens models. Next, we discuss two numerical methods suitable in the intermediate frequency range: 1) Regularized contour flow yields accurate answers in a fraction of a second for a broad range of frequencies. 2) Complex deformation is slower, but requires no knowledge of solutions to the geometric lens equation. Both methods are independent and complement each other. We verify sub-percent accuracy for several lens models, which should be sufficient for applications to GW astronomy in the near future. Apart from modelling lensed GWs, our method will also be applicable to the study of plasma lensing of radio waves and tests of gravity.