Papers for Tuesday, Oct 19 2021

Compact star-forming clumps observed in distant galaxies are often suggested to play a crucial role in galaxy assembly. In this paper, we use a novel approach of applying finite resolution deconvolution on ground-based images of the COSMOS field to resolve 20,185 star-forming galaxies (SFG) at 0.5<z<2 to an angular resolution of 0.3", and study their clumpy fractions. A comparison between the deconvolved and HST images across four different filters shows good agreement and validates the deconvolution. We model spectral energy distributions using the deconvolved 14-band images to provide resolved surface brightness and stellar mass density maps for these galaxies. We find that the fraction of clumpy galaxies decreases with increasing stellar masses, and with increasing redshift: from ~30% at z ~ 0.7 to ~50% at z ~ 1.7. Using abundance matching, we also trace the progenitors for galaxies at z ~ 0.7 and measure the fractional mass contribution of clumps toward their total mass budget. Clumps are observed to have a higher fractional mass contribution toward galaxies at higher redshift: increasing from ~1% at z ~ 0.7 to ~5% at z ~ 1.7. Finally, the majority of clumpy SFGs have higher specific star formation rates (sSFR) compared to the average SFGs at fixed stellar mass. We discuss the implication of this result to in-situ clump formation due to disk instability.

The central dark-matter fraction of galaxies is sensitive to feedback processes during galaxy formation. Strong gravitational lensing has been effective in the precise measurement of the dark-matter fraction inside massive early-type galaxies. Here, we compare the projected dark-matter fraction of early-type galaxies inferred from the SLACS strong-lens survey, with those obtained from the EAGLE, Illustris, and IllustrisTNG hydro-dynamical simulations. Previous comparisons with some simulations revealed a large discrepancy, with considerably higher inferred dark-matter fractions -- by factors 2-3 -- inside half of the effective radius in observed strong-lens galaxies as compared to simulated galaxies. Here, we report good agreement between EAGLE and SLACS for the dark-matter fractions inside both half of the effective radius and the effective radius as a function of the galaxy's stellar mass, effective radius, and total mass-density slope. However, for IllustrisTNG and Illustris, the dark-matter fractions are lower than observed. This work consistently assumes a Chabrier IMF, which suggests that a different IMF (although not excluded) is not necessary to resolve this mismatch. The differences in the stellar feedback model between EAGLE and Illustris and IllustrisTNG, are likely the dominant cause of the difference in their dark-matter fraction, and density slope.

The values of, and connection between, the cosmological observables of the primordial power spectrum tilt $n_s$ and the inflationary tensor to scalar ratio $r$ are key guideposts to the physics of inflation. Universality classes can be defined for the tilt from the scale free value proportional to $1/N$, where $N$ is the number of e-folds. We examine the consequences of a $\ln N$ next to leading order correction rather than an expansion in $1/N$, or introducing a new parameter. While nominally this can lower $r$ for some too-high $r$ simple inflation models (e.g. large field models), there is an interesting cancellation preventing such models from coming back into favor. On the other branch of the universality class, near Starobinsky inflation, $r$ can be raised, making it easier to detect.

We present a new, multi-mission catalogue of ultraluminous X-ray source (ULX) candidates, based on recent data releases from each of the XMM-Newton, Swift and Chandra observatories (the 4XMM-DR10, 2SXPS and CSC2 catalogues, respectively). This has been compiled by cross-correlating each of these X-ray archives with a large sample of galaxies primarily drawn from the HyperLEDA archive. Significant efforts have been made to clean the sample of known non-ULX contaminants (e.g. foreground stars, background AGN, supernovae), and also to identify ULX candidates that are common to the different X-ray catalogues utilized, allowing us to produce a combined 'master' list of unique sources. Our sample contains 1843 ULX candidates associated with 951 different host galaxies, making it the largest ULX catalogue compiled to date. Of these, 689 sources are catalogued as ULX candidates for the first time. Our primary motivation is to identify new sources of interest for detailed follow-up studies, and within our catalogue we have already found one new extreme ULX candidate that has high S/N data in the archive: NGC 3044 ULX1. This source has a peak luminosity of $L_{\rm{X,peak}} \sim 10^{40}$ erg s$^{-1}$, and the XMM-Newton spectrum of the source while at this peak flux is very similar to other, better-studied extreme ULXs that are now understood to be local examples of super-Eddington accretion. This likely indicates that NGC 3044 ULX1 is another source accreting at super-Eddington rates. We expect that this catalogue will be a valuable resource for planning future observations of ULXs - both with our current and future X-ray facilities - to further improve our understanding of this enigmatic population.

We obtained spectra of ASASSN-V J205543.90+240033.5 (J2055), a system that shows photometric variations similar to the white dwarf (WD) pulsar AR Scorpii (Kato et al. arXiv:2109.03979). Our spectra display a continuum rising steeply toward the blue as well as an array of emission lines. Resolved Balmer and Paschen lines are seen with H$\alpha$ and H$\beta$ having central absorption features. The strongest lines are unresolved CII, CIII, and NIII as well as doubly ionized helium. The spectra are similar to that of YY Hya (Kimeswenger et al. arXiv:2110.03935), and suggest that J2055 is a post-common envelope binary consisting of a hot compact star irradiating the face of a secondary of unknown spectral type. Velocity variations detected from the emission lines confirm the binary nature of J2055. The origin of the 10 minute photometric variation remains uncertain.

Standard sirens (SS) are the gravitational wave analog of the astronomical standard candles, and can provide powerful information about the dynamics of the Universe up to very high $z$ values. In this work, we generate three mock SS catalogs based on the merger of massive black hole binaries which are expected to be observed in the LISA operating frequency band. Then, we perform an analysis to test modifications of general relativity (GR) inspired by the No Slip Gravity framework. We find that in the best scenarios, we can constrain the free parameters which quantify deviations from GR to 21\% accuracy, while the Hubble parameter can be simultaneously fit to 6\% accuracy. In combination with CMB information, we find a 15\% accuracy on the modified gravity free parameters and 0.7\% accuracy on the Hubble parameter. Therefore, the SS events at very large cosmological distances to be observed in LISA band will provide a unique way to test nature of gravity.

The microlensing effect has developed into a powerful technique for a diverse range of applications including exoplanet discoveries, structure of the Milky Way, constraints on MAssive Compact Halo Objects, and measurements of the size and profile of quasar accretion discs. In this paper, we consider a special type of microlensing events where the sources are fast radio bursts with $\sim$milliseconds (ms) durations for which the relative motion between the lens and source is negligible. In this scenario, it is possible to temporally resolve the individual microimages. As a result, a method beyond the inverse ray shooting (IRS) method, which only evaluates the total magnification of all microimages, is needed. We therefore implement an algorithm for identifying individual microimages and computing their magnifications and relative time delays. We validate our algorithm by comparing to analytical predictions for a single microlens case and find excellent agreement. We show that the superposition of pulses from individual microimages produces a light curve that appears as multi-peaked FRBs. The relative time delays between pulses can reach 0.1--1 ms for stellar-mass lenses and hence can already be resolved temporally by current facilities. Although not yet discovered, microlensing of FRBs will become regular events and surpass the number of quasar microlensing events in the near future when $10^{4-5}$ FRBs are expected to be discovered on a daily basis. Our algorithm provides a way of generating the microlensing light curve that can be used for constraining stellar mass distribution in distant galaxies.

Quiescent filaments are usually affected by internal and/or external perturbations triggering oscillations of different kinds. In particular, external large-scale coronal waves can perturb remote quiescent filaments leading to large-amplitude oscillations. Observational reports have indicated that the activation time of oscillations coincides with the passage of a large-scale coronal wavefront through the filament, although the disturbing wave is not always easily detected. Aiming to contribute to understand how -- and to what extent -- coronal waves are able to excite filament oscillations, here we modelled with 2.5 MHD simulations a filament floating in a gravitationally stratified corona disturbed by a coronal shock wave. This simplified scenario results in a two-coupled oscillation pattern of the filament which is damped in a few cycles, enabling a detailed analysis. A parametric study was accomplished varying parameters of the scenario such as height, size and mass of the filament. An oscillatory analysis reveals a general tendency where periods of oscillations, amplitudes and damping times increase with height, whereas filaments of larger radius exhibit shorter periods and smaller amplitudes. The calculation of forces exerted on the filament shows that the main restoring force is the magnetic tension.

The REgolith X-ray Imaging Spectrometer (REXIS) instrument on board NASA's OSIRIS-REx mission to the asteroid Bennu is a Class-D student collaboration experiment designed to detect fluoresced X-rays from the asteroid's surface to measure elemental abundances. In July and November 2019 REXIS collected ~615 hours of integrated exposure time of Bennu's sun-illuminated surface from terminator orbits. As reported in Hoak et al. (2021), the REXIS data do not contain a clear signal of X-ray fluorescence from the asteroid, in part due to the low incident solar X-ray flux during periods of observation. To support the evaluation of the upper limits on the detectable X-ray signal that may provide insights for the properties of Bennu's regolith, we present an overview of the REXIS instrument, its operation, and details of its in-flight calibration on astrophysical X-ray sources. This calibration includes the serendipitous detection of the transient X-ray binary MAXI J0637-430 during Bennu observations, demonstrating the operational success of REXIS at the asteroid. We convey some lessons learned for future X-ray spectroscopy imaging investigations of asteroid surfaces.

We present the light curve for J1415+1320 blazar monitored at 327 MHz using the Ooty Radio Telescope from the period 1989 to 2018. The source displayed variability more than 5 $\sigma$ value above its mean flux density (2.70$\pm$0.02 Jy) at two epochs around 2008 and 2009. Here $\sigma$ is the RMS variation of the flux density over the epoch. The variability analysis techniques indicate weak variability but inadequate to quantify if it is due to an intrinsic or extrinsic mechanism.

Superclusters are the largest objects in the Universe, and they provide a unique opportunity to study how galaxy clusters are born at the junction of the cosmic web as well as the distribution of magnetic fields and relativistic particles beyond cluster volume. The field of radio astronomy is going through an exciting and important era of the Square Kilometer Array (SKA). We now have the most sensitive functional radio telescopes, such as the MeerKAT, which offers high angular resolution and sensitivity towards diffuse and faint radio sources. To study the radio environments around supercluster, we observed the (core part of) {\it Saraswati} supercluster with the MeerKAT. From our MeerKAT Observation of the {\it Saraswati} Supercluster (MOSS) project, the initial results of the pilot observations of two massive galaxy clusters, A2631 and ZwCl2341.1+0000, which are located around the dense central part of the {\it Saraswati} supercluster, were discussed. In this paper, we describe the observations and data analysis details, including direction-dependent calibration. In particular, we focus on the ZwCl2341.1+0000 galaxy cluster, which hosts double radio relics and puzzling diffuse radio source in the filamentary network. We have imaged these double radio relics in our high resolution and sensitive L-band MeerKAT observation and a puzzling radio source, located between relics, in the low-resolution image. We also derived the spectra of double radio relics using MeerKAT and archival GMRT observations. A following papers will focus on the formation of radio relics and halo, as well as radio galaxy properties in a supercluster core environment.

Low-metallicity, compact starburst galaxies referred to as Green Peas (GPs) provide a unique window to study galactic evolution across cosmic epochs. In this work, we present new deep optical spectra for three GPs from OSIRIS at the 10m Gran Telescopio Canarias (GTC), which are studied using a state-of-the-art methodology. A stellar population synthesis is conducted with 1098 spectral templates. The methodology succeeds at characterising stellar populations from 0.5 Myrs to 10 Gyrs. The light distribution shows a large red excess from a single population with $log\left(age\right) > 8.5yr$ in the GP sample analysed. This points towards an incomplete characterisation of the gas luminosity, whose continuum already accounts between $7.4\%$ and $27.6\%$ in the galaxy sample. The emission spectra are fitted with the largest Bayesian chemical model consisting of a electron temperature, a electron density, the logarithmic extinction coefficient and eleven ionic species under the direct method paradigm. Additionally, building on previous work, we propose a neural networks sampler to constrain the effective temperature and ionization parameter of each source from photoionization model grids. Finally, we combine both methodologies into a 16-dimensional model, which for the first time, simultaneously explores the direct method and photoionization parameter spaces. Both techniques consistently indicate a low metallicity gas, $7.76<12+log\left(\frac{O}{H}\right)<8.04$, ionized by strong radiation fields, in agreement with previous works.

We present a new approach to capturing the broad diversity of emission line and continuum properties in quasar spectra. We identify populations of spectrally similar quasars through pixel-level clustering on 12,968 high signal-to-noise spectra from the Sloan Digital Sky Survey (SDSS) in the redshift range of $1.57<z<2.4$. Our clustering analysis finds 396 quasar spectra that are not assigned to any population, 15 misclassified spectra, and six quasars with incorrect redshifts. We compress the quasar populations into a library of 684 high signal-to-noise composite spectra, anchored in redshift space by the Mg II emission line. Principal component analysis on the library results in an eigenspectrum basis spanning 1067 to 4007 $\r{A}$. We model independent samples of SDSS quasar spectra with the eigenbasis, allowing for a free redshift parameter. Our models achieve a median reduced chi-squared on non-BAL quasar spectra that is reduced by 8.5% relative to models using the eigenspectra from the SDSS spectroscopic pipeline. A significant contribution to the relative improvement is from the ability to reconstruct the range of emission line variation. The redshift estimates from our model are consistent with the Mg II emission line redshift with an average offset that displays 51.4% less redshift-dependent variation relative to the SDSS eigenspectra. Our method for developing quasar spectra models can improve automated classification and predict the intrinsic spectrum in regions affected by intervening absorbers such as Ly$\alpha$, C IV, and Mg II, thus benefiting studies of large-scale structure.

Three structural isomers of the C2H4O2 molecule, namely, methyl formate (MF; HCOOCH3), acetic acid (AA; CH3COOH), and glycol aldehyde (GA; HOCH2CHO), have attracted considerable attention as targets for understanding pathways towards molecular complexity in the interstellar medium (ISM). Among these isomers, MF is decisively abundant in various astronomical objects. For various formation pathways of MF, surface reactions on cosmic dust would play an important role. However, when compared to observations, the formation of MF has been found to be relatively inefficient in laboratory experiments in which methanol (CH3OH)-dominant ices were processed by ultraviolet (UV) photons and cosmic-ray analogues. Here, we show experimental results on the effective formation of MF by the photolysis of CH3OH on water ice at 10 K. We found that the key parameter leading to the efficient formation of MF is the supply of OH radicals by the photolysis of H2O, which significantly differs from CH3OH-rich experimental conditions. Moreover, using an ultra-high-sensitivity surface analysis method, we succeeded in constraining the decisive formation pathway of MF via the photolysis of methoxymethanol (MM; CH3OCH2OH), which would improve our current understanding of chemical evolution in the ISM.

Instabilities described by linear theory characterize an important form of wave-particle interaction in the solar wind. We diagnose unstable behavior of solar wind plasma between 0.3 and 1 au via the Nyquist criterion, applying it to fits of $\sim1.5$M proton and $\alpha$ particle Velocity Distribution Functions (VDFs) observed by \emph{Helios I} and \emph{II}. The variation of the fraction of unstable intervals with radial distance from the Sun is linear, signaling a gradual decline in the activity of unstable modes. When calculated as functions of the solar wind velocity and Coulomb number, we obtain more extreme, exponential trends in the regions where collisions appear to have a notable influence on the VDF. Instability growth rates demonstrate similar behavior, and significantly decrease with Coulomb number. We find that, for a non-negligible fraction of observations, the proton beam or secondary component might not be detected due to instrument resolution limitations, and demonstrate that the impact of this issue does not affect the main conclusions of this work.

Owing to lack of multiple components of prompt $\gamma$-ray emissions in short gamma-ray bursts (sGRBs), how these distinct components are correlated still keeps unclear. In this paper, we investigate the spectral and temporal properties of precursors, main peaks and extended emissions in 26 sGRBs including GRB 170817A. It is found that peak energies ($E_p$) in each pulse are uncorrelated with the pulse duration ($t_{dur}$). Meanwhile, we find that there is no obvious correlation between peak energy and energy fluence. Interestingly, there is no obvious spectral evolution from earlier precursors to later extended emissions in view of the correlations of $t_{dur}$ with either the $E_p$ or the low energy spectrum index, $\alpha$. A power-law correlation between the average flux ($F_{p}$) and the energy fluence ($S_\gamma$), $log F_p=(0.62\pm0.07) log S_\gamma + (0.27\pm0.07)$, is found to exist in the individual segments instead of mean peaks previously. Furthermore, we also find that the main peaks are on average brighter than the precursors or the extend emissions about one order of magnitude. On the basis of all the above analyses, one can conclude that three emissive components would share the same radiation mechanisms but they might be dominated by diverse physical processes.

Interplanetary scintillation (IPS) provides an approach for identifying the presence of sub-arcsec structures in radio sources, and very long baseline interferometry (VLBI) technique can help verify whether the IPS sources have fine structures on milli-arcsec (mas) scales. We searched the available VLBI archive for the 244 IPS sources detected by the Murchison Widefield Array at 162~MHz and found 63 cross-matches. We analysed the VLBI data of the 63 sources and characterised the compactness index in terms of the ratio $R$ of the VLBI-measured flux density at 4.3~GHz to the flux density estimated using the Very Large Array Sky Survey (VLASS) at 3~GHz and NRAO VLA Sky Survey (NVSS) at 1.4~GHz ($S_{\rm VLBI}/S_{\rm SA}$). Eleven sources are identified as blazars according to their flat spectra and strong variability. They show high compactness indices with $R>0.4$, compact core-jet structure, and a broad distribution of normalised scintillation index (NSI). Other sources show diverse morphologies (compact core, core and one-sided jet, core and two-sided jets), but there is a correlation between $R$ and NSI with a correlation coefficient $r=0.47$. A similar $R$--NSI correlation is found in sources showing single steep power-law or convex spectra. After excluding blazars (which are already known to be compact sources) from the VLBI-detected IPS sources, a strong correlation is found between the compactness and scintillation index of the remaining samples, indicating that stronger scintillating sources are more compact. This pilot study shows that IPS offers a convenient method to identify compact radio sources without the need to invoke high-resolution imaging observations, which often require significant observational time.

This work focuses on the study of the reconfiguration strategies available for uniformly distributed satellite constellations and slotting architectures. Particularly, this manuscript deals with the cases of reducing, maintaining, and also increasing the number of available positions for satellites in the space structure, and takes into account the potential minimum distances between spacecraft in the configuration to assure the safety of the system. To that end, several approaches to solve the reconfiguration problem are presented based on the properties of Flower Constellations, and more particularly, on the properties of uniformity and symmetries present in these uniform distributions.

This work focuses on the generation of non-self-intersecting relative trajectories, and their applications to satellite constellation design, slotting architectures, and Space Traffic Management. To that end, this paper introduces two theorems to determine when two spacecrafts share the same relative trajectory, and to identify the only conditions that allow the existence of non-self-intersecting relative trajectories. Then, these results are applied first to the estimation of the limits of the orbital capacity at a given altitude, and second, to the design of satellite constellations and slotting architectures that present no conjunctions between any element compliant with these space structures.

It has been proposed that the filament environment is closely connected to the pre-processing of galaxies, where their properties may have been changed by environmental effects in the filament before they fell into the galaxy cluster. We present the chemical properties of star-forming dwarf galaxies (SFDGs) in five filamentary structures (Virgo III, Leo Minor, Leo II A, Leo II B, and Canes Venatici) around the Virgo cluster using the Sloan Digital Sky Survey optical spectroscopic data and Galaxy Evolution Explorer ultraviolet photometric data. We investigate the relationship between stellar mass, gas-phase metallicity, and specific star formation rate (sSFR) of SFDGs in the Virgo filaments in comparison to those in the Virgo cluster and field. We find that, at a given stellar mass, SFDGs in the Virgo filaments show lower metallicity and higher sSFR than those in the Virgo cluster on average. We observe that SFDGs in the Virgo III filament show enhanced metallicities and suppressed star formation activities comparable to those in the Virgo cluster, whereas SFDGs in the other four filaments exhibit similar properties to the field counterparts. Moreover, about half of the galaxies in the Virgo III filament are found to be morphologically transitional dwarf galaxies that are supposed to be on the way to transforming into quiescent dwarf early-type galaxies. Based on the analysis of the galaxy perturbation parameter, we propose that the local environment represented by the galaxy interactions might be responsible for the contrasting features in "chemical pre-processing" found in the Virgo filaments.

Plasma loops are the elementary structures of solar flaring active regions which dominate the whole processes of flaring eruptions. The standard flare models are well explained the evolution and eruption after magnetic reconnection around the hot cusp-structure above the top of plasma loops, however, the early evolution of the plasma loops before the onset of magnetic reconnection has been poorly understood. Considering that magnetic-gradients are ubiquitous in solar plasma loops, this work applies the magnetic-gradient pumping (MGP) mechanism to study the early evolution of flaring plasma loops. The results indicate that the early evolution depend on the magnetic field distribution and the geometry of the plasma loops which dominates the balance between the accumulation and dissipation of energy around loop-tops. Driven by MGP process, both of the density and temperature as well as the plasma beta value around the looptop will increase in the early phase of the plasma loops evolution. In fact, the solar plasma loops will have two distinct evolutionary results: the low, initial dense plasma loops with relatively strong magnetic fields tend to be stable for their maximum beta value always smaller than the critical value, while the higher, initial dilute solar plasma loops with relatively weak magnetic fields tend to be unstable for their beta values exceeding the critical value at a time of about one hour after the formation of the solar magnetized plasma loop. The latter may produce ballooning instability and finally trigger the following magnetic reconnection and eruptions. These physical scenarios may provide us a new viewpoint to understand the nature and origin of solar flares.

The relation between warm absorber (WA) outflows of AGN and nuclear obscuration activities caused by optically-thick clouds (obscurers) crossing the line of sight is unclear. NGC 3227 is a suitable target to study the properties of both WAs and obscurers, because it matches the following selection criteria: WAs in both ultraviolet (UV) and X-rays, suitably variable, bright in UV and X-rays, good archival spectra for comparing with the obscured spectra. To investigate WAs and obscurers of NGC~3227, we used a broadband spectral-energy-distribution model built in our Paper I and the photoionization code of SPEX software to fit archival XMM-Newton and NuSTAR observations in 2006 and 2016. Using unobscured observations, we find four WAs with different ionization states (log$\xi$ [erg cm/s]~-1.0, 2.0, 2.5, 3.0). The highest-ionization WA has a higher hydrogen column density (~$10^{22}$/cm$^2$) than the other three WAs (~$10^{21}$/cm$^2$). Their outflow velocities range from 100 to 1300 km/s, and show a positive correlation with the ionization parameter. These WAs are estimated to be between the outer broad-line-region (BLR) and the narrow line region. Besides, we find an X-ray obscuration event in 2006, which was missed by previous studies. It can be explained by a single obscurer. We also study the previously published obscuration event in 2016, which needs two obscurers in the fit. A high-ionization obscurer (log$\xi$~2.80; covering factor $C_f$~30%) only appears in 2016, which has a high column density (~$10^{23}$/cm$^2$). A low-ionization obscurer (log$\xi$~1.0-1.9; $C_f$~20%-50%) exists in both 2006 and 2016, which has a lower column density (~$10^{22}$/cm$^2$). These obscurers are estimated to be in the BLR by their crossing time of transverse motions. The obscurers and WAs of NGC 3227 have different distances and number densities, which indicate that they might have different origins.

The Cherenkov Telescopic Array (CTA), the next-generation ground-based gamma-ray observatory, will have unprecedented sensitivity, providing answers to open questions in gamma-ray cosmology and fundamental physics. Using simulations of active galactic nuclei observations foreseen in the CTA Key Science Program, we find that CTA will measure gamma-ray absorption on the extragalactic background light with a statistical error below 15% up to the redshift of 2 and detect or establish limits on gamma halos induced by the intergalactic magnetic field of at least 0.3 pG. Extragalactic observations using CTA also demonstrate the potential for testing physics beyond the Standard Model. The best state-of-the-art constraints on the Lorentz invariance violation from astronomical gamma-ray observations will be improved at least two- to threefold. CTA will also probe the parameter space where axion-like particles can represent a significant proportion - if not all - of dark matter. Joint multiwavelength and multimessenger observations, carried out together with other future observatories, will further foster the growth of gamma-ray cosmology.

The taxonomic classification of asteroids has been mostly based on spectroscopic observations with wavelengths spanning from the VIS to the NIR. VIS-NIR spectra of $\sim$2500 asteroids have been obtained since the 1970s; the SDSS MOC 4 was released with $\sim$4 $\times$ 10$^{5}$ measurements of asteroid positions and colors in the early 2000s. A number of works then devised methods to classify these data within the framework of existing taxonomic systems. Some of these works, however, used 2D parameter space that displayed a continuous distribution of clouds of data points resulting in boundaries that were artificially defined. We introduce here a more advanced method to classify asteroids based on existing systems. This approach is simply represented by a triplet of SDSS colors. The distributions and memberships of each taxonomic type are determined by machine learning methods in the form of both unsupervised and semi-supervised learning. We apply our scheme to MOC 4 calibrated with VIS-NIR reflectance spectra. We successfully separate seven different taxonomy classifications with which we have a sufficient number of spectroscopic datasets. We found the overlapping regions of taxonomic types in a 2D plane were separated with relatively clear boundaries in the 3D space newly defined in this work. Our scheme explicitly discriminates between different taxonomic types, which is an improvement over existing systems. This new method for taxonomic classification has a great deal of scalability for asteroid research, such as space weathering in the S-complex, and the origin and evolution of asteroid families. We present the structure of the asteroid belt, and describe the orbital distribution based on our newly assigned taxonomic classifications. It is also possible to extend the methods presented here to other photometric systems, such as the Johnson-Cousins and LSST filter systems.

We analyse quasi-periodic pulsations (QPP) detected in the microwave and decimeter radio emission of the SOL2017-09-05T07:04 solar flare, using simultaneous observations by the Siberian Radioheliograph 48 (SRH-48, 4-8 GHz) and Mingantu Spectral Radioheliograph (MUSER-I, 0.4-2 GHz). The microwave emission was broadband with a typical gyrosynchrotron spectrum, while a quasi-periodic enhancement of the decimetric emission appeared in a narrow spectral band (500-700 MHz), consistent with the coherent plasma emission mechanism. The periodicity that we found in microwaves is about 30 s, coming from a compact loop-like source with a typical height of about 31 Mm. The decimetric emission demonstrated a periodicity about 6 s. We suggested a qualitative scenario linking the QPPs observed in both incoherent and coherent spectral bands and their generation mechanisms. The properties of the QPPs found in the microwave signal are typical for perturbations of the flare loop by the standing sausage mode of a fast magnetohydrodynamic (MHD) wave. Our analysis indicated that this sausage-oscillating flare loop was the primary source of oscillations in the discussed event. The suggested scenario is that a fundamental sausage harmonic is the dominant cause for the observed QPPs in the microwave emission. The initiation of oscillations in the decimetric emission is caused by the third sausage harmonic via periodic and nonlinear triggering of the acceleration processes in the current sheets, formed at the interface between the sausage-oscillating flare loop and the external coronal loop that extended to higher altitudes. Our results demonstrate the possible role of MHD wave processes in the release and transport of energy during solar flares, linking coherent and incoherent radio emission mechanisms.

We consider the issue of excessive TTV noise observed for the exoplanet HD 189733 b. Trying to explain it through the host star photospheric activity, we model the stellar surface brightness as a random field, then characterize statistical properties of the resulting transit signal perturbation and compute individual corrections to transit timings uncertainties. We find that possible effect of the photospheric brightness field can explain only a minor portion ($\sim 10$ s) of the observed ($\sim 70$ s) TTV excess of HD 189733, suggesting that the rest should be attributed to other sources. Regarding the photospheric pattern, we place an upper limit of $\sim 0.01$ on the combination $\varepsilon_{\rm cell} r_{\rm cell}$, where $\varepsilon_{\rm cell}$ is the relative magnitude of brightness variations, and $r_{\rm cell}$ is the geometric cellularity scale (relative to star radius).

We investigate the origins of cold sub-Saturns (CSS), an exoplanetary population inferred from microlensing surveys. If confirmed, these planets would rebut a theorised gap in planets' mass distribution between those of Neptune and Jupiter caused by the rapid runaway accretion of super-critical cores. In an attempt to resolve this theoretical-observational disparity, we examine the outcomes of giant collisions between sub-critical protoplanets. Due to the secular interaction among protoplanets, these events may occur in rapidly depleting discs. We show that impactors ~ 5% the mass of near-runaway envelopes around massive cores can efficiently remove these envelopes entirely via a thermally-driven super-Eddington wind emanating from the core itself, in contrast with the stellar Parker winds usually considered. After a brief cooling phase, the merged cores resume accretion. But, the evolution timescale of transitional discs is too brief for the cores to acquire sufficiently massive envelopes to undergo runaway accretion despite their large combined masses. Consequently, these events lead to the emergence of CSS without their transformation into gas giants. We show that these results are robust for a wide range of disc densities, grain opacities and silicate abundance in the envelope. Our fiducial case reproduces CSS with heavy (>= 30 M_Earth) cores and less massive (a few M_Earth) sub-critical envelopes. We also investigate the other limiting cases, where continuous mergers of comparable-mass cores yield CSS with wider ranges of core-to-envelope mass ratios and envelope opacities. Our results indicate that it is possible for CSS and Uranus and Neptune to emerge within the framework of well studied processes and they may be more common than previously postulated.

Studies have shown that remnants of destroyed planets and debris-disk planetesimals can survive the volatile evolution of their host stars into white dwarfs, but detection of intact planetary bodies around white dwarfs are few. Simulations predict that planets in Jupiter-like orbits around stars of $\lt 8 M_\odot$ avoid being destroyed by the strong tidal forces of their stellar host, but as yet there has been no observational confirmation of such a survivor. Here we report on the non-detection of a main-sequence lens star in the microlensing event MOA-2010-BLG-477Lb using near-infrared observations from the Keck Observatory. We determine this system contains a $0.53\pm0.11$ solar mass white dwarf host orbited by a $1.4 \pm 0.3$ Jupiter mass planet with a separation on the plane of the sky of $2.8\pm 0.5$ AU, which implies a semi-major axis larger than this. This system is evidence that planets around white dwarfs can survive the giant and asymptotic giant phases of their host's evolution, and supports the prediction that over half of white dwarfs are predicted to have Jovian planetary companions. Located at approximately 2.0 kpc toward the center of our Galaxy, it likely represents an analog to the end stages of the Sun and Jupiter in our own Solar System.

To investigate star-forming activities in early-type galaxies, we select a sample of 52 star-forming S0 galaxies (SFS0s) from the SDSS-IV MaNGA survey. We find that SFS0s have smaller stellar mass compared to normal S0s in MaNGA. After matching the stellar mass to select the control sample, we find that the mean S\'{e}rsic index of SFS0s' bulges (1.76$\pm$0.21) is significantly smaller than that of the control sample (2.57$\pm$0.20), suggesting the existence of a pseudo bulge in SFS0s. After introducing the environmental information, SFS0s show smaller spin parameters in the field than in groups, while the control sample has no obvious difference in different environments, which may suggest different dynamical processes in SFS0s. Furthermore, with derived N/O and O/H abundance ratios, SFS0s in the field show nitrogen enrichment, providing evidence for the accretion of metal-poor gas in the field environment. To study the star formation relation, we show that the slope of the spatially resolved star formation main sequence is nearly 1.0 with MaNGA IFU data, confirming the self-regulation of star formation activities at the kpc scales.

Recent observations of the gamma-ray source HAWC J2227+610 by Tibet AS+MD and LHAASO confirm the special interest of this source as a galactic PeVatron candidate in the northern hemisphere. HAWC J2227+610 emits Very High Energy (VHE) gamma-rays up to 500 TeV, from a region coincident with molecular clouds and significantly displaced from the nearby pulsar J2229+6114. Even if this morphology favours an hadronic origin, both leptonic or hadronic models can describe the current VHE gamma-ray emission. The morphology of the source is not well constrained by the present measurements and a better characterisation would greatly help the understanding of the underlying particle acceleration mechanisms. The Cherenkov Telescope Array (CTA) will be the future most sensitive Imaging Atmospheric Cherenkov Telescope and, thanks to its unprecedented angular resolution, could contribute to better constrain the nature of this source. The present work investigates the potentiality of CTA to study the morphology and the spectrum of HAWC J2227+610. For this aim, the source is simulated assuming the hadronic model proposed by the Tibet AS+MD collaboration, recently fitted on multi-wavelength data, and two spatial templates associated to the source nearby molecular clouds. Different CTA layouts and observation times are considered. A 3D map based analysis shows that CTA is able to significantly detect the extension of the source and to attribute higher detection significance to the simulated molecular cloud template compared to the alternative one. CTA data does not allow to disentangle the hadronic and the leptonic emission models. However, it permits to correctly reproduce the simulated parent proton spectrum characterized by a $\sim$ 500 TeV cutoff.

Single-line spectroscopic binaries recently contribute to the stellar-mass black hole discovery, independently of the X-ray transient method. We report the identification of a single-line binary system LTD064402+245919, with an orbital period of 14.50 days. The observed component is a subgiant with a mass of 2.77$\pm$0.68M$_{\odot}$, radius 15.5$\pm$2.5R$_{\odot}$, effective temperature $T_{\rm eff}$ 4500$\pm$200K, and surface gravity log\emph{g} 2.5$\pm$0.25dex. The discovery makes use of the LAMOST time-domain (LAMOST-TD) and ZTF survey. Our general-purpose software pipeline applies the Lomb-Scargle periodogram to determine the orbital period and uses machine-learning to classify the variable type from the folded light curves. We apply a combined model to estimate the orbital parameters from both the light and radial velocity curves, taking constraints on the primary star mass, mass function, and detection limit of secondary luminosity into consideration. We obtain a radial velocity semi-amplitude of 44.6$\pm$1.5 km s$^{-1}$, mass ratio of 0.73$\pm$0.07, and an undetected component mass of 2.02$\pm$0.49M$_{\odot}$ when the type of the undetected component is not set. We conclude that the inclination is not well constrained, and that the secondary mass is larger than 1M$_{\odot}$ when the undetected component is modelled as a compact object. According to our investigations using an MCMC simulation, increasing the spectra SNR by a factor of 3 would enable the secondary light to be distinguished (if present). The algorithm and software in this work are able to serve as general-purpose tools for the identification of compact objects quiescent in X-rays.

The search for Galactic pevatrons is now a well-identified key science project of all instruments operating in the very-high-energy domain. Indeed, in this energy range, the detection of gamma rays clearly indicates that efficient particle acceleration is taking place, and observations can thus help identify which astrophysical sources can energize particles up to the $\sim$PeV range, thus being $pevatrons$. In the search for the origin of Galactic cosmic rays (CRs), the PeV range is an important milestone, since the sources of Galactic CRs are expected to accelerate PeV particles. This is how the central scientific goal that is 'solving the mystery of the origin of CRs' has often been distorted into 'finding (a) pevatron(s)'. Since supernova remnants (SNRs) are often cited as the most likely candidates for the origin of CRs, 'finding (a) pevatron(s)' has often become 'confirming that SNRs are pevatrons'. Pleasingly, the first detection(s) of pevatron(s) were not associated to SNRs. Moreover, all clearly detected SNRs have yet revealed to not be pevatrons, and the detection from VHE gamma rays from regions unassociated with SNRs, are reminding us that other astrophysical sites might well be pevatrons. This short review aims at highlighting a few important results on the search for Galactic pevatrons.

In this article we present a formalism to incorporate the partial-sky maps to the Gibbs ILC algorithm to estimate the joint posterior density of the Cosmic Microwave Background (CMB) signal and the theoretical CMB angular power spectrum given the observed CMB maps. In order to generate the partial-sky maps we mask all the observed CMB maps provided by WMAP and Planck satellite mission using a Gaussian smoothed mask formed on the basis of thermal dust emissions in Planck 353 GHz map. The central galactic region from all the input maps is removed after the application of the smoothed mask. While implementing the Gibbs ILC method on the partial-sky maps, we convert the partial-sky cleaned angular power spectrum to the full-sky angular power spectrum using the mode-mode coupling matrix estimated from the smoothed mask. The main products of our analysis are partial-sky cleaned best-fit CMB map and an estimate of the underlying full-sky theoretical CMB angular power spectrum along with their error estimates. We validate the methodology by performing detailed Monte Carlo simulations after using realistic models of foregrounds and detector noise consistent with the WMAP and Planck frequency channels used in our analysis. We can estimate the posterior density and full-sky theoretical CMB angular power spectrum, without any need to explicitly model the foreground components, from partial-sky maps using our method. Another important feature of this method is that the power spectrum results along with the error estimates can be directly used for cosmological parameter estimations.

The polarization modulator unit for the low-frequency telescope in LiteBIRD employs an achromatic half-wave plate (AHWP). It consists of five layers of a-cut sapphire plate, which are stacked based on a Pancharatnam recipe. In this way, the retardance of the AHWP is a half-wave over a bandwidth of 34-161 GHz. The diameter of a single sapphire plate is about 500 mm and the thickness is 5 mm. When a large diameter AHWP is used for a space mission, it is important for the AHWP to survive launch vibration. A preliminary study indicates that the five-layer stacked HWP has a risk of breakage at the launch unless the five layers are glued together and mechanically treated as one disk. We report our investigation using a sodium silicate solution that can bond the sapphire plates. This technique has been previously investigated as a candidate of cryogenic bonding for a mirror material, including sapphire, of the gravitational wave experiments: LIGO, VIRGO, and KAGRA. We experimentally studied the mechanical strength of the bonded interface for two different surface conditions: polished and unpolished. We demonstrated that the tensile and shear strength > 20 MPa for samples with a polished surface. This satisfied the requirement of 5.5 MPa derived from the mechanical simulation assuming a launch load of 30G. We identified that samples glued on a polished surface exhibit higher strength than unpolished ones by a factor of 2 for tensile and 18 for shear strength. We measured the millimeter-wave transmittance between 90 and 140 GHz using sapphire plates with a diameter of 50 mm before and after bonding. We did not find any optical effects caused by the bonded interface within 2% error in transmittance, which originates from the measurement system.

We showcase a tool suite that enables the fitting of soft X-ray spectra in active galactic nuclei (AGN), without the need for specialist software, allowing access to AGN physics for school students. While these standardised Python tools were useful for measuring velocities (Note I), they offered significantly fewer capabilities for radiative recombination continua (RRC), and R and G ratios, utilised to obtain the internal plasma properties within the out owing wind seen in NGC 4151. Although further work is required for these tools to be used in outreach projects, we present findings of the plasma temperature and density in NGC 4151 spanning a 15 year period.

In recent years, the proposal that there is a large population of primordial black holes living in dense clusters has been gaining popularity. One natural consequence of these dense clusters will be that the black holes inside will gravitationally scatter off each other in hyperbolic encounters, emitting gravitational waves that can be observed by current detectors. In this paper we will derive how to compute the gravitational waves emitted by black holes in hyperbolic orbits, taking into account up to leading order spin effects. We will then study the signal these waves leave in the network of gravitational wave detectors currently on Earth. Using the properties of the signal, we will detail the data processing techniques that can be used to make it stand above the detector noise. Finally, we will look for these signals from hyperbolic encounters in the publicly available LIGO-Virgo data. For this purpose we will develop a two step trigger. The first step of the trigger will be based on looking for correlations between detectors in the time-frequency domain. The second step of the trigger will make use of a residual convolutional neural network, trained with the theoretical predictions for the signal, to look for hyperbolic encounters. With this trigger we find 8 hyperbolic encounter candidates in the 15.3 days of public data analyzed. Some of these candidates are promising, but the total number of candidates found is consistent with the number of false alarms expected from our trigger.

We construct Saturn equatorial neutral temperature and density profiles of H, H$_2$, He, and CH$_4$, between 10$^{-12}$ and 1 bar using measurements from Cassini's Ion Neutral Mass Spectrometer (INMS) taken during the spacecraft's final plunge into Saturn's atmosphere on 15 September 2017, combined with previous deeper atmospheric measurements from the Cassini Composite InfraRed Spectrometer (CIRS) and from the UltraViolet Imaging Spectrograph (UVIS). These neutral profiles are fed into an energy deposition model employing soft X-ray and Extreme UltraViolet (EUV) solar fluxes at a range of spectral resolutions ($\Delta\lambda=4\times10^{-3}$ nm to 1 nm) assembled from TIMED/SEE, from SOHO/SUMER, and from the Whole Heliosphere Interval (WHI) quiet Sun campaign. Our energy deposition model calculates ion production rate profiles through photo-ionisation and electron-impact ionisation processes, as well as rates of photo-dissociation of CH$_4$. The ion reaction rate profiles we determine are important to obtain accurate ion density profiles, meanwhile methane photo-dissociation is key to initiate complex organic chemical processes. We assess the importance of spectral resolution in the energy deposition model by using a high-resolution H$_2$ photo-absorption cross section, which has the effect of producing additional ionisation peaks near 800 km altitude. We find that these peaks are still formed when using low-resolution ($\Delta\lambda=1$ nm) or mid-resolution ($\Delta\lambda=0.1$ nm) solar spectra, as long as high-resolution cross sections are included in the model.

We use Gaia EDR3 astrometry to propose that a dynamical interaction between the multiple system $\theta^1$ Ori C and $\theta^1$ Ori F ejected the latter as a walkaway star ~1100 years ago (without deceleration) or somewhat later (with a more likely deceleration included). It is unclear whether the final 3-D velocity of $\theta^1$ Ori F will be large enough to escape the Orion nebula cluster.

Massive steam and CO$_2$ atmospheres have been proposed for magma ocean outgassing of Earth and terrestrial planets. Yet formation of such atmospheres depends on volatile exchange with the molten interior, governed by volatile solubilities and redox reactions. We determine the evolution of magma ocean-atmosphere systems for a range of oxygen fugacity, C/H ratio and hydrogen budget that include redox reactions for hydrogen (H$_2$-H$_2$O), carbon (CO-CO$_2$), methane (CH$_4$), and solubility laws for H$_2$O and CO$_2$ for silicate melts. We find that small initial budgets of hydrogen, high C/H ratios, and oxidising conditions, suppress outgassing of hydrogen until the late stages of magma ocean crystallisation. Hence early atmospheres in equilibrium with magma oceans are dominantly carbon-rich, and specifically CO-rich except at the most oxidising conditions. The high solubility of H$_2$O limits its outgassing to melt fractions below $\sim$30\%; the fraction at which the mantle transitions from vigorous- to sluggish convection with melt percolation. Inefficient melt percolation may prevent further crystallisation before a surface lid forms, which can delay or prevent most water ($>$75\%) from outgassing, thereby hindering the formation of surface oceans. Alternatively, with efficient melt percolation or mechanism to break the lid, a molten surface could be maintained, enabling the transition from a CO-rich to water-rich atmosphere during the late stage of crystallisation. We conclude that a large fraction of the water delivered to terrestrial planets over their accretion histories may be safely harboured in their interiors during the magma ocean phase and outgas only gradually over geological timescales.

Aims. We describe the first s-process post-processing models for asymptotic giant branch (AGB) stars of masses 3, 4 and 5 M at solar metallicity (Z=0.018) computed using the input from the stellar evolutionary code aton. Methods. The models are computed with the new code snuppat(S-process NUcleosynthesis Post-Processing code for aton), including an advective scheme for the convective overshoot that leads to the formation of the main neutron source, 13C. Each model is post-processed with 3 different values of the free overshoot parameter. Included in the code snuppat is the novel Patankar-Euler-Deflhard explicit numerical solver, that we use to solve the nuclear network system of differential equations. Results. The results are compared to those from other s-process nucleosynthesis codes (Monash,fruity, and NuGrid), as well as observations of s-process enhancement in AGB stars, planetary nebulae, and barium stars. This comparison shows that the relatively high abundance of12C in the He-rich intershell in aton results in as-process abundance pattern that favours the second over the first s-process peak for all the masses explored. Also, our choice of an advective as opposed to diffusive numerical scheme for the convective overshoot results in significants-process nucleosynthesis also for the 5 M models, which may be in contradiction with observations.

We revisit the moving groups (MGs) in the solar neighborhood with a sample of 91969 nearby stars constructed from LAMOST DR7. Using the wavelet technique and Monte Carlo simulations, five MGs together with a new candidate located at $V \simeq$ -130 km s$^{-1}$ are detected simultaneously in $V-\sqrt{U^2+2V^2}$ space. Taking into account the other known MGs, we conclude that MGs in the Galactic disk are spaced by approximately 15 $\sim$ 25 km s$^{-1}$ along $V$ velocity. The origin of detected MGs is analysed through the distributions of [Fe/H]$-$[Mg/Fe] and ages. Our results support attributing the origin to the continuous resonant mechanisms probably induced by the bar or spiral arms of the Milky Way.

We observed the predicted outburst of the Alpha Monocerotid (AMO) meteor shower on 2019 November 22 with our modernized video and photographic cameras. Due to the short duration and moderate intensity of the outburst, atmospheric trajectories and radiants were obtained for only ten meteors, seven of which included velocities, magnitudes, and orbits. In addition, one incomplete video spectrum was captured. The radiants and orbits were found to be compatible with that of the 1995 outburst. The spectrum confirmed that AMO meteoroids are deficient in sodium. Unlike any other meteor shower, meteor end heights were found to be distributed along a constant level of 90 km for all meteors with magnitudes between +4 and -2 and with atmospheric trajectory lengths up to 40 km. We propose that Alpha Monocerotids were formed from a devolatilized and fragile cometary crust composed from relatively large fundamental grains.

We present a model of a viscously evolving accretion disc around a magnetized neutron star. The model features the varying outer radius of the hot ionized part of the disc due to cooling and the varying inner radius of the disc due to interaction with the magnetosphere. It also includes hindering of accretion on the neutron star because of the centrifugal barrier and irradiation of the outer disc and companion star by X-rays from the neutron star and disc. When setting inner boundary conditions, we take into account that processes at the inner disc occur on a time scale much less than the viscous time scale of the whole disc. We consider three types of outflow from the disc inner edge: zero outflow, one based on MHD calculations, and a very efficient propeller mechanism. The light curves of an X-ray transient after the outburst peak can be calculated by a corresponding, publicly available code. We compare observed light curves of the 2013 burst of Aql X-1 in X-ray and optical bands with modeled ones. We find that the fast drop of the $0.3-10$ keV flux can be solely explained by a radial shrinking of the hot disc. At the same time, models with the neutron star magnetic field $>10^8$ G have better fits because the accretion efficiency behaviour emphasizes the 'knee' on the light curve. We also find that a plato emission can be produced by a disc-reservoir with stalled accretion.

We report on a NuSTAR observation of the young, energetic pulsar PSR J1617-5055. Parkes Observatory 3 GHz radio observations of the pulsar (taken about 7 years before the NuSTAR observations) are also reported here. NuSTAR detected pulsations at a frequency of $f\approx14.4$ Hz ($P\approx69.44$ ms) and, in addition, the observation was long enough to measure the source's frequency derivative, $\dot{f}\approx-2.8\times10^{-11}$ Hz s$^{-1}$. We find that the pulsar shows one peak per period at both hard X-ray and radio wavelengths, but that the hard X-ray pulse is broader (having a duty cycle of $\sim 0.7$), than the radio pulse (having a duty cycle of $\sim 0.08$). Additionally, the radio pulse is strongly linearly polarized. J1617's phase-integrated hard X-ray spectrum is well fit by an absorbed power-law model, with a photon index $\Gamma=1.59\pm 0.02$. The hard X-ray pulsations are well described by three Fourier harmonics, and have a pulsed fraction that increases with energy. We also fit the phase-resolved NuSTAR spectra with an absorbed power-law model in five phase bins and find that the photon index varies with phase from $\Gamma = 1.52\pm 0.03$ at phases around the flux maximum to $\Gamma=1.79\pm 0.06$ around the flux minimum. Lastly, we compare our results with other pulsars whose magnetospheric emission is detected at hard X-ray energies and find that, similar to previous studies, J1617's hard X-ray properties are more similar to the MeV pulsars than the GeV pulsars.

MESSENGER observations suggest a magma ocean formed on proto-Mercury, during which evaporation of metals and outgassing of C- and H-bearing volatiles produced an early atmosphere. Atmospheric escape subsequently occurred by plasma heating, photoevaporation, Jeans escape, and photoionization. To quantify atmospheric loss, we combine constraints on the lifetime of surficial melt, melt composition, and atmospheric composition. Consideration of two initial Mercury sizes and four magma ocean compositions determine the atmospheric speciation at a given surface temperature. A coupled interior-atmosphere model determines the cooling rate and therefore the lifetime of surficial melt. Combining the melt lifetime and escape flux calculations provide estimates for the total mass loss from early Mercury. Loss rates by Jeans escape are negligible. Plasma heating and photoionization are limited by homopause diffusion rates of $\sim10^{6}$ kg/s. Loss by photoevaporation depends on the timing of Mercury formation and assumed heating efficiency and ranges from $\sim10^{6.6}$ to $\sim10^{9.6}$ kg/s. The material for photoevaporation is sourced from below the homopause and is therefore energy-limited rather than diffusion-limited. The timescale for efficient interior-atmosphere chemical exchange is less than ten thousand years. Therefore, escape processes only account for an equivalent loss of less than 2.3 km of crust ($0.3\%$ of Mercury's mass). Accordingly, $\leq0.02\%$ of the total mass of H$_2$O and Na is lost. Therefore, cumulative loss cannot significantly modify Mercury's bulk mantle composition during the magma ocean stage. Mercury's high core:mantle ratio and volatile-rich surface may instead reflect chemical variations in its building blocks resulting from its solar-proximal accretion environment.

As more missions attempt to directly image and characterize exoplanets orbiting nearby sun-like stars, advance characterization of possible contaminating background sources becomes more important and can impact target selection. This paper describes an exploration of the Hubble Source catalog, Gaia catalog and Besancon galaxy simulations in order to determine the likelihood of having a contaminating source in the background of a set of high proper motion stars in the expected timeframe of observations for the Nancy Grace Roman Space Telescope Coronagraphic Instrument. The analysis shows that for most of the targets, there is a very low possibility of a star falling within the CGI field of view, but that at low galactic latitudes where there is a greater density of sources, faint stellar background sources could be a concern.

We report the discovery of emerged tidal tails around open cluster IC 4756 ($\sim$ 1 Gyr) based on 644 members identified from {\it Gaia} EDR3. Three-dimensional spatial positions, two-dimensional tangential velocities $\left( x, y, z, \kappa \cdot \mu_{\alpha}^{*}/\varpi, \kappa \cdot \mu_{\delta}/\varpi \right)$ are utilized to determine the co-moving member candidates of IC 4756. Using a Bayesian method, we correct the distance for each cluster member. Two tidal tails extend up to 180 pc and display a S-shape in $X^{\prime}Y^{\prime}$ space (Cartesian coordinates focused on cluster center). A clean sequence of our members in Color-Absolute-Magnitude Diagram (CAMD) indicates the coeval population and matches perfectly with the PARSEC isochrone with age from Bossini et al. (2019). Mass segregation is detected in this cluster as well. Finally, we derive the tidal radius and core radius of IC 4756 about $12.13$ pc and $4.33 \pm 0.75$ pc, respectively.

Context. Relativistic jets are ubiquitous in astrophysics. High Mass Microquasars (HMMQs) are useful labs to study such jets, as they are relatively close and evolve over observable time scales. The ambient medium into which the jet propagates is, however, far from homogeneous. Corresponding simulation studies to date consider various forms of a wind-shaped ambient medium but typically neglect radiative cooling and relativistic effects. Aims. We investigate the dynamical and structural effects of radiative losses and system parameters on relativistic jets in HMMQs, from the jet launch to its propagation over several tens of orbital separations. Methods. We use 3D relativistic hydrodynamical simulations including parameterized radiative cooling derived from relativistic thermal plasma distribution to carry out parameter studies around two fiducial cases inspired by Cygnus X-1 and Cygnus X-3. Results. Radiative losses are found to be more relevant in Cygnus X-3 than Cygnus X-1. Varying jet power, jet temperature, or the wind of the donor star tends to have a larger impact at early times, when the jet forms and instabilities initially develop, than at later times when the jet has reached a turbulent state. Conclusions. Radiative losses may be dynamically and structurally relevant at least in Cygnus X-3 case, and thus should be examined in more detail.

Chirality, or handedness, enters astrophysics in three distinct ways. Magnetic field and vortex lines tend to be helical and have a systematic twist in the northern and southern hemispheres of a star or a galaxy. Helicity is here driven by external factors. Chirality can also enter at the microphysical level and can then be traced back to the parity-breaking weak force. Finally, chirality can arise spontaneously, but this requires not only the presence of an instability, but also the action of nonlinearity. Examples can be found both in magnetohydrodynamics and in astrobiology, where homochirality among biomolecules probably got established at the origin of life. In this review, all three types of chirality production will be explored and compared.

Context. We aim to determine the capabilities of the MEGARA@GTC instrument integral-field unit to study stellar populations and exploit its combination of high spectral (R \sim 6,000, 12,000 and 20,000) and spatial (0.62") resolutions within its 12.5"x11.3" field of view. Aims. We pursue to establish a systematic method through which we can determine the properties of the stellar populations in the observations made with MEGARA, more specifically within the MEGADES legacy project and, for this paper, those of the stellar populations of NGC 7025. Methods. We use MEGARA observations of galaxy NGC 7025. We apply different approaches to estimate the properties of the stellar populations with the highest possible certainty. Numerous tests were also performed to check the reliability of the study. Results. All the studies we conduct (both full spectral fitting and absorption line indices) on the stellar populations of NGC 7025 indicate that the stars that form its bulge have supersolar metallicity and considerably old ages (\sim 10 Gyr), in general. We determined that the bulge of NGC 7025 has a mild negative mass-weighted age gradient using three different combinations of MEGARA spectral setups. Regarding its more detailed star formation history, our results indicate that, besides a rather constant star formation at early epochs, a peak in formation history of the stars in the bulge is also found 3.5-4.5 Gyr ago, partly explaining the mass-weighted age gradients measured. Conclusions. The scenario presented in NGC 7025 is that of an isolated galaxy under secular evolution that about 3.5-4.5 Gyr ago likely experimented a minor merger (mass ratio 1/10) that induced an increase in star formation and also perturbed the morphology of its outer disc. We report on different lessons learned for the ongoing exploitation of the MEGADES survey with GTC.

Spectroscopic eclipse observations, like those possible with the James Webb Space Telescope, should enable 3D mapping of exoplanet daysides. However, fully-flexible 3D planet models are overly complex for the data and computationally infeasible for data-fitting purposes. Here, we present ThERESA, a method to retrieve the 3D thermal structure of an exoplanet from eclipse observations by first retrieving 2D thermal maps at each wavelength and then placing them vertically in the atmosphere. This approach allows the 3D model to include complex thermal structures with a manageable number of parameters, hastening fit convergence and limiting overfitting. An analysis runs in a matter of days. We enforce consistency of the 3D model by comparing vertical placement of the 2D maps with their corresponding contribution functions. To test this approach, we generated a synthetic JWST NIRISS-like observation of a single hot-Jupiter eclipse using a global circulation model of WASP-76b and retrieved its 3D thermal structure. We find that a model which places the 2D maps at different depths depending on latitude and longitude is preferred over a model with a single pressure for each 2D map, indicating that ThERESA is able to retrieve 3D atmospheric structure from JWST observations. We successfully recover the temperatures of the planet's dayside, the eastward shift of its hotspot, and the thermal inversion. ThERESA is open-source and publicly available as a tool for the community.

It is well established that bars evolve significantly after they form in galaxy discs, often changing shape both in and out of the disc plane. In some cases they may bend or buckle out of the disc plane resulting in the formation of boxy/peanut/x-shape bulges. In this paper we show that the dark matter halo shape affects bar formation and buckling. We have performed N-body simulations of bar buckling in non-spherical dark matter halos and traced bar evolution for 8 Gyr. We find that bar formation is delayed in oblate halos, resulting in delayed buckling whereas bars form earlier in prolate halos leading to earlier buckling. However, the duration of first buckling remains almost comparable. All the models show two buckling events but the most extreme prolate halo exhibits three distinct buckling features. Bars in prolate halos also show buckling signatures for the longest duration compared to spherical and oblate halos. Since ongoing buckling events are rarely observed, our study suggests that most barred galaxies may have more oblate or spherical halos rather than prolate halos. Our measurement of BPX structures also shows that prolate halos promote bar thickening and disc heating more than oblate and spherical halos.

We explore how chromatic RFI flags affect 21-cm power spectrum measurements. We particularly study flags that are coarser than the analysis resolution. We find that such RFI flags produce excess power in the EoR window in much the same way as residual RFI. We use simulations to explain this as a result of chromatic disruptions in the interferometric sampling function of the array. We also show that without modifying current flagging strategies or implementing extremely accurate foreground subtraction, 21-cm EoR experiments will fail to make a significant detection.

We present BlendHunter, a proof-of-concept for a deep transfer learning based approach for the automated and robust identification of blended sources in galaxy survey data. We take the VGG-16 network with pre-trained convolutional layers and train the fully connected layers on parametric models of COSMOS images. We test the efficacy of the transfer learning by taking the weights learned on the parametric models and using them to identify blends in more realistic CFIS-like images. We compare the performance of this method to SEP (a Python implementation of SExtractor) as function of noise level and the separation between sources. We find that BlendHunter outperforms SEP by $\sim 15\%$ in terms of classification accuracy for close blends ($<10$ pixel separation between sources) regardless of the noise level used for training. Additionally, the method provides consistent results to SEP for distant blends ($\geq10$ pixel separation between sources) provided the network is trained on data with a relatively close noise standard deviation to the target images. The code and data have been made publicly available to ensure the reproducibility of the results.

Working in an idealised framework in which a series of phases of evolution defined by the second slow-roll parameter $\eta$ are matched together, we calculate the reduced bispectrum, $f_{\rm NL}$, for models of inflation with a large peak in their primordial power spectra. We find $f_{\rm NL}$ is typically approximately constant over scales at which the peak is located, and provide an analytic approximation for this value. This allows us to identify the conditions under which $f_{\rm NL}$ is large enough to have a significant impact on the resulting production of primordial black holes (PBHs) and scalar induced gravitational waves (SIGWs). Together with analytic formulae for the gradient of the rise and fall in the power spectrum, this provides a toolkit for designing or quickly analysing inflationary models that produce PBHs and SIGWs.

Innovative developments in data processing, archiving, analysis, and visualization are nowadays unavoidable to deal with the data deluge expected in next-generation facilities for radio astronomy, such as the Square Kilometre Array (SKA) and its precursors. In this context, the integration of source extraction and analysis algorithms into data visualization tools could significantly improve and speed up the cataloguing process of large area surveys, boosting astronomer productivity and shortening publication time. To this aim, we are developing a visual analytic platform (CIRASA) for advanced source finding and classification, integrating state-of-the-art tools, such as the CAESAR source finder, the ViaLactea Visual Analytic (VLVA) and Knowledge Base (VLKB). In this work, we present the project objectives and the platform architecture, focusing on the implemented source finding services.

The presence of radiatively driven outflows is well established in ultraluminous X-ray sources (ULXs). These outflows are optically thick and can reprocess a significant fraction of the accretion luminosity. Assuming isotropic emission, escaping radiation from the outflow's photosphere has the potential to irradiate the outer disc. Here, we explore how the atmosphere of the outer disc would respond to such irradiation, and specifically whether unstable heating may lead to significant mass loss via thermally-driven winds. We find that, for a range of physically relevant system parameters, this mass loss may actually switch off the inflow entirely and potentially drive limit-cycle behaviour (likely modulated on the timescale of the outer disc). In ULXs harbouring neutron stars, magnetic fields tend to have a slight destabilizing effect; for the strongest magnetic fields and highest accretion rates, this can push otherwise stable systems into the unstable regime. We explore the prevalence of the instability in a simulated sample of ULXs obtained from a binary population synthesis calculation. We find that almost all neutron star and black hole ULXs with Eddington-scaled accretion rates of $\dot{m}_0 < 100$ should be able to drive powerful outflows from their outer discs. Several known ULXs are expected to lie in this regime; the persistence of accretion in these sources implies the irradiation may be anisotropic which can be reconciled with the inferred reprocessed (optical) emission if some of this originates in the wind photosphere or irradiation of the secondary star.

Strongly non-geodesic, or rapidly turning, trajectories in multifield inflation have attracted much interest recently from both theoretical and phenomenological perspectives. For example, such trajectories can evade the eta-problem in fat inflation, in which all inflatons are heavier than the Hubble scale. Most models with large turning rates in the literature are formulated as effective field theories. In this paper we investigate rapid-turn inflation in supergravity as a first step toward understanding them in string theory. We find that large turning rates can be generated in a wide class of models, at the cost of high field space curvature. However, in these models one field is always tachyonic along the inflationary trajectory, implying that fat inflation may not be viable in supergravity. The high curvatures underscore the difficulty of obtaining rapid-turn inflation in realistic string-theoretical models.

Classical fields are ubiquitous in theoretical physics. They find applications in almost all areas of physics, from condensed matter and particle physics to cosmology and astrophysics. Scalar fields, in particular, can give rise to confined structures, such as boson stars, oscillatons or Q-balls. These objects are interesting hypothetical new "dark matter stars", but also good descriptions of dark matter cores when the fields are ultralight. In this thesis, we study the dynamical response of such confined bosonic structures when excited by external matter (stars, planets or black holes) in their vicinities. Such perturbers can either be piercing through the bosonic configuration or undergoing periodic motion at its center (e.g., binaries). Our setup can also efficiently describe the interaction between a moving massive black hole and the surrounding environment. It also depicts dark matter depletion as a reaction to an inspiralling binary within a dark matter core. Our results provide a complete picture of the interaction between black holes or stars and the ultralight dark matter core environment where they may live in. This thesis also deals with several classical field environmental effects on the motion (or, ultimately, the survival) of compact objects, like black holes.

Magnetars have already been a potential candidate as gravitational wave sources that could be detected by current and future terrestrial as well as ground based gravitational wave detectors. In this article, we focus on the gravitational wave emission from the distorted rotating neutron stars. The deformation is assumed to be symmetric around an axis that is perpendicular to the rotation axis. The form is applied in the context of a neutron star whose magnetic field has been deformed on its own. By introducing the effects from all magnetars in the Universe, based on various proposed magnetic field configurations, such as poloidal, toroidal; stochastic gravitational wave background (SGWB) can be generated. We choose to figure out exactly how the observations of the stochastic gravitational-wave background should be used to understand much more about physics correlated with the magnetar behaviour, based on the restriction on the eccentricity of the magnetar.

High energy solar flares and coronal mass ejections have the potential to destroy Earth's ground and satellite infrastructures, causing trillions of dollars in damage and mass human suffering. Destruction of these critical systems would disable power grids and satellites, crippling communications and transportation. This would lead to food shortages and an inability to respond to emergencies. A solution to this impending problem is proposed herein using satellites in solar orbit that continuously monitor the Sun, use artificial intelligence and machine learning to calculate the probability of massive solar explosions from this sensed data, and then signal defense mechanisms that will mitigate the threat. With modern technology there may be only safeguards that can be implemented with enough warning, which is why the best algorithm must be identified and continuously trained with existing and new data to maximize true positive rates while minimizing false negatives. This paper conducts a survey of current machine learning models using open source solar flare prediction data. The rise of edge computing allows machine learning hardware to be placed on the same satellites as the sensor arrays, saving critical time by not having to transmit remote sensing data across the vast distances of space. A system of systems approach will allow enough warning for safety measures to be put into place mitigating the risk of disaster.

This letter reports the first direct measurements of the permanent space-time strain component of the gravitational-wave memory effect, predicted by general relativity to accompany black hole merger events. By ignoring the details of how the memory effect develops over time, this approach circumvents the need for precise modeling of the non-memory, linear memory, and nonlinear memory signals. Applied to a selection 64 observations of black hole merger events in the LIGO/Virgo Gravitational Wave Transient Catalog, this analysis yields a mixture of probable detections and upper limits. These results are supported by both individual and ensemble statistical significance analyses, based on analyzing the merger events -- and a large number of time intervals shifted away from them -- in exactly the same way.

When a physicist says that a theory is fine-tuned, they mean that it must make a suspiciously precise assumption in order to explain a certain observation. This is evidence that the theory is deficient or incomplete. One particular case of fine-tuning is particularly striking. The data in question are not the precise measurements of cosmology or particle physics, but a more general feature of our universe: it supports the existence of life. This chapter reviews this Fine-Tuning of the Universe for Life.

We develop a systematic approach to realising linear detectors that saturate the Heisenberg limit. First, we consider the general constraints on the input-output transfer matrix of a linear detector. We then derive the physical realization of the most general transfer matrix using the quantum network synthesis technique, which allows for the inference of the physical setup directly from the input-output transfer matrix. By exploring the minimal realization which has the minimum number of internal modes, we show that such detectors that saturate the Heisenberg limit are internal squeezing schemes. Then, investigating the non-minimal realization, which is motivated by the parity-time symmetric systems, we arrive at the general quantum non-demolition measurement.

Most microwave readout architectures in quantum computing or sensing rely on a semiconductor amplifier at 4 K, typically a high-electron mobility transistor (HEMT). Despite its remarkable noise performance, a conventional HEMT dissipates several milliwatts of power, posing a practical challenge to scale up the number of qubits or sensors addressed in these architectures. As an alternative, we present an amplification chain consisting of a kinetic-inductance traveling-wave parametric amplifier (KI-TWPA) placed at 4 K, followed by a HEMT placed at 70 K, and demonstrate a chain-added noise $T_\Sigma = 6.3\pm0.5$ K between 3.5 and 5.5 GHz. While, in principle, any parametric amplifier can be quantum limited even at 4 K, in practice we find the KI-TWPA's performance limited by the temperature of its inputs, and by an excess of noise $T_\mathrm{ex} = 1.9$ K. The dissipation of the KI-TWPA's rf pump constitutes the main power load at 4 K and is about one percent that of a HEMT. These combined noise and power dissipation values pave the way for the KI-TWPA's use as a replacement for semiconductor amplifiers.