Abstract visibility
Abstract text size

Papers for Wednesday, Jun 28 2023

Papers with local authors

Ore Gottlieb, Danat Issa, Jonatan Jacquemin-Ide, Matthew Liska, Alexander Tchekhovskoy, Francois Foucart, Daniel Kasen, Rosalba Perna, Eliot Quataert, Brian D. Metzger

For movies of the simulations, see this https URL

0 votes
Paper 3 — arXiv:2306.14946
0 votes
Paper 3 — arXiv:2306.14946

The ongoing LIGO-Virgo-KAGRA observing run O4 provides an opportunity to discover new multi-messenger events, including binary neutron star (BNS) mergers such as GW170817, and the highly anticipated first detection of a multi-messenger black hole-neutron star (BH-NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly-relativistic outflows, it has remained uncertain whether the BH-NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly-relativistic outflows in BH-NS mergers. Adopting optimal binary properties, mass ratio of $q=2$ and rapidly rotating BH, we utilize numerical-relativity and general-relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary's evolution from the pre-merger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows, and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of $-14$ to $-15$ a few hours after the merger, longer than previous estimates, which did not consider the first-principles-based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution, or detect its peak up to distances $\gtrsim1$ Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within $ \sim 200 $ Mpc. BH-NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals.

Ore Gottlieb, Danat Issa, Jonatan Jacquemin-Ide, Matthew Liska, Francois Foucart, Alexander Tchekhovskoy, Brian D. Metzger, Eliot Quataert, Rosalba Perna, Daniel Kasen, Matthew D. Duez, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheeli

For movies of the simulations, see this https URL

0 votes
Paper 4 — arXiv:2306.14947
0 votes
Paper 4 — arXiv:2306.14947

We present the first numerical simulations that track the evolution of a black hole-neutron star (BH-NS) merger from pre-merger to $r\gtrsim10^{11}\,{\rm cm}$. The disk that forms after a merger of mass ratio $q=2$ ejects massive disk winds ($3-5\times10^{-2}\,M_\odot$). We introduce various post-merger magnetic configurations, and find that initial poloidal fields lead to jet launching shortly after the merger. The jet maintains a constant power due to the constancy of the large-scale BH magnetic flux, until the disk becomes magnetically arrested (MAD), where the jet power falls off as $L_j\sim t^{-2}$. All jets inevitably exhibit either excessive luminosity due to rapid MAD activation when accretion rate is high, or excessive duration due to delayed MAD activation, compared to typical short gamma-ray burst (sGRBs). This provides a natural explanation to long sGRBs such as GRB 211211A, but also raises a fundamental challenge to our understanding of jet formation in binary mergers. One possible implication being the necessity of higher binary mass ratios or moderate BH spins to launch typical sGRB jets. For post-merger disks with a toroidal magnetic field, dynamo processes delay jet launching such that the jets break out of the disk winds after several seconds. We show for the first time that sGRB jets with initial magnetization $\sigma_0>100$ retain significant magnetization ($\sigma\gg1$) at $r>10^{10}\,{\rm cm}$, emphasizing the importance of magnetic processes in the prompt emission. The jet-wind interaction leads to a power-law angular energy distribution by inflating an energetic cocoon, whose emission is studied in a companion paper.

Zoltán Haiman, Chengcheng Xin, Tamara Bogdanović, Pau Amaro Seoane, Matteo Bonetti, J. Andrew Casey-Clyde, Maria Charisi, Monica Colpi, Jordy Davelaar, Alessandra De Rosa, Daniel J. D'Orazio, Kate Futrowsky, Poshak Gandhi, Alister W. Graham, Jenny E. Greene, Melanie Habouzit, Daryl Haggard, Kelly Holley-Bockelmann, Xin Liu, Alberto Mangiagli, Alessandra Mastrobuono-Battisti, Sean McGee, Chiara M. F. Mingarelli, Rodrigo Nemmen, Antonella Palmese, Delphine Porquet, Alberto Sesana, Aaron Stemo, Alejandro Torres-Orjuela, Jonathan Zrake

White Paper for the Nancy Grace Roman Space Telescope's Core Community Surveys ( this https URL )

0 votes
Paper 15 — arXiv:2306.14990
0 votes
Paper 15 — arXiv:2306.14990

With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12 day orbital periods, which emit copious gravitational radiation and will inevitably merge on timescales of $10^{3-5}$ years. During the last few months of their merger, such binaries are observable with the Laser Interferometer Space Antenna (LISA), a joint ESA/NASA gravitational wave mission set to launch in the mid-2030s. Roman can thus find LISA precursors, provide uniquely robust constraints on the LISA source population, help identify the host galaxies of LISA mergers, and unlock the potential of multi-messenger astrophysics with massive black hole binaries.

Zhou Fan, Gang Zhao, Wei Wang, Jie Zheng, Jingkun Zhao, Chun Li, Yuqin Chen, Haibo Yuan, Haining Li, Kefeng Tan, Yihan Song, Fang Zuo, Yang Huang, Ali Luo, Ali Esamdin, Lu Ma, Bin Li, Nan Song, Frank Grupp, Haibin Zhao, Shuhrat A. Ehgamberdiev, Otabek A. Burkhonov, Guojie Feng, Chunhai Bai, Xuan Zhang, Hubiao Niu, Alisher S. Khodjaev, Bakhodir M. Khafizov, Ildar M. Asfandiyarov, Asadulla M. Shaymanov, Rivkat G. Karimov, Qudratillo Yuldashev, Hao Lu, Getu Zhaori, Renquan Hong, Longfei Hu, Yujuan Liu, Zhijian Xu

49 pages, 21 figures, 5 table, accepted for publication in ApJS

0 votes
Paper 67 — arXiv:2306.15611
0 votes
Paper 67 — arXiv:2306.15611

The Stellar Abundances and Galactic Evolution Survey (SAGES) of the northern sky is a specifically-designed multi-band photometric survey aiming to provide reliable stellar parameters with accuracy comparable to those from low-resolution optical spectra. It was carried out with the 2.3-m Bok telescope of Steward Observatory and three other telescopes. The observations in the $u_s$ and $v_s$ passband produced over 36,092 frames of images in total, covering a sky area of $\sim9960$ degree$^2$. The median survey completeness of all observing fields for the two bands are of $u_{\rm s}=20.4$ mag and $v_s=20.3$ mag, respectively, while the limiting magnitudes with signal-to-noise ratio (S/N) of 100 are $u_s\sim17$ mag and $v_s\sim18$ mag, correspondingly. We combined our catalog with the data release 1 (DR1) of the first of Panoramic Survey Telescope And Rapid Response System (Pan-STARRS1, PS1) catalog, and obtained a total of 48,553,987 sources which have at least one photometric measurement in each of the SAGES $u_s$ and $v_s$ and PS1 $grizy$ passbands, which is the DR1 of SAGES and it will be released in our paper. We compare our $gri$ point-source photometry with those of PS1 and found an RMS scatter of $\sim2$% in difference of PS1 and SAGES for the same band. We estimated an internal photometric precision of SAGES to be on the order of $\sim1$%. Astrometric precision is better than $0^{\prime\prime}.2$ based on comparison with the DR1 of Gaia mission. In this paper, we also describe the final end-user database, and provide some science applications.

All other papers

Neutron star mergers are believed to be a major cosmological source of rapid neutron-capture elements. The kilonovae associated with neutron star mergers have to date yielded only a single well-identified spectral signature: the P Cygni line of Sr$^+$ at about 1$\mu$m in the spectra of the optical transient, AT2017gfo. Such P Cygni lines are important, because they provide significant information not just potentially on the elemental composition of the merger ejecta, but also on the velocity, geometry, and abundance stratification of the explosion. In this paper we show evidence for a previously unrecognised P Cygni line in the spectra of AT2017gfo that emerges several days after the explosion, located at $\lambda \approx 760\,$nm. We show that the feature is well-reproduced by 4d$^2$-4d5p transitions of Y$^+$, which have a weighted mean wavelength around 760-770 nm, with the most prominent line at 788.19 nm. While the observed line is weaker than the Sr$^+$ feature, the velocity stratification of the new line provides an independent constraint on the expansion rate of the ejecta similar to the constraints from Sr$^+$.

Samantha C. Berek, Gwendolyn M. Eadie, Joshua S. Speagle, William E. Harris

16 pages, 9 figures, submitted to ApJ

Galaxy stellar mass is known to be monotonically related to the size of the galaxy's globular cluster (GC) population for Milky Way sized and larger galaxies. However, the relation becomes ambiguous for dwarf galaxies, where there is some evidence for a downturn in GC population size at low galaxy masses. Smaller dwarfs are increasingly likely to have no GCs, and these zeros cannot be easily incorporated into linear models. We introduce the Hierarchical ERrors-in-variables Bayesian lognormAL hurdle (HERBAL) model to represent the relationship between dwarf galaxies and their GC populations, and apply it to the sample of Local Group galaxies where the luminosity range coverage is maximal. This bimodal model accurately represents the two populations of dwarf galaxies: those that have GCs and those that do not. Our model thoroughly accounts for all uncertainties, including measurement uncertainty, uncertainty in luminosity to stellar mass conversions, and intrinsic scatter. The hierarchical nature of our Bayesian model also allows us to estimate galaxy masses and individual mass-to-light ratios from luminosity data within the model. We find that 50% of galaxies are expected to host globular cluster populations at a stellar mass of $\log_{10}(M_*)=6.996$, and that the expected mass of GC populations remains linear down to the smallest galaxies. Our hierarchical model recovers an accurate estimate of the Milky Way stellar mass. Under our assumed error model, we find a non-zero intrinsic scatter of $0.59_{-0.21}^{+0.30}$ (95% credible interval) that should be accounted for in future models.

M71E is a spider pulsar (i.e., a millisecond pulsar with a tight binary companion) with the shortest known orbital period of P=53.3 min discovered by Pan et al. (2023). Their favored evolutionary model suggests that it bridges between two types of spider pulsars, namely, it descended from a "redback" and will become a "black widow". Using Hubble Space Telescope (HST) archival imaging data, we report the first optical identification of its companion COM-M71E. The HST and pulsar timing coordinates are in excellent agreement (within ~10 mas). If M71E is associated with the globular cluster M71, our measured brightness of COM-M71E (m_F606W ~ 25.3) is broadly consistent with the expectation from Pan et al. (2023)'s preferred binary evolutionary model of a stripped dwarf companion.

Naoki Isobe, Hiroshi Nagai, Motoki Kino, Shunsuke Baba, Takao Nakagawa, Yuji Sunada, Makoto Tashiro

10 pages,4 figures, 3 tables, accepted for ApJ

In order to investigate the far-infrared excess detected from the west hot spot of the radio galaxy Pictor A with the Herschel observatory, a submillimeter photometry is performed with the Atacama Compact Array (ACA) of the Atacama Large Millimeter/submillimeter Array at Band 8 with the reference frequency of 405 GHz. A submillimeter source is discovered at the radio peak of the hot spot. Because the 405 GHz flux density of the source, $80.7\pm3.1$ mJy, agrees with the extrapolation of the synchrotron radio spectrum, the far-infrared excess is suggested to exhibit no major contribution at the ACA band. In contrast, by subtracting the power-law spectrum tightly constrained by the radio and ACA data, the significance of the excess in the Herschel band is well confirmed. No diffuse submillimeter emission is detected within the ACA field of view, and thus, the excess is ascribed to the west hot spot itself. In comparison to the previous estimate based on the Herschel data, the relative contribution of the far-infrared excess is reduced by a factor of $\sim 1.5$. The spectrum of the excess below the far-infrared band is determined to be harder than that of the diffusive shock acceleration. This strengthens the previous interpretation that the excess originates via the magnetic turbulence in the substructures within the hot spot. The ACA data are utilized to evaluate the magnetic field strength of the excess and of diffuse radio structure associated to the hot spot.

Thomas Kupfer, Camilla Danielski, Poshak Gandhi, Thomas J. Maccarone, Gijs Nelemans, Valeriya Korol, Liliana Rivera Sandoval

5 pages, 1 figure; Submitted to the NASA Roman Core Community Surveys White Paper Call

Galactic binaries with orbital periods less than 1 hour are strong gravitational wave sources in the mHz regime, ideal for the Laser Interferometer Space Antenna (LISA). At least several hundred, maybe up to a thousand of those binaries are predicted to be sufficiently bright in electromagnetic wavebands to allow detection in both the electromagnetic and the gravitational bands allowing us to perform multi-messenger studies on a statistically significant sample. Theory predicts that a large number of these sources will be located in the Galactic Plane and in particular towards the Galactic Bulge region. Some of these tight binaries may host sub-stellar tertiaries. In this white paper we propose an observing strategy for the Galactic Bulge Time Domain Survey which would use the unique observing capabilities of the Nancy Grace Roman Space telescope to discover and study several 10s of new strong LISA gravitational sources as well as exoplanet candidates around compact white dwarf binaries and other short period variables such as flaring stars, compact pulsators and rotators.

Himanish Ganjoo, M. Sten Delos

33 pages, 16 figures. Comments welcome

In cosmologies with an early matter-dominated era (EMDE) prior to Big Bang nucleosynthesis, the boosted growth of small-scale matter perturbations during the EMDE leads to microhalo formation long before halos would otherwise begin to form. For a range of models, halos can even form during the EMDE itself. These halos would dissipate at the end of the EMDE, releasing their gravitationally heated dark matter and thereby imprinting a free-streaming cut-off on the matter power spectrum. We conduct the first cosmological $N$-body simulations of the formation and evaporation of halos during and after an EMDE. We show that in these scenarios, the free-streaming cut-off after the EMDE can be predicted accurately from the linear matter power spectrum. Although the free streaming can erase much of the EMDE-driven boost to density perturbations, we use our findings to show that the (re-)formation of halos after the EMDE nevertheless proceeds before redshift $\sim 1000$. Early-forming microhalos are a key observational signature of an EMDE, and our prescription for the impact of gravitational heating will allow studies of the observational status and prospects of EMDE scenarios to cover a much wider range of parameters.

D. Benisty, A.-C. Davis, N.W. Evans (Cambridge)

ApJL, submitted

This Letter develops a method to constrain the Cosmological Constant $\Lambda$ from binary galaxies, focusing on the Milky Way and Andromeda. We provide an analytical solution to the two-body problem with $\Lambda$ and show that the ratio between the Keplerian period and $T_\Lambda = 2\pi/(c \sqrt{\Lambda}) \approx 63.2$ Gyr controls the importance of effects from the Cosmological Constant. The Andromeda-Milky Way orbit has a period of $\sim 20$ Gyr and so Dark Energy has to be taken into account. Using the current best mass estimates of the Milky Way and Andromeda galaxies, we find the Cosmological Constant value based only on the Local Group dynamics to be lower then $5.44$ times the value obtained by Planck. With future astrometric measurements, the bound on the Cosmological Constant can be reduced to $\left(1.67 \pm 0.79\right) \Lambda_{\rm PL}$. Our results offer the prospects of constraints on $\Lambda$ over very different scales than previously. The Local Group provides also a completely novel platform to test alternative theories of gravity. We illustrate this by deriving bounds on scalar-tensor theories of gravity over Megaparsec scales.

Kris Pardo, Tzu-Ching Chang, Olivier Doré, Yijun Wang

White paper. 9 pages, 2 figures. Submitted to the Roman Core Community Surveys white paper call

Gravitational waves (GWs) are a new avenue of observing our Universe. So far, we have seen them in the ~10-100 Hz range, and there are hints that we might soon detect them in the nanohertz regime. Multiple efforts are underway to access GWs across the frequency spectrum; however, parts of the frequency space are currently not covered by any planned or future observatories. Photometric surveys can bridge the microhertz gap in the spectrum between LISA and Pulsar Timing Arrays (PTAs) through relative astrometric measurements. Similar to PTA measurements, these astrometric measurements rely on the correlated spacetime distortions produced by gravitational waves at Earth, which induce coherent, apparent stellar position changes on the sky. To detect microhertz GWs with an imaging survey, a combination of high relative astrometric precision, a large number of observed stars, and a high cadence of exposures are needed. Roman's proposed core community survey, the Galactic Bulge Time Domain Survey (RGBTDS), would have all of these components. RGBTDS would be sensitive to GWs with frequencies ranging from $7.7\times 10^{-8}$ Hz to $5.6\times 10^{-4}$ Hz, which opens up a unique GW observing window for supermassive black hole binaries and their waveform evolution. We note that small changes to the survey could enhance Roman's sensitivity to GWs, making it possible to observe the GW background signal that PTAs have recently hinted at with an SNR $\sim$ 70.

Binay Rai, Mohammed Tobrej, Manoj Ghising, Ruchi Tamang, Bikash Chandra Paul

8 pages, 8 figures, accepted for publication in MNRAS

We study the recently discovered Be/X-ray pulsar MAXI J0655--013 using the 2022 \emph{NuSTAR} observations. The paper is the first detailed study of the timing and spectral properties of the source. The pulse profiles of the pulsar vary with energy. The pulsed fraction is found to increase monotonically with energy. In between the two \emph{NuSTAR} observations, a large spin-up rate of $\sim$ -1.23 s d$^{-1}$ is observed, which can be due to large spin-up torque acting on the pulsar during an outburst. Such a large spin-up rate is observed for the first time in an X-ray pulsar during an outburst. The variation of the spin period with time can be employed to obtain the orbital parameters of the binary system, and we found the orbital period to be $\sim$ 27.9 d. The second \emph{NuSTAR} observation is done in a low luminosity state ($L_{X} \sim$3.9$\times$10$^{34}$ \unilum). We have detected the pulsation of the pulsar in such a low luminosity state. In such a low luminosity state, the pulsar MAXI J0655--013 might be accreting from the cold disk.

Akshay Ghosh, L. C. Gallo, A. G. Gonzalez

14 pages, 14 figures, Accepted for publication to MNRAS

In this work, we examine the application of the wavelet transform to the X-ray timing analyses of active galactic nuclei (AGN) and quasi-periodic eruption sources (QPEs). Several scenarios are simulated to test the effectiveness of the wavelet analysis to stationary and non-stationary data. We find that the power spectral density (PSD) slope and the nature of the periodic signal can influence the ability to identify important features in the wavelet power spectrum. In general, weak and transient features can be discerned, which make the wavelet spectrum an important tool in examining AGN light curves. We carried out a wavelet analysis to four unique objects: Ark 120, IRAS 13224-3809, RE J1034+396, and the QPE GSN 069. The well-known quasi-periodic oscillation (QPO) in RE J1034+396 is significantly detected in the wavelet power spectrum. In IRAS 13224-3809, significant transient features appear during a flare at frequencies coincident with previously detected reverberation signals. Finally, the wavelet power spectrum of the QPE GSN 069 significantly reveals four persistent signals that exhibit a 3:2 ratio in oscillation frequencies, consistent with high-frequency QPOs in stellar mass X-ray binaries, but we cannot rule out the possibility this is an artefact of the calculation.

Alex Lazarian, Siyao Xu

10 pages, 2 figures, submitted to ApJ

In a weakly compressible high-$\beta$ medium, pitch-angle scattering and the associated scattering acceleration of cosmic rays (CRs) by anisotropic Alfv\'{e}n and slow modes of magnetohydrodynamic (MHD) turbulence is inefficient. To tap the energy from magnetic compressions for efficient particle acceleration, a diffusion mechanism that can effectively confine particles in space without causing their trapping or pitch-angle isotropization is needed. We find that the mirror diffusion in MHD turbulence recently identified in Lazarian and Xu (2021) satisfies all the above conditions and serves as a promising diffusion mechanism for efficient acceleration of CRs via their stochastic non-resonant interactions with magnetic compressions/expansions. The resulting mirror acceleration is dominated by the slow-mode eddies with their lifetime comparable to the mirror diffusion time of CRs. Consequently, we find that the acceleration time of mirror acceleration is independent of the spatial diffusion coefficient of CRs. The mirror acceleration brings new life for the particle acceleration in a weakly compressible/incompressible medium and has important implications for studying CR re-acceleration in the high-$\beta$ intracluster medium.

Sylvia Biscoveanu, Eric Burns, Philippe Landry, Salvatore Vitale

Submitted to RNAAS

Compact-object binary mergers consisting of one neutron star and one black hole (NSBHs) have long been considered promising progenitors for gamma-ray bursts, whose central engine remains poorly understood. Using gravitational-wave constraints on the population-level NSBH mass and spin distributions we find that at most $20~\mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$ of gamma-ray bursts in the local universe can have NSBH progenitors.

Gokce Ustunisik, Denton S. Ebel, David Walker, Roger L. Nielsen, Marina E. Gemma

23 pages, 15 figures, 5 tables

We determined the mineral-melt partition coefficients (Di's) and the compositional and/or temperature dependency between grossite, melilite, hibonite, olivine and Ca-, Al-inclusion (CAI)-type liquids for a number of light (LE), high field strength (HFSE), large ion lithophile (LILE), and rare earth (REE) elements including Li, Be, B, Sr, Zr, Nb, Ba, La, Ce, Eu, Dy, Ho, Yb, Hf, Ta, Th. A series of isothermal crystallization experiments was conducted at 5 kbar pressure and IW+1 in graphite capsules. The starting compositions were selected based on the calculated and experimentally confirmed phase relations during condensation in CI dust-enriched systems (Ebel and Grossman, 2000; Ebel, 2006; Ustunisik et al., 2014). Partition coefficients between melt and gehlenite, hibonite, and grossite show that the trace element budget of igneous CAIs is controlled by these three major Al-bearing phases in addition to pyroxene. In general, LE, LILE, REE, and HFSE partition coefficients (by mass) decrease in the order of Di(Gehlenite-Melt) > Di(Hibonite-Melt) > Di(Grossite-Melt). Results suggest that Di(Gehlenite-Melt) vary by a factor of 2-3 in different melt compositions at the same T (~1500 C). Increased melt Al and Ca, relative to earlier work, increases the compatibility of Di(Gehlenite-Melt), and also the compatibility of Di(Hibonite-Melt), especially for La and Ce. Olivine partitioning experiments confirm that olivine contribution to the trace element budget of CAIs is small due to the low Di(Olivine-Melt) at a range of temperatures while D-Eu, Yb(Olivine-Melt) are sensitive to changes in T and oxygen fugacity. The development of a predictive model for partitioning in CAI-type systems would require more experimental data and the use of analytical instruments capable of obtaining single phase analyses for crystals < 5 micron.

Andreas Wagner, Emilia K. J. Kilpua, Ranadeep Sarkar, Daniel J. Price, Anshu Kumari, Farhad Daei, Jens Pomoell, Stefaan Poedts

Accepted for publication in Astronomy & Astrophysics

Investigating the early-stage evolution of an erupting flux rope from the Sun is important to understand the mechanisms of how it looses its stability and its space weather impacts. Our aim is to develop an efficient scheme for tracking the early dynamics of erupting solar flux ropes and use the algorithm to analyse its early-stage properties. The algorithm is tested on a data-driven simulation of an eruption that took place in active region AR12473. We investigate the modelled flux rope's footpoint movement and magnetic flux evolution and compare with observational data from the Solar Dynamics Observatory's Atmospheric Imaging Assembly in the 211 $\unicode{x212B}$ and 1600 $\unicode{x212B}$ channels. To carry out our analysis, we use the time-dependent data-driven magnetofrictional model (TMFM). We also perform another modelling run, where we stop the driving of the TMFM midway through the flux rope's rise through the simulation domain and evolve it instead with a zero-beta magnetohydrodynamic (MHD) approach. The developed algorithm successfully extracts a flux rope and its ascend through the simulation domain. We find that the movement of the modelled flux rope footpoints showcases similar trends in both TMFM and relaxation MHD run: they recede from their respective central location as the eruption progresses and the positive polarity footpoint region exhibits a more dynamic behaviour. The ultraviolet brightenings and extreme ultraviolet dimmings agree well with the models in terms of their dynamics. According to our modelling results, the toroidal magnetic flux in the flux rope first rises and then decreases. In our observational analysis, we capture the descending phase of toroidal flux. In conclusion, the extraction algorithm enables us to effectively study the flux rope's early dynamics and derive some of its key properties such as footpoint movement and toroidal magnetic flux.

Jennifer C. Yee, Andrew Gould

10 pages, submitted to Roman Core Community Survey white paper call

Maximizing the scientific return of Roman requires focusing on the scientific discovery space opened up by Roman relative to the ground: i.e., planets in wide orbits (log s > 0.4), the smallest mass-ratio planets (log q < -4.5), and free-floating planet candidates (especially those with thetaE < 1 uas). However, capitalizing on that leverage requires not just detecting such planets but characterizing them sufficiently that they can be used in a statistical analysis. In particular, the signals from all three categories are all prone to light curve degeneracies that may lead to ambiguities in the planet mass-ratio q, separation s, and the size of the source rho (used to measure thetaE and constrain the host mass). Bound planets may also have light curves that are degenerate with models that include a second source rather than a planet. The most immediate need for designing the Roman Galactic Bulge Time Domain Survey is a detailed simulation of wide-orbit and small planetary perturbations to investigate how well the planet perturbations will be characterized. These investigations and related trade-studies must be done in order to maximize Roman's ability to take advantage of new parameter space.

Terrence S. Tricco, Daniel J. Price

7 pages, 5 figures, Proceedings of the "17th SPHERIC International Workshop", Rhodes, Greece, Jun 27-29, 2023

A new implementation for the time evolution of the magnetic vector potential is obtained for smoothed particle magnetohydrodynamics by considering the induction equation in integral form. Galilean invariance is achieved through proper gauge choice. This new discretisation is tested using the Orszag-Tang MHD vortex in a 3D configuration. The corresponding conservative equations of motion are derived, but are not found to solve the MHD equations in the continuum limit. Tests are performed using a hybrid approach instead, whereby the equations of motion based on the magnetic field instead of vector potential are used. Test results experience the same numerical instability as with the Price (2010) formulation. We conclude that this new formulation is non-viable.

Denton S. Ebel

43 pages, 4 color plates, 2 figures

Volatility-dependent fractionation of the rock-forming elements at high temperatures is an early, widespread process during formation of the earliest solids in protoplanetary disks. Equilibrium condensation calculations allow prediction of the identities and compositions of mineral and liquid phases coexisting with gas under presumed bulk chemical, pressure and temperature conditions. A graphical survey of such results is presented for systems of solar and non-solar bulk composition. Chemical equilibrium was approached to varying degrees in the local regions where meteoritic chondrules, Ca-Al-rich inclusions, matrix and other components formed. Early, repeated vapor-solid cycling and homogenization, followed by hierarchical accretion in dust-rich regions, is hypothesized for meteoritic inclusions. Disequilibrium chemical effects appear to have been common at all temperatures, but increasingly so in less refractory meteoritic components. Work is needed to better model high-temperature solid solutions, indicators of these processes.

Travis C. Fischer, Megan C. Johnson, Nathan J. Secrest, D. Michael Crenshaw, Steven B. Kraemer

13 pages, 4 figures, accepted for publication in ApJ

We present recent Very Long Baseline Array (VLBA) 5 GHz radio observations of the nearby, luminous Seyfert 2 galaxy NGC 1068 for comparison to similar VLBA observations made on 1997 April 26. By cross-correlating the positions of emitting regions across both epochs, we find that spatially-resolved extra-nuclear radio knots in this system have sub-relativistic transverse speeds (v < 0.1c). We discuss sources of the observed knots and how the radio emission relates to additional phases of gas in the central ~150 pcs of this system. We suggest that the most likely explanation for the observed emission is synchrotron radiation formed by shocked host media via interactions between AGN winds and the host environment.

J. H. Gillanders, S. A. Sim, S. J. Smartt, S. Goriely, A. Bauswein

Main text: 19 pages, 7 figures, 6 tables. Appendices: 10 pages, 3 tables. Submitted to MNRAS. Comments welcome

Binary neutron star mergers are the first confirmed site of rapid neutron capture (r-process) element nucleosynthesis. The kilonova AT2017gfo is the only electromagnetic counterpart of a neutron star merger spectroscopically observed. We analyse the entire spectral sequence of AT2017gfo (from merger to +10.4 days) and identify seven emission-like features. We confirm that the prominent 1.08um feature can be explained by the Sr II near-infrared triplet evolving from a P-Cygni profile through to pure emission. We calculate the expected strength of the [Sr II] doublet and show that its absence requires highly clumped ejecta. Near-infrared features at 1.58 and 2.07um emerge after three days and become more prominent as the spectra evolve. We model these as optically thick P-Cygni profiles and alternatively as pure emission features (with FWHM = 35600 +/- 6600 km/s), and favour the latter interpretation. The profile of the strong 2.07um emission feature is best reproduced with two lines, centred at 2.059 and 2.135um. We search for candidate ions for all prominent features in the spectra. Strong, permitted transitions of La III, Ce III, Gd III, Ra II and Ac I are plausible candidates for the emission features. If any of these features are produced by intrinsically weak, forbidden transitions, we highlight candidate ions spanning the three r-process peaks. The second r-process peak elements Te and I have plausible matches to multiple features. We highlight the need for more detailed and quantitative atomic line transition data.

Matilde Mingozzi, Bethan L. James, Danielle Berg, Karla Z. Arellano-Córdova, Adele Plat, Claudia Scarlata, Alessandra Aloisi, Ricardo O. Amorín, Jarle Brinchmann, Stéphane Charlot, John Chisholm, Anna Feltre, Simon Gazagnes, Matthew Hayes, Timothy Heckman, Svea Hernandez, Lisa J. Kewley, Nimisha Kumari, Claus Leitherer, Crystal L. Martin, Michael Maseda, Themiya Nanayakkara, Swara Ravindranath, Jane R. Rigby, Peter Senchyna, Evan D. Skillman, Yuma Sugahara, Stephen M. Wilkins, Aida Wofford, Xinfeng Xu

31 pages, submitted to ApJ, comments welcome

In the current JWST era, rest-frame UV spectra play a crucial role in enhancing our understanding of the interstellar medium (ISM) and stellar properties of the first galaxies in the epoch of reionization (EoR, $z>6$). Here, we compare well-known and reliable optical diagrams sensitive to the main ionization source (i.e., star formation, SF; active galactic nuclei, AGN; shocks) to UV counterparts proposed in the literature - the so-called ``UV-BPT diagrams'' - using the HST COS Legacy Archive Spectroscopic SurveY (CLASSY), the largest high-quality, high-resolution and broad-wavelength range atlas of far-UV spectra for 45 local star-forming galaxies. In particular, we explore where CLASSY UV line ratios are located in the different UV diagnostic plots, taking into account state-of-the-art photoionization and shock models and, for the first time, the measured ISM and stellar properties (e.g., gas-phase metallicity, ionization parameter, carbon abundance, stellar age). We find that the combination of C III] $\lambda\lambda$1907,9 He II $\lambda1640$ and O III] $\lambda$1666 can be a powerful tool to separate between SF, shocks and AGN at sub-solar metallicities. We also confirm that alternative diagrams without O III] $\lambda$1666 still allow us to define a SF-locus with some caveats. Diagrams including C IV $\lambda\lambda$1548,51 should be taken with caution given the complexity of this doublet profile. Finally, we present a discussion detailing the ISM conditions required to detect UV emission lines, visible only in low gas-phase metallicity (12+log(O/H) $\lesssim8.3$) and high ionization parameter (log($U$) $\gtrsim-2.5$) environments. Overall, CLASSY and our UV toolkit will be crucial in interpreting the spectra of the earliest galaxies that JWST is currently revealing.

Stefano Gariazzo, Olga Mena

Invited chapter for the edited book Hubble Constant Tension (Editors E. Di Valentino and D. Brout, Springer Singapore, expected in 2024)

Dark radiation, parameterized in terms of $N_{\rm eff}$, has been considered many times in the literature as a possible remedy in alleviating the Hubble constant ($H_0$) tension. We review here the effect of such an extra dark radiation component in the different cosmological observables, focusing mostly on $H_0$. While a larger value of $N_{\rm eff}$ automatically implies a larger value of the Hubble constant, and one would naively expect that such a simple scenario provides a decent solution, more elaborated models are required. Light sterile neutrinos or neutrino asymmetries are among the first-order corrections to the most economical (tree-level) massless dark radiation scenario. However, they are not fully satisfactory in solving the $H_0$ issue. We devote here special attention to second-order corrections: some interacting scenarios, such as those with new dark radiation degrees of freedom that exhibit a non-free streaming nature are highly satisfactory alternative cosmologies where to solve the Hubble constant tension. Models with self-interacting sterile neutrinos and/or majorons, both well-motivated beyond the Standard Model particles, will be discussed along our assessment.

Steph Sallum, Josh Eisner, Andy Skemer, Ruth Murray-Clay

24 pages, 11 figures, accepted for publication in ApJ

We present the highest angular resolution infrared monitoring of LkCa 15, a young solar analog hosting a transition disk. This system has been the subject of a number of direct imaging studies from the millimeter through the optical, which have revealed multiple protoplanetary disk rings as well as three orbiting protoplanet candidates detected in infrared continuum (one of which was simultaneously seen at H$\alpha$). We use high-angular-resolution infrared imaging from 2014-2020 to systematically monitor these infrared signals and determine their physical origin. We find that three self-luminous protoplanets cannot explain the positional evolution of the infrared sources, since the longer time baseline images lack the coherent orbital motion that would be expected for companions. However, the data still strongly prefer a time-variable morphology that cannot be reproduced by static scattered-light disk models. The multi-epoch observations suggest the presence of complex and dynamic substructures moving through the forward-scattering side of the disk at $\sim20$ AU, or quickly-varying shadowing by closer-in material. We explore whether the previous H$\alpha$ detection of one candidate would be inconsistent with this scenario, and in the process develop an analytical signal-to-noise penalty for H$\alpha$ excesses detected near forward-scattered light. Under these new noise considerations, the H$\alpha$ detection is not strongly inconsistent with forward scattering, making the dynamic LkCa 15 disk a natural explanation for both the infrared and H$\alpha$ data.

Warrick H. Ball

Author-produced version of article accepted in Astronomy & Geophysics

Warrick Ball highlights some recent discoveries in the context of the past, present and future of asteroseismology.

Steven Christe, Meriem Alaoui, Joel Allred, Marina Battaglia, Wayne Baumgartner, Juan Camilo Buitrago-Casas, Amir Caspi, Bin Chen, Thomas Chen, Brian Dennis, James Drake, Lindsay Glesener, Iain Hannah, Laura A. Hayes, Hugh Hudson, Andrew Inglis, Jack Ireland, James Klimchuk, Adam Kowalski, Säm Krucker, Anna Maria Massone, Sophie Musset, Michele Piana, Daniel Ryan, Albert Y. Shih, Astrid Veronig, Nicole Vilmer, Alexander Warmuth, Stephen White

White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 14 pages, 4 figures, 1 table

FOXSI is a direct-imaging, hard X-ray (HXR) telescope optimized for solar flare observations. It detects hot plasma and energetic electrons in and near energy release sites in the solar corona via bremsstrahlung emission, measuring both spatial structure and particle energy distributions. It provides two orders of magnitude faster imaging spectroscopy than previously available, probing physically relevant timescales (<1s) never before accessible to address fundamental questions of energy release and efficient particle acceleration that have importance far beyond their solar application (e.g., planetary magnetospheres, flaring stars, accretion disks). FOXSI measures not only the bright chromospheric X-ray emission where electrons lose most of their energy, but also simultaneous emission from electrons as they are accelerated in the corona and propagate along magnetic field lines. FOXSI detects emission from high in the tenuous corona, where previous instruments have been blinded by nearby bright features and will fully characterizes the accelerated electrons and hottest plasmas as they evolve in energy, space, and time to solve the mystery of how impulsive energy release leads to solar eruptions, the primary drivers of space weather at Earth, and how those eruptions are energized and evolve.

N. Werner, J. Řípa, C. Thöne, F. Münz, P. Kurfürst, M. Jelínek, F. Hroch, J. Benáček, M. Topinka, G. Lukes-Gerakopoulos, M. Zajaček, M. Labaj, M. Prišegen, J. Krtička, J. Merc, A. Pál, O. Pejcha, V. Dániel, J. Jon, R. Šošovička, J. Gromeš, J. Václavík, L. Steiger, J. Segiňák, E. Behar, S. Tarem, J.Salh, O. Reich, S. Ben-Ami, M. F. Barschke, D. Berge, A. Tohuvavohu, S. Sivanandam, M. Bulla, S. Popov, Hsiang-Kuang Chang

Submitted to Space Science Reviews

This is the first in a collection of three papers introducing the science with an ultra-violet (UV) space telescope on an approximately 100 kg small satellite with a moderately fast re-pointing capability and a real-time alert communication system that is being studied for a Czech national space mission. The mission, called Quick Ultra-Violet Kilonova surveyor - QUVIK, will provide key follow-up capabilities to increase the discovery potential of gravitational wave observatories and future wide-field multi-wavelength surveys. The primary objective of the mission is the measurement of the UV brightness evolution of kilonovae, resulting from mergers of neutron stars, to distinguish between different explosion scenarios. The mission, which is designed to be complementary to the Ultraviolet Transient Astronomy Satellite - ULTRASAT, will also provide unique follow-up capabilities for other transients both in the near- and far-UV bands. Between the observations of transients, the satellite will target other objects described in this collection of papers, which demonstrates that a small and relatively affordable dedicated UV-space telescope can be transformative for many fields of astrophysics.

J. Krtička, J. Benáček, J. Budaj, D. Korčáková, A. Pál, M. Piecka, M. Zejda, V. Bakış, M. Brož, Hsiang-Kuang Chang, N. Faltová, R. Gális, D. Jadlovský, J. Janík, J. Kára, J. Kolář, I. Krtičková, J. Kubát, B. Kubátová, P. Kurfürst, M. Labaj, J. Merc, Z. Mikulášek, F. Münz, E. Paunzen, M. Prišegen, T. Ramezani, T. Rievajová, J. Řípa, L. Schmidtobreick, M. Skarka, G. Szász, W. Weiss, N. Werner

Submitted to Space Science Reviews

We outline the impact of a small two-band UV-photometry satellite mission on the field of stellar physics, magnetospheres of stars, binaries, stellar clusters, interstellar matter, and exoplanets. On specific examples of different types of stars and stellar systems, we discuss particular requirements for such satellite missions in terms of specific mission parameters such as bandpass, precision, cadence, and mission duration. We show that such a mission may provide crucial data not only for hot stars that emit most of their light in UV, but also for cool stars, where UV traces their activity. This is important, for instance, for exoplanetary studies, because the level of stellar activity influences habitability. While the main asset of the two-band UV mission rests in time-domain astronomy, an example of open clusters proves that such a mission would be important also for the study of stellar populations. Properties of the interstellar dust are best explored when combining optical and IR information with observations in UV. It is well known that dust absorbs UV radiation efficiently. Consequently, we outline how such a UV mission can be used to detect eclipses of sufficiently hot stars by various dusty objects and study disks, rings, clouds, disintegrating exoplanets or exoasteroids. Furthermore, UV radiation can be used to study the cooling of neutron stars providing information about the extreme states of matter in the interiors of neutron stars and used for mapping heated spots on their surfaces.

M. Zajaček, B. Czerny, V. K. Jaiswal, M. Štolc, V. Karas, D. R. Pasham, M. Śniegowska, V. Witzany, P. Suková, F. Münz, N. Werner, J. Řípa, J. Merc, M. Labaj, P. Kurfürst

Submitted to Space Science Reviews

In this review (the third in the series focused on a small two-band UV-photometry mission), we assess possibilities for a small UV two-band photometry mission in studying accreting supermassive black holes (SMBHs; mass range $\sim 10^6$-$10^{10}\,M_{\odot}$). We focus on the following observational concepts: (i) dedicated monitoring of selected type-I Active Galactic Nuclei (AGN) in order to measure the time delay between the far-UV, the near-UV, and other wavebands (X-ray and optical), (ii) nuclear transients including (partial) tidal disruption events and repetitive nuclear transients, and (iii) the study of peculiar sources, such as changing-look AGN, hollows and gaps in accretion disks, low-luminosity AGN, and candidates for Intermediate-Mass Black Holes (IMBHs; mass range $\sim 10^2$-$10^5\,M_{\odot}$) in galactic nuclei. For tidal disruption events (TDEs), high-cadence UV monitoring is crucial for distinguishing among different scenarios for the origin of the UV emission. The small two-band UV space telescope will also provide the information about the near- and far-UV continuum variability for rare transients, such as repetitive partial TDEs and jetted TDEs. We also discuss the possibilities to study and analyze sources with non-standard accretion flows, such as AGN with gappy disks, low-luminosity active galactic nuclei with intermittent accretion, and SMBH binaries potentially involving intermediate-mass black holes.

K. Barkaoui, M. Timmermans, A. Soubkiou, B.V. Rackham, A. J. Burgasser, J. Chouqar, F.J. Pozuelos, K.A. Collins, S.B. Howell, R. Simcoe, C. Melis, K.G. Stassun, J. Tregloan-Reed, M. Cointepas, M. Gillon, X. Bonfils, E. Furlan, C.L. Gnilka, J.M. Almenara, R. Alonso, Z. Benkhaldoun, M. Bonavita, F. Bouchy, A. Burdanov, P. Chinchilla, F. Davoudi, L. Delrez, O. Demangeon, M. Dominik, B.-O. Demory, J. de Wit, G. Dransfield, E. Ducrot, A. Fukui, T. C. Hinse, M.J. Hooton, E. Jehin, J. M. Jenkins, U. G. Jørgensen, D. W. Latham, L. Garcia, S. Carrazco-Gaxiola, M. Ghachoui, Y. Gómez Maqueo Chew, M.N. Günther, J. McCormac, F. Murgas, C. A. Murray, N. Narita, P. Niraula, P. P. Pedersen, D. Queloz, R. Rebolo-López, G. Ricker, L. Sabin, S. Sajadian, N. Schanche, R. P. Schwarz, S. Seager, et al. (11 additional authors not shown)

Accepted for publication in A&A

We present the discovery and validation of two TESS exoplanets orbiting nearby M dwarfs: TOI-2084b, and TOI-4184b. We characterized the host stars by combining spectra from Shane/Kast and Magellan/FIRE, SED (Spectral Energy Distribution) analysis, and stellar evolutionary models. In addition, we used Gemini-South/Zorro & -North/Alopeke high-resolution imaging, archival science images, and statistical validation packages to support the planetary interpretation. We performed a global analysis of multi-colour photometric data from TESS and ground-based facilities in order to derive the stellar and planetary physical parameters for each system. We find that TOI-2084b and TOI-4184b are sub-Neptune-sized planets with radii of Rp = 2.47 +/- 0.13R_Earth and Rp = 2.43 +/- 0.21R_Earth, respectively. TOI-2084b completes an orbit around its host star every 6.08 days, has an equilibrium temperature of T_eq = 527 +/- 8K and an irradiation of S_p = 12.8 +/- 0.8 S_Earth. Its host star is a dwarf of spectral M2.0 +/- 0.5 at a distance of 114pc with an effective temperature of T_eff = 3550 +/- 50 K, and has a wide, co-moving M8 companion at a projected separation of 1400 au. TOI-4184b orbits around an M5.0 +/- 0.5 type dwarf star (Kmag = 11.87) each 4.9 days, and has an equilibrium temperature of T_eq = 412 +/- 8 K and an irradiation of S_p = 4.8 +/- 0.4 S_Earth. TOI-4184 is a metal poor star ([Fe/H] = -0.27 +/- 0.09 dex) at a distance of 69 pc with an effective temperature of T_eff = 3225 +/- 75 K. Both planets are located at the edge of the sub-Jovian desert in the radius-period plane. The combination of the small size and the large infrared brightness of their host stars make these new planets promising targets for future atmospheric exploration with JWST.

Isabella Soh Xiao Si, Michael D. Rhodes, Edwin Budding, Timothy Banks

15 pages, 6 figures, 3 tables. Accepted by JAAVSO

Estimates of orbital parameters were made using a Bayesian optimization technique on astrometric data for 25 visual binary systems catalogued a century ago by the ninth Astronomer Royal, Sir Frank Dyson. An advantage of this method is that it provides reliable, unbiased uncertainty estimates for the optimized parameters. Reasonable agreement is found for the short period (< 100 yr) systems between the current study and Dyson, with superior estimation for the longer systems through the inclusion of an additional century of data. Dynamical masses are presented for the systems through the inclusion of parallax measurements.

Shabbir Bawaji, Ujjaini Alam, Surajit Mondal, Divya Oberoi, Ayan Biswas

Accepted for publication in the Astrophysical Journal

The solar corona is extremely dynamic. Every leap in observational capabilities has been accompanied by unexpected revelations of complex dynamic processes. The ever more sensitive instruments now allow us to probe events with increasingly weaker energetics. A recent leap in the low-frequency radio solar imaging ability has led to the discovery of a new class of emissions, namely Weak Impulsive Narrowband Quiet Sun Emissions \citep[WINQSEs;][]{mondal2020}. They are hypothesized to be the radio signatures of coronal nanoflares and could potentially have a bearing on the long standing coronal heating problem. In view of the significance of this discovery, this work has been followed up by multiple independent studies. These include detecting WINQSEs in multiple datasets, using independent detection techniques and software pipelines, and looking for their counterparts at other wavelengths. This work focuses on investigating morphological properties of WINQSEs and also improves upon the methodology used for detecting WINQSEs in earlier works. We present a machine learning based algorithm to detect WINQSEs, classify them based on their morphology and model the isolated ones using 2D Gaussians. We subject multiple datasets to this algorithm to test its veracity. Interestingly, despite the expectations of their arising from intrinsically compact sources, WINQSEs tend to be resolved in our observations. We propose that this angular broadening arises due to coronal scattering. WINQSEs can, hence, provide ubiquitous and ever-present diagnostic of coronal scattering (and, in turn, coronal turbulence) in the quiet sun regions, which has not been possible till date.

Adolfo S. Carvalho, Lynne A. Hillenbrand, Franz-Josef Hambsch, Shawn Dvorak, Michael Sitko, Ray W. Russell, Victoria Hammond, Michael Connelley, Michael C.B. Ashley, Matthew J. Hankins

15 pages, 13 figures, 2 tables, Accepted to ApJ

We study the evolution of the FU Ori object V960 Mon since its outburst, using available multi-wavelength photometric time series over 8 years, complemented by several epochs of moderate-dispersion spectrophotometry. We find that the source fading can be well-described by a decrease in the temperature of the inner disk, which results from a combination of decreasing accretion rate and increasing inner disk radius. We model the system with a disk atmosphere model that produces the observed variations in multi-band photometry (this paper) and high resolution spectral lines (a companion paper).

Ying Wang, Wei Sun, Ji-lin Zhou, Ming Yang, Fu-yao Liu

16 pages, 5 figures, 3 tables

A possible polar-ring debris disc, the dynamics of which can be described by the outer hierarchical restricted three-body problem, has been detected in 99 Herculis. An empirical formula on the minimum radius beyond which test particles in polar orbits can keep stable within ${10^7}$ binary periods is provided through the numerical fitting, applying to the binary eccentricity $e_{1} \in \left[ {0,0.8} \right)$ and the mass ratio of binary $ \lambda \in \left[ {0.1,1} \right]$, where $ \lambda = m_0/m_1$ (${m_0}$ and ${m_1}$ represent the masses of the two binary stars). The polar planetary disc has the lowerest statistical accretion efficiency and moderate impact frequency of collisions among planetesimals (with a radius of 1-10km) compared to that in the circumbinary coplanar disc and the standard disc around the single host star. Colliding timescales in the circumbinary disk (both polar and coplanar configuration) are longer than $10^7$ yr exceeding the dissipation timescales of the gas disc. The stochastic simulations show that successive collisions cannot make planetesimal grow up which may explain the formation of the debris disc observed in 99 Herculis.

John P. Blakeslee, Michele Cantiello, Michael J. Hudson, Laura Ferrarese, Nandini Hazra, Joseph B. Jensen, Eric W. Peng, Gabriella Raimondo

5 pages; submitted to the call for Roman Core Community Survey white papers

The extragalactic distance scale is fundamental to our understanding of astrophysics and cosmology. In recent years, the surface brightness fluctuation (SBF) method, applied in the near-IR, has proven especially powerful for measuring galaxy distances, first with HST and now with a new JWST program to calibrate the method directly from the tip of the red giant branch (TRGB). So far, however, the distances from space have been gathered slowly, one or two at a time. With the Roman Space Telescope, we have the opportunity to measure uniformly high-quality SBF distances to thousands of galaxies out to hundreds of Mpc. The impact of these data on cosmology and galaxy studies depends on the specifics of the survey, including the filter selection, exposure depth, and (especially) the sky coverage. While the baseline HLWAS survey in four filters plus the grism would yield useful data, the impact would be limited by the relatively small area. A more optimal approach would concentrate on the most efficient passband (F146), adopt an exposure time sufficient to measure good quality distances well out into the Hubble flow, and then maximize the sky coverage within the total time constraints. Grism observations over the same area can provide the needed information on redshifts and spectral energy distributions for compact sources, while colors for larger objects can be obtained from lower resolution surveys. The proposed plan will enable accurate determination of the physical properties of thousands of nearby galaxies, an independent measure of the Hubble constant $H_0$ with negligible statistical error, and competitive constraints on $S_8{\,=\,}\sigma_8(\Omega_m/0.3)^{0.5}$. The resulting data set will be a phenomenal resource for a wide range of studies in astrophysics and cosmology.

Weak gravitational lensing induces flux dependent fluctuations in the observed galaxy number density distribution. This cosmic magnification (magnification bias) effect in principle enables lensing reconstruction alternative to cosmic shear and CMB lensing. However, the intrinsic galaxy clustering, which otherwise overwhelms the signal, has hindered its application. Through a scaling relation found by principal component analysis of the galaxy clustering in multi-band photometry space, we design a minimum variance linear estimator to suppress the intrinsic galaxy clustering and to reconstruct the lensing convergence map. In combination of the CosmoDC2 galaxy mock and the CosmicGrowth simulation, we test this proposal for a LSST-like galaxy survey with $ugrizY$ photometry bands. The scaling relation holds excellently at multipole $\ell<10^3$, and remains reasonably well to $\ell\sim 3000$. The linear estimator efficiently suppresses the galaxy intrinsic clustering, by a factor of $\sim 10^2$. For galaxies in the photo-z range $0.8<z_\kappa<1.2$, the reconstructed convergence map is cosmic variance limited per $\ell$ mode at $\ell<10^2$, and shot noise limited at $\ell>= 200$. Its cross-correlation with cosmic shear of galaxies can achieve $S/N >= 200$. When the source redshift of cosmic shear galaxies $z_\gamma<z_\kappa$, the systematic error is negligible at all investigated scales ($\ell<3000$). When $z_\gamma\geq z_\kappa$, the systematic error caused by the residual intrinsic galaxy clustering becomes non-negligible. We discuss possible mitigation of the residual intrinsic galaxy clustering required for accurate measurement at $\ell>10^3$. This work further demonstrates the potential of lensing measurement through cosmic magnification to enhance the weak lensing cosmology.

Surjit S. Wadhwa, Bojan Arbutina, Nick F.H. Tothill, Miroslav D. Filipovic, Ain Y. De Horta, Jelena Petrovic, Gojko Djurasevic

11 Pages, 2 Figures Accepted Publications of the Astronomical Society of the Pacific (PASP)

Photometric observations and analysis of twelve previously poorly studied contact binary systems is presented. All show total eclipses and have extremely low mass ratios ranging from 0.072 to 0.15. Also, all show characteristics of orbital instability with mass ratios within the theoretical orbital instability range. Although none demonstrate a significant O'Connell effect at least nine of the systems have other indicators of increased chromospheric and magnetic activity.

Kovi Rose, Joshua Pritchard, Tara Murphy, Manisha Caleb, Dougal Dobie, Laura Driessen, Stefan W. Duchesne, David L. Kaplan, Emil Lenc, Ziteng Wang

Accepted for publication in ApJ Letters; 11 pages, 3 figures and 2 tables

We present the detection of rotationally modulated, circularly polarized radio emission from the T8 brown dwarf WISE J062309.94-045624.6 between 0.9 and 2.0 GHz. We detected this high proper motion ultracool dwarf with the Australian SKA Pathfinder in $1.36$ GHz imaging data from the Rapid ASKAP Continuum Survey. We observed WISE J062309.94-045624.6 to have a time and frequency averaged Stokes I flux density of $4.17\pm0.41$ mJy beam$^{-1}$, with an absolute circular polarization fraction of $66.3\pm9.0\%$, and calculated a specific radio luminosity of $L_{\nu}\sim10^{14.8}$ erg s$^{-1}$ Hz$^{-1}$. In follow-up observations with the Australian Telescope Compact Array and MeerKAT we identified a multi-peaked pulse structure, used dynamic spectra to place a lower limit of $B>0.71$ kG on the dwarf's magnetic field, and measured a $P=1.912\pm0.005$ h periodicity which we concluded to be due to rotational modulation. The luminosity and period we measured are comparable to those of other ultracool dwarfs observed at radio wavelengths. This implies that future megahertz to gigahertz surveys, with increased cadence and improved sensitivity, are likely to detect similar or later-type dwarfs. Our detection of WISE J062309.94-045624.6 makes this dwarf the coolest and latest-type star observed to produce radio emission.

Tetiana Kozynets, Anatoli Fedynitch, D. Jason Koskinen

25 pages, 17 figures

The atmospheric lepton fluxes play a crucial role in many particle and astroparticle physics experiments, e.g. in establishing the neutrino signal and the muon background for neutrino oscillation measurements, or the atmospheric background for astrophysical neutrino searches. The Matrix Cascade Equations (MCEq) code is a numerical tool used to model the atmospheric lepton fluxes by solving a system of coupled differential equations for particle production, interaction, and decay at extremely low computational costs. Previously, the MCEq framework only accommodated longitudinal development of air showers, an approximation that works well for neutrino and muon fluxes at high energies (O(10 GeV) and above). However, for accurate calculations of atmospheric lepton angular distributions at lower energies, the lateral component of hadronic cascades becomes significant, necessitating three-dimensional calculation schemes. We introduce "2D MCEq", an efficient numerical approach for combined longitudinal and angular evolution of air showers that retains the low computational complexity. The accuracy of the "2D MCEq" is affirmed by its benchmark comparison with the standard Monte Carlo code CORSIKA. This study paves the way for efficient three-dimensional calculations of atmospheric neutrino fluxes.

Noam Soker (Technion, Israel)

Will be submitted in two days to allow for comments

I present the effervescent zone model to account for the compact dense circumstellar material (CSM) around the progenitor of the core collapse supernova (CCSN) SN 2023ixf. The effervescent zone is composed of bound dense clumps that are lifted by stellar pulsation and envelope convection to distances of tens AUs, and then fall back. The dense clumps provide most of the compact CSM mass and exist alongside the regular (escaping) wind. I crudely estimate that for a compact CSM within ~30 AU that contains ~0.01 Mo, the density of each clump is >3000 times the density of the regular wind at the same radius and that the total volume filling factor of the clumps is several percent. The clumps might cover only a small fraction of the CCSN photosphere in the first days post-explosion, accounting for the lack of strong narrow absorption lines. The long-lived effervescent zone is compatible with no evidence for outbursts in the years prior to SN 2023ixf explosion and the large-amplitude pulsations of its progenitor, and it is an alternative to the CSM scenario of several-years-long high mass loss rate wind.

Vlad D. Berlea, Steven Worm, Nirmal Kaipachery, Shrinivasrao R. Kulkarni, Shashank Kumar, Merlin F. Barschke, David Berge, Adi Birman, Shay Alfassi, Amos Fenigstein

Part of the conference: Frontier Detectors for Frontier Physics: 15th Pisa Meeting on Advanced Detectors, La Biodola - Isola d'Elba Published in: Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment Available online 15 June 2023, 168463. In Press, Journal Pre-proof

ULTRASAT (ULtraviolet TRansient Astronomy SATellite) is a wide-angle space telescope that will perform deep time-resolved surveys in the near-ultraviolet spectrum. ULTRASAT is a space mission led by the Weizmann Institute of Science and the Israel Space Agency and is planned for launch in 2025. The camera implements backside-illuminated, stitched pixel sensors. The pixel has a dual-conversion-gain 4T architecture, with a pitch of $9.5$ $\mu m$ and is produced in a $180$ $nm$ process by Tower Semiconductor. Before the final sensor was available for testing, test sensors provided by Tower were used to gain first insights into the pixel's radiation tolerance. One of the main contributions to sensor degradation due to radiation for the ULTRASAT mission is Total Ionizing Dose (TID). TID measurements on the test sensors have been performed with a Co-60 gamma source at Helmholz Zentrum Berlin and CC-60 facility at CERN and preliminary results are presented.

A. Asensio Ramos, M. C. M. Cheung, I. Chifu, R. Gafeira

100 pages, 13 figures, 286 references, accepted for publication as a Living Review in Solar Physics (LRSP)

The application of machine learning in solar physics has the potential to greatly enhance our understanding of the complex processes that take place in the atmosphere of the Sun. By using techniques such as deep learning, we are now in the position to analyze large amounts of data from solar observations and identify patterns and trends that may not have been apparent using traditional methods. This can help us improve our understanding of explosive events like solar flares, which can have a strong effect on the Earth environment. Predicting hazardous events on Earth becomes crucial for our technological society. Machine learning can also improve our understanding of the inner workings of the sun itself by allowing us to go deeper into the data and to propose more complex models to explain them. Additionally, the use of machine learning can help to automate the analysis of solar data, reducing the need for manual labor and increasing the efficiency of research in this field.

Zhenxin Lei, Ruijie He, Peter Nemeth, Xuan Zou, Huaping Xiao, Yong Yang, Jingkun Zhao

38 pages, 13 figures, 3 tables, accepted for publication in ApJ

Masses for 664 single-lined hot subdwarf stars identified in LAMOST were calculated by comparing synthetic fluxes from spectral energy distribution (SED) with observed fluxes from virtual observatory service. Three groups of hot subdwarf stars were selected from the whole sample according to their parallax precision to study the mass distributions. We found, that He-poor sdB/sdOB stars present a wide mass distribution from 0.1 to 1.0 $\mathrm{M}_{\odot}$ with a sharp mass peak around at 0.46 $\rm{M}_{\odot}$, which is consistent with canonical binary model prediction. He-rich sdB/sdOB/sdO stars present a much flatter mass distribution than He-poor sdB/sdOB stars and with a mass peak around 0.42 $\mathrm{M}_{\odot}$. By comparing the observed mass distributions to the predictions of different formation scenarios, we concluded that the binary merger channel, including two helium white dwarfs (He-WDs) and He-WD + main sequence (MS) merger, cannot be the only main formation channel for He-rich hot subdwarfs, and other formation channels such as the surviving companions from type Ia supernovae (SNe Ia) could also make impacts on producing this special population, especially for He-rich hot subdwarfs with masses less than 0.44 $\mathrm{M}_{\odot}$. He-poor sdO stars also present a flatter mass distribution with an inconspicuous peak mass at 0.18 $\mathrm{M}_{\odot}$. The similar mass - $\Delta RV_\mathrm{max}$ distribution between He-poor sdB/sdOB and sdO stars supports the scenario that He-poor sdO stars could be the subsequent evolution stage of He-poor sdB/sdOB stars.

Somnath Dutta, Chin-Fei Lee, Doug Johnstone, Jeong-Eun Lee, Naomi Hirano, James Di Francesco, Anthony Moraghan, Tie Liu, Dipen Sahu, Sheng-Yuan Liu, Kenichi Tatematsu, Chang Won Lee, Shanghuo Li, David Eden, Mika Juvela, Leonardo Bronfman, Shih-Ying Hsu, Kee-Tae Kim, Woojin Kwon, Patricio Sanhueza, Jesus Alejandro Lopez-Vazquez, Qiuyi Luo, Hee-Weon Yi

Submitted to Journal; 27 pages, 15 Figures and additional Appendix material

Protostellar outflows and jets are almost ubiquitous characteristics during the mass accretion phase, and encode the history of stellar accretion, complex-organic molecule (COM) formation, and planet formation. Episodic jets are likely connected to episodic accretion through the disk. Despite the importance, there is a lack of studies of a statistically significant sample of protostars via high-sensitivity and high-resolution observations. To explore episodic accretion mechanisms and the chronologies of episodic events, we investigated 42 fields containing protostars with ALMA observations of CO, SiO, and 1.3\,mm continuum emission. We detected SiO emission in 21 fields, where 19 sources are driving confirmed molecular jets with high abundances of SiO. Jet velocities, mass-loss rates, mass-accretion rates, and periods of accretion events are found to be dependent on the driving forces of the jet (e.g., bolometric luminosity, envelope mass). Next, velocities and mass-loss rates are positively correlated with the surrounding envelope mass, suggesting that the presence of high mass around protostars increases the ejection-accretion activity. We determine mean periods of ejection events of 20$-$175 years for our sample, which could be associated with perturbation zones of $\sim$ 2$-$25\,au extent around the protostars. Also, mean ejection periods are anti-correlated with the envelope mass, where high-accretion rates may trigger more frequent ejection events. The observed periods of outburst/ejection are much shorter than the freeze-out time scale of the simplest COMs like CH$_3$OH, suggesting that episodic events largely maintain the ice-gas balance inside and around the snowline.

S. Koc, T. Yontan

7 pages, 6 figures and 3 tables, accepted for publication in Bitlis Eren University Journal of Science

In this study, the structural and basic astrophysical parameters of the poorly studied open cluster Berkeley 6 are calculated. Analyses of the cluster are carried out using the third photometric, spectroscopic, and astrometric data release of Gaia (Gaia DR3). The membership probabilities of stars located in the direction of the cluster region are calculated by considering their astrometric data. Thus, we identified 119 physical members for Berkeley 6. The colour excess, distance, and age of the cluster are determined simultaneously on the colour-magnitude diagram. We fitted solar metallicity PARSEC isochrones to the colour-magnitude diagram by considering the most probable member stars and obtained $E(G_{\rm BP}-G_{\rm RP})$ colour excess as 0.918$\pm$0.145 mag. The distance and age of the cluster are determined as $d=2625\pm337$ pc and $t=350\pm50$ Myr, respectively.

S. Bellotti, R. Fares, A. A. Vidotto, J. Morin, P. Petit, G. A. J. Hussain, V. Bourrier, J.F. Donati, C. Moutou, E. Hebrard

9 pages, 8 figures, 3 tables

The space environment in which planets are embedded depends mainly on the host star and impacts the evolution of the planetary atmosphere. The quiet M dwarf GJ 436 hosts a close-in hot Neptune which is known to feature a comet-like tail of hydrogen atoms escaped from its atmosphere due to energetic stellar irradiation. Understanding such star-planet interactions is essential to shed more light on planet formation and evolution theories, in particular the scarcity of Neptune-size planets below 3 d orbital period, also known as ``Neptune desert''. We aimed at characterising the stellar environment around GJ 436, which requires an accurate knowledge of the stellar magnetic field. The latter is studied efficiently with spectropolarimetry, since it is possible to recover the geometry of the large-scale magnetic field by applying tomographic inversion on time series of circularly polarised spectra. We used spectropolarimetric data collected in the optical domain with Narval in 2016 to compute the longitudinal magnetic field, examine its periodic content via Lomb-Scargle periodogram and Gaussian Process Regression analysis, and finally reconstruct the large-scale field configuration by means of Zeeman-Doppler Imaging. We found an average longitudinal field of -12 G and a stellar rotation period of 46.6 d using a Gaussian Process model and 40.1 d using Zeeman-Doppler Imaging, both consistent with the literature. The Lomb-Scargle analysis did not reveal any significant periodicity. The reconstructed large-scale magnetic field is predominantly poloidal, dipolar and axisymmetric, with a mean strength of 16 G. This is in agreement with magnetic topologies seen for other stars of similar spectral type and rotation rate.

A. Tolamatti, K. K. Singh, K. K. Yadav

16 Pages, 10 Figures, Accepted for Publication in MNRAS

Machine learning based approaches are emerging as very powerful tools for many applications including source classification in astrophysics research due to the availability of huge high quality data from different surveys in observational astronomy. The Large Area Telescope on board \emph{Fermi} satellite (\emph{Fermi}-LAT) has discovered more than 6500 high energy gamma-ray sources in the sky from its survey over a decade. A significant fraction of sources observed by the \emph{Fermi}-LAT either remains unassociated or has been identified as \emph{Blazar Candidates of Uncertain type} (BCUs). We explore the potential of eXtreme Gradient Boosting (XGBoost)- a supervised machine learning algorithm to identify the blazar subclasses among a sample of 112 BCUs of the 4FGL catalog whose X-ray counterparts are available within 95$\%$ uncertainty regions of the \emph{Fermi}-LAT observations. We have used information from the multi-wavelength observations in IR, optical, UV, X-ray and $\gamma$-ray wavebands along with the redshift measurements reported in the literature for classification. Among the 112 uncertain type blazars, 62 are classified as BL Lacertae objects (BL Lacs) and 6 have been classified as Flat Spectrum Radio Quasars (FSRQs). This indicates a significant improvement with respect to the multi-perceptron neural network based classification reported in the literature. Our study suggests that the gamma-ray spectral index, and IR color indices are the most important features for identifying the blazar subclasses using the \emph{XGBoost} classifier. We also explore the importance of redshift in the classification BCU candidates.

Vikrant V. Jadhav (Uni. Bonn), Annapurni Subramaniam (IIA), Ram Sagar (IIA)

8 pages, 2 figures, catalogue available in CDS, accepted in A&A

NGC 6791 is one of the richest old open clusters in the Milky Way. Its position above the Galactic plane and the number density makes it an interesting middle ground between Galactic open and globular clusters. We aim to detect the UV bright population of NGC 6791 using \textit{AstroSat}/UVIT images in near-UV and far-UV filters and characterise the known post mass transfer systems such as blue straggler stars (BSSs). We identified 20 members with large UV flux (out of 91 cluster members among 1180 detections), suggestive of binarity, interactions or stellar activity using multi-wavelength spectral energy distribution analysis. We characterised 62 isolated cluster members, including five hot subdwarfs (sdA/sdB). Additionally, we detected ten sdA/sdB/extremely low mass white dwarf (ELM) type candidates hidden alongside other cluster members. Additionally, we report the discovery of four candidate blue lurkers, which are main sequence stars with mass accretion history. We report that this cluster has a variety of stellar (pre-)remnants, such as sdBs, sdAs, and ELM white dwarfs, that are by-products of binary evolution. The above are likely to be post mass transfer binaries found throughout the evolutionary phases from the main sequence to the post horizontal branch. Therefore, this dynamically old open cluster is unique, making it an ideal testbed for dynamical studies.

Camilla Giusti, Mario Cadelano, Francesco R. Ferraro, Barbara Lanzoni, Silvia Leanza, Cristina Pallanca, Enrico Vesperini, Emanuele Dalessandro, Alessio Mucciarelli

Accepted for publication in ApJ

In the context of a project aimed at characterizing the dynamical evolution of old globular clusters in the Large Magellanic Cloud, we have secured deep HST/WFC3 images of the massive cluster NGC 1835. In the field of view of the acquired images, at a projected angular separation of approximately 2 arcmin from the cluster, we detected the small stellar system KMK88-10. The observations provided the deepest color-magnitude diagram ever obtained for this cluster, revealing that it hosts a young stellar population with an age of 600-1000 Myr. The cluster surface brightness profile is nicely reproduced by a King model with a core radius rc = 4 arcsec (0.97 pc), an half-mass radius rhm = 12 arcsec (2.9 pc), and a concentration parameter c~1.3 corresponding to a truncation radius rt~81 arcsec (19.5 pc). We also derived its integrated absolute magnitude (MV=-0.71) and total mass (M~80-160 Msun). The most intriguing feature emerging from this analysis is that KMK88-10 presents a structure elongated in the direction of NGC 1835, with an intracluster over-density that suggests the presence of a tidal bridge between the two systems. If confirmed, this would be the first evidence of a tidal capture of a small star cluster by a massive globular.

Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Patrick Sheehan, Yuri Aikawa, Zhi-Yun Li, Leslie W. Looney, Jonathan P. Willians, Yusuke Aso, Rajeeb Sharma, Jinshi Sai (Insa Choi), Yoshihide Yamato, Jeong-Eun Lee, Kengo Tomida, Hsi-Wei Yen, Frankie J Encalada, Christian Flores, Sacha Gavino, Miyu Kido, Ilseung Han, Zhe-Yu Daniel Lin, Suchitra Narayanan, Nguyen Thi Phuong, Alejandro Santamaría-Miranda, Travis J. Thieme, Merel L. R. van 't Hoff, Itziar de Gregorio-Monsalvo, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Adele Plunkett, Kazuya Saigo, Shingo Hirano, Ka Ho Lam, Shoji Mori

This is a publication of a series of eDisk ALMA large program first-look papers

We present an overview of the Large Program, ``Early Planet Formation in Embedded Disks (eDisk)'', conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The ubiquitous detections of substructures, particularly rings and gaps, in protoplanetary disks around T Tauri stars raise the possibility that at least some planet formation may have already started during the embedded stages of star formation. In order to address exactly how and when planet formation is initiated, the program focuses on searching for substructures in disks around 12 Class 0 and 7 Class I protostars in nearby ($< $200 pc) star-forming regions through 1.3 mm continuum observations at a resolution of $\sim7$ au (0.04"). The initial results show that the continuum emission, mostly arising from dust disks around the sample protostars, has relatively few distinctive substructures, such as rings and spirals, in marked contrast to Class II disks. The dramatic difference may suggest that substructures quickly develop in disks when the systems evolve from protostars to Class II sources or alternatively that high optical depth of the continuum emission could obscure internal structures. Kinematic information obtained through CO isotopologue lines and other lines reveals the presence of Keplerian disks around protostars, providing us with crucial physical parameters, in particular, the dynamical mass of the central protostars. We describe the background of the eDisk program, the sample selection and their ALMA observations, the data reduction, and also highlight representative first-look results.

Merel L.R. van 't Hoff, John J. Tobin, Zhi-Yun Li, Nagayoshi Ohashi, Jes K. Jørgensen, Zhe-Yu Daniel Lin, Yuri Aikawa, Yusuke Aso, Itziar de Gregorio-Monsalvo, Sacha Gavino, Ilseung Han, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Leslie W. Looney, Suchitra Narayanan, Adele Plunkett, Jinshi Sai (Insa Choi), Alejandro Santamaría-Miranda, Rajeeb Sharma, Patrick D. Sheehan, Shigehisa Takakuwa, Travis J. Thieme, Jonathan P. Williams, Shih-Ping Lai, Nguyen Thi Phuong, Hsi-Wei Yen

27 pages, 16 figures, 2 tables, 10 pages appendix with 12 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

Studying the physical and chemical conditions of young embedded disks is crucial to constrain the initial conditions for planet formation. Here, we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum at $\sim$0.06" (8 au) resolution and molecular line emission at $\sim$0.17" (24 au) resolution toward the Class 0 protostar L1527 IRS from the Large Program eDisk (Early Planet Formation in Embedded Disks). The continuum emission is smooth without substructures, but asymmetric along both the major and minor axes of the disk as previously observed. The detected lines of $^{12}$CO, $^{13}$CO, C$^{18}$O, H$_2$CO, c-C$_3$H$_2$, SO, SiO, and DCN trace different components of the protostellar system, with a disk wind potentially visible in $^{12}$CO. The $^{13}$CO brightness temperature and the H$_2$CO line ratio confirm that the disk is too warm for CO freeze out, with the snowline located at $\sim$350 au in the envelope. Both molecules show potential evidence of a temperature increase around the disk-envelope interface. SO seems to originate predominantly in UV-irradiated regions such as the disk surface and the outflow cavity walls rather than at the disk-envelope interface as previously suggested. Finally, the continuum asymmetry along the minor axis is consistent with the inclination derived from the large-scale (100" or 14,000 au) outflow, but opposite to that based on the molecular jet and envelope emission, suggesting a misalignment in the system. Overall, these results highlight the importance of observing multiple molecular species in multiple transitions to characterize the physical and chemical environment of young disks.

Yoshihide Yamato, Yuri Aikawa, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Yusuke Aso, Jinshi Sai (Insa Choi), Christian Flores, Itziar de Gregorio-Monsalvo, Shingo Hirano, Ilseung Han, Miyu Kido, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Shoji Mori, Suchitra Narayanan, Nguyen Thi Phuong, Kazuya Saigo, Alejandro Santamaría-Miranda, Rajeeb Sharma, Travis J. Thieme, Kengo Tomida, Merel L. R. van 't Hoff, Hsi-Wei Yen

24 pages, 12 figures. Accepted for publication in The Astrophysical Journal as one of the first-look papers of the eDisk ALMA Large Program

Constraining the physical and chemical structure of young embedded disks is crucial to understanding the earliest stages of planet formation. As part of the Early Planet Formation in Embedded Disks Atacama Large Millimeter/submillimeter Array Large Program, we present high spatial resolution ($\sim$0$.\!\!^{\prime\prime}$1 or $\sim$15 au) observations of the 1.3 mm continuum and $^{13}$CO $J=$ 2-1, C$^{18}$O $J=$ 2-1, and SO $J_N=$ $6_5$-$5_4$ molecular lines toward the disk around the Class I protostar L1489 IRS. The continuum emission shows a ring-like structure at 56 au from the central protostar and a tenuous, optically thin emission extending beyond $\sim$300 au. The $^{13}$CO emission traces the warm disk surface, while the C$^{18}$O emission originates from near the disk midplane. The coincidence of the radial emission peak of C$^{18}$O with the dust ring may indicate a gap-ring structure in the gaseous disk as well. The SO emission shows a highly complex distribution, including a compact, prominent component at $\lesssim$30 au, which is likely to originate from thermally sublimated SO molecules. The compact SO emission also shows a velocity gradient along a slightly ($\sim15^\circ$) tilted direction with respect to the major axis of the dust disk, which we interpret as an inner warped disk in addition to the warp around $\sim$200 au suggested by previous work. These warped structures may be formed by a planet or companion with an inclined orbit, or by a gradual change in the angular momentum axis during gas infall.

Zhe-Yu Daniel Lin, Zhi-Yun Li, John J. Tobin, Nagayoshi Ohashi, Jes Kristian Jørgensen, Leslie W. Looney, Yusuke Aso, Shigehisa Takakuwa, Yuri Aikawa, Merel L. R. van 't Hoff, Itziar de Gregorio-Monsalvo, Frankie J. Encalada, Christian Flores, Sacha Gavino, Ilseung Han, Miyu Kido, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Nguyen Thi Phuong, Jinshi Sai (Insa Choi), Rajeeb Sharma, Patrick Sheehan, Travis J. Thieme, Jonathan P. Williams, Yoshihide Yamato, Hsi-Wei Yen

33 pages, 21 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

While dust disks around optically visible, Class II protostars are found to be vertically thin, when and how dust settles to the midplane are unclear. As part of the Atacama Large Millimeter/submillimeter Array (ALMA) large program, Early Planet Formation in Embedded Disks, we analyze the edge-on, embedded, Class I protostar IRAS 04302+2247, also nicknamed the ``Butterfly Star." With a resolution of 0.05" (8~au), the 1.3 mm continuum shows an asymmetry along the minor axis which is evidence of an optically thick and geometrically thick disk viewed nearly edge-on. There is no evidence of rings and gaps, which could be due to the lack of radial substructure or the highly inclined and optically thick view. With 0.1" (16~au) resolution, we resolve the 2D snow surfaces, i.e., the boundary region between freeze-out and sublimation, for $^{12}$CO $J$=2--1, $^{13}$CO $J$=2--1, C$^{18}$O $J$=2--1, $H_{2}$CO $J$=$3_{0,3}$--$2_{0,2}$, and SO $J$=$6_{5}$--$5_{4}$, and constrain the CO midplane snow line to $\sim 130$ au. We find Keplerian rotation around a protostar of $1.6 \pm 0.4 M_{\odot}$ using C$^{18}$O. Through forward ray-tracing using RADMC-3D, we find that the dust scale height is $\sim 6$ au at a radius of 100~au from the central star and is comparable to the gas pressure scale height. The results suggest that the dust of this Class~I source has yet to vertically settle significantly.

Miyu Kido, Shigehisa Takakuwa, Kazuya Saigo, Nagayoshi Ohashi, John J. Tobin, Jes K, Jørgensen, Yuri Aikawa, Yusuke Aso, Frankie J. Encalada, Christian Flores, Sacha Gavino, Itziar de Gregorio-Monsalvo, Ilseung Han, Shingo Hirano, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W.Looney, Shoji Mori, Suchitra Narayanan, Adele L. Plunkett, Nguyen Thi Phuong, Jinshi Sai (Insa Choi), Alejandro Santamarîa-Miranda, Rajeeb Sharma, Patrick Sheehan, Travis J. Thieme, Kengo Tomida, Merel L.R. van't Hoff, Jonathan P. Williams, Yoshihide Yamato, Hsi-Wei Yen

30 pages, 24 figures, accepted for publication in The Astrophysical Journal as one of the first-look papers of the eDisk ALMA Large Program

We present observations of the Class 0 protostar IRAS 16544-1604 in CB 68 from the ''Early Planet Formation in Embedded Disks (eDisk)'' ALMA Large program. The ALMA observations target continuum and lines at 1.3-mm with an angular resolution of $\sim$5 au. The continuum image reveals a dusty protostellar disk with a radius of $\sim$30 au seen close to edge-on, and asymmetric structures both along the major and minor axes. While the asymmetry along the minor axis can be interpreted as the effect of the dust flaring, the asymmetry along the major axis comes from a real non-axisymmetric structure. The C$^{18}$O image cubes clearly show the gas in the disk that follows a Keplerian rotation pattern around a $\sim$0.14 $M_{\odot}$ central protostar. Furthermore, there are $\sim$1500 au-scale streamer-like features of gas connecting from North-East, North-North-West, and North-West to the disk, as well as the bending outflow as seen in the $^{12}$CO (2-1) emission. At the apparent landing point of NE streamer, there are SO (6$_5$-5$_4$) and SiO (5-4) emission detected. The spatial and velocity structure of NE streamer can be interpreted as a free-falling gas with a conserved specific angular momentum, and the detection of the SO and SiO emission at the tip of the streamer implies presence of accretion shocks. Our eDisk observations have unveiled that the Class 0 protostar in CB 68 has a Keplerian rotating disk with flaring and non-axisymmetric structure associated with accretion streamers and outflows.

Yoshiaki Sofue, Jun Kataoka, Ryoji Iwashita

To appear in MNRAS, 8 pages, 8 figures

The ratio of the H$\alpha$ and radio continuum intensities in the North Polar Spur (NPS) is measured to be $\lesssim 50$, two orders of magnitude smaller than the values observed in the typical shell-type old supernova remnants (SNRs), Cygnus Loop and S147, of $\sim 10^4$.The extremely low} H$\alpha$-to-radio intensity ratio favors the GC explosion model}, which postulates a giant shock wave in the hot and low-density Galactic halo with low hydrogen recombination rate, over the local supernova(e) remnant model.

Jose Schiappacasse-Ulloa, Mario Pasquato, Sara Lucatello

10 pages, 6 figures. Submitted to the International Conference on Machine Learning 2023

The upcoming facilities like the Vera C. Rubin Observatory will provide extremely deep photometry of thousands of star clusters to the edge of the Galaxy and beyond, which will require adequate tools for automatic analysis, capable of performing tasks such as the characterization of a star cluster through the analysis of color-magnitude diagrams (CMDs). The latter are essentially point clouds in N-dimensional space, with the number of dimensions corresponding to the photometric bands employed. In this context, machine learning techniques suitable for tabular data are not immediately applicable to CMDs because the number of stars included in a given CMD is variable, and equivariance for permutations is required. To address this issue without introducing ad-hoc manipulations that would require human oversight, here we present a new CMD featurization procedure that summarizes a CMD by means of a quadtree-like structure through iterative partitions of the color-magnitude plane, extracting a fixed number of meaningful features of the relevant subregion from any given CMD. The present approach is robust to photometric noise and contamination and it shows that a simple linear regression on our features predicts distance modulus (metallicity) with a scatter of 0.33 dex (0.16 dex) in cross-validation.

O. S. Salafia, M. E. Ravasio, G. Ghirlanda, I. Mandel

30 pages, 23 figures, submitted to A&A. Comments are welcome!

We describe a model of the short gamma-ray burst (SGRB) population under a `quasi-universal jet' scenario in which jets can differ in their on-axis peak prompt emission luminosity $L_c$, but share a universal angular luminosity profile $\ell(\theta_v)=L(\theta_v)/L_c$ as a function of the viewing angle $\theta_v$. The model is fitted, through a Bayesian hierarchical approach inspired by gravitational wave (GW) population analyses, to 3 observed SGRB samples simultaneously: the Fermi/GBM sample of SGRBs with spectral information in the catalogue (367 events); a flux-complete sample of 16 Swift/BAT SGRBs also detected by GBM, with a measured redshift; and a sample of SGRBs with a binary neutron star (BNS) merger counterpart, which only includes GRB~170817A at present. The results favour a narrow jet core with half-opening angle $\theta_c=2.1_{-1.4}^{+2.4}$ deg (90\% credible intervals from our fiducial `full sample' analysis) whose on-axis peak luminosity is distributed as $p(L_c) \propto L_c^{-A}$ with $A=3.2_{-0.4}^{+0.7}$ above a minimum luminosity $L_c^\star = 5_{-2}^{+11}\times 10^{51}$ erg s$^{-1}$. For $\theta_v>\theta_c$, the luminosity scales as a power law $\ell\propto \theta_v^{-\alpha_L}$ with $\alpha_L=4.7_{-1.4}^{+1.2}$, with no evidence for a break. While the model implies an intrinsic `Yonetoku' correlation between $L$ and the peak photon energy $E_p$, its slope is somewhat shallower $E_p\propto L^{0.4\pm 0.2}$ than the apparent one, and the normalization is offset towards larger $E_p$, due to selection effects. The implied local rate density of SGRBs is between about 100 up to several thousands of events per Gpc$^{3}$ yr, in line with the BNS merger rate density inferred from GW observations. Based on the model, we predict 0.2 to 1.3 joint GW+SGRB detections per year by the Advanced GW detector network and Fermi/GBM during the O4 observing run.

W. S. Oh, T. Nordlander, G. S. Da Costa, M. S. Bessell, A. D. Mackey

Accepted for publication in MNRAS

We present results of a search for extremely metal-poor (EMP) stars in the Large Magellanic Cloud, which can provide crucial information about the properties of the first stars as well as on the formation conditions prevalent during the earliest stages of star formation in dwarf galaxies. Our search utilised SkyMapper photometry, together with parallax and proper motion cuts (from Gaia), colour-magnitude cuts (by selecting the red giant branch region) and finally a metallicity-sensitive cut. Low-resolution spectra of a sample of photometric candidates were taken using the ANU 2.3m telescope/WiFeS spectrograph, from which 7 stars with [Fe/H] $\leq$ -2.75 were identified, two of which have [Fe/H] $\leq$ -3. Radial velocities, derived from the CaII triplet lines, closely match the outer rotation curve of the LMC for the majority of the candidates in our sample. Therefore, our targets are robustly members of the LMC based on their 6D phase-space information (coordinates, spectrophotometric distance, proper motions and radial velocities), and they constitute the most metal-poor stars so far discovered in this galaxy.

F. Kirsten (Chalmers), O. Ould-Boukattine (ASTRON), W. Herrmann (Astropeiler Stockert), M. Gawronski (NCU, Torun), J. Hessels (University of Amsterdam, ASTRON), W. Lu (UC Berkeley), M. Snelders (ASTRON), P. Chawla (University of Amsterdam), J. Yang (Chalmers), R. Blaauw (ASTRON), K. Nimmo (MIT), W. Puchalska (NCU, Torun), P. Wolak (NCU, Torun), R. van Ruiten (University of Amsterdam)

25 pages, 9 figures. Submitted, comments welcome

Fast radio bursts (FRBs) are extremely energetic, millisecond-duration radio flashes that reach Earth from extragalactic distances. Broadly speaking, FRBs can be classified as repeating or (apparently) non-repeating. It is still unclear, however, whether the two types share a common physical origin, differing only in their activity rate. Here we report on an unprecedented observing campaign that targeted one hyperactive repeating source, FRB 20201124A, for more than $2000~\mathrm{hr}$ using four $25-32\mathrm{-m}$ class radio telescopes. In total, we detect $46$ high-energy bursts, many more than one would expect given previous observations of lower-energy bursts using larger radio telescopes. We find a high-energy burst distribution that resembles that of the non-repeating FRB population, suggesting that apparently non-repeating FRB sources may simply be the rarest bursts from repeating sources. We also discuss how FRB 20201124A contributes strongly to the all-sky FRB rate and how similar sources would be observable even at very high redshift.

Jin Beniyama, Shigeyuki Sako, Katsuhito Ohtsuka, Tomohiko Sekiguchi, Masateru Ishiguro, Daisuke Kuroda, Seitaro Urakawa, Fumi Yoshida, Asami Takumi, Natsuho Maeda, Jun Takahashi, Seiko Takagi, Hiroaki Saito, Tatsuya Nakaoka, Tomoki Saito, Tomohiro Ohshima, Ryo Imazawa, Masato Kagitani, Satoshi Takita

Resubmitted to AAS Journals. Any comments are welcome

Asteroid systems such as binaries and pairs are indicative of physical properties and dynamical histories of the Small Solar System Bodies. Although numerous observational and theoretical studies have been carried out, the formation mechanism of asteroid pairs is still unclear, especially for near-Earth asteroid (NEA) pairs. We conducted a series of optical photometric and polarimetric observations of a small NEA 2010 XC$_{15}$ in 2022 December to investigate its surface properties. The rotation period of 2010 XC$_{15}$ is possibly a few to several dozen hours and color indices of 2010 XC$_{15}$ are derived as $g-r=0.435\pm0.008$, $r-i=0.158\pm0.017$, and $r-z=0.186\pm0.009$ in the Pan-STARRS system. The linear polarization degrees of 2010 XC$_{15}$ are a few percent at the phase angle range of 58$^{\circ}$ to 114$^{\circ}$. We found that 2010 XC$_{15}$ is a rare E-type NEA on the basis of its photometric and polarimetric properties. Taking the similarity of not only physical properties but also dynamical integrals and the rarity of E-type NEAs into account, we suppose that 2010 XC$_{15}$ and 1998 WT$_{24}$ are of common origin (i.e., asteroid pair). These two NEAs are the sixth NEA pair and first E-type NEA pair ever confirmed, possibly formed by rotational fission. We conjecture that the parent body of 2010 XC$_{15}$ and 1998 WT$_{24}$ was transported from the main-belt through the $\nu_6$ resonance or Hungaria region.

J. R. Eskilt, K. Lee, D. J. Watts, V. Anshul, R. Aurlien, A. Basyrov, M. Bersanelli, L. P. L. Colombo, H. K. Eriksen, K. S. F. Fornazier, U. Fuskeland, M. Galloway, E. Gjerløw, L. T. Hergt, H. T. Ihle, J. G. S. Lunde, A. Marins, S. K. Nerval, S. Paradiso, F. Rahman, M. San, N.-O. Stutzer, I. K. Wehus

10 pages, 8 figures. Submitted to A&A

We implement support for a cosmological parameter estimation algorithm as proposed by Racine et al. (2016) in Commander, and quantify its computational efficiency and cost. For a semi-realistic simulation similar to Planck LFI 70 GHz, we find that the computational cost of producing one single sample is about 60 CPU-hours and that the typical Markov chain correlation length is $\sim$100 samples. The net effective cost per independent sample is $\sim$6 000 CPU-hours, in comparison with all low-level processing costs of 812 CPU-hours for Planck LFI and WMAP in Cosmoglobe Data Release 1. Thus, although technically possible to run already in its current state, future work should aim to reduce the effective cost per independent sample by at least one order of magnitude to avoid excessive runtimes, for instance through multi-grid preconditioners and/or derivative-based Markov chain sampling schemes. This work demonstrates the computational feasibility of true Bayesian cosmological parameter estimation with end-to-end error propagation for high-precision CMB experiments without likelihood approximations, but it also highlights the need for additional optimizations before it is ready for full production-level analysis.

Yue Shen, J. Andrew Casey-Clyde, Yu-Ching Chen, Arran Gross, Melanie Habouzit, Hsiang-Chih Hwang, Yuzo Ishikawa, Jun-Yao Li, Xin Liu, Chiara M. F. Mingarelli, D. Porquet, Aaron Stemo, Ming-Yang Zhuang

Roman Core Community Survey White Paper, focusing on the High Latitude Wide Area Survey

The hierarchical structure formation paradigm predicts the formation of pairs of supermassive black holes in merging galaxies. When both (or one) members of the SMBH pair are unobscured AGNs, the system can be identified as a dual (or offset) AGN. Quantifying the abundance of these AGN pairs as functions of separation, redshift and host properties is crucial to understanding SMBH formation and AGN fueling in the broad context of galaxy formation. The High Latitude Wide Area Survey with Roman, with its unprecedented combination of sensitivity, spatial resolution, area and NIR wavelength coverage, will revolutionize the study of galactic-scale environments of SMBH pairs. This white paper summarizes the science opportunities and technical requirements on the discovery and characterization of SMBH pairs down to galactic scales (i.e., less than tens of kpc) over broad ranges of redshift (1<z<7) and luminosity (Lbol>1E42 erg/s).

R. Fedriani, A. Caratti o Garatti, R. Cesaroni, J. C. Tan, B. Stecklum, L. Moscadelli, M. Koutoulaki, G. Cosentino, M. Whittle

Accepted for publication in Astronomy and Astrophysics. 10 pages, 8 figures

Massive stars have an impact on their surroundings from early in their formation until the end of their lives. However, very little is known about their formation. Episodic accretion may play a crucial role, but observations of these events have only been reported towards a handful of massive protostars. We aim to investigate the outburst event from the high-mass star-forming region S255IR where recently the protostar NIRS3 underwent an accretion outburst. We follow the evolution of this source both in photometry and morphology of its surroundings. Methods: We perform near-infrared adaptive optics observations on the S255IR central region using the Large Binocular Telescope in the K$_{\rm s}$ broad-band and the H$_2$ and Br$\gamma$ narrow-band filters with an angular resolution of $\sim0\farcs06$, close to the diffraction limit. We discover a new near-infrared knot north-east from NIRS3 that we interpret as a jet knot that was ejected during the last accretion outburst and observed in the radio regime as part of a follow-up after the outburst. We measure a mean tangential velocity for this knot of $450\pm50\,\mathrm{km\,s^{-1}}$. We analyse the continuum-subtracted images from H$_2$ which traces jet shocked emission, and Br$\gamma$ which traces scattered light from a combination of accretion activity and UV radiation from the central massive protostar. We observe a significant decrease in flux at the location of NIRS3, with K=13.48\,mag being the absolute minimum in the historic series. Our observations strongly suggest a scenario where the episodic accretion is followed by an episodic ejection response in the near-infrared, as it was seen in the earlier radio follow-up. The 30 years of $\sim2\,\mu{\rm m}$ photometry suggests that NIRS3 might have undergone another outburst in the late 1980s, being the first massive protostar with such evidence observed in the near-infrared.

Barbara Perri, Brigitte Schmieder, Pascal Démoulin, Stefaan Poedts, Florian Regnault

25 pages, 16 figures, accepted in ApJ

The propagation of Interplanetary Coronal Mass Ejections (ICMEs) in the heliosphere is influenced by many physical phenomena, related to the internal structure of the ICME and its interaction with the ambient solar wind and magnetic field. As the solar magnetic field is modulated by the 11-year dynamo cycle, our goal is to perform a theoretical exploratory study to assess the difference of propagation of an ICME in typical minimum and maximum activity backgrounds. We define a median representative CME at 0.1~au, using both observations and numerical simulations, and describe it using a spheromak model. We use the heliospheric propagator European Heliospheric FORecasting Information Asset (EUHFORIA) to inject the same ICME in two different background wind environments. We then study how the environment and the internal CME structure impact the propagation of the ICME towards Earth, by comparison with an unmagnetized CME. At minimum of activity, the structure of the heliosphere around the ecliptic causes the ICME to slow down, creating a delay with the polar parts of the ejecta. This delay is more important if the ICME is faster. At maximum of activity, a southern coronal hole causes a northward deflection. For these cases, we always find that the ICME at maximum of activity arrives first, while the ICME at minimum of activity is actually more geo-effective. The helicity sign of the ICME is also a crucial parameter but at minimum of activity only, since it affects the magnetic profile and the arrival time of up to 8 hours.

Nagaraj Vernekar, Annapurni Subramaniam, Vikrant V. Jadhav, Dominic M. Bowman

15 pages, 10 figures. Accepted for publication in MNRAS

Blue straggler stars (BSSs) are formed through mass transfer or mergers in binaries. The recent detections of white dwarf (WD) companions to BSSs in M67 suggested a mass transfer pathway of formation. In search of a close companion to five BSSs in M67 that are known to be spectroscopic binaries, we study the light curves from K2 and TESS data. We use PHOEBE to analyse the light curves and estimate the properties of the companions. We detect variability in WOCS 1007, and the light curve is dominated by ellipsoidal variation. Using the light curve and radial velocity measurements, we estimate its orbital period to be 4.212$\pm$0.041 d and $e$ = 0.206$\pm$002. The mass of the companion is estimated to be 0.22$\pm$0.05 M$_{\odot}$ with a radius of 0.078$\pm$0.027 R$_{\odot}$, confirming it to be a low mass WD with T$_{\rm eff}$ = 14300$\pm$1100 K. The estimated mass of the BSS, 1.95$\pm$0.26 M$_{\odot}$, is similar to that estimated from isochrones. The BSS in WOCS 1007 shows $\delta$ Scuti pulsations, although it is slightly deformed and likely to be formed through an efficient mass transfer. Though we detect a light curve for WOCS 4003 showing grazing eclipse with ellipsoidal variation, the estimated parameters are inconclusive. Apart from the 0.44 d period, we found smaller eclipses with a period of 1.1 d, suggesting a compact triple system. In the case of WOCS 4003, WOCS 5005, and WOCS 1025, no eclipses or pulsations are detected, confirming the absence of any short-period inner binary with high inclination in these BSSs.

Chengchao Yuan, Walter Winter

17 pages, 6 figures, 1 table. Submitted to ApJ. Comments are welcome

The potential association between Tidal Disruption Events (TDEs) and high-energy astrophysical neutrinos implies the acceleration of cosmic rays. These accelerated particles will initiate electromagnetic (EM) cascades spanning from keV to GeV energies by the processes related to neutrino production. We model the EM cascade and neutrino emissions by numerically solving the time-dependent transport equations and discuss the implications for AT2019dsg and AT2019fdr in the X-ray and $\gamma$-ray bands. We show that the $\gamma$-ray constraints from \emph{Fermi} can constrain the size of the radiation zone and the maximum energy of injected protons, and that the corresponding expected neutrino event numbers in follow-up searches are limited to be less than about 0.1. Depending on the efficiency of $p\gamma$ interactions, the X-ray and $\gamma$-ray signals can be expected closer to the peak of the optical-ultraviolet (OUV) luminosity, or to the time of the neutrino production.

Zhi-Wei Wang, Samuel L. Braunstein

Publication in Nature Astronomy (2023). arXiv admin note: substantial text overlap with arXiv:2109.10241

Dennis Sciama argued that the existence of life depended on many quantities, the fundamental constants, so in a random universe life should be highly unlikely. However, without full knowledge of these constants, his argument implies a universe that would appear to be `intelligently designed.'

Cyrus E. Dreyer, Rouven Essig, Marivi Fernandez-Serra, Aman Singal, Cheng Zhen

15 pages, plus references and appendices, 12 figures

We calculate target-material responses for dark matter--electron scattering at the \textit{ab-initio} all-electron level using atom-centered gaussian basis sets. The all-electron effects enhance the material response at high momentum transfers from dark matter to electrons, $q\gtrsim \mathcal{O}\left({10\ \alpha m_e}\right)$, compared to calculations using conventional plane wave methods, including those used in QEDark; this enhances the expected event rates at energy transfers $E \gtrsim 10$~eV, especially when scattering through heavy mediators. We carefully test a range of systematic uncertainties in the theory calculation, including those arising from the choice of basis set, exchange-correlation functional, number of unit cells in the Bloch sum, $\mathbf{k}$-mesh, and neglect of scatters with very high momentum transfers. We provide state-of-the-art crystal form factors, focusing on silicon and germanium. Our code and results are made publicly available as a new tool, called Quantum Chemistry Dark (``QCDark'').

Will Barker, Sebastian Zell

13 pages, 5 figures, 3 appendices

We construct a theory of gravity in which a propagating massive vector field arises from a quadratic curvature invariant. The Einstein-Cartan formulation and a partial suppression of torsion ensure the absence of ghost and strong-coupling problems, as we prove with nonlinear Lagrangian and Hamiltonian analysis. Augmenting General Relativity with a propagating torsion vector, our theory provides a purely gravitational origin of Einstein-Proca models and constrains their parameter space. As an outlook to phenomenology, we discuss the gravitational production of fermionic dark matter.

Mostafizur Rahman, Shailesh Kumar, Arpan Bhattacharyya

35 pages, 9 figures, 1 table

The discovery of gravitational waves and black holes has started a new era of gravitational wave astronomy that allows us to probe the underpinning features of gravity and astrophysics in extreme environments of the universe. In this article, we investigate one such study with an extreme mass-ratio inspiral system where the primary object is a spherically symmetric static black hole immersed in a dark matter halo governed by the Hernquist density distribution. We consider the eccentric equatorial orbital motion of the steller-mass object orbiting around the primary and compute measurable effects. We examine the behaviour of dark matter mass and halo radius in generated gravitational wave fluxes and the evolution of eccentric orbital parameters -- eccentricity and semi-latus rectum. We further provide an estimate of gravitational wave dephasing and find the seminal role of low-frequency detectors in the observational prospects of such an astrophysical environment.

Nirmal Raj

8 pages revtex4 + references, 3 figures, 1 table

Cosmologies in which dark matter clumps strongly on small scales are unfavorable to terrestrial detectors that are as yet unexposed to the clumps. I show that sub-hectometer clumps could trigger thermonuclear runaways by scattering on nuclei in white dwarf cores (carbon and oxygen) and neutron star oceans (carbon), setting off Type Ia-like supernovae and x-ray superbursts respectively. I consider two scenarios: ``dark clusters" that are essentially microhalos, and ``long-range dark nuggets", essentially macroscopic composites, with long-range Yukawa baryonic interactions that source the energy for igniting explosions. I constrain dark clusters weighing between the Planck mass and asteroid masses, and long-range dark nuggets over a wider mass range spanning forty orders of magnitude. These limits greatly complement searches I had co-proposed in 2109.04582 for scattering interactions of dark clumps in neutron stars, cosmic rays, and pre-historic minerals.

The title theory is formulated. It entails a quantum-coherent variant of the Fermi-Dirac distribution and casts new light on neutrino oscillations. It might enable the incorporation of neutrino mixing into the modeling of core-collapse supernovae and neutron-star mergers.

Separating signals from an additive mixture may be an unnecessarily hard problem when one is only interested in specific properties of a given signal. In this work, we tackle simpler "statistical component separation" problems that focus on recovering a predefined set of statistical descriptors of a target signal from a noisy mixture. Assuming access to samples of the noise process, we investigate a method devised to match the statistics of the solution candidate corrupted by noise samples with those of the observed mixture. We first analyze the behavior of this method using simple examples with analytically tractable calculations. Then, we apply it in an image denoising context employing 1) wavelet-based descriptors, 2) ConvNet-based descriptors on astrophysics and ImageNet data. In the case of 1), we show that our method better recovers the descriptors of the target data than a standard denoising method in most situations. Additionally, despite not constructed for this purpose, it performs surprisingly well in terms of peak signal-to-noise ratio on full signal reconstruction. In comparison, representation 2) appears less suitable for image denoising. Finally, we extend this method by introducing a diffusive stepwise algorithm which gives a new perspective to the initial method and leads to promising results for image denoising under specific circumstances.

Sergei D. Odintsov, Simone D'Onofrio, Tanmoy Paul

"Physics of the Dark Universe" Accepted

The growing cosmological interest of different entropy functions (like the Tsallis entropy, the R\'{e}nyi entropy, the Barrow entropy, the Sharma-Mittal entropy, the Kaniadakis entropy and the Loop Quantum gravity entropy) naturally raises an important question: "Does there exist a generalized entropy that can bring all the known entropies proposed so far within a single umbrella?" In spirit of this, recently a four parameter generalized entropy has been formulated that reduces to different known entropies for suitable limits of the parameters. Based on such four parameter generalized entropy (symbolized by $S_\mathrm{g}$), in the present paper, we examine the universe's evolution during its early phase, particularly from inflation to reheating, in the context of entropic cosmology where the entropic energy density acts as the inflaton. It turns out that the entropic energy successfully drives an early inflationary phase with a graceful exit, and moreover, the theoretical expectations of the observable indices get consistent with the recent Planck data for suitable ranges of the entropic parameters. After the inflation ends, the universe enters to a reheating stage when the entropic energy decays to relativistic particles with a certain decay rate. Actually the presence of the entropic parameters in the $S_\mathrm{g}$ ensures a continuous evolution of the Hubble parameter from a quasi de-Sitter phase during the inflation to a power law phase during the reheating stage dominated by a constant EoS parameter. Consequently we investigate the reheating phenomenology, and scan the entropic parameters from both the inflation and reheating requirements. We further address the possibility of instantaneous reheating in the present context of generalized entropy.

The equations for quintessential $\alpha$-attractor inflation with a single scalar field, radiation and matter in a spatially flat FLRW spacetime are recast into a regular dynamical system on a compact state space. This enables a complete description of the solution space of these models. The inflationary attractor solution is shown to correspond to the unstable center manifold of a de Sitter fixed point, and we describe connections between slow-roll and dynamical systems approximations for this solution, including Pad\'e approximants. We also introduce a new method for systematically obtaining initial data for quintessence evolution by using dynamical systems properties; in particular, this method exploits that there exists a radiation dominated line of fixed points with an unstable quintessence attractor submanifold, which plays a role that is reminiscent of that of the inflationary attractor solution for inflation.

Jonas P. Pereira, Jorge A. Rueda

10 pages, 4 figures. Accepted for publication in Universe (Feature Papers 2023-Gravitation)

We investigate within the Darmois-Israel thin shell formalism the match of neutral and asymptotically flat, slowly rotating spacetimes (up to the second order in the rotation parameter) when their boundaries are dynamic. It has several important applications in general relativistic systems, such as black holes and neutron stars, which we exemplify. We mostly focus on stability aspects of slowly rotating thin shells in equilibrium and surface degrees of freedom on the hypersurfaces splitting the matched slowly rotating spacetimes, e.g., surface energy density and surface tension. We show that the stability upon perturbations in the spherically symmetric case automatically implies stability in the slow rotation case. In addition, we show that when matching slowly rotating Kerr spacetimes through thin shells in equilibrium, surface degrees of freedom can decrease compared to their Schwarzschild counterparts, meaning that energy conditions could be weakened. Frame-dragging aspects of the match of slowly rotating spacetimes are also briefly discussed.

The small field inflation (SFI) of Coleman-Weinberg (CW) type suffers from precise tuning of the initial inflaton field value to be away from the true vacuum one. We propose a dynamical trapping mechanism to solve this problem: an ultra-supercooling caused by an almost scale-invariant CW potential traps the inflaton at the false vacuum, far away from the true vacuum dominantly created by the quantum scale anomaly, and allows the inflaton to dynamically start the slow-roll down due to a classical explicit-scale breaking effect. To be concrete, we employ a successful CW-SFI model and show that the proposed mechanism works consistently with the observed bounds on the inflation parameters. The proposed new mechanism thus provides new insights for developing small field inflation models.

Iason Baldes, Maximilian Dichtl, Yann Gouttenoire, Filippo Sala

5 pages plus references, 5 figures

In cosmological first-order phase transitions (PT) with relativistic bubble walls, high-energy shells of particles generically form on the inner and outer sides of the walls. Shells from different bubbles can then collide with energies much larger than the PT or inflation scales, and with sizeable rates, realising a `bubbletron'. As an application, we calculate the maximal dark matter mass $M_{DM}$ that can be produced from shell collisions in a U(1) gauge PT, for scales of the PT $v_\varphi$ from MeV to $10^{16}$ GeV. We find for example $M_{DM} \sim 10^6/10^{11}/10^{15}$ GeV for $v_\varphi \sim 10^{-2}/10^3/10^8$ GeV. The gravity wave signal sourced at the PT then links Pulsar Timing Arrays with the PeV scale, LISA with the ZeV one, and the Einstein Telescope with grand unification.

Angela D. V. Di Virgilio (on behalf of GINGER Collaboration)

6 pages, 5 figures

Large frame Ring laser gyroscopes, based on the Sagnac effect, are top sensitivity instrumentation to measure angular velocity with respect to the fixed stars. GINGER (Gyroscopes IN GEneral Relativity) project foresees the construction of an array of three large dimension ring laser gyroscopes, rigidly connected to the Earth. GINGER has the potentiality to measure general relativity effects and Lorentz Violation in the gravity sector, once a sensitivity of $10^{-9}$, or better, of the Earth rotation rate is obtained. Being attached to the Earth crust, the array will also provide useful data for geophysical investigation. For this purpose, it is at present under construction as part of the multi-components observatory called Underground Geophysics at Gran Sasso (UGSS). Sensitivity is the key point to determine the relevance of this instrument for fundamental science. The most recent progress in the sensitivity measurement, obtained on a ring laser prototype called GINGERINO, indicates that GINGER should reach the level of 1 part in $10^{11}$ of the Earth rotation rate.

We study the impact of out-of-equilibrium, dissipative effects on the dynamics of inspiraling neutron stars. We find that modeling dissipative processes (such as those from the stars internal effective fluid viscosity) requires that one introduce a new tidal deformability parameter--the dissipative tidal deformability--which modifies the phase of gravitational waves emitted during the inspiral phase of a neutron star binary. We show that the dissipative tidal deformability corrects the gravitational-wave phase at 4 post-Newtonian order for quasi-circular binaries. This correction receives a large finite-size enhancement by the stellar compactness, analogous to the case of the tidal deformability. Moreover, the correction is not degenerate with the time of coalescence, which also enters at 4PN order, because it contains a logarithmic frequency-dependent contribution. Using a simple Fisher analysis, we show that physically allowed values for the dissipative tidal deformability may be constrained by measurements of the phase of emitted gravitational waves to roughly the same extent as the (electric-type, quadrupolar) tidal deformability. Finally, we show that there are no out-of-equilibrium, dissipative corrections to the tidal deformability itself. We conclude that there are at least two relevant tidal deformability parameters that can be constrained with gravitational-wave phase measurements during the late inspiral of a neutron star binary: one which characterizes the adiabatic tidal response of the star, and another which characterizes the leading-order out-of-equilibrium, dissipative tidal response. These findings open a window to probe dissipative processes in the interior of neutron stars with gravitational waves.

We investigate the radiation of optical angular momentum by a dipole gas under uniform magnetic field with an unpolarized source at its center. Conservation of angular momentum implies that the radiation of angular momentum results in a torque on both the source and the surrounding environment. Moreover, we study the spin and orbital contributions to the radiated angular momentum.