Papers for Monday, Jan 08 2024

We aim to investigate the nature of time-variable X-ray sources detected in the {\it XMM-Newton} serendipitous survey. The X-ray light curves of objects in the {\it XMM-Newton} serendipitous survey were searched for variability and coincident serendipitous sources observed by {\it Chandra} were also investigated. Subsequent infrared spectroscopy of the counterparts to the X-ray objects that were identified using UKIDSS was carried out using {\it ISAAC} on the VLT. We found that the object 4XMM~J182531.5--144036 detected in the {\it XMM-Newton} serendipitous survey in April 2008 was also detected by {\it Chandra} as CXOU~J182531.4--144036 in July 2004. Both observations reveal a hard X-ray source displaying a coherent X-ray pulsation at a period of 781~s. The source position is coincident with a $K=14$ mag infrared object whose spectrum exhibits strong HeI and Br$\gamma$ emission lines and an infrared excess above that of early B-type dwarf or giant stars. We conclude that 4XMM~J182531.5--144036 is a Be/X-ray binary pulsar exhibiting persistent X-ray emission and is likely in a long period, low eccentricity orbit, similar to X Per.

Milky Way globular clusters (GCs) display chemical enrichment in a phenomenon called multiple stellar populations (MSPs). While the enrichment mechanism is not fully understood, there is a correlation between a cluster's mass and the fraction of enriched stars found therein. However, present-day GC masses are often smaller than their masses at the time of formation due to dynamical mass loss. In this work, we explore the relationship between mass and MSPs using the stellar stream 300S. We present the chemical abundances of eight red giant branch member stars in 300S with high-resolution spectroscopy from Magellan/MIKE. We identify one enriched star characteristic of MSPs and no detectable metallicity dispersion, confirming that the progenitor of 300S was a globular cluster. The fraction of enriched stars (12.5\%) observed in our 300S stars is less than the 50\% of stars found enriched in Milky Way GCs of comparable present-day mass ($\sim10^{4.5}$\msun). We calculate the mass of 300S's progenitor and compare it to the initial masses of intact GCs, finding that 300S aligns well with the trend between the system mass at formation and enrichment. 300S's progenitor may straddle the critical mass threshold for the formation of MSPs and can therefore serve as a benchmark for the stellar enrichment process. Additionally, we identify a CH star, with high abundances of \textit{s}-process elements, probably accreted from a binary companion. The rarity of such binaries in intact GCs may imply stellar streams permit the survival of binaries that would otherwise be disrupted.

We present the first results from a JWST program studying the role played by powerful radio jets in the evolution of the most massive galaxies at the onset of Cosmic Noon. Using NIRSpec integral field spectroscopy, we detect 24 rest-frame optical emission lines from the $z=3.5892$ radio galaxy 4C+19.71. 4C+19.71 contains one of the most energetic radio jets known, making it perfect for testing radio-mode feedback on the interstellar medium (ISM) of a $M_{\star}\sim10^{11}\,\rm M_{\odot}$ galaxy. The rich spectrum enables line ratio diagnostics showing that the radiation from the active galactic nucleus (AGN) dominates the ionization of the entire ISM out to at least $25\,$kpc, the edge of the detection. Sub-kpc resolution reveals filamentary structures and emission blobs in the warm ionized ISM distributed on scales of $\sim5$ to $\sim20\,$kpc. A large fraction of the extended gaseous nebula is located near the systemic velocity. This nebula may thus be the patchy ISM which is illuminated by the AGN after the passage of the jet. A radiatively-driven outflow is observed within $\sim5\,$kpc from the nucleus. The inefficient coupling ($\lesssim 10^{-4}$) between this outflow and the quasar and the lack of extreme gas motions on galactic scales are inconsistent with other high-$z$ powerful quasars. Combining our data with ground-based studies, we conclude that only a minor fraction of the feedback processes is happening on $<25\,$kpc scales.

The Galactic Center Excess (GCE) remains an enduring mystery, with leading explanations being annihilating dark matter or an unresolved population of millisecond pulsars. Analyzing the morphology of the GCE provides critical clues to identify its exact origin. We investigate the robustness of the inferred GCE morphology against the effects of masking, an important step in the analysis where the gamma-ray emission from point sources and the galactic disk are excluded. Using different masks constructed from Fermi point source catalogs and a wavelet method, we find that the GCE morphology, particularly its ellipticity and cuspiness, is relatively independent of the choice of mask for energies above 2-3 GeV. The GCE morphology systematically favors an approximately spherical shape, as expected for dark matter annihilation. Compared to various stellar bulge profiles, a spherical dark matter annihilation profile better fits the data across different masks and galactic diffuse emission backgrounds, except for the stellar bulge profile from Coleman et al. (2020), which provides a similar fit to the data. Modeling the GCE with two components, one from dark matter annihilation and one tracing the Coleman Bulge, we find this two-component model outperforms any single component or combinations of dark matter annihilation and other stellar bulge profiles. Uncertainty remains about the exact fraction contributed by each component across different background models and masks. However, when the Coleman Bulge dominates, its corresponding spectrum lacks characteristics typically associated with millisecond pulsars, suggesting that it mostly models the emission from other sources instead of the GCE that is still present and spherically symmetric.

Insights from JWST observations suggest that AGN feedback evolved from a short-lived, high redshift phase in which radiatively cooled turbulence and/or momentum-conserving outflows stimulated vigorous early star formation (``positive'' feedback), to late, energy-conserving outflows that depleted halo gas reservoirs and quenched star formation. The transition between these two regimes occurred at $z\sim 6$, independently of galaxy mass, for simple assumptions about the outflows and star formation process. Observational predictions provide circumstantial evidence for the prevalence of massive black holes at the highest redshifts hitherto observed, and we discuss their origins.

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys.

Henize 2-10 is a dwarf starburst galaxy hosting a $\sim10^{6}~M_{\odot}$ black hole (BH) that is driving an ionized outflow and triggering star formation within the central $\sim100$ pc of the galaxy. Here we present ALMA continuum observations from 99 to 340 GHz, as well as spectral line observations of the molecules CO (1-0, 3-2), HCN (1-0, 3-2), and HCO$^{+}$ (1-0, 3-2), with a focus on the BH and its vicinity. Incorporating cm-wave radio measurements from the literature, we show that the spectral energy distribution of the BH is dominated by synchrotron emission from 1.4 to~340 GHz with a spectral index of $\alpha\approx-0.5$. We analyze the spectral line data and identify an elongated molecular gas structure around the BH with a velocity distinct from the surrounding regions. The physical extent of this molecular gas structure is $\approx130~{\rm pc}\times30$ pc and the molecular gas mass is $\sim10^{6}~M_{\odot}$. Despite an abundance of molecular gas in this general region, the position of the BH is significantly offset from the peak intensity, which may explain why the BH is radiating at a very low Eddington ratio. Our analysis of the spatially-resolved line ratio between CO J=3-2 and J=1-0 implies that the CO gas in the vicinity of the BH is highly excited, particularly at the interface between the BH outflow and the regions of triggered star formation. This suggests that the cold molecular gas is being shocked by the bipolar outflow from the BH, supporting the case for positive BH feedback.

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) is designed to detect and measure the redshifts of more than one million Ly$\alpha$ emitting galaxies (LAEs) between $1.88 < z < 3.52$. In addition to its cosmological measurements, these data enable studies of Ly$\alpha$ spectral profiles and the underlying radiative transfer. Using the roughly half a million LAEs in the HETDEX Data Release 3, we stack various subsets to obtain the typical Ly$\alpha$ profile for the $z \sim 2-3$ epoch and to understand their physical properties. We find clear absorption wings around Ly$\alpha$ emission, which extend $\sim 2000$ km $\mathrm{s}^{-1}$ both redward and blueward of the central line. Using far-UV spectra of nearby ($0.002 < z < 0.182$) LAEs in the CLASSY treasury and optical/near-IR spectra of $2.8 < z < 6.7$ LAEs in the MUSE-Wide survey, we observe absorption profiles in both redshift regimes. Dividing the sample by volume density shows that the troughs increase in higher density regions. This trend suggests that the depth of the absorption is dependent on the local density of objects near the LAE, a geometry that is similar to damped Lyman-$\alpha$ systems. Simple simulations of Ly$\alpha$ radiative transfer can produce similar troughs due to absorption of light from background sources by HI gas surrounding the LAEs.

We report on the detection of radio bursts from the Galactic bulge using the real-time transient detection and localization system, realfast. The pulses were detected commensally on the Karl G. Jansky Very Large Array during a survey of unidentified Fermi $\gamma$-ray sources. The bursts were localized to subarcsecond precision using realfast fast-sampled imaging. Follow-up observations with the Green Bank Telescope detected additional bursts from the same source. The bursts do not exhibit periodicity in a search up to periods of 480s, assuming a duty cycle of < 20%. The pulses are nearly 100% linearly polarized, show circular polarization up to 12%, have a steep radio spectral index of -2.7, and exhibit variable scattering on timescales of months. The arcsecond-level realfast localization links the source confidently with the Fermi $\gamma$-ray source and places it nearby (though not coincident with) an XMM-Newton X-ray source. Based on the source's overall properties, we discuss various options for the nature of this object and propose that it could be a young pulsar, magnetar, or a binary pulsar system.

Context. Coronal mass ejections (CMEs) are the main driver of solar wind disturbances near Earth. When directed towards us, the internal magnetic field of the CME can interact with the Earth`s magnetic field and cause geomagnetic storms. In order to better predict and avoid damage coming from such events, the optimized heliospheric model Icarus has been implemented. Aims. The impact of a CME at Earth is greatly affected by its internal magnetic field structure. The aim of this work is to enable modelling the evolution of the magnetic field configuration of the CME throughout its propagation in Icarus. The focus of the study is on the global magnetic structure of the CME and its evolution and interaction with the solar wind. Methods. The magnetized CME model that is implemented in Icarus is the Linear Force-Free Spheromak and is imported from EUHFORIA. Advanced techniques, such as grid stretching and AMR are applied. Different AMR levels are applied in order to obtain high resolution locally, where needed. The results of all the simulations are compared in detail and the wall-clock times of the simulations are provided. Results. The results from the performed simulations are analyzed. The arrival time is better approximated by the EUHFORIA simulation, with the CME shock arriving 1.6 and 1.09 hours later than in the AMR level 4 and 5 simulations, respectively. The profile features and variable strengths are best modelled by Icarus simulations with AMR level 4 and 5. Conclusions. The arrival time is closer to the observed time in the EUHFORIA simulation, but the profiles of the different variables show more features and details in the Icarus simulations. Considering the small difference in the modelled results, and the large difference in computational resources, the AMR level 4 simulation is considered to have performed the best.

We here present the results of an analysis of the optical spectroscopy of 42 globular cluster (GC) candidates in the nearby spiral galaxy M81 (3.61~Mpc). The spectra were obtained using the long-slit and MOS modes of the OSIRIS instrument at the 10.4~m Gran Telescopio Canarias (GTC) at a spectral resolution of $\sim$1000. We used the classical H$\beta$ vs [MgFe]$'$ index diagram to separate genuine old GCs from clusters younger than 3 Gyr. Of the 30 spectra with continuum signal-to-noise ratio $>10$, we confirm 17 objects to be classical GCs (age $>10$~Gyr, $-1.4<$[Fe/H]$<-$0.4), with the remaining 13 being intermediate-age clusters (1-7.5~Gyr). We combined age and metallicity data of other nearby spiral galaxies ($\lesssim18$~Mpc) obtained using similar methodology like the one we have used here to understand the origin of GCs in spiral galaxies in the cosmological context. We find that the metal-poor ([Fe/H]<$-$1) GCs continued to form up to 6~Gyr after the first GCs were formed, with all younger systems (age $<8$~Gyr) being metal-rich.

We present an overview of recent numerical advances in the theoretical characterization of massive binary black hole (MBBH) mergers in astrophysical environments. These systems are among the loudest sources of gravitational waves (GWs) in the universe and particularly promising candidates for multimessenger astronomy. Coincident detection of GWs and electromagnetic (EM) signals from merging MBBHs is at the frontier of contemporary astrophysics. One major challenge in observational efforts searching for these systems is the scarcity of strong predictions for EM signals arising before, during, and after merger. Therefore, a great effort in theoretical work to-date has been to characterize EM counterparts emerging from MBBHs concurrently to the GW signal, aiming to determine distinctive observational features that will guide and assist EM observations. To produce sharp EM predictions of MBBH mergers it is key to model the binary inspiral down to coalescence in a full general relativistic fashion by solving Einstein's field equations coupled with the magnetohydrodynamics equations that govern the evolution of the accreting plasma in strong-gravity. We review the general relativistic numerical investigations that have explored the astrophysical manifestations of MBBH mergers in different environments and focused on predicting potentially observable smoking-gun EM signatures that accompany the gravitational signal.

Observed protostellar outflows exhibit a variety of asymmetrical features, including remarkable unipolar outflows and bending outflows. Revealing the formation and early evolution of such asymmetrical protostellar outflows, especially the unipolar outflows, is essential for a better understanding of the star and planet formation because they can dramatically change the mass accretion and angular momentum transport to the protostars and protoplanetary disks. Here, we perform the three-dimensional non-ideal magnetohydrodynamics simulations to investigate the formation and early evolution of the asymmetrical protostellar outflows in magnetized turbulent isolated molecular cloud cores. We find, for the first time to our knowledge, that the unipolar outflow forms even in the single low-mass protostellar system. The results show that the unipolar outflow is driven in the weakly magnetized cloud cores with the dimensionless mass-to-flux ratios of $\mu=8$ and $16$. Furthermore, we find the $\textit{protostellar rocket effect}$ of the unipolar outflow, which is similar to the launch and propulsion of a rocket. The unipolar outflow ejects the protostellar system from the central dense region to the outer region of the parent cloud core, and the ram pressure caused by its ejection suppresses the driving of additional new outflows. In contrast, the bending bipolar outflow is driven in the moderately magnetized cloud core with $\mu=4$. The ratio of the magnetic to turbulent energies of a parent cloud core may play a key role in the formation of asymmetrical protostellar outflows.

This paper explores core-collapse supernovae as crucial targets for neutrino telescopes, addressing uncertainties in their simulation results. We comprehensively analyze eighteen modern simulations and discriminate among supernova models using realistic detectors and interactions. A significant correlation between the total neutrino energy and cumulative counts, driven by massive lepton neutrinos and oscillations, is identified, particularly noticeable with the DUNE detector. Bayesian techniques indicate strong potential for model differentiation during a Galactic supernova event, with HK excelling in distinguishing models based on equation of state, progenitor mass, and mixing scheme.

The intracluster light (ICL) fraction is a well-known indicator of the dynamical activity in intermediate-redshift clusters. Merging clusters in the redshift interval $0.18<z<0.56$ have a distinctive peak in the ICL fractions measured between $\sim 3800-4800$ \AA. In this work, we analyze two higher-redshift, clearly merging clusters, ACT-CLJ0102-49151 and CL J0152.7-1357, at $z>0.8$, using the HST optical and infrared images obtained by the RELICS survey. We report the presence of a similar peak in the ICL fractions, although wider and redshifted to the wavelength interval $\sim 5200-7300$ \AA. The fact that this excess in the ICL fractions is found at longer wavelengths can be explained by an assorted mixture of stellar populations in the ICL, direct inheritance of an ICL that was mainly formed by major galaxy mergers with the BCG at $z>1$ and whose production is instantaneously burst by the merging event. The ubiquity of the ICL fraction merging signature across cosmic time enhances the ICL as a highly reliable and powerful probe to determine the dynamical stage of galaxy clusters, which is crucial for cluster-based cosmological inferences that require relaxation of the sample.

More than 36 years have passed since the discovery of the infrared excess from circumstellar dust orbiting the white dwarf G29-38, which at 17.5 pc it is the nearest and brightest of its class. The precise morphology of the orbiting dust remains only marginally constrained by existing data, subject to model-dependent inferences, and thus fundamental questions of its dynamical origin and evolution persist. This study presents a means to constrain the geometric distribution of the emitting dust using stellar pulsations measured at optical wavelengths as a variable illumination source of the dust, which re-radiates primarily in the infrared. By combining optical photometry from the Whole Earth Telescope with 0.7-2.5 micron spectroscopy obtained with SpeX at NASA's Infrared Telescope Facility, we detect luminosity variations at all observed wavelengths, with variations at most wavelengths corresponding to the behavior of the pulsating stellar photosphere, but towards the longest wavelengths the light curves probe the corresponding time-variability of the circumstellar dust. In addition to developing methodology, we find pulsation amplitudes decrease with increasing wavelength for principal pulsation modes, yet increase beyond approximately 2 microns for nonlinear combination frequencies. We interpret these results as combination modes deriving from principal modes of identical l values and discuss the implications for the morphology of the warm dust. We also draw attention to some discrepancies between our findings and theoretical expectations for the results of the non-linearity imposed by the surface convection zone on mode--mode interactions and on the behavior of the first harmonic of the highest-amplitude pulsation mode.

Ly$\alpha$ emission is an exceptionally informative tracer of the life cycle of evolving galaxies and the escape of ionising photons. However, theoretical studies of Ly$\alpha$ emission are often limited by insufficient numerical resolution, incomplete sets of physical models, and poor line-of-sight (LOS) statistics. To overcome such limitations, we utilize here the novel PANDORA suite of high-resolution dwarf galaxy simulations that include a comprehensive set of state-of-the-art physical models for ionizing radiation, magnetic fields, supernova feedback and cosmic rays. We post-process the simulations with the radiative transfer code \textsc{RASCAS} to generate synthetic observations and compare to observed properties of Ly$\alpha$ emitters. Our simulated Ly$\alpha$ haloes are more extended than the spatial region from which the intrinsic emission emanates and our spatially resolved maps of spectral parameters of the Ly$\alpha$ emission are very sensitive to the underlying spatial distribution and kinematics of neutral hydrogen. Ly$\alpha$ and LyC emission display strongly varying signatures along different LOS depending on how each LOS intersects low-density channels generated by stellar feedback. Comparing galaxies simulated with different physics, we find the Ly$\alpha$ signatures to exhibit systematic offsets determined by the different levels of feedback strength and the clumpiness of the neutral gas. Despite this variance, and regardless of the different physics included in each model, we find universal correlations between Ly$\alpha$ observables and LyC escape fraction, demonstrating a robust connection between Ly$\alpha$ and LyC emission. Ly$\alpha$ observations from a large sample of dwarf galaxies should thus give strong constraints on their stellar feedback-regulated LyC escape and confirm their important role for the reionization of the Universe.

Stellar deformations play a significant role in the dynamical evolution of stars in binary systems, impacting the tidal dissipation and the outcomes of mass transfer processes. The prevalent method for modelling the deformations and tidal interactions of celestial bodies solely relies on the perturbative approach, which assumes that stellar deformations are minor perturbations to the spherical symmetry. An observable consequence of stellar deformations is the apsidal motion in eccentric systems. Our objective is to assert the reliability of the perturbative approach when applied to close and strongly deformed binary systems. We have developed a non-perturbative 3D modelling method designed to account for high stellar deformations to explore the limitations of the perturbative models. Our research highlights that the perturbative model becomes imprecise and underestimates the tidal force and rate of apsidal motion at a short orbital separation. This discrepancy primarily results from the first-order treatment in the perturbative approach, and cannot be rectified using straightforward mathematical corrections due to the strong non-linearity and numerous parameters of the problem. We have determined that our methodology affects the modelling of approximately 42% of observed binary systems with measured apsidal motion, introducing a discrepancy greater than 2% when the normalised orbital separation verifies q^(-1/5)a(1-e^2)/R1 < 6.5. The perturbative approach underestimates tidal interactions between bodies up to ~40% for close low-mass binaries. All the subsequent modelling is impacted by our findings, in particular, the tidal dissipation is significantly underestimated. As a result, all binary stellar models are imprecise when applied to systems with a low orbital separation, and the outcomes of these models are also affected by these inaccuracies.

The lifetime of mm size dust grains, such as chondrules, in the nominal solar nebula model is limited to $\sim 10^{5}$ yr due to an inward drift driven by gas drag. However, isotopic and petrological studies on primitive meteorites indicate a discrepancy of $\gtrsim 10^{6}$ yr between the formation time of chondrules and that of chondritic parent bodies. Therefore chondrules should survive for $\gtrsim 10^{6}$ yr in the solar nebula against the inward drift without subsequent growth (i.e., planetesimal formation). Here we investigate the conditions of the solar nebula that are suitable for the long lifetime of chondrule-sized dust particles. We take the turbulent strength, the radial pressure gradient force, and the disk metallicity of the solar nebula as free parameters. For 1 mm-radius-chondrules to survive and keep their size for $\gtrsim 10^{6}$ yr, the suitable condition is a weak turbulence ($\alpha \sim 10^{-6}$), a flat radial profile ($\eta \lesssim 10^{-3}$), and a high metallicity ($Z\sim 0.1$). This condition is qualitatively consistent with the characteristics of protoplanetary disks suggested by recent observations. We eventually propose that planetesimal formation may be induced by the disk evolution, e.g., the inside-out dispersal of the gas component due to the disk wind.

Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral scan survey from the ALMA Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1".6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH$^+$ and the first interferometric images of C$_3$H$^+$, NO, HCS$^+$. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our dataset is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as i) young starburst tracers characterized by high-excitation transitions of HC$_3$N and complex organic molecules (COMs) versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing PDRs, ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH$_3$OH, HNCO, HOCO$^+$, and OCS) versus shocks from star-formation-induced outflows (high-excitation transitions of SiO), as well as iii) outflows showing emission from HOC$^+$, CCH, H$_3$O$^+$, CO isotopologues, HCN, HCO$^+$, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback.

We present a study of high-resolution spectra of RS Canum Venaticorum (RS CVn), a prototype of active binary systems. Our data were obtained from 1998 to 2017 using different telescopes. We analyze the chromospheric activity indicators Ca II IRT, H$_{\alpha}$, Na I D$_{1}$, D$_{2}$ doublet, He I D$_{3}$, and H$_{\beta}$ using a spectral subtraction technique. The chromospheric emission stems mainly from the K2 IV primary star, while the F5 V secondary star only shows weak emission features in a few of our spectra. We find excess absorption features in the subtracted H$_{\alpha}$ lines and other activity indicators from spectra taken near primary eclipse, which we ascribe to prominence-like material associated with the primary star. We estimate size limits of these tentative prominences based on the geometry of the binary system, and investigate the physical properties of the strongest prominence. An optical flare, characterized by He I D$_{3}$ line emission, together with stronger emission in other activity lines, was detected. The flare energy is roughly comparable to strong flares observed on other RS CVn-type stars. The chromospherically active longitudes of RS CVn most frequently appear near the two quadratures of the system and display changes between observing runs, which indicates an ongoing evolution of its active regions.

We present the study on continuous high-resolution spectroscopic observations of two long-period single-lined RS Canum Venaticorum (RS CVn) binary stars IM Pegasi (IM Peg) and $\sigma$ Geminorum ($\sigma$ Gem), obtained with the Hertzsprung SONG telescope during the 2015-2016 season. Chromospheric activity indicators H$_{\alpha}$, Na I D$_{1}$, D$_{2}$ doublet, He I D$_{3}$, and H$_{\beta}$ lines have been analyzed by using the spectral subtraction technique. The expected chromospheric emission features in the H$_{\alpha}$, Na I D$_{1}$, D$_{2}$ doublet, and H$_{\beta}$ lines confirm that both of two stars are very active systems. In the spectra, the He I D$_{3}$ line had been always detected in absorption feature. Although the behavior of chromospheric activity indicators is very similar for both stars, the activity level of IM Peg is much stronger than that of $\sigma$ Gem. Moreover, the EW variations of the H$_{\alpha}$, He I D$_{3}$, and H$_{\beta}$ line subtractions correlate well and show different behavior among different orbital cycles, which indicates the presence and evolution of activity longitudes over the surface of two stars. Furthermore, the subtracted H$_{\alpha}$ line profile is usually asymmetric. The red-shifted excess absorption features could be interpreted as a strong down-flow of cool absorbing material, while the blue-shifted emission component is probably caused by up-flow of hot materials through microflare events.

In the forthcoming era of big astronomical data, it is a burden to find out target sources from ground-based and space-based telescopes. Although Machine Learning (ML) methods have been extensively utilized to address this issue, the incorporation of in-depth data analysis can significantly enhance the efficiency of identifying target sources when dealing with massive volumes of astronomical data. In this work, we focused on the task of finding AGN candidates and identifying BL Lac/FSRQ candidates from the 4FGL DR3 uncertain sources. We studied the correlations among the attributes of the 4FGL DR3 catalogue and proposed a novel method, named FDIDWT, to transform the original data. The transformed dataset is characterized as low-dimensional and feature-highlighted, with the estimation of correlation features by Fractal Dimension (FD) theory and the multi-resolution analysis by Inverse Discrete Wavelet Transform (IDWT). Combining the FDIDWT method with an improved lightweight MatchboxConv1D model, we accomplished two missions: (1) to distinguish the Active Galactic Nuclei (AGNs) from others (Non-AGNs) in the 4FGL DR3 uncertain sources with an accuracy of 96.65%, namely, Mission A; (2) to classify blazar candidates of uncertain type (BCUs) into BL Lacertae objects (BL Lacs) or Flat Spectrum Radio Quasars (FSRQs) with an accuracy of 92.03%, namely, Mission B. There are 1354 AGN candidates in Mission A, 482 BL Lacs candidates and 128 FSRQ candidates in Mission B were found. The results show a high consistency of greater than 98% with the results in previous works. In addition, our method has the advantage of finding less variable and relatively faint sources than ordinary methods.

It is currently unknown how matter behaves at the extreme densities found within the cores of neutron stars. Measurements of the neutron star equation of state probe nuclear physics that is otherwise inaccessible in a laboratory setting. Gravitational waves from binary neutron star mergers encode details about this physics, allowing the equation of state to be inferred. Planned third-generation gravitational-wave observatories, having vastly improved sensitivity, are expected to provide tight constraints on the neutron star equation of state. We combine simulated observations of binary neutron star mergers by the third-generation observatories Cosmic Explorer and Einstein Telescope to determine future constraints on the equation of state across a plausible neutron star mass range. In one year of operation, a network consisting of one Cosmic Explorer and the Einstein Telescope is expected to detect $\gtrsim 3\times 10^5$ binary neutron star mergers. By considering only the 75 loudest events, we show that such a network will be able to constrain the neutron star radius to at least $\lesssim 200$ m (90% credibility) in the mass range $1-1.97$ $M_{\odot}$ -- about ten times better than current constraints from LIGO-Virgo-KAGRA and NICER. The constraint is $\lesssim 75$ m (90% credibility) near $1.4-1.6$ $M_{\odot}$ where we assume the the binary neutron star mass distribution is peaked. This constraint is driven primarily from the loudest $\sim 20$ events.

Determining the size distribution of asteroids is key for understanding the collisional history and evolution of the inner Solar System. We aim at improving our knowledge on the size distribution of small asteroids in the Main Belt by determining the parallaxes of newly detected asteroids in the Hubble Space Telescope (HST) Archive and hence their absolute magnitudes and sizes. Asteroids appear as curved trails in HST images due to the parallax induced by the fast orbital motion of the spacecraft. The parallax effect can be computed to obtain the distance to the asteroids by fitting simulated trajectories to the observed trails. Using distance, we can obtain the object's absolute magnitude and size estimation assuming an albedo value, along with some boundaries for its orbital parameters. In this work we analyse a set of 632 serendipitously imaged asteroids found in the ESA HST Archive. An object-detection machine learning algorithm was used to perform this task during previous work. Our raw data consists of 1,031 asteroids trails from unknown objects (not matching any entries in the MPC database). We also found 670 trails from known objects (objects featuring matching entries in the MPC). After an accuracy assessment and filtering process, our analysed HST set consists of 454 unknown objects and 178 known objects. We obtain a sample dominated by potential Main Belt objects featuring absolute magnitudes (H) mostly between 15 and 22 mag. The absolute magnitude cumulative distribution confirms the previously reported slope change for 15 < H < 18, from 0.56 to 0.26, maintained in our case down to absolute magnitudes around H = 20, hence expanding the previous results by approximately two magnitudes. HST archival observations can be used as an asteroid survey since the telescope pointings are statistically randomly oriented in the sky and they cover long periods of time.

Under the hypothesis of gravitational redshift induced by the central supermassive black hole, and based on line widths and shifts of redward shifted H$\beta$ and H$\alpha$ broad emission lines for more than 8000 SDSS DR7 AGNs, we measure the virial factor in determining supermassive black hole masses. The virial factor had been believed to be independent of accretion radiation pressure on gas clouds in broad-line region (BLR), and only dependent on inclination effects of BLR. The virial factor measured spans a very large range. For the vast majority of AGNs ($>$96%) in our samples, the virial factor is larger than $f=1$ usually used in literatures. The $f$ correction makes the percent of high-accreting AGNs decrease by about 100 times. There are positive correlations of $f$ with the dimensionless accretion rate and Eddington ratio. The redward shifts of H$\beta$ and H$\alpha$ are mainly the gravitational origin, confirmed by a negative correlation between the redward shift and the dimensionless radius of BLR. Our results show that radiation pressure force is a significant contributor to the measured virial factor, containing the inclination effects of BLR. The usually used values of $f$ should be corrected for high-accreting AGNs, especially high redshift quasars. The $f$ correction increases their masses by one--two orders of magnitude, which will make it more challenging to explain the formation and growth of supermassive black holes at high redshifts.

The full-disk chromosphere was routinely monitored in the He I 10830\AA\, line at the National Solar Observatory/Kitt Peak from 2004 Nov. to 2013 March, and thereby, synoptic maps of He I line intensity from Carrington rotations 2032 to 2135 were acquired. They are utilized to investigate the differential rotation of the chromosphere and the quiet chromosphere during the one falling (descending part of solar cycle 23) and the one rising (ascending part of solar cycle 24) period of a solar cycle. Both the quiet chromosphere and the chromosphere are found to rotate slower and have a more prominent differential rotation, in the rising period of solar cycle 24 than in the falling period of solar cycle 23, and an illustration is offered.

The corona is a structure possessed by stars, including the sun. The abnormal heating of the solar corona and chromosphere is one of the greatest mysteries in modern astronomy. While state-of-the-art observations have identified some candidates of magnetic activity events that could be responsible for this abnormal heating, and theoretical studies have proposed various heating modes, a complete physical picture of how they are heated as a whole remains elusive. In this study, the characteristics of the heated corona and chromosphere are investigated, and for the first time, the question of how they are abnormally heated is explicitly answered by analyzing the long-term observations of the global chromosphere in the Ca II K line and the global corona in the coronal green line. The findings reveal that both the quiet chromosphere and corona are in anti-phase with the solar cycle, whereas the active chromosphere and corona are in phase with it. Different parts of the solar corona and chromosphere exhibit significantly different variation characteristics, and are found to be heated by different magnetic categories and probably in different modes. This study posits that unraveling the heating mystery is best approached through the lens of magnetic categories, rather than magnetic activity events.

Many sub-Neptune exoplanets have been believed to be composed of a thick hydrogen-dominated atmosphere and a high-temperature heavier-element-dominant core. From an assumption that there is no chemical reaction between hydrogen and silicates/metals at the atmosphere-interior boundary, the cores of sub-Neptunes have been modeled with molten silicates and metals (magma) in previous studies. In large sub-Neptunes, pressure at the atmosphere-magma boundary can reach tens of gigapascals where hydrogen is a dense liquid. A recent experiment showed that hydrogen can induce the reduction of Fe$^{2+}$ in (Mg,Fe)O to Fe$^0$ metal at the pressure-temperature conditions relevant to the atmosphere-interior boundary. However, it is unclear if Mg, one of the abundant heavy elements in the planetary interiors, remains oxidized or can be reduced by H. Our experiments in the laser-heated diamond-anvil cell found that heating of MgO + Fe to 3500-4900 K (close to or above their melting temperatures) in a H medium leads to the formation of Mg$_2$FeH$_6$ and H$_2$O at 8-13 GPa. At 26-29 GPa, the behavior of the system changes, and Mg-H in an H fluid and H$_2$O were detected with separate FeH$_x$. The observations indicate the dissociation of the Mg-O bond by H and subsequent production of hydride and water. Therefore, the atmosphere-magma interaction can lead to a fundamentally different mineralogy for sub-Neptune exoplanets compared with rocky planets. The change in the chemical reaction at the higher pressures can also affect the size demographics (i.e., "radius cliff") and the atmosphere chemistry of sub-Neptune exoplanets.

Atmospheric photochemistry on Titan continuously transforms methane and nitrogen gases into various organic compounds. This study explores the fate of these molecules when they land on Titan's surface. Our analytical exploration reveals that most simple organics found in Titan's atmosphere, including all nitriles, triple-bonded hydrocarbons, and benzene, land as solids. Only a few compounds are in the liquid phase, while only ethylene remains gaseous. For the simple organics that land as solids, we further examine their interactions with Titan's lake liquids. Utilizing principles of buoyancy, we found that flotation can be achieved via porosity-induced (25-60% porosity) or capillary force-induced buoyancy for HCN ices on ethane-rich lakes. Otherwise, these ices would sink and become lakebed sediments. By evaluating the timescale of flotation, our findings suggest that porosity-induced flotation of millimeter-sized and larger sediments is the only plausible mechanism for floating solids to explain the transient "magic islands" phenomena on Titan's lakes.

30$\%$--50$\%$ of the luminous pulsating red-giant stars show light variations of a longer period than the pulsation periods. Those periods are called long secondary periods (LSP). There has been debated for many years but the origin of the LSP is still unknown. To explain the LSP variations, there have been many approaches in not only observations but also theoretical studies. However, the invariance of the length of the LSPs has been investigated little. Thus, we studied the temporal variations of the period by performing the weighted wavelet-Z-transform analysis. Using the OGLE-III database, the $I$-band light curves of 6904 and 1945 LSP candidates in the Large/Small Magellanic Clouds, respectively, were analyzed. Most of our sample stars indicated that the period corresponding to the LSP was constant during the observation term. However, 101 and 44 LSP stars in the LMC and SMC, respectively, showed the signature of the temporal variations of the LSPs. There were diversities of the period modulations i.e. monotonic increase or decrease, or constant until the middle and then increase, etc. The comet-like companion is one of the possible explanations for the LSP variations, but this hypothesis cannot explain the period modulations because the LSP is determined by the orbital period.

Context. In September 2022, the transient neutron star low-mass X-ray binary XTE J1701-462 went into a new outburst. Aims. The objective of this work is to examine the evolution of the accretion geometry of XTE J1701-462 by studying the spectro-polarimetric properties along the Z track of this source. The simultaneous observations archived by the Insight-Hard X-ray Modulation Telescope (HXMT) and the Imaging X-ray Polarimetry Explorer (IXPE) give us the opportunity. Methods. We present a comprehensive X-ray spectro-polarimetric analysis of XTE J1701-462, using simultaneous observations from IXPE, Insight-HXMT and NuSTAR. For IXPE observations, two methods are employed to measure the polarization: a model-independent measurement with PCUBE and a model-dependent polarization-spectral analysis with XSPEC. The corresponding spectra from Insight-HXMT and NuSTAR are studied with two configurations that correspond to a slab-like corona and a spherical shell-like corona, respectively. Results. Significant polarization characteristics are detected in XTE J1701-462. The polarization degree shows a decreasing trend along the Z track, reducing from (4.84 $\pm$ 0.37)% to (3.76 $\pm$ 0.43)% on the horizontal branch and jumping to less than 1% on the normal branch. The simultaneous spectral analysis from Insight-HXMT and NuSTAR suggests that the redistribution between the thermal and Comptonized emission could be the reason for the PD evolution along the Z track. Based on the correlated spectro-polarimetric properties, we propose that this source likely has a slab coronal geometry and the size/thickness of the corona decreases along the Z track.

The colliding-wind region in binary systems made of massive stars allows us to investigate various aspects of shock physics, including particle acceleration. Particle accelerators of this kind are tagged as Particle-Accelerating Colliding-Wind Binaries, and are mainly identified thanks to their synchrotron radio emission. Our objective is first to validate the idea that obtaining snapshot high-resolution radio images of massive binaries constitutes a relevant approach to unambiguously identify particle accelerators. Second, we intend to exploit these images to characterize the synchrotron emission of two specific targets, HD167971 and HD168112, known as particle accelerators. We traced the radio emission from the two targets at 1.6 GHz with the European Very Long Baseline Interferometry Network, with an angular resolution of a few milli-arcseconds. Our measurements allowed us to obtain images for both targets. For HD167971, our observation occurs close to apastron, at an orbital phase where the synchrotron emission is minimum. For HD168112, we resolved for the very first time the synchrotron emission region. The emission region appears slightly elongated, in agreement with expectation for a colliding-wind region. In both cases the measured emission is significantly stronger than the expected thermal emission from the stellar winds, lending strong support for a non-thermal nature. Our study brings a significant contribution to the still poorly addressed question of high angular resolution radio imaging of colliding-wind binaries. We show that snapshot Very Long Baseline Interferometry measurements constitute an efficient approach to investigate these objects, with promising results in terms of identification of additional particle accelerators, on top of being promising as well to reveal long period binaries.

The Enhanced X-ray Timing and Polarimetry (eXTP) mission is a space mission to be launched in the late 2020s that is currently in development led by China in international collaboration with European partners. Here we provide a progress report on the Czech contribution to the eXTP science. We report on our simulation results performed in Opava (Institute of Physics of the Silesian University in Opava) and Prague (Astronomical Institute of the Czech Academy of Sciences), where the advanced timing capabilities of the satellite have been assessed for bright X-ray binaries that contain an accreting neutron star (NS) and exhibit the quasi-periodic oscillations. Measurements of X-ray variability originating in oscillations of fluid in the innermost parts of the accretion region determined by general relativity, such as the radial or Lense-Thirring precession, can serve for sensitive tests enabling us to distinguish between the signatures of different viable dense matter equations of state. We have developed formulae describing non-geodesic oscillations of accreted fluid and their simplified practical forms that allow for an expeditious application of the universal relations determining the NS properties. These relations, along with our software tools for studying the propagation of light in strong gravity and neutron star models, can be used for precise modeling of the X-ray variability while focusing on properties of the intended Large Area Detector (LAD). We update the status of our program and set up an electronic repository that will provide simulation results and gradual updates as the mission specifications progress toward their final formulation.

We analyze JWST NIRSpec$+$MIRI/MRS observations of the infrared (IR) gas-phase molecular bands of the most enshrouded source (D1) within the interacting system and luminous IR galaxy II Zw 096. We report the detection of rovibrational lines of H$_2$O $\nu_2$=1-0 ($\sim$5.3-7.2 $\mu$m) and $^{12}$CO $\nu$=1-0 ($\sim$4.45-4.95 $\mu$m) in D1. The CO band shows the R- and P-branches in emission and the spectrum of the H$_2$O band shows the P-branch in emission and the R-branch in absorption. The H$_2$O R-branch in absorption unveils an IR-bright embedded compact source in D1 and the CO broad component features a highly turbulent environment. From both bands, we also identified extended intense star-forming (SF) activity associated with circumnuclear photodissociation regions (PDRs), consistent with the strong emission of the ionised 7.7 $\mu$m polycyclic aromatic hydrocarbon band in this source. By including the 4.5-7.0 $\mu$m continuum information derived from the H$_2$O and CO analysis, we modelled the IR emission of D1 with a dusty torus and SF component. The torus is very compact (diameter of $\sim$3 pc at 5 $\mu$m) and characterised by warm dust ($\sim$ 370 K), giving an IR surface brightness of $\sim$3.6$\times$10$^{8}$ L$_{\rm sun}$/pc$^2$. This result suggests the presence of a dust-obscured active galactic nucleus (AGN) in D1, which has an exceptionally high covering factor that prevents the direct detection of AGN emission. Our results open a new way to investigate the physical conditions of inner dusty tori via modelling the observed IR molecular bands.

(Abridged). We explore the fraction of radio loud quasars in the eHAQ+GAIA23 sample, which contains quasars from the High A(V) Quasar (HAQ) Survey, the Extended High A(V) Quasar (eHAQ) Survey, and the Gaia quasar survey. All quasars in this sample have been found using a near-infrared color selection of target candidates that have otherwise been missed by the Sloan Digital Sky Survey (SDSS). We implemented a redshift-dependent color cut in g-i to select red quasars in the sample and divided them into redshift bins, while using a nearest-neighbors algorithm to control for luminosity and redshift differences between our red quasar sample and a selected blue sample from the SDSS. Within each bin, we cross-matched the quasars to the Faint Images of the Radio Sky at Twenty centimeters (FIRST) survey and determined the radio-detection fraction. We find similar radio-detection fractions for red and blue quasars within 1 sigma, independent of redshift. This disagrees with what has been found in the literature for red quasars in SDSS. It should be noted that the fraction of broad absorption line (BAL) quasars in red SDSS quasars is about five times lower than in our sample. BAL quasars have been observed to be more frequently radio quiet than other quasars, therefore the difference in BAL fractions could explain the difference in radio-detection fraction. The observed higher proportion of BAL quasars in our dataset relative to the SDSS sample, along with the higher rate of radio detections, indicates an association of the redness of quasars and the inherent BAL fraction within the overall quasar population. This finding highlights the need to explore the underlying factors contributing to both the redness and the frequency of BAL quasars, as they appear to be interconnected phenomena.

We aim to understand the formation of cloud condensation nuclei in oxygen-rich substellar atmospheres by calculating fundamental properties of the energetically most favorable vanadium oxide molecules and clusters. A hierarchical optimization approach is applied in order to find the most favorable structures for clusters of (VO)$_{N}$ and (VO$_2$)$_{N}$ for N=1-10, and (V$_2$O$_5$)$_{N}$ for N=1-4 and to calculate their thermodynamical potentials. The candidate geometries are initially optimized applying classical interatomic potentials and then refined at the B3LYP/cc-pVTZ level of theory to obtain accurate zero-point energies and thermochemical quantities. We present previously unreported vanadium oxide cluster structures as lowest-energy isomers. We report revised cluster energies and their thermochemical properties. Chemical equilibrium calculations are used to asses the impact of the updated and newly derived thermodynamic potentials on the gas-phase abundances of vanadium-bearing species. In chemical equilibrium, larger clusters from different stoichiometric families are found to be the most abundant vanadium-bearing species for temperatures below ~1000 K, while molecular VO is the most abundant between ~1000 K and ~2000 K. We determine the nucleation rates of each stoichiometric family for a given (T$_{gas}$, p$_{gas}$) profile of a brown dwarf using classical and non-classical nucleation theory. Small differences in the revised Gibbs free energies of the clusters have a large impact on the abundances of vanadium bearing species in chemical equilibrium at temperatures below ~1000 K, which subsequently has an impact on the nucleation rates of each stoichiometric family. We find that with the revised and more accurate cluster data non-classical nucleation rates are up to 15 orders of magnitude higher than classical nucleation rates.

Relic galaxies are massive, compact, quiescent objects observed in the local Universe that have not experienced any significant interaction episodes or merger events since about $z = 2$, remaining relatively unaltered since their formation. On the other hand, massive and compact Early Type Galaxies (cETGs) in the local Universe appear to show similar properties to Relic galaxies, despite having substantial accretion history. Relic galaxies, with frozen history, can provide important clues about the intrinsic processes related to the evolutionary pathways of ETGs and the role that mergers play in their evolution. Using the high-resolution cosmological simulation TNG50-1 from the Illustris Project, we investigate the assembly history of a sample of massive, compact, old, and quiescent subhalos split by satellite accretion fraction. We compare the evolutionary pathways at three cosmic epochs: $z = 2$, $z = 1.5$, and $z = 0$, using the orbital decomposition numerical method to investigate the stellar dynamics of each galactic kinematical component and their environmental correlations. Our results point to a steady pathway across time that is not strongly dependent on the environment. Relics and cETGs do not show a clear preference for high or low-density environments within the volume explored at TNG50. However, progenitors of Relic galaxies are shown to be located in high density since $z = 2$. The merger history can be recovered from the hot inner stellar halo imprints in the local Universe. In the current scenario, the mergers that drive the growth of cETGs do not give rise to a new and distinct evolutionary pathway when compared to Relics. This is despite the reported effects on the age and metallicity of the kinematic components.

The diffusion of water molecules through mesoporous dust of amorphous carbon (a-C) is a key process in the evolution of prestellar, protostellar, and protoplanetary dust, as well as in that of comets. It also plays a role in the formation of planets. Given the absence of data on this process, we experimentally studied the isothermal diffusion of water molecules desorbing from water ice buried at the bottom of a mesoporous layer of aggregated a-C nanoparticles, a material analogous to protostellar and cometary dust. We used infrared spectroscopy to monitor diffusion in low temperature (160 to 170 K) and pressure (6 $\times$ 10$^{-5}$ to 8 $\times$ 10$^{-4}$ Pa) conditions. Fick's first law of diffusion allowed us to derive diffusivity values on the order of 10$^{-2}$ cm$^2$ s$^{-1}$, which we linked to Knudsen diffusion. Water vapor molecular fluxes ranged from 5 $\times$ 10$^{12}$ to 3 $\times$ 10$^{14}$ cm$^{-2}$ s$^{-1}$ for thicknesses of the ice-free porous layer ranging from 60 to 1900 nm. Assimilating the layers of nanoparticles to assemblies of spheres, we attributed to this cosmic dust analog of porosity 0.80-0.90 a geometry correction factor, similar to the tortuosity factor of tubular pore systems, between 0.94 and 2.85. Applying the method to ices and refractory particles of other compositions will provides us with other useful data.

MWC 758 is a young star hosting a spiral protoplanetary disk. The spirals are likely companion-driven, and two previously-identified candidate companions have been identified -- one at the end the Southern spiral arm at ~0.6 arcsec, and one interior to the gap at ~0.1 arcsec. With JWST/NIRCam, we provide new images of the disk and constraints on planets exterior to ~1". We detect the two-armed spiral disk, a known background star, and a spatially resolved background galaxy, but no clear companions. The candidates that have been reported are at separations that are not probed by our data with sensitivity sufficient to detect them -- nevertheless, these observations place new limits on companions down to ~2 Jupiter-masses at ~150 au and ~0.5 Jupiter masses at ~600 au. Owing to the unprecedented sensitivity of JWST and youth of the target, these are among the deepest mass-detection limits yet obtained through direct imaging observations, and provide new insights into the system's dynamical nature.

We present JWST/NIRCam F187N, F200W, F405N and F410M direct imaging data of the disk surrounding SAO 206462. Previous images show a very structured disk, with a pair of spiral arms thought to be launched by one or more external perturbers. The spiral features are visible in three of the four filters, with the non-detection in F410M due to the large detector saturation radius. We detect with a signal-to-noise ratio of 4.4 a companion candidate (CC1) that, if on a coplanar circular orbit, would orbit SAO 206462 at a separation of $\sim300$ au, $2.25\sigma$ away from the predicted separation for the driver of the eastern spiral. According to the BEX models, CC1 has a mass of $M_\mathrm{CC1}=0.8\pm0.3~M_\mathrm{J}$. No other companion candidates were detected. At the location predicted by simulations of both spirals generated by a single massive companion, the NIRCam data exclude objects more massive than $\sim2.2~M_\mathrm{J}$ assuming the BEX evolutionary models. In terms of temperatures, the data are sensitive to objects with $T_{\text{eff}}\sim650-850$ K, when assuming planets emit like blackbodies ($R_\mathrm{p}$ between 1 and $3 R_\mathrm{J}$). From these results, we conclude that if the spirals are driven by gas giants, these must be either cold or embedded in circumplanetary material. In addition, the NIRCam data provide tight constraints on ongoing accretion processes. In the low extinction scenario we are sensitive to mass accretion rates of the order $\dot{M}\sim10^{-9} M_\mathrm{J}$ yr$^{-1}$. Thanks to the longer wavelengths used to search for emission lines, we reach unprecedented sensitivities to processes with $\dot{M}\sim10^{-7} M_\mathrm{J}$ yr$^{-1}$ even towards highly extincted environments ($A_\mathrm{V}\approx50$~mag).

Stellar population synthesis (SPS) is essential for understanding galaxy formation and evolution. However, the recent discovery of rotation-driven phenomena in star clusters warrants a review of uncertainties in SPS models caused by overlooked factors, including stellar rotation. In this study, we investigate the impact of rotation on SPS specifically using the PARSEC V2.0 rotation model and its implications for high redshift galaxies with the JWST. Rotation enhances the ultraviolet (UV) flux for up to $\sim 400$ Myr after the starburst, with the slope of UV increasing as the population gets faster rotating and more metal-poor. Using the Prospector tool, we construct simulated galaxies and deduce their properties associated with dust and star formation. Our results suggest that rapid rotation models result in a gradual UV slope up to 0.1 dex higher and an approximately 50\% increase in dust attenuation for identical wide-band spectral energy distributions. Furthermore, we investigate biases if the stellar population should be characterized by rapid rotation and demonstrate that accurate estimation can be achieved for rotation rates up to $\omega_\text{i}=0.6$. Accounting for the bias in the case of rapid rotation aligns specific star formation rates more closely with predictions from theoretical models. Notably, this also implies a slightly higher level of dust attenuation than previously anticipated, while still allowing for a `dust-free' interpretation of the galaxy. The impact of rapid rotation SPS models on the rest-UV luminosity function is found to be minimal. Overall, our findings have potentially important implications for comprehending dust attenuation and mass assembly history in the high-redshift Universe.

Many mechanisms have been proposed to alleviate the magnetic catastrophe, which prevents the Keplerian disk from forming inside a collapsing magnetized core. Such propositions include inclined field and non-ideal magnetohydrodynamics effects, and have been supported with numerical experiments. Models have been formulated for typical disk sizes when a field threads the rotating disk, parallel to the rotation axis, while observations at the core scales do not seem to show evident correlation between the directions of angular momentum and the magnetic field. In the present study, we propose a new model that considers both vertical and horizontal fields and discuss their effects on the protoplanetary disk size.

We derive a new criterion for estimating characteristic dynamical timescales in N-body simulations. The criterion uses the second, third, and fourth derivatives of particle positions: acceleration, jerk, and snap. It can be used for choosing timesteps in integrators with adaptive step size control. For any two-body problem the criterion is guaranteed to determine the orbital period and pericenter timescale regardless of eccentricity. We discuss why our criterion is the simplest derivative-based expression for choosing adaptive timesteps with the above properties and show its superior performance over existing criteria in numerical tests. Because our criterion uses lower order derivatives, it is less susceptible to rounding errors caused by finite floating point precision. This significantly decreases the volume of phase space where an adaptive integrator fails or gets stuck due to unphysical timestep estimates. For example, our new criterion can accurately estimate timesteps for orbits around a 50m sized Solar System object located at 40AU from the coordinate origin when using double floating point precision. Previous methods where limited to objects larger than 10km. We implement our new criterion in the high order IAS15 integrator which is part of the freely available N-body package REBOUND.

Using a 12 ks archival Chandra X-ray Observatory ACIS-S observation on the massive globular cluster (GC) M14, we detect a total of 7 faint X-ray sources within its half-light radius at a 0.5-7 keV depth of $2.5\times 10^{31}\,\mathrm{erg~s^{-1}}$. We cross-match the X-ray source positions with a catalogue of the Very Large Array radio point sources and a Hubble Space Telescope (HST) UV/optical/near-IR photometry catalogue, revealing radio counterparts to 2 and HST counterparts to 6 of the X-ray sources. In addition, we also identify a radio source with the recently discovered millisecond pulsar PSR 1737-0314A. The brightest X-ray source, CX1, appears to be consistent with the nominal position of the classic nova Ophiuchi 1938 (Oph 1938), and both Oph 1938 and CX1 are consistent with a UV-bright variable HST counterpart, which we argue to be the source of the nova eruption in 1938. This makes Oph 1938 the second classic nova recovered in a Galactic GC since Nova T Scorpii in M80. CX2 is consistent with the steep-spectrum radio source VLA8, which unambiguously matches a faint blue source; the steepness of VLA8 is suggestive of a pulsar nature, possibly a transitional millisecond pulsar with a late K dwarf companion, though an active galactic nucleus (AGN) cannot be ruled out. The other counterparts to the X-ray sources are all suggestive of chromospherically active binaries or background AGNs, so their nature requires further membership information.

We present an analysis of the X-ray properties 10 luminous, dust-reddened quasars from the FIRST-2MASS (F2M) survey based on new and archival Chandra observations. These systems are interpreted to be young, transitional objects predicted by merger-driven models of quasar/galaxy co-evolution. The sources have been well-studied from the optical through mid-infrared, have Eddington ratios above 0.1, and possess high-resolution imaging, most of which shows disturbed morphologies indicative of a recent or ongoing merger. When combined with previous X-ray studies of five other F2M red quasars, we find that the sources, especially those hosted by mergers, have moderate to high column densities ($N_H \simeq 10^{22.5-23.5}$ cm$^{-2}$) and Eddington ratios high enough to enable radiation pressure to blow out the obscuring material. We confirm previous findings that red quasars have dust-to-gas ratios that are significantly lower than the value for the Milky Way's interstellar medium, especially when hosted by a merger. The dust-to-gas ratio for two red quasars that lack evidence for merging morphology is consistent with the Milky Way and they do not meet the radiative feedback conditions for blowout. These findings support the picture of quasar/galaxy co-evolution in which a merger results in feeding of and feedback from an AGN. We compare the F2M red quasars to other obscured and reddened quasar populations in the literature, finding that, although morphological information is lacking, nearly all such samples meet blowout conditions and exhibit outflow signatures suggestive of winds and feedback.

Dust attenuation in star-forming galaxies (SFGs), as parameterized by the infrared excess (IRX $\equiv L_{\rm IR}/L_{\rm UV}$), is found to be tightly correlated with star formation rate (SFR), metallicity and galaxy size, following a universal IRX relation up to $z=3$. This scaling relation can provide a fundamental constraint for theoretical models to reconcile galaxy star formation, chemical enrichment, and structural evolution across cosmic time. We attempt to reproduce the universal IRX relation over $0.1\leq z\leq 2.5$ using the EAGLE hydrodynamical simulations and examine sensitive parameters in determining galaxy dust attenuation. Our findings show that while the predicted universal IRX relation from EAGLE approximately aligns with observations at $z\leq 0.5$, noticeable disparities arise at different stellar masses and higher redshifts. Specifically, we investigate how modifying various galaxy parameters can affect the predicted universal IRX relation in comparison to the observed data. We demonstrate that the simulated gas-phase metallicity is the critical quantity for the shape of the predicted universal IRX relation. We find that the influence of the infrared luminosity and infrared excess is less important while galaxy size has virtually no significant effect. Overall, the EAGLE simulations are not able to replicate some of the observed characteristics between IRX and galaxy parameters of SFGs, emphasizing the need for further investigation and testing for our current state-of-the-art theoretical models.

Isolated galaxies are the ideal reference sample to study the galaxy structure minimising potential environmental effects. We selected a complete sample of 14 nearby, late-type, highly inclined ($i\geq80^{\circ}$), isolated galaxies from the Catalogue of Isolated Galaxies (CIG) which offers a vertical view of their disc structure. We aim to study extraplanar Diffuse Ionized Gas (eDIG) by comparing the old and young disc components traced by near-infrared (NIR) and Ultraviolet (UV) imaging with the H$\alpha$ emission structure. We obtained H$\alpha$ monochromatic maps from the Fabry-Perot (FP) interferometry, while the old and young discs structures are obtained from the photometric analysis of the 2MASS K$_{s}$-band, and GALEX NUV and FUV images, thereby identifying the stellar disc and whether the eDIG is present. The H$\alpha$ morphology is peculiar in CIG 71, CIG 183, CIG 593 showing clear asymmetries. In general, geometric parameters (isophotal position angle, peak light distribution, inclination) measured from H$\alpha$, UV and NIR show minimal differences (e.g. $\Delta i\leq\pm$10$^{\circ}$), suggesting that interaction does not play a significant role in shaping the morphology, as expected in isolated galaxies. From H$\alpha$ maps, the eDIG was detected vertically in 11 out of 14 galaxies. Although the fraction of eDIG is high, the comparison between our sample and a generic sample of inclined spirals suggests that the phenomenon is uncorrelated to the galaxy environment. As suggested by the extraplanar UV emission found in 13 out of 14 galaxies the star formation extends well beyond the disc defined by the H$\alpha$ map.

The isolated globule B335 contains a single, low luminosity Class 0 protostar associated with a bipolar nebula and outflow system seen nearly perpendicular to its axis. We observed the innermost regions of this outflow as part of JWST/NIRCam GTO program 1187, primarily intended for wide-field slitless spectroscopy of background stars behind the globule. We find a system of expanding shock fronts with kinematic ages of only a few decades emerging symmetrically from the position of the embedded protostar, which is not directly detected at NIRCam wavelengths. The innermost and youngest of the shock fronts studied here shows strong emission from CO. The next older shock front shows less CO and the third shock front shows only H_2 emission in our data. This third and most distant of these inner shock fronts shows substantial evolution of its shape since it was last observed with high spatial resolution in 1996 with Keck/NIRC. This may be evidence of a faster internal shock catching up with a slower one and of the two shocks merging.

We report dynamical mass measurements of the individual stars in the most luminous and massive stellar member of the nearby Ophiuchus star-forming region, the young tight binary system S1. We combine 28 archival datasets with seven recent, proprietary VLBA observations obtained as part of the \textit{Dynamical Masses of Young Stellar Multiple Systems with the VLBA} project (DYNAMO--VLBA), to constrain the astrometric and orbital parameters of the system, and recover high accuracy dynamical masses. The primary component, S1A, is found to have a mass of 4.11$\pm$0.10~M$_\odot$, significantly less than the typical value, $\sim$~6~M$_\odot$ previously reported in the literature. We show that the spectral energy distribution of S1A can be reproduced by a reddened blackbody with a temperature between roughly 14,000~K and 17,000~K. According to evolutionary models, this temperature range corresponds to stellar masses between 4~M$_\odot$ and 6~M$_\odot$ so the SED is not a priori inconsistent with the dynamical mass of S1A. The luminosity of S1 derived from SED-fitting, however, is only consistent with models for stellar masses above 5~M$_\odot$. Thus, we cannot reconcile the evolutionary models with the dynamical mass measurement of S1A: the models consistent with the location of S1A in the HR diagram correspond to masses at least 25\% higher than the dynamical mass. For the secondary component, S1B, a mass of 0.831~$\pm$~0.014~M$_\odot $ is determined, consistent with a low-mass young star. While the radio flux of S1A remains roughly constant throughout the orbit, the flux of S1B is found to be higher near the apastron.

The HERMES (High Energy Rapid Modular Ensemble of Satellites) Pathfinder mission aims to develop a constellation of nanosatellites to study astronomical transient sources, such as gamma-ray bursts, in the X and soft $\gamma$ energy range, exploiting a novel inorganic scintillator. This study presents the results obtained describing, with an empirical model, the unusually intense and long-lasting residual emission of the GAGG:Ce scintillating crystal after irradiating it with high energy protons (70 MeV) and ultraviolet light ($\sim$ 300 nm). From the model so derived, the consequences of this residual luminescence for the detector performance in operational conditions has been analyzed. It was demonstrated that the current generated by the residual emission peaks at 1-2 pA, thus ascertaining the complete compatibility of this detector with the HERMES Pathfinder nanosatellites.

Aims: Ab initio global particle-in-cell (PIC) simulations of compact neutron star magnetospheres in the align rotator configuration to investigate the role of GR and plasma supply on the polar cap particle acceleration efficiency - precursor of coherent radio emission. Methods: A new module for the PIC code OSIRIS to model plasma dynamics around compact objects with full self-consistent GR effects is presented. A detailed description of the extensions and implementation methods are provided for the main sub-algorithms of the PIC loop, including the field solver, particle pusher and charge-conserving current-deposit scheme. Results: Leptons are efficiently accelerated in the polar caps of neutron stars with force-free magnetospheres. This solution supports strong poloidal currents, which are easily turned to spacelike at any stellar compactness due to the GR frame-dragging effect. Charge-separated magnetospheric solutions, in opposition, depend on the plasma supply and compactness to activate the polar cap. The furthest away the solution is from force-free, the highest the compactness is required for the polar cap activation. Conclusions: GR effects are crucial to explain the pulsar mechanism for low obliquity rotators. Focusing on the aligned rotator, we show that GR relaxes the minimum required poloidal magnetospheric current for the transition of the polar cap to the accelerator regime, thus justifying the observation of weak pulsars beyond the expected death line. Also, we demonstrate how the interplay between the polar cap and outer gaps might explain the intermittent behaviour of the measured spin-down luminosity and the existence of radio sub-pulse nullings for older pulsars.

We report on a detailed spatial and spectral analysis of the large-scale X-ray emission from the merging cluster Cygnus A. We use 2.2 Ms Chandra and 40 ks XMM-Newton archival datasets to determine the thermodynamic properties of the intracluster gas in the merger region between the two sub-clusters in the system. These profiles exhibit temperature enhancements that imply significant heating along the merger axis. Possible sources for this heating include the shock from the ongoing merger, past activity of the powerful AGN in the core, or a combination of both. To distinguish between these scenarios, we compare the observed X-ray properties of Cygnus A with simple, spherical cluster models. These models are constructed using azimuthally averaged density and temperature profiles determined from the undisturbed regions of the cluster and folded through MARX to produce simulated Chandra observations. The thermodynamic properties in the merger region from these simulated X-ray observations were used as a baseline for comparison with the actual observations. We identify two distinct components in the temperature structure along the merger axis, a smooth, large-scale temperature excess we attribute to the ongoing merger, and a series of peaks where the temperatures are enhanced by 0.5-2.5 keV. If these peaks are attributable to the central AGN, the location and strength of these features imply that Cygnus A has been active for the past 300 Myr injecting a total of $\sim$10$^{62}$ erg into the merger region. This corresponds to $\sim$10% of the energy deposited by the merger shock.

Flux excesses in the early time light curves of Type Ia supernovae (SNe\,Ia) are predicted by multiple theoretical models and have been observed in a number of nearby SNe\,Ia over the last decade. However, the astrophysical processes that cause these excesses may affect their use as standardizable candles for cosmological parameter measurements. In this paper, we perform a systematic search for early-time excesses in SNe\,Ia observed by the Zwicky Transient Facility (ZTF) to study whether SNe\,Ia with these excesses yield systematically different Hubble residuals. We analyze two compilations of ZTF SN\,Ia light curves from its first year of operations: 127 high-cadence light curves from \citet{Yao19} and 305 light curves from the ZTF cosmology data release of \citet{Dhawan22}. We detect significant early-time excesses for 17 SNe\,Ia in these samples and find that the excesses have an average $g-r$ color of $0.06\pm0.09$~mag; we do not find a clear preference for blue excesses as predicted by several models. Using the SALT3 model, we measure Hubble residuals for these two samples and find that excess-having SNe\,Ia may have lower Hubble residuals (HR) after correcting for shape, color, and host-galaxy mass, at $\sim$2-3$\sigma$ significance; our baseline result is $\Delta HR = -0.056 \pm 0.026$~mag ($2.2 \sigma$). We compare the host-galaxy masses of excess-having and no-excess SNe\,Ia and find they are consistent, though at marginal significance excess-having SNe\,Ia may prefer lower-mass hosts. Additional discoveries of early excess SNe\,Ia will be a powerful way to understand potential biases in SN\,Ia cosmology and probe the physics of SN\,Ia progenitors.

We have carried out a SEM-EPMA-TEM study to determine the textures and compositions of relict primary iron sulfides and their alteration products in a suite of moderately to heavily-altered CM1 carbonaceous chondrites. We observed four textural groups of altered primary iron sulfides: 1) pentlandite+phyllosilicate (2P) grains, characterized by pentlandite with submicron lenses of phyllosilicates, 2) pyrrhotite+pentlandite+magnetite (PPM) grains, characterized by pyrrhotite-pentlandite exsolution textures with magnetite veining and secondary pentlandite, 3) pentlandite+serpentine (PS) grains, characterized by relict pentlandite exsolution, serpentine, and secondary pentlandite, and 4) pyrrhotite+pentlandite+magnetite+serpentine (PPMS) grains, characterized by features of both the PPM and PS grains. We have determined that all four groups were initially primary iron sulfides, which formed from crystallization of immiscible sulfide melts within silicate chondrules in the solar nebula. The fact that such different alteration products could result from the same precursor sulfides within even the same meteorite sample further underscores the complexity of the aqueous alteration environment for the CM chondrites. The different alteration reactions for each textural group place constraints on the mechanisms and conditions of alteration with evidence for acidic environments, oxidizing environments, and changing fluid compositions (Ni-bearing and Si-Mg-bearing).

Constraining the physical conditions of the ionized media in the vicinity of an active supermassive black hole (SMBH) is crucial to understanding how these complex systems operate. Metal emission lines such as iron (Fe) are useful probes to trace the gaseous media's abundance, activity, and evolution in these accreting systems. Among these, the FeII emission has been the focus of many prior studies to investigate the energetics, kinematics, and composition of the broad-emission line region (BELR) from where these emission lines are produced. In this work, we present the first simultaneous FeII modeling in the optical and near-infrared (NIR) regions. We use CLOUDY photoionization code to simulate both spectral regions in the wavelength interval 4000-12000 Angstroms. We compare our model predictions with the observed line flux ratios for IZw1 - a prototypical strong FeII-emitting active galactic nuclei (AGN). This allows putting constraints on the BLR cloud density and metal content that is optimal for the production of the FeII emission, which can be extended to IZw1-like sources, by examining a broad parameter space. We demonstrate the salient and distinct features of the FeII pseudo-continuum in the optical and NIR, giving special attention to the effect of micro-turbulence on the intensity of the FeII emission.

We present the full Hubble diagram of photometrically-classified Type Ia supernovae (SNe Ia) from the Dark Energy Survey supernova program (DES-SN). DES-SN discovered more than 20,000 SN candidates and obtained spectroscopic redshifts of 7,000 host galaxies. Based on the light-curve quality, we select 1635 photometrically-identified SNe Ia with spectroscopic redshift 0.10$< z <$1.13, which is the largest sample of supernovae from any single survey and increases the number of known $z>0.5$ supernovae by a factor of five. In a companion paper, we present cosmological results of the DES-SN sample combined with 194 spectroscopically-classified SNe Ia at low redshift as an anchor for cosmological fits. Here we present extensive modeling of this combined sample and validate the entire analysis pipeline used to derive distances. We show that the statistical and systematic uncertainties on cosmological parameters are $\sigma_{\Omega_M,{\rm stat+sys}}^{\Lambda{\rm CDM}}=$0.017 in a flat $\Lambda$CDM model, and $(\sigma_{\Omega_M},\sigma_w)_{\rm stat+sys}^{w{\rm CDM}}=$(0.082, 0.152) in a flat $w$CDM model. Combining the DES SN data with the highly complementary CMB measurements by Planck Collaboration (2020) reduces uncertainties on cosmological parameters by a factor of 4. In all cases, statistical uncertainties dominate over systematics. We show that uncertainties due to photometric classification make up less than 10% of the total systematic uncertainty budget. This result sets the stage for the next generation of SN cosmology surveys such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time.

The paradigm-changing possibility of collective neutrino-antineutrino oscillations was recently advanced in analogy to collective flavor oscillations. However, the amplitude for the backward scattering process $\nu_{\mathbf{p}_1}\overline\nu_{\mathbf{p}_2}\to\nu_{\mathbf{p}_2}\overline\nu_{\mathbf{p}_1}$ is helicity-suppressed and vanishes for massless neutrinos, implying that there is no off-diagonal refractive index between $\nu$ and $\overline\nu$ of a single flavor of massless neutrinos. For a nonvanishing mass, collective helicity oscillations are possible, representing de-facto $\nu$--$\overline\nu$ oscillations in the Majorana case. However, such phenomena are suppressed by the smallness of neutrino masses as discussed in the previous literature.

Polynomial inflation is a very simple and well motivated scenario. A potential with a concave ``almost'' saddle point at field value $\phi = \phi_0$ fits well the cosmic microwave background (CMB) data and makes testable predictions for the running of the spectral index and the tensor to scalar ratio. In this work we analyze leptogenesis in the polynomial inflation framework. We delineate the allowed parameter space giving rise to the correct baryon asymmetry as well as being consistent with data on neutrino oscillations. To that end we consider two different reheating scenarios. $(i)$ If the inflaton decays into two bosons, the reheating temperature can be as high as $T_\text{rh} \sim 10^{14}$ GeV without spoiling the flatness of the potential, allowing vanilla $N_1$ thermal leptogenesis to work if $T_\text{rh}> M_1$ where $N_1$ is the lightest right--handed neutrino and $M_1$ its mass. Moreover, if the dominant decay of the inflaton is into Higgs bosons of the Standard Model, we find that rare three--body inflaton decays into a Higgs boson plus one light and one heavy neutrino allow leptogenesis even for $T_\text{rh} < M_1$ if the inflaton mass is of order $10^{12}$ GeV or higher; in the polynomial inflation scenario this requires $\phi_0 \gtrsim 2.5~M_P$. This novel mechanism of non--thermal leptogenesis is quite generic, since the coupling leading to the three--body final state is required in the type I see--saw mechanism. $(ii)$ If the inflaton decays into two fermions, the flatness of the potential implies a lower reheating temperature. In this case inflaton decay to two $N_1$ still allows successful non--thermal leptogenesis if $\phi_0 \gtrsim 0.1~M_P$ and $T_\text{rh} \gtrsim 10^{6}$ GeV.

We explore the prospects for identifying differences in simulated gravitational-wave signals of binary neutron star (BNS) mergers associated with the way thermal effects are incorporated in the numerical-relativity modelling. We consider a hybrid approach in which the equation of state (EoS) comprises a cold, zero temperature, piecewise-polytropic part and a thermal part described by an ideal gas, and a tabulated approach based on self-consistent, microphysical, finite-temperature EoS. We use time-domain waveforms corresponding to BNS merger simulations with four different EoS. Those are injected into Gaussian noise given by the sensitivity of the third-generation detector Einstein Telescope and reconstructed using BayesWave, a Bayesian data-analysis algorithm that recovers the signals through a model-agnostic approach. The two representations of thermal effects result in frequency shifts of the dominant peaks in the spectra of the post-merger signals, for both the quadrupole fundamental mode and the late-time inertial modes. For some of the EoS investigated those differences are large enough to be told apart, especially in the early post-merger phase when the signal amplitude is the loudest. These frequency shifts may result in differences in the inferred tidal deformability, which might be resolved by third-generation detectors up to distances of about tens of Mpc at most.

The Milky Way galaxy is estimated to be home to ten million to a billion stellar-mass black holes (BHs). Accurately determining this number and distribution of BH masses can provide crucial information about the processes involved in BH formation, the possibility of the existence of primordial BHs, and interpreting gravitational wave (GW) signals detected in LIGO-VIRGO-KAGRA. Sahu et al. recently confirmed one isolated stellar-mass BH in our galaxy using astrometric microlensing. This work proposes a novel method to identify such BHs using the gravitational analog of the Gertsenshtein-Zel'dovich (GZ) effect. We explicitly demonstrate the generation of GWs when a kilohertz(kHz) electromagnetic (EM) pulse from a pulsar is intervened by a spherically symmetric compact object situated between the pulsar and Earth. Specifically, we show that the curvature of spacetime acts as the catalyst, akin to the magnetic field in the GZ effect. Using the covariant semi-tetrad formalism, we quantify the GW generated from the EM pulse through the Regge-Wheeler tensor and express the amplitude of the generated GW in terms of the EM energy and flux. We demonstrate how GW detectors can detect stellar-mass BHs by considering known pulsars within our galaxy. This approach has a distinct advantage in detecting stellar mass BHs at larger distances since the GW amplitude falls as $1/r$.

We investigate quantum decoherence in a class of models which interpolates between expanding (inflation) and contracting (ekpyrosis) scenarios. For the cases which result in a scale-invariant power spectrum, we find that ekpyrotic universes lead to complete decoherence of the curvature perturbation before the bounce. This is in stark contrast to the inflationary case, where recoherence has been previously observed in some situations. Although the purity can be computed for couplings of all sizes, we also study the purity perturbatively and observe that late-time (secular growth) breakdown of perturbation theory often occurs in these cases. Instead, we establish a simple yet powerful late-time purity resummation which captures the exact evolution to a remarkable level, while maintaining analytical control. We conclude that the cosmological background plays a crucial role in the decoupling of the heavy fields during inflation and alternatives.

In this work we present a novel approach to the study of cosmological particle production in asymptotically Minkowski spacetimes. We emphasize that it is possible to determine the amount of particle production by focusing on the mathematical properties of the mode function equations,i.e. their singularities and monodromies, sidestepping the need to solve those equations. We consider in detail creation of scalar and spin 1/2 particles in four dimensional asymptotically Minkowski flat FLRW spacetimes. We explain that when the mode function equation for scalar fields has only regular singular points, the corresponding scale factors are asymptotically Minkowski. For Dirac spin 1/2 fields, the requirement of mode function equations with only regular points is more restrictive, and picks up a subset of the aforementioned scale factors. For the scalar case, we argue that there are two different regimes of particle production; while most of the literature has focused on only one of these regimes, the other regime presents enhanced particle production. On the other hand, for Dirac fermions we find a single regime of particle production. Finally, we very briefly comment on the possibility of studying particle production in spacetimes that don't asymptote to Minkowski, by considering mode function equations with irregular singular points.

The quanta of phantom dark energy models are negative energy particles, whose maximum magnitude of energy $|E|$ must be less than a cutoff $\Lambda\lesssim 20\,$MeV. They are produced by spontaneous decay of the vacuum into phantoms plus normal particles. I review general cosmological constraints that have been derived from the effects of such phantom fluid production, and a possible application: the generation of boosted dark matter or radiation that could be directly detected. Recent excess events from the DAMIC experiment can be well-fit by such processes.