Search

Stars in close binaries are tidally distorted, and this has a strong effect on their pulsation modes. We compute the mode frequencies and geometries of tidally distorted stars using perturbation theory, accounting for the effects of the Coriolis force and the coupling between different azimuthal orders $m$ of a multiplet induced by the tidal distortion. For tidally coupled dipole pressure modes, the tidal coupling dominates over the Coriolis force and the resulting pulsations are ``triaxial", with each of the three modes in a multiplet ``tidally tilted" to be aligned with the one of the three principal axes of the star. The observed amplitudes and phases of the dipole modes aligned orthogonal to the spin axis are modulated throughout the orbit, producing doublets in the power spectrum that are spaced by exactly twice the orbital frequency. Quadrupole modes have similar but slightly more complex behavior. This amplitude modulation allows for mode identification which can potentially enable detailed asteroseismic analyses of tidally tilted pulsators. Pressure modes should exhibit this behavior in stellar binaries close enough to be tidally synchronized, while gravity modes should remain aligned with the star's spin axis. We discuss applications to various types of pulsating stars, and the relationship between tidal tilting of pulsations and the ``single-sided" pulsations sometimes observed in very tidally distorted stars., Comment: Submitted to ApJ, comments welcome!

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 origin of magnetic fields in white dwarfs (WDs) remains mysterious. Magnetic WDs are traditionally associated with field strengths $\gtrsim1\,\mathrm{MG}$, set by the sensitivity of typical spectroscopic magnetic field measurements. Informed by recent developments in red giant magnetoasteroseismology, we revisit the use of WD pulsations as a seismic magnetometer. WD pulsations primarily probe near-surface magnetic fields, whose effect on oscillation mode frequencies is to asymmetrize rotational multiplets and, if strong enough, suppress gravity-mode propagation altogether. The sensitivity of seismology to magnetic fields increases strongly with mode period and decreases quickly with the depth of the partial ionization-driven surface convective zone. We place upper limits for magnetic fields in $24$ pulsating WDs: $20$ hydrogen-atmosphere (DAV) and three helium-atmosphere (DBV) carbon-oxygen WDs, and one extremely low-mass (helium-core) pulsator. These bounds are typically $\sim1$-$10\,\mathrm{kG}$, although they can reach down to $\sim10$-$100\,\mathrm{G}$ for DAVs and helium-core WDs in which lower-frequency modes are excited. Seismic magnetometry may enable new insights into the formation and evolution of WD magnetism., Comment: 21 pages, 9 figures, submitted to ApJ

The Dice Similarity Coefficient (DSC) is the current de facto standard to determine agreement between a reference segmentation and one generated by manual / auto-contouring approaches. This metric is useful for non-spatially important images; however, radiation therapy requires consideration of nearby Organs-at-Risk (OARs) and their radiosensitivity which are currently unaccounted for with the traditional DSC. In this work, we introduce the OAR-DSC which accounts for nearby OARs and their radiosensitivity when computing the DSC. We illustrate the importance of this through cases where two proposed contours have similar DSC, but lower OAR-DSC when one contour expands closer to the surrounding OARs. This work is important because the OAR-DSC may be used by deep learning auto-contouring algorithms in a radiation therapy specific loss function, thereby progressing on the current disregard for the importance of these differences on the final radiation dose plan generation, delivery, and risks of patient toxicity., Comment: Proceedings of the XX-th international conference on the use of computers in radiation therapy (ICCR)

CMOS-based microelectronics are limited to ~150{\deg}C and therefore not suitable for the extreme high temperatures in aerospace, energy, and space applications. While wide bandgap semiconductors can provide high-temperature logic, nonvolatile memory devices at high temperatures have been challenging. In this work, we develop a nonvolatile electrochemical memory cell that stores and retains analog and digital information at temperatures as high as 600 {\deg}C. Through correlative electron microscopy, we show that this high-temperature information retention is a result of composition phase separation between the oxidized and reduced forms of amorphous tantalum oxide. This result demonstrates a memory concept that is resilient at extreme temperatures and reveals phase separation as the principal mechanism that enables nonvolatile information storage in these electrochemical memory cells., Comment: 19 pages, 4 figures

The recently discovered Pa 30 nebula, the putative type Iax supernova remnant associated with the historical supernova of 1181 AD, shows puzzling characteristics that make it unique among known supernova remnants. In particular, Pa 30 exhibits a complex morphology, with a unique radial and filamentary structure, and it hosts a hot stellar remnant at its center, which displays oxygen-dominated, ultra-fast winds. Because of the surviving stellar remnant and the lack of hydrogen and helium in its filaments, it has been suggested that Pa 30 is the product of a failed thermonuclear explosion in a near- or super-Chandrasekhar white dwarf, which created a sub-luminous transient, a rare sub-type of the Ia class of supernovae called type Iax. We here present a detailed study of the 3D structure and velocities of a full radial section of the remnant. The Integral Field Unit (IFU) observations, obtained with the new red channel of the Keck Cosmic Web Imager spectrograph, reveal that the ejecta are consistent with being ballistic, with velocities close to the free-expansion velocity. Additionally, we detect a large cavity inside the supernova remnant and a sharp inner edge to the filamentary structure, which coincides with the outer edge of a bright ring detected in infrared images. Finally, we detect a strong asymmetry in the amount of ejecta along the line of sight, which might hint to an asymmetric explosion. Our analysis provides strong confirmation that the explosion originated from SN 1181., Comment: 13 pages, 8 figures. Accepted for publication in ApJL

Failed supernovae (SNe), which are likely the main channel for forming stellar-mass black holes, are predicted to accompany mass ejections much weaker than typical core-collapse SNe. We conduct a grid of one-dimensional radiation hydrodynamical simulations to explore the emission of failed SNe from red supergiant progenitors, leveraging recent understanding of the weak explosion and the dense circumstellar matter (CSM) surrounding these stars. We find from these simulations and semi-analytical modeling that diffusion in the CSM prolongs the early emission powered by shock breakout/cooling. The early emission has peak luminosities of $\sim 10^7$-$10^8~L_\odot$ in optical and UV, and durations of days to weeks. The presence of dense CSM aids detection of the early bright peak from these events via near-future wide-field surveys such as Rubin Observatory, ULTRASAT and UVEX., Comment: 17 pages, 8 figures, submitted to ApJ. Model light curves available here: https://github.com/DTsuna/failedSNeLCs

We show that strong interactions violating the conservation of lepton numbers in the neutrino sector could significantly alter the standard low entropy picture for the pre-supernova collapsing core of a massive star. A rapid neutrino-antineutrino equilibration leads to entropy generation and enhanced electron capture and, hence, a lower electron fraction than in the standard model. This would affect the downstream core evolution, the prospects for a supernova explosion, and the emergent neutrino signal. If realized by lepton-number-violating neutrino self-interactions (LNV $\nu$SI), the relevant mediator mass and coupling ranges can be probed by future accelerator-based experiments., Comment: 5+8 pages, 2+5 figures. Comments are welcome!

We analyze JWST MIRI/MRS IFU observations of three Seyferts and showcase the intriguing polycyclic aromatic hydrocarbon (PAH) emission characteristics in regions of $\sim 500\,\rm pc$ scales over or around their active galactic nuclei (AGN). Combining the model predictions and the measurements of PAH features and other infrared emission lines, we find that the central regions containing a high fraction of neutral PAHs with small sizes, e.g., those in ESO137-G034, are in highly heated environments, due to collisional shock heating, with hard and moderately intense radiation fields. Such environments are proposed to be associated with inhibited growth or preferential erosion of PAHs, decreasing the average PAH size and the overall abundance of PAHs. We additionally find that the central regions containing a high fraction of ionized PAHs with large sizes, e.g., those in MCG-05-23-016, are likely experiencing severe photo-ionization because of the radiative effects from the radiative shock precursor besides the AGN. The severe photo-ionization can contribute to the ionization of all PAHs and further destruction of small PAHs. Overall, different Seyferts, even different regions in the same galaxy, e.g., those in NGC\,3081, can contain PAH populations of different properties. Specifically, Seyferts that exhibit similar PAH characteristics to ESO137-G034 and MCG-05-23-016 also tend to have similar emission line properties to them, suggesting that the explanations for PAH characteristics of ESO137-G034 and MCG-05-23-016 may also apply generally. These results have promising application in the era of JWST, especially in diagnosing different (i.e., radiative, and kinetic) AGN feedback modes., Comment: ApJL accepted on September 26th, title slightly modified in accordance with ApJL standards, Fig. 2 updated with additional labels

Utilizing JWST MIRI/MRS IFU observations of the kiloparsec scale central regions, we showcase the diversity of ionized gas distributions and kinematics in six nearby Seyfert galaxies included in the GATOS survey. Specifically, we present spatially resolved flux distribution and velocity field maps of six ionized emission lines covering a large range of ionization potentials ($15.8-97.1$ eV). Based on these maps, we showcase the evidence of ionized gas outflows in the six targets, and find some highly disturbed regions in NGC\,5728, NGC\,5506, and ESO137-G034. We propose AGN-driven radio jets plausibly play an important role in triggering these highly disturbed regions. With the outflow rates estimated based on [Ne~{\footnotesize V}] emission, we find the six targets tend to have ionized outflow rates converged to a narrower range than previous finding. These results have important implication for the outflow properties in AGN of comparable luminosity., Comment: 34 pages (11 pages in the appendix), 18 figures in the main text, ApJ in press (accepted on July 26th)

Since the invention of the laser nonlinear effects such as filamentation, Rabi-cycling and collective emission have been explored in the optical regime leading to a wide range of scientific and industrial applications. X-ray free electron lasers (XFELs) have led to the extension of many optical techniques to X-rays for their advantages of angstrom scale spatial resolution and elemental specificity. One such example is XFEL driven population inversion of 1s core hole states resulting in inner-shell K${\alpha}$ (2p to 1s) X-ray lasing in elements ranging from neon to copper, which has been utilized for nonlinear spectroscopy and development of next generation X-ray laser sources. Here we show that strong lasing effects, similar to those observed in the optical regime, can occur at 1.5 to 2.1 angstrom wavelengths during high intensity (> ${10^{19}}$ W/cm${^{2}}$) XFEL driven inner-shell lasing and superfluorescence of copper and manganese. Depending on the temporal substructure of the XFEL pump pulses, the resulting inner-shell X-ray laser pulses can exhibit strong spatial inhomogeneities as well as spectral inhomogeneities and broadening. Through 3D Maxwell Bloch theory we show that the observed spatial inhomogeneities result from X-ray filamentation, and that the spectral broadening is driven by Rabi cycling with sub-femtosecond periods. These findings indicate that we have generated Angstrom-wavelength x-ray pulses (containing ${10^{6}}$ - ${10^{8}}$ photons) in the strong lasing regime, some of them with pulse lengths of less than 100 attoseconds.

The main treatment modality for oropharyngeal cancer (OPC) is radiotherapy, where accurate segmentation of the primary gross tumor volume (GTVp) is essential. However, accurate GTVp segmentation is challenging due to significant interobserver variability and the time-consuming nature of manual annotation, while fully automated methods can occasionally fail. An interactive deep learning (DL) model offers the advantage of automatic high-performance segmentation with the flexibility for user correction when necessary. In this study, we examine interactive DL for GTVp segmentation in OPC. We implement state-of-the-art algorithms and propose a novel two-stage Interactive Click Refinement (2S-ICR) framework. Using the 2021 HEad and neCK TumOR (HECKTOR) dataset for development and an external dataset from The University of Texas MD Anderson Cancer Center for evaluation, the 2S-ICR framework achieves a Dice similarity coefficient of 0.713 $\pm$ 0.152 without user interaction and 0.824 $\pm$ 0.099 after five interactions, outperforming existing methods in both cases.

We analyze JWST MIRI/MRS observations of the infrared PAH bands in the nuclear and circumnuclear regions of local AGN from the GATOS Survey. In this work, we examine the PAH properties in the circumnuclear regions of AGN and AGN-outflows, and compare them to those in star-forming regions and the innermost regions of AGN. This study employs 4.9-28.1 micron sub-arcsecond angular resolution data to investigate the properties of PAH in three nearby sources (DL~30-40 Mpc). Our findings align with previous JWST studies, showing that the central regions of AGN show a larger fraction of neutral PAH molecules (i.e. elevated 11.3/6.2 and 11.3/7.7 PAH ratios) compared to star-forming galaxies. We find that the AGN might affect not only the PAH population in the innermost region but also in the extended regions up to ~kpc scales. By comparing our observations to PAH diagnostic diagrams, we find that, in general, regions located in the projected direction of the AGN-outflow occupy similar positions on the PAH diagnostic diagrams as those of the innermost regions of AGN. Star-forming regions that are not affected by the AGN in these galaxies share the same part of the diagram as Star-forming galaxies. We examine the potential of the PAH-H2 diagram to disentangle AGN versus star-forming activity. Our results suggest that in Sy-like AGN, illumination and feedback from the AGN might affect the PAH population at nuclear and kpc scales, in particular, the ionization state of the PAH grains. However, PAH sizes are rather similar. The carriers of the ionized PAH bands (6.2 and 7.7 micron) are less resilience than those of neutral PAH bands (11.3 micron), which might be particularly important for strongly AGN-host coupled systems. Therefore, caution must be applied when using PAH bands as star-formation rate indicators in these systems even at kpc scales, with the ionized ones being more affected by the AGN., Comment: Accepted for publication in A&A. 21 pages, 13 Figures

Filaments are believed to play a key role in high-mass star formation. We present a systematic study of the filaments and their hosting clumps in the G35 molecular complex using JCMT SCUBA-2 850 $\micron$ continuum data. We identified five clouds in the complex and 91 filaments within them, some of which form 10 hub-filament systems (HFSs), each with at least 3 hub-composing filaments. We also compiled a catalogue of 350 dense clumps, 183 of which are associated with the filaments. We investigated the physical properties of the filaments and clumps, such as mass, density, and size, and their relation to star formation. We find that the global mass-length trend of the filaments is consistent with a turbulent origin, while the hub-composing filaments of high line masses ($m_{\rm l}\,>$\,230\,$\mathrm{M_{\odot}~pc^{-1}}$) in HFSs deviate from this relation, possibly due to feedback from massive star formation. We also find that the most massive and densest clumps (R\,$>$\,0.2\,pc, M\,$>35\,\mathrm{M_{\odot}}$, $\mathrm{\Sigma}>\,0.05\,\mathrm{g~cm^{-2}}$) are located in the filaments and in the hubs of HFS with the latter bearing a higher probability of occurrence of high-mass star-forming signatures, highlighting the preferential sites of HFSs for high-mass star formation. We do not find significant variation in the clump mass surface density across different evolutionary environments of the clouds, which may reflect the balance between mass accretion and stellar feedback., Comment: 34 pages, 17 figures. Accepted for publication in ApJ

The Transiting Exoplanet Survey Satellite (TESS) has enabled the discovery of numerous tidally tilted pulsators (TTPs), which are pulsating stars in close binaries where the presence of a tidal bulge has the effect of tilting the primary star's pulsation axes into the orbital plane. Recently, the modeling framework developed to analyze TTPs has been applied to the emerging class of tri-axial pulsators, which exhibit nonradial pulsations about three perpendicular axes. In this work, we report on the identification of the second-ever discovered tri-axial pulsator, with sixteen robustly-detected pulsation multiplets, of which fourteen are dipole doublets separated by 2$\nu_{\rm orb}$. We jointly fit the spectral energy distribution (SED) and TESS light curve of the star, and find that the primary is slightly evolved off the zero-age main sequence, while the less massive secondary still lies on the zero-age main sequence. Of the fourteen doublets, we associate eight with $Y_{10x}$ modes and six with novel $Y_{10y}$ modes. We exclude the existence of $Y_{11x}$ modes in this star and show that the observed pulsation modes must be $Y_{10y}$. We also present a toy model for the tri-axial pulsation framework in the context of this star. The techniques presented here can be utilized to rapidly analyze and confirm future tri-axial pulsator candidates., Comment: 19 pages, 11 figures, 4 tables. Accepted for publication in ApJ

Understanding how stars form, evolve and impact molecular clouds is key to understanding why star formation is such an inefficient process globally. In this paper, we use the infrared bright fraction, $f_\text{IRB}$ (the fraction of a given molecular cloud that appears bright against the 8 $\mu$m Milky Way background) as a proxy for time evolution to test how cloud properties change as star formation evolves. We apply this metric to 12,000 high-mass star-forming molecular clouds we identify using the Herschel-Hi-GAL survey between $|l|<70^\circ$ on the Milky Way plane. We find clouds are not static while forming stars. Instead, molecular clouds continuously gain mass while star formation progresses. By performing principal component analysis on the cloud properties, we find that they evolve down two paths distinguished by their mass gain. Most clouds (80%) gain four times more mass as a function of $f_\text{IRB}$. The remaining 20% experience an extreme period of growth, growing in mass by a factor of 150 on average and during this period, they initially gain mass fast enough to outpace their star formation. For all clouds, it is only after half their area becomes star forming that mass loss occurs. We expect stellar feedback and potentially galactic shear is responsible. By analysing cloud positions, we suggest that the rate of mass growth may be linked to the larger galactic environment. Altogether, these results have strong implications on how we assess star forming ability on cloud scales when assuming molecular cloud masses are fixed in time., Comment: 20 pages, 24 figures and 4 tables. Submitted to MNRAS

Hot subdwarf B (sdB) stars are stripped helium-burning stars that are often found in close binaries, where they experience strong tidal interactions. The dissipation of tidally excited gravity waves alters their rotational evolution throughout the sdB lifetime. While many sdB binaries have well-measured rotational and orbital frequencies, there have been few theoretical efforts to accurately calculate the tidal torque produced by gravity waves. In this work, we directly calculate the tidal excitation of internal gravity waves in realistic sdB stellar models and integrate the coupled spin-orbit evolution of sdB binaries. We find that for canonical sdB ($M_\mathrm{sdB}=0.47\,M_\odot$) binaries, the transitional orbital period below which they could reach tidal synchronization in the sdB lifetime is $\sim \! 0.2\;\mathrm{days}$, with weak dependence on the companion masses. For low-mass sdBs ($M_\mathrm{sdB}=0.37\,M_\odot$) formed from more massive progenitor stars, the transitional orbital period becomes $\sim \! 0.15\;\mathrm{days}$. These values are very similar to the tidal synchronization boundary ($\sim \! 0.2\;\mathrm{days}$) evident from observations. We discuss the dependence of tidal torques on stellar radii, and we make predictions for the rapidly rotating white dwarfs formed from synchronized sdB binaries., Comment: Accepted to ApJ

In this Roadmap, we present a vision for the future of submillimetre and millimetre astronomy in the United Kingdom over the next decade and beyond. This Roadmap has been developed in response to the recommendation of the Astronomy Advisory Panel (AAP) of the STFC in the AAP Astronomy Roadmap 2022. In order to develop our stragetic priorities and recommendations, we surveyed the UK submillimetre and millimetre community to determine their key priorities for both the near-term and long-term future of the field. We further performed detailed reviews of UK leadership in submillimetre/millimetre science and instrumentation. Our key strategic priorities are as follows: 1. The UK must be a key partner in the forthcoming AtLAST telescope, for which it is essential that the UK remains a key partner in the JCMT in the intermediate term. 2. The UK must maintain, and if possible enhance, access to ALMA and aim to lead parts of instrument development for ALMA2040. Our strategic priorities complement one another: AtLAST (a 50m single-dish telescope) and an upgraded ALMA (a large configurable interferometric array) would be in synergy, not competition, with one another. Both have identified and are working towards the same overarching science goals, and both are required in order to fully address these goals., Comment: 91 pages plus cover, 38 figures. Submitted to the Science and Technology Facilities Council, August 2024. One figure corrected (v2); new appendix with STFC Q&A; corrected SMA access statement; updated references, acronyms & author list (v3)

We use data from the ALMA Evolutionary Study of High Mass Protocluster Formation in the Galaxy (ALMAGAL) survey to study 100 ALMAGAL regions at $\sim$ 1 arsecond resolution located between $\sim$ 2 and 6 kpc distance. Using ALMAGAL $\sim$ 1.3mm line and continuum data we estimate flow rates onto individual cores. We focus specifically on flow rates along filamentary structures associated with these cores. Our primary analysis is centered around position velocity cuts in H$_2$CO (3$_{0,3}$ - 2$_{0,2}$) which allow us to measure the velocity fields, surrounding these cores. Combining this work with column density estimates we derive the flow rates along the extended filamentary structures associated with cores in these regions. We select a sample of 100 ALMAGAL regions covering four evolutionary stages from quiescent to protostellar, Young Stellar Objects (YSOs), and HII regions (25 each). Using dendrogram and line analysis, we identify a final sample of 182 cores in 87 regions. In this paper, we present 728 flow rates for our sample (4 per core), analysed in the context of evolutionary stage, distance from the core, and core mass. On average, for the whole sample, we derive flow rates on the order of $\sim$10$^{-4}$ M$_{sun}$yr$^{-1}$ with estimated uncertainties of $\pm$50%. We see increasing differences in the values among evolutionary stages, most notably between the less evolved (quiescent/protostellar) and more evolved (YSO/HII region) sources. We also see an increasing trend as we move further away from the centre of these cores. We also find a clear relationship between the flow rates and core masses $\sim$M$^{2/3}$ which is in line with the result expected from the tidal-lobe accretion mechanism. Overall, we see increasing trends in the relationships between the flow rate and the three investigated parameters; evolutionary stage, distance from the core, and core mass., Comment: 11 pages, 11 figures, accepted for publication in A&A

The intrinsic frustrated nature of a kagome lattice is amenable to the realization of exotic phases of matter, such as quantum spin liquids or spin ices, and more recently the multiple-$\mathrm{\textbf{q}}$ charge density waves (CDW) in the kagome metals. Despite intense efforts to understand the mechanism driving the electronic modulations, its origin is still unknown and hindered by competing interactions and intertwined orders. Here, we identify a dimerization-driven 2D hexagonal charge-diffuse precursor in the antiferromagnetic kagome metal FeGe and demonstrate that the fraction of dimerized/undimerized states is the relevant order parameter of the multiple-$\mathrm{\textbf{q}}$ CDW of a continuous phase transition. The pretransitional charge fluctuations with propagation vector $\mathrm{\textbf{q}=\textbf{q}_M}$ at T$_{\mathrm{CDW}}$$<$T$<$T$^*$(125 K) are anisotropic, hence holding a quasi-long-range bond-orientational order. The broken translational symmetry emerges from the anisotropic diffuse precursor, akin to the Ising scenario of antiferromagnetic triangular lattices. The temperature and momentum dependence of the critical scattering show parallels to the stacked hexatic $\mathrm{B}$-phases reported in liquid crystals and transient states of CDWs and highlight the key role of the topological defect-mediated melting of the CDW in FeGe.

Tensor networks and quantum computation are two of the most powerful tools for the simulation of quantum many-body systems. Rather than viewing them as competing approaches, here we consider how these two methods can work in tandem. We introduce a novel algorithm that combines tensor networks and quantum computation to produce results that are more accurate than what could be achieved by either method used in isolation. Our algorithm is based on multiproduct formulas (MPF) - a technique that linearly combines Trotter product formulas to reduce algorithmic error. Our algorithm uses a quantum computer to calculate the expectation values and tensor networks to calculate the coefficients used in the linear combination. We present a detailed error analysis of the algorithm and demonstrate the full workflow on a one-dimensional quantum simulation problem on $50$ qubits using two IBM quantum computers: $ibm\_torino$ and $ibm\_kyiv$.

We present observations of the type-2 Seyfert NGC7172 obtained with the medium-resolution spectrometer (MRS) of the Mid-Infrared Instrument (MIRI) on board of the JWST. This galaxy hosts one of the lowest ionised gas mass outflow rates (Mout~0.005 M/yr) in a sample of six AGN with similar bolometric luminosities (log Lbol~44erg/s) within the Galactic Activity, Torus and Outflow Survey (GATOS). We aim to understand the properties of the ionised gas outflow, mainly using the emission lines from the neon transitions, that cover a broad range of ionisation potentials (IP) from ~20 eV to ~130 eV. We applied parametric and non-parametric methods to characterise the line emission and kinematics. The low excitation lines (IP<25eV, e.g.[NeII]) trace the rotating disc emission. The high excitation lines (IP>90eV, e.g.[NeV]), which are likely photoionised exclusively by the AGN, are expanding in the direction nearly perpendicular to the disc of the galaxy, with maximum projected velocities of ~350-500 km/s. In particular, [NeV] and [NeVI] lines reveal a biconical ionised gas outflow emerging N-S from the nuclear region, extending at least ~2.5"N and 3.8"S (projected distance of ~450 and 680 pc). Most of the emission arising in the northern part of the cone was not previously detected due to obscuration. Given the almost face-on orientation of the outflow and the almost edge-on orientation of the galaxy, NGC7172 may be a case of weak coupling. Nevertheless, we found evidence for positive feedback in two distinct outflowing clumps at projected distances of 3.1" and 4.3" (i.e. ~560 and 780 pc) SW from the AGN. We estimated a star formation rate in these regions using the [NeII] and [NeIII] luminosities of 0.08 M/yr, that is ~10% of that found in the circumnuclear ring. The star formation activity might have been triggered by the interaction between the ionised gas outflow and the ISM of the galaxy., Comment: 20 pages, 6 figures. Accepted for publication in A&A. Abstract adapted for the arxiv version

The estimation of low energies of many-body systems is a cornerstone of computational quantum sciences. Variational quantum algorithms can be used to prepare ground states on pre-fault-tolerant quantum processors, but their lack of convergence guarantees and impractical number of cost function estimations prevent systematic scaling of experiments to large systems. Alternatives to variational approaches are needed for large-scale experiments on pre-fault-tolerant devices. Here, we use a superconducting quantum processor to compute eigenenergies of quantum many-body systems on two-dimensional lattices of up to 56 sites, using the Krylov quantum diagonalization algorithm, an analog of the well-known classical diagonalization technique. We construct subspaces of the many-body Hilbert space using Trotterized unitary evolutions executed on the quantum processor, and classically diagonalize many-body interacting Hamiltonians within those subspaces. These experiments show that quantum diagonalization algorithms are poised to complement their classical counterpart at the foundation of computational methods for quantum systems., Comment: 25 pages, 13 figures

The expanding molecular ring (EMR) manifests itself as a parallelogram in the position-velocity diagram of spectral line emission from the Central Molecular Zone (CMZ) surrounding the Galacic centre (GC). Using multiwavelength data, we investigate the gas kinematics, star formation activity, and the presence of shocked gas in a 200 pc long high velocity gas stream (V~ +150 km/s) with a double helix morphology named the helix stream, that is located 15-55 pc above the CMZ and is kinematically associated with the EMR/parallelogram. We carried out molecular line observations using the IRAM 30m, Yebes 40m, and APEX 12m telescopes. The detection of four rotational transitions of the SiO molecule indicate the presence of shocks. We derived the SiO column densities and abundances in different regions of the helix stream. The presence of protostellar clumps and a candidate HII region signify the ongoing star formation activity within the helix stream. The cloud is massive (2.5x10^6 M_sun) and highly turbulent. We find evidence of cloud-cloud collisions towards the eastern edge (l~1.3{\deg}), suggesting a dynamic interaction with the CMZ. An expanding shell is detected within the cloud with radius of 6.7 pc and an expansion velocity of 35 km/s. The shell might be powered by several supernovae or a single hypernova. The SiO abundance within the helix stream implies extensive shock processes occurring on large scales. The helical or cork-screw velocity structure of the helix stream indicates twisting and turning motions within the cloud. We propose that the helix stream is the continuation of the near side bar lane, that is overshooting after brushing the CMZ. Our findings carry profound implications for understanding star formation in extreme conditions and elucidate the intricate properties of gas and dust associated with nuclear inflows in barred spiral galaxies., Comment: 20 pages, 20 figures, 4 tables, accepted for publication in A&A

Digital twin models are of high interest to Head and Neck Cancer (HNC) oncologists, who have to navigate a series of complex treatment decisions that weigh the efficacy of tumor control against toxicity and mortality risks. Evaluating individual risk profiles necessitates a deeper understanding of the interplay between different factors such as patient health, spatial tumor location and spread, and risk of subsequent toxicities that can not be adequately captured through simple heuristics. To support clinicians in better understanding tradeoffs when deciding on treatment courses, we developed DITTO, a digital-twin and visual computing system that allows clinicians to analyze detailed risk profiles for each patient, and decide on a treatment plan. DITTO relies on a sequential Deep Reinforcement Learning digital twin (DT) to deliver personalized risk of both long-term and short-term disease outcome and toxicity risk for HNC patients. Based on a participatory collaborative design alongside oncologists, we also implement several visual explainability methods to promote clinical trust and encourage healthy skepticism when using our system. We evaluate the efficacy of DITTO through quantitative evaluation of performance and case studies with qualitative feedback. Finally, we discuss design lessons for developing clinical visual XAI applications for clinical end users.

We developed a state-space age-structured catch-at-length (ACL) assessment model for redfish in NAFO Divisions 3LN. The model was developed to address limitations in the surplus production model that was previously used to assess this stock. The ACL model included temporal variations in recruitment, growth, and mortality rates, which were limitations identified for the surplus production model. Our ACL model revealed some important discrepancies in survey and fishery length compositions. Our model also required large population dynamics process errors to achieve good fits to survey indices and catch estimates, which also demonstrated that additional understanding of these data and other model assumptions is required. As such, we do not propose the ACL model to provide management advice for 3LN redfish, but we do provide research recommendations that should provide a better basis to model the 3LN redfish stock dynamics. Recommendations include implementing sampling programs to determine redfish species/ecotypes in commercial and research survey catches and improving biological sampling for maturity and age., Comment: 27 pages including references, tables, and figures. In addition 12 pages figures in an Appendix

This paper presents a study of the global medium-scale (scales<620 km) gravity wave (GW) activity (in terms of zonal wind variance) and its seasonal, local time and longitudinal variations by employing the enhanced-resolution (~50 km) Whole Atmosphere Model (WAMT254) and space-based observations for geomagnetically quiet conditions. It is found that the GW hotspots produced by WAMT254 in the troposphere and stratosphere agree well with previously well-studied orographic and non-orographic sources. In the ionosphere-thermosphere (IT) region, GWs spread out forming latitudinal band-like hotspots. During solstices, a primary maximum in GW activity is observed in WAMT254 and GOCE over winter mid-high latitudes, likely associated with higher-order waves with primary sources in polar night jet, fronts and polar vortex. During all the seasons, the enhancement of GWs around the geomagnetic poles as observed by GOCE (at ~250 km) is well captured by simulations. WAMT254 GWs in the IT region also show dependence on local time due to their interaction with migrating tides leading to diurnal and semidiurnal variations. The GWs are more likely to propagate up from the MLT region during westward/weakly-eastward phase of thermospheric tides, signifying the dominance of eastward GW momentum flux in the MLT. Additionally, as a novel finding, a wavenumber-4 signature in GW activity is predicted by WAMT254 between 6-12 LT in the tropics at ~250 km, which propagates eastward with local time. This behavior is likely associated with the modulation of GWs by wave-4 signal of non-migrating tides in the lower thermospheric zonal winds.

The mechanical complexity of flapping wings, their unsteady aerodynamic flow, and challenge of making measurements at the scale of a sub-gram flapping-wing flying insect robot (FIR) make its behavior hard to predict. Knowing the precise mapping from voltage input to torque output, however, can be used to improve their mechanical and flight controller design. To address this challenge, we created a sensitive force-torque sensor based on a flexured gimbal that only requires a standard motion capture system or accelerometer for readout. Our device precisely and accurately measures pitch and roll torques simultaneously, as well as thrust, on a tethered flapping-wing FIR in response to changing voltage input signals. With it, we were able to measure cross-axis coupling of both torque and thrust input commands on a 180 mg FIR, the UW Robofly. We validated these measurements using free-flight experiments. Our results showed that roll and pitch have maximum cross-axis coupling errors of 8.58% and 17.24%, respectively, relative to the range of torque that is possible. Similarly, varying the pitch and roll commands resulted in up to a 5.78% deviation from the commanded thrust, across the entire commanded torque range. Our system, the first to measure two torque axes simultaneously, shows that torque commands have a negligible cross-axis coupling on both torque and thrust., Comment: This work has been submitted to the IEEE for possible publication

Flying insects can perform rapid, sophisticated maneuvers like backflips, sharp banked turns, and in-flight collision recovery. To emulate these in aerial robots weighing less than a gram, known as flying insect robots (FIRs), a fast and responsive control system is essential. To date, these have largely been, at their core, elaborations of proportional-integral-derivative (PID)-type feedback control. Without exception, their gains have been painstakingly tuned by hand. Aggressive maneuvers have further required task-specific tuning. Optimal control has the potential to mitigate these issues, but has to date only been demonstrated using approxiate models and receding horizon controllers (RHC) that are too computationally demanding to be carried out onboard the robot. Here we used a more accurate stroke-averaged model of forces and torques to implement the first demonstration of optimal control on an FIR that is computationally efficient enough to be performed by a microprocessor carried onboard. We took force and torque measurements from a 150 mg FIR, the UW Robofly, using a custom-built sensitive force-torque sensor, and validated them using motion capture data in free flight. We demonstrated stable hovering (RMS error of about 4 cm) and trajectory tracking maneuvers at translational velocities up to 25 cm/s using an optimal linear quadratic regulator (LQR). These results were enabled by a more accurate model and lay the foundation for future work that uses our improved model and optimal controller in conjunction with recent advances in low-power receding horizon control to perform accurate aggressive maneuvers without iterative, task-specific tuning., Comment: The video of the results can be accessed using www.youtube.com/watch?v=0o7j1nS2KHA

The ionisation cones of NGC5728 have a deficit of molecular gas based on millimetre observations of CO(2-1) emission. Although photoionisation from the active nucleus may lead to suppression of this transition, warm molecular gas can still be present. We report the detection of eight mid-infrared rotational H$_2$ lines throughout the central kiloparsec, including the ionisation cones, using integral field spectroscopic observations with JWST/MIRI MRS. The H$_2$ line ratios, characteristic of a power-law temperature distribution, indicate that the gas is warmest where it enters the ionisation cone through disk rotation, suggestive of shock excitation. In the nucleus, where the data can be combined with an additional seven ro-vibrational H$_2$ transitions, we find that moderate velocity (30 km s$^{-1}$) shocks in dense ($10^5$ cm$^{-3}$) gas, irradiated by an external UV field ($G_0 = 10^3$), do provide a good match to the full set. The warm molecular gas in the ionisation cone that is traced by the H$_2$ rotational lines has been heated to temperatures $>200$ K. Outside of the ionisation cone the molecular gas kinematics are undisturbed. However, within the ionisation cone, the kinematics are substantially perturbed, indicative of a radial flow, but one that is quantitatively different from the ionised lines. We argue that this outflow is in the plane of the disk, implying a short 50 pc acceleration zone up to speeds of about 400 km s$^{-1}$ followed by an extended deceleration over $\sim$700 pc where it terminates. The deceleration is due to both the radially increasing galaxy mass, and mass-loading as ambient gas in the disk is swept up., Comment: A&A accepted; 16 pages

We calculate the extent to which collisionless dark matter impacts the stability of supermassive stars $(M\gtrsim10^4\,M_\odot)$. We find that, depending on the star's mass, a dark matter content in excess of ${\sim}1\%$ by mass throughout the entire star can raise the critical central density for the onset general relativistic instability, in some cases by orders of magnitude. We consider implications of this effect for the onset of nuclear burning and significant neutrino energy losses., Comment: 13 pages, 4 figures V3: added references, added section of new analysis with two new figures, overall results unaltered, added journal reference

Binary mass transfer can occur at high rates due to rapid expansion of the donor's envelope. In the case where mass transfer is unstable, the binary can rapidly shrink its orbit and lead to a merger. In this work we consider the appearance of the system preceding merger, specifically for the case of a low-mass ($\approx 2.5$-$3~M_\odot$) helium star with a neutron star (NS) companion. Modeling the mass transfer history as well as the wind launched by super-Eddington accretion onto the NS, we find that such systems can power slowly rising transients with timescales as long as years, and luminosities of $\sim 10^{40}$-$10^{41}$ erg s$^{-1}$ from optical to UV. The final explosion following the merger (or core-collapse of the helium star in some cases) leads to an interaction-powered transient with properties resembling Type Ibn supernovae (SNe), possibly with a bright early peak powered by shock cooling emission for merger-powered explosions. We apply our model to the Type Ibn SN 2023fyq, that displayed a long-term precursor activity from years before the terminal explosion., Comment: Accepted to the Open Journal of Astrophysics. A typographic error in equation 34 has been corrected, and correspondingly Figure 5 has been slightly updated (other results are unaffected). Light curve source code available here: https://github.com/DTsuna/merger_precursor.git

The mass loss mechanism of red supergiant stars is not well understood, even though it has crucial consequences for their stellar evolution and the appearance of supernovae that occur upon core-collapse. We argue that outgoing shock waves launched near the photosphere can support a dense chromosphere between the star's surface and the dust formation radius at several stellar radii. We derive analytic expressions for the time-averaged density profile of the chromosphere, and we use these to estimate mass loss rates due to winds launched by radiation pressure at the dust formation radius. These mass loss rates are similar to recent observations, possibly explaining the upward kink in mass loss rates of luminous red supergiants. Our models predict that low-mass red supergiants lose less mass than commonly assumed, while high-mass red supergiants lose more. The chromospheric mass of our models is $\sim$0.01 solar masses, most of which lies within a few stellar radii. This can help explain the early light curves and spectra of type-II P supernovae without requiring extreme pre-supernova mass loss. We discuss implications for stellar evolution, type II-P supernovae, SN 2023ixf, and Betelgeuse., Comment: Accepted to the Open Journal of Astrophysics

High-resolution images offer more information about scenes that can improve model accuracy. However, the dominant model architecture in computer vision, the vision transformer (ViT), cannot effectively leverage larger images without finetuning -- ViTs poorly extrapolate to more patches at test time, although transformers offer sequence length flexibility. We attribute this shortcoming to the current patch position encoding methods, which create a distribution shift when extrapolating. We propose a drop-in replacement for the position encoding of plain ViTs that restricts attention heads to fixed fields of view, pointed in different directions, using 2D attention masks. Our novel method, called LookHere, provides translation-equivariance, ensures attention head diversity, and limits the distribution shift that attention heads face when extrapolating. We demonstrate that LookHere improves performance on classification (avg. 1.6%), against adversarial attack (avg. 5.4%), and decreases calibration error (avg. 1.5%) -- on ImageNet without extrapolation. With extrapolation, LookHere outperforms the current SoTA position encoding method, 2D-RoPE, by 21.7% on ImageNet when trained at $224^2$ px and tested at $1024^2$ px. Additionally, we release a high-resolution test set to improve the evaluation of high-resolution image classifiers, called ImageNet-HR., Comment: NeurIPS 2024 Camera Ready

To study the impact of the initial effects of metallicity (i.e., the abundance of elements heavier than helium) on star formation and the formation of different molecular species, we searched for hot molecular cores in the sub-solar metallicity environment of the Large Magellanic Cloud (LMC). We conducted an ALMA Band 6 observations of 20 fields centered on young stellar objects (YSOs) distributed over the LMC in order to search for hot molecular cores in this galaxy. We detected a total of 65 compact 1.2 mm continuum cores in the 20 ALMA fields and analyzed their spectra with XCLASS software. The main temperature tracers are CH3OH and SO2, with more than two transitions detected in the observed frequency ranges. Other molecular lines with high detection rates in our sample are CS , SO, H13CO+, H13CN, HC15 N, and SiO. More complex molecules, such as HNCO, HDCO, HC3N, CH3CN, and NH2CHO, and multiple transitions of SO and SO2 isotopologues showed tentative or definite detection toward a small subset of the cores. According to the chemical richness of the cores and high temperatures from the XCLASS fitting, we report the detection of four hot cores and one hot core candidate. With one new hot core detection in this study, the number of detected hot cores in the LMC increases to seven. Six out of seven hot cores detected in the LMC to date are located in the stellar bar region of this galaxy. These six hot cores show emission from complex organic molecules (COMs), such as CH3OH, CH3CN, CH3OCHO, and CH3OCH3. The only known hot core in the LMC with no detection of COMs is located outside the bar region. The metallicity in the LMC presents a shallow gradient increasing from outer regions toward the bar. We suggest that the formation of hot molecular cores containing COMs ensues from the new generation of stars forming in the more metal-rich environment of the LMC bar.

A patient's digital twin is a computational model that describes the evolution of their health over time. Digital twins have the potential to revolutionize medicine by enabling individual-level computer simulations of human health, which can be used to conduct more efficient clinical trials or to recommend personalized treatment options. Due to the overwhelming complexity of human biology, machine learning approaches that leverage large datasets of historical patients' longitudinal health records to generate patients' digital twins are more tractable than potential mechanistic models. In this manuscript, we describe a neural network architecture that can learn conditional generative models of clinical trajectories, which we call Digital Twin Generators (DTGs), that can create digital twins of individual patients. We show that the same neural network architecture can be trained to generate accurate digital twins for patients across 13 different indications simply by changing the training set and tuning hyperparameters. By introducing a general purpose architecture, we aim to unlock the ability to scale machine learning approaches to larger datasets and across more indications so that a digital twin could be created for any patient in the world.