Search

The hierarchical interplay among gravity, magnetic fields, and turbulent gas flows in delivering the necessary material to form massive protostellar clusters remains enigmatic. We have performed high-resolution (beam size $\sim$14 arcsec $\simeq$ 0.05 pc at a distance 725 pc) 850 $\mu$m dust polarization and C$^{18}$O molecular line observations of Cepheus A (Cep A), the second closest massive star-forming region, using the 15-meter James Clerk Maxwell Telescope (JCMT) along with the SCUBA-2/POL-2 and HARP instruments. Our key analyses reveal that (i) morphologically, all three fields--gravitational (G), magnetic (B), and kinetic (K) fields--are aligned with each other, and (ii) energetically, they exhibit a hierarchical relationship with gravitational ($E_{\mathrm{G}}$) > magnetic ($E_{\mathrm{B}}$) > kinetic ($E_{\mathrm{K}}$). Gravity dominates in Cep A clump and, as a primary active player, dictates the other two agents. Consequently, gravity plays two active roles: (i) induces gas flows and (ii) drags B-field lines toward the gravitational trough. Since magnetic energy dominates kinetic energy, $E_{\mathrm{B}}$ > $E_{\mathrm{K}}$, the "dragged-in" B-field as a secondary active player can mitigate turbulence and instabilities, thereby regularizing gas flows into a more ordered configuration. At the $\sim$0.60 pc clump scale, these flows deliver material at a substantially high rate of $\sim$ 2.1$\times$10$^{-4}$ M$_{\odot}$ yr$^{-1}$ toward the cluster center. Our study, for the first time, presents new insights into how B-fields and turbulent gas flows passively assist the active role of gravity in the formation of a protostellar cluster, contrasting with the standard notion that these agents primarily oppose gravitational collapse., Comment: 33 pages, 6 figures in the main text and 8 figures in the supplementary material, Submitted, and comments are welcome

Many active galaxies harbor powerful relativistic jets, however, the detailed mechanisms of their formation and acceleration remain poorly understood. To investigate the area of jet acceleration and collimation with the highest available angular resolution, we study the innermost region of the bipolar jet in the nearby low-ionization nuclear emission-line region (LINER) galaxy NGC 1052. We combined observations of NGC 1052 taken with VLBA, GMVA, and EHT over one week in the spring of 2017. For the first time, NGC 1052 was detected with the EHT, providing a size of the central region in-between both jet bases of 250 RS (Schwarzschild radii) perpendicular to the jet axes. This size estimate supports previous studies of the jets expansion profile which suggest two breaks of the profile at around 300 RS and 10000 RS distances to the core. Furthermore, we estimated the magnetic field to be 1.25 Gauss at a distance of 22 {\mu}as from the central engine by fitting a synchrotron-self absorption spectrum to the innermost emission feature, which shows a spectral turn-over at about 130 GHz. Assuming a purely poloidal magnetic field, this implies an upper limit on the magnetic field strength at the event horizon of 26000 Gauss, which is consistent with previous measurements. The complex, low-brightness, double-sided jet structure in NGC 1052 makes it a challenge to detect the source at millimeter (mm) wavelengths. However, our first EHT observations have demonstrated that detection is possible up to at least 230 GHz. This study offers a glimpse through the dense surrounding torus and into the innermost central region, where the jets are formed. This has enabled us to finally resolve this region and provide improved constraints on its expansion and magnetic field strength., Comment: 22 pages, 10 figures, published in A&A

The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed: the main science targets, Sgr A* and M87 along with various calibrators. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-K\"onigl (BK) jet model. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2-86 GHz range. Then, we studied the dependences of the VLBI core flux density, size, and brightness temperature on the frequency measured in the AGN host frame. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Our results indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations.

We present a study of very low-mass stars and brown dwarfs in the rich star-forming core of the Rho Ophiuchi cloud complex. The selection of the sample relies on detecting the inherent water absorption characteristic in young substellar objects. Of the 22 water-bearing candidates selected, 15 have a spectral type of M6 or later. Brown dwarf candidates too faint for membership determination by Gaia have their proper motions derived by deep-infrared images spanning six years. Astrometric analysis confirms 21/22 sources as members, one identified as a contaminant. Infrared colors and the spectral energy distribution of each water-bearing candidate are used to diagnose the mass, age, and possible existence of circumstellar dust. 15 sources exhibit evidence of disks in their spectral energy distributions, as late as in M8-type objects. Spectroscopy for bright candidates has confirmed one as an M8 member and verified two sources (with disks) exhibiting signatures of magnetospheric accretion., Comment: 25 pages, 35 figures, 4 tables. Accepted for publication in The Astrophysical Journal

Using the Event Horizon Telescope (EHT), the gravitationally lensed rings around the supermassive black holes (SMBHs) in Messier 87 (M87) and Sagittarius A* (Sgr A*) have now been successfully imaged at a resolution under 10 gravitational radii (R$_{\rm g}$ $ = \rm{GM/c^2}$). To expand studies beyond M87 and Sgr A*, we have constructed the Event Horizon and Environs (ETHER) sample, a comprehensive database encompassing approximately 3.15 million SMBH mass estimates, $\sim$ 20,000 Very-Long Baseline Interferometry (VLBI) radio flux densities, and $\sim$ 36,000 hard X-ray flux densities. This database is designed to identify and optimize target selection for the EHT and its upgrades on the ground and in space. We have identified a Gold Sample (GS) of nearby low-luminosity Active Galactic Nuclei (AGNs) within it that are ideal for studying jet bases and potentially imaging black hole shadows. We observed 27 of these AGNs using the EHT from 2022 to 2024, providing an opportunity to resolve and image accretion flows and jets at resolutions of $\leq$ 100 R$_{\rm g}$. Only a few SMBHs have sufficiently high enough flux density to be imaged at scales of $\leq$ 50 R$_{\rm g}$ with the present EHT. Among these are M87, Sgr A*, NGC4594 (Sombrero/M104), NGC4261, and NGC4374 (Messier 84/M84). Of these, NGC4261, Sombrero, and M84 have been observed and/or are scheduled for deep imaging with EHT+ALMA from 2023 to 2025. Sombrero, NGC4261, M84, NGC4278, and NGC5232 are clearly detected in our EHT+ALMA observations in 2022, indicating that the 230 GHz flux density from the accretion flows is significantly high. Ongoing imaging of the ETHER GS will enable measurements of black hole mass and spin, help constrain General Relativity, and enrich our understanding of jet launching and accretion inflows across a broad multi-parameter space, including black hole mass, spin, accretion rate, and orientation., Comment: 9 pages, 6 figures, 1 table, published in Proceedings of the 16th EVN Symposium, Ed. E. Ros, P. Benke, S.A. Dzib, I. Rottmann, & J.A. Zensus, Bonn: Max-Planck-Institut f\"ur Radioastronomie, 2024, pages 75-84, https://cloud.mpifr-bonn.mpg.de/index.php/s/BkX2CC2Xjn2aKR4

As part of the BISTRO survey, we present JCMT 850 $\mu$m polarimetric observations towards the Orion Integral-Shaped Filament (ISF) that covers three portions known as OMC-1, OMC-2, and OMC-3. The magnetic field threading the ISF seen in the JCMT POL-2 map appears as a tale of three: pinched for OMC-1, twisted for OMC-2, and nearly uniform for OMC-3. A multi-scale analysis shows that the magnetic field structure in OMC-3 is very consistent at all the scales, whereas the field structure in OMC-2 shows no correlation across different scales. In OMC-1, the field retains its mean orientation from large to small scales, but shows some deviations at small scales. Histograms of relative orientations between the magnetic field and filaments reveal a bimodal distribution for OMC-1, a relatively random distribution for OMC-2, and a distribution with a predominant peak at 90$^\circ$ for OMC-3. Furthermore, the magnetic fields in OMC-1 and OMC-3 both appear to be aligned perpendicular to the fibers, which are denser structures within the filament, but the field in OMC-2 is aligned along with the fibers. All these suggest that gravity, turbulence, and magnetic field are each playing a leading role in OMC-1, 2, and 3, respectively. While OMC-2 and 3 have almost the same gas mass, density, and non-thermal velocity dispersion, there are on average younger and fewer young stellar objects in OMC-3, providing evidence that a stronger magnetic field will induce slower and less efficient star formation in molecular clouds., Comment: published in the ApJ Letters

Asteroseismology has emerged as a powerful tool to unravel the intricate relationships between evolved stars and their planetary systems. In this study, we leverage this technique to investigate the evolutionary stages of five exoplanet host stars, each exhibiting solar-like oscillations. Building on our previous work that identified two host stars as red clump and red giant branch (RGB) stars, this study focuses on a new and broader sample. By precisely measuring asteroseismic parameters such as the period spacing of dipole gravity modes ($\Delta\Pi_{1}$), we provide definitive confirmation of these stars' evolutionary states as subgiants or RGB stars. These results are not only crucial for understanding the internal structures of evolved stars but also for predicting the eventual fate of their planetary companions, which may face engulfment as their host stars expand. This research highlights the profound role of asteroseismology in advancing our knowledge of planetary system evolution and opens new pathways for exploring how stellar evolution impacts planetary survival. Our findings set the stage for future studies on the dynamic fates of exoplanets, providing key insights into the intricate processes of stellar and planetary evolution.

Eclipsing Algol-type systems containing a $\delta$ Scuti (hereafter $\delta$ Sct) star enable precise determination of physical parameters and the investigation of stellar internal structure and evolution. We present the absolute parameters of CZ Aquarius (hereafter CZ Aqr) based on TESS data. CZ Aqr has an orbital period of 0.86275209 d, a mass ratio of 0.489 (6), and the secondary component nearly fills its Roche lobe. $O-C$ analysis reveals a downward parabolic trend and a cyclical variation with a period of 88.2 yr. The downward parabola suggests a long-term decrease in the orbital period with $\dot{P}$ = -3.09$\times$$10^{-8}$ d $\textrm{yr}^{-1}$. The mass loss rate is estimated to be 4.54$\times$$10^{-9}$ M$_{\odot}$ $\textrm{yr}^{-1}$, which possibly due to magnetic stellar wind or hot spot. The cyclical variation might be caused by the light travel time effect via the presence of a third body with a minimum mass of $M_{3min}$ = 0.312 (21) M$_{\odot}$. Additionally, there are two possible celestial bodies in a 2:7 resonance orbit around CZ Aqr. The asymmetric light curve is explained by adding a hot spot on the surface of the primary star. After removing the binary model, 26 frequencies were extracted from TESS data. Two radial modes were newly identified among three possible independent frequencies. Our results show that the eclipsing Algol-type system is composed of a $\delta$ Sct primary star and a subgiant star in a quadruple system.

The massive semi-detached binary with a long-term decreasing orbital period may involve a rapid mass-transfer phase in Case A, and thus they are good astrophysical laboratories for investigating the evolution of massive binary stars. In this work, by using the long-term observational light curves from the OGLE project and other data in the low-metallicity LMC, four semi-detached massive binaries with long-term decreases in the orbital periods are detected from 165 EB-type close binaries. It is found that the more massive component in S07798 is filling its Roche lobe where the period decrease is caused by mass transfer from the primary to the secondary. However, the other three (S03065, S12631, S16873) are semi-detached binaries with a lobe-filling secondary where the mass transfer between the components should cause the period to increase if the angular momentum is conservative. The long-term period decreases in these three systems may be caused by the angular momentum loss. Additionally, the orbital periods of three systems (S03065, S07798, S16873) are detected to show cyclic variation with periods shorter than 11 years, which can be plausibly explained by the presence of close-in third bodies in these massive binaries. Based on all of these results, it is suggested that the detected four semi-detached binaries almost have multiplicity. The companion stars are crucial for the origin and evolution of these massive close binaries.

We present 850 $\mu$m polarization observations of the IC 348 star-forming region in the Perseus molecular cloud as part of the B-fields In STar-forming Region Observation (BISTRO) survey. We study the magnetic properties of two cores (HH 211 MMS and IC 348 MMS) and a filamentary structure of IC 348. We find that the overall field tends to be more perpendicular than parallel to the filamentary structure of the region. The polarization fraction decreases with intensity, and we estimate the trend by power-law and the mean of the Rice distribution fittings. The power indices for the cores are much smaller than 1, indicative of possible grain growth to micron size in the cores. We also measure the magnetic field strengths of the two cores and the filamentary area separately by applying the Davis-Chandrasekhar-Fermi method and its alternative version for compressed medium. The estimated mass-to-flux ratios are 0.45-2.20 and 0.63-2.76 for HH 211 MMS and IC 348 MMS, respectively, while the ratios for the filament is 0.33-1.50. This result may suggest that the transition from subcritical to supercritical conditions occurs at the core scale ($\sim$ 0.05 pc) in the region. In addition, we study the energy balance of the cores and find that the relative strength of turbulence to the magnetic field tends to be stronger for IC 348 MMS than HH 211 MMS. The result could potentially explain the different configurations inside the two cores: a single protostellar system in HH 211 MMS and multiple protostars in IC 348 MMS., Comment: Accepted for publication in ApJ. 21 pages, 12 figures

The first very long baseline interferometry (VLBI) detections at 870$\mu$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescopes in Chile, Hawaii, and Spain, obtained during observations in October 2018. The longest-baseline detections approach 11$\,$G$\lambda$ corresponding to an angular resolution, or fringe spacing, of 19$\mu$as. The Allan deviation of the visibility phase at 870$\mu$m is comparable to that at 1.3$\,$mm on the relevant integration time scales between 2 and 100$\,$s. The detections confirm that the sensitivity and signal chain stability of stations in the Event Horizon Telescope (EHT) array are suitable for VLBI observations at 870$\mu$m. Operation at this short wavelength, combined with anticipated enhancements of the EHT, will lead to a unique high angular resolution instrument for black hole studies, capable of resolving the event horizons of supermassive black holes in both space and time., Comment: Corresponding author: S. Doeleman

Heartbeat stars (HBSs) with tidally excited oscillations (TEOs) are ideal astrophysical laboratories for studying the internal properties of the systems. In this paper, five new HBSs exhibiting TEOs are discovered using TESS photometric data. The orbital parameters are derived using a corrected version of Kumar et al.'s model based on the Markov Chain Monte Carlo (MCMC) method. The TEOs in these objects are examined, and their pulsation phases and modes are identified. The pulsation phases of the TEOs in TIC 266809405, TIC 266894805, and TIC 412881444 are consistent with the dominant $l=2$, $m=0$, or $\pm2$ spherical harmonic. For TIC 11619404, although the TEO phase is close to the $m=+2$ mode, the $m = 0$ mode cannot be excluded because of the low inclination in this system. The TEO phase in TIC 447927324 shows a large deviation ($>2\sigma$) from the adiabatic expectations, suggesting that it is expected to be a traveling wave rather than a standing wave. In addition, these TEOs occur at relatively low orbital harmonics, and we cautiously suggest that this may be an observational bias. These objects are valuable sources for studying the structure and evolution of eccentricity orbit binaries and extending the TESS HBS catalog with TEOs., Comment: 9pages, 5figures

Heartbeat stars (HBSs) are ideal astrophysical laboratories to study the formation and evolution of binary stars in eccentric orbits and the internal structural changes of their components under strong tidal action. We discover 23 new HBSs based on TESS photometric data. The orbital parameters, including orbital period, eccentricity, orbital inclination, argument of periastron, and epoch of periastron passage of these HBSs are derived by using a corrected version of Kumar et al.'s model based on the Markov Chain Monte Carlo (MCMC) method. The preliminary results show that these HBSs have orbital periods in the range from 2.7 to 20 days and eccentricities in the range from 0.08 to 0.70. The eccentricity-period relation of these objects shows a positive correlation between eccentricity and period, and also shows the existence of orbital circularization. The Hertzsprung-Russell diagram shows that the HBSs are not all located in a particular area. The distribution of the derived parameters suggests a selection bias within the TESS survey towards massive HBSs with shorter orbital periods, higher temperatures and luminosities. These objects are a very useful source to study the structure and evolution of eccentricity orbit binaries and to extend the TESS HBS catalog., Comment: 10 pages, 7 figures

Variability can be the pathway to understanding the physical processes in astrophysical jets, however, the high-cadence observations required to test particle acceleration models are still missing. Here we report on the first attempt to produce continuous, >24 hour polarization light curves of blazars using telescopes distributed across the globe and the rotation of the Earth to avoid the rising Sun. Our campaign involved 16 telescopes in Asia, Europe, and North America. We observed BL Lacertae and CGRaBS J0211+1051 for a combined 685 telescope hours. We find large variations in the polarization degree and angle for both sources in sub-hour timescales as well as a ~180 degree rotation of the polarization angle in CGRaBS J0211+1051 in less than two days. We compared our high-cadence observations to Particle-In-Cell magnetic reconnection and turbulent plasma simulations. We find that although the state of the art simulation frameworks can produce a large fraction of the polarization properties, they do not account for the entirety of the observed polarization behavior in blazar jets., Comment: 20 pages, 15 figures, 2 tables, accepted for publication in A&A. The data used in the paper are available here: https://doi.org/10.7910/DVN/IETSXS

We reported a cyclic variation of $O-C$ diagram with a semi-amplitude of 0.0033 days and a period of 1.05 years for the pulsating eclipsing binary HZ Dra. The cyclic variation can be explained by the light travel-time effect via the presence of a close-in third body orbiting around HZ Dra in an elliptical orbit with a maximum semi-major axis of 0.92 au. Based on the W-D code, the contribution of the third light to the total system is determined to be 29 $\%$, which is in agreement with the estimated value. Our light curve modelling indicates an evolving hot and cool spot on the surface of the primary and secondary components, respectively. Their positions are roughly symmetrical to the inner Lagrangian point L1, which could be used to explain the variation in the O$^{'}$Connell effect. Our frequency analysis detects 1 radial p-mode, 7 non-radial p-modes and 1 non-radial g-mode. In addition, a total of 6 multiplets are identified, spaced by the orbital frequency, which can be explained as a tidally split mode caused by the equilibrium tides of the close binary system with a circular orbit. These pulsating features suggest that the primary of HZ Dra is a $\delta$ Scuti star, pulsating in both p- and g-mode and influenced by tidal forces., Comment: 15 pages, 11 figures. The first version is Under Review in ApJ

We model the TESS light curve of the binary system RX Dra, and also first calculate a lot of theoretical models to fit the g-mode frequencies previously detected from the TESS series of this system. The mass ratio is determined to be $q$=0.9026$^{+0.0032}_{-0.0032}$. We newly found that there are 16 frequencies (F1-F7, F11-F20) identified as dipole g-modes, two frequencies (F21, F22) identified as quadrupole g-modes, and another two frequencies (F23, F24) identified as g-sextupole modes, based on these model fits. The primary star is newly determined to be a $\gamma$ Dor pulsator in the main-sequence stage with a rotation period of about 5.7$^{+0.7}_{-0.2}$ days, rotating slower than the orbital motion. The fundamental parameters of two components are firstly estimated as follows: $M_1$=1.53$^{+0.00}_{-0.17}$ M $_{\odot}$, $M_2$= 1.38$^{+0.18}_{-0.00}$ M $_{\odot}$, $T_1$=7240$^{+490}_{-44}$ K, $T_2$=6747$^{+201}_{-221}$ K, $R_1$=1.8288$^{+0.0260}_{-0.0959}$ R $_{\odot}$, $R_2$= 1.3075$^{+0.0450}_{-0.2543}$ R $_{\odot}$, $L_1$=8.2830$^{+1.8015}_{-0.6036}$ L $_{\odot}$ and $L_2$= 3.4145$^{+0.1320}_{-0.1843}$ L $_{\odot}$. Our results show that the secondary star is in the solar-like pulsator region of the H-R diagram, indicating that it could be a pulsating star of this type. Finally, the radius of the convective core of the primary star is estimated to be about 0.1403$^{+0.0206}_{-0.0000}$ R$_{\odot}$., Comment: 13 pages,7 figures

Asteroseismology offers a profound window into stellar interiors and has emerged as a pivotal technique in exoplanetary research. This study harnesses the Transiting Exoplanet Survey Satellite (TESS) observations to reveal, for the first time, the asteroseismic oscillations of four exoplanet-hosting stars. Through meticulous analysis, we extracted their asteroseismic signatures, enabling the precise determination of stellar masses, radii, luminosities, and surface gravities. These parameters exhibit markedly reduced uncertainties compared to those derived from spectroscopic methods. Crucially, the exoplanets orbiting these stars were initially identified via radial velocity measurements. The refinement of host stellar masses necessitates a corresponding adjustment in planetary characteristics. Employing asteroseismology, we recalibrated the exoplanets' minimum masses and semi-major axes - a novel approach in the field. For instance, the exoplanet HD 5608 b's minimum mass, denoted as Msini, was ascertained to be 1.421 +/- 0.091 M_J through the integration of asteroseismic and radial velocity data. Similarly, two planets within the 7 CMa system yielded Msini values of 1.940 +/- 0.064 MJ and 0.912 +/- 0.067 M_J, respectively. Two planets in HD 33844 system presented Msini figures of 1.726 +/- 0.145 M_J and 1.541 +/- 0.182 M_J, while the HIP 67851 system's planets registered Msini at 1.243 +/- 0.139 M_J and a notably higher 5.387 +/- 0.699 M_J. This investigation extends beyond mere parameter refinement; it underscores the synergy between asteroseismology and exoplanetology, yielding unprecedented precision in system metrics. Focusing on quartet of K-type giants in advanced evolutionary phases, our work positions these systems as invaluable astrophysical laboratories, offering insights into the potential trajectory of our own solar system's fate.

We present the eclipsing Heartbeat Star KIC 9704906 with tidally excited oscillations (TEOs) and gamma Doradus/delta Scuti hybrid pulsations. The derived parameters show that it has an orbital period of $P$=8.7529108(1) days, a high eccentricity of $e$=0.467(3), and a high inclination of $i$=78$^{\circ}$.81(6). The mass ratio $q$=0.981(5), the relative radii (radius divided by semi-major axis) $r_1$=0.0639(2), and $r_2$=0.0715(4) indicate that the secondary component has a less mass and a larger radius, and may have evolved off the main sequence. The eight derived TEO candidates, $n$ = 3, 4, 5, 6, 7, 12, 40, and 44 harmonics, are consistent with or close to the dominant spherical harmonic $l=2$, $m=0$, or $\pm2$, assuming that the spin and orbital axes are aligned, and the pulsations are adiabatic and standing waves. We also identify ten independent frequency candidates, but one of them, $\mathit{f_{7}}$, is more like a modulation of a quasi-periodic signal and the orbits. According to the g-mode frequencies, we find that the rotation period of one component is 11.52(29) days. Although the masses and radii cannot be further constrained due to the lack of sufficient high-precision spectra, the fascinating phenomena in the Fourier spectra are evident and valuable in this system., Comment: 11 pages, 7 figures

V583 Lyr is an extremely low mass ratio Algol-type binary with an orbital period of 11.2580 days. We determined an effective temperature of T_{eff1} = 9000 \pm 350 K from newly observed spectra, which might be an underestimate due to binary mass transfer. The binary mass ratio q = 0.1 \pm 0.004 and the orbital inclination i = 85.5{\deg} are determined based on the assumption that the secondary fills its Roche lobe and rotates synchronously. The radial velocity curve is obtained from time series spectra, allowing for improved estimation of stellar masses and radii: M1 = 3.56 \pm 0.5 Msun, R1 = 2.4 \pm 0.2 Rsun; and M2 = 0.36 \pm 0.02 Msun, R2= 6.9 \pm 0.4 Rsun. The variations in the double-peaked H_{\alpha} emission indicate the formation of a stable disk during mass transfer. V583 Lyr appears to be a post-mass-reversal system, according to the estimated mass transfer using O-C period analysis. Its orbital period is slowly increasing, from which the rate of mass accretion by the primary star is estimated to be dM1/dt = 3.384 \times10^{-8} Msun/yr. The pulsation analysis was conducted on the residuals of the light curve. The primary component was found to be a g-mode pulsating star with 26 frequencies extracted lower than 9 d^{-1}. The frequency groups and rotational splitting properties of the g-mode were studied in detail. This study provides compelling evidence for an accretion disk surrounding the g-mode pulsating primary.

We report 850 $\mu$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity., Comment: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

Tidally excited oscillations (TEOs) in Heartbeat Stars (HBSs) are an essential probe of the internal properties of the systems, but their potential has yet to be fully exploited. Based on the orbital parameters of TEO candidates from our previous works, we identify the pulsation phases and amplitudes of TEOs in fourteen Kepler HBSs. Most pulsation phases of most systems can be explained by the dominant being $l=2$, $m=0$, or $\pm2$ spherical harmonic, assuming that the spin and orbital axes are aligned, and the pulsations are adiabatic and standing waves. The largest deviation ($>6\sigma$) occurs in KIC 8459354, which can be explained by the spin-orbit misalignment, and KIC 5877364 has a similar scenario. For KIC 11122789, almost half of the harmonics show large deviations; we cautiously suggest that these harmonics may not be considered TEO candidates. A similar scenario also exists in KIC 6290740. This phases and mode identification approach can also be used inversely to verify the TEO candidates derived by the Fourier analysis. Furthermore, the harmonics with large deviations ($>2\sigma$) in KIC 4377638, KIC 5090937, and KIC 11403032 can be expected to be travelling waves rather than standing waves. In addition, we also suggest that the apsidal motion could cause large deviations in TEO phases from theoretical values., Comment: published on MNRAS

TV Col is a long-period eclipsing intermediate polar (IPs) prototype star for the negative superhump (NSH) system. We investigate the eclipse minima, eclipse depth, and NSH amplitude based on TESS photometry. Using the same analytical method as SDSS J081256.85+191157.8, we find periodic variations of the O-C for eclipse minima and NSH amplitudes with periods of 3.939(25) d and 3.907(30) d, respectively. The periodic variation of the NSH amplitude of TV Col confirms that periodic NSH amplitude changes in response to the tilted disk precession may be universal, which is another evidence that the origin of the NSHs is related to the tilted disk precession. We suggest that the NSH amplitude variation may be similar to the superorbital signal, coming from the periodic change in visual brightness of the energy released by streams touching the tilted disk with tilted disk precession. Finally, we find for the first time that the eclipse depth exhibits bi-periodic variations with periods of P1 = 3.905(11) d and P2 = 1.953(4) d, respectively. P2 is about half of P1 and the disk precession period (P1~ Pprec ~ 2 * P2). We suggest that P1 may come from the periodic change in the brightness of the eclipse center due to tilted disk precession, while P2 may come from two accretion curtains precessing together with the tilted disk, but more verification and discussion are necessary. The discovery of bi-periodic variations in eclipse depth provides a new window for studying IPs and tilted disk precession., Comment: 18 pages, 7 figures, 3 tables; Accepted for publication in The Astrophysical Journal

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., Comment: 79 pages, 3 figures (main text) + 7 figures (extended data) + 2 figures (supplementary information). Published online in Nature on 15 November 2023

We uncover a hub-filament system associated with massive star formation in the G083.097$+$03.270. Diagnosed with simultaneous $^{12}$CO, $^{13}$CO, and C$^{18}$O line observations, the region is found to host two distinct and elongated filaments having separate velocity components, interacting spatially and kinematically, that appear to have seeded the formation of a dense hub at the intersection. A large velocity spread at the hub in addition to clear bridging feature connecting the filaments in velocity are indicating merging of filaments. Along the filaments axis, the velocity gradient reveals a global gas motion with an increasing velocity dispersion inward to the hub signifying turbulence. Altogether, the clustering of Class I sources, a high excitation temperature, a high column density, and presence of a massive outflow at the central hub suggest enhanced star formation. We propose that merging of large-scale filaments and velocity gradients along filaments are the driving factors in the mass accumulation process at the hub that have sequentially led to the massive star formation. With two giant filaments merging to coincide with a hub therein with ongoing star formation, this site serves as a benchmark for the `filaments to clusters' star-forming paradigm., Comment: 12 pages, 7 figures, 1 table, accepted for publication in ApJ

Negative superhumps (NSHs) are signals a few percent shorter than the orbital period of a binary star and are considered to originate from the reverse precession of the tilted disk. Based on TESS photometry, we find nine new cataclysmic variable stars (CVs) with NSHs. Three (ASAS J1420, TZ Per, and V392 Hya) of these stars similar to AH Her still have NSHs during dwarf nova outbursts, and the NSH amplitude varies with the outburst. The variation in the radius of the accretion disk partially explains this phenomenon. However, it does not explain the rebound of the NSH amplitude after the peak of the outburst and the fact that the NSH amplitude of the quiescence is sometimes not the largest, and it is necessary to combine the disk instability model (DIM) and add other ingredients. Therefore, we suggest that the variation of NSH amplitude with outburst can be an essential basis for studying the origin of NSHs and improving the DIM. The six ( ASASSN-V J1137, ASASSN-V J0611, 2MASS J0715, LAMOST J0925, ASASSN-17qj, and ZTF18acakuxo) remaining stars have been poorly studied, and for the first time we determine their orbital periods, NSHs and Superorbital signal (SOR) periods. The NSH periods and amplitudes of ASASSN-V J1137 and ASASSN-17qj vary with the SOR, and based on the comparison of the observations with the theory, we suggest that a single change in tilted disk angle does not explain the observations of the SOR and that other ingredients need to be considered as well., Comment: 21 pages, 12 figures, accepted for publication in the Astrophysical Journal

Heartbeat stars (HBSs) with tidally excited oscillations (TEOs) are ideal laboratories for studying the effect of equilibrium and dynamical tides. However, studies of TEOs in Kepler HBSs are rare due to the need for better modeling of the equilibrium tide in light curves. We revisit the HBSs in our previous work and study the TEOs in these HBSs based on the derived orbital parameters that could express the equilibrium tide. We also compile a set of analytic procedures to examine the harmonic and anharmonic TEOs in their Fourier spectra. The TEOs of 21 HBSs have been newly analyzed and presented. Twelve of these HBSs show prominent TEOs (signal-to-noise ratio of the harmonics $S/N \ge 10$). The relation between the orbital eccentricities and the harmonic number of the TEOs shows a positive correlation. The relation between the orbital periods and the harmonic number also shows a positive correlation. Furthermore, the distribution of HBSs with TEOs in the Hertzsprung-Russell (H-R) diagram shows that TEOs are more visible in hot stars with surface effective temperatures $T$ $\gtrsim$ 6500 K. These samples may also be valuable targets for future studies of the effect of tidal action in eccentric orbits., Comment: Published on ApJ

The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($\lambda$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region., Comment: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

This study conducts a theoretical investigation into the signal amplification arising from multiple Rabi sideband coherence within a one-dimensional waveguide quantum electrodynamics system. We utilize a semi-infinite waveguide to drive an anharmonic multi-level transmon with a strong coherent microwave field, examining the scattering behavior by introducing a probe signal. Our findings reveal signal amplification under specific resonant conditions, presenting spectra that reveal finer details than previously documented in the literature. To elucidate the mechanisms behind this amplification, we develop a model that explicitly accounts for multiple dressed sidebands in the presence of a strong driving field. From this model, we derive the reflection amplitude of the probe signal. Notably, our results indicate that amplification can occur due to either population inversion or, in some instances, through the constructive interference of multiple sidebands even in the absence of population inversion. Additionally, we explore how qubit dephasing impacts the amplification process.

Using membership of 85 open clusters from previous studies (Pang et al. 2021a,b, 2022b; Li et al. 2021) based on Gaia DR3 data, we identify binary candidates in the color-magnitude diagram, for systems with mass ratio q > 0.4. The binary fraction is corrected for incompleteness at different distances due to the Gaia angular resolution limit. We find a decreasing binary fraction with increasing cluster age, with substantial scatter. For clusters with a total mass > 200$M_\odot$, the binary fraction is independent of cluster mass. The binary fraction depends strongly on stellar density. Among four types of cluster environments, the lowest-density filamentary and fractal stellar groups have the highest mean binary fraction: 23.6% and 23.2%, respectively. The mean binary fraction in tidal-tail clusters is 20.8%, and is lowest in the densest halo-type clusters: 14.8%. We find clear evidence of early disruptions of binary stars in the cluster sample. The radial binary fraction depends strongly on the cluster-centric distance across all four types of environments, with the smallest binary fraction within the half-mass radius $r_h$, and increasing towards a few $r_h$. Only hints of mass segregation is found in the target clusters. The observed amount of mass segregation is not significant to generate a global effect inside the target clusters. We evaluate the bias of unresolved binary systems (assuming a primary mass of 1$M_\odot$) in 1D tangential velocity, which is 0.1-1$\,\rm km\,s^{-1}$. Further studies are required to characterize the internal star cluster kinematics using Gaia proper motions.

We present observations of polarized dust emission at 850 $\mu$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $\mu$G in the main starless core and up to $\sim$90$\pm$40 $\mu$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfv\'enic. We also present a new method of data reduction for these denser but fainter objects like starless cores., Comment: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

The nearby radio galaxy M87 is a prime target for studying black hole accretion and jet formation^{1,2}. Event Horizon Telescope observations of M87 in 2017, at a wavelength of 1.3 mm, revealed a ring-like structure, which was interpreted as gravitationally lensed emission around a central black hole^3. Here we report images of M87 obtained in 2018, at a wavelength of 3.5 mm, showing that the compact radio core is spatially resolved. High-resolution imaging shows a ring-like structure of 8.4_{-1.1}^{+0.5} Schwarzschild radii in diameter, approximately 50% larger than that seen at 1.3 mm. The outer edge at 3.5 mm is also larger than that at 1.3 mm. This larger and thicker ring indicates a substantial contribution from the accretion flow with absorption effects in addition to the gravitationally lensed ring-like emission. The images show that the edge-brightened jet connects to the accretion flow of the black hole. Close to the black hole, the emission profile of the jet-launching region is wider than the expected profile of a black-hole-driven jet, suggesting the possible presence of a wind associated with the accretion flow., Comment: 50 pages, 18 figures, 3 tables, author's version of the paper published in Nature