Papers for Thursday, Apr 14 2022

Fast radio bursts (FRBs) represent an exciting frontier in the study of gravitational lensing, due to their brightness, extragalactic nature, and the compact, coherent characteristics of their emission. In a companion work [Kader, Leung+2022], we use a novel interferometric method to search for gravitationally lensed FRBs in the time domain using bursts detected by CHIME/FRB. There, we dechannelize and autocorrelate electric field data at a time resolution of 1.25 ns. This enables a search for FRBs whose emission is coherently deflected by gravitational lensing around a foreground compact object such as a primordial black hole (PBH). Here, we use our non-detection of lensed FRBs to place novel constraints on the PBH abundance outside the Local Group. We use a novel two-screen model to take into account decoherence from scattering screens in our constraints. Our constraints are subject to a single astrophysical model parameter -- the effective distance between an FRB source and the scattering screen, for which we adopt a fiducial distance of 1 parsec. We find that coherent FRB lensing is a sensitive probe of sub-solar mass compact objects. Having observed no lenses in $172$ bursts from $114$ independent sightlines through the cosmic web, we constrain the fraction of dark matter made of compact objects, such as PBHs, to be $f \lesssim 0.8$, if their masses are $\sim 10^{-3} M_{\odot}$.

We apply the theory of nonlinear resonance capture to the problem of a black hole binary (BHB) orbiting a supermassive black hole (SMBH) while embedded in the accretion disk of an active galactic nucleus (AGN). If successful, resonance capture can trigger dramatic growth in the BHB eccentricity, with important consequences for the BHB merger timescale, as well as for the gravitational wave (GW) signature of such an eccentric merger. This resonance capture may occur when the orbital period around the SMBH (the "outer binary") and the apsidal precession of the BHB (the "inner binary") are in a 1:1 commensurability. This effect is analogous to the phenomenon of lunar evection resonance in the early Sun-Earth-Moon system, with the distinction that in the present case, the BHB apsidal precession is due to general relativity, rather than rotationally-induced distortion. In contrast to the case of lunar evection, however, the inner binary undergoes orbital decay driven by GW emission, rather than orbital expansion driven by tidal dissipation. This distinction fundamentally alters the three-body dynamics, forbidding resonance capture, and limiting eccentricity growth. However, if the BHB migrates through of a gaseous AGN disk, the change in the outer binary can counterbalance the suppressing effect of BHB decay, permitting evection resonance capture and the production of eccentric BHB mergers. We compute the likelihood of resonance capture assuming an agnostic distribution of parameters for the three bodies involved and for the properties of the AGN disk. We find that intermediate-mass ratio BHBs (involving an intermediate-mass black hole and a stellar-mass black hole) are the most likely to be captured into evection resonance and thus undergo an eccentric merger. We also compute the GW signature of these mergers, showing that they can enter the LISA band while eccentric.

The Radcliffe Wave is a $\sim3$ kpc long coherent gas structure containing most of the star-forming complexes near the Sun. In this Letter we aim to find a Galactic context for the Radcliffe Wave by looking into a possible relationship between the gas structure and the Orion (Local) Arm. We use catalogs of massive stars and young open clusters based on \textit{Gaia} EDR3 astrometry, in conjunction with kiloparsec-scale 3D dust maps, to investigate the Galactic \textit{XY} spatial distributions of gas and young stars. We find a quasi-parallel offset between the luminous blue stars and the Radcliffe Wave, in that massive stars and clusters are found essentially inside and downstream from the Radcliffe Wave. We examine this offset in the context of color gradients observed in the spiral arms of external galaxies, where the interplay between density wave theory, spiral shocks, and triggered star formation has been used to interpret this particular arrangement of gas/dust and OB stars, and outline other potential explanations as well. We hypothesize that the Radcliffe Wave constitutes the gas reservoir of the Orion (Local) Arm, and presents itself as a prime laboratory to study the interface between Galactic structure, the formation of molecular clouds in the Milky Way, and star formation.

The ultra-relativistic outflows powering gamma-ray bursts (GRBs) acquire angular structure through their interaction with stellar material (dynamical ejecta) in long-soft (short-hard) GRBs. They can often be characterized by a compact, nearly uniform narrow core (with half-opening angle $\theta_{c,\{\epsilon,\Gamma\}}$) surrounded by material with energy per unit solid angle ($\epsilon=\epsilon_c\Theta_{\epsilon}^{-a}$, where $\Theta_{\{\epsilon,\Gamma\}}=[1+\theta^2/\theta_{c,\{\epsilon,\Gamma\}}^2]^{1/2}$) and initial specific kinetic energy ($\Gamma_0-1=[\Gamma_c-1]\Theta_\Gamma^{-b}$) declining as power laws. Multi-wavelength afterglow lightcurves of off-axis jets (with viewing angle $\theta_{obs} >\theta_c$) offer robust ways to constrain a, b and the external density radial profile ($\rho\propto R^{-k}$), even while other burst parameters may remain highly degenerate. We extend our work BGG2020 on such afterglows to include more realistic angular structure profiles derived from three-dimensional hydrodynamic simulations of GRBs (addressing also jets with shallow angular energy profiles, whose emission exhibits unique evolution). We present afterglow lightcurves based on our parameterized power-law jet angular profiles for different viewing angles $\theta_{obs}$ and k=0,1,2. We identify a unique evolutionary power-law phase of the characteristic synchrotron frequencies ($\nu_m$ and $\nu_c$) that manifests when the lightcurve is dominated by emission sensitive to the angular structure of the outflow. We calculate the criterion for obtaining single or double peaked light curves in the general case when $\theta_{c,\Gamma}\neq\theta_{c,\epsilon}$. Following our earlier work, we emphasize how the shape of the lightcurve and the temporal evolution of $\nu_m$ and $\nu_c$ can be used to constrain the outflow structure and potentially distinguish between magnetic and hydrodynamic jets.

Solar flares and coronal mass ejections are the largest energy release phenomena in the current solar system. They cause drastic enhancements of electromagnetic waves of various wavelengths and sometimes eject coronal material into the interplanetary space, disturbing the magnetic surroundings of orbiting planets including the Earth. It is generally accepted that solar flares are a phenomenon in which magnetic energy stored in the solar atmosphere above an active region is suddenly released through magnetic reconnection. Therefore, to elucidate the nature of solar flares, it is critical to estimate the complexity of the magnetic field and track its evolution. Magnetic helicity, a measure of the twist of coronal magnetic structures, is thus used to quantify and characterize the complexity of flare-productive active regions. This chapter provides an overview of solar flares and discusses how the different concepts of magnetic helicity are used to understand and predict solar flares.

We present methods to (i) graphically identify robust redshifts using emission lines in the (sub)mm regime, (ii) evaluate the capabilities of different (sub)mm practices for measuring spectroscopic redshifts, and (iii) optimise future (sub)mm observations towards increasing the fraction of robust redshifts. Using this publicly-available code (https://github.com/tjlcbakx/redshift-search-graphs), we discuss scenarios where robust redshifts can be identified using both single- and multiple-line detections, as well as scenarios where the redshift remains ambiguous, even after the detection of multiple lines. Using the redshift distribution of (sub)mm samples, we quantify the efficiencies of various practices for measuring spectroscopic redshifts, including interferometers, as well as existing and future instruments specifically designed for redshift searches. Finally, we provide a method to optimise the observation strategy for future (sub)mm spectroscopic redshift searches with the Atacama Large Millimetre/submillimetre Array, where 2 mm proves indispensable for robust redshifts in the z = 2 - 4 region.

It is believed that millisecond pulsars attain their fast spins by accreting matter and angular momentum from companion stars. Theoretical modelling of the accretion process suggests a spin-up line in the period-period derivative ($P$-$\dot{P}$) diagram of millisecond pulsars, which plays an important role in population studies of radio millisecond pulsars and accreting neutron stars in X-ray binaries. Here we present observational evidence for such a spin-up line using a sample of 46 radio millisecond pulsars with $P< 10$ ms and intrinsic spin-down rate $\gtrsim 10^{-20}$. We also find that PSRs J1823$-$3021A and J1824$-$2452A, located near the classic spin-up line, are consistent with the broad population of millisecond pulsars. Finally, we show that our approach of Bayesian inference can probe accretion physics, allowing constraints to be placed on the accretion rate and the disk-magnetosphere interaction.

We report an Atacama Large Millimeter/submillimeter Array 0.88 mm (Band 7) continuum detection of the accretion disk around SR 12 c, an $\sim$11 $M_{\rm Jup}$ planetary-mass companion (PMC) orbiting its host binary at 980 au. This is the first submillimeter detection of a circumplanetary disk around a wide PMC. The disk has a flux density of $127 \pm14~\mu$Jy and is not resolved by the $\sim$0.1" beam, so the dust disk radius is likely less than 5 au and can be much smaller if the dust continuum is optically thick. If, however, the dust emission is optically thin, then the SR 12 c disk has a comparable dust mass to the circumplanetary disk around PDS 70 c but is about five times lower than that of the $\sim$12 $M_{\rm Jup}$ free-floating OTS 44. This suggests that disks around bound and unbound planetary-mass objects can span a wide range of masses. The gas mass estimated with an accretion rate of $10^{-11}~M_\odot$ yr$^{-1}$ implies a gas-to-dust ratio higher than 100. If cloud absorption is not significant, a nondetection of ${}^{12}$CO(3-2) implies a compact gas disk around SR 12 c. Future sensitive observations may detect more PMC disks at 0.88 mm flux densities of $\lesssim$100 $\mu$Jy.

The gravitational field of compact objects, such as primordial black holes, can create multiple images of background sources. For transients such as fast radio bursts (FRBs), these multiple images can be resolved in the time domain. Under certain circumstances, these images not only have similar burst morphologies but are also phase-coherent at the electric field level. With a novel dechannelization algorithm and a matched filtering technique, we search for repeated copies of the same electric field waveform in observations of FRBs detected by the FRB backend of the Canadian Hydrogen Mapping Intensity Experiment (CHIME). An interference fringe from a coherent gravitational lensing signal will appear in the time-lag domain as a statistically-significant peak in the time-lag autocorrelation function. We calibrate our statistical significance using telescope data containing no FRB signal. Our dataset consists of $\sim$100-ms long recordings of voltage data from 172 FRB events, dechannelized to 1.25-ns time resolution. This coherent search algorithm allows us to search for gravitational lensing signatures from compact objects in the mass range of $10^{-4}-10^{4} ~\mathrm{M_{\odot}}$. After ruling out an anomalous candidate due to diffractive scintillation, we find no significant detections of gravitational lensing in the 172 FRB events that have been analyzed. In a companion work [Leung, Kader+2022], we interpret the constraints on dark matter from this search.

The Transiting Exoplanet Survey Satellite (TESS) has an exceptionally large plate scale of 21"/px, causing most TESS light curves to record the blended light of multiple stars. This creates a danger of misattributing variability observed by TESS to the wrong source, which would invalidate any analysis. We develop a method that can localize the origin of variability on the sky to better than one fifth of a pixel. Given measured frequencies of observed variability (e.g., from periodogram analysis), we show that the corresponding best-fit sinusoid amplitudes to raw light curves extracted from each pixel are distributed the same as light from the variable source. The primary assumption of this method is that other nearby stars are not variable at the same frequencies. Essentially, we are using the high frequency resolution of TESS to overcome limitations from its low spatial resolution. We have implemented our method in an open source Python package, TESS Localize (github.com/Higgins00/TESS-Localize), that determines the location of a variable source on the sky given TESS pixel data and a set of observed frequencies of variability. Our method utilizes the TESS Pixel Response Function models, and we characterize systematics in the residuals of fitting these models to data. Given the ubiquity of source blending in TESS light curves, verifying the source of observed variability should be a standard step in TESS analyses.

Motivated by the desire for wide-field images with well-defined statistical properties for 21cm cosmology, we implement an optimal mapping pipeline that computes a maximum likelihood estimator for the sky using the interferometric measurement equation. We demonstrate this direct optimal mapping with data from the Hydrogen Epoch of Reionization (HERA) Phase I observations. After validating the pipeline with simulated data, we develop a maximum likelihood figure-of-merit for comparing four sky models at 166MHz with a bandwidth of 100kHz. The HERA data agree with the GLEAM catalogs to $<$10%. After subtracting the GLEAM point sources, the HERA data discriminate between the different continuum sky models, providing most support for the model from Byrne et al. 2021. We report the computation cost for mapping the HERA Phase I data and project the computation for the HERA 320-antenna data, both are feasible with a modern server. The algorithm is broadly applicable to other interferometers, and is valid for wide-field and non-coplanar arrays.

We present new results of a 12CO(J=1-0) imaging survey using the Atacama Compact Array (ACA) for 31 HI detected galaxies in the IC 1459 and NGC 4636 groups. This is the first CO imaging survey for loose galaxy groups. We obtained well-resolved CO data (~0.7-1.5 kpc) for a total of 16 galaxies in two environments. By comparing our ACA CO data with the HI and UV data, we probe the impacts of the group environment on the cold gas components (CO and HI gas) and star formation activity. We find that CO and/or HI morphologies are disturbed in our group members, some of which show highly asymmetric CO distributions (e.g., IC 5264, NGC 7421, and NGC 7418). In comparison with isolated galaxies in the xCOLD GASS sample, our group galaxies tend to have low star formation rates and low H2 gas fractions. Our findings suggest that the group environment can change the distribution of cold gas components, including the molecular gas, and star formation properties of galaxies. This is supporting evidence that pre-processing in the group-like environment can play an important role in galaxy evolution.

Line-of-sight extinction estimates to well-studied young T Tauri and Herbig Ae/Be stars are based on many different measurements and analysis methods. This has resulted in wide scatter among the published $A_V$ values for the same star. In this work, we discuss the challenges in measuring extinction to actively accreting and especially outbursting young stellar objects (YSOs), and explore a method not previously applied to young stars, utilizing diffuse interstellar bands (DIBs). In early-type stars, narrow correlations exist between DIB equivalent widths and the column density of interstellar material, and therefore the line-of-sight extinction. Here, we measure equivalent widths of the 5780 \AA\ and 6614 \AA\ DIB features in a sample of actively accreting YSOs, and apply established DIB--reddening calibrations to derive reddening, and subsequently extinction. We also compare the DIBs-inferred optical line-of-sight extinction values with previous extinction estimates for our sample stars.

Accretion and ejection of matter in active galactic nuclei (AGN) are tightly connected phenomena and represent fundamental mechanisms regulating the growth of the central supermassive black hole and the evolution of the host galaxy. However, the exact physical processes involved are not yet fully understood. We present a high-resolution spectral analysis of a simultaneous \xmm\ and \nustar\ observation of the narrow line Seyfert 1 (NLS1) AGN 1H 1934-063, during which the X-ray flux dropped by a factor of $\sim6$ and subsequently recovered within 140 kiloseconds. By means of the time-resolved and flux-resolved X-ray spectroscopy, we discover a potentially variable warm absorber and a relatively stable ultra-fast outflow (UFO, $v_\mathrm{UFO}\sim-0.075\,c$) with a mild ionization state ($\log(\xi/\mathrm{erg\,cm\,s^{-1})}\sim1.6$). The detected emission lines (especially a strong and broad feature around 1\,keV) are of unknown origin and cannot be explained with emission from plasmas in photo- or collisional-ionization equilibrium. Such emission lines could be well described by a strongly blueshifted ($z\sim-0.3$) secondary reflection off the base of the equatorial outflows, which may reveal the link between the reprocessing of the inner accretion flow photons and the ejection. However, this scenario although being very promising is only tentative and will be tested with future observations.

We use results of a survey for low-surface-gravity stars in Galactic (and LMC) globular clusters to show that "yellow" post-asymptotic-giant-branch (yPAGB) stars are likely to be excellent extragalactic standard candles, capable of producing distances to early-type galaxies that are accurate to a few percent. We show that the mean bolometric magnitude of the 10 known yPAGB stars in globular clusters is <Mbol> = -3.38 +/- 0.03, a value that is ~0.2 mag brighter than that predicted from the latest post-horizontal-branch evolutionary tracks. More importantly, we show that the observed dispersion in the distribution is only 0.10 mag, i.e., smaller than the scatter for individual Cepheids. We describe the physics that can produce such a small dispersion, and show that, if one restricts surveys to the color range 0 < (B-V)0 < 0.5, then samples of non-variable yPAGB stars can be identified quite easily with a minimum of contamination. The bright absolute V magnitudes of these stars (<Mv> = -3.37) make them, by far, the visually brightest objects in old stellar populations and ideal Population II standard candles for measurements out to ~10 Mpc with current instrumentation. A Hubble Space Telescope survey in the halos of galaxies in the M81 and Sculptor groups could therefore serve as an effective cross-check on both the Cepheid and TRGB distance scales.

Old, digitized astronomical images taken before the human spacefaring age offer a unique view of the sky devoid of known artificial satellites. In this paper, we have carried out the first optical searches ever for non-terrestrial artifacts near the Earth following the method proposed in Villarroel et al. (2022). We use images contained in the First Palomar Sky Survey to search for simultaneous (during a plate exposure time) transients that in addition to being point-like, are aligned. We provide a shortlist of the most promising candidates of aligned transients, that must be examined with the help of a microscope to separate celestial sources from plate defects with coincidentally star-like brightness profiles. We further explore one possible, but not unique, interpretation in terms of fast reflections off high-albedo objects in geosynchronous orbits around Earth. If a future study rules out each multiple transient candidate, the estimated surface density becomes an upper limit of $<10^{-9}$ objects km$^{-2}$ non-terrestrial artifacts in geosynchronous orbits around Earth. Finally, we conclude that observations and analysis of multiple, simultaneously appearing and vanishing light sources on the sky merit serious further attention, regardless of their origin.

We present a combined radio/X-ray study of six massive galaxy clusters, aimed at determining the potential for heating of the intra-cluster medium (ICM) by non-central radio galaxies. Since X-ray cavities associated with the radio lobes of non-central galaxies are generally not detectable, we use Giant Metrewave Radio Telescope 610~MHz observations to identify jet sources and estimate their size, and Chandra data to estimate the pressure of the surrounding ICM. In the radio, we detect 4.5% of galaxies above the spectroscopic survey limit (M*K+2.0) of the Arizona cluster redshift survey (ACReS) which covers five of our six clusters. Approximately one tenth of these are extended radio sources. Using star formation (SF) rates determined from mid-infrared data, we estimate the expected contribution to radio luminosity from the stellar population of each galaxy, and find that most of the unresolved or poorly-resolved radio sources are likely star formation dominated. The relatively low frequency and good spatial resolution of our radio data allows us to trace star formation emission down to galaxies of stellar mass ~10^9.5 Msol. We estimate the enthalpy of the (AGN dominated) jet/lobe and tailed sources, and place limits on the energy available from unresolved radio jets. We find jet powers in the range ~10^43-10^46 erg/s, comparable to those of brightest cluster galaxies. Our results suggest that while cluster-central sources are the dominant factor balancing ICM cooling over the long term, non-central sources may have a significant impact, and that further investigation is possible and warranted.

Observational evidence for the accelerated expansion of the universe requires dark energy for its explanation if general relativity is an accurate model of gravity. However, dark energy is a mysterious quantity and we do not know much about its nature so understanding dark energy is an exciting scientific challenge. Cosmological dark energy models are fairly well tested in the low and high redshift parts of the universe. The highest of the low redshift, $z\sim2.3$, region is probed by baryon acoustic oscillation (BAO) measurements and the only high redshift probe is the cosmic microwave background anisotropy which probes the $z\sim1100$ part of redshift space. In the intermediate redshift range $2.3 < z < 1100$ there are only a handful of observational probes and cosmological models are poorly tested in this region. In this thesis we constrain three pairs of general relativistic cosmological dark energy models using observational data which reach beyond the current BAO limit. We use quasar X-ray and UV flux measurements, the current version of these data span $0.009 \leq z \leq 7.5413$. We have discovered that most of these data cannot be standardized using the proposed method. However, the lower redshift part, $z \lesssim 1.5-1.7$, of these data are standardizable and can be used to derive lower-$z$ cosmological constraints. Another data set we use are gamma-ray burst measurements which span $0.3399 \leq z \leq 8.2$. Cosmological constraints derived from these data are significantly weaker than, but consistent with, those obtained from better-established cosmological probes. We also study and standardize 78 reverberation-measured Mg II time-lag quasars in the redshift range $0.0033 \leq z \leq 1.89$ by using their radius-luminosity relation. We also study 118 reverberation-measured H$\beta$ time-lag quasars which span $0.0023 \leq z \leq 0.89$.

Black widows (BWs) are a type of eclipsing millisecond pulsars (MSPs) with companion masses $\lesssim 0.05\,\rm M_\odot$, which can be used to study the accretion history and the radiation of pulsars, as well as the origin of isolated MSPs. Recent observations indicate that there are two sub-types of BWs. One is the BWs with companion masses $M_2$ $\lesssim 0.01\,\rm M_\odot$, whereas another with $M_2$ $\sim$ $0.01-0.05\,\rm M_\odot$. However, the origin of the former is still highly uncertain. In this paper, we investigated the formation of BWs with $M_2$ $\lesssim 0.01\,\rm M_\odot$ through ultra-compact X-ray binaries (UCXBs) with He star companions, in which a neutron star (NS) accretes material from a He star through Roche-lobe overflow. By considering different He star masses and evaporation efficiencies with the stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA), we evolved a series of NS+He star systems that can undergo UCXB stage. We found that this channel can explain the formation of BWs with $M_2$ $\lesssim 0.01\,\rm M_\odot$ within the Hubble time, especially three widely studied BWs, i.e. PSRs J1719-1438, J2322-2650 and J1311-3430. We also found that X-ray irradiation feedback does not affect the evolutionary tracks of evaporation process. Our simulations indicate that the origin of BWs with $M_2$ $\lesssim 0.01\,\rm M_\odot$ is different with another sub-type of BWs, and that the UCXB channel with He star companions are the potential progenitors of isolated MSPs. In addition, the present work suggests that the BWs with $M_2$ $\lesssim 0.01\,\rm M_\odot$ may not be produced by redback systems.

Terzan 5 is a rich globular cluster within the galactic bulge that contains 39 known millisecond pulsars, the largest known population of any globular cluster. The Terzan 5 pulsars are faint, so that individual observations of most of the pulsars have too little signal-to-noise (S/N) to measure reliable flux density or polarization information. We combined over 5.2\,days of archival data, at each of 1.5\,GHz and 2.0\,GHz, taken with the Green Bank Telescope over the past 11\,years. We created high S/N profiles for 32 of the pulsars and determined precise rotation measures (RMs) for 28 of them. We used the RMs, and the known pulsar positions and dispersion measures (DMs), to map the projected parallel component of the Galactic magnetic field toward the cluster. The $\langle B_{||}\rangle$ shows a rough gradient of $\sim$6\,nG/arcsec ($\sim$160\,nG/parsec), or fractionally, a change of $\sim$20$\%$ in the right ascension direction across the cluster, implying Galactic magnetic field variability at sub-parsec scales. We also measured average flux densities $S_\nu$ for the pulsars, ranging from $\sim$10\,$\mu$Jy to $\sim$2\,mJy, and an average spectral index $\alpha = -1.35$, where $S_\nu \propto \nu^{\alpha})$. This spectral index is flatter than most known pulsars, likely a selection effect due to the high frequencies used in pulsar searches to mitigate dispersion and scattering. The inferred pulsar luminosity function is roughly power-law, with slope $(d\log N)/(d\log L) = -1$ at the high-luminosity end. At the low-luminosity end, there are incompleteness effects implying that Terzan 5 contains many more pulsars to be found.

We investigated the detection and localization of binary neutron star (BNS) populations with decihertz gravitational-wave observatories in a realistic detecting strategy, including real-time observations and early warnings. Assuming 4 years' operation of B-DECIGO, we found that the detected BNSs can be divided into three categories: (a) sources that merge within 1 year, which could be localized with an uncertainty of $\Delta\Omega \sim 10^{0}$ deg$^2$; (b) sources that merge in 1-4 years, which take up three quarters of the total events and yield the most precise angular resolution with $\Delta \Omega\sim 10^{-2}$ deg$^2$ and time-of-merger accuracy with $\Delta t_c\sim 10^{-1}$ s; and (c) sources that do not merge during the 4-yr mission window, which enable possible early warnings, with $\Delta \Omega\sim 10^{-1}$ deg$^2$ and $\Delta t_c\sim 10^{0}$ s. Furthermore, we compared the pros and cons of B-DECIGO with the Einstein Telescope, and explored the prospects of detections using 3 other decihertz observatories and 4 BNS population models. In realistic observing scenarios, we found that decihertz detectors could even provide early-warning alerts to a source decades before its merger while their localizations are still more accurate than ground-based facilities. Finally we found a decrease of events when considering the confusion noise, but this could be partially solved by a proper noise subtraction.

We compare the 230 GHz near-horizon emission from Sagittarius A* to simulations representing three classes of accretion flows. Using the structure function to capture the variability statistics of the light curve, we find a noticeable discrepancy between the observations and models based on torus-fed accretion disks, whether those disks bring in a small or large amount of net magnetic flux. On the other hand, the simulations that are fed more realistically by stellar winds match the observed structure function very well. We describe the differences between models, arguing that feeding by stellar winds may be a critical component in constructing theoretical models for accretion in the Galactic Center.

During the transition phase from prestellar to protostellar cloud core, one or several protostars can form within a single gas core. The detailed physical processes of this transition, however, still remain unclear. We present 1.3 mm dust continuum and molecular line observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward 43 protostellar cores in the Orion molecular cloud complex ({\lambda}Orionis, Orion B, and Orion A) with an angular resolution of \sim 0. 35 (\sim 140 au). In total, we detect 13 binary/multiple systems. We derive an overall multiplicity frequency (MF) of 28% \pm 4% and a companion star fraction (CSF) of 51% \pm 6%, over a separation range of 300-8900 au. The median separation of companions is about 2100 au. The occurrence of stellar multiplicity may depend on the physical characteristics of the dense cores. Notably, those containing binary/multiple systems tend to show higher gas density and Mach number than cores forming single stars. The integral-shaped filament (ISF) of Orion A GMC, which has the highest gas density and hosts high-mass star formation in its central region (Orion Nebula Cluster or ONC), shows the highest MF and CSF among Orion GMCs. In contrast, the {\lambda} Orionis Giant Molecular Cloud (GMC) has a lower MF and CSF than the Orion B and Orion A GMCs, indicating that feedback from Hii regions may suppress the formation of multiple systems. We also find that the protostars comprising binary/multiple systems are usually at different evolutionary stages.

In this paper, we report on 606 S-type stars identified from Data Release 9 of the LAMOST medium-resolution spectroscopic (MRS) survey, and 539 of them are reported for the first time. The discovery of these stars is a three-step process, i.e., selecting with the ZrO band indices greater than 0.25, excluding non-S-type stars with the iterative Support Vector Machine method, and finally retaining stars with absolute bolometric magnitude larger than -7.1. The 606 stars are consistent with the distribution of known S-type stars in the color-magnitude diagram. We estimated the C/Os using the [C/Fe] and [O/Fe] provided by APOGEE and the MARCS model for S-type stars, respectively, and the results of the two methods show that C/Os of all stars are larger than 0.5. Both the locations on the color-magnitude diagram and C/Os further verify the nature of our S-type sample. Investigating the effect of TiO and atmospheric parameters on ZrO with the sample, we found that log g has a more significant impact on ZrO than Teff and [Fe/H], and both TiO and log g may negatively correlate with ZrO. According to the criterion of Tian et al. (2020), a total of 238 binary candidates were found by the zero-point-calibrated radial velocities from the officially released catalog of LAMOST MRS and the catalog of Zhang et al. (2021). A catalog of these 606 S-type stars is available from the following link https://doi.org/10.12149/101097.

The dynamics of the partially ionized solar atmosphere is controlled by the frequent collision and charge exchange between the predominant neutral Hydrogen atoms and charged ions. At signal frequencies below or of the order of either of the collision or charge exchange frequencies the magnetic stress is {\it felt} by both the charged and neutral particles simultaneously. The resulting neutral-mass loading of the ions leads to the rescaling of the effective ion-cyclotron frequency-it becomes the Hall frequency, and the resultant effective Larmor radius becomes of the order of few kms. Thus the finite Larmor radius (FLR) effect which manifests as the ion and neutral pressure stress tensors operates over macroscopic scales. Whereas parallel and perpendicular (with respect to the magnetic field) viscous momentum transport competes with the Ohm and Hall diffusion of the magnetic field in the photosphere-chromosphre, the gyroviscous effect becomes important only in the transition region between the chromosphere and corona, where it competes with the ambipolar diffusion. The wave propagation in the gyroviscous effect dominated medium depends on the plasma $\beta$ (a ratio of the thermal and magnetic energies). The abundance of free energy makes gyro waves unstable with the onset condition exactly opposite of the Hall instability. However, the maximum growth rate is identical to the Hall instability. For a flow gradient $\sim 0.1 \,\mbox{s}^{-1}$ the instability growth time is one minute. Thus, the transition region may become subject to this fast growing, gyroviscous instability.

The Boltzmann-Gibbs thermodynamic equilibrium state of charged particles pitch-angle scattered by weak plasma waves is discussed. Degrees of freedom of these waves play a fundamental role in constructing the grand canonical ensemble. Via the gyro-resonance condition, fast particles have an inverse break power-law spectrum for $ \varepsilon -\mu \ll T $, where $ \varepsilon $ is the particle energy, $ \mu $ is the chemical potential, $ T $ is the temperature. The break energies are the rest energy and $ -\mu $. For $ \varepsilon \ll -\mu \ll T $, the energy spectral index $ \alpha $ is $ \delta /2+1 $ and $ \delta +1 $ for non- and ultra-relativistic particles, respectively, with $ \delta $ an effective fractal dimension of background magnetic field lines. The spectral index for $ -\mu \ll \varepsilon \ll T $ is $ \alpha +1 $. This thermal equilibrium scenario, combined with the leaky-box model and cosmic-ray observations, seems to suggest that the Galactic magnetic field is super-diffusive with $ \delta \approx 1.4 $.

UNIONS is an ongoing deep photometric multi-band survey of the Northern sky. As part of UNIONS, CFIS provides r-band data which we use to study weak-lensing peak counts for cosmological inference. We assess systematic effects for weak-lensing peak counts and their impact on cosmological parameters for the UNIONS survey. In particular, we present results on local calibration, metacalibration shear bias, baryonic feedback, the source galaxy redshift estimate, intrinsic alignment, and the cluster member dilution. For each uncertainty and systematic effect, we describe our mitigation scheme and the impact on cosmological parameter constraints. We obtain constraints on cosmological parameters from MCMC using CFIS data and MassiveNuS N-body simulations as a model for peak counts statistics. Depending on the calibration (local versus global, and the inclusion of the residual multiplicative shear bias), the mean matter density parameter $\Omega_m$ can shift up to $-0.024$ ($-0.5\sigma$). We also see that including baryonic corrections can shift $\Omega_m$ by $+0.027$ ($+0.5 \sigma$) with respect to the DM-only simulations. Reducing the impact of the intrinsic alignment and cluster member dilution through signal-to-noise cuts can lead to a shift in $\Omega_m$ of $+0.027$ ($+0.5 \sigma$). Finally, with a mean redshift uncertainty of $\Delta \bar{z} = 0.03$, we see that the shift of $\Omega_m$ ($+0.001$ which corresponds to $+0.02 \sigma$) is not significant. This paper investigates for the first time with UNIONS weak-lensing data and peak counts the impact of systematic effects. The value of $\Omega_m$ is the most impacted and can shift up to $\sim 0.03$ which corresponds to $0.5\sigma$ depending on the choices for each systematics. We expect constraints to become more reliable with future (larger) data catalogues, for which the current pipeline will provide a starting point.

We present the first unbiased survey of neutral hydrogen (HI) absorption in the Small Magellanic Cloud (SMC). The survey utilises pilot HI observations with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope as part of the Galactic ASKAP HI (GASKAP-HI) project whose dataset has been processed with the GASKAP-HI absorption pipeline, also described here. This dataset provides absorption spectra towards 229 continuum sources, a 275% increase in the number of continuum sources previously published in the SMC region, as well as an improvement in the quality of absorption spectra over previous surveys of the SMC. Our unbiased view, combined with the closely matched beam size between emission and absorption, reveals a lower cold gas faction (11%) than the 2019 ATCA survey of the SMC and is more representative of the SMC as a whole. We also find that the optical depth varies greatly between the SMC's bar and wing regions. In the bar we find that the optical depth is generally low (correction factor to the optically thin column density assumption of $\mathcal{R}_{\rm HI} \sim 1.04$) but increases linearly with column density. In the wing however, there is a wide scatter in optical depth despite a tighter range of column densities.

The accuracy of the estimated stellar atmospheric parameter decreases evidently with the decreasing of spectral signal-to-noise ratio (SNR) and there are a huge amount of this kind observations, especially in case of SNR$<$30. Therefore, it is helpful to improve the parameter estimation performance for these spectra and this work studied the ($T_\texttt{eff}, \log~g$, [Fe/H]) estimation problem for LAMOST DR8 low-resolution spectra with 20$\leq$SNR$<$30. We proposed a data-driven method based on machine learning techniques. Firstly, this scheme detected stellar atmospheric parameter-sensitive features from spectra by the Least Absolute Shrinkage and Selection Operator (LASSO), rejected ineffective data components and irrelevant data. Secondly, a Multi-layer Perceptron (MLP) method was used to estimate stellar atmospheric parameters from the LASSO features. Finally, the performance of the LASSO-MLP was evaluated by computing and analyzing the consistency between its estimation and the reference from the APOGEE (Apache Point Observatory Galactic Evolution Experiment) high-resolution spectra. Experiments show that the Mean Absolute Errors (MAE) of $T_\texttt{eff}, \log~g$, [Fe/H] are reduced from the LASP (137.6 K, 0.195 dex, 0.091 dex) to LASSO-MLP (84.32 K, 0.137 dex, 0.063 dex), which indicate evident improvements on stellar atmospheric parameter estimation. In addition, this work estimated the stellar atmospheric parameters for 1,162,760 low-resolution spectra with 20$\leq$SNR$<$30 from LAMOST DR8 using LASSO-MLP, and released the estimation catalog, learned model, experimental code, trained model, training data and test data for scientific exploration and algorithm study.

We update and extend our previous CMB anisotropy constraints on primordial magnetic fields through their dissipation by ambipolar diffusion and MHD decaying turbulence effects on the post-recombination ionization history. We derive the constraints using the latest Planck 2018 data release which improves on the E-mode polarization leading to overall tighter constraints with respect to Planck 2015. We also use the low-multipole E-mode polarization likelihood obtained by the SROLL2 map making algorithm and we note how it is compatible with larger magnetic field amplitudes than the Planck 2018 baseline, especially for positive spectral indices. The 95% CL constraints on the amplitude of the magnetic fields from the combination of the effects is $\sqrt{\langle B^2 \rangle} <0.69 (<0.72)$ nG for Planck 2018 (SROLL2) by marginalizing on the magnetic spectral index. We also investigate the impact of a damping scale allowed to vary and the interplay between the magnetic field effects and the lensing amplitude parameter.

This is a short summary of a project to construct a first principles cosmology of the Standard Model epoch, the period starting shortly before the electro-weak transition. The cosmology is derived from a simple initial state entirely within the SM and General Relativity. The initial state is semi-classical -- concentrated near a classical solution of the SM equations of motion -- and is precisely specified by a few simple conditions. The dark matter is a classical effect, a coherent state of the SU(2)-weak gauge field and the Higgs field. The leading order, classical universe contains only the dark matter. Ordinary matter is a correction due to the fluctuations of the SM fields around the classical trajectory. The initial state produces a homogeneous, isotropic, flat universe. There are no adjustable parameters. No physics beyond the SM is invoked. Only the classical calculations have been done so far. The time evolution of the fluctuations remains to be calculated.

We use cosmological hydrodynamic zoom-in simulations and semi-analytical models to study the effects of primordial black holes (PBHs) on first star formation. Our models self-consistently combine two competing effects: initial (isocurvature) perturbations induced by PBHs and BH accretion feedback. Focusing on PBHs with masses $\sim 30\ \rm M_{\odot}$, we find that the standard picture of first star formation in molecular-cooling minihaloes is not changed by PBHs, as the simulated star-forming gas clouds in the central parsec are very similar to those in the $\rm \Lambda CDM$ case when PBHs make up $f_{\rm PBH}\sim 10^{-4}-0.1$ of dark matter. With a dynamical friction timescale of $\sim 2-10\ \rm Myr$ when the central gas density reaches $10^{5}\ \rm cm^{-3}$, it is also unlikely that PBHs can sink into star-forming discs and affect the evolution of protostars, although they may interact with the stars during the main-sequence stage. At larger scales, PBHs tend to shift star formation to more massive haloes, and accelerate structure formation. The latter effect is stronger in regions with higher initial overdensities. For $f_{\rm PBH}\sim 10^{-4}-0.01$ (allowed by observational constraints), the collapsed mass fraction of haloes hosting Population III stars is similar (within a factor of $\sim2$ at $z\lesssim 30$) to that in $\rm \Lambda CDM$, implying that the impact of stellar-mass PBHs on the cosmic star formation history at $z\gtrsim 10$ is small. We also find that the Lyman-Werner photons from PBH accretion in atomic-cooling haloes may facilitate the formation of direct-collapse BHs.

UCHII regions are an important phase in the formation and early evolution of massive stars and a key component of the interstellar medium. The main objectives of this work are to study the young stellar objects (YSOs) associated with the G45.07+0.13 (G45.07) and G45.12+0.13 (G45.12) UCHII regions, as well as the interstellar medium in which they are embedded. We determined the distribution of the hydrogen column density (N(H2)) and dust temperature (Td) in the molecular cloud using Herschel images. We used infrared photometric data to identify and classify the YSOs. We also constructed a colour-magnitude diagram and K luminosity functions (KLFs) to compare the parameters of YSOs with the results of the radiative transfer models. We found that N(H2) varies from about 3.0 to 5.5x10^(23)cm^(-2) within the G45.07 and G45.12 regions, respectively. The maximum Td value is 35K in G45.12 and 42K in G45.07. Td reaches about 18-20K at distances of 2.6 and 3.7pc from IRAS19110+1045 (G45.07) and IRAS19111+1048 (G45.12), respectively. The gas-dust mass value included in G45.12 is 3.4x10^5Msun and 1.7x10^5Msun in G45.07. The UCHII regions are connected through a cold (Td=19K) bridge. The density distribution of 518 YSOs exhibits dense clusters around both IRAS sources. The parameters of YSOs in the IRAS clusters (124 objects) and 394 non-cluster objects surrounding them show some differences. About 75% of the YSOs belonging to the IRAS clusters have an evolutionary age greater than 10^6 years. Their slope alpha of the KLF agrees well with a Salpeter-type Initial Mass Function (IMF) (gamma=1.35) for a high mass range (O-F stars, beta=2) at 1Myr. The non-cluster objects are uniformly distributed in the molecular cloud, 80% of which are located to the right of the 0.1Myr isochrone. Based on the small age spread of the stellar objects, we suggest that the clusters originate from a single triggering shock.

We present LyMAS2, an improved version of the "Lyman-{\alpha} Mass Association Scheme" aiming at predicting the large-scale 3d clustering statistics of the Lyman-{\alpha} forest (Ly-{\alpha}) from moderate resolution simulations of the dark matter (DM) distribution, with prior calibrations from high resolution hydrodynamical simulations of smaller volumes. In this study, calibrations are derived from the Horizon-AGN suite simulations, (100 Mpc/h)^3 comoving volume, using Wiener filtering, combining information from dark matter density and velocity fields (i.e. velocity dispersion, vorticity, line of sight 1d-divergence and 3d-divergence). All new predictions have been done at z=2.5 in redshift-space, while considering the spectral resolution of the SDSS-III BOSS Survey and different dark matter smoothing (0.3, 0.5 and 1.0 Mpc/h comoving). We have tried different combinations of dark matter fields and found that LyMAS2, applied to the Horizon-noAGN dark matter fields, significantly improves the predictions of the Ly-{\alpha} 3d clustering statistics, especially when the DM overdensity is associated with the velocity dispersion or the vorticity fields. Compared to the hydrodynamical simulation trends, the 2-point correlation functions of pseudo-spectra generated with LyMAS2 can be recovered with relative differences of ~5% even for high angles, the flux 1d power spectrum (along the light of sight) with ~2% and the flux 1d probability distribution function exactly. Finally, we have produced several large mock BOSS spectra (1.0 and 1.5 Gpc/h) expected to lead to much more reliable and accurate theoretical predictions.

The dust in protoplanetary disks is an important ingredient in planet formation and can be investigated with model simulations and quantitative imaging polarimetry of the scattered stellar light. This study explores circumstellar disks with calculations for the intensity and polarization of the reflected light. The photon scattering and absorption by the disk are calculated with a Monte Carlo method for a grid of simple, rotationally symmetric models approximated at each point by a plane-parallel dusty atmosphere. The results of our simple disk models reproduce well the measurements for the intensity $I/I_{\star}$, azimuthal polarization $Q_{\varphi}/I_{\star}$, and fractional polarization $p_{\varphi}$ reported for well-observed transition disks. They describe the dependencies of the scattered radiation on the disk geometry and the dust scattering parameters in detail. Particularly strong constraints on disk properties can be obtained from certain diagnostic quantities. Therefore, they can be used as a diagnostic tool for the analysis of quantitative measurements, specifically in constraining or even determining the scattering properties of the dust particles in disks. Collecting and comparing such information for many systems is required to understand the nature of the scattering dust in planet-forming disks.

We investigated the presence of planetary companions around the nearby (7.6 pc) and bright ($V=9$ mag) early-type M dwarf Gl 514, analysing 540 radial velocities collected over nearly 25 years with the HIRES, HARPS, and CARMENES spectrographs. The data are affected by time-correlated signals at the level of 2-3 ms$^{-1}$ due to stellar activity, that we filtered out testing three different models based on Gaussian process regression. As a sanity cross-check, we repeated the analyses using HARPS radial velocities extracted with three different algorithms. We used HIRES radial velocities and Hipparcos-Gaia astrometry to put constraints on the presence of long-period companions, and we analysed TESS photometric data. We found strong evidence that Gl 514 hosts a super-Earth on a likely eccentric orbit, residing in the conservative habitable zone for nearly $34\%$ of its orbital period. The planet Gl 514 b has minimum mass $m_b\sin i_b=5.2\pm0.9$ $M_{\rm Earth}$, orbital period $P_b=140.43\pm0.41$ days, and eccentricity $e_b=0.45^{+0.15}_{-0.14}$. No evidence for transits is found in the TESS light curve. There is no evidence for a longer period companion in the radial velocities and, based on astrometry, we can rule out a $\sim0.2$ $M_{\rm Jup}$ planet at a distance of $\sim3-10$ au, and massive giant planets/brown dwarfs out to several tens of au. We discuss the possible presence of a second low-mass companion at a shorter distance from the host than Gl 514 b. Gl 514 b represents an interesting science case to study the habitability of planets on eccentric orbits. We advocate for additional spectroscopic follow-up to get more accurate and precise planetary parameters. Further follow-up is also needed to investigate sub \ms and shorter period signals.

The flattening rotation velocity $v(r)\to {\rm constant}$ found by Vera Rubin and collaborators and very apparent in the SPARC galaxy-rotation data coincides with Kepler's law in one less dimension. Thus, it is naturally reproduced by elongated dark matter distributions with the axis of prolateness perpendicular to the galactic plane. This theoretical understanding is borne out by the detailed fits to the rotation data that we here report: for equal dark matter profile, elongated distributions provide smaller $\chi^2$ than purely spherical ones. We also propose to use the geometric mean of the individual halo ellipticities, as opposed to their arithmetic average, because $s=c/a\in (0,\infty)$ corresponds to spherical haloes for $s=1$, so that the usually reported average is skewed towards oblateness and fails to reveal the large majority of prolate haloes. Several independently coded fitting exercises concur in yielding $s<1$ for most of the database entries and the oblate exceptions are understood and classified. This likely prolateness is of consequence for the estimated dark matter density near Earth.

The cosmic origin of fluorine is still under debate. Asymptotic giant branch (AGB) stars are among the few suggested candidates to efficiently synthesis F in our Galaxy, however their relative contribution is not clear. In this paper, we briefly review the theoretical studies from stellar yield models of the F synthesis and chemical equilibrium models of the F-containing molecules in the outflow around AGB stars. Previous detections of the F-bearing species towards AGB and post-AGB stars are also highlighted. We suggest that high-resolution ALMA observations of the AlF, one of the two main carriers of F in the outflow of AGB stars, can provide a reliable tracer of the F-budget in AGB stars. This will be helpful to quantify the role of AGB stars in the Galactic F budget.

Understanding how super-massive black holes form and grow in the early Universe has become a major challenge since the discovery of luminous quasars only 700 million years after the Big Bang. Simulations indicate an evolutionary sequence of dust-reddened quasars emerging from heavily dust-obscured starbursts that then transition to unobscured luminous quasars by expelling gas and dust. Although the last phase has been identified out to a redshift of 7.6, a transitioning quasar has not been found at similar redshifts owing to their faintness at optical and near-infrared wavelengths. Here we report observations of an ultraviolet compact object, GNz7q, associated with a dust-enshrouded starburst at a redshift of z=7.1899+/-0.0005. The host galaxy is more luminous in dust emission than any other known object at this epoch, forming 1,600 solar masses of stars per year within a central radius of 480 parsec. A red point source in the far-ultraviolet is identified in deep, high-resolution imaging and slitless spectroscopy. GNz7q is extremely faint in X-rays, which indicates the emergence of a uniquely ultraviolet compact star-forming region or a Compton-thick super-Eddington black-hole accretion disk at the dusty starburst core. In the latter case, the observed properties are consistent with predictions from cosmological simulations and suggest that GNz7q is an antecedent to unobscured luminous quasars at later epochs.

We consider a sample of Galactic classical Cepheids with highly accurate estimates of their distances taken from Skowron et al., where they were determined based on the period-luminosity relation. We have refined the geometric characteristics of two spiral arm segments-the Carina-Sagittarius and Outer ones. For this purpose, we have selected 269 Cepheids belonging to the Carina-Sagittarius arm with ages in the range 80-120 Myr. From them we have estimated the pitch angle of the spiral pattern $i=-11.9\pm0.2^\circ$ and the position of this arm $a_0=7.32\pm0.05$ kpc for the adopted $R_0=8.1\pm0.1$ kpc. In the Outer arm we have selected 343 Cepheids with ages in the range 120-300 Myr. From them we have found $i=-11.5\pm0.5^\circ$ and $a_0=12.89\pm0.06$ kpc. Adhering to the model of a grand-design spiral pattern in the Galaxy with one pitch angle for all arms, we can conclude that this angle is close to $-12^\circ$.

Her X-1/HZ Her is one of the best studied accreting X-ray pulsars. In addition to the pulsating and orbital periods, the X-ray and optical light curves of the source exhibit an almost periodic 35-day variability caused by a precessing accretion disk. The nature of the observed long-term stability of the 35-day cycle has been debatable. The X-ray pulse frequency of Her X-1 measured by the Fermi/GBM demonstrates periodical variations with X-ray flux at the Main-on state of the source. We explain the observed periodic sub-microsecond pulse frequency changes by the free precession of a triaxial neutron star with parameters previously inferred from an independent analysis of the X-ray pulse evolution over the 35-day cycle. In the Fermi/GBM data, we identified several time intervals with a duration of half a year or longer where the neutron star precession period describing the pulse frequency variations does not change. We found that the NS precession period varies within one per cent in different intervals. Such variations in the free precession period on a year time scale can be explained by <1% changes in the fractional difference between the triaxial neutron star's moments of inertia due to the accreted mass readjustment or variable internal coupling of the neutron star crust with the core.

We examine the light curve parameters of 97 nearby Type Ia supernovae in the ultraviolet and optical using observations from the Swift Ultra-Violet/Optical Telescope. Our light curve models used a linear combinations of templates, which were based on Functional Principal Component Analysis of the optical photometry of SNe Ia. The weights for each principal component used in the fit capture certain aspects of the light curve properties. We find that there is little difference in the overall variability of principal component scores across filters. We also find correlations between certain filters and PC components, indicating that the UV and optical properties seem to be tied together. This is a promising step in UV SNe Ia analysis, and suggests that UV light curves may be used to infer optical properties, paving the way for future cosmological studies.

Combining large segmented space telescopes, coronagraphy and wavefront control methods is a promising solution to produce a dark hole (DH) region in the coronagraphic image of an observed star and study planetary companions. The thermal and mechanical evolution of such a high-contrast facility leads to wavefront drifts that degrade the DH contrast during the observing time, thus limiting the ability to retrieve planetary signals. Lyot-style coronagraphs are starlight suppression systems that remove the central part of the image for an unresolved observed star, the point spread function, with an opaque focal plane mask (FPM). When implemented with a flat mirror containing an etched pinhole, the mask rejects part of the starlight through the pinhole which can be used to retrieve information about low-order aberrations. We propose an active control scheme using a Zernike wavefront sensor (ZWFS) to analyze the light rejected by the FPM, control low-order aberrations, and stabilize the DH contrast. The concept formalism is first presented before characterizing the sensor behavior in simulations and in laboratory. We then perform experimental tests to validate a wavefront control loop using a ZWFS on the HiCAT testbed. By controlling the first 11 Zernike modes, we show a decrease in wavefront error standard deviation by a factor of up to 9 between open- and closed-loop operations using the ZWFS. In the presence of wavefront perturbations, we show the ability of this control loop to stabilize a DH contrast around 7x10^-8 with a standard deviation of 7x10^-9. Active control with a ZWFS proves a promising solution in Lyot coronagraphs with an FPM-filtered beam to control and stabilize low-order wavefront aberrations and DH contrast for exoplanet imaging with future space missions.

Dust attenuation varies substantially from galaxy to galaxy and as of yet cannot be reproduced from first principles in theoretical models. In Nagaraj et al. (2022), we developed the first Bayesian population model of dust attenuation as a function of stellar population properties and projected galaxy shape, built on spectral energy distribution (SED) fits of nearly 30,000 galaxies in the 3D-HST grism survey with broadband photometric coverage from the rest-frame UV to IR. In this paper, we apply the model to galaxies from the large-volume cosmological simulation TNG100. We produce a UVJ diagram and compare it with one obtained in previous work by applying approximate radiative transfer to the simulated galaxies. We find that the UVJ diagram based on our empirical model is in better agreement with observations than the previous effort, especially in the number density of dusty star forming galaxies. We also construct the intrinsic dust-free UVJ diagram for TNG and 3D-HST galaxies at z ~ 1, finding qualitative agreement but residual differences at the 10-20% level. These differences can be largely attributed to the finding that TNG galaxies have, on average, 29% younger stellar populations and 0.28 dex higher metallicities than observed galaxies.

The possible time variation of the fundamental constants of nature has been an active subject of research since the large-number hypothesis was proposed by Dirac. In this paper, we propose a new method to investigate a possible time variation of the speed of light ($c$) along with the fine-structure constant ($\alpha$) using Strong Gravitational Lensing (SGL) and Type Ia Supernovae (SNe Ia) observations. We assume a general approach to describe the mass distribution of lens-type galaxies, the one in favor of the power-law index model (PLAW). We also consider the runaway dilaton model to describe a possible time-variation of $\alpha$. In order to explore the results deeply, we split the SGL sample into five sub-samples according to the lens stellar velocity dispersion and three sub-samples according to lens redshift. The results suggest that it is reasonable to treat the systems separately, but no strong indication of varying $c$ was found.

The initial mass-final mass relationship (IFMR) of white dwarfs (WD) represents a crucial benchmark for stellar evolution models, especially for the efficiency of mixing episodes and mass loss during the asymptotic giant branch (AGB) phase. In this study, we argue that such relation offers the opportunity of constraining the third dredge-up (3DU), with important consequences for chemical yields. The results are discussed in light of recent studies that have identified a kink in the IFMR for initial masses close to $2\,M_{\odot}$. Adopting a physically-sound approach in which the efficiency $\lambda$ of the 3DU varies as a function of core and envelope masses, we calibrate $\lambda$ in solar-metallicity TP-AGB models in order to reproduce the final masses of their WD progeny, over a the range of initial masses $0.9 \le M_{\rm i}/M_{\odot} \le 6$. In particular, we find that in low-mass stars with $1.4 \le M_{\rm i}/M_{\odot} \le 2.0$ the efficiency is small, $\lambda \le 0.3$, it steeply rises to about $\lambda \simeq 0.65$ in intermediate-mass stars with $2.0 \le M_{\rm i}/M_{\odot} \le 4.0$, and then it drops in massive TP-AGB stars with $4.0 \le M_{\rm i}/M_{\odot} \le 6.0$. Our study also suggests that a second kink may show up in the IFMR at the transition between the most massive carbon stars and those that are dominated by hot-bottom burning.

Galactic outflows driven by star formation and active galactic nuclei blow bubbles into their local environments, causing galactic magnetic fields to be carried into intergalactic space. We explore the redshift-dependent effect of these magnetized bubbles on the Faraday Rotation Measure (RM) of extragalactic radio sources. Using the IllustrisTNG cosmological simulations, we separate the contribution from magnetic bubbles from that of the volume-filling magnetic component expected to be due to the seed field originating in the Early Universe. We use this separation to extract the redshift dependence of each component and to compare TNG model predictions with observation measurements of the NRAO VLA Sky Survey (NVSS). We find that magnetized bubbles provide a sizeable contribution to the extragalactic RM, with redshift-independent $\langle |{\rm RM}| \rangle \simeq 13$ rad/m$^2$ for sources at redshifts $z\ge 2$. This is close to the mean residual RM of $16$ rad/m$^2$ found from NVSS data in this redshift range. Using the IllustrisTNG simulations, we also evaluate a simple model for the contribution to residual RM from individual host galaxies and show that this contribution is negligible at high-redshift. While the contribution from magnetic bubbles in the IllustrisTNG model is currently compatible with observational measurements of residual RM, the next-generation RM sky surveys, which will be free from the wrapping uncertainty, have larger statistics and better sensitivity should be able to observe predicted flat contribution from magnetic bubbles at large redshifts. This should allow to experimentally probe magnetic bubbles and check models of galaxy feedback in cosmological simulations.

It is critical for James Webb Space Telescope (JWST) science that instrumental units are converted to physical units. We detail the design of the JWST absolute flux calibration program that has the core goal of ensuring a robust flux calibration internal to and between all the science instruments for both point and extended source science. This program will observe a sample of calibration stars that have been extensively vetted based mainly on Hubble Space Telescope, Spitzer Space Telescope, and Transiting Exoplanet Survey Satellite observations. The program uses multiple stars of three different, well understood types (hot stars, A dwarfs, and solar analogs) to allow for the statistical (within a type) and systematic (between types) uncertainties to be quantified. The program explicitly includes observations to calibrate every instrument mode, further vet the set of calibration stars, measure the instrumental repeatability, measure the relative calibration between subarrays and full frame, and check the relative calibration between faint and bright stars. For photometry, we have set up our calibration to directly support both the convention based on the band average flux density and the convention based on the flux density at a fixed wavelength.

Observations of time-resolved thermal emission from tidally locked exoplanets can tell us about their atmospheric temperature structure. Telescopes such as JWST and ARIEL will improve the quality and availability of these measurements. This motivates an improved understanding of the processes that determine atmospheric temperature structure, particularly atmospheric circulation. The circulation is important in determining atmospheric temperatures, not only through its ability to transport heat, but also because any circulation pattern needs to be balanced by horizontal pressure contrasts, therefore implying a particular temperature structure. In this work, we show how the global temperature field on a tidally locked planet can be decomposed into contributions that are balanced by different components of the atmospheric circulation. These are the superrotating jet, stationary Rossby waves, and the divergent circulation. To achieve this, we partition the geopotential field into components balanced by the divergent circulation and the rotational circulation, with the latter comprising the jet and Rossby waves. The partitioned geopotential then implies a corresponding partitioning of the temperature via the hydrostatic relation. We apply these diagnostics to idealised general circulation model simulations, to show how the separate rotational and divergent circulations together make up the total three-dimensional atmospheric temperature structure. We also show how each component contributes distinct signatures to the thermal phase curve of a tidally locked planet. We conclude that this decomposition is a physically meaningful separation of the temperature field that explains its global structure, and can be used to fit observations of thermal emission.

Using a new method to estimate total galaxy mass (M$_{\rm T}$) and two samples of low luminosity galaxies containing measurements of the number of globular clusters (GCs) per galaxy (N$_{\rm GC}$), we revisit the N$_{\rm GC}-$M$_{\rm T}$ relation using a total of 203 galaxies, 157 of which have M$_{\rm T}$ $\ \le 10^{10}$ M$_\odot$. We find that the relation is nearly linear, N$_{\rm GC} \propto$ M$_{\rm T}^{0.92\pm0.08}$ down to at least M$_{\rm T} \sim 10^{8.75}$ M$_\odot$. Because the relationship extends to galaxies that average less than one GC per galaxy and to a mass range in which mergers are relatively rare, the relationship cannot be solely an emergent property of hierarchical galaxy formation. The character of the radial GC distribution in low mass galaxies, and the lack of mergers at these galaxy masses, also appears to challenge models in which the GCs form in central, dissipatively concentrated high-density, high-pressure regions and are then scattered to large radius. The slight difference between the fitted power-law exponent and a value of one, leaves room for a shallow M$_{\rm T}$-dependent variation in the mean mass per GC that would allow the relation between total mass in GCs and M$_{\rm T}$ to be linear.

In recent work, we examined how different modes in the ringdown phase of a binary coalescence are excited as a function of the final plunge geometry. At least in the large mass ratio limit, we found a clean mapping between angles describing the plunge and the amplitude of different quasi-normal modes (QNMs) which constitute the ringdown. In this study, we use that mapping to construct a waveform model expressed as a sum of QNMs where the mode amplitudes and phases are determined by the source plunge parameters. We first generate a large number of calibration waveforms and interpolate between fits of each mode amplitude and phase up to $\ell \leq 8$ and $\ell - |m| \leq 4$. The density of our calibration data allows us to resolve important features such as phase transition discontinuities at large misalignments. Using our ringdown waveform model, we then perform Bayesian parameter estimation with added white Gaussian noise to demonstrate that, in principle, the mode amplitudes can be measured and used to constrain the plunge geometry. We find that inferences are substantially improved by incorporating prior information constraining mode excitation, which motivates work to understand and characterize how the QNM excitation depends on the coalescence geometry. These results are part of a broader effort to map the mode excitation from arbitrary masses and spins, which will be useful for characterizing ringdown waves in upcoming gravitational-wave measurements.

We consider a decaying scalar dark matter (DM) with mass $m_\chi$ in the range 10 GeV - 10 TeV and vary the branching ratios of all possible two-body SM final states (excluding and including $\nu\bar{\nu}$) in the range $0\%-100\%$ to derive constraints on the total decay width $\Gamma$ by combining the data collected by several astrophysical and cosmological observations. We find that, $\Gamma \lesssim 10^{-26} - 10^{-27}\,{\rm s}^{-1}$ (excluding $\nu\bar{\nu}$) and $\Gamma \lesssim 10^{-24} - 10^{-26}\,{\rm s}^{-1}$ (including $\nu\bar{\nu}$) are allowed, depending on the values of $m_\chi$, which are most robust upper limits on $\Gamma$ for a generic decaying scalar DM. We then investigate the prospect of the upcoming Square Kilometer Array (SKA) radio telescope in detecting the DM decay induced radio signals originating inside the dwarf spheroidal (dSph) galaxies. We have classified the DM parameter space, allowed by the existing observations, independently of the branching ratio in each individual two-body SM final state, based on the detectability at the SKA. Excluding the $\nu\bar{\nu}$ decay mode, we find that, throughout the DM mass range 10 GeV - 10 TeV, $\Gamma \gtrsim 10^{-30}\,{\rm s}^{-1} - 10^{-29}\,{\rm s}^{-1}$ is detectable for all possible branching ratio combinations at the SKA (for 100 hours of observation time), with conservative choices for the relevant astrophysical parameters. On the other hand, when arbitrary branching ratios are allowed also for the $\nu\bar{\nu}$ decay mode, DM decays can be probed model-independently for $\Gamma \gtrsim 2 \times 10^{-29}\,{\rm s}^{-1}$ provided DM masses are greater than a few hundreds of GeV.

The development of high-resolution, large-baseline optical interferometers would revolutionize astronomical imaging. However, classical techniques are hindered by physical limitations including loss, noise, and the fact that the received light is generally quantum in nature. We show how to overcome these issues using quantum communication techniques. We present a general framework for using quantum error correction codes for protecting and imaging starlight received at distant telescope sites. In our scheme, the quantum state of light is coherently captured into a non-radiative atomic state via Stimulated Raman Adiabatic Passage, which is then imprinted into a quantum error correction code. The code protects the signal during subsequent potentially noisy operations necessary to extract the image parameters. We show that even a small quantum error correction code can offer significant protection against noise. For large codes, we find noise thresholds below which the information can be preserved. Our scheme represents an application for near-term quantum devices that can increase imaging resolution beyond what is feasible using classical techniques.

In this work we shall calculate in detail the effect of an GW170817-compatible Einstein-Gauss-Bonnet theory on the energy spectrum of the primordial gravitational waves. The spectrum is affected by two characteristics, the overall amplification/damping factor caused by the GW170817-compatible Einstein-Gauss-Bonnet theory and by the tensor spectral index and the tensor-to-scalar ratio. We shall present the formalism for studying the inflationary dynamics and post-inflationary dynamics of GW170817-compatible Einstein-Gauss-Bonnet theories for all redshifts starting from the radiation era up to the dark energy era. We exemplify our formalism by using two characteristic models, which produce viable inflationary and dark energy eras. As we demonstrate, remarkably the overall damping/amplification factor is of the order of unity, thus the GW170817-compatible Einstein-Gauss-Bonnet models affect the primordial gravitational waves energy spectrum only via their tensor spectral index and the tensor-to-scalar ratio. Both models have a blue tilted tensor spectrum, and therefore the predicted energy spectrum of the primordial gravity waves can be detectable by most of the future gravitational waves experiments, for various reheating temperatures.

In this study, we present the phase-space analysis of Quintessence models specified by the choice of two potentials, namely the Recliner potential and what we call the broken exponential-law potential, which is a new proposal. Using a dynamical system analysis we provide a systematic study of the cosmological evolution of the two models and their properties. We find new scaling solutions characterised by a constant ratio between the energy density of the scalar field and that of the matter component. These solutions are of high interest in light of the possibility to alleviate the coincidence problem. Additionally, the models also show attractor solutions. We finally construct concrete models built using a double potential according to which one potential realises the early-time scaling regime and the second one allows to exit this regime and to enter in the epoch of cosmic acceleration driven by a scalar-field dominated attractor point.

Following our earlier work on the 3-dimensional effective velocity distribution of Galactic WIMPs (not only impinging on our detectors but also) scattering off target nuclei, in this paper, we demonstrate the normal and a "reverse" annual modulations of elastic WIMP-nucleus scattering signals, which could be observed in direct Dark Matter detection experiments. Our simulations show that, once the WIMP mass is as light as only a few tens GeV, the event number and the accumulated recoil energy of WIMP-induced scattering events off both of light and heavy target nuclei would indeed be maximal (minimal) in summer (winter). However, once the WIMP mass is as heavy as a few hundreds GeV, the event number and the accumulated recoil energy of WIMP scattering events off heavy nuclei would inversely be minimal in summer. Understandably, for an intermediate WIMP mass, the event number and the accumulated recoil energy of scattering events off some middle-mass nuclei would show an approximately uniform time dependence.

In this paper we address the hierarchy in the neutron-anti-neutron ($n-\bar{n}$) oscillation through clockwork mechanism (CW) with a set of mirror neutrons acting as clockwork gears. This is achieved by coupling Standard Model (SM) neutrons to the zeroth gear while having an explicit mirror baryon number ($\bar{B}$) breaking term at the $N^{th}$ one. The explicit $\bar{B}$ breaking term induces neutron oscillation in the mirror world which propagates to the SM baryon number ($B$) violation through the CW gears. This mirror world is a Standard Model singlet and possess the $\bar{B}$ symmetry that prohibits its decay to SM matter fields. This symmetry is broken to the combination, $B-\bar{B}$, by a scalar field to generate the clockwork. Without introducing large hierarchies among the mirror neutrons, we predict the time period of neutron-anti-neutron oscillation in our world of $\tau_{n\bar{n}} \gtrsim 10^8 s$, which is within the reach of current experiments and the upcoming ones. This mechanism also predicts neutron to mirror neutron transition along with the emission of SM particles, given the phase space is favourable. This could be searched for at neutron stars and will place bounds on the symmetry breaking mechanism in the mirror world. The current constraint on neutron-mirror neutron oscillation, coming from neutron stars, is shown to be satisfied, in the parameter space which predicts a detectable neutron-anti-neutron oscillation in our world. This scenario also predicts double neutron annihilation to produce two mirror neutrons which could also be tested via the luminosity measurements of neutron stars.

The origin of life might be sparked by the polymerization of the first RNA molecules in Darwinian ponds during wet-dry cycles. The key life-building block ribose was found in carbonaceous chondrites. Its exogenous delivery onto the Hadean Earth could be a crucial step toward the emergence of the RNA world. Here, we investigate the formation of ribose through a simplified version of the formose reaction inside carbonaceous chondrite parent bodies. Following up on our previous studies regarding nucleobases with the same coupled physico-chemical model, we calculate the abundance of ribose within planetesimals of different sizes and heating histories. We perform laboratory experiments using catalysts present in carbonaceous chondrites to infer the yield of ribose among all pentoses (5Cs) forming during the formose reaction. These laboratory yields are used to tune our theoretical model that can only predict the total abundance of 5Cs. We found that the calculated abundances of ribose were similar to the ones measured in carbonaceous chondrites. We discuss the possibilities of chemical decomposition and preservation of ribose and derived constraints on time and location in planetesimals. In conclusion, the aqueous formose reaction might produce most of the ribose in carbonaceous chondrites. Together with our previous studies on nucleobases, we found that life-building blocks of the RNA world could be synthesized inside parent bodies and later delivered onto the early Earth.

General Relativity (GR) was proven via the direct detection of gravitational waves from the mergers of the binary black holes and binary neutron stars by the Advanced LIGO and Advanced Virgo detectors. These detections confirmed the prediction of GR and provided the first direct evidence of the existence of stellar-mass black holes (BHs). However, the occurrence of singularities at the centers of BHs suggests that GR is inapplicable because of the breakdown of the equivalence principle at the singularities. The fact that these singularities exist indicates that GR cannot be a universal theory of space-time. In the low-energy limit, the theoretical and observational challenges faced by the $\Lambda$CDM model also indicate that we might have to look beyond GR as the underlying theory of gravity. Unlike GR, whose field equations contain only up to second-order derivatives, the modified theories with higher derivative Ricci/Riemann tensor gravity models include higher derivatives. Therefore, one expects significant differences between GR and modified theories. Since there are many ways of modifying GR in the strong-gravity and cosmological distances, each model has unique features. This leads to the following crucial question: Are there a set of unique signatures that distinguish GR from modified gravity (MG) theories? This review discusses three aspects of MG theories: (1) Why do we need to consider MG theories? (2) How to modify GR? and (3) What are the observational consequences? The review is written in a pedagogical style with the expectation that it will serve as a useful reference for theorists and observers and those interested in bridging the divide between theory and observations.