Papers for Friday, Mar 10 2023

We show that dynamical Chern-Simons (dCS) gravity imprints a parity-violating signal in primordial scalar perturbations. Specifically, we find that, after dCS amplifies one graviton helicity due to a tachyonic instability, the graviton-mediated correlation between two pairs of scalars develops a parity-odd component. This correlation, the primordial scalar trispectrum, is then transferred to the corresponding curvature correlator and thus is imprinted in both LSS and the CMB. We find that the parity-odd piece has roughly the same amplitude as its parity-even counterpart, scaled linearly by the degree of gravitational circular polarization $\Pi_{\rm circ} \sim \sqrt{\varepsilon}[H^2/(M_{\rm Pl} f)] \leq 1$, with $\varepsilon$ the slow-roll parameter, $H$ the inflationary Hubble scale, $f$ the dCS decay constant, and the upper bound saturated for purely circularly-polarized gravitons. We also find that, in the collapsed limit, the ratio of the two trispectra contains direct information about the graviton's spin. In models beyond standard inflationary dCS, e.g. those with multiple scalar fields or superluminal scalar sound speed, there can be a large enhancement factor $F \gtrsim 10^6$ to the trispectrum. We find that an LSS survey that contains $N_{\rm modes}$ linear modes would place an $n\sigma$ constraint on $\Pi_{\rm circ}r$ of $\sim 0.04\ (n/3)(10^6/F)(10^6/N_{\rm modes})^{1/2}$ from the parity-odd galaxy trispectrum, for tensor-to-scalar ratio $r$. We also forecast for several spectroscopic and 21-cm surveys. This constraint implies that, for high-scale single-field inflation parameters, LSS can probe very large dCS decay constants $f \lesssim 4\times 10^9\ {\rm GeV}(3/n)(F/10^6)\left(N_{\rm modes}/10^6\right)^{1/2}$. Our result is the first example of a massless particle yielding a parity-odd scalar trispectrum through spin-exchange.

We measure the extinction law in the 30 Dor star formation region in the Large Magellanic Cloud using Early Release Observations taken with Near-Infrared Camera (NIRCam) onboard the JWST, thereby extending previous studies with the Hubble Space Telescope to the infrared. We use red clump stars to derive the direction of the reddening vector in twelve bands and present the extinction law in this massive star forming region from $0.3$ to $4.7\,\mu$m. At wavelengths longer than 1 $\mu$m, we find a ratio of total and selective extinction twice as high as in the diffuse Milky Way interstellar medium and a change in the relative slope from the optical to the infrared domain. Additionally, we derive an infrared extinction map and find that extinction closely follows the highly embedded regions of 30 Dor.

The investigation of the feedback cycle in galaxy clusters has historically been performed for systems where feedback is ongoing ("mature-feedback" clusters), that is where the central radio galaxy has inflated radio lobes, pushing aside the intracluster medium (ICM). In this pilot study we present results from "pre-feedback" clusters, where the central newly active radio galaxies (age $<10^{3}$ yr) may not yet have had time to alter the thermodynamic state of the ICM. We analyze $Chandra$ and MUSE observations of two such systems, evaluating the hot gas entropy and cooling time profiles, and characterizing the morphology and kinematics of the warm gas. Based on our exploratory study of these two sources, we find that the hot gas meets the expectations for an as-yet unheated ICM. Specifically, the entropy and cooling time of pre-feedback clusters within 20 kpc from the center fall below those of mature-feedback clusters by a factor $\sim$2. We speculate that with an estimated mechanical power of $\sim10^{44} - 10^{45}$ erg s$^{-1}$, the two young radio galaxies may restore the entropy levels in a few tens of Myr, which are typical values of power outbursts and lifetimes for radio galaxies in clusters. Conversely, the properties of the gas at $\sim10^{4}$ K seem to remain invariant between the two feedback stages, possibly suggesting that the warm gas reservoir accumulates over long periods ($10^{7}$ - $10^{8}$ yr) during the growth of the radio galaxy. We conclude that the two cluster-central young radio galaxies we analyzed provide an interesting case to study the onset of AGN feedback, but the investigation of other similar sources is essential to confirm our results.

We use the ASTRID cosmological hydrodynamic simulation to investigate the properties and evolution of triple and quadruple Massive Black Hole (MBH) systems at $z = 2-3$. Only a handful of MBH tuple systems have been detected to date. In ASTRID, we find $4\%$ of the $M_{\rm BH}>10^7\,M_\odot$ are in tuples with $\Delta r_{\rm max} < 200\,{\rm kpc}$. The tuple systems span a range of separations with the majority of the observable AGN systems at $\Delta r \sim 50-100$ kpc. They include some of the most massive BHs (up to $10^{10} \,M_\odot$) but with at least one of the components of $M_{\rm BH} \sim 10^7 \,M_\odot$. Tuples' host galaxies are typically massive with $M_* \sim 10^{10-11} \,M_\odot$. We find that $>10\%$ massive halos with $M_{\rm halo} > 10^{13} M_\odot$ host MBH tuples. Following the subsequent interactions between MBHs in tuples, we found that in $\sim 5\%$ of the triplets all three MBHs merge within a Gyr, and $15\%$ go through one merger. As a by-product of the complex multi-galaxy interaction of these systems, we also find that up to $\sim 5\%$ of tuples lead to runaway MBHs. In ASTRID, virtually all of the ultramassive black holes ($>10^{10} \,M_\odot $) have undergone a triple quasar phase while for BHs with $M_{\rm BH} \sim 10^9 \,M_\odot$ this fraction drops to $50\%$.

JWST observations reveal a surprising excess of luminous galaxies at $z\sim 10$, consistent with efficient conversion of the accreted gas into stars, unlike the suppression of star formation by feedback at later times. We show that the high densities and low metallicities at this epoch guarantee a high star-formation efficiency in the most massive dark-matter haloes. Feedback-free starbursts (FFBs) occur when the free-fall time is shorter than $\sim 1$ Myr, below the time for low-metallicity massive stars to develop winds and supernovae. This corresponds to a characteristic density of $\sim 3\times 10^3$cm$^{-3}$. A comparable threshold density permits a starburst by allowing cooling to star-forming temperatures in a free-fall time. The galaxies within $\sim 10^{11} M_\odot$ haloes at $z \sim 10$ are expected to have FFB densities. The halo masses allow efficient gas supply by cold streams in a halo crossing time $\sim 80$ Myr. The FFBs gradually turn all the accreted gas into stars in clusters of $\sim 10^{4-7.5} M_\odot$ within galaxies that are rotating discs or shells. The starbursting clouds are shielded against feedback from earlier stars. We predict high star-formation efficiency above thresholds in redshift and halo mass, where the density is $10^{3-4}$cm$^{-3}$. The $z\sim 10$ haloes of $\sim 10^{10.8} M_\odot$ are predicted to host galaxies of $\sim 10^{10} M_\odot$ with SFR $\sim 65 M_\odot$ yr$^{-1}$ and sub-kpc sizes. The metallicity is $\leq 0.1 Z_\odot$ with little gas, dust, outflows and hot circumgalactic gas, allowing a top-heavy IMF but not requiring it. The post-FFB evolution of compact galaxies with thousands of young clusters may have implications on black-hole growth and globular clusters at later times.

Many barred galaxies exhibit upturns (shoulders) in their bar major-axis density profile. Simulation studies have suggested that shoulders are supported by looped x1 orbits, are present in growing bars, and can appear after bar buckling. Here we investigate the orbital support and evolution of shoulders via frequency analyses of orbits in simulations. We confirm that looped orbits are important to shoulders, and can remain so, to a lesser extent, after they have been vertically thickened. We show that looped orbits are found at the Inner Lindblad Resonance (ILR) with typical ratios of vertical to radial frequencies $1 \lesssim\Omega_z/\Omega_R \lesssim 3/2$ ($\textit{warm}$ ILR). $\textit{Cool}$ ILR orbits (those with $\Omega_z/\Omega_R > 3/2$), which are vertically thin and have no loops, contribute negligibly to shoulders. As bars slow and thicken, either secularly or by buckling, they populate warm ILR orbits. Further thickening carries these orbits towards the vertical ILR (vILR) at $\Omega_z/\Omega_R=1$, where they convert in-plane motion to vertical motion, become chaotic, kinematically hotter and less shoulder-supporting. Because of this heating by the vILR, persistent shoulders require bars to trap new orbits, consistent with the need for a growing bar. Since buckling speeds up the trapping of stars on warm ILR orbits, it can be followed by the formation of shoulders, as seen in simulations. This evolutionary sequence supports the recent observational finding that shoulder formation likely precedes the emergence of BP-bulges. The python module used for the frequency analysis, $\texttt{naif}$, is made publicly available.

We analyse the molecular and atomic emission lines of 71 bright Herschel-selected galaxies between redshifts 1.4 to 4.6 detected by the Atacama Large Millimetre/submillimetre Array. These lines include a total of 156 CO, [C I], and H2O emission lines. For 46 galaxies, we detect two transitions of CO lines, and for these galaxies we find gas properties similar to those of other dusty star-forming galaxy (DSFG) samples. A comparison to photo-dissociation models suggests that most of Herschel-selected galaxies have similar interstellar medium conditions as local infrared-luminous galaxies and high-redshift DSFGs, although with denser gas and more intense far-ultraviolet radiation fields than normal star-forming galaxies. The line luminosities agree with the luminosity scaling relations across five orders of magnitude, although the star-formation and gas surface density distributions (i.e., Schmidt-Kennicutt relation) suggest a different star-formation phase in our galaxies (and other DSFGs) compared to local and low-redshift gas-rich, normal star-forming systems. The gas-to-dust ratios of these galaxies are similar to Milky Way values, with no apparent redshift evolution. Four of 46 sources appear to have CO line ratios in excess of the expected maximum (thermalized) profile, suggesting a rare phase in the evolution of DSFGs. Finally, we create a deep stacked spectrum over a wide rest-frame frequency (220-890 GHz) that reveals faint transitions from HCN and CH, in line with previous stacking experiments.

Drawing connections between heliospheric spacecraft and solar wind sources is a vital step in understanding the evolution of the solar corona into the solar wind and contextualizing \textit{in situ} timeseries. Furthermore, making advanced predictions of this linkage for ongoing heliospheric missions, such as Parker Solar Probe (PSP), is necessary for achieving useful coordinated remote observations and maximizing scientific return. The general procedure for estimating such connectivity is straightforward (i.e. magnetic field line tracing in a coronal model) but validating the resulting estimates difficult due to the lack of an independent ground truth and limited model constraints. In its most recent orbits, PSP has reached perihelia of 13.3$R_\odot$ and moreover travels extremely fast prograde relative to the solar surface, covering over 120 degrees longitude in three days. Here we present footpoint predictions and subsequent validation efforts for PSP Encounter 10, the first of the 13.3$R_\odot$ orbits, which occurred in November 2021. We show that the longitudinal dependence of \textit{in situ} plasma data from these novel orbits provides a powerful method of footpoint validation. With reference to other encounters, we also illustrate that the conditions under which source mapping is most accurate for near-ecliptic spacecraft (such as PSP) occur when solar activity is low, but also requires that the heliospheric current sheet is strongly warped by mid-latitude or equatorial coronal holes. Lastly, we comment on the large-scale coronal structure implied by the Encounter 10 mapping, highlighting an empirical equatorial cut of the Alfv\`{e}n surface consisting of localized protrusions above unipolar magnetic separatrices.

Brown dwarf spectra offer vital testbeds for our understanding of the chemical and physical processes that sculpt substellar atmospheres. Recently, atmospheric retrieval approaches have been applied to a number of low-resolution (R~100) spectra of brown dwarfs, yielding constraints on the abundances of chemical species and temperature structures of these atmospheres. Medium-resolution (R~1e3) spectra of brown dwarfs offer significant additional insight, as molecular features are more easily disentangled from one another and the thermal structure of the upper atmosphere is more readily probed. We present results from a GPU-based retrieval analysis of a high signal-to-noise, medium-resolution (R~6000) FIRE spectrum from 0.85-2.5 microns of a T9 dwarf. At 60x higher spectral resolution than previous brown dwarf retrievals, a number of novel challenges arise. We examine the strong effect of different opacity sources on our retrieved constraints, in particular for CH4. Furthermore, we find that flaws in the data such as errors from order stitching can greatly bias our results. We compare these results to those obtained for a R~100 spectrum of the same object, revealing how constraints on atmospheric abundances and temperatures improve by an order of magnitude or more (depending on the species) with increased spectral resolution. In particular, we precisely constrain the abundance of H2S, which is undetectable at lower spectral resolution. While these medium-resolution retrievals offer the potential of precise, stellar-like constraints on atmospheric abundances (~0.02 dex), our retrieved radius is unphysically small (R~0.50 R$_{Jup}$), indicating lingering shortcomings with our modeling framework. This work is an initial investigation into brown dwarf retrievals at medium spectral resolution, offering guidance for future ground-based studies and JWST observations of substellar objects.

Two-body scatterings under the potential of a massive object are very common in astrophysics. If the massive body is far enough away that the two small bodies are in their own gravitational sphere of influence, the gravity of the massive body can be temporarily ignored. However, this requires the scattering process to be fast enough that the small objects do not spend too much time at distances near the surface of the sphere of influence. In this paper, we derive the validation criteria for effective two-body scattering and establish a simple analytical solution for this process, which we verify through numerical scattering experiments. We use this solution to study star-black hole scatterings in the disks of Active Galactic Nuclei and planet-planet scatterings in planetary systems, and calculate their one-dimensional cross-section analytically. Our solution will be valuable in reducing computational time when treating two-body scatterings under the potential of a much more massive third body, provided that the problem settings are in the valid parameter space region identified by our study.

T CrB is one of the most-famous and brightest novae known, and is a recurrent nova with prior eruptions in 1866 and 1946 that peak at $V$=2.0. I have constructed light curves spanning 1842--2022 with 213,730 magnitudes, where the $B$ and $V$ magnitudes are fully corrected to the Johnson system. These light curves first reveal a unique complex high-state (with 20$\times$ higher accretion rate than the normal low-state) stretching from -10 to +9 years after eruption, punctuated with a deep pre-eruption dip (apparently from dust formation in a slow mass ejection) and a unique enigmatic secondary eruption (with 10 per cent of the energy of the primary eruption), with the light curves identical for the 1866 and 1946 eruptions. Starting in 2015, T CrB entered the high-state, like in 1936, so a third eruption in upcoming years has been widely anticipated. With the pre-1946 light curve as a template, I predict a date of 2025.5$\pm$1.3 for the upcoming eruption, with the primary uncertainty arising from a possible lengthening of the pre-eruption high-state. I use the large-amplitude ellipsoidal modulation to track the orbital phase of the binary from 1867--2022. I measure that the orbital period increased abruptly by $+$0.185$\pm$0.056 days across the 1946 eruption, the 1947--2022 years had a steady period decrease of ($-$8.9$\pm$1.6)$\times$10$^{-6}$ days-per-day, and the 1867--1946 years had a steady period change consistent with zero, at ($+$1.75$\pm$4.5)$\times$10$^{-6}$ days-per-day. These large period changes cannot be explained by any published mechanism.

A comprehensive study of the nucleus and western hotspot of Pictor A is carried out using AstroSat observations, 13 years of Fermi, and archival Swift observations, along with other published data. We report the first detection of the western hotspot of Pictor A in the far-UV band using observations from AstroSat-UVIT. The broad-band SED of the western hotspot is explained by a multizone emission scenario where X-ray emission is caused by synchrotron emission process in the substructures embedded in the diffuse region, while the emission in radio to optical is caused by synchrotron emission process in the diffuse region. We do not notice any excess in the IR band and an additional zone (beyond 2-zone) is not required to account for the X-ray emission. Our broad-band spectro-temporal study and associated modelling of the core and hotspot of Pictor A suggests that (a) gamma-rays originate in the nuclear jet and not from the hotspot (b) X-ray emission from the core of Pictor A has nuclear jet-origin instead of previously reported disk-origin.

We introduce a practical methodology for investigating the star formation and chemical evolution history of a galaxy: age-divided mean stellar populations (ADPs) from full spectrum fitting. In this method, the mass-weighted mean stellar populations and mass fractions (f_mass) of young and old stellar components in a galaxy are separately estimated, which are divided with an age cut (selected to be 10^9.5 yr ~ 3.2 Gyr in this paper). To examine the statistical reliability of ADPs, we generate 10,000 artificial galaxy spectra, each of which consists of five random simple stellar population components. Using the Penalized PiXel-Fitting (pPXF) package, we conduct full spectrum fitting to the artificial spectra with noise as a function of wavelength, imitating the real noise of Sydney-Australian Astronomical Observatory Multi-object Integral field spectrograph (SAMI) galaxies. As a result, the \Delta (= output - input) of age and metallicity appears to significantly depend on not only signal-to-noise ratio (S/N), but also luminosity fractions (f_lum) of young and old components. At given S/N and f_lum, \Delta of young components tends to be larger than \Delta of old components; e.g., \sigma(\Delta [M/H]) ~ 0.40 versus 0.23 at S/N = 30 and f_lum = 50 per cent. The age-metallicity degeneracy appears to be insignificant, but \Delta log(age/yr) shows an obvious correlation with \Delta f_mass for young stellar components (R ~ 0.6). The impact of dust attenuation and emission lines appears to be mostly insignificant. We discuss how this methodology can be applied to spectroscopic studies of the formation histories of galaxies, with a few examples of SAMI galaxies.

Geometric optical distortion is a significant contributor to the astrometric error budget in large telescopes using adaptive optics. To increase astrometric precision, optical distortion calibration is necessary. Modern distortion calibration systems use back-illuminated pinhole masks to image a regular grid of point sources through a system's optics. The optical distortion of the instrument is then calibrated to an astrometric precision limited by the source placement error on the mask. Because of this, pinhole masks require an extreme level of precision, making their manufacture challenging and expensive. Additionally, they must be designed for a particular system, rendering them incompatible between instruments. We investigate using smartphone OLED screens as astrometric calibrators. Smartphones are low cost, have stable illumination, and can be quickly reconfigured to probe different spatial frequencies in the optical system. In this work, we characterize the astrometric accuracy of a Samsung S20 smartphone, with a view towards providing large format, flexible astrometric calibrators for the next generation of astronomical instruments. Investigating the placement error of 826 green OLED pixels, the non-linear deviations are measured to be 189 nm +/- 15 nm RMS. At this level of error, milli-arcsecond astrometric accuracy can be obtained on modern astronomical instruments.

We present Spitzer, WISE, and Herschel observations of the young supernova remnant (SNR) N132D in the LMC, including 3-40 microns Spitzer IRS mapping, 12 microns WISE and 70, 100, 160, 250, 350, and 500 microns Herschel images. The high-velocity lines of [Ne II] at 12.8 microns, [Ne III] at 15.5 microns, and [O IV] 26 microns reveal infrared ejecta concentrated in a central ring and coincide the optical and X-ray ejecta. Herschel images reveal far-IR emission coinciding with the central ejecta, which suggests that the IR emission is freshly formed, cold dust in the SN-ejecta. The infrared spectra are remarkably similar to those of another young SNR of 1E0102 with Ne and O lines. Shock modeling of the Ne ejecta emission suggests a gas temperature of 300 - 600 K and densities in the range 1000-20,000 cm^{-3} in the post-shock photoionized region. The IR continuum from the ejecta shows an 18 microns-peak dust feature. We performed spectral fitting to the IRS dust continuum and Herschel photometry. The dust mass associated with the central ejecta is 1.25+-0.65 Msun, while the 18 microns dust feature requires forsterite grains. The dust mass of the central ejecta region in N132D is higher than those of other young SNRs, which is likely associated with its higher progenitor mass. We discuss the dust productivity in the ejecta of N132D and infer its plausible implications for the dust in the early Universe.

Observationally, two main spectral states, i.e., the low/hard state and the high/soft state, are identified in black hole X-ray binaries (BH-XRBs). Meanwhile, the transitions between the two states are often observed. In this paper, we re-investigate the transition luminosities in the framework of the self-similar solution of the advection-dominated accretion flow (ADAF). Specifically, we search for the critical mass accretion rate $\dot m_{\rm crit}$ of ADAF for different radii $r$ respectively. It is found that $\dot m_{\rm crit}$ decreases with decreasing $r$. By testing the effects of BH mass $m$, the magnetic parameter $\beta$ and the viscosity parameter $\alpha$, it is found that only $\alpha$ has significant effects on $\dot m_{\rm crit}-r$ relation. We define the minimum $\dot m_{\rm crit}$ (roughly at the innermost stable circular orbit) as the hard-to-soft transition rate $\dot m_{\rm tr: H\rightarrow S}$, above which BH will gradually transit from the low/hard state to the high/soft state, and $\dot m_{\rm crit}$ at $30$ Schwarzschild radii as the soft-to-hard transition rate $\dot m_{\rm tr: S\rightarrow H}$, below which BH will gradually transit from the high/soft state to the low/hard state. We derive fitting formulae of $\dot m_{\rm tr: H\rightarrow S}$ and $\dot m_{\rm tr: S\rightarrow H}$ as functions of $\alpha$ respectively. By comparing with observations, it is found that the mean value of $\alpha$ are $\alpha \sim 0.85$ and $\alpha \sim 0.33$ for the hard-to-soft transition and the soft-to-hard transition respectively, which indicates that two classes of $\alpha$ are needed for explaining the hysteresis effect during the state transition. Finally, we argue that such a constrained $\alpha$ may provide valuable clues for further exploring the accretion physics in BH-XRBs.

The progenitor system of Type Ia supernovae (SNe Ia) is expected to be a close binary system of a carbon/oxygen white dwarf (WD) and a non-degenerate star or another WD. Here, we present results from a high-cadence monitoring observation of SN 2021hpr in a spiral galaxy, NGC 3147, and constraints on the progenitor system based on its early multi-color light curve data. First, we classify SN 2021hpr as a normal SN Ia from its long-term photometric and spectroscopic data. More interestingly, we found a significant "early excess" in the light curve over a simple power-law $\sim t^{2}$ evolution. The early light curve evolves from blue to red and blue during the first week. To explain this, we fitted the early part of $BVRI$-band light curves with a two-component model of the ejecta-companion interaction and a simple power-law model. The early excess and its color can be explained by shock cooling emission due to a companion star having a radius of $8.84\pm0.58$$R_{\odot}$. We also examined HST pre-explosion images with no detection of a progenitor candidate, consistent with the above result. However, we could not detect signs of a significant amount of the stripped mass from a non-degenerate companion star ($\lesssim0.003\,M_{\odot}$ for H$\alpha$ emission). The early excess light in the multi-band light curve supports a non-degenerate companion in the progenitor system of SN 2021hpr. At the same time, the non-detection of emission lines opens a door for other methods to explain this event.

Gradients in the mass-to-light ratio of distant galaxies impede our ability to characterize their size and compactness. The long-wavelength filters of $JWST$'s NIRCam offer a significant step forward. For galaxies at Cosmic Noon ($z\sim2$), this regime corresponds to the rest-frame near-infrared, which is less biased towards young stars and captures emission from the bulk of a galaxy's stellar population. We present an initial analysis of an extraordinary lensed dusty star-forming galaxy (DSFG) at $z=2.3$ behind the $El~Gordo$ cluster ($z=0.87$), named $El~Anzuelo$ ("The Fishhook") after its partial Einstein-ring morphology. The FUV-NIR SED suggests an intrinsic star formation rate of $81^{+7}_{-2}~M_\odot~{\rm yr}^{-1}$ and dust attenuation $A_V\approx 1.6$, in line with other DSFGs on the star-forming main sequence. We develop a parametric lens model to reconstruct the source-plane structure of dust imaged by the Atacama Large Millimeter/submillimeter Array, far-UV to optical light from $Hubble$, and near-IR imaging with 8 filters of $JWST$/NIRCam, as part of the Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS) program. The source-plane half-light radius is remarkably consistent from $\sim 1-4.5~\mu$m, despite a clear color gradient where the inferred galaxy center is redder than the outskirts. We interpret this to be the result of both a radially-decreasing gradient in attenuation and substantial spatial offsets between UV- and IR-emitting components. A spatial decomposition of the SED reveals modestly suppressed star formation in the inner kiloparsec, which suggests that we are witnessing the early stages of inside-out quenching.

Crescent-like asymmetric dust structures discovered in protoplanetary disks indicate dust aggregations. Thus, the research on them helps us understand the planet formation process. Here we analyze the ALMA data of the protoplanetary disk around the T-Tauri star SR 21, which has asymmetric structures detected in previous sub-millimeter observations. Imaged at ALMA Band 6 (1.3 mm) with a spatial resolution of about 0.$\arcsec$04, the disk is found to consist of two rings and three asymmetric structures, with two of the asymmetric structures being in the same ring. Compared to the Band 6 image, the Band 3 (2.7 mm) image also shows the three asymmetric structures but with some clumps. The elongated asymmetric structures in the outer ring could be due to the interactions of a growing planet. By fitting the Band 3 and Band 6 dust continuum data, two branches of solutions of maximum dust size in the disk are suggested: one is larger than 1 mm, and the other is smaller than 300 $\mu m$. High-resolution continuum observations at longer wavelengths as well as polarization observations can help break the degeneracy. We also suggest that the prominent spiral previously identified in VLT/SPHERE observations to the south of the star at 0.$\arcsec$25 may be the scattered light counterpart of the Inner Arc, and the structure is a dust-trapping vortex in nature. The discovered features in SR 21 make it a good target for studying the evolution of asymmetric structures and planet formation.

Although multidimensional simulations have investigated the processes of double WD mergers, post-merger evolution only focused on the carbon-oxygen (CO) WD or helium (He) WD merger remnants. In this work, we investigate for the first time the evolution of the remnants stemmed from the merger of oxygen-neon (ONe) WDs with CO WDs. Our simulation results indicate that the merger remnants can evolve to hydrogen- and helium-deficient giants with maximum radius of about 300Rsun. Our models show evidence that merger remnants more massive than 1.95Msun can ignite Ne before significant mass-loss ensues, and they thus would become electron-capture supernovae (ECSNe). However, remnants with initial masses less than 1.90Msun will experience further core contraction and longer evolutionary time before reaching at the conditions for Ne-burning. Therefore their fates are more dependent on mass-loss rates due to stellar winds, and thus more uncertain. Relatively high mass-loss rates would cause such remnants to end their lives as ONe WDs. Our evolutionary models can naturally explain the observational properties of the double WD merger remnant IRAS 00500+6713 (J005311). As previously suggested in the literature, we propose and justify that J005311 may be the remnant from the coalescence of an ONe WD and an CO WD. We deduce that the final outcome of J005311 would be a massive ONe WD rather than a supernova explosion. Our investigations may be able to provide possible constraints on the wind mass-loss properties of the giants which have CO-dominant envelopes.

In the framework of the Visible NEAs Observations Survey (ViNOS) that uses several telescopes at the Canary Islands observatories since 2018, we observed two super fast rotator NEAs, 2021 NY$_1$ and 2022 AB. We obtained photometry and spectrophotometry of both targets and visible spectroscopy of 2022 AB. Light curves of 2021 NY$_1$ obtained in 4 different nights between Sept. 30 and Oct. 16, 2021 return a rotation period $P=13.3449\pm0.0013$ minutes and a light curve amplitude $A = 1.00$ mag. We found that 2021 NY$_1$ is a very elongated super fast rotator with an axis ratio $a/b \ge 3.6$. We also report colours $(g-r) = 0.664 \pm 0.013$, $(r-i) = 0.186 \pm 0.013$, and $(i-z_s) = -0.117 \pm 0.012$ mag. These are compatible with an S-type asteroid. The light curves of 2022 AB obtained on Jan. 5 and Jan. 8, 2021 show a rotation period $P=3.0304\pm0.0008$ minutes, with amplitudes $A = 0.52$ and $A =0.54$ mag. 2022 AB is also an elongated object with axis ratio $a/b \ge 1.6$. The obtained colours are $(g-r) = 0.400 \pm 0.017$, $(r-i) = 0.133 \pm 0.017$, and $(i-z_s) = 0.093 \pm 0.016$. These colours are similar to those of the X-types, but with an unusually high $(g-r)$ value. Spectra obtained on Jan. 12 and Jan. 14, 2022, are consistent with the reported colours. The spectral upturn over the 0.4 - 0.6 $\mu m $ region of 2022 AB does not fit with any known asteroid taxonomical class or meteorite spectrum, confirming its unusual surface properties.

We derive the oblateness parameter q of the dark matter halo of a sample of gas rich, face-on disk galaxies. We have assumed that the halos are triaxial in shape but their axes in the disk plane (a and b) are equal, so that q=c/a measures the halo flattening. We have used the HI velocity dispersion, derived from the stacked HI emission lines and the disk surface density to determine the disk potential and the halo shape at the R25 and 1.5R25 radii. We have applied our model to 20 nearby galaxies, of which 6 are large disk galaxies with M(stellar)>10^10 solar mass, 8 have moderate stellar masses and 6 are low surface brightness dwarf galaxies. Our most important result is that gas rich galaxies that have M(gas)/M(baryons)>0.5 have oblate halos (q < 0.55), whereas stellar dominated galaxies have a range of q values, ranging from 0.21+-0.07 in NGC4190 to 1.27+-0.61 in NGC5194. Our results also suggest a positive correlation between the stellar mass and the halo oblateness q, which indicates that galaxies with massive stellar disks have a higher probability of having halos that are spherical or slightly prolate, whereas low mass galaxies have oblate halos (q < 0.55).

Neutron stars (NSs) can be used to constrain dark matter (DM) since a NS can transform into a black hole (BH) if it captures sufficient DM particles and exceeds the Chandrasekhar limit. We extend earlier work and for the first time take into account the Galactic motion of individual NSs, which changes the amount of the captured DM by as large as one to two orders of magnitude. We systematically apply the analysis to 414 NSs in the Milky Way, and constrain the DM particle mass and its interaction with nucleon simultaneously. We find that the most stringent bound is placed by a few NSs and the bound becomes stronger after considering the Galactic motion. The survival of observed NSs already excludes a cross section $\sigma_{nX}\gtrsim 10^{-45} \, {\rm cm}^2$ for DM particles with mass from $100\, {\rm MeV}$ to $10^3 \, {\rm GeV}$. Especially for a mass around $10 \, {\rm GeV}$, the constraint on the cross section is as stringent as $\sigma_{nX}\lesssim 10^{-49} \, {\rm cm}^2$.

{We simulate the evolution of relativistic electrons injected into the intracluster medium by five radio galaxies. We study the spatial transport and the emission properties of the injected radio plasma over a $\sim 5$ Gyr period, and the sequence of cooling and re-acceleration events experienced by electrons, using a Lagrangian approach joint with a numerical method to model the evolution of momentum spectra of relativistic electrons. When compared with electrons injected by shock waves, electrons injected by radio galaxies (here limited to a single injection event) in our tests are unable to fuel large, $\sim \rm ~Mpc$ sized radio relics with fossil electrons, as required by current theoretical models, while electrons previously seeded by other shocks can do this. On the other hand, the combination of seeding from radio galaxies, and of re-acceleration events from plasma perturbation, can produce detectable, small scale and filamentary emissions in the proximity ($\leq 100-200$ kpc) of radio galaxies.

In recent years, large radio surveys of Active Galactic Nuclei (AGNs), comprising millions of sources, have become available where one could investigate dipole asymmetries, assumedly arising due to a peculiar motion of the Solar system. Investigations of such dipoles have yielded in past much larger amplitudes (a factor of 2 to 20) than that of the CMB dipole, though their inferred directions in sky seem to lie, within statistical uncertainties, close to the CMB dipole. Here we investigate dipole asymmetries in two recent large radio surveys, Very Large Array Sky Survey (VLASS) containing 1.9 million sources, carried out at 3 GHz with the Very Large Array, covering the sky north of -40 deg dec, and the Rapid ASKAP Continuum Survey (RACS) containing 2.1 million sources, covering the sky south of 30 deg dec, carried out at 887.5 MHz. We find dipoles determined from the VLASS and RACS surveys to be significantly larger, about 4 and 7 times respectively, than the CMB dipole. Though directions of the VLASS and RACS dipoles differ significantly from each other, however, along with a number of other previously determined dipoles, they all appear to be pointing in a rather narrow sky region about the CMB dipole, with a joint probability of occurrence by a random chance to be less than $10^{-7}$, which argues for the various dipoles to be related somehow. Nonetheless, significant differences in their derived peculiar velocities, including that of the CMB, cannot be explained by a peculiar motion of the Solar system, which would necessarily be a single value. Instead, their discordant peculiar velocities may be indicating that different cosmic reference frames are moving relative to each other or that the matter distribution on cosmic scales is not homogeneous and isotropic, either scenario being in contravention of what expected from the Cosmological Principle (CP).

The formation process of fossil groups (FGs) is still under debate, and large samples of such objects are still missing. The aim of this paper is to increase the sample of known FGs, and to analyse the properties of their brightest group galaxies (BGG) and compare them with a control sample of non-FG BGGs. Based on the Tinker spectroscopic catalogue of haloes and galaxies, we extract 87 FG and 100 non-FG candidates. For all the objects with data available in UNIONS in the u and r bands, and/or in an extra r-band processed to preserve all low surface brightness features (rLSB), we made a 2D photometric fit of the BGG with GALFIT with one or two Sersic components and analysed how the subtraction of intracluster light contribution modifies the BGG properties. From the SDSS spectra available for the BGGs of 65 FGs and 82 non-FGs, we extracted the properties of their stellar populations with Firefly. We also investigated the origin of the emission lines in a nearby FG, NGC 4104, that has an AGN. A single Sersic profile can fit most objects in the u band, while two Sersics are needed in the r and rLSB bands, both for FGs and non-FGs. Non-FG BGGs cover a larger range of Sersic index. FG BGGs follow the Kormendy relation derived for almost one thousand brightest cluster galaxies (BCGs) by Chu et al. (2022) while non-FGs BGGs are mostly located below this relation, suggesting that FG BGGs have evolved similarly to BCGs, while non-FG BGGs have evolved differently. The above properties can be strongly modified by the subtraction of intracluster light contribution. The stellar populations of FG and non-FG BGGs do not differ significantly. Our results suggest FG and non-FG BGGs have had different formation histories, but it is not possible to trace differences in their stellar populations or large scale distributions.

Henri Deslandres initiated imaging spectroscopy of the solar atmosphere in 1892 at Paris observatory. He invented, concurrently with George Hale in Kenwood (USA) but quite independently, the spectroheliograph designed for monochromatic imagery of the Sun. Deslandres developed two kinds of spectrographs: the ''spectroh{\'e}liographe des formes'', i.e. the narrow bandpass instrument to reveal chromospheric structures (such as filaments, prominences, plages and active regions); and the ''spectroh{\'e}liographe des vitesses'', i.e. the ''section'' spectroheliograph to record line profiles of cross sections of the Sun, in order to measure the Dopplershifts of dynamic features. Deslandres moved to Meudon in 1898 with his instruments and improved the spectral and spatial resolutions, leading to the large quadruple spectroheliograph in 1908 developed with Lucien d'Azambuja. CaII K systematic observations started at this date and were followed in 1909 by H$\alpha$ with two dedicated 3-metres spectroheliographs. The observing service was organized by d'Azambuja who also intensively used the large 7-metres spectroheliograph for his research and thesis (1930). This paper summarizes fifty years of research by Mr and Mrs d'Azambuja, who explored various photospheric and chromospheric lines, performing special spectroheliograms with the high dispersion 7-metres instrument. They also observed intensively filaments and prominences (memoir published in 1948) and recorded rare solar activity events with the two 3-metres spectroheliographs, during the first half of the twentieth century.

Inhomogeneous reionization enhances the 1D Lyman-$\alpha$ forest power spectrum on large scales at redshifts $z\geq4$. This is due to coherent fluctuations in the ionized hydrogen fraction that arise from large-scale variations in the post-reionization gas temperature, which fade as the gas cools. It is therefore possible to use these relic fluctuations to constrain inhomogeneous reionization with the power spectrum at wavenumbers $\log_{10}(k/{\rm km^{-1}\,s})\lesssim -1.5$. We use the Sherwood-Relics suite of hybrid radiation hydrodynamical simulations to perform a first analysis of new Lyman-$\alpha$ forest power spectrum measurements at $4.0\leq z \leq 4.6$. These data extend to wavenumbers $\log_{10}(k/{\rm km^{-1}\,s})\simeq -3$, with a relative uncertainty of $10$--$20$ per cent in each wavenumber bin. Our analysis returns a $2.7\sigma$ preference for an enhancement in the Lyman-$\alpha$ forest power spectrum at large scales, in excess of that expected for a spatially uniform ultraviolet background. This large-scale enhancement could be a signature of inhomogeneous reionization, although the statistical precision of these data is not yet sufficient for obtaining a robust detection of the relic post-reionization fluctuations. We show that future power spectrum measurements with relative uncertainties of $\lesssim 2.5$ per cent should provide unambiguous evidence for an enhancement in the power spectrum on large scales.

We study the prospects of Machine Learning algorithms like Gaussian processes (GP) as a tool to reconstruct the Hubble parameter $H(z)$ with two upcoming gravitational wave missions, namely the evolved Laser Interferometer Space Antenna (eLISA) and the Einstein Telescope (ET). We perform non-parametric reconstructions of $H(z)$ with GP using realistically generated catalogues, assuming various background cosmological models, for each mission. We also take into account the effect of early-time and late-time priors separately on the reconstruction, and hence on the Hubble constant ($H_0$). Our analysis reveals that GPs are quite robust in reconstructing the expansion history of the Universe within the observational window of the specific mission under study. We further confirm that both eLISA and ET would be able to constrain $H(z)$ and $H_0$ to a much higher precision than possible today, and also find out their possible role in addressing the Hubble tension for each model, on a case-by-case basis.

Reactions of SiO molecules have been postulated to initiate efficient formation of silicate dust particles in outflows around dying (AGB) stars. Both OH radicals and H$_2$O molecules can be present in these environments and their reactions with SiO and the smallest SiO cluster, Si$_2$O$_2$, affect the efficiency of eventual dust formation. Rate coefficients of gas-phase oxidation and clustering reactions of SiO, Si$_2$O$_2$ and Si$_2$O$_3$ have been calculated using master equation calculations based on density functional theory calculations. The calculations show that the reactions involving OH are fast. Reactions involving H$_2$O are not efficient routes to oxidation but may under the right conditions lead to hydroxylated species. The reaction of Si$_2$O$_2$ with H$_2$O, which has been suggested as efficient producing Si$_2$O$_3$, is therefore not as efficient as previously thought. If H$_2$O molecules dissociate to form OH radicals, oxidation of SiO and dust formation could be accelerated. Kinetics simulations of oxygen-rich circumstellar environments using our proposed reaction scheme suggest that under typical conditions only small amounts of SiO$_2$ and Si$_2$O$_2$ are formed and that most of the silicon remains as molecular SiO.

The capture of the free-floating planets and primordial black holes into a collapsing protostellar cloud is considered. Although the last stage of rapid contraction leading to the star formation lasts for a relatively short time $\sim 10^5$ years, during this time there is a strong change in the gravitational potential created by the movement of the entire cloud's mass ($\sim M_\odot$). As a result, the probability of capturing an object into a contracting cloud is comparable to the probability of capturing into an already formed planetary system. Taking into account the collapse of the cloud increases by 70% the full probability of the planets capture at the orbits with large semi-axis $a<10^3$ au. Capture in the cloud can explain the wide inclined orbit of the supposed 9th planet in the solar system. At the same time, the probability of primordial black holes capturing from the galactic halo into a contracting cloud is extremely small.

Primordial black holes (PBHs) are black holes that might have formed in high density regions in the early universe. The presence of local-type non-Gaussianity can lead to large-scale fluctuations in the PBH formation rate. If PBHs make up a non-negligible fraction of dark matter, these fluctuations can appear as isocurvature modes, and be used to constrain the amplitude of non-Gaussianity. We build upon the results of previous work by extending the calculation to include peaks theory and making use of the compaction $C$ for the formation criteria, accounting for non-linearities between $C$ and the curvature perturbation $\zeta$. For quadratic models of non-Gaussianity, our updated calculation gives constraints that are largely unaltered compared to those previously found, while for cubic models the constraints worsen significantly. In case all of the DM is made up of PBHs, the parameters of non-Gaussianity are $-2.9\cdot10^{-4}<f<3.8\cdot10^{-4}$ and $-1.5\cdot10^{-3}<g<1.9\cdot10^{-3}$ for quadratic and cubic models respectively.

The main propose of this work is to investigate under which circumstances the tidal response of a stratified body can be approximated by that of a homogeneous body. We show that any multilayered body can be approximated by a homogeneous body, with the same dissipation of tidal energy as a function of the excitation frequency, as long as the rheology of the homogeneous model is sufficiently complex. Moreover we provide a straightforward recipe to associate the aforementioned homogeneous rheology to a given stratified body. These results highlight the fact that the two models cannot be distinguished from each other only by the measurement of the second degree tidal Love number and quality factor, and that we do not need the complexity of the multilayer planet model in order to estimate its tidal dissipation.

We suggest a common envelope jets supernova (CEJSN) origin to the supernova remnant (SNR) W49B where jets launched by a neutron star (NS) that collapsed to a black hole (BH) together with a thermonuclear outburst of the disrupted red super giant's (RGS's) core powered and shaped the ejecta. The jets account for the highly non-spherical morphology of W49B and the thermonuclear outburst to its high iron abundance. CEJSNe are violent events powered by jets that a NS or a BH launch as they orbit inside a red supergiant star and accrete mass from its envelope and then from its core. We classify the CEJSN process to either a case where the NS/BH enters the core to form a common envelope evolution (CEE) inside the core or to a case where the NS/BH tidally disrupts the core. In the later case the core material forms an accretion disk around the NS that might experience a thermonuclear outburst, leading to an energetic event powered by both jets and thermonuclear burning. We term this scenario thermonuclear CEJSN. We find that the maximum core mass that leads to this scenario with a NS is $2 M_{\rm \odot} \lesssim M_{\rm core} \lesssim 3.5M_{\rm \odot}$. We estimate the event rates of CEJSN that go through tidal disruption of the core by a NS to be 5 per 1000 core collapse supernovae.

While an astronomer's job is typically to look out from Earth, the seriousness of the climate crisis has meant a shift in many astronomers' focus. Astronomers are starting to consider how our resource requirements may contribute to this crisis and how we may better conduct our research in a more environmentally sustainable fashion. Astronomers for Planet Earth is an international organisation (more than 1,700 members from over 70 countries as of November 2022) that seeks to answer the call for sustainability to be at the heart of astronomers' practices. In this article, we review the organisation's history, summarising the proactive, collaborative efforts and research into astronomy sustainability conducted by its members. We update the state of affairs with respect to the carbon footprint of astronomy research, noting an improvement in renewable energy powering supercomputing facilities in Australia, reducing that component of our footprint by a factor of 2--3. We discuss how, despite accelerated changes made throughout the pandemic, we still must address the format of our meetings. Using recent annual meetings of the Australian and European astronomical societies as examples, we demonstrate that the more online-focussed a meeting is, the greater its attendance and the lower its emissions.

MHD models and the observation of accretion streamers confirmed that protostars can undergo late accretion events after the initial collapse phase. To provide better constraints, we study the evolution of stellar masses in MHD simulations of a (4 pc)^3 molecular cloud. Tracer particles allow us to accurately follow the trajectory of accreting material for all protostars and thereby constrain the accretion reservoir of the stars. The diversity of the accretion process implies that stars in the solar mass regime can have vastly different accretion histories. Some stars accrete most of their mass during the initial collapse phase, while others gain >50 % of their final mass from late infall. The angular momentum budget of stars that experience substantial infall, so-called late accretors, is significantly higher than for stars without or with only little late accretion. As the probability of late infall increases with increasing final stellar mass, the specific angular momentum budget of higher mass stars is on average higher. The hypothetical centrifugal radius computed from the accreting particles at the time of formation is orders of magnitude higher than observed disk sizes, which emphasizes the importance of angular momentum transport during disk formation. Nevertheless, we find a correlation that the centrifugal radius is highest for stars with substantial infall, which suggests that very large disks are the result of recent infall events. There are also indications for a subtle trend of increasing centrifugal radius with increasing final stellar mass, which is in agreement with an observed marginal correlation of disk size and stellar mass. Finally, we show that late accretors become embedded again during infall. As a consequence, late accretors are (apparently) rejuvenated and would be classified as Class 0 objects according to their bolometric temperature despite being 1 Myr old.

Compact Hierarchical Triples (CHT) are systems with the tertiary star orbiting the inner binary in an orbit shorter than 1000 days. CHT with an eclipsing binary as its inner binary can help us extract a multitude of information about all three stars in the system. In this study, we use independent observational techniques to estimate the orbital, stellar, and atmospheric parameters of two triple-lined CHT: BD+442258 and KIC06525196. We find that the masses of stars in BD+442258 are $1.011\pm0.029 M_{\odot}$, $0.941\pm0.033 M_{\odot}$, and $0.907\pm0.065 M_{\odot}$ while in KIC06525196 the estimated masses are $1.0351\pm0.0055 M_{\odot}$, $0.9712\pm0.0039 M_{\odot}$, and $0.777\pm0.012 M_{\odot}$. Using spectral disentangling, we obtained individual spectra of all the stars and combined them with light curve modeling to obtain radii, metallicities, and temperatures. Using stellar evolution models from MESA, we constrain the log(age) of BD+442258 to be 9.89 and 9.49 for KIC06525196. Two stars in BD+442258 are found to be sub-giants while all three stars in KIC06525196 are main-sequence stars. We constrain the mutual inclinations to certain angles for BD+442258 and KIC06525196 using numerical integration. Integrating with tidal interaction schemes and stellar evolution models, we find that KIC06525196 is a stable system. But the inner binary of BD+442258 merges within 550 Myrs. The time of this merger is affected by the orientation of the tertiary, even rushing the collapse by 100 Myrs when the mutual inclination is close to 90 degrees.

In this work, using \textit{NICER} and \textit{Insight}-HXMT observations, we present a study of the broadband spectral and timing evolution of the source throughout the first reflare, which occurred about 4 months after the major outburst. Our findings suggest that during the reflare, below a critical luminosity $L_{\rm crit}\sim2.5\times10^{36}$ (D/2.2 kpc)$^{2}$ erg s$^{-1}$, the scale of the corona shrinks in the radial direction, whereas the inner radius of the disk does not change considerably; however, the inner radius of the disk starts to move inward when the source exceeds the critical luminosity. We conclude that at low luminosity the increase in accretion rate only heats up the inner zone of the accretion disc without the transfer of angular momentum which occurs above a certain luminosity.

We present a new simulation setup using the MURaM radiative MHD code that allows to study the formation of collisional polarity inversion lines (cPILs) in the photosphere and the coronal response including flares. In the setup we start with a bipolar sunspot configuration and set the spots on collision course by imposing the appropriate velocity field at the footpoints in the subphotospheric boundary. We vary parameters such as the initial spot separation, collision speed and collision distance. While all setups lead to the formation of a sigmoid structure, only the cases with a close passing of the spots cause flares and mass eruptions. The energy release is in the $1-2 \times 10^{31}$ ergs range, putting the simulated flares into the upper C to lower M-class range. While the case with the more distant passing of the spots does not lead to a flare, the corona is nonetheless substantially heated, suggesting non-eruptive energy release mechanisms. We focus our discussion on two setups that differ in spot coherence and resulting cPIL length. We find different timings in the transition from a sheared magnetic arcade (SMA) to magnetic flux rope (MFR); the setup with a short cPIL produces a MFR during the eruption, while the MFR is pre-existing in the setup with a longer cPIL. While both result in flares of comparable strength, only the setup with pre-existing MFR produces a CME.

We conduct a binary population synthesis study to investigate the formation of wind-fed high-mass X-ray binaries containing black holes (BH-HMXBs). We evolve multiple populations of high-mass binary stars and consider BH-HMXB formation rates, masses, spins and separations. We find that systems similar to Cygnus X-1 likely form after stable Case A mass transfer (MT) from the main sequence progenitors of black holes, provided such MT is characterised by low accretion efficiency, $\beta \lesssim 0.1$, with modest orbital angular momentum losses from the non-accreted material. Additionally, BH-HMXB formation relies on a new simple treatment for Case A MT that allows donors to retain larger core masses compared to traditional rapid population-synthesis procedures. At solar metallicity, our Preferred model yields $\mathcal{O}(1)$ observable BH-HMXBs in the Galaxy today, consistent with observations. In this simulation, $8\%$ of BH-HMXBs go on to merge as binary black holes or neutron star-black hole binaries within a Hubble time. Compact-object mergers that evolve via a BH-HMXB phase may contribute $\gtrsim20$~Gpc$^{-3}$~yr$^{-1}$ to the merger rate inferred from gravitational-wave observations. We also suggest that MT efficiency is higher during stable Case B MT than during Case A MT.

With $n$-body simulations we investigate the stability of tilted circumbinary planetary systems consisting of two nonzero mass planets. The planets are initially in circular orbits that are coplanar to each other, as would be expected if they form in a flat but tilted circumbinary gas disc and decouple from the disc within a time difference that is much less than the disc nodal precession period. We constrain the parameters of stable multiple planet circumbinary systems. Both planet-planet and planet-binary interactions can cause complex planet tilt oscillations which can destabilise the orbits of one or both planets. The system is considerably more unstable than the effects of these individual interactions would suggest, due to the interplay between these two interactions. The stability of the system is sensitive to the binary eccentricity, the orbital tilt and the semi-major axes of the two circumbinary planets. With an inner planet semi-major axis of $5\,a_{\rm b}$, where $a_{\rm b}$ is semi-major axis of the binary, the system is generally stable if the outer planet is located at $\gtrsim 8\,a_{\rm b}$, beyond the 2:1 mean motion resonance with the inner planet. For larger inner planet semi-major axis the system is less stable because the von-Zeipel--Kozai--Lidov mechanism plays a significant role, particularly for low binary-eccentricity cases. For the unstable cases, the most likely outcome is that one planet is ejected and the other remains bound on a highly eccentric orbit. Therefore we suggest that this instability is an efficient mechanism for producing free-floating planets.

Radiative magnetohydrodynamic simulation includes sufficiently realistic physics to allow for the synthesis of remote sensing observables that can be quantitatively compared with observations. We analyze the largest flare in a simulation of the emergence of large flare-productive active regions described by Chen et al. The flare is accompanied by a spectacular coronal mass ejection and reaches M2 class, as measured from synthetic soft X-ray flux. The eruption reproduces many key features of observed solar eruptions. A pre-existing magnetic flux rope is formed along the highly sheared polarity inversion line between a sunspot pair and is covered by an overlying multi-pole magnetic field. During the eruption, the progenitor flux rope actively reconnects with the canopy field and evolves to the large-scale multi-thermal flux rope that is observed in the corona. Meanwhile, the magnetic energy released via reconnection is channeled down to the lower atmosphere and gives rise to bright soft X-ray post-flare loops and flare ribbons that reproduce the morphology and dynamic evolution of observed flares. The model helps to shed light on questions of where and when the a flux rope may form and how the magnetic structures in an eruption are related to observable emission properties.

We estimate the carbon footprint of astronomical research infrastructures, including space telescopes and probes and ground-based observatories. Our analysis suggests annual greenhouse gas emissions of $1.2\pm0.2$ MtCO$_2$e yr$^{-1}$ due to construction and operation of the world-fleet of astronomical observatories, corresponding to a carbon footprint of 36.6$\pm$14.0 tCO$_2$e per year and average astronomer. We show that decarbonising astronomical facilities is compromised by the continuous deployment of new facilities, suggesting that a significant reduction in the deployment pace of new facilities is needed to reduce the carbon footprint of astronomy. We propose measures that would bring astronomical activities more in line with the imperative to reduce the carbon footprint of all human activities.

A field with particularly exciting results over the past few years is the study of the interaction of cosmic rays with interstellar matter. For star formation to take place, gas and dust need to be sufficiently cold for gravity to overcome thermal pressure, and the ionisation fraction must be low enough to enable substantial decoupling between the gas and the Galactic magnetic field. As soon as the visual extinction is of the order of 3-4 magnitudes, the ultraviolet photon flux from the interstellar radiation field is fully quenched, thus the only source of ionisation and heating is provided by low-energy cosmic rays. We will briefly focus on the Galactic and local origin of cosmic rays and on their effects on medium ionisation.

Solar Energetic Particle (SEP) events and their major subclass, Solar Proton Events (SPEs), can result in unfavorable consequences to numerous aspects of life and technology, making them one of the most prevalent and harmful effects of solar activity. Garnering knowledge leading up to such events by studying proton and soft X-ray (SXR) flux data to alleviate the burdens they cause is therefore critical for their forecasting. Our previous SEP prediction study, Sadykov et al. 2021 indicated that it may be sufficient to utilize only proton and SXR parameters for SPE forecasts considering a limited data set from Solar Cycle (SC) 24. In this work we report the completion of a catalog of $\geq$ 10 MeV $\geq$10 particle flux unit (pfu) SPEs observed by Geostationary Operational Environmental Satellite (GOES) detectors operated by the National Oceanic and Atmospheric Administration (NOAA), with records of their properties spanning through SCs 22-24. We report an additional catalog of daily proton and SXR flux statistics. We use these catalogs to test the application of machine learning (ML) for the prediction of SPEs using a Support Vector Machine (SVM) algorithm. We explore how previous SCs can train and test on each other using both earlier and longer data sets during the training phase, evaluating how transferable an algorithm is across different time periods. Validation against the effects of cross-cycle transferability is an understudied area in SEP research, but should be considered for verifying the cross-cycle robustness of an ML-driven forecast.

Clouds more massive than about $10^5$ M$_\odot$ are potential sites of massive cluster formation. Studying the properties of such clouds in the early stages of their evolution offers an opportunity to test various cluster formation processes. We make use of CO, Herschel, and UKIDSS observations to study one such cloud, G148.24+00.41. Our results show the cloud to be of high mass ($\sim$ $1.1\times10^5$ M$_\odot$), low dust temperature ($\sim$ 14.5 K), nearly circular (projected radius $\sim$ 26 pc), and gravitationally bound with a dense gas fraction of $\sim 18$% and a density profile with a power-law index of $\sim -1.5$. Comparing its properties with those of nearby molecular clouds, we find that G148.24+00.41 is comparable to the Orion-A molecular cloud in terms of mass, size, and dense gas fraction. From our analyses, we find that the central area of the cloud is actively forming protostars and is moderately fractal with a Q-value of $\sim$ 0.66. We also find evidence of global mass-segregation in the cloud, with a degree of mass-segregation ($\Lambda_{MSR}) \approx3.2$. We discuss these results along with the structure and compactness of the cloud, the spatial and temporal distribution of embedded stellar population, and their correlation with the cold dust distribution, in the context of high-mass cluster formation. Comparing our results with models of star cluster formation, we conclude that the cloud has the potential to form a cluster in the mass range $\sim$ 2000--3000 M$_\odot$ through dynamical hierarchical collapse and assembly of both gas and stars.

We present a detailed spectro-polarimetric study of Black hole X-ray binary 4U $1630-47$ during its $2022$ outburst with IXPE and NICER observations. The source is observed in disk dominated thermal state (kT$_{in}\approx1.4$ keV) with clear detection of absorption features at $6.69\pm0.01$ keV and $6.97\pm0.01$ keV from both NICER as well as IXPE spectra, likely indicating a coupling of disk-wind. A significant degree of polarization(PD) $= 8.33\pm0.17\%$ and polarization angle(PA) $=17.78^{\circ}\pm0.60^{\circ}$ in the energy range of $2-8$ keV are measured with IXPE. PD is found to be an increasing function of energy whereas PA remains roughly same within the energy range. Simultaneous energy spectra from NICER in the range of $0.5-12$ keV are modelled to study the spectral properties. Furthermore, the spin parameter of the black hole is estimated with spectro-polarimetric data as a$_{\ast}=0.927\pm0.001\,(1\sigma)$ which is corroborated by NICER observations. Finally, we discuss the implications of our findings.

We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications. As in previous works, the turbulent mass entrainment rate extrapolated to nominal heating is unrealistically high ($1.58\times 10^{-4}\ \mathrm{M_\odot/yr}$), which is discussed in terms of the non-equilibrium response of the simulations to the initial stratification. We measure quantitatively the effect of mixing due to internal gravity waves excited by core convection interacting with the boundary in our simulations. The wave power spectral density as a function of frequency and wavelength agrees well with the GYRE eigenmode predictions based on the 1D spherically averaged radial profile. A diffusion coefficient profile that reproduces the spherically averaged abundance distribution evolution is determined for each simulation. Through a combination of eigenmode analysis and scaling relations it is shown that in the $N^2$-peak region, mixing is due to internal gravity waves and follows the scaling relation $D_\mathrm{IGW-hydro} \propto L^{4/3}$ over a $\gtrapprox 2\ \mathrm{dex}$ range of heating factors. Different extrapolations of the mixing efficiency down to nominal heating are discussed. If internal gravity wave mixing is due to thermally-enhanced shear mixing, an upper limit is $D_\mathrm{IGW} \lessapprox 2$ to $3\times 10^4\ \mathrm{cm^2/s}$ at nominal heating in the $N^2$-peak region above the convective core.

We propose structures of size between $\sim 1$ meter to 100 meters that drastically alter the local distribution of the Cosmic Neutrino Background ($C\nu B$). These structures have a shape reminiscent of a sea urchin: They consist of rods of width $w$ and length $L \gg w$ periodically arranged on the surface of sphere of radius $R\sim L$. Such a structure functions as a diffraction phase grating and produces a region around its center where the fractional neutrino-antineutrino asymmetry is $\sim k\delta_\nu L$, where $k$ is the neutrino momentum, and $\delta_\nu$ the deviation of the neutrino index of refraction from unity. The asymmetry has a gradient set by the rod width. We find that the local neutrino asymmetry can be enhanced by $\mathcal{O}(\text{few}\times 10^6)$ relative to the naive Standard Model expectation, for reasonably sized structures. This results in a force $\mathcal{O}(10^3)$ times bigger than the one we recently pointed out due to the neutrinos' reflection on the surface of the Earth. While in this paper we do not propose a concrete detection setup, we estimate that the $\mathcal{O}(G_F)$ force on a test mass can be close to the Standard Quantum Limit of a torsion balance or a low frequency harmonic oscillator. Finally, we show that this $C \nu B$ diffractor can be used as a Dark Matter diffractor. For example, the QCD axion Dark Matter with decay constant $f_a$ around $10^9$ GeV can be sufficiently diffracted to produce a gradient force that is up to $\mathcal{O}(10^2)$ times larger than the one from the $C \nu B$. This is the first setup of this kind and the simplicity of our design suggests that there could be significant improvements that escape us.

In this paper we use the effective Schr\"{o}dinger-Poisson and square-root Klein-Gordon-Poisson models to study the quantum and relativistic quantum energy band structure of finite temperature electron gas in a neutralizing charge background. Based the plasmon band gap appearing above the Fermi level, new definitions on plasmonic excitations and plasma parameters in a wide electron temperature-density regime is suggested. The new equation of state (EoS) for excited electrons to the plasmon band leads to novel aspects of relativistic collective quantum excitations such as the plasmon black-out and quantum pressure collapse which are studied using both non-relativistic and relativistic quantum models. The plasmon black-out effect may be used to explain why metallic elements do not show collective behavior at low temperatures. The model can be used to predict phases of matter in which the plasmonic activities is shut down, hence, it may behave like a mysterious dark matter. On the other hand, the energy band structure model predicts the plasmon pressure collapse in temperature-density coordinates matching that of a white dwarf star. The prediction of energy band structure of collective quantum excitations may have direct implications for the inertial confinement fusion (ICF), the EoS of warm dense matter (WDM) and evolution of stellar and other unknown cosmological structures. It is found that predictions of non-relativistic and relativistic quantum excitation models closely match up to temperature-density of degenerate stars which confirms the relevance of non-relativistic plasmon models used in the warm and dense matter regime. The effect of positron on band structure of collective quantum excitations is also studied.

In this work, I consider an inflation model with a quadratic potential and a negative cosmological constant. An analytical solution of the equation of motion for the inflaton field is found without slow-roll approximation. The result is that the inflation field is rolling at a constant speed. The prediction for cosmological perturbation is calculated.

A brief history of the development of surface detectors for the study of the high-energy cosmic rays is presented. The paper is based on an invited talk given at UHECR2022 held in LAquila, October 2022. In a complementary talk, P Sokolsky discussed the development of the fluorescence technique for air-shower detection.

We show that the general vacuum states that respect the de Sitter symmetry, known as the {\alpha}-vacua, can introduce non-vanishing parity-violating tensor non-Gaussianities. This is due to the mixing by the Bogoliubov transformation of the positive and negative frequency modes of the Bunch-Davies vacuum. We calculate explicitly the bispectra of tensor perturbations and show that the amplitude can be exponentially enhanced for certain choices of the squeezing parameter {\alpha} and the phase {\phi} of the {\alpha}-vacua. We find a new shape for the parity-violating tensor bispectrum which peaks in the flattened configuration

This paper presents and investigates non-Gaussian perturbations for the warm k-inflation model that is driven by pure kinetic energy. The two complementary components of the overall non-Gaussianity are the three-point and four-point correlations. The intrinsic non-Gaussian component, denoted as the nonlinear parameter f^{int}_{NL}, is rooted in the three-point correlation for the inflation field. Meanwhile, the {\delta}N part non-Gaussianity, denoted as f^{{\delta}N}_{NL}, is the contribution attributed to the four-point correlation function of the inflation field. In this paper, the above two components in warm k-inflation are individually computed and analyzed. Then, comparisons and discussions between them are conducted, and the non-Gaussian theoretical results are compared with experimental observations to determine the range of model parameters within the allowable range of observation.

We consider a possible crystalline equation of state (EOS) that follows from nucleons that are chiral solitons with deep relativistic quark bound states. In this model conventional quark matter does not occur till a high density threshold at which the quark bound states in the nucleons get compressed and merge with the continuum. Once the barrier at this threshold is overcome, roughly when, $n_B \sim 1/fm^3$, we expect the crystalline nuclear matter to make a sudden transition into quark matter when the EOS becomes soft through a decompression. The sudden increase in density triggers a collapse releasing large amount of gravitational potential energy that can generate a spherical outburst (or kilonova) of ejected matter.

Dark photons could be produced resonantly by the oscillating axion field in the early universe. This resonant production mechanism has been used in various contexts, including dark photon dark matter and primordial magnetic field production. However, for this resonant production to work in an expanding universe, a large axion-dark photon coupling is required, which is not easy to realize in terms of model building and requires the introduction of many charged fermions and/or the complex clockwork mechanism. In this paper, we present a new scenario that efficiently produces dark photons from the axion with a much smaller coupling. This is possible by modifying the dynamics of axion and significantly delaying the onset of oscillations, as in the so-called trapped misalignment mechanism. As a specific example, we consider models in which dark photon production occurs efficiently despite the small axion-dark photon coupling by temporally trapping an axion in a wrong minimum and releasing it after the Hubble parameter becomes much smaller than the axion mass. In this scenario, it is expected that the polarization asymmetry of dark photons and gravitational waves generated from dark photons will be significantly reduced.