Papers for Thursday, Aug 24 2023

We calculated empirical reddening and extinction coefficients with respect to the dust reddening map of Schlegel et al. for the Swift-UVOT passbands, using the 'star pair' method and photometric data from the UVOT Serendipitous Source Catalogue and spectroscopic data from LAMOST Data Release 7 and 2MASS. The reddening coefficients for the UVW2-UVM2, UVM2-UVW1, UVW1-U, U-B, and B-V colors are -1.39, 2.08, 0.78, 0.72, and 0.84, respectively.The extinction coefficients for the UVW2, UVM2, UVW1, U, B, and V bands are 5.60, 6.99, 4.91, 4.13, 3.41 and 2.57, respectively. The numbers are consistent with predictions by the Fitzpatrick's extinction law of R(V)=3.0. Temperature-dependent variations of the coefficients are found and discussed, particularly for the ultraviolet passbands. We recommend using the new reddening and extinction coefficients in future when dereddening the Swift-UVOT data.

Neutron star surfaces and atmospheres are unique environments that sustain the largest-known magnetic fields in the universe. Our knowledge of neutron star material properties, including the composition and equation of state, remains highly unconstrained. Electron-atom collisions are integral to theoretical thermal conduction and spectral emission models that describe neutron star surfaces. The theory of scattering in magnetic fields was developed in the 1970s, but focused only on bare nuclei scattering. In this work, we present a quantum treatment of atom-electron collisions in magnetic fields; of significant importance is the inclusion of Pauli repulsion arising from two interacting electrons. We find strange behaviors not seen in collisions without a magnetic field. In high magnetic fields, Pauli repulsion can lead to orders of magnitude enhancements of collision cross sections. Additionally, the elastic collision cross sections that involve the ground state become comparable to those involving excited states, and states with large orbits have the largest contribution to the collisions. We anticipate significant changes to transport properties and spectral line broadening in neutron star surfaces and atmospheres, which will aid in spectral diagnostics of these extreme environments.

Theoretical models and observations suggest that the abundances of molecular ions in protoplanetary disks should be highly sensitive to the variable ionization conditions set by the young central star. We present a search for temporal flux variability of HCO+ J=1-0, which was observed as a part of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We split out and imaged the line and continuum data for each individual day the five sources were observed (HD 163296, AS 209, GM Aur, MWC 480, and IM Lup, with between 3 to 6 unique visits per source). Significant enhancement (>3\sigma) was not observed, but we find variations in the spectral profiles in all five disks. Variations in AS 209, GM Aur, and HD 163296 are tentatively attributed to variations in HCO+ flux, while variations in IM Lup and MWC 480 are most likely introduced by differences in the \textit{uv} coverage, which impact the amount of recovered flux during imaging. The tentative detections and low degree of variability are consistent with expectations of X-ray flare driven HCO+ variability, which requires relatively large flares to enhance the HCO+ rotational emission at significant (>20%) levels. These findings also demonstrate the need for dedicated monitoring campaigns with high signal to noise ratios to fully characterize X-ray flare driven chemistry.

Ly$\alpha$ damping wings in the spectra of bright objects at high redshift are a useful probe of the ionization state of the intergalactic medium during the reionization epoch. It has recently been noted that, despite the inhomogeneous nature of reionization, these damping wings have a characteristic shape which is a strong function of the volume-weighted average neutral hydrogen fraction of the intergalactic medium. We present here a closer examination of this finding using a simulation of patchy reionization from the Sherwood-Relics simulation suite. We show that the characteristic shape and scatter of the damping wings are determined by the average neutral hydrogen density along the line of sight, weighted by its contribution to the optical depth producing the damping wing. We find that there is a redshift dependence in the characteristic shape due to the expansion of the Universe. Finally, we show that it is possible to differentiate between the shapes of damping wings in galaxies and young (or faint) quasars at different points in the reionization history at large velocity offsets from the point where the transmission first reaches zero.

We report on 85-101 GHz light curves of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), observed in April 2017 with the Atacama Large Millimeter/submillimeter Array (ALMA). This study of high-cadence full-Stokes data provides new measurements of the fractional linear polarization at a 1-2% level resolved in 4 s time segments, and stringent upper limits on the fractional circular polarization at 0.3%. We compare these findings to ALMA light curves of Sgr A* at 212-230 GHz observed three days later, characterizing a steep depolarization of the source at frequencies below about 150 GHz. We obtain time-dependent rotation measure (RM) measurements, with the mean RM at 85-101 GHz being a factor of two lower than that at 212-230 GHz. Together with the rapid temporal variability of the RM and its different statistical characteristics in both frequency bands, these results indicate that the Faraday screen in Sgr A* is largely of internal character, with about half of the Faraday rotation taking place inside the inner 10 gravitational radii, contrary to the common external Faraday screen assumption. We then demonstrate how this observation can be reconciled with theoretical models of radiatively inefficient accretion flows for a reasonable set of physical parameters. Comparisons with numerical general relativistic magnetohydrodynamic simulations suggest that the innermost part of the accretion flow in Sgr A* is much less variable than what these models predict, in particular, the observed magnetic field structure appears to be coherent and persistent.

The detection of emission lines associated with accretion processes is a direct method for studying how and where gas giant planets form, how young planets interact with their natal protoplanetary disk and how volatile delivery to their atmosphere takes place. H$\alpha$ ($\lambda=0.656\,\mu$m) is expected to be the strongest accretion line observable from the ground with adaptive optics systems, and is therefore the target of specific high-contrast imaging campaigns. We present MagAO-X and HST data obtained to search for H$\alpha$ emission from the previously detected protoplanet candidate orbiting AS209, identified through ALMA observations. No signal was detected at the location of the candidate, and we provide limits on its accretion. Our data would have detected an H$\alpha$ emission with $F_\mathrm{H\alpha}>2.5\pm0.3 \times10^{-16}$ erg s$^{-1}$ cm$^{-2}$, a factor 6.5 lower than the HST flux measured for PDS70b (Zhou et al., 2021). The flux limit indicates that if the protoplanet is currently accreting it is likely that local extinction from circumstellar and circumplanetary material strongly attenuates its emission at optical wavelengths. In addition, the data reveal the first image of the jet north of the star as expected from previous detections of forbidden lines. Finally, this work demonstrates that current ground-based observations with extreme adaptive optics systems can be more sensitive than space-based observations, paving the way to the hunt for small planets in reflected light with extremely large telescopes.

The "tuning-fork" (TF) analysis of CO and Halpha emission has been used to estimate the lifetimes of molecular clouds in nearby galaxies. With simple model calculations, we show that this analysis does not necessarily estimate cloud lifetimes, but instead captures a duration of the cloud evolutionary cycle, from dormant to star forming, and then back to a dormant phase. We adopt a hypothetical setup in which molecular clouds (e.g., traced in CO) live forever and form stars (e.g., HII regions) at some frequency, which then drift away from the clouds. The TF analysis still returns a timescale for the immortal clouds. This model requires drifting motion to separate the newborn stars from the clouds, and we discuss its origin. We also discuss the physical origin of the characteristic spatial separation term in the TF analysis and a bias due to systematic error in the determination of the reference timescale.

We model the interstellar dust content of the reionization era with a suite of cosmological, fluid-dynamical simulations of galaxies with stellar masses ranging from $\sim 10^5 - 10^9 M_{\odot}$ in the first $1.2$ billion years of the universe. We use a post-processing method that accounts for dust creation and destruction processes, allowing us to systematically vary the parameters of these processes to test whether dust-dependent observable quantities of galaxies at these epochs could be useful for placing constraints on dust physics. We then forward model observable properties of these galaxies to compare to existing data. We find that we are unable to simultaneously match existing observational constraints with any one set of model parameters. Specifically, the models which predict the largest dust masses $D/Z \gtrsim 0.1$ at $z = 5$ -- because of high assumed production yields and/or efficient growth via accretion in the interstellar medium -- are preferred by constraints on total dust mass and infrared luminosities, but these models produce far too much extinction in the ultraviolet, preventing them from matching observations of $\beta_{\rm UV}$. To investigate this discrepancy, we analyze the relative spatial distribution of stars and dust as probed by infrared (IR) and ultraviolet (UV) emission, which appear to exhibit overly symmetric morphologies compared to existing data, likely due to the limitations of the stellar feedback model used in the simulations. Our results indicate that the observable properties of the dust distribution in high redshift galaxies are a particularly strong test of stellar feedback.

We conduct a comprehensive spectro-temporal analysis of repeating Fast Radio Bursts (FRBs) utilizing nine distinct sources, the largest sample to date. Our data set includes 175 sub-bursts and 31 multi-component bursts from 11 data sets, with centre frequencies ranging from 149-7144 MHz and durations spanning from 73 \mu s-13 ms. Our findings are consistent with the predictions of the Triggered Relativistic Dynamical Model (TRDM) of FRB emission. We affirm the predicted quadratic relationship between sub-burst slope and frequency, as well as a linear dependence of bandwidth on frequency that is consistent with mildly-relativistic Doppler broadening of narrow-band emission. Most importantly, we confirm the sub-burst slope law, a predicted inverse relationship between sub-burst slope and duration, to hold consistently across different sources. Remarkably, we also discover that the drift rates of multi-component bursts follow the same law as the sub-burst slopes, an unexplained result that warrants further investigation. These findings not only support the TRDM as a viable framework for explaining several aspects of FRB emission, but also provide new insights into the complex spectro-temporal properties of FRBs.

We report on the close similarity of coronal mass ejection (CME) properties in ground level enhancement (GLE) in solar energetic particle (SEP) events and sustained gamma ray emission (SGRE) from the Sun as indicated by low frequency type III radio bursts observed in the interplanetary medium. The complex type III bursts have an average 1 MHz duration of 36 and 34 min in the SGRE and GLE events, respectively. Similarly, the CMEs underlying SGRE and GLE have average space speeds of 1866 and 2084 km/s, respectively. These are larger than the corresponding values (32 min, 1407 km/s) for a control sample of type III bursts associated with frontside halo CMEs with sky plane speed exceeding 800 km/s. These results are consistent with the idea that energetic CME driven shocks accelerate particles to very high energies that are responsible for GLE and SGRE events.

The Origins, Spectral Interpretation, Resource Identification, and Security Regolith Explorer (OSIRIS-REx) spacecraft arrived at its target, near-Earth asteroid 101955 Bennu, in December 2018. After one year of operating in proximity, the team selected a primary site for sample collection. In October 2020, the spacecraft descended to the surface of Bennu and collected a sample. The spacecraft departed Bennu in May 2021 and will return the sample to Earth in September 2023. The analysis of the returned sample will produce key data to determine the history of this B-type asteroid and that of its components and precursor objects. The main goal of the OSIRIS-REx Sample Analysis Plan is to provide a framework for the Sample Analysis Team to meet the Level 1 mission requirement to analyze the returned sample to determine presolar history, formation age, nebular and parent-body alteration history, relation to known meteorites, organic history, space weathering, resurfacing history, and energy balance in the regolith of Bennu. To achieve this goal, this plan establishes a hypothesis-driven framework for coordinated sample analyses, defines the analytical instrumentation and techniques to be applied to the returned sample, provides guidance on the analysis strategy for baseline, overguide, and threshold amounts of returned sample, including a rare or unique lithology, describes the data storage, management, retrieval, and archiving system, establishes a protocol for the implementation of a micro-geographical information system to facilitate co-registration and coordinated analysis of sample science data, outlines the plans for Sample Analysis Readiness Testing, and provides guidance for the transfer of samples from curation to the Sample Analysis Team.

WR 137 (HD 192641) is a binary system consisting of a carbon-rich Wolf-Rayet star and an Oe companion star in a 13-year orbit. Near periastron, the winds of the two stars collide and form carbonaceous dust. We obtained three mid-infrared grism spectra of the system with SOFIA and FORCAST during the last year of SOFIA's operations in July 2021, February 2021, and May 2022 (Cycle 9). Within these spectra, we have identified several wind lines from He I, He II, C III, and C IV that are emitted from the Wolf-Rayet wind as well as a weak emission feature around 6.3-6.4 $\mu$m that may have shifted its peak flux from 6.29 to 6.41$\mu$m through this time period. The weak feature grew as the continuum dust emission grew while the WR emission appeared to decline due to lower contrast with the continuum. Furthermore, we observe that the peak of the feature shifts to redder wavelengths during the observations. We compare this feature to the UIR feature and other emission lines identified in dusty WC binaries. For WR 137, we speculate that mixing of the winds in the system with the Oe star's disk is important for starting the dust formation and that it is less important as dust formation continues. Previous infrared photometry shows "mini-eruptions" of dust production which could then be explained with variations of the Oe star disk.

Context. The detection of Earth-like planets with the radial-velocity (RV) method is extremely challenging today due to the presence of non-Doppler signatures such as stellar activity and instrumental signals that mimic and hide the signals of exoplanets. In a previous paper, we presented the YARARA pipeline, which implements corrections for telluric absorption, stellar activity and instrumental systematics at the spectral level, then extracts line-by-line (LBL) RVs with significantly better precision than standard pipelines. Aims. In this paper, we demonstrate that further gains in RVs precision can be achieved by performing Principal Component Analysis (PCA) decomposition on the LBL RVs. Methods. The mean-insensitive nature of PCA means that it is unaffected by true Doppler shifts, and thus can be used to isolate and correct nuisance signals other than planets. Results. We analysed the data of 20 intensively observed HARPS targets by applying our PCA approach on the LBL RVs obtained by YARARA. The first principal components show similarities across most of the stars and correspond to newly identified instrumental systematics, which we can now correct for. For several targets, this results in an unprecedented RV root-mean-square of around 90 cm/s over the full lifetime of HARPS. We use the corrected RVs to confirm a previously published 120-day signal around 61Vir, and to detect a Super-Earth candidate (K = 60 +/- 6 cm/s, m sin i = 6.6 +/- 0.7 Earth mass) around the G6V star HD20794, which spends part of its 600-day orbit within the habitable zone of the host star. Conclusions. This study highlights the potential of LBL PCA to identify and correct hitherto unknown, long-term instrumental effects and thereby extend the sensitivity of existing and future instruments towards the Earth analogue regime.

We examined the grain size in the dust ring encircling the 0.19~$M_\sun$ T Tauri star CIDA 1 using the Karl G. Jansky Very Large Array (JVLA) at multiple centimeter wavelengths, with a spatial resolution of 0$\farcs$2--0$\farcs$9. We detected distinct partial-ring structures at these wavelengths around CIDA~1. Based on spatial distributions and spectral indexes, we determined that these centimeter emissions originated from dust, rather than free-free or synchrotron emissions. To estimate the maximum grain size ($a_{\rm max}$) within the ring, we compared the observed spectral energy distribution (SED) with SEDs calculated for different $a_{\rm max}$ values using radiative transfer calculations. Our findings indicate an $a_{\rm max}$ value of approximately 2.5~cm in the ring, assuming the dust opacity can be approximated by the DSHARP models. These results suggest that grain growth took place within the CIDA~1 ring, potentially facilitating more efficient planet formation through pebble accretion scenarios involving centimeter-sized pebbles.

In the process of producing the roughly three ionizing photons per atom required to reionize the IGM, the same massive stars explode and eject metals into their surroundings. While the overly sensitive Lya transition makes Gunn-Peterson absorption of background quasar light an ineffective probe of reionization at z > 6, strong low-ionization transitions like the MgII doublet will give rise to a detectable "metal-line forest", if metals pollute the neutral IGM. We measure the auto-correlation of the MgII forest transmission using a sample of ten ground based z >= 6.80 quasar spectra probing the redshift range 5.96 < z_MgII < 7.42 (z_MgII,median = 6.47). The correlation function exhibits strong small-scale clustering and a pronounced peak at the doublet velocity (768 km/s) arising from strong absorbers in the CGM of galaxies. After these strong absorbers are identified and masked the signal is consistent with noise. Our measurements are compared to a suite of models generated by combining a large hydrodynamical simulation with a semi-numerical reionization topology, assuming a simple uniform enrichment model. We obtain a 95% credibility upper limit of [Mg/H] < -3.73 at z_MgII,median = 6.47, assuming uninformative priors on [Mg/H] and the IGM neutral fraction x_HI. Splitting the data into low-z (5.96 < z_MgII < 6.47; z_MgII,median = 6.235) and high-z (6.47 < z_MgII < 7.42; z_MgII,median = 6.72) subsamples again yields null-detections and 95% upper limits of [Mg/H] < -3.75 and [Mg/H] < -3.45, respectively. These first measurements set the stage for an approved JWST Cycle 2 program (GO 3526) targeting a similar number of quasars that will be an order of magnitude more sensitive, making the Mgii forest an emerging powerful tool to deliver precision constraints on the reionization and enrichment history of the Universe.

EX Lup is the archetype for the class of young stars that undergoes repeated accretion outbursts of $\sim 5$ mag at optical wavelengths and that last for months. Despite extensive monitoring that dates back 130 years, the accretion history of EX Lup remains mostly qualitative and has large uncertainties. We assess historical accretion rates of EX Lup by applying correlations between optical brightness and accretion, developed on multi-band magnitude photometry of the $\sim 2$ mag optical burst in 2022. Two distinct classes of bursts occur: major outbursts ($\Delta V\sim5$ mag) have year-long durations, are rare, reach accretion rates of $\dot{M}_{\rm acc}\sim10^{-7}~M_\odot~{\rm yr^{-1}}$ at peak, and have a total accreted mass of around 0.1 Earth masses. The characteristic bursts ($\Delta V\sim2$ mag) have durations of $\sim 2-3$ months, are more common, reach accretion rates of $\dot{M}_{\rm acc}\sim10^{-8}~M_\odot~{\rm yr^{-1}}$ at peak, and have a total accreted mass of around $10^{-3}$ Earth masses. The distribution of total accreted mass in the full set of bursts is poorly described by a power law, which suggests different driving causes behind the major outburst and characteristic bursts. The total mass accreted during two classes of bursts is around two times the masses accreted during quiescence. Our analysis of the light curves reveals a color-dependent time lag in the 2022 post-burst light curve, attributed to the presence of both hot and cool spots on the stellar surface.

The study reports photometric and spectroscopic observations of two recently recognised contact binary systems. Both systems show total eclipses and analysis of the light curves indicate both have a very low mass ratios of less than 0.3. We derive absolute parameters from colour and distance based calibrations and show that although both have low mass ratios they are likely to be in a stable orbit and unlikely to merge. In other respects both systems have characteristics similar to other contact binaries with the secondary larger and brighter than main sequence counterparts and we also find that the secondary is considerably denser than the primary in both systems.

The formation of black hole-neutron star (BH-NS) or BH-BH systems may be accompanied with special supernova (SN) signals, due to the accretion feedback from the companion BH. The additional heating, which is mainly attributed to the Blandford-Payne mechanism, would disrupt the isotropic nature of the luminosity distribution on the surface of the SN ejecta, leading to the appearance of polarization. Here we develop a three dimensional (3D) Monte Carlo polarization simulation code (MCPSC) to conduct simulations for these special SNe. We find that the maximum polarization level of approximately \sim 2 occurs at the peak time of SN emission in the "close-binary" scenario, while in the "faraway-binary" case, maximum polarization (i.e. \sim 0.7) is observed at a considerably later time than the peak of the SN. The magnitude of polarization is dependent on the degree of unevenness in the luminosity distribution and the angle between the line of sight and the equatorial direction. When considering the geometric distortion of supernova ejecta at the same time, the magnitude of polarization may either increase (for a oblate ellipsoidal shape) or decrease (for a prolate ellipsoidal shape). The polarization signatures represent a promising auxiliary instrument to facilitate the identification of the companion-fed SNe. Moreover, by comparing the event rate of these special SNe with the event rate density of LIGO-Virgo detected BH-NS/BH systems could further help to distinguish the BH-NS/BH formation channel.

The kinematics of coronal mass ejections (CMEs) are essential for understanding their initiation mechanisms and predicting their planetary impact. Most acceleration and deceleration occur below 4 R$\odot$, which is crucial for initiation understanding. Furthermore, the kinematics of CMEs in the inner corona ($<$ 3 R$_\odot$) are closely related to their propagation in the outer corona and their eventual impact on Earth. Since the CME kinematics are mainly probed using coronagraph data, it is crucial to investigate how imaging cadence affects the precision of data analysis and conclusions drawn and also for determining the flexibility of designing observational campaigns with upcoming coronagraphs. We study ten CMEs observed by the K-Coronagraph of the MLSO. We manually track the CMEs using high cadence (15 s) white-light observations of K-Cor and vary the cadence as 30 s, 1 min, 2 min, and 5 min to study the impact of cadence on the kinematics. We also employed the bootstrapping method to estimate the fitting parameters. Our results indicate that the average velocity of the CMEs does not have a high dependence on the imaging cadence, while the average acceleration shows significant dependence on the same, with the confidence interval showing significant shifts for the average acceleration for different cadences. The decrease in cadence also influences the determination of acceleration onset time. We further find that it is difficult to find an optimum cadence to study all CMEs, as it is also influenced by the pixel resolution of the instrument and the speed of the CME. However, except for very slow CMEs (speeds less than 300 Kms$^{-1}$), our results indicate a cadence of 1 min to be reasonable for the study of their kinematics. The results of this work will be important in the planning of observational campaigns for the existing and upcoming missions that will observe the inner corona.

Ultra-high energy cosmic rays (UHECRs), whose energy are beyond $10^{18}~\mathrm{eV}$, are the most energetic particles we have ever detected. The latest results seem to indicate a heavier composition at the highest energies, complicating the search for their origins. Due to the limited number of UHECR events, we need to build an instrument with an order of magnitude larger effective-exposure to collect UHECRs in future decades. The Fluorescence detector Array of Single-pixel Telescopes (FAST) is a proposed low-cost, easily deployable UHECR detector suitable for a future ground array. It is essential to validate the telescope design and autonomous observational techniques using prototypes located in both hemispheres. Here we report on the current status of observations, recent performance results of prototypes, and developments towards a future mini-array.

The Fluorescence detector Array of Single-pixel Telescopes (FAST) is one of several proposed designs for a next-generation cosmic-ray detector. Such detectors will require enormous collecting areas whilst also needing to remain cost-efficient. To meet these demands, the FAST collaboration has designed a simplified, low-cost fluorescence telescope consisting of only four photomultiplier tubes (PMTs). Since standard air shower reconstruction techniques cannot be used with so few PMTs, FAST utilises an alternative two-step approach. In the first step, a neural network is used to provide a first estimate of the true shower parameters. This estimate is then used as the initial guess in a minimisation procedure where the measured PMT traces are compared to simulated ones, and the best-fit shower parameters are found. A detailed explanation of these steps is given, with the expected performance of FAST prototypes at the Telescope Array experiment acting as a demonstration of the technique.

The MACE (Major Atmospheric Cherenkov Experiment) telescope has started its regular gamma-ray observations at Hanle in India. Located at an altitude of $\sim$ 4.3 km above sea level and equipped with a 21 m diameter large quasi-parabolic reflector, it has the capability to explore the gamma-ray sky in the energy range above 20 GeV with very high sensitivity. In this work, we present the results from the feasibility studies for searching high-energy gamma-ray signals from dark matter interaction in potential astrophysical environments. We study the impact of MACE response function and other instrumental characteristics to probe the velocity average interaction cross-section ($<\sigma v>$) of the weakly interacting massive particles (WIMPs), expected from the thermal dark matter freeze-out during the decoupling era. We consider the presence of dark matter in the form of pure WIMPs in the mass range 200 GeV - 10 TeV to produce distinctive gamma-ray spectra through its self-annihilation into standard model particles using the Pythia simulation package. The convolution of gamma-ray spectra corresponding to different standard model channels with the MACE response function is used to estimate the upper limit on $<\sigma v>$ for 100 hours of expected MACE observation of Segue1 (a dwarf spheroidal galaxy) which is a potential site of dark matter.

Yttrium (Y), a key s-process element, is commonly used in nucleosynthesis studies and as a Galactic chemical clock when combined with magnesium (Mg). We study the applicability of the previously assumed LTE line formation assumption in Y abundance studies of main-sequence and red giant stars, and probe the impact of NLTE effects on the [Y/Mg] ratio, a proposed stellar age indicator. We derive stellar parameters, ages, and NLTE abundances of Fe, Mg, and Y for 48 solar analogue stars from high-resolution spectra acquired within the Gaia-ESO survey. For Y, we present a new NLTE atomic model. We determine a solar NLTE abundance of A(Y)$_{\rm NLTE}=2.12\pm0.04$ dex, $0.04$ dex higher than LTE. NLTE effects on Y abundance are modest for optical Y II lines, which are frequently used in Sun-like stars diagnostics. NLTE has a small impact on the [Y/Mg] ratio in such stars. For metal-poor red giants, NLTE effects on Y II lines are substantial, potentially exceeding $+0.5$ dex. For the Gaia/4MOST/WEAVE benchmark star, HD 122563, we find the NLTE abundance ratio of [Y/Fe]$_{\rm NLTE}=-0.55\pm0.04$ dex with consistent abundances obtained from different Y II lines. NLTE has a differential effect on Y abundance diagnostics in late-type stars. They notably affect Y II lines in red giants and very metal-poor stars, which are typical Galactic enrichment tracers of neutron-capture elements. For main-sequence stars, NLTE effects on optical diagnostic Y II lines remain minimal across metallicities. This affirms the [Y/Mg] ratio's reliability as a cosmochronometer for Sun-like stars.

We report on exceptionally bright bursts (>400 Jy ms) detected from the repeating fast radio burst source FRB 20220912A using the Nan\c{c}ay Radio Telescope (NRT), as part of the ECLAT (Extragalactic Coherent Light from Astrophysical Transients) monitoring campaign. These bursts exhibit extremely luminous, broadband, short-duration structures (~ 16 microseconds), which we term 'microshots' and which can be especially well studied in the NRT data given the excellent signal-to-noise and dynamic range (32-bit samples). The estimated peak flux density of the brightest microshot is 450 Jy. We show that the microshots are clustered into dense 'forests', by modelling them as Weibull distributions and obtaining Weibull shape parameters of approximately 0.5. Our polarimetric analysis reveals that the bursts are nearly 100% linearly polarised; have < 10% circular polarisation fractions; a near-zero average rotation measure of 0.10(6) rad/m^2; and varying polarisation position angles over the burst duration. For one of the bursts, we analyse raw voltage data from simultaneous observations with the Westerbork RT-1 single 25-m dish. These data allow us to measure the scintillation bandwidth, 0.30(3) MHz, and to probe the bursts on (sub-)microsecond timescales. Some important nuances related to dedispersion are also discussed. We propose that the emission mechanism for the broadband microshots is potentially different from the emission mechanism of the broader burst components which still show a residual drift of a few hundred MHz/ms after correcting for dispersion using the microshots. We discuss how the observed emission is phenomenologically analogous to different types of radio bursts from the Sun.

The Askaryan Radio Array (ARA) is an in-ice ultrahigh energy (UHE, $>10$ PeV) neutrino experiment at the South Pole that aims to detect radio emissions from neutrino-induced particle cascades. ARA has five independent stations which together have collected nearly 24 station-years of data. Each of these stations search for UHE neutrinos by burying in-ice clusters of antennas $\sim 200$ m deep in a roughly cubical lattice with side length $\sim 15$ m. Additionally, the fifth ARA station (A5) has a beamforming trigger, referred to as the Phased Array (PA), consisting of a trigger array of 7 tightly packed vertically-polarized antennas. In this proceeding, we will present a neutrino search with the data of this "hybrid" station, emphasizing its capabilities for improved analysis efficiencies, background rejection, and neutrino vertex reconstruction. This is enabled by combining the closely packed trigger antennas with the long-baselines of the outrigger antennas. We will also place the A5 analysis into the context of the broader five station analysis program, including efforts to characterize and calibrate the detector, model and constrain backgrounds, and reject noise across the entire array. We anticipate this full neutrino search to set world-leading limits above 100 PeV, and inform the next generation of neutrino detection experiments.

NGC 2071 is a star-forming region that overlaps with three $\gamma$-ray sources detected by the Fermi Space Telescope. We propose that strong flare activity in T Tauri stars could produce $\gamma$-ray emission in a way that makes them a counterpart to some unidentified sources detected by the Large Area Telescope aboard the Fermi satellite. We have performed a spectral and temporal analysis for two Fermi data sets: the first 2 yr and the entire 14 yr of observations. We have found that the $\gamma$-ray source is detectable at 3.2$\sigma$ above the background at energies above 100 GeV during the first 2 yr of observation. The analysis of the expected frequency of the highest energy flares occurring in T Tauri stars is consistent with our estimate. In addition, we have determined the minimum energy of the flare that would produce $\gamma$-ray emission, which is $\sim 5 \times 10^{37}$ erg. This agreement becomes a hard observational constraint supporting previous hypotheses about rare flares as the origin of unidentified $\gamma$-ray sources in star-forming regions.

A planet embedded in a protoplanetary disk produces a gap by disk-planet interaction. It also generates velocity perturbation of gas, which can also be observed as deviations from the Keplerian rotation in the channel map of molecular line emission, called kinematic planetary features. These observed signatures provide clues to determine the mass of the planet. We investigated the features induced by the planet with an inclined orbit through three-dimensional hydrodynamic simulations. We found that a smaller planet, with the inclination being $\sim 10^{\circ}$ -- $20^{\circ}$, can produce kinematic features as prominent as those induced by the massive coplanar planet. Despite the kinematic features being similar, the gap is shallower and narrower as compared with the case in which the kinematic features are formed by the coplanar planet. We also found that the kinematic features induced by the inclined planet were fainter for rarer CO isotopologues because the velocity perturbation is weaker at the position closer to the midplane, which was different in the case with a coplanar massive planet. This dependence on the isotopologues is distinguished if the planet has the inclined orbit. We discussed two observed kinematic features in the disk of HD 163296. We concluded that the kink observed at 220 au can be induced by the inclined planet, while the kink at 67 au is consistent to that induced by the coplanar planet.

Blue early-type galaxies with galaxy-scale ongoing star formation are interesting targets in order to understand the stellar mass buildup in elliptical and S0 galaxies in the local Universe. We study the star-forming population of blue early-type galaxies to understand the origin of star formation in these otherwise red and dead stellar systems. The legacy survey imaging data taken with the dark energy camera in the $g$, $r$, and $z$ bands for 55 star-forming blue early-type galaxies were examined, and $g-r$ color maps were created. We identified low surface brightness features near 37 galaxies, faint-level interaction signatures near 15 galaxies, and structures indicative of recent merger activity in the optical color maps of all 55 galaxies. These features are not visible in the shallow Sloan Digital Sky Survey imaging data in which these galaxies were originally identified. Low surface brightness features found around galaxies could be remnants of recent merger events. The star-forming population of blue early-type galaxies could be post-merger systems that are expected to be the pathway for the formation of elliptical galaxies. We hypothesize that the star-forming population of blue early-type galaxies is a stage in the evolution of early-type galaxies. The merger features will eventually disappear, fuel for star formation will cease, and the galaxy will move to the passive population of normal early-type galaxies.

Context. A characteristic feature of young stellar objects is their variability, which is caused by a variety of different physical processes. High-resolution interferometric observations in the near- and mid-infrared wavelength ranges spanning multiple epochs allow the detailed study of these processes. Aims. We aim at investigating the expected variations of the interferometric observables connected to changes in the measured photometric fluxes of a typical variable accreting central young stellar object with a circumstellar disk. Methods. We calculated visibilities and closure phases as well as the photometric flux of brightness distributions obtained using 3D Monte Carlo radiative transfer simulations for a model of a circumstellar disk with an accreting central star. Results. Changes in the accretion luminosity of the central object, that is, an accreting pre-main-sequence star, can lead to significant variations in the visibility and closure phase of the star-disk system measured with instruments at the Very Large Telescope Interferometer (VLTI) that can be related to changes in the photometric flux. Taking into account additional effects due to baseline variation, interferometric observations can provide valuable contributions to the understanding of the underlying processes. Additionally, we provide the web application VLTI B-VAR that allows the impact of the hour angle on the visibility and closure phase for customized intensity maps to be estimated.

Two dark matter searches performed with charge-coupled devices (CCDs) in the DAMIC cryostat at SNOLAB reported with high statistical significance the presence of an unidentified source of low-energy events in bulk silicon. The observed spectrum is consistent with nuclear recoils from the elastic scattering of weakly interacting massive particles (WIMPs) with masses between 2 and 3 GeV. In the standard scenario of spin-independent WIMP-nucleus scattering, the derived cross section is conclusively excluded by results in argon by the DarkSide-50 experiment. We identify isospin-violating and spin-dependent scenarios where interactions with $^{40}$Ar are strongly suppressed and the interpretation of the DAMIC excess as WIMP-nucleus elastic scattering remains viable.

The Be/X-ray binary system LS V +44 17 (RX J0440.9+4431) is a potential member of the rare class of gamma-ray binaries. The system is comprised of a Be star and a neutron star companion with an orbital period of 150 days. In December of 2022, MAXI detected an X-ray outburst from the source, which peaked in early January before declining and then re-brightening. During the second peak, the flux exceeded 1 Crab in the 15-50 keV range, and exhibited a pulsed emission component with a pulse period of 208 seconds. VERITAS observations were conducted close to the peak of the second outburst, from January 24 to January 27, 2023. We report here on the search for very high energy (VHE) gamma-ray emission in these data.

Around 70 $z>6.5$ luminous quasars have been discovered, strongly biased toward the bright end, thus not providing a comprehensive view on quasar abundance beyond cosmic dawn. We present the predicted results of Roman/Rubin high-redshift quasar survey, yielding 3 times more, $2-4$ magnitudes deeper quasar samples, probing high-redshift quasars across broad range of luminosities, especially faint quasars at $L_\mathrm{bol}\sim 10^{10}\;L_{\odot}$ or $M_\mathrm{1450} \sim-22$ that are currently poorly explored. We include high-$z$ quasars, galactic dwarfs and low-$z$ compact galaxies with similar colors as quasar candidates. We create mock catalogs based on population models to evaluate selection completeness and efficiency. We utilize classical color dropout method in $z$ and $Y$ bands to select primary quasar candidates, followed up with Bayesian selection method to identify quasars. We show that overall selection completeness $> 80\%$ and efficiency $\sim 10\%$ at $6.5<z<9$, with 180 quasars at $z>6.5$, 20 at $z > 7.5$ and 2 at $z > 8.5$. The quasar yields depend sensitively on the assumed quasar luminosity shape and redshift evolution. Brown dwarf rejection through proper motion up to 50$\%$ can be made for stars brighter than 25 mag, low-$z$ galaxies dominate at fainter magnitude. Our results show that Roman/Rubin are able to discover a statistical sample of the earliest and faintest quasars in the Universe. The new valuable datasets worth follow up studies with James Webb Space Telescope and Extremely Large Telescopes, to determine quasar luminosity function faint end slope and constraint the supermassive black holes growth in the early Universe.

Our knowledge of the circumgalactic medium (CGM) is mostly based on quasar absorption-line measurements. These have uncovered a multiphase medium that is likely highly turbulent, but constraints of this turbulence are limited to measurements of the non-thermal width of absorption-line components ($b_{turb}$) and the line-of-sight velocity dispersion between components ($\sigma_{LOS}$). Here we analyze a suite of CGM simulations to determine how well these indirect measures are related to the underlying CGM. Our simulations track the non-equilibrium evolution of all commonly observed ions, and consist of two main types: small-scale simulations of regions of homogeneous CGM turbulence and global simulations of inhomogeneous turbulence throughout a galactic halo. From each simulation, we generate mock spectra of Si II, Si IV, C IV, and O VI, which allow us to directly compare $b_{turb}$ and $\sigma_{LOS}$ to the true line-of-sight turbulence ($\sigma_{1D}$). In the small-scale simulations, $b_{turb}$ is only weakly correlated with $\sigma_{1D}$, likely because it measures random motions within individual warm CGM clouds, which do not sample the overall random motions. Meanwhile, $\sigma_{LOS}$ and $\sigma_{1D}$ are strongly correlated, with $\sigma_{1D}\approx\sigma_{LOS}+10$ km s$^{-1}$ in the densest regions we simulated, though, the strength of this correlation depended weakly on the gas phase being probed. Our large-scale simulations also indicate that $b_{turb}$ and $\sigma_{1D}$ are largely uncorrelated, and that $\sigma_{1D}\approx\sigma_{LOS}+10$ kms$^{-1}$ on average, although it varies along individual sightlines. Moreover, the $\sigma_\mathrm{LOS}$ distributions from our global simulations are similar to recent observations, suggesting that this quantity may provide useful constraints on circumgalactic turbulence regardless of the axis probed.

The annual GLEAM outdoor art exhibition features curated, large-scale light installations in the Olbrich Botanical Gardens in Madison, Wisconsin. With submissions from local, regional, and international artists, the annual two-month long event draws tens of thousands of visitors each fall. "Tidal Disruption", an art-science light sculpture representative of the form and behavior of a tidal disruption event, debuted at GLEAM 2022.The installation was a collaboration between artists Hosale and Abdu'Allah, and astrophysicist Madsen. The immersive light and sound exhibit stylistically depicted the death of a star as it spirals into a black hole in a roughly 4 minute sequence. Designed primarily by Hosale, the components were custom fabricated using clear PVC pipe at the Physical Sciences Laboratory at UW-Madison. In this poster, we will describe the process of developing, fabricating, and installing "Tidal Disruption" and the response from the viewers.

When the inflaton is coupled to the gluon Chern-Simons term for successful reheating, mixing between the inflaton and the QCD axion is generally expected given the solution of the strong CP problem by the QCD axion. This is particularly natural if the inflaton is a different, heavier axion. We propose a scenario in which the QCD axion plays the role of the inflaton by mixing with heavy axions. In particular, if the energy scale of inflation is lower than the QCD scale, a hybrid inflation is realized where the QCD axion plays the role of the inflaton in early stages. We perform detailed numerical calculations to take account of the mixing effects. Interestingly, the initial misalignment angle of the QCD axion, which is usually a free parameter, is determined by the inflaton dynamics. It is found to be close to $\pi$ in simple models. This is the realization of the pi-shift inflation proposed in previous literature, and it shows that QCD axion dark matter and inflation can be closely related. The heavy axion may be probed by future accelerator experiments.

It is believed that dark matter (DM) could accumulate inside neutron stars and significantly change their masses, radii and tidal properties. We study what effect bosonic dark matter, modelled as a massive and self-interacting scalar or vector field, has on neutron stars. We derive equations to compute the tidal deformability of the full Einstein-Hilbert-Klein-Gordon system self-consistently, and probe the influence of the scalar field mass and self-interaction strength on the total mass and tidal properties of the combined system, called fermion boson stars (FBS). We are the first to combine Proca stars with neutron stars to mixed systems of fermions and a vector field in Einstein-Proca theory, which we name fermion Proca stars (FPS). We construct equilibrium solutions of FPS, compute their masses, radii and analyse them regarding their stability and higher modes. We find that FPS tend to be more massive and geometrically larger than FBS for equal boson masses and self-interaction strengths. Both FBS and FPS admit DM core and DM cloud solutions and we find that they can produce degenerate results. Core solutions compactify the neutron star component and lower their tidal deformability, cloud solutions have the inverse effect. Electromagnetic observations of certain cloud-like configurations would appear to violate the Buchdahl limit. The self-interaction strength is found to significantly affect both mass and tidal deformability. We discuss observational constraints and the connection to anomalous detections. We also show how models with an effective equation of state compare to the self-consistent solution of FBS and find the self-interaction strength where both solutions converge sufficiently.

Strong lensing of gravitational waves can produce several detectable images as repeated events in the upcoming observing runs, which can be detected with the posterior overlap analysis (Bayes factor). The choice of the binary black hole population plays an important role in the analysis as two gravitational-wave events could be similar either because of lensing or astrophysical coincidence. In this study, we investigate the biases induced by different population models on the Bayes factor. We build up a mock catalogue of gravitational-wave events following a benchmark population and reconstruct it using both non-parametric and parametric methods. Using these reconstructions, we compute the Bayes factor for lensed pair events by utilizing both models and compare the results with a benchmark model. We show that the use of a non-parametric population model gives a smaller bias than parametric population models. Therefore, our study demonstrates the importance of choosing a sufficiently agnostic population model for strong lensing analyses.

Problems with uniform probabilities on an infinite support show up in contemporary cosmology. This paper focuses on the context of inflation theory, where it complicates the assignment of a probability measure over pocket universes. The measure problem in cosmology, whereby it seems impossible to pick out a uniquely well-motivated measure, is associated with a paradox that occurs in standard probability theory and crucially involves uniformity on an infinite sample space. This problem has been discussed by physicists, albeit without reference to earlier work on this topic. The aim of this article is both to introduce philosophers of probability to these recent discussions in cosmology and to familiarize physicists and philosophers working on cosmology with relevant foundational work by Kolmogorov, de Finetti, Jaynes, and other probabilists. As such, the main goal is not to solve the measure problem, but to clarify the exact origin of some of the current obstacles. The analysis of the assumptions going into the paradox indicates that there exist multiple ways of dealing consistently with uniform probabilities on infinite sample spaces. Taking a pluralist stance towards the mathematical methods used in cosmology shows there is some room for progress with assigning probabilities in cosmological theories.