Papers for Friday, Aug 19 2022

Context. The rapidly improving quality and resolution of both low surface brightness observations and cosmological simulations of galaxies enables an increasingly thorough investigation of the imprints of the formation history in the outer, unrelaxed regions of galaxies, and a direct comparison to another tracer of galaxy formation, the internal kinematics. Aims. Using the state-of-the-art hydrodynamical cosmological simulation Magneticum Pathfinder, we identify tidal tails, shells, streams, and satellite planes, and connect their existence to the amount of rotational support and the formation histories of the host galaxies. Methods. Tidal features are visually classified from a three-dimensional rendering of the simulated galaxies by several scientists independently. Only features that were identified by at least half of the participating individuals are considered as existing features. The results are compared to observations of the MATLAS survey. Results. Shells are preferentially found around kinematically slowly rotating galaxies in both simulations and observations, while streams can be found around all kind of galaxies with a slightly higher probability to be present around less rotationally supported galaxies. Tails and satellite planes, however, appear independently of the internal kinematics of the central galaxy, indicating that they are formed through processes that have not (yet) affected the internal kinematics. Conclusions. As shells are formed through radial merger events while streams are remnants of circular merger infall, this suggests that the orbital angular momentum of the merger event could play a more crucial role in transforming the host galaxy than previously anticipated. The existence of a shell around a given slow rotator can further be used to distinguish the radial merger formation scenario from other formation pathways of slow rotators.

Low- and intermediate-mass ($\rm 0.8~M_\odot < M < 8~M_\odot$) stars that evolve into planetary nebulae (PNe) play an important role in tracing and driving Galactic chemical evolution. Spectroscopy of PNe enables access to both the initial composition of their progenitor stars and products of their internal nucleosynthesis, but determining accurate ionic and elemental abundances of PNe requires high-quality optical spectra. We obtained new optical spectra of eight highly-extincted PNe with limited optical data in the literature using the Low Resolution Spectrograph 2 (LRS2) on the Hobby-Eberly Telescope (HET). Extinction coefficients, electron temperatures and densities, and ionic and elemental abundances of up to 11 elements (He, N, O, Ne, S, Cl, Ar, K, Fe, Kr, and Xe) are determined for each object in our sample. Where available, astrometric data from Gaia eDR3 is used to kinematically characterize the probability that each object belongs to the Milky Way's thin disk, thick disk, or halo. Four of the PNe show kinematic and chemical signs of thin disk membership, while two may be members of the thick disk. The remaining two targets lack Gaia data, but their solar O, Ar, and Cl abundances suggest thin disk membership. Additionally, we report the detection of broad emission features from the central star of M 3-35. Our results significantly improve the available information on the nebular parameters and chemical compositions of these objects, which can inform future analyses.

Future missions like Roman, HabEx, and LUVOIR will directly image exoplanets in reflected light. While current near infrared direct imaging searches are only sensitive to young, self-luminous planets whose brightness is independent of their orbital phase, reflected light imaging will reveal changes in planet brightness over the course of an orbit due to phase variations. One of the first objectives will be determining the planet's orbit via astrometry, the projected position of the planet with respect to its host star in the sky plane. We show that phase variations can significantly improve the accuracy and precision of orbital retrieval with two or three direct images. This would speed up the classification of exoplanets and improve the efficiency of subsequent spectroscopic characterization. We develop a forward model to generate synthetic observations of the two dimensional position of the planet with respect to its host star on the sky plane, and the planet/star flux ratio. Synthetic data are fitted with Keplerian orbits and Henyey-Greenstein phase variations to retrieve orbital and phase parameters. For astrometric uncertainties of 0.01 AU in projected separation and flux ratio uncertainties of 10^-12, using photometry in orbit retrieval improves the accuracy of semi-major axis by 47% for two epochs and 61% for three epochs if the phase curves have a known shape, but unknown amplitude. In the more realistic scenario where the shape and amplitude of the phase curve are a priori unknown, photometry still improves accuracy by 16% for two epochs and 50% for three.

VHS J1256-1257AB is an ultracool dwarf binary that hosts a wide-separation planetary-mass companion that is a key target of the JWST Exoplanet Early Release Science program. Using Keck adaptive optics imaging and aperture masking interferometry, we have determined the host binary's orbit ($a=1.96\pm0.03$ au, $P=7.31\pm0.02$ yr, $e=0.8826^{+0.0025}_{-0.0024}$) and measured its dynamical total mass ($0.141\pm0.008$ $M_{\odot}$). This total mass is consistent with VHS J1256-1257AB being a brown dwarf binary or pair of very low-mass stars. In addition, we measured the orbital motion of VHS J1256-1257 b with respect to the barycenter of VHS J1256-1257AB, finding that the wide companion's orbit is also eccentric ($e=0.73^{+0.09}_{-0.10}$), with a mutual inclination of $116^{\circ}\pm16^{\circ}$ with respect to the central binary. This orbital architecture is consistent with VHS J1256-1257 b attaining a significant mutual inclination through dynamical scattering and thereafter driving Kozai-Lidov cycles to pump the eccentricity of VHS J1256-1257AB. We derive a cooling age of $140\pm20$ Myr for VHS J1256-1257AB from low-mass stellar/substellar evolutionary models. At this age, the luminosity of VHS J1256-1257 b is consistent with both deuterium-inert and deuterium-fusing evolutionary tracks. We thus find a bimodal probability distribution for the mass of VHS J1256-1257 b, either $11.8\pm0.2$ $M_{\rm Jup}$ or $16\pm1$ $M_{\rm Jup}$, from Saumon & Marley (2008) hybrid models. Future spectroscopic data to measure isotopologues such as HDO and CH$_3$D could break this degeneracy and provide a strong test of substellar models at the deuterium-fusion mass boundary.

AR Scorpii (AR Sco) is the only radio-pulsing white dwarf known to date. It shows a broad-band spectrum extending from radio to X-rays whose luminosity cannot be explained by thermal emission from the system components alone, and is instead explained through synchrotron emission powered by the spin-down of the white dwarf. We analysed NTT/ULTRACAM, TNT/ULTRASPEC, and GTC/HiPERCAM high-speed photometric data for AR Sco spanning almost seven years and obtained a precise estimate of the spin frequency derivative, now confirmed with 50-sigma significance. Using archival photometry, we show that the spin down rate of P/Pdot = 5.6e6 years has remained constant since 2005. As well as employing the method of pulse-arrival time fitting used for previous estimates, we also found a consistent value via traditional Fourier analysis for the first time. In addition, we obtained optical time-resolved spectra with WHT/ISIS and VLT/X-shooter. We performed modulated Doppler tomography for the first time for the system, finding evidence of emission modulated on the orbital period. We have also estimated the projected rotational velocity of the M-dwarf as a function of orbital period and found that it must be close to Roche lobe filling. Our findings provide further constraints for modelling this unique system.

The streaming instability is a promising mechanism to induce the formation of planetesimals. Nonetheless, this process has been found in previous studies to require either a dust-to-gas surface density ratio or a dust size that is enhanced compared to observed values. Employing two-dimensional global simulations of protoplanetary disks, we show that the vertical shear instability and the streaming instability in concert can cause dust concentration that is sufficient for planetesimal formation for lower surface density ratios and smaller dust sizes than the streaming instability in isolation, and in particular under conditions that are consistent with observational constraints. This is because dust overdensities forming in pressure bumps induced by the vertical shear instability act as seeds for the streaming instability and are enhanced by it. While our two-dimensional model does not include self-gravity, we find that strong dust clumping and the formation (and dissolution) of gravitationally unstable overdensities can be robustly inferred from the evolution of the maximum or the mean dust-to-gas volume density ratio. The vertical shear instability puffs up the dust layer to an average mid-plane dust-to-gas density ratio that is significantly below unity. We therefore find that reaching a mid-plane density ratio of one is not necessary to trigger planetesimal formation via the streaming instability when it acts in unison with the vertical shear instability.

Actively accreting supermassive black holes significantly impact the evolution of their host galaxies, truncating further star formation by expelling large fractions of gas with wide-angle outflows. The X-ray band is key to understanding how these black hole winds affect their environment, as the outflows have high temperatures ($\sim$10$^{5-8}$K). We have developed a Bayesian framework for characterizing Active Galactic Nuclei (AGN) outflows with an improved ability to explore parameter space and perform robust model selection. We applied this framework to a new 700 ks and an archival 315 ks Chandra High Energy Transmission Gratings observation of the Seyfert galaxy NGC 4051. We have detected six absorbers intrinsic to NGC 4051. These wind components span velocities from 400 km s$^{-1}$ to 30,000 km s$^{-1}$. We have determined that the most statistically significant wind component is purely collisionally ionized, which is the first detection of such an absorber. This wind has $T\approx10^7$ K and $v\approx880$ km s$^{-1}$ and remains remarkably stable between the two epochs. Other slow components also remain stable across time. Fast outflow components change their properties between 2008 and 2016, suggesting either physical changes or clouds moving in and out of the line of sight. For one of the fast components we obtain one of the tightest wind density measurements to date, log $n/$[cm$^{-3}$]=13.0$^{+0.01}_{-0.02}$, and determine that it is located at $\sim$240 gravitational radii. The estimated total outflow power surpasses 5% of the bolometric luminosity (albeit with large uncertainties) making it important in the context of galaxy-black hole interactions.

We review the history of spacecraft encounters with comets, concentrating on those that took place in the recent past, since the publication of the Comets II book. This includes the flyby missions Stardust and Deep Impact, and their respective extended missions, the Rosetta rendezvous mission, and serendipitous encounters. While results from all of these missions can be found throughout this book, this chapter focuses on the questions that motivated each mission, the technologies that were required to answer these questions, and where each mission opened new areas to investigate. There remain a large number of questions that will require future technologies and space missions to answer; we also describe planned next steps and routes forward that may be pursued by missions that have yet to be selected, and eventually lead to cryogenic sample return of nucleus ices for laboratory study.

We present the details of a strong lens model of SMACS J0723.3-7327, which was made public as part of the data and high level science products (HLSP) release of the RELICS HST treasury program (Reionization Lensing Cluster Survey; GO-14096, PI: Coe). The model products were made available on the Mikulski Archive for Space Telescopes (MAST) via 10.17909/T9SP45 in 2017. Here, we provide the list of constraints that were used in the HST-based RELICS lens model, as well as other information related to our modeling choices, which were not published with the data and HLSP release. This model was computed with Lenstool, used multiple images of 8 sources, with no spectroscopic redshifts. The image plane RMS was 0".58.

We present new absolute proper motion measurements for the Arches and Quintuplet clusters, two young massive star clusters near the Galactic center. Using multi-epoch HST observations, we construct proper motion catalogs for the Arches ($\sim$35,000 stars) and Quintuplet ($\sim$40,000 stars) fields in ICRF coordinates established using stars in common with the Gaia EDR3 catalog. The bulk proper motions of the clusters are measured to be ($\mu_{\alpha*}$, $\mu_{\delta}$) = (-0.80 $\pm$ 0.032, -1.89 $\pm$ 0.021) mas/yr for the Arches and ($\mu_{\alpha*}$, $\mu_{\delta}$) = (-0.96 $\pm$ 0.032, -2.29 $\pm$ 0.023) mas/yr for the Quintuplet, achieving $\sim$5x higher precision than past measurements. We place the first constraints on the properties of the cluster orbits that incorporate the uncertainty in their current line-of-sight distances. The clusters will not approach closer than $\sim$25 pc to SgrA*, making it unlikely that they will inspiral into the Nuclear Star Cluster within their lifetime. Further, the cluster orbits are not consistent with being circular; the average value of r$_{apo}$ / r$_{peri}$ is $\sim$1.9 (equivalent to eccentricity of $\sim$0.31) for both clusters. Lastly, we find that the clusters do not share a common orbit, challenging one proposed formation scenario in which the clusters formed from molecular clouds on the open stream orbit derived by Kruijssen et al. (2015). Meanwhile, our constraints on the birth location and velocity of the clusters offer mild support for a scenario in which the clusters formed via collisions between gas clouds on the x1 and x2 bar orbit families.

The Dark Energy Spectroscopic Instrument (DESI) survey will measure large-scale structures using quasars as direct tracers of dark matter in the redshift range 0.9<z<2.1 and using Ly-alpha forests in quasar spectra at z>2.1. We present several methods to select candidate quasars for DESI, using input photometric imaging in three optical bands (g, r, z) from the DESI Legacy Imaging Surveys and two infrared bands (W1, W2) from the Wide-field Infrared Explorer (WISE). These methods were extensively tested during the Survey Validation of DESI. In this paper, we report on the results obtained with the different methods and present the selection we optimized for the DESI main survey. The final quasar target selection is based on a Random Forest algorithm and selects quasars in the magnitude range 16.5<r<23. Visual selection of ultra-deep observations indicates that the main selection consists of 71% quasars, 16% galaxies, 6% stars and 7% inconclusive spectra. Using the spectra based on this selection, we build an automated quasar catalog that achieves a >99% purity for a nominal effective exposure time of ~1000s. With a 310 per sq. deg. target density, the main selection allows DESI to select more than 200 QSOs per sq. deg. (including 60 quasars with z>2.1), exceeding the project requirements by 20%. The redshift distribution of the selected quasars is in excellent agreement with quasar luminosity function predictions.

The Dark Energy Spectroscopic Instrument (DESI) will precisely constrain cosmic expansion and the growth of structure by collecting $\sim$40 million extra-galactic redshifts across $\sim$80\% of cosmic history and one third of the sky. The Emission Line Galaxy (ELG) sample, which will comprise about one-third of all DESI tracers, will be used to probe the Universe over the $0.6 < z < 1.6$ range, which includes the $1.1<z<1.6$ range, expected to provide the tightest constraints. We present the target selection of the DESI SV1 Survey Validation and Main Survey ELG samples, which relies on the Legacy Surveys imaging. The Main ELG selection consists of a $g$-band magnitude cut and a $(g-r)$ vs.\ $(r-z)$ color box, while the SV1 selection explores extensions of the Main selection boundaries. The Main ELG sample is composed of two disjoint subsamples, which have target densities of about 1940 deg$^{-2}$ and 460 deg$^{-2}$, respectively. We first characterize their photometric properties and density variations across the footprint. Then we analyze the DESI spectroscopic data obtained since December 2020 during the Survey Validation and the Main Survey up to December 2021. We establish a preliminary criterion to select reliable redshifts, based on the \oii~flux measurement, and assess its performance. Using that criterion, we are able to present the spectroscopic efficiency of the Main ELG selection, along with its redshift distribution. We thus demonstrate that the the main selection with higher target density sample should provide more than 400 deg$^{-2}$ reliable redshifts in both the $0.6<z<1.1$ and the $1.1<z<1.6$ ranges.

We describe the Milky Way Survey (MWS) that will be undertaken with the Dark Energy Spectroscopic Instrument (DESI) on the Mayall 4m Telescope at the Kitt Peak National Observatory. Over the next 5 years DESI MWS will observe approximately 7 million stars at Galactic latitudes |b|>20 deg, with an inclusive target selection scheme focused on the thick disk and stellar halo. MWS will also include several high-completeness samples of rare stellar types, including white dwarfs, low-mass stars within 100pc of the Sun, and horizontal branch stars. We summarize the potential of DESI to advance understanding of Galactic structure and stellar evolution. We introduce the final definitions of the main MWS target classes and estimate the number of stars in each class that will be observed. We describe our pipelines to derive radial velocities, atmospheric parameters and chemical abundances. We use ~500,000 spectra of unique stellar targets from the DESI Survey Validation program (SV) to demonstrate that our pipelines can measure radial velocities to approximately 1 km/s and [Fe/H] accurate to approximately 0.2 dex for typical stars in our main sample. We find stellar parameter distributions from 100 sq. deg. of SV observations with >90% completeness on our main sample are in good agreement with expectations from mock catalogues and previous surveys.

The Dark Energy Spectroscopic Instrument (DESI) is carrying out a 5-year survey that aims to measure the redshifts of tens of millions of galaxies and quasars, including 8 million luminous red galaxies (LRGs) in the redshift range of $0.4<z<{\sim}\,1.0$. Here we present the selection of the DESI LRG sample and assess its spectroscopic performance using data from Survey Validation (SV) and the first 2 months of the Main Survey. The DESI LRG sample, selected using $g$, $r$, $z$ and $W1$ photometry from the DESI Legacy Imaging Surveys, is highly robust against imaging systematics. The sample has a target density of 605 deg$^{-2}$ and a comoving number density of $5\times10^{-4}\ h^3\mathrm{Mpc}^{-3}$ in $0.4<z<0.8$; this is significantly higher density than previous LRG surveys (such as SDSS, BOSS and eBOSS) while also extends to much higher redshifts. After applying a bright star veto mask developed for the sample, $98.9\%$ of the observed LRG targets yield confident redshifts (with a catastrophic failure rate of $0.2\%$ in the confident redshifts), and only $0.5\%$ of the LRG targets are stellar contamination. The LRG redshift efficiency varies with source brightness and effective exposure time, and we present a simple model that accurately characterizes this dependence. In the appendices, we describe the extended LRG samples observed during SV.

A key component of the Dark Energy Spectroscopic Instrument (DESI) survey validation (SV) is a detailed visual inspection (VI) of the optical spectroscopic data to quantify key survey metrics. In this paper we present results from VI of the quasar survey using deep coadded SV spectra. We show that the majority (~70%) of the main-survey targets are spectroscopically confirmed as quasars, with ~16% galaxies, ~6% stars, and ~8% low-quality spectra lacking reliable features. A non-negligible fraction of the quasars are misidentified by the standard DESI spectroscopic pipeline but we show that the majority can be recovered using post-pipeline "afterburner" quasar-identification approaches. We combine these "afterburners" with our standard pipeline to create a modified pipeline to improve the overall quasar completeness. At the depth of the main DESI survey both pipelines achieve a good-redshift purity (reliable redshifts measured within 3000 km/s) of ~99%; however, the modified pipeline recovers ~94% of the visually inspected quasars, as compared to just ~86% from the standard pipeline. We demonstrate that both pipelines achieve an overall redshift precision and accuracy of ~100 km/s and ~70 km/s, respectively. We constructed composite spectra to investigate why some quasars are missed by the standard spectroscopic pipeline and find that they are more host-galaxy dominated and/or dust reddened than the standard-pipeline quasars. We also show example spectra to demonstrate the overall diversity of the DESI quasar sample and provide strong-lensing candidates where two targets contribute to a single DESI spectrum.

Commonly used wavefront sensors, the Shack Hartmann wavefront sensor and the pyramid wavefront sensor, for example, have large dynamic range or high sensitivity, trading one regime for the other. A new type of wavefront sensor is being developed and is currently undergoing testing at the University of Arizona's Center for Astronomical Adaptive Optics. This sensor builds on linear optical differentiation theory by using linear, spatially varying halfwave plates in an intermediate focal plane. These filters, along with the polarizing beam splitters, divide the beam into four pupil images, similar to those produced by the pyramid wavefront sensor. The wavefront is then reconstructed from the local wavefront slope information contained in these images. The ODWFS is ideally suited for wavefront sensing on extended objects because of its large dynamic range and because it operates in a pupil plane which allows for on chip resampling even for arbitrarily shaped sources. We have assembled the ODWFS on a testbed using 32 by 32 square 1000 actuator deformable mirror to introduce aberration into a simulated telescope beam. We are currently testing the system's spatial frequency response and are comparing the resulting data to numerical simulations. This paper presents the results of these initial experiments.

We present Atacama Large Millimeter/submillimeter Array observations at $\approx100$ GHz with $0.05$ arcsec (3 pc) resolution of the kiloparsec-scale jet seen in the nearby Seyfert galaxy NGC 1068, and we report the presence of parsec-scale blobs at the head of the jet. The combination of the detected radio flux ($\approx0.8$ mJy), spectral index ($\approx0.5$), and the blob size ($\approx10$ pc) suggests a strong magnetic field of $B\approx240\,\mu$G. Such a strong magnetic field most likely implies magnetic field amplification by streaming cosmic rays. The estimated cosmic-ray power by the jet may exceed the limit set by the star formation activity in this galaxy. This result suggests that even modest-power jets can increase the galactic cosmic-ray content while propagating through the galactic bulge.

The relative abundances of singly-deuterated methanol isotopologues, [CH$_{2}$DOH]/[CH$_{3}$OD], in star-forming regions deviate from the statistically expected ratio of 3. In Orion KL, the nearest high-mass star-forming region to Earth, the singly-deuterated methanol ratio is about 1, and the cause for this observation has been explored through theory for nearly three decades. We present high-angular resolution observations of Orion KL using the Atacama Large Millimeter/submillimeter Array to map small-scale changes in CH$_{3}$OD column density across the nebula, which provide a new avenue to examine the deuterium chemistry during star and planet formation. By considering how CH$_{3}$OD column densities vary with temperature, we find evidence of chemical processes that can significantly alter the observed column densities. The astronomical data are compared with existing theoretical work and support D-H exchange between CH$_{3}$OH and heavy water (i.e., HDO and D$_{2}$O) at methanol's hydroxyl site in the icy mantles of dust grains. The enhanced CH$_{3}$OD column densities are localized to the Hot Core-SW region, a pattern that may be linked to the coupled evolution of ice mantel chemistry and star formation in giant molecular clouds. This work provides new perspectives on deuterated methanol chemistry in Orion KL and informs considerations that may guide future theoretical, experimental, and observational work.

The Astro2020 decadal survey recommended an infrared, optical, ultra-violet (IR/O/UV) telescope with a $\sim$6~m inscribed diameter and equipped with a coronagraph instrument to directly image exoEarths in the habitable zone of their host star. A telescope of such size may need to be segmented to be folded and then carried by current launch vehicles. However, a segmented primary mirror introduces the potential for additional mid spatial frequency optical wavefront instabilities during the science operations that would degrade the coronagraph performance. A coronagraph instrument with a wavefront sensing and control (WS\&C) system can stabilize the wavefront with a picometer precision at high temporal frequencies ($>$1Hz). In this work, we study a realistic set of aberrations based on a finite element model of a slightly larger (8m circumscribed, 6.7m inscribed diameter) segmented telescope with its payload. We model an adaptive optics (AO) system numerically to compute the post-AO residuals. The residuals then feed an end-to-end model of a vortex coronagraph instrument. We report the long exposure contrast and discuss the overall benefits of the adaptive optics system in the flagship mission success.

We present the design and science goals of SPT-3G+, a new camera for the South Pole Telescope, which will consist of a dense array of 34100 kinetic inductance detectors measuring the cosmic microwave background (CMB) at 220 GHz, 285 GHz, and 345 GHz. The SPT-3G+ dataset will enable new constraints on the process of reionization, including measurements of the patchy kinematic Sunyaev-Zeldovich effect and improved constraints on the optical depth due to reionization. At the same time, it will serve as a pathfinder for the detection of Rayleigh scattering, which could allow future CMB surveys to constrain cosmological parameters better than from the primary CMB alone. In addition, the combined, multi-band SPT-3G and SPT-3G+ survey data will have several synergies that enhance the original SPT-3G survey, including: extending the redshift-reach of SZ cluster surveys to $z > 2$; understanding the relationship between magnetic fields and star formation in our Galaxy; improved characterization of the impact of dust on inflationary B-mode searches; and characterizing astrophysical transients at the boundary between mm and sub-mm wavelengths. Finally, the modular design of the SPT-3G+ camera allows it to serve as an on-sky demonstrator for new detector technologies employing microwave readout, such as the on-chip spectrometers that we expect to deploy during the SPT-3G+ survey. In this paper, we describe the science goals of the project and the key technology developments that enable its powerful yet compact design.

We report object 282P/(323137) 2003 BM80 is undergoing a sustained activity outburst, lasting over 15 months thus far. These findings stem in part from our NASA Partner Citizen Science project Active Asteroids (this http URL), which we introduce here. We acquired new observations of 282P via our observing campaign (Vatican Advanced Technology Telescope, Lowell Discovery Telescope, and the Gemini South telescope), confirming 282P was active on UT 2022 June 7, some 15 months after 2021 March images showed activity in the 2021/2022 epoch. We classify 282P as a member of the Quasi-Hilda Objects, a group of dynamically unstable objects found in an orbital region similar to, but distinct in their dynamical characteristics to, the Hilda asteroids (objects in 3:2 resonance with Jupiter). Our dynamical simulations show 282P has undergone at least five close encounters with Jupiter and one with Saturn over the last 180 years. 282P was most likely a Centaur or Jupiter Family Comet (JFC) 250 years ago. In 350 years, following some 15 strong Jovian interactions, 282P will most likely migrate to become a JFC or, less likely, a main-belt asteroid. These migrations highlight a dynamical pathway connecting Centaurs and JFC with Quasi-Hildas and, potentially, active asteroids. Synthesizing these results with our thermodynamical modeling and new activity observations, we find volatile sublimation is the primary activity mechanism. Observations of a quiescent 282P, which we anticipate will be possible in 2023, will help confirm our hypothesis by measuring a rotation period and ascertaining spectral type.

We report the discovery of substellar objects in the young star cluster IC 348 and the neighboring Barnard 5 dark cloud, both at the eastern end of the Perseus star-forming complex. The substellar candidates are selected using narrowband imaging, i.e., on and off photometric technique with a filter centered around the water absorption feature at 1.45 microns, a technique proven to be efficient in detecting water-bearing substellar objects. Our spectroscopic observations confirm three brown dwarfs in IC 348. In addition, the source WBIS 03492858+3258064, reported in this work, is the first confirmed brown dwarf discovered toward Barnard 5. Together with the young stellar population selected via near- and mid-infrared colors using the Two Micron All Sky Survey and the Wide-field Infrared Survey Explorer, we diagnose the relation between stellar versus substellar objects with the associated molecular clouds. Analyzed by Gaia EDR3 parallaxes and kinematics of the cloud members across the Perseus region, we propose the star formation scenario of the complex under influence of the nearby OB association.

We present results from an HI counterpart search using the HI Parkes All Sky Survey (HIPASS) for a sample of low surface brightness galaxies (LSBGs) and ultradiffuse galaxies (UDGs) identified from the Dark Energy Survey (DES). We aimed to establish the redshifts of the DES LSBGs to determine the UDG fraction and understand their properties. Out of 409 galaxies investigated, none were unambiguously detected in HI. Our study was significantly hampered by the high spectral rms of HIPASS and thus in this paper we do not make any strong conclusive claims but discuss the main trends and possible scenarios our results reflect. The overwhelming number of non-detections suggest that: (A) Either all the LSBGs in the groups, blue or red, have undergone environment aided pre-processing and are HI deficient or the majority of them are distant galaxies, beyond the HIPASS detection threshold. (B) The sample investigated is most likely dominated by galaxies with HI masses typical of dwarf galaxies. Had there been Milky Way (MW) size (R_e) galaxies in our sample, with proportionate HI content, they would have been detected, even with the limitations imposed by the HIPASS spectral quality. This leads us to infer that if some of the LSBGs have MW size optical diameters, their HI content is possibly in the dwarf range. More sensitive observations using the SKA precursors in future may resolve these questions.

We study the behavior of dust temperature and its infrared emission of FirstLight1 simulated galaxies at the redshift of 6 and 8, by using POLARIS2 as a Monte Carlo photon transport simulator. To calculate the dust temperature ($T_{dust}$) of the Interstellar medium (ISM) of galaxies, POLARIS requires three essential parameters as an input - (1) The physical characteristics of galaxies such as the spatial distribution of stars and dust, which are taken from FirstLight galaxies. (2) The intrinsic properties of dust grains that are derived from theDiscrete Dipole Approximation Code (DDSCAT) model. (3) The optical properties of star-particles are in the form of their spectral energy distributions (SEDs) which are extracted from the Binary Population and Spectral Synthesis (BPASS) model. Our simulations produced the 3D maps of the equilibrium dust temperature along with the sight-line infrared emission maps of galaxies. Our results show the importance of excess heating of dust by the Cosmic Microwave Background (CMB) radiations at high redshifts that results in increased Mid and Far infrared (M-FIR) dust emission. The different evaluations of dust temperature models relate diversely to the optical and intrinsic properties of galaxies

X-ray mirror fabrication for astronomy is challenging; this is due to the Wolter I optical geometry and the tight tolerances on roughness and form error to enable accurate and efficient X-ray reflection. The performance of an X-ray mirror, and ultimately that of the telescope, is linked to the processes and technologies used to create it. The goal of this chapter is to provider the reader with an overview of the different technologies and processes used to create the mirrors for X-ray telescopes. The objective is to present this diverse field in the framework of the manufacturing methodologies (subtractive, formative, fabricative & additive) and how these methodologies influence the telescope attributes (angular resolution and effective area). The emphasis is placed upon processes and technologies employed in recent X-ray space telescopes and those that are being actively investigated for future missions such as Athena and concepts such as Lynx. Speculative processes and technologies relating to Industry 4.0 are introduced to imagine how X-ray mirror fabrication may develop in the future.

The Galactic source G2.4$+$1.4 is an optical and radio nebula containing an extreme Wolf--Rayet star. At one time this source was regarded as a supernova remnant, because of its apparent non-thermal radio spectrum, although this was based on limited observations. Subsequent observations instead supported a flat, optically thin thermal radio spectrum for G2.4$+$1.4, and it was identified as a photoionized, mass-loss bubble, not a supernova remnant. Recently, however, it has been claimed that this source has a non-thermal integrated radio spectrum. I discuss the integrated radio flux densities available for G2.4$+$1.4 from a variety of surveys, and show that it has a flat spectrum at gigahertz frequencies (with a spectral index $\alpha$ of $0.02 \pm 0.08$, where flux density $S$ scales with frequency $\nu$ as $S \propto \nu^{-\alpha}$).

Recently, some fast radio burst (FRB) repeaters were reported to exhibit complex, diverse variations of Faraday rotation measures (RMs), which implies that they are surrounded by an inhomogeneous, dynamically evolving, magnetized environment. We systematically investigate some possible astrophysical processes that may cause RM variations of an FRB repeater. The processes include (1) a supernova remnant (SNR) with a fluctuating medium; (2) a binary system with stellar winds from a massive/giant star companion or stellar flares from a low-mass star companion; (3) a pair plasma medium from a neutron star (including pulsar winds, pulsar wind nebulae and magnetar flares); (4) outflows from a massive black hole. For the SNR scenario, a large relative RM variation during a few years requires that the SNR is young with a thin and local anisotropic shell, or the size of dens gas clouds in interstellar/circumstellar medium around the SNR is extremely small. If the RM variation is caused by the companion medium in a binary system, it is more likely from stellar winds of a massive/giant star companion. The RM variation contributed by stellar flares from a low-mass star is disfavored, because this scenario predicts an extremely large relative RM variation during a short period of time. The scenarios invoking a pair plasma from a neutron star can be ruled out due to their extremely low RM contributions. Outflows from a massive black hole could provide a large RM variation if the FRB source is in the vicinity of the black hole.

We focus on post-asymptotic giant branch (post-AGB) binaries and study the interaction between the different components of these complex systems. These components comprise the post-AGB primary, a main sequence secondary, a circumbinary disk, as well as a fast bipolar outflow (jet) launched by the companion. We obtained well-sampled time series of high resolution optical spectra over the last decade and these spectra provide the basis of our study. The jet is detected in absorption, at superior conjunction, when the line of sight towards the primary goes through the bipolar cone. Our spectral time series scan the jets during orbital motion. Our spatio-kinematic model is constrained by these dynamical spectra. We complement this with a radiative-transfer model in which the Balmer series are used to derive total mass-loss rates in the jets. The jets are found to be wide and display an angle-dependent density structure with a dense and slower outer region near the jet cone and a fast inner part along the jet symmetry axes. The deprojected outflow velocities confirm that the companions are main sequence companions. The total mass-loss rates are large (10^{-8} and 10^{-5}\,solar mass per year), from which we can infer that the mass-accretion rates onto the companion star must be high as well. The circumbinary disk is likely the main source for the accretion disk around the companion. All systems with full disks that start near the sublimation radius show jets, whereas for systems with evolved transition disks, this lowers to a detection rate of 50%. Objects without an infrared excess do not show jets. We conclude that jet creation in post-AGB binaries is a mainstream process. The interaction between the circumbinary disks and the central binary provide the needed accretion flow, but the presence of a circumbinary disk does not seem to be the only prerequisite to launch a jet.

[Abridged] The ionization of polycyclic aromatic hydrocarbons (PAHs), by ultraviolet (UV) photons from massive stars is expected to account for a large fraction of the heating of neutral gas in galaxies. Evaluation of this proposal, however, has been limited by our ability to directly compare observational diagnostics to the results of a molecular model describing PAH ionization. The objective of this article is to take advantage of the most recent values of molecular parameters derived from laboratory experiments and quantum chemical calculations on PAHs and provide a detailed comparison between modeled values and observational diagnostics for the PAH charge state and the heating efficiency for PAHs. Despite the use of a simple analytical model, we obtain a good agreement between model results and observational diagnostics over a wide range of radiation fields and physical conditions, in environments such as star-forming regions, galaxies, and protoplanetary disks. In addition, we found that the modeled photoelectric heating rates by PAHs are close to the observed cooling rates given by the gas emission. These results show that PAH ionization is the main source of neutral gas heating in these environments. The results of our photoelectric heating model by PAHs can thus be used to assess the contribution of UV radiative heating in galaxies (vs shocks, for instance). We provide the empirical formulas fitted to the model results, and the full python code itself, to calculate the heating rates and heating efficiencies for PAHs.

Recently, Kojima and co-authors have reported a record low oxygen abundance, 12+logO/H=6.90+/-0.03 in the low-mass star-forming galaxy HSC J1631+4426. This exceptionally low oxygen abundance was obtained by the direct method, using the [OIII]4363 emission line. However, using the strong-line method by Izotov et al. (2019b), these authors have derived a significantly higher metallicity 12+logO/H=7.175+/-0.005. To clarify the situation, we have obtained new observations of HSC J1631+4426 with the Large Binocular Telescope (LBT)/Multi-Object Dual Spectrograph (MODS). We have derived a higher oxygen abundance, 12+logO/H=7.14+/-0.03, using the direct method, a value similar to the oxygen abundance obtained by the strong-line method. Thus, HSC J1631+4426 has a metallicity close to that of the well known blue compact dwarf galaxy IZw18.

Acetaldehyde is one of the most common and abundant gaseous interstellar complex organic molecules, found in cold and hot regions of the molecular interstellar medium. Its presence in the gas-phase depends on the chemical formation and destruction routes, and its binding energy (BE) governs whether acetaldehyde remains frozen onto the interstellar dust grains or not. In this work, we report a combined study of the acetaldehyde BE obtained via laboratory TPD (Temperature Programmed Desorption) experiments and theoretical quantum chemical computations. BEs have been measured and computed as a pure acetaldehyde ice and as mixed with both polycrystalline and amorphous water ice. Both calculations and experiments found a BE distribution on amorphous solid water that covers the 4000--6000 K range, when a pre-exponential factor of $1.1\times 10^{18}s^{-1}$ is used for the interpretation of the experiments. We discuss in detail the importance of using a consistent couple of BE and pre-exponential factor values when comparing experiments and computations, as well as when introducing them in astrochemical models. Based on the comparison of the acetaldehyde BEs measured and computed in the present work with those of other species, we predict that acetaldehyde is less volatile than formaldehyde, but much more than water, methanol, ethanol, and formamide. We discuss the astrochemical implications of our findings and how recent astronomical high spatial resolution observations show a chemical differentiation involving acetaldehyde, which can easily explained as due to the different BEs of the observed molecules.

Recent advancements in low-cost astronomical equipment, including high-quality medium-aperture telescopes and low-noise CMOS detectors, have made the deployment of large optical telescope arrays both financially feasible and scientifically interesting. The Argus Optical Array is one such system, composed of 900 eight-inch telescopes, which is planned to cover the entire night sky in each exposure and is capable of being the deepest and fastest Northern Hemisphere sky survey. With this new class of telescope comes new challenges: determining optimal individual telescope pointings to achieve required sky coverage and overlaps for large numbers of telescopes, and realizing those pointings using either individual mounts, larger mounting structures containing telescope subarrays, or the full array on a single mount. In this paper, we describe a method for creating a pointing pattern, and an algorithm for rapidly evaluating sky coverage and overlaps given that pattern, and apply it to the Argus Array. Using this pattern, telescopes are placed into a hemispherical arrangement, which can be mounted as a single monolithic array or split into several smaller subarrays. These methods can be applied to other large arrays where sky packing is challenging and evenly spaced array subdivisions are necessary for mounting.

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ~ 1.5 compared to SOL. Our simulation IMFs agree well with the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF.

Recent observations indicate that magnetars commonly reside in merging compact binaries and at least part of fast radio bursts (FRBs) are sourced by magnetar activities. It is natural to speculate that a class of merging neutron star binaries may have FRB emitters. In this work, we study the observational aspects of these binaries - particularly those with FRB repeaters, which are promising multi-band and multi-messenger observation targets of radio telescopes and ground based gravitational wave detectors as the former telescopes can probe the systems at a much earlier stage in the inspiral than the latter. We show that observations of FRB repeaters in compact binaries have a significant advantage in pinning down the binary spin dynamics, constraining neutron star equation of state, probing FRB production mechanisms, and testing beyond standard physics. As a proof of principle, we investigate several mock observations of FRB pulses originating from pre-merger neutron star binaries, and we find that using the information of FRB arriving times alone, the intrinsic parameters of this system (including the stellar masses, spins, and quadrupole moments) can be measured with high precision, and the angular dependence of the FRB emission pattern can also be well reconstructed. The measurement of stellar masses (with an error of $\mathcal{O}(10^{-6}-10^{-5})$) and quadrupole moments (with an error of $\mathcal{O}(1\%-10\%)$) may be an unprecedented discriminator of nuclear equations of state in neutron stars. In addition, we find the multi-band and multi-messenger observations of this binary will be sensitive to alternative theories of gravity and beyond standard models, e.g., dynamical Chern-Simons gravity and axion field that is coupled to matter.

The pulsar search was started at the radio telescope LPA LPI at the frequency 111~MHz. The first results deals of a search for right ascension $0^h - 24^h$ and declinations $+21^{\circ} - +42^{\circ}$ are presented in paper. The data with sampling 100 ms and with 6 frequency channals was used. It were found 34 pulsars. Seventeen of them previously been observed at radio telescope LPA LPI, and ten known pulsars has not previously been observed. It were found 7 new pulsars.

In this paper we describe a new way to understand the formation of galaxies via the infall of baryonic matter (BM) and dark matter (DM) onto a pre-existing over density. Unlike BM, DM particles can fly through this area without being captured. In this case it is impossible to explain the existence of a dark halo around galaxies, which contains most of the mass of the galaxy. We propose a simple model for DM capture. If during the flight the mass of the galaxy has increased, then slow DM particles are captured by the galaxy, further increasing its mass, while faster particles slow down, transferring part of their energy to the galaxy. This model allows to estimate the minimum initial velocity of a particle required for a passage without capture through the center of the galaxy and derive an nonlinear equation describing the rate of galaxy mass increase. An analysis carried out using the ideas of catastrophes theory shows that for intensive capture of dark matter, an increase in the mass of galactic baryonic matter is necessary, exceeding a certain threshold value. It may be associated with the accretion of matter or the merger of galaxies. Additionally, the density of intergalactic DM must exceed some threshold value. Then the rate of increase in the mass of DM can be much higher than the one of baryonic matter. The capture sharply decreases after the DM density drops below the threshold value due to expansion.

Cubic gold-platinum free-falling test masses (TMs) constitute the mirrors of future LISA and LISA-like interferometers for low-frequency gravitational wave detection in space. High-energy particles of Galactic and solar origin charge the TMs and thus induce spurious electrostatic and magnetic forces that limit the sensitivity of these interferometers. Prelaunch Monte Carlo simulations of the TM charging were carried out for the LISA Pathfinder (LPF) mission, that was planned to test the LISA instrumentation. Measurements and simulations were compared during the mission operations. The measured net TM charging agreed with simulation estimates, while the charging noise was three to four times higher. We aim to bridge the gap between LPF TM charging noise simulations and observations. New Monte Carlo simulations of the LPF TM charging due to both Galactic and solar particles were carried out with the FLUKA/LEI toolkit. This allowed propagating low-energy electrons down to a few electronvolt. These improved FLUKA/LEI simulations agree with observations gathered during the mission operations within statistical and Monte Carlo errors. The charging noise induced by Galactic cosmic rays is about one thousand charges per second. This value increases to tens of thousands charges per second during solar energetic particle events. Similar results are expected for the LISA TM charging.

The spin evolution of young protostars, surrounded by an accretion disk, still poses problems for observations and theoretical models. In recent studies, the importance of the magnetic star-disk interaction for stellar spin evolution has been elaborated. The accretion disk in these studies, however, is only represented by a simplified model and important features are not considered. We combined the implicit hydrodynamic TAPIR disk code with a stellar spin evolution model. The influence of stellar magnetic fields on the disk dynamics, the radial position of the inner disk radius, as well as the influence of stellar rotation on the disk were calculated self-consistently. Within a defined parameter space, we can reproduce the majority of fast and slow rotating stars observed in young stellar clusters. Additionally, the back reaction of different stellar spin evolutionary tracks on the disk can be analyzed. Disks around fast rotating stars are located closer to the star. Consequently, the disk midplane temperature in the innermost disk region increases significantly compared to slow rotating stars. We can show the effects of stellar rotation on episodic accretion outbursts. The higher temperatures of disks around fast rotating stars result in more outbursts and a longer outbursting period over the disk lifetime. The combination of a long-term hydrodynamic disk and a stellar spin evolution model allows the inclusion of previously unconsidered effects such as the back-reaction of stellar rotation on the long-term disk evolution and the occurrence of accretion outbursts. However, a wider parameter range has to be studied to further investigate these effects. Additionally, a possible interaction between our model and a more realistic stellar evolution code (e.g., the MESA code) can improve our understanding of the stellar spin evolution and its effects on the pre-main sequence star.

We aim to prepare the machine-learning ground for the next generation of spectroscopic surveys, such as 4MOST and WEAVE. Our goal is to show that convolutional neural networks can predict accurate stellar labels from relevant spectral features in a physically meaningful way. We built a neural network and trained it on GIRAFFE spectra with associated stellar labels from the sixth internal Gaia-ESO data release. Our neural network predicts the atmospheric parameters Teff and log(g) as well as the chemical abundances [Mg/Fe], [Al/Fe], and [Fe/H] for 30115 stellar spectra. The scatter of predictions from eight slightly different network models shows a high internal precision of the network results: 24 K for Teff, 0.03 for log(g), 0.02 dex for [Mg/Fe], 0.03 dex for [Al/Fe], and 0.02 dex for [Fe/H]. The network gradients reveal that the network is inferring the labels in a physically meaningful way from spectral features. Validation with benchmark stars and several scientific applications confirm that our network predictions are accurate for individual stars and recover the properties of different stellar populations in the Milky Way galaxy. Such a study provides very good insights into the application of machine-learning for the spectral analysis of large-scale spectroscopic surveys, such as WEAVE and 4MIDABLE-LR and -HR (4MOST Milky Way disk and bulge low- and high-resolution). The community will have to put a substantial effort into building proactive training sets for machine-learning methods to minimize the possible systematics.

Originating from several sources (Big Bang, stars, cosmic rays) and being strongly depleted during stellar lifetime, the lithium element (Li) is of great interest as its chemical evolution in the Milky Way is not yet well understood. To help constrain stellar and galactic chemical evolution models, numerous and precise Li abundances are necessary for a large range of evolutionary stages, metallicities, and Galactic volume. In the age of industrial parametrization, spectroscopic surveys such as APOGEE, GALAH, RAVE, and LAMOST have used data-driven methods to rapidly and precisely infer stellar labels (atmospheric parameters and abundances). To prepare grounds for future spectroscopic surveys like 4MOST and WEAVE, we aim to apply machine-learning techniques for Li study/measurement. We train a Convolution Neural-Network (CNN) coupling Gaia-ESO Survey iDR6 stellar labels ($\mathrm{\textit{T}_{\text{eff}}}$, log(g), [Fe/H] and A(Li)) and GIRAFFE HR15N spectra, to infer the atmospheric parameters and Li abundances for $\sim$40000 stars. We show that the CNN properly learns the physics of the stellar labels, from relevant spectral features, over a large range of evolutionary stages and stellar parameters. The Li feature at 6707.8 A is successfully singled out by our CNN, among the thousands of lines in the GIRAFFE HR15N setup. Rare objects like Li-rich giants are found in our sample. Such performances are achieved thanks to a meticulously built high-quality and homogeneous training sample. The CNN approach is very well adapted for the next generations of spectroscopic surveys aiming at studying (among other elements) Li, such as the 4MIDABLE-LR/HR (4MOST Milky Way disk and bulge low- and high-resolution) surveys. In this context, the caveats of the machine-learning applications should be properly investigated along with realistic label uncertainties and upper limits for abundances.

Hi-5 is the L'-band (3.5-4.0 $\mu$m) high-contrast imager of Asgard, an instrument suite in preparation for the visitor focus of the VLTI. The system is optimized for high-contrast and high-sensitivity imaging within the diffraction limit of a single UT/AT telescope. It is designed as a double-Bracewell nulling instrument producing spectrally-dispersed (R=20, 400, or 2000) complementary nulling outputs and simultaneous photometric outputs for self-calibration purposes. In this paper, we present an update of the project with a particular focus on the overall architecture, opto-mechanical design of the warm and cold optics, injection system, and development of the photonic beam combiner. The key science projects are to survey (i) nearby young planetary systems near the snow line, where most giant planets are expected to be formed, and (ii) nearby main sequence stars near the habitable zone where exozodiacal dust that may hinder the detection of Earth-like planets. We present an update of the expected instrumental performance based on full end-to-end simulations using the new GRAVITY+ specifications of the VLTI and the latest planet formation models.

We investigate how the constraints on cosmological and astrophysical parameters ($\Omega_{\rm m}$, $\sigma_{8}$, $A_{\rm SN1}$, $A_{\rm SN2}$) vary when exploiting information from multiple fields in cosmology. We make use of a convolutional neural network to retrieve the salient features from different combinations of field maps from IllustrisTNG in the CAMELS project. The fields considered are neutral hydrogen (HI), gas density (Mgas), magnetic fields (B) and gas metallicity (Z). We estimate the predictive uncertainty on the predictions of our model by using Monte Carlo dropout, a Bayesian approximation. Results show that overall, the performance of the model improves on all parameters as the number of channels of its input is increased. As compared to previous works, our model is able to predict the astrophysical parameters with up to $5\%$ higher in accuracy. In the best setup which includes all fields (four channel input, Mgas-HI-B-Z) the model achieves $R^{2} > 0.96$ on all parameters. Similarly, we find that the total uncertainty, which is dominated by the aleatoric uncertainty, decreases as more fields are used to train the model in general. The uncertainties obtained by dropout variational inference are overestimated on all parameters in our case, in that the predictive uncertainty is much larger than the actual squared error. After calibration, which consists of a simple $\sigma$ scaling method, the average deviation of the total uncertainty from the actual error goes down to $25\%$ at most (on $A_{\rm SN1}$).

Hercules X-1 is a nearly edge-on accreting X-ray pulsar with a warped accretion disk, precessing with a period of about 35 days. The disk precession allows for unique and changing sightlines towards the X-ray source. To investigate the accretion flow at a variety of sightlines, we obtained a large observational campaign on Her X-1 with XMM-Newton (380 ks exposure) and Chandra (50 ks exposure) for a significant fraction of a single disk precession cycle, resulting in one of the best datasets taken to date on a neutron star X-ray binary. Here we present the spectral analysis of the High State high-resolution grating and CCD datasets, including the extensive archival data available for this famous system. The observations reveal a complex Fe K region structure, with three emission line components of different velocity widths. Similarly, the high-resolution soft X-ray spectra reveal a number of emission lines of various widths. We correct for the uncertain gain of the EPIC-pn Timing mode spectra, and track the evolution of these spectral components with Her X-1 precession phase and observed luminosity. We find evidence for three groups of emission lines: one originates in the outer accretion disk (10^5 RG from the neutron star). The second line group plausibly originates at the boundary between the inner disk and the pulsar magnetosphere (10^3 RG). The last group is too broad to arise in the magnetically-truncated disk and instead must originate very close to the neutron star surface, likely from X-ray reflection from the accretion curtain (~10^2 RG).

The relationship between solar eruption and sunspot rotation has been widely reported, and the underlying mechanism requires to be studied. Here we performed a full 3D MHD simulation of data-constrained approach to study the mechanism of flare eruptions in active region (AR) NOAA 10930, which is characterized by continuous sunspot rotation and homologous eruptions. We reconstructed the potential magnetic field from the magnetogram of Hinode/SOT as the initial condition and drove the MHD system by applying continuous sunspot rotation at the bottom boundary. The key magnetic structure before the major eruptions and the pre-formed current sheet were derived, which is responsible for the complex MHD evolution with multiple stages. The major eruptions were triggered directly by fast reconnection in the pre-formed current sheet above the main polarity inversion line between the two major magnetic polarities of the AR. Furthermore, our simulation shows the homologous eruption successfully. It has reasonable consistence with observations in relative strength, energy release, X-ray and H{\alpha} features and time interval of eruptions. In addition, the rotation angle of the sunspot before the first eruption in the simulation is also close to the observed value. Our simulation offers a scenario different from many previous studies based on ideal instabilities of twisted magnetic flux rope, and shows the importance of sunspot rotation and magnetic reconnection in efficiently producing homologous eruptions by continuous energy injection and impulsive energy release in a recurrent way.

Ultraluminous X-ray bursts (hereafter ULXBs) are ultraluminous X-ray flares with a fast rise ($\sim$ one minute) and a slow decay ($\sim$ an hour), which are commonly observed in extragalactic globular clusters. Most ULXBs are observational one-off bursts, whereas five flares from the same source in NGC 5128 were discovered by Irwin et al. (2016). In this Letter, we propose a neutron star (NS)-white dwarf (WD) binary model with super-Eddington accretion rates to explain the repeating behavior of the ULXB source in NGC 5128. With an eccentric orbit, the mass transfer occurs at the periastron where the WD fills its Roche lobe. The ultraluminous X-ray flares can be produced by the accretion column around the NS magnetic poles. On the other hand, some repeating fast radio bursts (hereafter FRBs) were also found in extragalactic globular clusters. Repeating ULXBs and repeating FRBs are the most violent bursts in the X-ray and radio bands, respectively. We propose a possible association between the repeating ULXBs and the repeating FRBs. Such an association is worth further investigation by follow-up observations on nearby extragalactic globular clusters.

We present a novel mechanism for thermal dark matter production, characterized by a "bounce": the dark matter equilibrium distribution transitions from the canonical exponentially falling abundance to an exponentially rising one, resulting in an enhancement of the freezeout abundance by several orders of magnitude. We discuss multiple qualitatively different realizations of bouncing dark matter. The bounce allows the present day dark matter annihilation cross section to be significantly larger than the canonical thermal target, improving the prospects for indirect detection signals.

The worldine effective field theory (EFT) gives a gauge-invariant definition of black hole conservative tidal responses (Love numbers), dissipation numbers, and their spin-0 and spin-1 analogs. In the first part of this paper we show how the EFT allows us to circumvent the source/response ambiguity without having to use the analytic continuation prescription. The source/response ambiguity appears if Post-Newtonian (PN)corrections to external sources overlap with the response. However, these PN corrections can be clearly identified and isolated using the EFT.We illustrate that by computing static one-point functions of various external fields perturbing the four-dimensional Schwarzschild geometry. Upon resumming all relevant Feynman diagrams, we find that the PN terms that may mimic the response actually vanish for static black holes. Thus, the extraction of Love numbers from matching the EFT and general relativity (GR) calculations is completely unambiguous, and it implies that the Love numbers vanish identically for all types of perturbations. We also study in detail another type of fine tuning in the EFT, the absence of Love numbers' running. We show that logarithmic corrections to Love numbers do stem from individual loop diagrams in generic gauges, but cancel after all diagrams are summed over.In the particular cases of spin-0 and spin-2 fields the logarithms are completely absent if one uses the Kaluza-Klein metric decomposition.In the second part of the paper we compute frequency-dependent dissipative response contributions to the one-point functions using the Schwinger-Keldysh formalism. We extract black hole dissipation numbers by comparing the one-point functions in the EFT and GR. Our results are in perfect agreement with those obtained from a manifestly gauge-invariant matching of absorption cross-sections.

In the standard picture of stellar evolution, pair-instability -- the energy loss in stellar cores due to electron-positron pair production -- is predicted to prevent the collapse of massive stars into black holes with mass in the range between approximately 50 and 130 solar masses -- a range known as the "{\em black hole mass gap}". LIGO detection of black hole binary mergers containing one or both black holes with masses in this {\em mass gap} thus challenges the standard picture, possibly pointing to an unexpected merger history, unanticipated or poorly understood astrophysical mechanisms, or new physics. Here, we entertain the possibility that a "dark sector" exists, consisting of dark electrons, dark protons, and electromagnetic-like interactions, but no nuclear forces. Dark stars would inevitably form given such dark sector constituents, possibly collapsing into black holes with masses within the mass gap. We study in detail the cooling processes necessary for successful stellar collapse in the dark sector and show that for suitable choices of the particle masses, we indeed predict populating the mass gap with dark sector black holes. In particular, we numerically find that the heavier of the two dark sector massive particles cannot be lighter than, approximately, the visible sector proton for the resulting dark sector black holes to have masses within the mass gap. We discuss constraints on this scenario and how to test it with future, larger black hole merger statistics.

Cosmic rays and air showers at ultra-high energy are unique tools to test the validity of Lorentz invariance. A brief overview is given on such tests focusing on isotropic, non-birefringent Lorentz violation (LV) in the photon sector. Based on the apparent absence of vacuum Cherenkov radiation and photon decay, the LV parameter $\kappa$ is bound to $-0.6 \cdot 10^{-20} < \kappa < 6 \cdot 10^{-20}$ (98\% CL). We report an updated limit from cosmic-ray photon observations and preliminary results on testing vacuum Cherenkov radiation in air showers.

The use of Deep Learning (DL) algorithms has improved the performance of vision-based space applications in recent years. However, generating large amounts of annotated data for training these DL algorithms has proven challenging. While synthetically generated images can be used, the DL models trained on synthetic data are often susceptible to performance degradation, when tested in real-world environments. In this context, the Interdisciplinary Center of Security, Reliability and Trust (SnT) at the University of Luxembourg has developed the 'SnT Zero-G Lab', for training and validating vision-based space algorithms in conditions emulating real-world space environments. An important aspect of the SnT Zero-G Lab development was the equipment selection. From the lessons learned during the lab development, this article presents a systematic approach combining market survey and experimental analyses for equipment selection. In particular, the article focus on the image acquisition equipment in a space lab: background materials, cameras and illumination lamps. The results from the experiment analyses show that the market survey complimented by experimental analyses is required for effective equipment selection in a space lab development project.