Papers for Friday, Feb 24 2023

Classification is a popular task in the field of Machine Learning (ML) and Artificial Intelligence (AI), and it happens when outputs are categorical variables. There are a wide variety of models that attempts to draw some conclusions from observed values, so classification algorithms predict categorical class labels and use them in classifying new data. Popular classification models including logistic regression, decision tree, random forest, Support Vector Machine (SVM), multilayer perceptron, Naive Bayes, and neural networks have proven to be efficient and accurate applied to many industrial and scientific problems. Particularly, the application of ML to astronomy has shown to be very useful for classification, clustering, and data cleaning. It is because after learning computers, these tasks can be done automatically by them in a more precise and more rapid way than human operators. In view of this, in this paper, we will review some of these popular classification algorithms, and then we apply some of them to the observational data of nonvariable and the RR Lyrae variable stars that come from the SDSS survey. For the sake of comparison, we calculate the accuracy and F1-score of the applied models.

The Hot Big Bang is often considered as the origin of all matter and radiation in the Universe. Primordial nucleosynthesis (BBN) provides strong evidence that the early Universe contained a hot plasma of photons and baryons with a temperature $T>\text{MeV}$. However, the earliest probes of dark matter originate from much later times around the epoch of structure formation. In this work we describe a scenario in which dark matter (and possibly dark radiation) can be formed around or even after BBN in a second Big Bang which we dub the ``Dark Big Bang''. The latter occurs through a phase transition in the dark sector which transforms dark vacuum energy into a hot dark plasma of particles; in this paper we focus on a first-order phase transition for the Dark Big Bang. The correct dark matter abundance can be set by dark matter cannibalism or by pair-annihilation within the dark sector followed by a thermal freeze-out. Alternatively ultra-heavy ``dark-zilla'' dark matter can originate directly from bubble collisions during the Dark Big Bang. We will show that the Dark Big Bang is consistent with constraints from structure formation and the Cosmic Microwave Background (CMB) if it occurred when the Universe was less than one month old, corresponding to a temperature in the visible sector above $\mathcal{O}$(keV). While the dark matter evades direct and indirect detection, the Dark Big Bang gives rise to striking gravity wave signatures to be tested at pulsar timing array experiments. Furthermore, the Dark Big Bang allows for realizations of self-interacting and/or warm dark matter which suggest exciting discovery potential in future small-scale structure observations.

We present a seven band (g, r, i, z, y, NB816, NB921) catalog derived from a Subaru Hyper Suprime-Cam (HSC) imaging survey of the North Ecliptic Pole (NEP). The survey, known as HEROES, consists of 44 sq. deg of contiguous imaging reaching median 5-sigma depths of g: 26.5, r: 26.2, i: 25.7, z: 25.1, y: 23.9, NB816: 24.4, NB921: 24.4 mag. We reduced these data with the HSC pipeline software hscPipe, and produced a resulting multiband catalog containing over 25 million objects. We provide the catalog in three formats: (1) a collection of hscPipe format forced photometry catalogs, (2) a single combined catalog containing every object in that dataset with selected useful columns, and (3) a smaller variation of the combined catalog with only essential columns for basic analysis or low memory machines. The catalog uses all the available HSC data on the NEP and may serve as the primary optical catalog for current and future NEP deep fields from instruments and observatories such as SCUBA-2, eROSITA, Spitzer, Euclid, and JWST.

We present reverberation mapping measurements for the prominent ultraviolet broad emission lines of the active galactic nucleus Mrk817 using 165 spectra obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Our ultraviolet observations are accompanied by X-ray, optical, and near-infrared observations as part of the AGN Space Telescope and Optical Reverberation Mapping Program 2 (AGN STORM 2). Using the cross-correlation lag analysis method, we find significant correlated variations in the continuum and emission-line light curves. We measure rest-frame delayed responses between the far-ultraviolet continuum at 1180 A and Ly$\alpha$ $\lambda1215$ A ($10.4_{-1.4}^{+1.6}$ days), N V $\lambda1240$ A ($15.5_{-4.8}^{+1.0}$days), SiIV + OIV] $\lambda1397$ A ($8.2_{-1.4}^{+1.4}$ days), CIV $\lambda1549$ A ($11.8_{-2.8}^{+3.0}$ days), and HeII $\lambda1640$ A ($9.0_{-1.9}^{+4.5}$ days) using segments of the emission-line profile that are unaffected by absorption and blending, which results in sampling different velocity ranges for each line. However, we find that the emission-line responses to continuum variations are more complex than a simple smoothed, shifted, and scaled version of the continuum light curve. We also measure velocity-resolved lags for the Ly$\alpha$, and CIV emission lines. The lag profile in the blue wing of Ly$\alpha$ is consistent with virial motion, with longer lags dominating at lower velocities, and shorter lags at higher velocities. The CIV lag profile shows the signature of a thick rotating disk, with the shortest lags in the wings, local peaks at $\pm$ 1500 $\rm km\,s^{-1}$, and a local minimum at line center. The other emission lines are dominated by broad absorption lines and blending with adjacent emission lines. These require detailed models, and will be presented in future work.

Solid-phase photo-processes involving icy dust grains greatly affect the chemical evolution of the interstellar medium by leading to the formation of complex organic molecules and by inducing photodesorption. So far, the focus of laboratory studies has been mainly on the impact of energetic ultraviolet (UV) photons on ices, but direct vibrational excitation by infrared (IR) photons is expected to influence the morphology and content of interstellar ices as well. However, little is still known about the mechanisms through which this excess vibrational energy is dissipated, and its implications on the structure and ice photochemistry. In this work, we present a systematic investigation of the behavior of interstellar relevant CO and CH3OH ice analogues upon resonant excitation of vibrational modes using tunable infrared radiation, leading to both the quantification of infrared-induced photodesorption and insights in the impact of vibrational energy dissipation on ice morphology. We utilize an ultrahigh vacuum setup at cryogenic temperatures to grow pure CO and CH3OH ices, as well as mixtures of the two. We expose the ices to intense, near-monochromatic mid-infrared free-electron-laser radiation to selectively excite the species. The dissipation of vibrational energy is observed to be highly dependent on the excited mode and the chemical environment of the ice. All amorphous ices undergo some degree of restructuring towards a more organized configuration upon on-resonance irradiation. Moreover, IR-induced photodesorption is observed to occur for both pure CO and CH3OH ices, with interstellar photodesorption efficiencies of the order of 10 molecules cm-2 s-1 (i.e., comparable to or higher than UV-induced counterparts). Indirect photodesorption of CO upon vibrational excitation of CH3OH in ice mixtures is also observed to occur, particularly in environments rich in methanol.

We present ALMA multiwavelength observations of the protoplanetary disk around the nearby (d$\sim$100 pc) young solar analog MP Mus (PDS 66). These observations at 0.89 mm, 1.3 mm, and 2.2 mm have angular resolutions of $\sim$ 1", 0.05", and 0.25", respectively, and probe the dust and gas in the system with unprecedented detail and sensitivity. The disk appears smooth down to the 4 au resolution of the 1.3 mm observations, in contrast with most disks observed at comparable spatial scales. The dust disk has a radius of 60$\pm$5 au, a dust mass of $0.14_{-0.06}^{+0.11} M_{\rm Jup}$, and a mm spectral index $<2$ in the inner 30 au, suggesting optically thick emission from grains with high albedo in this region. Several molecular gas lines are also detected extending up to 130$\pm$15 au, similar to small grains traced by scattered light observations. Comparing the fluxes of different CO isotopologues with previous models yields a gas mass of $0.1-1 M_{\rm Jup}$, implying a gas to dust ratio of 1-10. We also measure a dynamical stellar mass of $M_{\rm dyn}$=1.30$\pm$0.08 $M_\odot$ and derive an age of 7-10 Myr for the system. The survival of large grains in an evolved disk without gaps/rings is surprising, and it is possible that existing substructures remain undetected due to optically thick emission at 1.3 mm. Alternatively, small structures may still remain unresolved with the current observations. Based on simple scaling relations for gap-opening planets and gap widths, this lack of substructures places upper limits to the masses of planets in the disk as low as 2 $M_\oplus$-0.06 $M_{\rm Jup}$ at $r > 40$ au. The lack of mm emission at radii $r > 60$ au also suggests that the gap in scattered light between 30-80 au is likely not a gap in the disk density, but a shadow cast by a puffed-up inner disk.

We examine to what extent disk and outflow models can reproduce observations of H I gas within a few virial radii of galaxies in pairs and groups. Using highly-sensitive HST/COS and FOS spectra of the Q0107 quasar triplet covering Ly$\alpha$ for z$\lesssim$1, as well as a deep galaxy redshift survey including VIMOS, DEIMOS, GMOS and MUSE data, we test simple disk and outflow models against the H I absorption along three lines-of-sight (separated by 200-500 kpc) through nine galaxy groups in this field. These can be compared with our previous results in which these models can often be fit to the absorption around isolated galaxies. Our models can reproduce $\approx$ 75$\%$ of the 28 identified absorption components within 500 km/s of a group galaxy, so most of the H I around groups is consistent with a superposition of the CGM of the individual galaxies. Gas stripped in interactions between galaxies may be a plausible explanation for some of the remaining absorption, but neither the galaxy images nor the galaxy and absorber kinematics provide clear evidence of such stripped material, and these unexplained absorbers do not preferentially occur around close pairs of galaxies. We find H I column densities typically higher than at similar impact parameters around isolated galaxies ($\approx$ 2.5$\sigma$), as well as more frequent detections of O VI than around isolated galaxies (30$\%$ of sightlines to 7$\%$).

Fingering convection (also known as thermohaline convection) is a process that drives the radial transport of chemical elements in regions of stellar radiative zones where the mean molecular weight increases with radius. It is often accounted for in one-dimensional stellar evolution calculations by using simplified mixing prescriptions. Recently, Harrington & Garaud (2019) used three-dimensional direct numerical simulations to show that a vertical magnetic field can dramatically enhance the rate of chemical mixing by fingering convection. Furthermore, they proposed a so-called "parasitic saturation" theory to model this process. Here, we test their model over a broad range of parameter space, using a suite of direct numerical simulations of magnetized fingering convection varying the magnetic Prandtl number, magnetic field strength, and composition gradient. We find that the rate of chemical mixing is not always predicted accurately by this existing model, in particular when the magnetic diffusivity is large. We present a first attempt at extending the Harrington & Garaud (2019) model to solve this issue, but are forced to conclude that some of core assumptions of the parasitic saturation paradigm must be revisited in the magnetized case.

With about a hundred mergers of binary black holes (BBHs) detected via gravitational waves by the LIGO-Virgo-KAGRA (LVK) Collaboration, our understanding of the darkest objects in the Universe has taken unparalleled steps forward. While most of the events are expected to consist of BHs directly formed from the collapse of massive stars, some may contain the remnants of previous BBH mergers. In the most massive globular clusters and in nuclear star clusters, successive mergers can produce second- (2G) or higher-generation BHs, and even form intermediate-mass BHs. Overall, we predict that up to $\sim 10\%$, $\sim 1\%$ or $\sim 0.1\%$ of the BBH mergers have one component being a 2G, 3G, or 4G BH, respectively. Assuming that $\sim 500$ BBH mergers will be detected in O4 by LVK, this means that $\sim 50$, $\sim 5$, or $\sim 0.5$ events, respectively, will involve a 2G, 3G, or 4G BH, if most sources are produced dynamically in dense star clusters. With their distinctive signatures of higher masses and spins, such hierarchical mergers offer an unprecedented opportunity to learn about the BH populations in the densest stellar systems and to shed light on the elusive intermediate-mass BHs that may form therein.

Constraints on the linear growth rate, $f\sigma_8$, using small scale redshift space distortion measurements have a significant statistical advantage over those made on large scales. However, these measurements need to carefully disentangle the linear and non-linear information when interpreting redshift space distortions in terms of $f\sigma_8$. It is particularly important to do this given that some previous measurements found a significant deviation from the expectation based on the $\Lambda$CDM model constrained by Planck CMB data. We construct a new emulator-based model for small scale galaxy clustering with scaling parameters for both the linear and non-linear velocities of galaxies, allowing us to isolate the linear growth rate. We train the emulator using simulations from the AbacusCosmos suite, and apply it to data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) luminous red galaxy sample. We obtain a value of $f\sigma_8(z=0.737)=0.368\pm0.041$, in 2.3-$\sigma$ tension with the Planck 2018 $\Lambda$CDM expectation, and find less dependence on the minimum measurement scale than previous analyses.

The project MOMO (Multiwavelength Observations and Modelling of OJ 287) was set up to test predictions of binary supermassive black hole (SMBH) scenarios and to understand disk-jet physics of the blazar OJ 287. After a correction, the precessing binary (PB) SMBH model predicted the next main outburst of OJ 287 in 2022 October, no longer in July, making the outburst well observable and the model testable. We have densely covered this period in our ongoing multi-frequency radio, optical, UV, and X-ray monitoring. The predicted outburst was not detected. Instead, OJ 287 was at low optical-UV emission levels, declining further into November. The predicted thermal bremsstrahlung spectrum was not observed either, at any epoch. Further, applying scaling relations, we estimate a SMBH mass of OJ 287 of 10^8 M_sun. The latest in a sequence of deep low-states that recur every 1-2 yrs is used to determine an upper limit on the Eddington ratio and on the accretion-disk luminosity. This limit is at least a factor of 10 lower than required by the PB model with its massive primary SMBH of >10^{10} M_sun. All these results favor alternative binary SMBH models of OJ 287 that neither require strong orbital precession nor a very large mass of the primary SMBH.

Magnetohydrodynamical (MHD) turbulence is ubiquitous in magnetized astrophysical plasmas, and it radically changes a great variety of astrophysical processes. In this review, we give the concept of MHD turbulence and explain the origin of its scaling. We consider the implications of MHD turbulence to various problems: dynamo in different types of stars, flare activity, solar and stellar wind from different stars, propagation of cosmic rays, and star formation. We also discuss how the properties of MHD turbulence provide a new way of tracing magnetic fields in interstellar and intracluster media.

Recent X-ray polarimetric data on the prototypical black hole X-ray binary Cyg~X-1 from the Imaging X-ray Polarimetry Explorer present tight constraints on accretion geometry in the hard spectral state. Contrary to general expectations of a low, $\lesssim1\%$ polarization degree (PD), the observed average PD was found to be a factor of 4 higher. Aligned with the jet position angle on the sky, the observed polarization favors geometry of the X-ray emission region stretched in the accretion disk plane. The high PD is, however, difficult to reconcile with the low orbital inclination of the binary $i\approx30^\circ$. We suggest that this puzzle can be explained if the emitting plasma is outflowing with a mildly relativistic velocity $\gtrsim0.4\,c$. Our radiative transfer simulations show that Comptonization in the outflowing medium radiates X-rays with the degree and direction of polarization consistent with observations at $i\approx30^\circ$.

In this work, we investigate the influence of the phase transition and a stiffer fluid in neutron stars' cores on the static equilibrium configuration, dynamical stability, and tidal deformability. For this aim, it is taken into account that the fluid in the core and the envelope follow the relativistic polytropic equation of state. We find that the phase transition and a stiffer fluid in the core will reflect in the total mass, radius, speed of sound, core radius, radial stability with a slow and rapid conversion at the interface, and tidal deformability. We also investigate the dimensionless tidal deformability $\Lambda_1$ and $\Lambda_2$ for a binary neutron stars system with chirp mass equal to GW$170817$. Finally, we contrast our results with observational data to show the role that phase transition and a stiffer core fluid could play in the study of neutron stars.

The NASA/Dawn mission has acquired an unprecedented amount of data from the surface of the dwarf planet Ceres, providing a thorough characterization of its surface composition. The current favorite compositional model includes a mixture of ultra-carbonaceous material, phyllosilicates, carbonates, organic, Fe-oxides, and volatiles, as determined by the Dawn/VIR-IR imaging spectrometer, with the constraints on the elemental composition provided by the Dawn/GRAND instrument. The recent calibration of the VIS channel of the VIR spectrometer provided further constraints on the overall compositional model, suggesting another type of silicate not previously considered. This work reviews the VIR-VIS calibration process, and the resulting compositional model of Ceres surface considering both VIS and IR channels. Laboratory experiments are needed for a better understanding of the global (VIS-IR) shape of Ceres average surface, and in particular for the interpretation of the emergent absorption at 1{\mu}m. Current and future experiments are also discussed in this work.

We report observations of nine MBCs or candidate MBCs, most of which were obtained when the targets were apparently inactive. We find effective nucleus radii (assuming albedos of p_V=0.05+/-0.02) of r_n=(0.24+/-0.05) km for 238P/Read, r_n=(0.9+/-0.2) km for 313P/Gibbs, r_n=(0.6+/-0.1) km for 324P/La Sagra, r_n=(1.0+/-0.2) km for 426P/PANSTARRS, r_n=(0.5+/-0.1) km for 427P/ATLAS, r_n<(0.3+/-0.1) km for P/2016 J1-A (PANSTARRS), r_n<(0.17+/-0.04) km for P/2016 J1-B (PANSTARRS), r_n<(0.5+/-0.2) km for P/2017 S9 (PANSTARRS), and r_n=(0.4+/-0.1) km for P/2019 A3 (PANSTARRS). We identify evidence of activity in observations of 238P in 2021, and find similar inferred activity onset times and net initial mass loss rates for 238P during perihelion approaches in 2010, 2016, and 2021. P/2016 J1-A and P/2016 J1-B are also found to be active in 2021 and 2022, making them collectively the tenth MBC confirmed to be recurrently active near perihelion and therefore likely to be exhibiting sublimation-driven activity. The nucleus of 313P is found to have colors of g'-r'=0.52+/-0.05 and r'-i'=0.22+/-0.07, consistent with 313P being a Lixiaohua family member. We also report non-detections of P/2015 X6 (PANSTARRS), where we conclude that its current nucleus size is likely below our detection limits (r_n<0.3 km). Lastly, we find that of 17 MBCs or candidate MBCs for which nucleus sizes (or inferred parent body sizes) have been estimated, >80% have r_n<1.0 km, pointing to an apparent physical preference toward small MBCs, where we suggest that YORP spin-up may play a significant role in triggering and/or facilitating MBC activity.

Rapidly spinning magnetic grains can acquire large magnetic dipole moments due to the Barnett effect. Here we study the new effect of Barnett magnetic dipole-dipole interaction on grain-grain collisions and grain growth, assuming grains spun up by radiative torques. We find that the collision rate between grains having embedded iron inclusions can be significantly enhanced due to Barnett magnetic dipole-dipole interaction when grains rotate suprathermally by radiative torques. We discuss the implications of enhanced collision rate for grain growth and destruction in the circumstellar envelope of evolved stars, photodissociation regions, and protostellar environments. Our results first reveal the importance of the dust magnetic properties and the local radiation field on grain growth and destruction.

While solar-like oscillations in red giants have been observed at massive scale by the Kepler mission, few features of these oscillation mode frequencies, other than their global properties, have been exploited for stellar characterization. The signatures of acoustic glitches in mode frequencies have been used for studying main-sequence stars, but the validity of applying such techniques to evolved red giants, particularly pertaining to the inclusion of nonradial modes, has been less well-examined. Making use of new theoretical developments, we characterize glitches using the $\pi$ modes associated with red giant stellar models, and use our procedure to examine for the first time how properties of the He II acoustic glitch -- specifically its amplitude and associated acoustic depth -- vary over the course of evolution up the red giant branch, and with respect to other fundamental stellar properties. We find that the acoustic depths of these glitches, in conjunction with other spectroscopic information, discriminates between red giants in the first-ascent and core-helium-burning phases. We critically reexamine previous attempts to constrain acoustic glitches from nonradial (in particular dipole) modes in red giants. Finally, we apply our fitting procedure to Kepler data, to evaluate its robustness to noise and other observational systematics.

Solar filament eruptions often show complex and dramatic geometric deformation that is highly relevant to the underlying physical mechanism triggering the eruptions. It has been well known that the writhe of filament axes is a key parameter characterizing its global geometric deformation, but a quantitative investigation of the development of writhe during its eruption is still lacking. Here we introduce the Writhe Application Toolkit (WAT) which can be used to characterize accurately the topology of filament axes. This characterization is achieved based on the reconstruction and writhe number computation of three-dimensional paths of the filament axes from dual-perspective observations. We apply this toolkit to four dextral filaments located in the northern hemisphere with a counterclockwise (CCW) rotation during their eruptions. Initially, all these filaments possess a small writhe number (=<0.20) indicating a weak helical deformation of the axes. As the CCW rotation kicks in, their writhe numbers begin to decrease and reach large negative values. Combined with the extended C\u{a}lug\u{a}reanu theorem, the absolute value of twist is deduced to decrease during the rotation. Such a quantitative analysis strongly indicates a consequence of the conversion of twist into writhe for the studied events.

Searching for X-ray and gamma-ray bursts, including Gamma-ray bursts (GRBs), Soft Gamma-ray Repeaters (SGRs) and high energy transients associated with Gravitational wave (GW) events or Fast radio bursts (FRBs), is of great importance in the multi-messenger and multi-wavelength era. Although a coherent search based on the likelihood ratio and Gaussian statistics has been established and utilized in many studies, this Gaussian-based method could be problematic for weak and short bursts which usually have very few counts. To deal with all bursts including weak ones, here we propose the coherent search in Poisson statistics. We studied the difference between Poisson-based and Gaussian-based search methods by Monte Carlo simulations, and find that the Poisson-based search method has advantages compared to the Gaussian case especially for weak bursts. Our results show that, for very weak bursts with very low number of counts, the Poisson-based search can provide higher significance than the Gaussian-based search, and its likelihood ratio (for background fluctuation) still generally follows the chi2 distribution, making the significance estimation of searched bursts very convenient. Thus, we suggest that the coherent search based on Poisson likelihood ratio is more appropriate in the search for generic transients including very weak ones.

We report results from the TESS photometric data and new high-resolution spectra of the Algol system X Tri showing short-period pulsations. From the echelle spectra, the radial velocities of the eclipsing pair were measured, and the rotational rate and effective temperature of the primary star were obtained to be $v_1$$\sin$$i=84\pm6$ km s$^{-1}$ and $T_{\rm eff,1}=7900 \pm 110$ K, respectively. The synthetic modeling of these observations implies that X Tri is in synchronous rotation and is physically linked to a visual companion TIC 28391715 at a separation of about 6.5 arcsec. The absolute parameters of our target star were accurately and directly determined to be $M_1 = 2.137 \pm 0.018$ $M_\odot$, $M_2 = 1.101 \pm 0.010$ $M_\odot$, $R_1 = 1.664 \pm 0.010$ $R_\odot$, $R_2 = 1.972 \pm 0.010$ $R_\odot$, $L_1 = 9.67 \pm 0.55$ $L_\odot$, and $L_2 = 2.16 \pm 0.09$ $L_\odot$. The phase-binned mean light curve was used to remove the binary effect from the observed TESS data. Multifrequency analysis of the residuals revealed 16 significant frequencies, of which the high-frequency signals between 37 day$^{-1}$ and 48 day$^{-1}$ can be considered probable pulsation modes. Their oscillation periods of 0.021$-$0.027 days and pulsation constants of 0.014$-$0.018 days are typical values of $\delta$ Sct variables. The overall results demonstrate that X Tri is an oEA star system, consisting of a $\delta$ Sct primary and its lobe-filling companion in the semi-detached configuration.

(Sub) mm VLBI observations are strongly hindered by limited sensitivity, with the fast tropospheric fluctuations being the dominant culprit. We predict great benefits from applying next-generation frequency phase transfer calibration techniques for the next generation Event Horizon Telescope, using simultaneous multi-frequency observations. We present comparative simulation studies to characterise its performance, the optimum configurations, and highlight the benefits of including observations at 85\,GHz along with the 230 and 340\,GHz bands. The results show a transformational impact on the ngEHT array capabilities, with orders of magnitude improved sensitivity, observations routinely possible over the whole year, and ability to carry out micro-arcsecond astrometry measurements at the highest frequencies, amongst others. This will enable the addressing of a host of innovative open scientific questions in astrophysics. We present a solution for highly scatter-broadened sources such as SgrA*, a prime ngEHT target. We conclude that adding the 85\,GHz band provides a pathway to an optimum and robust performance for ngEHT in sub-millimeter VLBI, and strongly recommmend its inclusion in the simultaneous multi-frequency receiver design.

Visible and near-infrared (V+NIR) emission lines were the first to be discovered in the corona, during total solar eclipses, and they continue to offer unique opportunities to study the physical properties of the corona. The most commonly observed coronal emission lines today are in the extreme ultraviolet, which are dominated by collisionally excited emission. V+NIR lines on the other hand, are radiatively excited out to high helioprojective distances. Indeed, recent eclipse observations have demonstrated the diagnostic potential of V+NIR lines, which are still observable out to at least 3.4 Rs. V+NIR lines can be used to infer key plasma parameters such as: the electron and ion temperatures, magnetic field strength and morphologies, the ionic freeze-in distances, Doppler motions of coronal plasmas, and the dynamics of coronal mass ejections through time variations of these parameters. Current and planned space-based coronagraphs, such as Solar Orbiter and Proba 3, will have some filters for V+NIR lines, but will only have an exceptionally small selection. They will thus be limited in their ability to infer electron or ion temperatures, as well as other crucial physical properties of the corona. The ground-based DKIST and UCoMP will soon offer V+NIR line observations, but they will be limited to a maximum helioprojective distance of about 1.5 Rs. To better explore the middle corona, and to understand the formation of the solar wind and space weather events, it is essential that we deploy additional space-based assets to measure a wide selection of V+NIR emission lines at helioprojective distances beyond 1.5 Rs. Occulting of the solar disk could be achieved by a conventional coronagraph, by novel methods such as an external occulter, by lunar occultations in situ in orbit around the Moon, or by lunar based observations of lunar eclipses (i.e., total solar eclipse on the Moon).

There has been an unfortunate gap in coronal emission line observations from space in the visible and near IR (V+NIR). Their distinct scientific advantage stems from the dominance of radiative excitation in their formation, whereby their emission can be detected out to several solar radii above the limb. V+NIR emission lines can thus yield the only inferences of the physical properties of the coronal plasma, such as species temperatures, densities, elemental abundances, and speeds along and perpendicular to the line of sight in this critical spatial span. These diagnostics have been demonstrated with decades of unsurpassed high-resolution imaging and spectroscopic observations during total solar eclipses. This white paper calls for dedicated funding for a Total Solar Eclipse Earth-Based Mission for ground, airborne and seaborne observations of the corona during totality for the next decade starting in 2024. The proposed Mission capitalizes on the unique diagnostic potential offered by the V+NIR coronal emission lines for the inference of key plasma parameters over a distance range of at least 5 Rs from the solar surface. This critical coronal space is currently missing from existing and to-be launched coronagraphic instrumentation in the proposed time frame. Multi-site observing platforms for each eclipse would further capture the temporal variability of coronal plasmas over a time span of at least 1 hour, with a temporal resolution of a fraction of a minute. Furthermore, this Mission offers unsurpassed opportunities for the exploration of new technologies for future implementation with coronagraphs. This Mission has a unique significant broader impact for outreach opportunities to engage the public and the younger generations in heliospheric science from an awe-inspiring cosmic event.

We present updated measurements of the X-ray properties of the pulsar wind nebula associated with the TeV $\gamma$-ray source HESS J1640-465 derived from Chandra & NuSTAR data. We report a high $N_{H}$ value along line of sight, consistent with previous work, which led us to incorporate effects of dust scattering in our spectral analysis. Due to uncertainties in the dust scattering, we report a range of values for the PWN properties (photon index and un-absorbed flux). In addition, we fit the broadband spectrum of this source and found evidence for spectral softening and decreasing unabsorbed flux as we go to higher photon energies. We then used a one zone time dependent evolutionary model to reproduce the dynamical and multi-wavelength spectral properties of our source. Our model suggests a short spin-down time scale, a relatively higher than average magnetized pulsar wind, a strong PWN magnetic field and maximum electron energy up to PeV, suggesting HESS J1640-465 could be a PeVatron candidate.

A fast $N$-body code has been developed for simulating a stellar disk embedded in a live dark matter halo. In generating its Poisson solver, a self-consistent field (SCF) code which inherently possesses perfect scalability is incorporated into a tree code which is parallelized using a library termed Framework for Developing Particle Simulators (FDPS). Thus, the code developed here is called SCF-FDPS. This code has realized the speedup of a conventional tree code by applying an SCF method not only to the calculation of the self-gravity of the halo but also to that of the gravitational interactions between the disk and halo particles. Consequently, in the SCF-FDPS code, a tree algorithm is applied only to calculate the self-gravity of the disk. On a many-core parallel computer, the SCF-FDPS code has performed at least three times, in some case nearly an order of magnitude, faster than an extremely-tuned tree code on it, if the numbers of disk and halo particles are, respectively, fixed for both codes. In addition, the SCF-FDPS code shows that the cpu cost scales almost linearly with the total number of particles and almost inversely with the number of cores. We find that the time evolution of a disk-halo system simulated with the SCF-FDPS code is, in large measure, similar to that obtained using the tree code. We suggest how the present code can be extended to cope with a wide variety of disk-galaxy simulations.

In systems where the standard $\alpha$ effect is inoperative, one often explains the existence of mean magnetic fields by invoking the `incoherent $\alpha$ effect', which appeals to fluctuations of the mean kinetic helicity. Previous studies, while considering fluctuations in the mean kinetic helicity, treated the mean turbulent kinetic energy as a constant, despite the fact that both these quantities involve second-order velocity correlations. The mean turbulent kinetic energy causes both turbulent diffusion and diamagnetic pumping of the mean magnetic field. In this work, we use a double-averaging procedure to analytically show that fluctuations of the mean turbulent kinetic energy (giving rise to $\eta$-fluctuations, where $\eta$ is the turbulent diffusivity) can lead to the growth of a large-scale magnetic field even when the kinetic helicity is zero pointwise. Constraints on the operation of such a dynamo are expressed in terms of dynamo numbers that depend on the correlation length, correlation time, and strength of these fluctuations. In the white-noise limit, we find that these fluctuations reduce the overall turbulent diffusion, while also contributing a drift term which does not affect the growth of the field. We also study the effects of nonzero correlation time and anisotropy. Diamagnetic pumping, which arises due to inhomogeneities in the turbulent kinetic energy, leads to growing mean field solutions even when the $\eta$-fluctuations are isotropic. Our results suggest that fluctuations of the turbulent kinetic energy may be relevant in astrophysical contexts.

We have studied the temporal and spectral properties of SXP 15.3 observed by \textit{NuSTAR} in hard energy range 3-79 keV during late 2018. The timing analysis of \textit{NuSTAR} observation predicts coherent pulsation at $15.2388\;\pm\;0.0002\;s$. The pulse profiles in different energy bands demonstrates energy dependence. The shape of the pulse profile was generally suggestive of a fan-beamed dominated pattern, which when combined with the measured luminosity predicts that the source may be accreting in the super-critical regime. Non-monotonic increase of Pulse Fraction was observed with energy. The \textit{NuSTAR} observation finds that the pulse period of the source have spun-up at a rate of --0.0176 s $yr^{-1}$ when compared to the previous analysis by the same observatory more than 1 year ago. The source flux in the present \textit{NuSTAR} study in 3-79 keV energy range is $\sim\;1.36\;\times\;10^{-10}\;erg\;cm^{-2}\;s^{-1}$ which corresponds to a luminosity of $\sim\;6\;\times\;10^{37}\;erg\;s^{-1}$. The cyclotron line energy of the source is detected at $\sim$5 keV. Pulse phase resolved spectroscopy shows that the cyclotron line energy varies significantly with pulse phase and the photon index becomes softer with increasing flux. In addition, we have studied the evolution of luminosity with time using 2017 and 2018 \textit{Swift/XRT} observations. The analysis of \textit{Swift/XRT} data reveals that the photon index are positively correlated with the source luminosity which is a characteristic of super-critical accretion phenomena.

The abundances of low First Ionisation Potential (FIP) elements are three to four times higher (FIP bias) in the closed loop active corona than in the photosphere, known as the FIP effect. Observations suggest that the abundances vary in different coronal structures. Here, we use the soft X-ray spectroscopic measurements from the Solar X-ray Monitor (XSM) on board the Chandrayaan-2 orbiter to study the FIP effect in multiple A-class flares observed during the minimum of solar cycle 24. Using time-integrated spectral analysis, we derive the average temperature, emission measure, and the abundances of four elements - Mg, Al, Si, and S. We find that the temperature and emission measure scales with the flares sub-class while the measured abundances show an intermediate FIP bias for the lower A-flares (e.g., A1), while for the higher A-flares, the FIP bias is near unity. To investigate it further, we perform a time-resolved spectral analysis for a sample of the A-class flares and examine the evolution of temperature, emission measure, and abundances. We find that the abundances drop from the coronal values towards their photospheric values in the impulsive phase of the flares, and after the impulsive phase, they quickly return to the usual coronal values. The transition of the abundances from the coronal to photospheric values in the impulsive phase of the flares indicates the injection of fresh unfractionated material from the lower solar atmosphere to the corona due to chromospheric evaporation. However, explaining the quick recovery of the abundances from the photospheric to coronal values in the decay phase of the flare is challenging.

The Cherenkov Telescope Array (CTA) is the future observatory for ground-based imaging atmospheric Cherenkov telescopes. Each telescope will provide a snapshot of gamma-ray induced particle showers by capturing the induced Cherenkov emission at ground level. The simulation of such events provides camera images that can be used as training data for convolutional neural networks (CNNs) to differentiate signals from background events and to determine the energy of the initial gamma-ray events. Pattern spectra are commonly used tools for image classification and provide the distributions of the sizes and shapes of features comprising an image. The application of pattern spectra on a CNN allows the selection of relevant combinations of features within an image. In this work, we generate pattern spectra from simulated gamma-ray images to train a CNN for signal-background separation and energy reconstruction for CTA. We compare our results to a CNN trained with CTA images and find that the pattern spectra-based analysis is computationally less expensive but not competitive with the purely CTA images-based analysis. Thus, we conclude that the CNN must rely on additional features in the CTA images not captured by the pattern spectra.

The MeerKAT Fornax Survey maps the distribution and kinematics of atomic neutral hydrogen gas (HI) in the nearby Fornax galaxy cluster using the MeerKAT telescope. The 12 deg^2 survey footprint covers the central region of the cluster out to ~ Rvir and stretches out to ~ 2 Rvir towards south west to include the NGC 1316 galaxy group. The HI column density sensitivity (3 sigma over 25 km/s) ranges from 5e+19/cm^2 at a resolution of ~ 10" (~ 1 kpc at the 20 Mpc distance of Fornax) down to ~ 1e+18/cm^2 at ~ 1' (~ 6 kpc), and slightly below this level at the lowest resolution of ~ 100" (~ 10 kpc). The HI mass sensitivity (3 sigma over 50 km/s) is 6e+5 Msun. The HI velocity resolution is 1.4 km/s. In this paper we describe the survey design and HI data processing, and we present a sample of six galaxies with long, one-sided, star-less HI tails (of which only one was previously known) radially oriented within the cluster and with measurable internal velocity gradients. We argue that the joint properties of the HI tails represent the first unambiguous evidence of ram pressure shaping the distribution of HI in the Fornax cluster. The disturbed optical morphology of all host galaxies supports the idea that the tails consist of HI initially pulled out of the galaxies' stellar body by tidal forces. Ram pressure was then able to further displace the weakly bound HI and give the tails their present direction, length and velocity gradient.

We use the recently discovered simple photometric parameter Central Intensity Ratio (CIR, Aswathy & Ravikumar 2018) determined for a sample of 57 nearby (z < 0.02) Seyfert galaxies to explore the central features of galaxies and their possible connection with galaxy evolution. The sample of galaxies shows strong anti-correlation between CIR and mass of their central supermassive black holes (SMBH). The SMBH masses of ellipticals are systematically higher for a given CIR value than that for lenticulars and spirals in the sample. However, the correlation between CIR and central velocity dispersion is weak. CIR appears less influenced by the excess flux produced by the central engine in these galaxies, when compared to spectroscopic parameters like velocity dispersion and OIV flux, and proves a fast and reliable tool for estimating central SMBH mass.

The nova super-remnant (NSR) surrounding M31N 2008-12a (12a), the annually erupting recurrent nova (RN), is the only known example of this phenomenon. As this structure has grown as a result of frequent eruptions from 12a, we might expect to see NSRs around other RNe; this would confirm the RN--NSR association and strengthen the connection between novae and type Ia supernovae (SN Ia) as NSRs centered on SN Ia provide a lasting, unequivocal signpost to the single degenerate progenitor type of that explosion. The only previous NSR simulation used identical eruptions from a static white dwarf (WD). In this Paper, we simulate the growth of NSRs alongside the natural growth/erosion of the central WD, within a range of environments, accretion rates, WD temperatures, and initial WD masses. The subsequent evolving eruptions create dynamic NSRs tens of parsecs in radius comprising a low-density cavity, bordered by a hot ejecta pile-up region, and surrounded by a cool high-density, thin, shell. Higher density environments restrict NSR size, as do higher accretion rates, whereas the WD temperature and initial mass have less impact. NSRs form around growing or eroding WDs, indicating that NSRs also exist around old novae with low-mass WDs. Observables such as X-ray and H$\alpha$ emission from the modelled NSRs are derived to aid searches for more examples; only NSRs around high accretion rate novae will currently be observable. The observed properties of the 12a NSR can be reproduced when considering both the dynamically grown NSR and photoionisation by the nova system.

*Context: The optimisation of new multiplex spectrographs (resolution, wavelength range,...), their associated surveys (choice of setup), or their parameterisation pipelines require methods that estimate which wavelengths contain useful information. *Aim: We propose a method that establishes the usefulness (purity & detectability) of an atomic line. We show two applications: a) optimising an instrument, by comparing the number of useful lines at a given setup, and b) optimising the line-list for a given setup by choosing the least blended lines detectable at different signal-to-noise ratios. *Method: The method compares pre-computed synthetic stellar spectra containing all of the elements and molecules with spectra containing the lines of specific elements alone. Then, the flux ratios between the full spectrum and the element spectrum are computed to estimate the line purities. The method identifies automatically (i) the line's central wavelength, (ii) its detectability based on its depth and a given S/N threshold and (iii) its usefulness based on the purity ratio. *Results: We compare the three WEAVE high-resolution setups (Blue: 404-465nm, Green: 473-545nm, Red: 595-685nm), and find that the Green+Red setup both allows one to measure more elements and contains more useful lines. However, there is a disparity in terms of which elements are detected, which we characterise. We also study the performances of R~20 000 and R~6000 spectra covering the entire optical range. Assuming a purity threshold of 60%, we find that the HR setup contains a much wealthier selection of lines, for any of the considered elements, whereas the LR has a "loss" of 50 to 90% of the lines even for higher S/N. *Conclusions: The method provides a diagnostic of where to focus to get the most out of a spectrograph, and is easy to implement for future instruments, or for pipelines that require line masks.

We present an overview of the Southern-sky MWA Rapid Two-metre (SMART) pulsar survey that exploits the MWA's large field of view and voltage capture system to survey the sky south of 30 degree in declination for pulsars and fast transients in the 140-170 MHz band. The survey is enabled by the advent of the Phase II MWA's compact configuration, which offers an enormous efficiency in beam-forming and processing costs, thereby making an all-sky survey of this magnitude tractable with the MWA. Even with the long dwell times of the survey (4800 s), data collection can be completed in < 100 hours of telescope time, while still retaining the ability to reach a limiting sensitivity of ~2-3 mJy. Each observation is processed to generate ~5000-8000 tied-array beams that tessellate the full ~610 square degree field of view, which are then processed to search for pulsars. The voltage-capture recording allows a multitude of post hoc processing options including the reprocessing of data for higher time resolution. Due to the substantial computational cost in pulsar searches at low frequencies, processing is undertaken in multiple passes: in the first pass, a shallow survey is performed, where 10 minutes of each observation is processed, reaching about one-third of the full search sensitivity. Here we present the system overview and initial results. Further details including first pulsar discoveries and a census of low-frequency detections are presented in a companion paper. Future plans include deeper searches to reach the full sensitivity and acceleration searches to target binary and millisecond pulsars. Simulation analysis forecasts ~300 new pulsars upon the completion of full processing. The SMART survey will also generate a complete digital record of the low-frequency sky, which will serve as a valuable reference for future pulsar searches planned with the low-frequency Square Kilometre Array.

In Paper I, we presented an overview of the Southern-sky MWA Rapid Two-metre (SMART) survey, including the survey design and search pipeline. While the combination of MWA's large field-of-view and the voltage capture system brings a survey speed of ~450 square degrees per hour, the survey progression relies on the availability of compact configuration of the Phase II array. Over the past few years, by taking advantage of multiple windows of opportunity when the compact configuration was available, we have advanced the survey to 75% completion. To date, about 10% of the data collected thus far have been processed for a first-pass search, where 10 minutes of observation is processed for dispersion measures out to 250 ${\rm pc\,cm^{-3}}$, to realise a shallow survey for long-period pulsars. The ongoing analysis has led to two new pulsar discoveries, as well as an independent discovery and a rediscovery of a previously incorrectly characterised pulsar, all from ~3% of the data for which candidate scrutiny is completed. Here we describe the strategies for further detailed follow-up including improved sky localisation and convergence to timing solution, and illustrate them using example pulsar discoveries. The processing has also led to re-detection of 120 pulsars in the SMART observing band, bringing the total number of pulsars detected to date with the MWA to 180, and these are used to assess the search sensitivity of current processing pipelines. The planned second-pass (deep survey) processing is expected to yield a three-fold increase in sensitivity for long-period pulsars, and a substantial improvement to millisecond pulsars by adopting optimal de-dispersion plans. The SMART survey will complement the highly successful Parkes High Time Resolution Universe survey at 1.2-1.5 GHz, and inform future large survey efforts such as those planned with the low-frequency Square Kilometre Array.

In this work we present an observational technique and a detailed analysis of the stellar interferograms produced by three bright stars: Betelgeuse, Rigel and Sirius. It is shown that the atmospheric turbulence is responsible for the reduction of the long-exposure fringe visibility of the obtained interference patterns. By using different baselines in our interferometer, we are able to distinguish the decay of the visibility with the baseline, how different parameters such us the diameter of the holes in our interferometer or their distribution affects the pattern, and to measure the turbulence with the estimation of the Fried parameter r0. The work and methodology are presented as a method for postgraduate students that targets practical learning of optical interferometry in astronomy and how it is affected by several causes, such as the atmospheric turbulence.

J1004+4112 is a lensed quasar for which the first broad emission line profile deformations due to microlensing were identified. Detailed interpretations of these features have nevertheless remained controversial. Based on 15 spectra obtained from 2003 to 2018, we revisit the microlensing effect that distorts the CIV broad emission line profile. We show that the microlensing-induced line profile distortions in image A, although variable, are remarkably similar over a period of 15 years. They are characterized by a strong magnification of the blue part of the line profile, a strong demagnification of the red part of the line profile, and a small-to-negligible demagnification of the line core. We used the microlensing effect to constrain the broad emission-line region (BLR) size, geometry, and kinematics. For this purpose, we modeled the deformation of the emission lines considering three simple, representative BLR models: a Keplerian disk, an equatorial wind, and a biconical polar wind, with various inclinations with respect to the line of sight. We find that the observed magnification profile of the CIV emission line can be reproduced with the simple BLR models we considered, without the need for more complex BLR features. The magnification appears dominated by the position of the BLR with respect to the caustic network -- and not by the velocity-dependent size of the BLR. The favored models for the CIV BLR are either the Keplerian disk or the equatorial wind, depending on the orientation of the BLR axis with respect to the caustic network. We also find that the polar wind model can be discarded. We measured the CIV BLR half-light radius as $r_{1/2} = 2.8^{+2.0}_{-1.7}$ light-days. This value is smaller than the BLR radius expected from the radius-luminosity relation derived from reverberation mapping, but it is still in reasonable agreement given the large uncertainties.

It has been suggested recently that the appparent accelerated expansion of the universe could be explained by a bias in the SNIa measurements. Such events indeed occur mainly in overdense regions, where matter is located, and whose dynamics can perhaps not been considered as representative of the one of the universe. In this article, we develop a model to investigate in more detail the effect of this bias. This model depends on one single parameter, related to the void fraction of space, and leads to simple analytical relations. We in particular determine the average metric tensor in overdense regions, and deduce that the scale factor and the rate at which time progresses in such regions differ significantly from the corresponding values expected on the average space. We then quantitatively deduce how redshift and luminosity distance measurements are affected by the bias, taking into account the perturbation of the metric tensor. Using a value for the void fraction corresponding to the order of magnitude found in the literature, we show that the model is able to predict a distance modulus versus redshift relation being in excellent aggreement with the one corresponding to a universe characterized by $\Omega_{m,0} = 0.3$ and $\Omega_{\Lambda,0} = 0.7$.

The Be stars display variable optical emission lines originating in the circumstellar disc. Here we analyse high resolution spectroscopic observations of Be stars and the distance between the peaks of H-alpha, H-beta, and H-gamma emission lines ($\Delta V_\alpha$, $\Delta V_\beta$, and $\Delta V_\gamma$ respectively). Combining published data, spectra from the ELODIE archive (obtained in the period 1998 -- 2003) and Rozhen spectra (obtained 2015 -- 2023) of 93 Be stars, we find a set of relations connecting $\Delta V_\alpha$, $\Delta V_\beta$ and $\Delta V_\gamma$. They are effective for $30 \le \Delta V_\alpha \le 500$ km s$^{-1}$, $80 \le \Delta V_\beta \le 600$ km s$^{-1}$, and $40 \le \Delta V_\gamma \le 300$ km s$^{-1}$. The new equations are in the form $y=ax + b$ and are valid for a wider velocity range than in previous studies.

In the nearby universe jets from AGN are observed to have a dramatic impact on their surrounding extragalactic environment. Their effect at the `cosmic noon' (z>1.5), the epoch when star formation and AGN activity peak, is instead much less constrained. Here we present a study of the giant (750 kpc) radio galaxy 103025+052430 located at the centre of a protocluster at redshift z=1.7, with a focus on its interaction with the external medium. We present new LOFAR observations at 144 MHz, which we combine with VLA 1.4 GHz and 0.5-7 keV Chandra archival data. The new map at 144 MHz confirms that the source has a complex morphology, possibly consistent with the `hybrid morphology' classification. The large size of the source gave us the possibility to perform a resolved radio spectral index analysis, a very unique opportunity for a source at such high redshift. This reveals a tentative flattening of the radio spectral index at the edge of the backflow in the Western lobe, which might be indicating plasma compression. The spatial coincidence between this region and the thermal X-ray bubble C suggests a causal connection between the two. Contrary to previous estimates for the bright X-ray component A, we find that inverse Compton scattering between the radio-emitting plasma of the Eastern lobe and the CMB photons can account for a large fraction (~45%-80%) of its total 0.5-7 keV measured flux. Finally, the X-ray bubble C, which is consistent with a thermal origin, is found to be significantly overpressurised with respect to the ambient medium. This suggests that it will tend to expand and release its energy in the surroundings, contributing to the overall intracluster medium heating. Overall, 103025+052430 gives us the chance to investigate the interaction between AGN jets and the surrounding gas in a system that is likely the predecessor of the rich galaxy clusters we all well know at z=0.

Detailed knowledge of the different classes of stellar orbits that can be accommodated in a given galactic potential is a prerequisite when building self-consistent models using for instance the Schwazschild technique. Furthermore, observational properties of galaxies depend on what these classes of orbits are and on the presence of chaos in the systems. In the realistic case in which the starting point for modeling is not a gravitational potential, but an observed density distribution, we will require a gravitational theory to make the connection between the stars that we see and the movement these stars may be having. The argument can be turned upside down: understanding what orbits may be allowed by each gravitational theory may give us a greater insight on what these theories are and on how we can test them. Our aim is thus to understand novel properties of orbits that are predicted by the latest extension of the MOND phenomenology into the relativistic world. We integrated orbits numerically in a fixed density distribution. The potential required for such integration was obtained also numerically by assuming different gravitational models. We find that thanks to the presence of a mass term in the field equations, the theory can allocate new classes of orbits that do not exist in Newtonian gravity nor standard MOND. We discuss consequences that these new families of orbits can have in non-linear cosmological structure formation as well as explore a possible alternative model for galactic structure based on them.

The Tucana-Horologium association is one of the closest young stellar groups to the Sun and despite the close proximity its age is still debated in the literature. We take advantage of the state-of-the-art astrometry delivered by the third data release of the Gaia space mission combined with precise radial velocity measurements obtained from high-resolution spectroscopy to investigate the dynamical age of the association. We perform an extensive traceback analysis using a combination of different samples of cluster members, metrics to evaluate the minimum size of the association in the past and models for the galactic potential to integrate the stellar orbits back in time. The dynamical age of $38.5^{+1.6}_{-8.0}$ Myr that we derive in this paper is consistent with the various age estimates obtained from isochrone fitting in the literature (30-50 Myr) and reconciles, for the first time, the dynamical age of the Tucana-Horologium association with the age obtained from lithium depletion ($\sim40$ Myr). Our results are independent from stellar models and represent one more step towards constructing a self-consistent age scale for the young stellar groups of the Solar neighbourhood based on the 3D space motion of the stars.

We study the impact of pions in simulations of neutron star mergers and explore the impact on gravitational-wave observables. We model charged and neutral pions as a non-interacting Boson gas with a chosen, constant effective mass. We add the contributions of pions, which can occur as a condensate or as a thermal population, to existing temperature and composition dependent equations of state. Compared to the models without pions, the presence of a pion condensate decreases the characteristic properties of cold, non-rotating neutron stars such as the maximum mass, the radius and the tidal deformability. We conduct relativistic hydrodynamical simulations of neutron star mergers for these modified equations of state models and compare to the original models, which ignore pions. Generally, the inclusion of pions leads to a softening of the equation of state, which is more pronounced for smaller effective pion masses. We find a shift of the dominant postmerger gravitational-wave frequency by up to 150~Hz to higher frequencies and a reduction of the threshold binary mass for prompt black-hole formation by up to 0.07~$M_\odot$. We evaluate empirical relations between the threshold mass or the dominant postmerger gravitational-wave frequency and stellar parameters of nonrotating neutron stars. These relations are constructed to extract these stellar properties from merger observations and are built based on large sets of equation of state models which do not include pions. Comparing to our calculations with pions, we find that these empirical relations remain valid to good accuracy, which justifies their use although they neglect a possible impact of pions. We also address the mass ejection by neutron star mergers and observe a moderate enhancement of the ejecta mass by a few ten per cent. (abridged)

We present BISTRO Survey 850 {\mu}m dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965).

Context. We are operating an elastic LIDAR for the monitoring of atmospheric conditions during regular observations of the MAGIC Telescopes. Aims. We present and evaluate methods to convert aerosol extinction profiles, obtained with the LIDAR, into corrections of the reconstructed gamma-ray event energy and Instrument Response Functions of Imaging Atmospheric Cherenkov Telescopes. Methods. We assess the performance of these correction schemes with almost seven years of Crab Nebula data taken by the MAGIC Telescopes under various zenith angles and different aerosol extinction scenarios of Cherenkov light. Results. The methods enable the reconstruction of data taken under non-optimal atmospheric conditions with aerosol transmissions down to around 0.65 with systematic uncertainties comparable to those for data taken under optimal conditions. For the first time, the correction of data affected by clouds has been included in the assessment. The data can also be corrected when the transmission is lower than 0.65, but the results are less accurate and suffer from larger systematics.

The evolution of many astrophysical systems is dominated by the interaction between matter and radiation such as photons or neutrinos. The dynamics can be described by the evolution equations of radiation hydrodynamics in which reactions between matter particles and radiation quanta couples the hydrodynamic equations to those of radiative transfer, see Munier & Weaver (1986a), Munier & Weaver (1986b).. The numerical treatment has to account for their potential stiffness (e.g., in optically thick environments). In this article, we will present a new method to numerically integrate these equations in a stable way by using minimally implicit Runge-Kutta methods. With these methods, the inversion of the implicit operator can be done analytically. We also take into account the physical behavior of the evolved variables in the limit of the stiff regime. We will show the results of applying this method to the reactions between neutrinos and matter in core-collapse supernovae simulations.

The mid-infrared (IR) regime is well suited to directly detect the thermal signatures of exoplanets in our solar neighborhood. The NEAR experiment: demonstration of high-contrast imaging (HCI) capability at ten microns, can reach sub-mJy detection sensitivity in a few hours of observation time, which is sufficient to detect a few Jupiter mass planets in nearby systems. One of the big limitations for HCI in the mid-IR is thermal sky-background. In this work, we show that precipitate water vapor (PWV) is the principal contributor to thermal sky background and science PSF quality. In the presence of high PWV, the HCI performance is significantly degraded in the background limited regime.

For deep-space mission design, the gravity of the Sun and the Moon can be first considered and utilized. Their gravity can provide the energy change for launching spacecraft and retrieving spacecraft as well as asteroids. Regarding an asteroid retrieval mission, it can lead to mitigation of asteroid hazards and easy exploration and exploitation of the asteroid. This paper discusses the application of the Sun-driven lunar swingby sequence for asteroid missions. To characterize the capacity of this technique is not only interesting in terms of dynamic insights but also non-trivial for trajectory design. This paper elucidates the capacity of a Sun-driven lunar swingby sequence with the help of the "Swingby-Jacobi" graph. The capacity can be represented by a range of Jacobi integral that encloses around 660 asteroids currently cataloged. To facilitate trajectory design, a database of Sun-perturbed Moon-to-Moon transfers including multi-revolution cases is generated and employed. Massive trajectory options for spacecraft launch and asteroid capture can then be explored and optimized. Finally, a number of asteroid flyby, rendezvous, sample-return, and retrieval mission options enabled by the proposed technique are obtained.

The KAGRA cryogenic gravitational-wave observatory has begun joint observation with the worldwide gravitational waves detector network. Precise calibration of the detector response is essential for the parameter estimation of gravitational wave sources. The photon calibrator is the main calibrator in LIGO, Virgo and KAGRA, and we used this calibrator in joint observation 3 on 2020 April with GEO600 in Germany. KAGRA improved the system for joint observation 3 with three unique points: high laser power, power stabilization system, and remote beam position control. KAGRA employs the 20 W laser and divides it into two beams injected on the mirror surface. By using a high-power laser, we can calibrate the response at the kHz region. To control the power of each laser independently, we also installed an optical follower servo for each beam power stabilization. By controlling the optical path of the photon calibrater beam positions with pico-motors, we were able to characterize the rotation response of the detector. We also installed a telephoto camera and QPD to monitor beam position and controlled beam position to optimize mirror response. In this paper, we discussed the statistical error with the result of the relative power noise measurement. We also discussed systematic errors about the power calibration model of photon calibrator and simulation of elastic deformation effect with the finite element analysis.

In this paper, we present data from the Blind Ultra-Deep HI Environmental Survey (BUDHIES), which is a blind 21-cm HI spectral line imaging survey undertaken with the Westerbork Synthesis Radio Telescope (WSRT). Two volumes were surveyed, each with a single pointing and covering a redshift range of 0.164 < z < 0.224. Within these two volumes, this survey targeted the clusters Abell 963 and Abell 2192, which are dynamically different and offer unique environments to study the process of galaxy evolution within clusters. With an integration time of 117x12h on Abell 963 and 72x12h on Abell 2192, a total of 166 galaxies were detected and imaged in HI. While the clusters themselves occupy only 4 per cent of the 73,400 Mpc$^3$ surveyed by BUDHIES, most of the volume consists of large-scale structures in which the clusters are embedded, including foreground and background overdensities and voids. We present the data processing and source detection techniques and counterpart identification based on a wide-field optical imaging survey using the Isaac Newton Telescope (INT) and deep ultra-violet GALEX imaging. Finally, we present HI and optical catalogues of the detected sources as well as atlases of their global HI properties, which include integrated column density maps, position-velocity diagrams, global HI profiles, and optical and UV images of the HI sources.

We utilize high-resolution cosmological simulations to reveal that high-redshift galaxies tend to undergo a robust `wet compaction' event when near a `golden' stellar mass of $\sim 10^{10} M_{\odot}$. This is a gaseous shrinkage to a compact star-forming phase, a `blue nugget' (BN), followed by central quenching of star formation to a compact passive stellar bulge, a `red nugget' (RN), and a buildup of an extended gaseous disc and ring. Such nuggets are observed at cosmic noon and seed today's early-type galaxies. The compaction is triggered by a drastic loss of angular momentum due to, e.g., wet mergers, counter-rotating cold streams, or violent disc instability. The BN phase marks drastic transitions in the galaxy structural, compositional and kinematic properties. The transitions are from star-forming to quenched inside-out, from diffuse to compact with an extended disc-ring and a stellar envelope, from dark matter to baryon central dominance, from prolate to oblate stellar shape, from pressure to rotation support, from low to high metallicity, and from supernova to AGN feedback. The central black hole growth, first suppressed by supernova feedback when below the golden mass, is boosted by the compaction, and the black hole keeps growing once the halo is massive enough to lock in the supernova ejecta.

Scientific articles published prior to the "age of digitization" in the late 1990s contain figures which are "trapped" within their scanned pages. While progress to extract figures and their captions has been made, there is currently no robust method for this process. We present a YOLO-based method for use on scanned pages, after they have been processed with Optical Character Recognition (OCR), which uses both grayscale and OCR-features. We focus our efforts on translating the intersection-over-union (IOU) metric from the field of object detection to document layout analysis and quantify "high localization" levels as an IOU of 0.9. When applied to the astrophysics literature holdings of the NASA Astrophysics Data System (ADS), we find F1 scores of 90.9% (92.2%) for figures (figure captions) with the IOU cut-off of 0.9 which is a significant improvement over other state-of-the-art methods.

There may exist extended configurations in the dark matter sector that are analogues of structures in the visible sector. In this work, we explore non-topological solitonic configurations, specifically Q-balls, and study when they may form macroscopic astrophysical structures and what their distinct characteristics might be. We study in some detail theoretical bounds on their sizes and constraints on the underlying parameters, based on criteria for an astrophysical Q-ball's existence, gravitational stability and viability of solutions. Following this path, one is able to obtain novel limits on astrophysical Q-ball sizes and their underlying parameters. We also explore the gravitational lensing features of different astrophysical Q-ball profiles, which are more general than the simple thin-wall limit. It is seen that the magnification characteristics may be very distinct, depending on the actual details of the solution, even for astrophysical Q-balls having the same size and mass. Assuming that such astrophysical Q-balls may form a small component of the dark matter in the universe, we place limits on this fraction from the gravitational microlensing surveys EROS-2, OGLE-IV and the proposed future survey WFIRST. Exploring various astrophysical Q-ball profiles and sizes, it is found that while for smaller masses, the dark matter fraction comprising astrophysical Q-balls is at most sub-percent, for larger masses, it may be significantly higher.

High temporal and high spatial resolution geoelectric field models of two M\"ants\"al\"a, Finnish pipeline GIC intervals that occurred within the 7-8 September, 2017 geomagnetic storm have been made. The geomagnetic measurements with 10 s sampling rate of 28 IMAGE ground magnetometers distributed over the north Europe (from $52.07^\circ$ to $69.76^\circ$ latitude) are the bases for the study. A GeoElectric Dynamic Mapping (GEDMap) code was developed for this task. GEDMap considers 4 different methods of interpolation and allows a grid of $0.05^\circ$ (lat.)$\times 0.2^\circ$ (lon.) spatial scale resolution. The geoelectric field dynamic mapping output gives both spatial and temporal variations of the magnitude and direction of fields. The GEDMap results show very rapid and strong variability of geoelectric field and the extremely localized peak enhancements. The magnitude of geoelectric fields over M\"ants\"al\"a at the time of the two GIC peaks were 279.7 mV/km and 336.9 mV/km. The comparison of the GIC measurements in M\"ants\"al\"a and our modeling results show very good agreement with a correlation coefficient higher than 0.8. It is found that the auroral electrojet geoelectric field has very rapid changes in both magnitude and orientation causing the GICs. It is also shown that the electrojet is not simply oriented in the east-west direction. It is possible that even higher time resolution base magnetometer data of 1 s will yield even more structure, so this will be our next effort.

Pulsar timing arrays (PTA) are a promising probe to the cosmologically novel nanohertz gravitational wave (GW) regime through the stochastic GW background. In this work, we consider subluminal GW modes as a possible source of correlations in a PTA, utilizing the public code PTAfast and the 12.5 years correlations data by NANOGrav. Our results show no evidence in support of tensor- or vector-induced GW correlations in the data, and that vector correlations are disfavored. This places an upper bound to the graviton mass, $m_{\rm g} \lesssim 10^{-22}$ eV, characteristic of the PTA GW energy scale.

In this work we analyze the observational properties of incompressible relativistic fluid spheres with and without thin-shells, when surrounded by thin accretion disks. We consider a set of six configurations with different combinations of the star radius $R$ and the thin-shell radius $r_\Sigma$ to produce solutions with neither thin-shells nor light-rings, with either of those features, and with both. Furthermore, we consider three different models for the intensity profile of the accretion disk, based on the Gralla-Lupsasca-Marrone (GLM) disk model, for which the peaks of intensity occur at the Innermost Stable Circular Orbit (ISCO), the Light-Ring (LR), and the center of the star. The observed images and intensity profiles for an asymptotic observer are produced using a Mathematica-based ray-tracing code. Our results indicate that, in the absence of a light-ring, the presence of a thin-shell produces a negligible effect in the observational properties of the stars. However, when the spacetime features a light-ring, the portion of the mass of the star that is stored in the thin-shell has a strong effect on its observational properties, particularly in the magnitude of the central gravitational redshift effect responsible for causing a central shadow-like dimming in the observed images. A comparison with the Schwarzschild spacetime is also provided and the most compact configurations are shown to produce observational imprints similar to those of black-hole solutions, with subtle qualitative differences, most notably extra secondary image components that decrease the radius of the shadow and are potentially observable.

In this study, cosmological models are considered, where dark matter and dark energy are coupled and may exchange energy through non-gravitational interactions with one other. These interacting dark energy (IDE) models have previously been introduced to address problems with the standard $\Lambda$CDM model of cosmology (which include the coincidence problem, Hubble tension and $S_8$ discrepancy). However, conditions ensuring positive energy densities have often been overlooked. Assuming two different linear dark energy couplings, $Q = \delta H \rho_{\rm{de}}$ and $Q = \delta H \rho_{\rm{dm}}$, we find that negative energy densities are inevitable if energy flows from dark matter to dark energy (iDMDE regime) and that consequently, we should only seriously consider models where energy flows from dark energy to dark matter (iDEDM regime). To additionally ensure that these models are free from early time instabilities, we need to require that dark energy is in the `phantom' ($\omega<-1$) regime. This has the consequence that model $Q=\delta H \rho_{\rm{dm}}$ will end with a future big rip singularity, while $Q = \delta H \rho_{\rm{de}}$ may avoid this fate with the right choice of cosmological parameters.

Gravitational wave science is a new and rapidly expanding field of observational astronomy. Multimessenger observations of the binary neutron star merger GW170817 have provided some iconic results including the first gravitational-wave standard-siren measurement of the Hubble constant, opening up a new way to probe cosmology. The majority of the compact binary sources observed in gravitational waves are however without bright electromagnetic counterparts. In these cases, one can fall back on the ``dark standard siren'' approach to include information statistically from potential host galaxies. For such a measurement, we need to be cautious about all possible sources of systematic errors. In this paper, we begin to study the possible errors coming from the galaxy catalogue sector, and in particular, look into the effect of galaxy redshift uncertainties for the cases where these are photometry-based. We recalculate the dark standard siren Hubble constant using the latest GWTC-3 events and associated galaxy catalogues, with different galaxy redshift uncertainty models, namely, the standard Gaussian, a modified Lorentzian, and no uncertainty at all. We find that not using redshift uncertainties at all can lead to a potential bias comparable with other potential systematic effects previously considered for the GWTC-3 $H_0$ measurement (however still small compared to the overall statistical error in this measurement). The difference between different uncertainty models leads to small differences in the results for the current data; their impact is much smaller than the current statistical errors and other potential sources of systematic errors which have been considered in previous robustness studies.

We present a comprehensive simulation of the spatial acquisition of optical links for the LISA mission in the in-field pointing architecture, where a fast pointing mirror is used to move the field-of-view of the optical transceiver, which was studied as an alternative scheme to the baselined telescope pointing architecture. The simulation includes a representative model of the far-field intensity distribution and the beam detection process using a realistic detector model, and a model of the expected platform jitter for two alternative control modes with different associated jitter spectra. For optimally adjusted detector settings and accounting for the actual far-field beam profile, we investigate the dependency of acquisition performance on the jitter spectrum and the track-width of the search spiral, while scan speed and detector integration time are varied over several orders of magnitude. Results show a strong dependency of the probability for acquisition failure on the width of the auto-correlation function of the jitter spectrum, which we compare to predictions of analytical models. Depending on the choice of scan speed, three different regimes may be entered which differ in failure probability by several orders of magnitude. We then use these results to optimize the acquisition architecture for the given jitter spectra with respect to failure rate and overall duration, concluding that the full constellation could be acquired on average in less than one minute. Our method and findings can be applied to any other space mission using a fine-steering mirror for link acquisition.

We propose a modification of no-scale supergravity models which incorporates sgoldstino stabilization and supersymmetry (SUSY) breaking with a tunable cosmological constant by introducing a Kahler potential which yields a kinetic pole of order one. The resulting scalar potential may develop an inflection point, close to which an inflationary period can be realized for subplanckian field values consistently with the observational data. For central value of the spectral index ns, the necessary tuning is of the order of 10^-6, the tensor-to-scalar ratio is tiny whereas the running of ns is around -3x10^-3. Our proposal is compatible with high-scale SUSY and the results of LHC on the Higgs boson mass.