Showing 1–50 of 294 results for author: Rodriguez, C

Once a Triple, Not Always a Triple: The Evolution of Hierarchical Triples that Yield Merged Inner Binaries

Authors: Cheyanne Shariat, Smadar Naoz, Kareem El-Badry, Antonio C. Rodriguez, Bradley M. S. Hansen, Isabel Angelo, Alexander P. Stephan

Abstract: More than half of all main-sequence (MS) stars have one or more companions, and many of those with initial masses <8 M$_\odot$ are born in hierarchical triples. These systems feature two stars in a close orbit (the inner binary) while a tertiary star orbits them on a wider orbit (the outer binary). In hierarchical triples, three-body dynamics combined with stellar evolution drives interactions and… ▽ More More than half of all main-sequence (MS) stars have one or more companions, and many of those with initial masses <8 M$_\odot$ are born in hierarchical triples. These systems feature two stars in a close orbit (the inner binary) while a tertiary star orbits them on a wider orbit (the outer binary). In hierarchical triples, three-body dynamics combined with stellar evolution drives interactions and, in many cases, merges the inner binary entirely to create a renovated `Post-Merger Binary' (PMB). By leveraging dynamical simulations and tracking binary interactions, we explore the outcomes of merged triples and investigate whether PMBs preserve signatures of their three-body history. Our findings indicate that in 26-54% of wide double WD binaries (s>100 au), the more massive white dwarf (WD) is a merger product, implying that these DWD binaries were previously triples. Overall, we estimate that $44\pm14\%$ of observed wide DWDs originated in triple star systems and thereby have rich dynamical histories. Additionally, our results suggest that the separations of inner and outer binaries are uncorrelated at birth, providing insights into stellar formation processes. We also examine MS+MS and MS+Red Giant mergers manifesting as Blue Straggler stars (BSSs). These PMBs have orbital configurations and ages similar to most observed BSS binaries. While the triple+merger formation channel can explain the observed chemical abundances, moderate eccentricities, and companion masses in BSS binaries, it likely only accounts for $\sim$20-25% of BSSs. Meanwhile, we predict that the majority of observed single BSSs formed as collisions in triples and harbor long-period (>10 yr) companions. Furthermore, both BSS binaries and DWDs exhibit signatures of WD birth kicks. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ. Comments are welcome

Isotropy of cosmic rays beyond $10^{20}$ eV favors their heavy mass composition

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the resul… ▽ More We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the results are consistent with a relatively heavy injected composition at E ~ 10 EeV that becomes lighter up to E ~ 100 EeV, while the composition at E > 100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extra-galactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 8 pages, 3 figures, accepted for publication in PRL

Mass composition of ultra-high energy cosmic rays from distribution of their arrival directions with the Telescope Array

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale struc… ▽ More We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale structure (LSS) of the Universe. As we report in the companion letter, the TA data show large deflections with respect to the LSS which can be explained, assuming small extra-galactic magnetic fields (EGMF), by an intermediate composition changing to a heavy one (iron) in the highest energy bin. Here we show that these results are robust to uncertainties in UHECR injection spectra, the energy scale of the experiment and galactic magnetic fields (GMF). The assumption of weak EGMF, however, strongly affects this interpretation at all but the highest energies E > 100 EeV, where the remarkable isotropy of the data implies a heavy injected composition even in the case of strong EGMF. This result also holds if UHECR sources are as rare as $2 \times 10^{-5}$ Mpc$^{-3}$, that is the conservative lower limit for the source number density. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 18 pages, 11 figures, accepted for publication in PRD

Jacobian-Free Newton-Krylov method for multilevel NLTE radiative transfer problems

Authors: D. Arramy, J. de la Cruz Rodríguez, J. Leenaarts

Abstract: The calculation of the emerging radiation from a model atmosphere requires knowledge of the emissivity and absorption coefficients, which are proportional to the atomic level population densities of the levels involved in each transition. Due to the intricate interdependency of the radiation field and the physical state of the atoms, iterative methods are required in order to calculate the atomic… ▽ More The calculation of the emerging radiation from a model atmosphere requires knowledge of the emissivity and absorption coefficients, which are proportional to the atomic level population densities of the levels involved in each transition. Due to the intricate interdependency of the radiation field and the physical state of the atoms, iterative methods are required in order to calculate the atomic level population densities. A variety of different methods have been proposed to solve this problem, which is known as the Non-Local Thermodynamical Equilibrium (NLTE) problem. In this study we have developed a Jacobian-Free Newton-Krylov method (JFNK) to solve multi-level NLTE radiative transfer problems. Using the Rybicki & Hummer (1992) method as a reference (Rybicki, G. B. & Hummer, D. G. 1992, A&A, 262, 209), our results show that our JFNK solver can achieve up to a factor two speed up when using local approximate operators / preconditioners, while also achieving a lower residual error in the statistical equilibrium equations. Another advantage of this method is that the addition of charge conservation and partial redistribution effects should be straight forward. Our method can help accelerating the calculation of the emerging spectra from numerical models and also the reconstruction of chromospheric datasets through NLTE inversions. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: A&A revised version

From Young Massive Clusters to Old Globular Clusters: Density Profile Evolution and IMBH Formation

Authors: Kuldeep Sharma, Carl L. Rodriguez

Abstract: The surface brightness profiles of globular clusters are conventionally described with the well-known King profile. However, observations of young massive clusters (YMCs) in the local Universe suggest that they are better fit by simple models with flat central cores and simple power-law densities in their outer regions (such as the Elson-Fall-Freeman, or EFF, profile). Depending on their initial c… ▽ More The surface brightness profiles of globular clusters are conventionally described with the well-known King profile. However, observations of young massive clusters (YMCs) in the local Universe suggest that they are better fit by simple models with flat central cores and simple power-law densities in their outer regions (such as the Elson-Fall-Freeman, or EFF, profile). Depending on their initial central density, these YMCs may also facilitate large numbers of stellar collisions, potentially creating very massive stars that will directly collapse to intermediate-mass black holes (IMBHs). Using Monte Carlo $N$-body models of YMCs, we show that EFF-profile clusters transform to Wilson or King profiles through natural dynamical evolution, but that their final $W_0$ parameters do not strongly correlate to their initial concentrations. The most centrally-dense YMCs can produce runaway stellar mergers as massive as $4000\,M_{\odot}$ (the largest resolved mass in our simulations) which can collapse to produce IMBHs of similar masses. In doing so, these runaway collisions also deplete the clusters of their primordial massive stars, reducing the number of stellar-mass BHs by as much as $\sim$ 40\%. This depletion will accelerate the core collapse of clusters, suggesting that the process of IMBH formation itself may produce the high densities observed in some core-collapsed clusters. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 22 pages, 9 figures, 3 tables

One-dimensional, geometrically stratified semi-empirical models of the quiet-Sun photosphere and lower chromosphere

Authors: J. M. Borrero, I. Milic, A. Pastor Yabar, A. J. Kaithakkal, J. de la Cruz Rodriguez

Abstract: One-dimensional, semi-empirical models of the solar atmosphere are widely employed in numerous contexts within solar physics, ranging from the determination of element abundances and atomic parameters to studies of the solar irradiance and from Stokes inversions to coronal extrapolations. These models provide the physical parameters (i.e. temperature, gas pressure, etc.) in the solar atmosphere as… ▽ More One-dimensional, semi-empirical models of the solar atmosphere are widely employed in numerous contexts within solar physics, ranging from the determination of element abundances and atomic parameters to studies of the solar irradiance and from Stokes inversions to coronal extrapolations. These models provide the physical parameters (i.e. temperature, gas pressure, etc.) in the solar atmosphere as a function of the continuum optical depth $τ_{\rm c}$. The transformation to the geometrical $z$ scale (i.e. vertical coordinate) is provided via vertical hydrostatic equilibrium. Our aim is to provide updated, one-dimensional, semi-empirical models of the solar atmosphere as a function of $z,$ but employing the more general case of three-dimensional magneto-hydrostatic equilibrium (MHS) instead of vertical hydrostatic equilibrium (HE). We employed a recently developed Stokes inversion code that, along with non-local thermodynamic equilibrium effects, considers MHS instead of HE. This code is applied to spatially and temporally resolved spectropolarimetric observations of the quiet Sun obtained with the CRISP instrument attached to the Swedish Solar Telescope. We provide average models for granules, intergranules, dark magnetic elements, and overall quiet-Sun as a function of both $τ_{\rm c}$ and $z$ from the photosphere to the lower chromosphere. We demonstrate that, in these quiet-Sun models, the effect of considering MHS instead of HE is negligible. However, employing MHS increases the consistency of the inversion results before averaging. We surmise that in regions with stronger magnetic fields (i.e. pores, sunspots, network) the benefits of employing the magneto-hydrostatic approximation will be much more palpable. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: accepted for publication in Astronomy and Astrophysics; 17 pages; 10 figures

A Survey for Radio Emission from White Dwarfs in the VLA Sky Survey

Authors: Ingrid Pelisoli, Laura Chomiuk, Jay Strader, T. R. Marsh, Elias Aydi, Kristen C. Dage, Rebecca Kyer, Isabella Molina, Teresa Panurach, Ryan Urquhart, Thomas J. Maccarone, R. Michael Rich, Antonio C. Rodriguez, E. Breedt, A. J. Brown, V. S. Dhillon, M. J. Dyer, Boris. T. Gaensicke, J. A. Garbutt, M. J. Green, M. R. Kennedy, P. Kerry, S. P. Littlefair, James Munday, S. G. Parsons

Abstract: Radio emission has been detected from tens of white dwarfs, in particular in accreting systems. Additionally, radio emission has been predicted as a possible outcome of a planetary system around a white dwarf. We searched for 3 GHz radio continuum emission in 846,000 candidate white dwarfs previously identified in Gaia using the Very Large Array Sky Survey (VLASS) Epoch 1 Quick Look Catalogue. We… ▽ More Radio emission has been detected from tens of white dwarfs, in particular in accreting systems. Additionally, radio emission has been predicted as a possible outcome of a planetary system around a white dwarf. We searched for 3 GHz radio continuum emission in 846,000 candidate white dwarfs previously identified in Gaia using the Very Large Array Sky Survey (VLASS) Epoch 1 Quick Look Catalogue. We identified 13 candidate white dwarfs with a counterpart in VLASS within 2". Five of those were found not to be white dwarfs in follow-up or archival spectroscopy, whereas seven others were found to be chance alignments with a background source in higher-resolution optical or radio images. The remaining source, WDJ204259.71+152108.06, is found to be a white dwarf and M-dwarf binary with an orbital period of 4.1 days and long-term stochastic optical variability, as well as luminous radio and X-ray emission. For this binary, we find no direct evidence of a background contaminant, and a chance alignment probability of only ~2 per cent. However, other evidence points to the possibility of an unfortunate chance alignment with a background radio and X-ray emitting quasar, including an unusually poor Gaia DR3 astrometric solution for this source. With at most one possible radio emitting white dwarf found, we conclude that strong (> 1-3 mJy) radio emission from white dwarfs in the 3 GHz band is virtually nonexistent outside of interacting binaries. △ Less

Submitted 14 May, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: 18 pages, 17 figures. Updated to match version accepted for publication in MNRAS

Improved reconstruction of solar magnetic fields from imaging spectropolarimetry through spatio-temporal regularisation

Authors: Jaime de la Cruz Rodríguez, Jorrit Leenaarts

Abstract: Determination of solar magnetic fields with a spatial resolution set by the diffraction limit of a telescope is difficult because the time required to measure the Stokes vector with sufficient signal-to-noise is long compared to the solar evolution timescale. This difficulty gets worse with increasing telescope size as the photon flux per diffraction-limited resolution element remains constant but… ▽ More Determination of solar magnetic fields with a spatial resolution set by the diffraction limit of a telescope is difficult because the time required to measure the Stokes vector with sufficient signal-to-noise is long compared to the solar evolution timescale. This difficulty gets worse with increasing telescope size as the photon flux per diffraction-limited resolution element remains constant but the evolution timescale decreases linearly with the diffraction-limited resolution. We aim to improve magnetic field reconstruction at the diffraction limit without averaging the observations in time or space, and without applying noise filtering. The magnetic field vector tends to evolve slower than the temperature, velocity and microturbulence. We exploit this by adding temporal regularisation terms for the magnetic field to the linear least-squares fitting used in the weak-field approximation, as well as to the Levenberg-Marquardt algorithm used in inversions. The other model parameters are allowed to change in time without constraints. We infer the chromospheric magnetic field from Ca II 854.2 nm observations using the weak field approximation and the photospheric magnetic field from Fe I 617.3 nm observations using Milne-Eddington inversions, both with and without temporal regularisation. Temporal regularisation reduce the noise in the reconstructed maps of the magnetic field and provides a better coherency in time in both the weak-field approximation and Milne-Eddington inversions. Temporal regularisation markedly improves magnetic field determination from spatially and temporally resolved observations. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 8 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 685, A85 (2024)

Dramatic rebrightening of the type-changing stripped-envelope supernova SN 2023aew

Authors: Yashvi Sharma, Jesper Sollerman, Shrinivas R. Kulkarni, Takashi J. Moriya, Steve Schulze, Stan Barmentloo, Michael Fausnaugh, Avishay Gal-Yam, Anders Jerkstrand, Tomás Ahumada, Eric C. Bellm, Kaustav K. Das, Andrew Drake, Christoffer Fremling, Saarah Hall, K. R. Hinds, Theophile Jegou du Laz, Viraj Karambelkar, Mansi M. Kasliwal, Frank J. Masci, Adam A. Miller, Guy Nir, Daniel A. Perley, Josiah N. Purdum, Yu-Jing Qin , et al. (10 additional authors not shown)

Abstract: Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and lo… ▽ More Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and long (~100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew's spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak's photospheric phase, before H$α$ resurfaces again during the nebular phase. The nebular lines ([O I], [Ca II], Mg I], H$α$) exhibit a double-peaked structure which hints towards a clumpy or non-spherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than normal SESNe as well as requiring a very high $^{56}$Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak which also powers the second peak of the light curve through SN-CSM interaction. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 22 pages, 11 figures, 5 tables

Spatial resolution effects on the solar open flux estimates

Authors: Ivan Milic, Rebecca Centeno, Xudong Sun, Matthias Rempel, Jaime de la Cruz Rodriguez

Abstract: Spectropolarimetric observations used to infer the solar magnetic fields are obtained with a limited spatial resolution. The effects of this limited resolution on the inference of the open flux over the observed region have not been extensively studied. We aim to characterize the biases that arise in the inference of the mean flux density by performing an end-to-end study that involves the generat… ▽ More Spectropolarimetric observations used to infer the solar magnetic fields are obtained with a limited spatial resolution. The effects of this limited resolution on the inference of the open flux over the observed region have not been extensively studied. We aim to characterize the biases that arise in the inference of the mean flux density by performing an end-to-end study that involves the generation of synthetic data, its interpretation (inversion), and a comparison of the results with the original model. We synthesized polarized spectra of the two magnetically sensitive lines of neutral iron around 630\,nm from a state-of-the-art numerical simulation of the solar photosphere. We then performed data degradation to simulate the effect of the telescope with a limited angular resolution and interpreted (inverted) the data using a Milne-Eddington spectropolarimetric inversion code. We then studied the dependence of the inferred parameters on the telescope resolution. The results show a significant decrease in the mean magnetic flux density -- related to the open flux observed at the disk center -- with decreasing telescope resolution. The original net magnetic field flux is fully resolved by a 1m telescope, but a 20\,cm aperture telescope yields a 30\% smaller value. Even in the fully resolved case, the result is still biased due to the corrugation of the photospheric surface. Even the spatially averaged quantities, such as the open magnetic flux in the observed region, are underestimated when the magnetic structures are unresolved. The reason for this is the presence of nonlinearities in the magnetic field inference process. This effect might have implications for the modeling of large-scale solar magnetic fields; for example, those corresponding to the coronal holes, or the polar magnetic fields, which are relevant to our understanding of the solar cycle. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: To be published in A&A

Magnetic braking below the cataclysmic variable period gap and the observed dearth of period bouncers

Authors: Arnab Sarkar, Antonio C. Rodriguez, Sivan Ginzburg, Lev Yungelson, Christopher A. Tout

Abstract: Period bouncers are cataclysmic variables (CVs) that have evolved past their orbital period minimum. The strong disagreement between theory and observations of the relative fraction of period bouncers is a severe shortcoming in the understanding of CV evolution. We test the implications of the hypothesis that magnetic braking (MB), which is suggested to be an additional angular momentum loss (AML)… ▽ More Period bouncers are cataclysmic variables (CVs) that have evolved past their orbital period minimum. The strong disagreement between theory and observations of the relative fraction of period bouncers is a severe shortcoming in the understanding of CV evolution. We test the implications of the hypothesis that magnetic braking (MB), which is suggested to be an additional angular momentum loss (AML) mechanism for CVs below the period gap ($P_\mathrm{orb}\lesssim 120$ min), weakens around their period minimum. We compute the evolution of CV donors below the period gap using the MESA code, assuming that the evolution of the system is driven by AML by gravitational wave radiation (GWR) and MB. We parametrize the MB strength as $\mathrm{AML_{MB}}=κ\mathrm{AML_{GWR}}$. We compute two qualitatively different sets of models, one where $κ$ is a constant and the other where $κ$ depends on stellar parameters. We find that two crucial effects drive the latter set of models. (1) A decrease in $κ$ as CVs approach the period minimum stalls their evolution so that they spend a long time in the observed period minimum spike ($80\lesssim P_\mathrm{orb}/\,\mathrm{min}\lesssim 86$). Here, they become difficult to distinguish from pre-bounce systems in the spike. (2) A strong decrease in the mass-transfer rate makes them virtually undetectable as they evolve further. So, the CV stalls around the period minimum and then `disappears'. This reduces the number of detectable bouncers. Physical processes, such as dynamo action, white dwarf magnetism, and dead zones, may cause such a weakening of MB at short orbital periods. The weakening magnetic braking formalism provides a possible solution to the problem of the lack of period bouncers in CV observational surveys. △ Less

Submitted 19 June, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A Letters. Significant revision from the original version

From Active Stars to Black Holes: A Tool for the SRG/eROSITA X-ray Survey and New Discoveries as Proof of Concept

Authors: Antonio C. Rodriguez

Abstract: Galactic X-ray sources are diverse, ranging from active M dwarfs to compact object binaries, and everything in between. The X-ray landscape of today is rich, with point source catalogs such as those from XMM-Newton, Chandra, and Swift, each with $\gtrsim10^5$ sources and growing. Furthermore, X-ray astronomy is on the verge of being transformed through data releases from the all-sky SRG/eROSITA su… ▽ More Galactic X-ray sources are diverse, ranging from active M dwarfs to compact object binaries, and everything in between. The X-ray landscape of today is rich, with point source catalogs such as those from XMM-Newton, Chandra, and Swift, each with $\gtrsim10^5$ sources and growing. Furthermore, X-ray astronomy is on the verge of being transformed through data releases from the all-sky SRG/eROSITA survey. Many X-ray sources can be associated with an optical counterpart, which in the era of Gaia, can be determined to be Galactic or extragalactic through parallax and proper motion information. Here, I present a simple diagram -- the ``X-ray Main Sequence", which distinguishes between compact objects and active stars based on their optical color and X-ray-to-optical flux ratio ($F_X/F_\textrm{opt}$). As a proof of concept, I present optical spectroscopy of six exotic accreting WDs discovered using the X-ray Main Sequence as applied to the XMM-Newton catalog. Looking ahead to surveys of the near future, I additionally present SDSS-V optical spectroscopy of new systems discovered using the X-ray Main Sequence as applied to the SRG/eROSITA eFEDS catalog. △ Less

Submitted 17 January, 2024; originally announced January 2024.

Comments: Submitted to PASP

A Joint SRG/eROSITA + ZTF Search: Discovery of a 97-min Period Eclipsing Cataclysmic Variable with Evidence of a Brown Dwarf Secondary

Authors: Ilkham Galiullin, Antonio C. Rodriguez, Shrinivas R. Kulkarni, Rashid Sunyaev, Marat Gilfanov, Ilfan Bikmaev, Lev Yungelson, Jan van Roestel, Boris T. Gänsicke, Irek Khamitov, Paula Szkody, Kareem El-Badry, Mikhail Suslikov, Thomas A. Prince, Mikhail Buntov, Ilaria Caiazzo, Mark Gorbachev, Matthew J. Graham, Rustam Gumerov, Eldar Irtuganov, Russ R. Laher, Pavel Medvedev, Reed Riddle, Ben Rusholme, Nail Sakhibullin , et al. (2 additional authors not shown)

Abstract: Cataclysmic variables (CVs) that have evolved past the period minimum during their lifetimes are predicted to be systems with a brown dwarf donor. While population synthesis models predict that around $\approx 40-70\%$ of the Galactic CVs are post-period minimum systems referred to as "period bouncers", only a few dozen confirmed systems are known. We report the study and characterisation of a new… ▽ More Cataclysmic variables (CVs) that have evolved past the period minimum during their lifetimes are predicted to be systems with a brown dwarf donor. While population synthesis models predict that around $\approx 40-70\%$ of the Galactic CVs are post-period minimum systems referred to as "period bouncers", only a few dozen confirmed systems are known. We report the study and characterisation of a new eclipsing CV, SRGeJ041130.3+685350 (SRGeJ0411), discovered from a joint SRG/eROSITA and ZTF program. The optical spectrum of SRGeJ0411 shows prominent hydrogen and helium emission lines, typical for CVs. We obtained optical high-speed photometry to confirm the eclipse of SRGeJ0411 and determine the orbital period to be $P_\textrm{orb} \approx 97.530$ minutes. The spectral energy distribution suggests that the donor has an effective temperature of $\lesssim 1,800$ K. We constrain the donor mass with the period--density relationship for Roche-lobe-filling stars and find that $M_\textrm{donor} \lesssim 0.04\ M_\odot$. The binary parameters are consistent with evolutionary models for post-period minimum CVs, suggesting that SRGeJ0411 is a new period bouncer. The optical emission lines of SRGeJ0411 are single-peaked despite the system being eclipsing, which is typically only seen due to stream-fed accretion in polars. X-ray spectroscopy hints that the white dwarf in SRGeJ0411 could be magnetic, but verifying the magnetic nature of SRGeJ0411 requires further investigation. The lack of optical outbursts has made SRGeJ0411 elusive in previous surveys, and joint X-ray and optical surveys highlight the potential for discovering similar systems in the near future. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 17 pages, 16 figures, accepted to MNRAS

Microlensing Events in Five Years of Photometry from the Zwicky Transient Facility

Authors: Ruocheng Zhai, Antonio C. Rodriguez, Casey Y. Lam, Eric C. Bellm, Josiah Purdum, Frank J. Masci, Avery Wold

Abstract: Microlensing has a unique advantage for detecting dark objects in the Milky Way, such as free floating planets, neutron stars, and stellar-mass black holes. Most microlensing surveys focus towards the Galactic bulge, where higher stellar density leads to a higher event rate. However, microlensing events in the Galactic plane are closer, and take place over longer timescales. This enables a better… ▽ More Microlensing has a unique advantage for detecting dark objects in the Milky Way, such as free floating planets, neutron stars, and stellar-mass black holes. Most microlensing surveys focus towards the Galactic bulge, where higher stellar density leads to a higher event rate. However, microlensing events in the Galactic plane are closer, and take place over longer timescales. This enables a better measurement of the microlensing parallax, which serves as an independent constraint on the mass of the dark lens. In this work, we systematically searched for microlensing events in Zwicky Transient Facility (ZTF) Data Release 17 from 2018--2023 in the Galactic plane region $|b| < 20^\circ$. We find 124 high-confidence microlensing events and 54 possible events. In the event selection, we use the efficient \texttt{EventFinder} algorithm to detect microlensing signals, which could be used for large datasets such as future ZTF data releases or data from the Rubin Observatory Legacy Survey of Space and Time (LSST). With detection efficiencies of ZTF fields from catalog-level simulations, we calculate the mean Einstein timescale to be $\langle t_\mathrm{E}\rangle = 51.7 \pm 3.3$ days, smaller than previous results of the Galactic plane to within 1.5-$σ$. We calculate optical depths and event rates, which we interpret with caution due to the use of visual inspection in creating our final sample. With two years of additional ZTF data in DR17, we have more than doubled the amount of microlensing events (60) found in the three-year DR5 search and found events with longer Einstein timescales than before. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Comments: 9 figures, 3 tables. Submitted to ApJ

Stellar Black Holes and Compact Stellar Remnants

Authors: Guglielmo Costa, Martyna Chruślińska, Jakub Klencki, Floor S. Broekgaarden, Carl L. Rodriguez, Tana D. Joseph, Sara Saracino

Abstract: The recent observations of gravitational waves (GWs) by the LIGO-Virgo-KAGRA collaboration (LVK) have provided a new opportunity for studying our Universe. By detecting several merging events of black holes (BHs), LVK has spurred the astronomical community to improve theoretical models of single, binary, and multiple star evolution in order to better understand the formation of binary black hole (… ▽ More The recent observations of gravitational waves (GWs) by the LIGO-Virgo-KAGRA collaboration (LVK) have provided a new opportunity for studying our Universe. By detecting several merging events of black holes (BHs), LVK has spurred the astronomical community to improve theoretical models of single, binary, and multiple star evolution in order to better understand the formation of binary black hole (BBH) systems and interpret their observed properties. The final BBH system configuration before the merger depends on several processes, including those related to the evolution of the inner stellar structure and those due to the interaction with the companion and the environment (such as in stellar clusters). This chapter provides a summary of the formation scenarios of stellar BHs in single, binary, and multiple systems. We review all the important physical processes that affect the formation of BHs and discuss the methodologies used to detect these elusive objects and constrain their properties. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: To appear in Chapter 1 in the book Black Holes in the Era of Gravitational Wave Astronomy, ed. Arca Sedda, Bortolas, Spera, pub. Elsevier. All authors equally contributed to the chapter. Figures from other publications have been reproduced with permission

Great Balls of FIRE III: Modeling Black Hole Mergers from Massive Star Clusters in Simulations of Galaxies

Authors: Tristan Bruel, Carl L. Rodriguez, Astrid Lamberts, Michael Y. Grudic, Zachary Hafen, Robert Feldmann

Abstract: After the nearly hundred gravitational-wave detections reported by the LIGO-Virgo-KAGRA Collaboration, the question of the cosmological origin of merging binary black holes (BBHs) remains open. The two main formation channels generally considered are from isolated field binaries or via dynamical assembly in dense star clusters. Here, we focus on understanding the dynamical formation of merging BBH… ▽ More After the nearly hundred gravitational-wave detections reported by the LIGO-Virgo-KAGRA Collaboration, the question of the cosmological origin of merging binary black holes (BBHs) remains open. The two main formation channels generally considered are from isolated field binaries or via dynamical assembly in dense star clusters. Here, we focus on understanding the dynamical formation of merging BBHs within massive clusters in galaxies of different masses. To this end, we apply a new framework to consistently model the formation and evolution of massive star clusters in zoom-in cosmological simulations of galaxies. Each simulation, taken from the FIRE project, provides a realistic star formation environment with a unique star formation history and hosts realistic giant molecular clouds that constitute the birthplace of star clusters. Combined with the code for star cluster evolution CMC, we are able to produce populations of dynamically formed merging BBHs across cosmic time in different environments. As the most massive star clusters preferentially form in dense massive clouds of gas, we find that, despite their low metallicities favourable to the creation of black holes, low-mass galaxies contain few massive clusters and therefore have a limited contribution to the global production of dynamically formed merging BBHs. Furthermore, we find that massive clusters can host hierarchical BBH mergers with clear identifiable physical properties. Looking at the evolution of the BBH merger rate in different galaxies, we find strong correlations between BBH mergers and the most extreme episodes of star formation. Finally, we discuss the implications for future LIGO-Virgo-KAGRA gravitational wave observations. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 14 pages, 9 figures, 3 tables

Journal ref: A&A 686, A106 (2024)

Minutes-duration Optical Flares with Supernova Luminosities

Authors: Anna Y. Q. Ho, Daniel A. Perley, Ping Chen, Steve Schulze, Vik Dhillon, Harsh Kumar, Aswin Suresh, Vishwajeet Swain, Michael Bremer, Stephen J. Smartt, Joseph P. Anderson, G. C. Anupama, Supachai Awiphan, Sudhanshu Barway, Eric C. Bellm, Sagi Ben-Ami, Varun Bhalerao, Thomas de Boer, Thomas G. Brink, Rick Burruss, Poonam Chandra, Ting-Wan Chen, Wen-Ping Chen, Jeff Cooke, Michael W. Coughlin , et al. (52 additional authors not shown)

Abstract: In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Seve… ▽ More In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source, such as X-ray variability, prolonged ultraviolet emission, a tentative X-ray quasiperiodic oscillation, and large energies coupled to fast (but subrelativistic) radio-emitting ejecta. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the "Tasmanian Devil"). The flares occur over a period of months, are highly energetic, and are likely nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that in some AT2018cow-like transients the embedded energy source is a compact object, either a magnetar or an accreting black hole. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 79 pages, 3 figures (main text) + 7 figures (extended data) + 2 figures (supplementary information). Published online in Nature on 15 November 2023

No X-Rays or Radio from the Nearest Black Holes and Implications for Future Searches

Authors: Antonio C. Rodriguez, Yvette Cendes, Kareem El-Badry, Edo Berger

Abstract: Astrometry from the Gaia mission was recently used to discover the two nearest known stellar-mass black holes (BHs), Gaia BH1 and Gaia BH2. Both systems contain $\sim 1\,M_{\odot}$ stars in wide orbits ($a\approx$1.4 AU, 4.96 AU) around $\sim9\,M_{\odot}$ BHs. These objects are among the first stellar-mass BHs not discovered via X-rays or gravitational waves. The companion stars -- a solar-type ma… ▽ More Astrometry from the Gaia mission was recently used to discover the two nearest known stellar-mass black holes (BHs), Gaia BH1 and Gaia BH2. Both systems contain $\sim 1\,M_{\odot}$ stars in wide orbits ($a\approx$1.4 AU, 4.96 AU) around $\sim9\,M_{\odot}$ BHs. These objects are among the first stellar-mass BHs not discovered via X-rays or gravitational waves. The companion stars -- a solar-type main sequence star in Gaia BH1 and a low-luminosity red giant in Gaia BH2 -- are well within their Roche lobes. However, the BHs are still expected to accrete stellar winds, leading to potentially detectable X-ray or radio emission. Here, we report observations of both systems with the Chandra X-ray Observatory and radio observations with the Very Large Array (for Gaia BH1) and MeerKAT (for Gaia BH2). We did not detect either system, leading to X-ray upper limits of $L_X < 10^{29.4}$ and $L_X < 10^{30.1}\,\rm erg\,s^{-1}$ and radio upper limits of $L_r < 10^{25.2}$ and $L_r < 10^{25.9}\,\rm erg\,s^{-1}$. For Gaia BH2, the non-detection implies that the the accretion rate near the horizon is much lower than the Bondi rate, consistent with recent models for hot accretion flows. We discuss implications of these non-detections for broader BH searches, concluding that it is unlikely that isolated BHs will be detected via ISM accretion in the near future. We also calculate evolutionary models for the binaries' future evolution using Modules for Experiments in Stellar Astrophysics (MESA). We find that Gaia BH1 will be X-ray bright for 5--50 Myr when the star is a red giant, including 5 Myr of stable Roche lobe overflow. Since no symbiotic BH X-ray binaries are known, this implies either that fewer than $\sim 10^4$ Gaia BH1-like binaries exist in the Milky Way, or that they are common but have evaded detection, perhaps due to very long outburst recurrence timescales. △ Less

Submitted 15 November, 2023; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: Submitted to PASP

The High Energy X-ray Probe (HEX-P): Probing Accretion onto Stellar Mass Black Holes

Authors: Riley Connors, John Tomsick, Paul Draghis, Benjamin Coughenour, Aarran Shaw, Javier Garcia, Dominic Walton, Kristin Madsen, Daniel Stern, Nicole Cavero Rodriguez, Thomas Dauser, Melania Del Santo, Jiachen Jiang, Henric Krawczynski, Honghui Liu, Joseph Neilsen, Michael Nowak, Sean Pike, Andrea Santangelo, Navin Sridhar, Andrew West, Joern Wilms, the HEX-P Team

Abstract: Accretion is a universal astrophysical process that plays a key role in cosmic history, from the epoch of reionization to galaxy and stellar formation and evolution. Accreting stellar-mass black holes in X-ray binaries are one of the best laboratories to study the accretion process and probe strong gravity -- and most importantly, to measure the angular momentum, or spin, of black holes, and its r… ▽ More Accretion is a universal astrophysical process that plays a key role in cosmic history, from the epoch of reionization to galaxy and stellar formation and evolution. Accreting stellar-mass black holes in X-ray binaries are one of the best laboratories to study the accretion process and probe strong gravity -- and most importantly, to measure the angular momentum, or spin, of black holes, and its role as a powering mechanism for relativistic astrophysical phenomena. Comprehensive characterization of the disk-corona system of accreting black holes, and their co-evolution, is fundamental to measurements of black hole spin. Here, we use simulated data to demonstrate how key unanswered questions in the study of accreting stellar-mass black holes will be addressed by the {\it High Energy X-ray Probe} (\hexp). \hexp\ is a probe-class mission concept that will combine high spatial resolution X-ray imaging and broad spectral coverage ($0.2\mbox{--}80$keV) with a sensitivity superior to current facilities (including \xmm\ and \nustar) to enable revolutionary new insights into a variety of important astrophysical problems. We illustrate the capability of \hexp\ to: 1) measure the evolving structures of black hole binary accretion flows down to low ($\lesssim0.1\%$) Eddington-scaled luminosities via detailed X-ray reflection spectroscopy; 2) provide unprecedented spectral observations of the coronal plasma, probing its elusive geometry and energetics; 3) perform detailed broadband studies of stellar mass black holes in nearby galaxies, thus expanding the repertoire of sources we can use to study accretion physics and determine the fundamental nature of black holes; and 4) act as a complementary observatory to a range of future ground and space-based astronomical observatories, thus providing key spectral measurements of the multi-component emission from the inner accretion flows of BH-XRBs. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 31 pages, 7 figures

Measuring the CMB primordial B-modes with Bolometric Interferometry

Authors: A. Mennella, P. Ade, A. Almela, G. Amico, L. H. Arnaldi, J. Aumont, S. Banfi, E. S. Battistelli, B. Bélier, L. Bergé, J. -Ph. Bernard, P. de Bernardis, M. Bersanelli, J. Bonaparte, J. D. Bonilla, E. Bunn, D. Buzi, F. Cacciotti, D. Camilieri, F. Cavaliere, P. Chanial, C. Chapron, L. Colombo, F. Columbro, A. Coppolecchia , et al. (89 additional authors not shown)

Abstract: The Q&U Bolometric Interferometer for Cosmology (QUBIC) is the first bolometric interferometer designed to measure the primordial B-mode polarization of the Cosmic Microwave Background (CMB). Bolometric interferometry is a novel technique that combines the sensitivity of bolometric detectors with the control of systematic effects that is typical of interferometry, both key features in the quest fo… ▽ More The Q&U Bolometric Interferometer for Cosmology (QUBIC) is the first bolometric interferometer designed to measure the primordial B-mode polarization of the Cosmic Microwave Background (CMB). Bolometric interferometry is a novel technique that combines the sensitivity of bolometric detectors with the control of systematic effects that is typical of interferometry, both key features in the quest for the faint signal of the primordial B-modes. A unique feature is the so-called "spectral imaging", i.e., the ability to recover the sky signal in several sub-bands within the physical band during data analysis. This feature provides an in-band spectral resolution of Δν/ν \sim 0.04 that is unattainable by a traditional imager. This is a key tool for controlling the Galactic foregrounds contamination. In this paper, we describe the principles of bolometric interferometry, the current status of the QUBIC experiment and future prospects. △ Less

Submitted 5 November, 2023; originally announced November 2023.

Comments: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Characterization of diamond-turned optics for SCALES

Authors: Isabel J. Kain, Phil Hinz, Marius Doetz, Benjamin Bulla, Renate Kupke, Daren Dillon, Andrew Skemer, Deno Stelter, Michael Gonzales, Nicholas MacDonald, Aditi Gangadharan, Cristian Rodriguez, Christopher Ratliff, Mackenzie R. Lach, Steph Sallum

Abstract: High-contrast imaging has been used to discover and characterize dozens of exoplanets to date. The primary limiting performance factor for these instruments is contrast, the ratio of exoplanet to host star brightness that an instrument can successfully resolve. Contrast is largely determined by wavefront error, consisting of uncorrected atmospheric turbulence and optical aberrations downstream of… ▽ More High-contrast imaging has been used to discover and characterize dozens of exoplanets to date. The primary limiting performance factor for these instruments is contrast, the ratio of exoplanet to host star brightness that an instrument can successfully resolve. Contrast is largely determined by wavefront error, consisting of uncorrected atmospheric turbulence and optical aberrations downstream of AO correction. Single-point diamond turning allows for high-precision optics to be manufactured for use in astronomical instrumentation, presenting a cheaper and more versatile alternative to conventional glass polishing. This work presents measurements of wavefront error for diamond-turned aluminum optics in the Slicer Combined with an Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument, a 2-5 micron coronagraphic integral field spectrograph under construction for Keck Observatory. Wavefront error measurements for these optics are used to simulate SCALES' point spread function using physical optics propagation software poppy, showing that SCALES' contrast performance is not limited by wavefront error from internal instrument optics. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: Techniques and Instrumentation for Detection of Exoplanets XI

Journal ref: SPIE Proceedings Volume 12680, Techniques and Instrumentation for Detection of Exoplanets XI; 1268020 (2023)

An Extremely Massive White Dwarf Escaped From the Hyades Star Cluster

Authors: David R. Miller, Ilaria Caiazzo, Jeremy Heyl, Harvey B. Richer, Kareem El-Badry, Antonio C. Rodriguez, Zachary P. Vanderbosch, Jan van Roestel

Abstract: We searched the Gaia DR3 database for ultramassive white dwarfs with kinematics consistent with having escaped the nearby Hyades open cluster, identifying three such candidates. Two of these candidates have masses estimated from Gaia photometry of approximately 1.1 solar masses; their status as products of single stellar evolution that have escaped the cluster was deemed too questionable for immed… ▽ More We searched the Gaia DR3 database for ultramassive white dwarfs with kinematics consistent with having escaped the nearby Hyades open cluster, identifying three such candidates. Two of these candidates have masses estimated from Gaia photometry of approximately 1.1 solar masses; their status as products of single stellar evolution that have escaped the cluster was deemed too questionable for immediate follow-up analysis. The remaining candidate has an expected mass >1.3 solar masses, significantly reducing the probability of it being an interloper. Analysis of follow-up Gemini GMOS spectroscopy for this source reveals a non-magnetized hydrogen atmosphere white dwarf with a mass and age consistent with having formed from a single star. Assuming a single-stellar evolution formation channel, we estimate a 97.8% chance that the candidate is a true escapee from the Hyades. With a determined mass of 1.317 solar masses, this is potentially the most massive known single-evolution white dwarf and is by far the most massive with a strong association with an open cluster. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 13 pages, 4 figures, 1 table. Accepted for publication in the Astrophysical Journal Letters

Solar Atmospheric Heating Due to Small-scale Events in an Emerging Flux Region

Authors: Rahul Yadav, Maria D. Kazachenko, Andrey N. Afanasyev, Jaime de la Cruz Rodríguez, Jorrit Leenaarts

Abstract: We investigate the thermal, kinematic and magnetic structure of small-scale heating events in an emerging flux region (EFR). We use high-resolution multi-line observations (including Ca II 8542~Å, Ca II K, and Fe I 6301~Åline pair) of an EFR located close to the disk center from the CRISP and CHROMIS instruments at the Swedish 1-m Solar Telescope. We perform non-LTE inversions of multiple spectral… ▽ More We investigate the thermal, kinematic and magnetic structure of small-scale heating events in an emerging flux region (EFR). We use high-resolution multi-line observations (including Ca II 8542~Å, Ca II K, and Fe I 6301~Åline pair) of an EFR located close to the disk center from the CRISP and CHROMIS instruments at the Swedish 1-m Solar Telescope. We perform non-LTE inversions of multiple spectral lines to infer the temperature, velocity, and magnetic field structure of the heating events. Additionally, we use the data-driven Coronal Global Evolutionary Model to simulate the evolution of the 3D magnetic field configuration above the events and understand their dynamics. Furthermore, we analyze the differential emission measure to gain insights into the heating of the coronal plasma in the EFR. Our analysis reveals the presence of numerous small-scale heating events in the EFR, primarily located at polarity inversion lines of bipolar structures. These events not only heat the lower atmosphere but also significantly heat the corona. The data-driven simulations, along with the observed enhancement of currents and Poynting flux, suggest that magnetic reconnection in the lower atmosphere is likely responsible for the observed heating at these sites. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: 20 pages, 15 figures, accepted for publication in the ApJ

A rotating white dwarf shows different compositions on its opposite faces

Authors: Ilaria Caiazzo, Kevin B. Burdge, Pier-Emmanuel Tremblay, James Fuller, Lilia Ferrario, Boris T. Gaensicke, J. J. Hermes, Jeremy Heyl, Adela Kawka, S. R. Kulkarni, Thomas R. Marsh, Przemek Mroz, Thomas A. Prince, Harvey B. Richer, Antonio C. Rodriguez, Jan van Roestel, Zachary P. Vanderbosch, Stephane Vennes, Dayal Wickramasinghe, Vikram S. Dhillon, Stuart P. Littlefair, James Munday, Ingrid Pelisoli, Daniel Perley, Eric C. Bellm , et al. (13 additional authors not shown)

Abstract: White dwarfs, the extremely dense remnants left behind by most stars after their death, are characterised by a mass comparable to that of the Sun compressed into the size of an Earth-like planet. In the resulting strong gravity, heavy elements sink toward the centre and the upper layer of the atmosphere contains only the lightest element present, usually hydrogen or helium. Several mechanisms comp… ▽ More White dwarfs, the extremely dense remnants left behind by most stars after their death, are characterised by a mass comparable to that of the Sun compressed into the size of an Earth-like planet. In the resulting strong gravity, heavy elements sink toward the centre and the upper layer of the atmosphere contains only the lightest element present, usually hydrogen or helium. Several mechanisms compete with gravitational settling to change a white dwarf's surface composition as it cools, and the fraction of white dwarfs with helium atmospheres is known to increase by a factor ~2.5 below a temperature of about 30,000 K; therefore, some white dwarfs that appear to have hydrogen-dominated atmospheres above 30,000 K are bound to transition to be helium-dominated as they cool below it. Here we report observations of ZTF J203349.8+322901.1, a transitioning white dwarf with two faces: one side of its atmosphere is dominated by hydrogen and the other one by helium. This peculiar nature is likely caused by the presence of a small magnetic field, which creates an inhomogeneity in temperature, pressure or mixing strength over the surface. ZTF J203349.8+322901.1 might be the most extreme member of a class of magnetic, transitioning white dwarfs -- together with GD 323, a white dwarf that shows similar but much more subtle variations. This new class could help shed light on the physical mechanisms behind white dwarf spectral evolution. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Comments: 45 pages, 11 figures

Journal ref: Nature 620, 61-66 (2023)

A transiting brown dwarf in a 2 hour orbit

Authors: Kareem El-Badry, Kevin B. Burdge, Jan van Roestel, Antonio C. Rodriguez

Abstract: We report the discovery of ZTF J2020+5033, a high-mass brown dwarf (BD) transiting a low-mass star with an orbital period of 1.90 hours. Phase-resolved spectroscopy, optical and infrared light curves, and precise astrometry from Gaia allow us to constrain the masses, radii, and temperatures of both components with few-percent precision. We infer a BD mass of $M_{\rm BD} = 80.1\pm 1.6\,M_{\rm J}$,… ▽ More We report the discovery of ZTF J2020+5033, a high-mass brown dwarf (BD) transiting a low-mass star with an orbital period of 1.90 hours. Phase-resolved spectroscopy, optical and infrared light curves, and precise astrometry from Gaia allow us to constrain the masses, radii, and temperatures of both components with few-percent precision. We infer a BD mass of $M_{\rm BD} = 80.1\pm 1.6\,M_{\rm J}$, almost exactly at the stellar/substellar boundary, and a moderately inflated radius, $R_{\rm BD} = 1.05\pm 0.02\,R_{\rm J}$. The transiting object's temperature, $T_{\rm eff}\approx 1700\,\rm K$, is well-constrained by the depth of the infrared secondary eclipse and strongly suggests it is a BD. The system's high tangential velocity ($v_\perp = 98\,\rm km\,s^{-1}$) and thick disk-like Galactic orbit imply the binary is old; its close distance ($d\approx 140$ pc) suggests that BDs in short-period orbits are relatively common. ZTF J2020+5033 is the shortest-period known transiting BD by more than a factor of 7. Today, the entire binary would comfortably fit inside the Sun. However, both components must have been considerably larger in youth, implying that the orbit has shrunk by at least a factor of $\sim 5$ since formation. The simplest explanation is that magnetic braking continues to operate efficiently in at least some low-mass stars and BDs. △ Less

Submitted 13 September, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 15 pages, 9 figures. Accepted to OJA

SRGeJ045359.9+622444: A 55-min Period Eclipsing AM CVn Discovered from a Joint SRG/eROSITA + ZTF Search

Authors: Antonio C. Rodriguez, Ilkham Galiullin, Marat Gilfanov, Shrinivas R. Kulkarni, Irek Khamitov, Ilfan Bikmaev, Jan van Roestel, Lev Yungelson, Kareem El-Badry, Rashid Sunayev, Thomas A. Prince, Mikhail Buntov, Ilaria Caiazzo, Andrew Drake, Mark Gorbachev, Matthew J. Graham, Rustam Gumerov, Eldar Irtuganov, Russ R. Laher, Frank J. Masci, Pavel Medvedev, Josiah Purdum, Nail Sakhibullin, Alexander Sklyanov, Roger Smith , et al. (2 additional authors not shown)

Abstract: AM CVn systems are ultra-compact binaries where a white dwarf accretes from a helium-rich degenerate or semi-degenerate donor. Some AM CVn systems will be among the loudest sources of gravitational waves for the upcoming Laser Interferometer Space Antenna (LISA), yet the formation channel of AM CVns remains uncertain. We report the study and characterisation of a new eclipsing AM CVn, SRGeJ045359.… ▽ More AM CVn systems are ultra-compact binaries where a white dwarf accretes from a helium-rich degenerate or semi-degenerate donor. Some AM CVn systems will be among the loudest sources of gravitational waves for the upcoming Laser Interferometer Space Antenna (LISA), yet the formation channel of AM CVns remains uncertain. We report the study and characterisation of a new eclipsing AM CVn, SRGeJ045359.9+622444 (hereafter SRGeJ0453), discovered from a joint SRG/eROSITA and ZTF program to identify cataclysmic variables (CVs). We obtained optical photometry to confirm the eclipse of SRGeJ0453 and determine the orbital period to be $P_\textrm{orb} = 55.0802 \pm 0.0003$ min. We constrain the binary parameters by modeling the high-speed photometry and radial velocity curves and find $M_\textrm{donor} = 0.044 \pm0.024 M_{\odot}$ and $R_\textrm{donor}=0.078 \pm 0.012 R_{\odot}$. The X-ray spectrum is approximated by a power-law model with an unusually flat photon index of $Γ\sim 1$ previously seen in magnetic CVs with SRG/eROSITA, but verifying the magnetic nature of SRGeJ0453 requires further investigation. Optical spectroscopy suggests that the donor star of SRGeJ0453 could have initially been a He star or a He white dwarf. SRGeJ0453 is the ninth eclipsing AM CVn system published to date, and its lack of optical outbursts have made it elusive in previous surveys. The discovery of SRGeJ0453 using joint X-ray and optical surveys highlights the potential for discovering similar systems in the near future. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: Submitted to ApJ

Young Star Clusters Dominate the Production of Detached Black Hole-Star Binaries

Authors: Ugo Niccolò Di Carlo, Poojan Agrawal, Carl L. Rodriguez, Katelyn Breivik

Abstract: The recent discovery of two detached black hole-star (BH-star) binaries from Gaia's third data release has sparkled interest in understanding the formation mechanisms of these systems. We investigate the formation of these systems by dynamical processes in young open star clusters (SCs) and via isolated binary (IB) evolution, using a combination of direct $N$-body models and population synthesis s… ▽ More The recent discovery of two detached black hole-star (BH-star) binaries from Gaia's third data release has sparkled interest in understanding the formation mechanisms of these systems. We investigate the formation of these systems by dynamical processes in young open star clusters (SCs) and via isolated binary (IB) evolution, using a combination of direct $N$-body models and population synthesis simulations. By comparing dynamical and isolated systems created using the same model of binary stellar evolution, we find that dynamical formation in SCs is nearly 40 times as efficient per unit of star formation at producing BH-star binaries compared to IB evolution. We expand this analysis to the full Milky Way (MW) using a FIRE-2 hydrodynamical simulation of a MW-mass galaxy. Even assuming that only $10\%$ of star formation produces SCs with masses $> 1000\,\mathrm{M_{\odot}}$, we find that the MW contains $\sim 2 \times 10^5$ BH-star systems, with approximately 4 out of every 5 systems being formed dynamically. Many of these dynamically-formed systems have larger orbital periods, eccentricities, and black hole masses than their isolated counterparts. For binaries older than 100 Myr, we show that any detectable system with $e\gtrsim0.5$ or $M_{\rm BH}\gtrsim 10\,\mathrm{M_{\odot}}$ can only be formed through dynamical processes. Our MW model predicts between 61 and 210 such detections from the complete DR4 Gaia catalog, with the majority of systems being dynamically formed in massive and metal-rich SCs. Finally, we compare our populations to the recently discovered Gaia BH1 and Gaia BH2, and conclude that the dynamical scenario is the most favorable formation pathway for both systems. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 14 pages, 8 figures, 2 tables. Submitted to ApJ, comments welcome

The fastest stars in the Galaxy

Authors: Kareem El-Badry, Ken J. Shen, Vedant Chandra, Evan Bauer, Jim Fuller, Jay Strader, Laura Chomiuk, Rohan Naidu, Ilaria Caiazzo, Antonio C. Rodriguez, Pranav Nagarajan, Natsuko Yamaguchi, Zachary P. Vanderbosch, Benjamin R. Roulston, Jan van Roestel, Boris Gänsicke, Jiwon Jesse Han, Kevin B. Burdge, Alexei V. Filippenko, Thomas G. Brink, WeiKang Zheng

Abstract: We report a spectroscopic search for hypervelocity white dwarfs (WDs) that are runaways from Type Ia supernovae (SNe Ia) and related thermonuclear explosions. Candidates are selected from Gaia data with high tangential velocities and blue colors. We find six new runaways, including four stars with radial velocities (RVs) $>1000\,\rm km\,s^{-1}$ and total space velocities… ▽ More We report a spectroscopic search for hypervelocity white dwarfs (WDs) that are runaways from Type Ia supernovae (SNe Ia) and related thermonuclear explosions. Candidates are selected from Gaia data with high tangential velocities and blue colors. We find six new runaways, including four stars with radial velocities (RVs) $>1000\,\rm km\,s^{-1}$ and total space velocities $\gtrsim 1300\,\rm km\,s^{-1}$. These are most likely the surviving donors from double-degenerate binaries in which the other WD exploded. The other two objects have lower minimum velocities, $\gtrsim 600\,\rm km\,s^{-1}$, and may have formed through a different mechanism, such as pure deflagration of a WD in a Type Iax supernova. The four fastest stars are hotter and smaller than the previously known "D$^6$ stars," with effective temperatures ranging from $\sim$20,000 to $\sim$130,000 K and radii of $\sim 0.02-0.10\,R_{\odot}$. Three of these have carbon-dominated atmospheres, and one has a helium-dominated atmosphere. Two stars have RVs of $-1694$ and $-2285\rm \,km\,s^{-1}$ -- the fastest systemic stellar RVs ever measured. Their inferred birth velocities, $\sim 2200-2500\,\rm km\,s^{-1}$, imply that both WDs in the progenitor binary had masses $>1.0\,M_{\odot}$. The high observed velocities suggest that a dominant fraction of the observed hypervelocity WD population comes from double-degenerate binaries whose total mass significantly exceeds the Chandrasekhar limit. However, the two nearest and faintest D$^6$ stars have the lowest velocities and masses, suggesting that observational selection effects favor rarer, higher-mass stars. A significant population of fainter low-mass runaways may still await discovery. We infer a birth rate of D$^6$ stars that is consistent with the SN Ia rate. The birth rate is poorly constrained, however, because the luminosities and lifetimes of $\rm D^6$ stars are uncertain. △ Less

Submitted 25 July, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: 26 pages, 17 figures. Accepted to OJA

Asteroseismology using quadrupolar f-modes revisited: breaking of universal relationships in the slow hadron-quark conversion scenario

Authors: Ignacio F. Ranea-Sandoval, Mauro Mariani, Marcos O. Celi, M. Camila Rodríguez, Lucas Tonetto

Abstract: In this work, we consider polar perturbations and we calculate the frequency and damping time of the quadrupolar fundamental f -mode of compact objects, constructed using a wide range of model-independent hybrid equations of state that include quark matter. We give special attention to the impact of the hadron-quark conversion speed that, in the slow case, gives rise to a branch of slow stable hyb… ▽ More In this work, we consider polar perturbations and we calculate the frequency and damping time of the quadrupolar fundamental f -mode of compact objects, constructed using a wide range of model-independent hybrid equations of state that include quark matter. We give special attention to the impact of the hadron-quark conversion speed that, in the slow case, gives rise to a branch of slow stable hybrid stars. Moreover, we study the validity of universal relationships proposed in the literature and find out that none of them remains valid when slow stable hybrid stars are taken into account. This fact could constrain the applicability of asteroseismology methods with fundamental modes designed to estimate the properties of pulsating compact objects. We hope that this result could be tested with the start up of the third-generation gravitational wave observatories, which might shed some light on the f -mode emission from compact objects. △ Less

Submitted 28 June, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: 13 pages, 13 figures. Comments are welcome

Journal ref: Phys. Rev. D 107, 123028 2023

The geometry of the hot corona in MCG-05-23-16 constrained by X-ray polarimetry

Authors: D. Tagliacozzo, A. Marinucci, F. Ursini, G. Matt, S. Bianchi, L. Baldini, T. Barnouin, N. Cavero Rodriguez, A. De Rosa, L. Di Gesu, M. Dovciak, D. Harper, A. Ingram, V. Karas, D. E. Kim, H. Krawczynski, G. Madejski, F. Marin, R. Middei, H. L. Marshall, F. Muleri, C. Panagiotou, P. O. Petrucci, J. Podgorny, J. Poutanen , et al. (83 additional authors not shown)

Abstract: We report on the second observation of the radio-quiet active galactic nucleus (AGN) MCG-05-23-16 performed with the Imaging X-ray Polarimetry Explorer (IXPE). The observation started on 2022 November 6 for a net observing time of 640 ks, and was partly simultaneous with NuSTAR (86 ks). After combining these data with those obtained in the first IXPE pointing on May 2022 (simultaneous with XMM-New… ▽ More We report on the second observation of the radio-quiet active galactic nucleus (AGN) MCG-05-23-16 performed with the Imaging X-ray Polarimetry Explorer (IXPE). The observation started on 2022 November 6 for a net observing time of 640 ks, and was partly simultaneous with NuSTAR (86 ks). After combining these data with those obtained in the first IXPE pointing on May 2022 (simultaneous with XMM-Newton and NuSTAR) we find a 2-8 keV polarization degree $Π$ = 1.6 $\pm$ 0.7 (at 68 per cent confidence level), which corresponds to an upper limit $Π$ = 3.2 per cent (at 99 per cent confidence level). We then compare the polarization results with Monte Carlo simulations obtained with the MONK code, with which different coronal geometries have been explored (spherical lamppost, conical, slab and wedge). Furthermore, the allowed range of inclination angles is found for each geometry. If the best fit inclination value from a spectroscopic analysis is considered, a cone-shaped corona along the disc axis is disfavoured. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: 9 pages, 8 figures, 3 table. Submitted to MNRAS

Large-scale detector testing for the GAPS Si(Li) Tracker

Authors: Mengjiao Xiao, Achim Stoessl, Brandon Roach, Cory Gerrity, Ian Bouche, Gabriel Bridges, Philip von Doetinchem, Charles J. Hailey, Derik Kraych, Anika Katt, Michael Law, Alexander Lowell, Evan Martinez, Kerstin Perez, Maggie Reed, Chelsea Rodriguez, Nathan Saffold, Ceaser Stringfield, Hershel Weiner, Kelsey Yee

Abstract: Lithium-drifted silicon [Si(Li)] has been used for decades as an ionizing radiation detector in nuclear, particle, and astrophysical experiments, though such detectors have frequently been limited to small sizes (few cm$^2$) and cryogenic operating temperatures. The 10-cm-diameter Si(Li) detectors developed for the General Antiparticle Spectrometer (GAPS) balloon-borne dark matter experiment are n… ▽ More Lithium-drifted silicon [Si(Li)] has been used for decades as an ionizing radiation detector in nuclear, particle, and astrophysical experiments, though such detectors have frequently been limited to small sizes (few cm$^2$) and cryogenic operating temperatures. The 10-cm-diameter Si(Li) detectors developed for the General Antiparticle Spectrometer (GAPS) balloon-borne dark matter experiment are novel particularly for their requirements of low cost, large sensitive area (~10 m$^2$ for the full 1440-detector array), high temperatures (near -40$\,^\circ$C), and energy resolution below 4 keV FWHM for 20--100-keV x-rays. Previous works have discussed the manufacturing, passivation, and small-scale testing of prototype GAPS Si(Li) detectors. Here we show for the first time the results from detailed characterization of over 1100 flight detectors, illustrating the consistent intrinsic low-noise performance of a large sample of GAPS detectors. This work demonstrates the feasibility of large-area and low-cost Si(Li) detector arrays for next-generation astrophysics and nuclear physics applications. △ Less

Submitted 7 September, 2023; v1 submitted 29 April, 2023; originally announced May 2023.

Comments: Updated to version accepted in IEEE Trans Nucl Sci. Minor changes to text, fixed plotting error on Fig. 5. Conclusions unchanged

Journal ref: IEEE Trans. Nucl. Sci. 70 (8), 2023, pp. 2125-2133

Sodium enhancement in evolved cataclysmic variables

Authors: Natsuko Yamaguchi, Kareem El-Badry, Antonio C. Rodriguez, Maude Gull, Benjamin R. Roulston, Zachary P. Vanderbosch

Abstract: We present follow-up spectroscopy of 21 cataclysmic variables (CVs) with evolved secondaries and ongoing or recently-terminated mass transfer. Evolutionary models predict that the secondaries should have anomalous surface abundances owing to nuclear burning in their cores during their main-sequence evolution and subsequent envelope stripping by their companion white dwarfs. To test these models, w… ▽ More We present follow-up spectroscopy of 21 cataclysmic variables (CVs) with evolved secondaries and ongoing or recently-terminated mass transfer. Evolutionary models predict that the secondaries should have anomalous surface abundances owing to nuclear burning in their cores during their main-sequence evolution and subsequent envelope stripping by their companion white dwarfs. To test these models, we measure sodium (Na) abundances of the donors from the Fraunhofer "D" doublet. Accounting for interstellar absorption, we find that {\it all} objects in our sample have enhanced Na abundances. We measure 0.3 $\lesssim$ [Na/H] $\lesssim$ 1.5 dex across the sample, with a median [Na/H] = 0.956 dex, i.e., about an order of magnitude enhancement over solar values. To interpret these values, we run MESA binary evolution models of CVs in which mass transfer begins just as the donor leaves the main sequence. These generically predict Na enhancement in donors with initial donor masses $\gtrsim 1\,M_{\odot}$, consistent with our observations. In the models, Na enrichment occurs in the donors' cores via the NeNa cycle near the end of their main-sequence evolution. Na-enhanced material is exposed when the binaries reach orbital periods of a few hours. Donors with higher initial masses are predicted to have higher Na abundances at fixed orbital period owing to their higher core temperatures during main-sequence evolution. The observed [Na/H] values are on average $\approx$0.3 dex higher than predicted by the models. Surface abundances of evolved CV donors provide a unique opportunity to study nuclear burning products in the cores of intermediate-mass stars. △ Less

Submitted 20 June, 2023; v1 submitted 26 April, 2023; originally announced April 2023.

Comments: accepted to MNRAS

Spectroscopic follow-up of black hole and neutron star candidates in ellipsoidal variables from Gaia DR3

Authors: Pranav Nagarajan, Kareem El-Badry, Antonio C. Rodriguez, Jan van Roestel, Benjamin Roulston

Abstract: We present multi-epoch spectroscopic follow-up of a sample of ellipsoidal variables selected from Gaia DR3 as candidates for hosting quiescent black holes (BHs) and neutron stars (NSs). Our targets were identified as BH/NS candidates because their optical light curves -- when interpreted with models that attribute variability to tidal distortion of a star by a companion that contributes negligible… ▽ More We present multi-epoch spectroscopic follow-up of a sample of ellipsoidal variables selected from Gaia DR3 as candidates for hosting quiescent black holes (BHs) and neutron stars (NSs). Our targets were identified as BH/NS candidates because their optical light curves -- when interpreted with models that attribute variability to tidal distortion of a star by a companion that contributes negligible light -- suggest that the companions are compact objects. From the likely BH/NS candidates identified in recent work accompanying Gaia DR3, we select 14 of the most promising targets for follow-up. We obtained spectra for each object at 2-10 epochs, strategically observing near conjunction to best-constrain the radial velocity semi-amplitude. From the measured semi-amplitudes of the radial velocity curves, we derive minimum companion masses of $M_{2,\min} \leq 0.5 ~ M_{\odot}$ in all cases. Assuming random inclinations, the typical inferred companion mass is $M_2 \sim 0.15 ~ M_{\odot}$. This makes it unlikely that any of these systems contain a BH or NS, and we consider alternative explanations for the observed variability. We can best reproduce the observed light curves and radial velocities with models for unequal-mass contact binaries with starspots. Some of the objects in our sample may also be detached main-sequence binaries, or even single stars with pulsations or starspot variability masquerading as ellipsoidal variation. We provide recommendations for future spectroscopic efforts to further characterize this sample and more generally to search for compact object companions in close binaries. △ Less

Submitted 27 July, 2023; v1 submitted 14 April, 2023; originally announced April 2023.

Comments: 18 pages, 12 figures, Published in MNRAS

Three approaches for the classification of protoneutron star oscillation modes

Authors: M. C. Rodriguez, Ignacio F. Ranea-Sandoval, C. Chirenti, D. Radice

Abstract: The future detection of gravitational waves (GWs) from a galactic core-collapse supernova will provide information on the physics inside protoneutron stars (PNS). In this work, we apply three different classification methods for the PNS non-radial oscillation modes: Cowling classification, Generalized Cowling Nomenclature (GCN), and a Classification Based on Modal Properties (CBMP). Using PNS mode… ▽ More The future detection of gravitational waves (GWs) from a galactic core-collapse supernova will provide information on the physics inside protoneutron stars (PNS). In this work, we apply three different classification methods for the PNS non-radial oscillation modes: Cowling classification, Generalized Cowling Nomenclature (GCN), and a Classification Based on Modal Properties (CBMP). Using PNS models from $3$D simulations of core-collapse supernovae, we find that in the early stages of the PNS evolution, typically before $0.4$ seconds after the bounce, the Cowling classification is inconsistent, but the GCN and the CBMP provide complementary information that helps to understand the evolution of the modes. In the GCN, we note several avoided crossings as the mode frequencies evolve at early times, while the CBMP tracks the modes across the avoided crossings. We verify that the strongest emission of GWs by the PNS corresponds to the $f$-mode in the GCN, indicating that the mode trapping region alternates between the core and the envelope at each avoided crossing. At later times, approximately $0.4$ seconds after the bounce, the three classification methods present a similar description of the mode spectrum. We use our results to test universal relations for the PNS modes according to their classification and find that the behaviour of the universal relations for $f$- and $p$-modes is remarkably simple in the CBMP. △ Less

Submitted 12 May, 2023; v1 submitted 31 March, 2023; originally announced April 2023.

Comments: 11 pages, 8 figures. Matches the version accepted on MNRAS

Designing wavelength sampling for Fabry-Pérot observations. Information-based spectral sampling

Authors: C. J. Díaz Baso, L. Rouppe van der Voort, J. de la Cruz Rodríguez, J. Leenaarts

Abstract: Fabry-Pérot interferometers (FPIs) have become very popular in solar observations because they offer a balance between cadence, spatial resolution, and spectral resolution through a careful design of the spectral sampling scheme according to the observational requirements of a given target. However, an efficient balance requires knowledge of the expected target conditions, the properties of the ch… ▽ More Fabry-Pérot interferometers (FPIs) have become very popular in solar observations because they offer a balance between cadence, spatial resolution, and spectral resolution through a careful design of the spectral sampling scheme according to the observational requirements of a given target. However, an efficient balance requires knowledge of the expected target conditions, the properties of the chosen spectral line, and the instrumental characteristics. Our aim is to find a method that allows finding the optimal spectral sampling of FPI observations in a given spectral region. In this study, we propose a technique based on a sequential selection approach where a neural network is used to predict the spectrum (or physical quantities, if the model is known) from the information at a few points. Only those points that contain relevant information and improve the model prediction are included in the sampling scheme. The method adapts the separation of the points according to the spectral resolution of the instrument, the typical broadening of the spectral shape, and the typical Doppler velocities. The experiments using the CaII 8542 A line show that the resulting wavelength scheme naturally places more points in the core than in the wings, consistent with the sensitivity of the spectral line at each wavelength interval. The method can also be used as an accurate interpolator, to improve the inference of the magnetic field when using the weak-field approximation. Overall, this method offers an objective approach for designing new instrumentation or observing proposals with customized configurations for specific targets. This is particularly relevant when studying highly dynamic events in the solar atmosphere with a cadence that preserves spectral coherence without sacrificing much information. △ Less

Submitted 24 March, 2023; originally announced March 2023.

Comments: 8 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 673, A35 (2023)

Modelling stellar evolution in mass-transferring binaries and gravitational-wave progenitors with METISSE

Authors: Poojan Agrawal, Jarrod Hurley, Simon Stevenson, Carl L. Rodriguez, Dorottya Szecsi, Alex Kemp

Abstract: Massive binaries are vital sources of various transient processes, including gravitational-wave mergers. However, large uncertainties in the evolution of massive stars, both physical and numerical, present a major challenge to the understanding of their binary evolution. In this paper, we upgrade our interpolation-based stellar evolution code METISSE to include the effects of mass changes, such as… ▽ More Massive binaries are vital sources of various transient processes, including gravitational-wave mergers. However, large uncertainties in the evolution of massive stars, both physical and numerical, present a major challenge to the understanding of their binary evolution. In this paper, we upgrade our interpolation-based stellar evolution code METISSE to include the effects of mass changes, such as binary mass transfer or wind-driven mass loss, not already included within the input stellar tracks. METISSE's implementation of mass loss (applied to tracks without mass loss) shows excellent agreement with the SSE fitting formulae and with detailed MESA tracks, except in cases where the mass transfer is too rapid for the star to maintain equilibrium. We use this updated version of METISSE within the binary population synthesis code BSE to demonstrate the impact of varying stellar evolution parameters, particularly core overshooting, on the evolution of a massive (25M$_\odot$ and 15M$_\odot$) binary system with an orbital period of 1800 days. Depending on the input tracks, we find that the binary system can form a binary black hole or a black hole-neutron star system, with primary(secondary) remnant masses ranging between 4.47(1.36)M$_\odot$ and 12.30(10.89)M$_\odot$, and orbital periods ranging from 6 days to the binary becoming unbound. Extending this analysis to a population of isolated binaries uniformly distributed in mass and orbital period, we show that the input stellar models play an important role in determining which regions of the binary parameter space can produce compact binary mergers, paving the way for predictions for current and future gravitational-wave observatories. △ Less

Submitted 28 July, 2023; v1 submitted 17 March, 2023; originally announced March 2023.

Comments: 19 pages, 14 figures, accepted for publication in MNRAS

A reconnection driven magnetic flux cancellation and a quiet Sun Ellerman bomb

Authors: Anjali. J. Kaithakkal, J. M. Borrero, A. Pastor Yabar, J. de la Cruz Rodríguez

Abstract: The focus of this investigation is to quantify the conversion of magnetic to thermal energy initiated by a quiet Sun cancellation event and to explore the resulting dynamics from the interaction of the opposite polarity magnetic features. We used imaging spectroscopy in the H$α$ line, along with spectropolarimetry in the \ion{Fe}{I} 6173~Å and \ion{Ca}{II} 8542~Å lines from the Swedish Solar Teles… ▽ More The focus of this investigation is to quantify the conversion of magnetic to thermal energy initiated by a quiet Sun cancellation event and to explore the resulting dynamics from the interaction of the opposite polarity magnetic features. We used imaging spectroscopy in the H$α$ line, along with spectropolarimetry in the \ion{Fe}{I} 6173~Å and \ion{Ca}{II} 8542~Å lines from the Swedish Solar Telescope (SST) to study a reconnection-related cancellation and the appearance of a quiet Sun Ellerman bomb (QSEB). We observed, for the first time, QSEB signature in both the wings and core of the \ion{Fe}{I} 6173~Å line. We also found that, at times, the \ion{Fe}{I} line-core intensity reaches higher values than the quiet Sun continuum intensity. From FIRTEZ-dz inversions of the Stokes profiles in \ion{Fe}{I} and \ion{Ca}{II} lines, we found enhanced temperature, with respect to the quiet Sun values, at the photospheric ($\logτ_c$ = -1.5; $\sim$1000 K) and lower chromospheric heights ($\logτ_c$ = -4.5; $\sim$360 K). From the calculation of total magnetic energy and thermal energy within these two layers it was confirmed that the magnetic energy released during the flux cancellation can support heating in the aforesaid height range. Further, the temperature stratification maps enabled us to identify cumulative effects of successive reconnection on temperature pattern, including recurring temperature enhancements. Similarly, Doppler velocity stratification maps revealed impacts on plasma flow pattern, such as a sudden change in the flow direction. △ Less

Submitted 11 March, 2023; originally announced March 2023.

Comments: 16 pages, 12 figures, accepted for publication in MNRAS

Decay of a photospheric transient filament at the boundary of a pore and the chromospheric response

Authors: Philip Lindner, Rolf Schlichenmaier, Nazaret Bello González, Jaime de la Cruz Rodríguez

Abstract: Intermediate stages between pores and sunspots are a rare phenomenon and can manifest with the formation of transient photospheric penumbral-like filaments. Although the magnetic field changes rapidly during the evolution of such filaments, they have not been shown to be connected to magnetic reconnection events yet. We analyzed observations of a pore in NOAA AR 12739 from the Swedish Solar Telesc… ▽ More Intermediate stages between pores and sunspots are a rare phenomenon and can manifest with the formation of transient photospheric penumbral-like filaments. Although the magnetic field changes rapidly during the evolution of such filaments, they have not been shown to be connected to magnetic reconnection events yet. We analyzed observations of a pore in NOAA AR 12739 from the Swedish Solar Telescope including spectropolarimetric data of the Fe I 6173 Å and the Ca II 8542 Å line and spectroscopic data of the Ca II K 3934 Å line. The VFISV Milne-Eddington inversion code and the multi-line Non-LTE inversion code STiC were utilized to obtain atmospheric parameters in the photosphere and the chromosphere. Multiple filamentary structures of inclined magnetic fields are found in photospheric inclination maps at the boundary of the pore, although the pore never developed a penumbra. One of the filaments shows a clear counterpart in continuum intensity maps in addition to photospheric blueshifts. During its decay, a brightening in the blue wing of the Ca II 8542 Å line is observed. The Ca II K 3934 Å and the Ca II 8542 Å lines show complex spectral profiles in this region. Depth-dependent STiC inversion results using data from all available lines yield a temperature increase (roughly 1000 Kelvin) and bidirectional flows (magnitudes up to 8 km/s) at log tau=-3.5. The temporal and spatial correlation of the decaying filament (observed in the photosphere) to the temperature increase and the bidirectional flows in the high photosphere/low chromosphere suggests that they are connected. We propose scenarios in which magnetic reconnection happens at the edge of a rising magnetic flux tube in the photosphere. This leads to both the decay of the filament in the photosphere and the observed temperature increase and the bidirectional flows in the high photosphere/low chromosphere. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: Paper accepted by "Astronomy & Astrophysics" on 02/03/2023

Journal ref: A&A 673, A65 (2023)

Does H$α$ Stokes~$V$ profiles probe the chromospheric magnetic field? An observational perspective

Authors: Harsh Mathur, K. Nagaraju, Jayant Joshi, Jaime de la Cruz Rodríguez

Abstract: We investigated the diagnostic potential of the Stokes $V$ profile of the H$α$ line to probe the chromospheric line-of-sight (LOS) magnetic field ($B_{\mathrm{LOS}}$) by comparing the $B_{\mathrm{LOS}}$ inferred from the weak field approximation (WFA) with that of inferred from the multi-line inversions of the Ca II 8542 Å, Si I 8536 Å and Fe I 8538 Å lines using the STiC inversion code. Simultane… ▽ More We investigated the diagnostic potential of the Stokes $V$ profile of the H$α$ line to probe the chromospheric line-of-sight (LOS) magnetic field ($B_{\mathrm{LOS}}$) by comparing the $B_{\mathrm{LOS}}$ inferred from the weak field approximation (WFA) with that of inferred from the multi-line inversions of the Ca II 8542 Å, Si I 8536 Å and Fe I 8538 Å lines using the STiC inversion code. Simultaneous spectropolarimetric observations of a pore in the Ca II 8542 Å and H$α$ spectral lines obtained from the SPINOR at the Dunn Solar Telescope on the 4th of December, 2008 are used in this study. The WFA was applied on the Stokes $I$ and $V$ profiles of H$α$ line over three wavelength ranges viz.: around line core ($Δλ=\pm0.35$ Å), line wings ($Δλ=[-1.5, -0.6]$ and $[+0.6, +1.5]$ Å) and full spectral range of the line ($Δλ=\pm1.5$ Å) to derive the $B_{\mathrm{LOS}}$. We found the maximum $B_{\mathrm{LOS}}$ strengths of $\sim+800$ and $\sim+600$ G at $\logτ_{\mathrm{500}}$ = $-$1 and $-$4.5, respectively in the pore. The morphological map of the $B_{\mathrm{LOS}}$ inferred from the H$α$ line core is similar to the $B_{\mathrm{LOS}}$ map at $\logτ_{\mathrm{500}}$ = $-$4.5 inferred from multi-line inversions. The $B_{\mathrm{LOS}}$ map inferred from the H$α$ line wings and full spectral range have a similar morphological structure to the $B_{\mathrm{LOS}}$ map inferred at $\logτ_{\mathrm{500}}$ = $-$1. The $B_{\mathrm{LOS}}$ estimated from H$α$ using WFA is weaker by a factor of $\approx 0.53$ than that of inferred from the multi-line inversions. △ Less

Submitted 25 February, 2023; originally announced February 2023.

Comments: Accepted to be published in the ApJ on 25th February, 2023

Constraints on the Cosmological Coupling of Black Holes from the Globular Cluster NGC 3201

Authors: Carl L. Rodriguez

Abstract: Globular clusters are among the oldest stellar populations in the Milky Way; consequently, they also host some of the oldest known stellar-mass black holes, providing insight into black hole formation and evolution in the early ($z\gtrsim 2$) Universe. Recent observations of supermassive black holes in elliptical galaxies have been invoked to suggest the possibility of a cosmological coupling betw… ▽ More Globular clusters are among the oldest stellar populations in the Milky Way; consequently, they also host some of the oldest known stellar-mass black holes, providing insight into black hole formation and evolution in the early ($z\gtrsim 2$) Universe. Recent observations of supermassive black holes in elliptical galaxies have been invoked to suggest the possibility of a cosmological coupling between astrophysical black holes and the surrounding expanding Universe, offering a mechanism for black holes to grow over cosmic time, and potentially explaining the origin of dark energy. In this paper, I show that the mass functions of the two radial velocity black hole candidates in NGC 3201 place strong constraints on the cosmologically-coupled growth of black holes. In particular, the amount of coupling required to explain the origin of dark energy would either require both NGC 3201 black holes to be nearly face on (a configuration with probability of at most $10^{-4}$) or one of the BHs would need to have formed with a mass below that of the most massive neutron stars ($2.2M_{\odot}$). This emphasizes that these and other detached black hole-star binaries can serve not only as laboratories for compact object and binary astrophysics, but as constraints on the long-term evolution of astrophysical black holes. △ Less

Submitted 3 April, 2023; v1 submitted 23 February, 2023; originally announced February 2023.

Comments: 9 Pages, 2 Figures, matches version accepted by ApJL

Ultra-high resolution observations of plasmoid-mediated magnetic reconnection in the deep solar atmosphere

Authors: L. Rouppe van der Voort, M. van Noort, J. de la Cruz Rodriguez

Abstract: Magnetic reconnection in the deep solar atmosphere can give rise to enhanced emission in the Balmer hydrogen lines, a phenomenon referred to as Ellerman bombs. To effectively trace magnetic reconnection below the canopy of chromospheric fibrils, we analyzed unique spectroscopic observations of Ellerman bombs in the H-alpha line. We analyzed a 10 min dataset of a young emerging active region observ… ▽ More Magnetic reconnection in the deep solar atmosphere can give rise to enhanced emission in the Balmer hydrogen lines, a phenomenon referred to as Ellerman bombs. To effectively trace magnetic reconnection below the canopy of chromospheric fibrils, we analyzed unique spectroscopic observations of Ellerman bombs in the H-alpha line. We analyzed a 10 min dataset of a young emerging active region observed with the prototype of the Microlensed Hyperspectral Imager (MiHI) at the Swedish 1-m Solar Telescope (SST). The MiHI instrument is an integral field spectrograph that is capable of achieving simultaneous ultra-high resolution in the spatial, temporal and spectral domains. With the combination of the SST adaptive optics system and image restoration techniques, MiHI can deliver diffraction limited observations if the atmospheric seeing conditions allow. The dataset samples the H-alpha line over 4.5 A with 10 mA/pix, with 0.065"/pix over a field of view of 8.6" x 7.7", and at a temporal cadence of 1.33s. This constitutes a hyperspectral data cube that measures 132 x 118 spatial pixels, 456 spectral pixels, and 455 time steps. There were multiple sites with Ellerman bomb activity associated with strong magnetic flux emergence. The Ellerman bomb activity is very dynamic, showing rapid variability and small-scale substructure. We found a number of plasmoid-like blobs with full-width-half-maximum sizes between 0.1" - 0.4" and moving with apparent velocities between 14 and 77 km/s. Some of these blobs have Ellerman bomb spectral profiles with a single peak at a Doppler offset between 47 and 57 km/s. Our observations support the idea that fast magnetic reconnection in Ellerman bombs is mediated by the formation of plasmoids. These MiHI observations demonstrate that a micro-lens based integral field spectrograph is capable of probing fundamental physical processes in the solar atmosphere. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: accepted for publication in A&A. Movies can be found at https://www.mn.uio.no/astro/english/people/aca/rouppe/movies/

Journal ref: A&A 673, A11 (2023)

A red giant orbiting a black hole

Authors: Kareem El-Badry, Hans-Walter Rix, Yvette Cendes, Antonio C. Rodriguez, Charlie Conroy, Eliot Quataert, Keith Hawkins, Eleonora Zari, Melissa Hobson, Katelyn Breivik, Arne Rau, Edo Berger, Sahar Shahaf, Rhys Seeburger, Kevin B. Burdge, David W. Latham, Lars A. Buchhave, Allyson Bieryla, Dolev Bashi, Tsevi Mazeh, Simchon Faigler

Abstract: We report spectroscopic and photometric follow-up of a dormant black hole (BH) candidate from Gaia DR3. The system, which we call Gaia BH2, contains a $\sim 1M_{\odot}$ red giant and a dark companion with mass $M_2 = 8.9\pm 0.3\,M_{\odot}$ that is very likely a BH. The orbital period, $P_{\rm orb} = 1277$ days, is much longer than that of any previously studied BH binary. Our radial velocity (RV)… ▽ More We report spectroscopic and photometric follow-up of a dormant black hole (BH) candidate from Gaia DR3. The system, which we call Gaia BH2, contains a $\sim 1M_{\odot}$ red giant and a dark companion with mass $M_2 = 8.9\pm 0.3\,M_{\odot}$ that is very likely a BH. The orbital period, $P_{\rm orb} = 1277$ days, is much longer than that of any previously studied BH binary. Our radial velocity (RV) follow-up over a 7-month period spans more than 90% of the orbit's dynamic range in RV and is in excellent agreement with predictions of the Gaia solution. UV imaging and high-resolution optical spectra rule out all plausible luminous companions that could explain the orbit. The star is a bright ($G=12.3$), slightly metal-poor ($\rm [Fe/H]=-0.22$) low-luminosity giant ($T_{\rm eff}=4600\,\rm K$; $R = 7.8\,R_{\odot}$; $\log\left[g/\left({\rm cm\,s^{-2}}\right)\right] = 2.6$). The binary's orbit is moderately eccentric ($e=0.52$). The giant is strongly enhanced in $α-$elements, with $\rm [α/Fe] = +0.26$, but the system's Galactocentric orbit is typical of the thin disk. We obtained X-ray and radio nondetections of the source near periastron, which support BH accretion models in which the net accretion rate at the horizon is much lower than the Bondi-Hoyle-Lyttleton rate. At a distance of 1.16 kpc, Gaia BH2 is the second-nearest known BH, after Gaia BH1. Its orbit -- like that of Gaia BH1 -- seems too wide to have formed through common envelope evolution. Gaia BH1 and BH2 have orbital periods at opposite edges of the Gaia DR3 sensitivity curve, perhaps hinting at a bimodal intrinsic period distribution for wide BH binaries. Dormant BH binaries like Gaia BH1 and Gaia BH2 likely significantly outnumber their close, X-ray bright cousins, but their formation pathways remain uncertain. △ Less

Submitted 19 March, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Comments: Accepted to MNRAS

Estimating the longitudinal magnetic field in the chromosphere of quiet-Sun magnetic concentrations

Authors: S. Esteban Pozuelo, A. Asensio Ramos, J. de la Cruz Rodríguez, J. Trujillo Bueno, M. J. Martínez González

Abstract: Details of the magnetic field in the quiet Sun chromosphere are key to our understanding of essential aspects of the solar atmosphere. We aim to determine the longitudinal magnetic field component (B_lon) of quiet Sun regions depending on their size. We estimated B_lon by applying the weak-field approximation (WFA) to high-spatial-resolution Ca II 854.2 nm data taken with the Swedish 1m Solar Tele… ▽ More Details of the magnetic field in the quiet Sun chromosphere are key to our understanding of essential aspects of the solar atmosphere. We aim to determine the longitudinal magnetic field component (B_lon) of quiet Sun regions depending on their size. We estimated B_lon by applying the weak-field approximation (WFA) to high-spatial-resolution Ca II 854.2 nm data taken with the Swedish 1m Solar Telescope. Specifically, we analyzed the estimates inferred for different spectral ranges using the data at the original cadence and temporally integrated signals. The longitudinal magnetic field in each considered plasma structure correlates with its size. Using a spectral range restricted to the line core leads to chromospheric longitudinal fields varying from 50 G at the edges to 150-500 G at the center of the structure. These values increase as the spectral range widens due to the photospheric contribution. However, the difference between this contribution and the chromospheric one is not uniform for all structures. Small and medium-sized concentrations show a steeper height gradient in B_lon compared to their chromospheric values, so estimates for wider ranges are less trustworthy. Signal addition does not alleviate this situation as the height gradients in B_lon are consistent with time. Finally, despite the amplified noise levels that deconvolving processes may cause, data restored with the destretching technique show similar results, though are affected by smearing. We obtained B_lon estimates similar to those previously found, except for large concentrations and wide spectral ranges. In addition, we report a correlation between the height variation of B_lon compared to the chromospheric estimates and the concentration size. This correlation affects the difference between the photospheric and chromospheric magnetic flux values and the reliability of the estimates for wider spectral ranges. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: 14 pages, 14 figures. Accepted for publication in Astronomy & Astrophysics; abstract has been abridged

Journal ref: A&A 672, A141 (2023)

Runaway and Hypervelocity Stars from Compact Object Encounters in Globular Clusters

Authors: Tomás Cabrera, Carl L. Rodriguez

Abstract: The dense environments in the cores of globular clusters (GCs) facilitate many strong dynamical encounters among stellar objects. These encounters have been shown capable of ejecting stars from the host GC, whereupon they become runaway stars, or hypervelocity stars if unbound to the galactic potential. We study high speed stellar ejecta originating from GCs by using Monte Carlo N-body models, in… ▽ More The dense environments in the cores of globular clusters (GCs) facilitate many strong dynamical encounters among stellar objects. These encounters have been shown capable of ejecting stars from the host GC, whereupon they become runaway stars, or hypervelocity stars if unbound to the galactic potential. We study high speed stellar ejecta originating from GCs by using Monte Carlo N-body models, in particular focusing on binary-single encounters involving compact objects. We pair our model-discriminated populations with observational catalogs of Milky Way GCs to compose a present-day galactic population of stellar ejecta. We find that these kinds of encounters can accelerate stars to velocities in excess of 2000 km/s, to speeds beyond the previously predicted limits for ejecta from star-only encounters and in the same regime of galactic center ejections. However, the same ejections can only account for 1.5-20% of the total population of stellar runaways, and only 0.0001-1% of hypervelocity stars, with similar relative rates found for runaway white dwarfs. We also provide credible regions for ejecta from 149 Milky Way GCs, which we hope will be useful as supplementary evidence when pairing runaway stars with origin GCs. △ Less

Submitted 15 February, 2024; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 16 pages, 11 figures, accepted for publication in ApJ. Updated Zenodo link

A very luminous jet from the disruption of a star by a massive black hole

Authors: Igor Andreoni, Michael W. Coughlin, Daniel A. Perley, Yuhan Yao, Wenbin Lu, S. Bradley Cenko, Harsh Kumar, Shreya Anand, Anna Y. Q. Ho, Mansi M. Kasliwal, Antonio de Ugarte Postigo, Ana Sagues-Carracedo, Steve Schulze, D. Alexander Kann, S. R. Kulkarni, Jesper Sollerman, Nial Tanvir, Armin Rest, Luca Izzo, Jean J. Somalwar, David L. Kaplan, Tomas Ahumada, G. C. Anupama, Katie Auchettl, Sudhanshu Barway , et al. (56 additional authors not shown)

Abstract: Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jett… ▽ More Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jetted TDE to date is Swift J1644+57, which was discovered in gamma-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical discovery of AT2022cmc, a rapidly fading source at cosmological distance (redshift z=1.19325) whose unique lightcurve transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-rays, sub-millimeter, and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron "afterglow", likely launched by a SMBH with spin $a \gtrsim 0.3$. Using 4 years of Zwicky Transient Facility (ZTF) survey data, we calculate a rate of $0.02 ^{+ 0.04 }_{- 0.01 }$ Gpc$^{-3}$ yr$^{-1}$ for on-axis jetted TDEs based on the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations. Correcting for the beaming angle effects, this rate confirms that about 1% of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in Nature

Chromospheric Heating from Local Magnetic Growth and Ambipolar Diffusion Under Non-Equilibrium Conditions

Authors: Juan Martínez Sykora, Jaime de la Cruz Rodríguez, Milan Gošić, Alberto Sainz Dalda, Viggo H. Hansteen, Bart De Pontieu

Abstract: The heating of the chromosphere in internetwork regions remains one of the foremost open questions in solar physics. In the present study we tackle this old problem by using a very high spatial-resolution simulation of quiet-Sun conditions performed with radiative MHD numerical models and IRIS observations. We have expanded a previously existing 3D radiative MHD numerical model of the solar atmosp… ▽ More The heating of the chromosphere in internetwork regions remains one of the foremost open questions in solar physics. In the present study we tackle this old problem by using a very high spatial-resolution simulation of quiet-Sun conditions performed with radiative MHD numerical models and IRIS observations. We have expanded a previously existing 3D radiative MHD numerical model of the solar atmosphere, which included self-consistently locally driven magnetic amplification in the chromosphere, by adding ambipolar diffusion and time-dependent non-equilibrium hydrogen ionization to the model. The energy of the magnetic field is dissipated in the upper chromosphere, providing a large temperature increase due to ambipolar diffusion and the non-equilibrium ionization (NEQI). At the same time, we find that adding the ambipolar diffusion and NEQI in the simulation has a minor impact on the local growth of the magnetic field in the lower chromosphere and its dynamics. Our comparison between synthesized Mg II profiles from these high spatial resolution models, with and without ambipolar diffusion and NEQI, and quiet Sun and coronal hole observations from IRIS now reveal a better correspondence. The intensity of profiles is increased and the line cores are slightly broader when ambipolar diffusion and NEQI effects are included. Therefore, the Mg II profiles are closer to those observed than in previous models, but some differences still remain. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: submitted to ApJL, 3 figures

Transverse oscillations in 3D along Ca II K bright fibrils in the Solar chromosphere

Authors: Sepideh Kianfar, Jorrit Leenaarts, Sara Esteban Pozuelo, João M. da Silva Santos, Jaime de la Cruz Rodríguez, Sanja Danilovic

Abstract: Fibrils in the solar chromosphere carry transverse oscillations as determined from non-spectroscopic imaging data. They are estimated to carry an energy flux of several kW m$^{-2}$, which is a significant fraction of the average chromospheric radiative energy losses. We aim to determine oscillation properties of fibrils not only in the plane-of-the-sky (horizontal) direction, but also along the li… ▽ More Fibrils in the solar chromosphere carry transverse oscillations as determined from non-spectroscopic imaging data. They are estimated to carry an energy flux of several kW m$^{-2}$, which is a significant fraction of the average chromospheric radiative energy losses. We aim to determine oscillation properties of fibrils not only in the plane-of-the-sky (horizontal) direction, but also along the line-of-sight (vertical) direction. We obtained imaging-spectroscopy data in Fe I 6173, Ca II 8542, and Ca II K obtained with the Swedish 1-m Solar Telescope. We created a sample of 120 bright Ca II K fibrils and measured their horizontal motions. Their vertical motion was determined through non-LTE inversion of the observed spectra. We determined the periods and velocity amplitudes of the fibril oscillations, as well as phase differences between vertical and horizontal oscillations in the fibrils. The bright Ca II K fibrils carry transverse waves with a mean period of $2.2\times10^2$ s, and a horizontal velocity amplitude of 2 km s$^{-1}$, consistent with earlier results. The mean vertical velocity amplitude is 1 km s$^{-1}$. We find that 118 out of the 120 fibrils carry waves in both the vertical and horizontal directions, and 55 of those have identical periods. For those 55, we find that all phase differences between $0$ and $2π$ occur, with a mild but significant preference for linearly polarized waves (phase difference of $0$ or $π$). The results are consistent with the scenario where transverse waves are excited by granular buffeting at the photospheric footpoints of the fibrils. Estimates of transverse wave flux based on imaging data only are too low because they ignore the contribution of the vertical velocity. △ Less

Submitted 25 October, 2022; originally announced October 2022.

Comments: 12 pages, 12 figures

Prospects and challenges of numerical modelling of the Sun at millimetre wavelengths

Authors: Sven Wedemeyer, Gregory Fleishman, Jaime de la Cruz Rodriguez, Stanislav Gunar, Joao M. da Silva Santos, Patrick Antolin, Juan Camilo Guevara Gomez, Mikolaj Szydlarski, Henrik Eklund

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) offers new diagnostic possibilities that complement other commonly used diagnostics for the study of our Sun. In particular, ALMA's ability to serve as an essentially linear thermometer of the chromospheric gas at unprecedented spatial resolution at millimetre wavelengths and future polarisation measurements have great diagnostic potential. S… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) offers new diagnostic possibilities that complement other commonly used diagnostics for the study of our Sun. In particular, ALMA's ability to serve as an essentially linear thermometer of the chromospheric gas at unprecedented spatial resolution at millimetre wavelengths and future polarisation measurements have great diagnostic potential. Solar ALMA observations are therefore expected to contribute significantly to answering long-standing questions about the structure, dynamics and energy balance of the outer layers of the solar atmosphere. In this regard, current and future ALMA data are also important for constraining and further developing numerical models of the solar atmosphere, which in turn are often vital for the interpretation of observations. The latter is particularly important given the Sun's highly intermittent and dynamic nature that involves a plethora of processes occurring over extended ranges in spatial and temporal scales. Realistic forward modelling of the Sun therefore requires time-dependent three-dimensional radiation magnetohydrodynamics that account for non-equilibrium effects and, typically as a separate step, detailed radiative transfer calculations, resulting in synthetic observables that can be compared to observations. Such artificial observations sometimes also account for instrumental and seeing effects, which, in addition to aiding the interpretation of observations, provide instructive tools for designing and optimising ALMA's solar observing modes. In the other direction, ALMA data in combination with other simultaneous observations enables the reconstruction of the solar atmospheric structure via data inversion techniques. This article highlights central aspects of the impact of ALMA for numerical modelling for the Sun, their potential and challenges, together with selected examples. △ Less

Submitted 25 October, 2022; originally announced October 2022.

Comments: 32 pages, 7 figures. Accepted for publication in Frontiers in Astronomy and Space Sciences (Stellar and Solar Physics)

Status of QUBIC, the Q&U Bolometer for Cosmology

Authors: L. Mousset, P. Ade, A. Almela, G. Amico, L. H. Arnaldi, J. Aumont, S. Banfi, E. S. Battistelli, B. Bélier, L. Bergé, J. -Ph. Bernard, P. de Bernardis, M. Bersanelli, J. Bonaparte, J. D. Bonilla, E. Bunn, D. Buzi, D. Camilieri, F. Cavaliere, P. Chanial, C. Chapron, S. Colombo, F. Columbro, A. Coppolecchia, B. Costanza , et al. (86 additional authors not shown)

Abstract: The Q&U Bolometric Interferometer for Cosmology (QUBIC) is a novel kind of polarimeter optimized for the measurement of the B-mode polarization of the Cosmic Microwave Back-ground (CMB), which is one of the major challenges of observational cosmology. The signal is expected to be of the order of a few tens of nK, prone to instrumental systematic effects and polluted by various astrophysical foregr… ▽ More The Q&U Bolometric Interferometer for Cosmology (QUBIC) is a novel kind of polarimeter optimized for the measurement of the B-mode polarization of the Cosmic Microwave Back-ground (CMB), which is one of the major challenges of observational cosmology. The signal is expected to be of the order of a few tens of nK, prone to instrumental systematic effects and polluted by various astrophysical foregrounds which can only be controlled through multichroic observations. QUBIC is designed to address these observational issues with a novel approach that combines the advantages of interferometry in terms of control of instrumental systematics with those of bolometric detectors in terms of wide-band, background-limited sensitivity. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: Contribution to the 2022 Cosmology session of the 33rd Rencontres de Blois. arXiv admin note: substantial text overlap with arXiv:2203.08947

A Sun-like star orbiting a black hole

Authors: Kareem El-Badry, Hans-Walter Rix, Eliot Quataert, Andrew W. Howard, Howard Isaacson, Jim Fuller, Keith Hawkins, Katelyn Breivik, Kaze W. K. Wong, Antonio C. Rodriguez, Charlie Conroy, Sahar Shahaf, Tsevi Mazeh, Frédéric Arenou, Kevin B. Burdge, Dolev Bashi, Simchon Faigler, Daniel R. Weisz, Rhys Seeburger, Silvia Almada Monter, Jennifer Wojno

Abstract: We report discovery of a bright, nearby ($G = 13.8;\,\,d = 480\,\rm pc$) Sun-like star orbiting a dark object. We identified the system as a black hole candidate via its astrometric orbital solution from the Gaia mission. Radial velocities validated and refined the Gaia solution, and spectroscopy ruled out significant light contributions from another star. Joint modeling of radial velocities and a… ▽ More We report discovery of a bright, nearby ($G = 13.8;\,\,d = 480\,\rm pc$) Sun-like star orbiting a dark object. We identified the system as a black hole candidate via its astrometric orbital solution from the Gaia mission. Radial velocities validated and refined the Gaia solution, and spectroscopy ruled out significant light contributions from another star. Joint modeling of radial velocities and astrometry constrains the companion mass to $M_2 = 9.62\pm 0.18\,M_{\odot}$. The spectroscopic orbit alone sets a minimum companion mass of $M_2>5\,M_{\odot}$; if the companion were a $5\,M_{\odot}$ star, it would be $500$ times more luminous than the entire system. These constraints are insensitive to the mass of the luminous star, which appears as a slowly-rotating G dwarf ($T_{\rm eff}=5850\,\rm K$, $\log g = 4.5$, $M=0.93\,M_{\odot}$), with near-solar metallicity ($\rm [Fe/H] = -0.2$) and an unremarkable abundance pattern. We find no plausible astrophysical scenario that can explain the orbit and does not involve a black hole. The orbital period, $P_{\rm orb}=185.6$ days, is longer than that of any known stellar-mass black hole binary. The system's modest eccentricity ($e=0.45$), high metallicity, and thin-disk Galactic orbit suggest that it was born in the Milky Way disk with at most a weak natal kick. How the system formed is uncertain. Common envelope evolution can only produce the system's wide orbit under extreme and likely unphysical assumptions. Formation models involving triples or dynamical assembly in an open cluster may be more promising. This is the nearest known black hole by a factor of 3, and its discovery suggests the existence of a sizable population of dormant black holes in binaries. Future Gaia releases will likely facilitate the discovery of dozens more. △ Less

Submitted 28 February, 2023; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: 30 pages, 15 figures, 1-2 black holes. Accepted to MNRAS