Papers for Wednesday, Jan 20 2021

We present the results of a search for stellar flares from stars neighbouring the target sources in the Kepler short cadence data. These flares have been discarded as contaminants in previous surveys and therefore provide an unexplored resource of flare events, in particular high energy events from faint stars. We have measured M dwarf flare energies up to 1.5$\times$10^35 erg, pushing the limit for flare energies measured using Kepler data. We have used our sample to study theflaring activity of wide binaries, finding that the lower mass counterpart in a wide binary flares more often at a given energy. Of the 4430 flares detected in our original search, 298 came from a neighbouring star, a rate of 6.7$\pm$0.4 per cent for the Kepler short cadence lightcurves. We have used our sample to estimate a 5.8$\pm$0.1 per cent rate of false positive flare events in studies using TESS short cadence data.

After the prediction of many sub- and super-Chandrasekhar (at least a dozen for the latter) limiting mass white dwarfs, hence apparently peculiar class of white dwarfs, from the observations of luminosity of type Ia supernovae, researchers have proposed various models to explain these two classes of white dwarfs separately. We earlier showed that these two peculiar classes of white dwarfs, along with the regular white dwarfs, can be explained by a single form of the f(R) gravity, whose effect is significant only in the high-density regime, and it almost vanishes in the low-density regime. However, since there is no direct detection of such white dwarfs, it is difficult to single out one specific theory from the zoo of modified theories of gravity. We discuss the possibility of direct detection of such white dwarfs in gravitational wave astronomy. It is well-known that in f(R) gravity, more than two polarization modes are present. We estimate the amplitudes of all the relevant modes for the peculiar as well as the regular white dwarfs. We further discuss the possibility of their detections through future-based gravitational wave detectors, such as LISA, ALIA, DECIGO, BBO, or Einstein Telescope, and thereby put constraints or rule out various modified theories of gravity. This exploration links the theory with possible observations through gravitational wave in f(R) gravity.

The spectra of active galactic nuclei usually exhibit wings in some emission lines, such as [OIII]$\lambda\lambda$5007,4959, with these wings generally being blueshifted and related to strong winds and outflows. The aim of this work was to analyse the [OIII] emission lines in broad line Seyfert 1 (BLS1) galaxies in order to detect the presence of wings, and to study the [OIII] line properties and their possible connection with the central engine. In addition, we attempted to compare the black hole mass distribution in both BLS1 galaxies with symmetric and blue-asymmetric [OIII] profiles. For this purpose, we carried out a spectroscopic study of a sample of 45 nearby southern BLS1 galaxies from the 6 Degree Field Galaxy survey. The [OIII] emission lines were well fitted using a single Gaussian function in 23 galaxies, while 22 objects presented a wing component and required a double-Gaussian decomposition. By computing the radial velocity difference between the wing and core centroids (i.e. $\Delta$v), we found 18 galaxies exhibiting blueshifted wings, 2 objects presenting red wings and 2 galaxies showing symmetric wings ($\Delta$v$= 0$). Moreover, $\Delta$v was slightly correlated with the black hole mass. In addition, we computed the radial velocity difference of the blue-side full extension of the wing relative to the centroid of the core component through the \emph{blue emission} parameter, which revealed a correlation with black hole mass, in agreement with previous results reported for narrow line galaxies. Finally, in our sample, similar black hole mass distributions were observed in both BLS1 galaxies with symmetric and blueshifted asymmetric [OIII] profiles.

H$\alpha$ blobs are off-galaxy emission-line regions with weak or no optical counterparts. They are mostly visible in H$\alpha$ line, appearing as concentrated blobs. Such unusual objects have been rarely observed and studied, and their physical origin is still unclear. We have identified 13 H$\alpha$ blobs in the public data of MaNGA survey, by visually inspecting both the SDSS optical images and the spatially resolved maps of H$\alpha$ line for $\sim 4600$ galaxy systems. Among the 13 H$\alpha$ blobs, 2 were reported in previously MaNGA-based studies and 11 are newly discovered. This sample, though still small in size, is by far the largest sample with both deep imaging and integral field spectroscopy. Therefore, for the first time we are able to perform statistical studies to investigate the physical origin of H$\alpha$ blobs. We examine the physical properties of these H$\alpha$ blobs and their associated galaxies, including their morphology, environments, gas-phase metallicity, kinematics of ionized gas, and ionizing sources. We find that the H$\alpha$ blobs in our sample can be broadly divided into two groups. One is associated with interacting/merging galaxy systems, of which the ionization is dominated by shocks or diffuse ionized gas. It is likely that these H$\alpha$ blobs used to be part of their nearby galaxies, but were stripped away at some point due to tidal interactions. The other group is found in gas-rich systems, appearing as low-metallicity star-forming regions that are visually detached from the main galaxy.

Mergers of black hole (BH) and neutron star (NS) binaries are of interest since the emission of gravitational waves (GWs) can be followed by an electromagnetic (EM) counterpart, which could power short gamma-ray bursts. Until now, LIGO/Virgo has only observed a candidate BH-NS event, GW190426\_152155, which was not followed by any EM counterpart. We show how the presence (absence) of a remnant disk, which powers the EM counterpart, can be used along with spin measurements by LIGO/Virgo to derive a lower (upper) limit on the radius of the NS. For the case of GW190426\_152155, large measurement errors on the spin prevent from placing an upper limit on the NS radius. Our proposed method works best when the aligned component of the BH spin (with respect to the orbital angular momentum) is the largest, and can be used to complement the information that can be extracted from the GW signal to derive valuable information on the NS equation of state.

We present Space-VLBI RadioAstron observations at 1.6 GHz and 4.8 GHz of the flat spectrum radio quasar 3C 273, with detections on baselines up to 4.5 and 3.3 Earth Diameters, respectively. Achieving the best angular resolution at 1.6 GHz to date, we have imaged limb-brightening in the jet, not previously detected in this source. In contrast, at 4.8 GHz, we detected emission from a central stream of plasma, with a spatial distribution complementary to the limb-brightened emission, indicating an origin in the spine of the jet. While a stratification across the jet width in the flow density, internal energy, magnetic field, or bulk flow velocity are usually invoked to explain the limb-brightening, the different jet structure detected at the two frequencies probably requires a stratification in the emitting electron energy distribution. Future dedicated numerical simulations will allow the determination of which combination of physical parameters are needed to reproduce the spine/sheath structure observed by RadioAstron in 3C 273.

The study of stripped-envelope core-collapse supernovae (SNe), with evidence for strong interaction of SN ejecta with the circumstellar medium (CSM), provides insights into the pre-supernova progenitor, and a fast-forwarded view of the progenitor mass-loss history. In this context, we present late-time radio observations of SN2004dk, a type Ibc supernova located in the galaxy, NGC 6118, at a distance of $d_L \approx 23$ Mpc. About 15 years after explosion, SN2004dk has shown evidence for H$\alpha$ emission, possibly linked to the SN ejecta interacting with an H-rich CSM. Using data from the VLA Low Band Ionosphere and Transient Experiment (VLITE), we confirm the presence of a late-time radio re-brightening accompanying the observed H$\alpha$ emission. We model the SN2004dk radio light curves within the (spherically symmetric) synchrotron-self-absorption (SSA) model. Within this model, our VLITE observations combined with previously collected VLA data favor an interpretation of SN2004dk as a strongly CSM-interacting radio SN going through a complex environment shaped by a non-steady mass-loss from the SN progenitor.

The correct evaluation of gradients is at the cornerstone of the smoothed particle hydrodynamics (SPH) technique. Using an integral approach to estimate gradients has proven to enhance accuracy substantially. Such approach retains the Lagrangian structure of SPH equations and is fully conservative. In this paper we study, among other things, the connection between the choice of the volume elements (VEs), which enters in the SPH summations, and the accuracy in the gradient estimation within the integral approach scheme (ISPH). A new kind of VEs are proposed which improve the partition of unit and are fully compatible with the Lagrangian formulation of SPH, including the grad-h corrections. Using analytic considerations, simple static toy models in 1D, and a few full 3D test cases, we show that any improvement in the partition of unit also leads to a better calculation of gradients when the integral approach is used jointly. Additionally, we propose a simple-to-implement variant of the ISPH scheme which is more adequate to handle sharp density contrasts.

The major programme for observing young, non-recycled pulsars with the Parkes telescope has transitioned from a narrow-band system to an ultra-wideband system capable of observing between 704 and 4032 MHz. We report here on the initial two years of observations with this receiver. Results include dispersion measure (DM) and Faraday rotation measure (RM) variability with time, determined with higher precision than hitherto, flux density measurements and the discovery of several nulling and mode changing pulsars. PSR J1703-4851 is shown to be one of a small subclass of pulsars that has a weak and a strong mode which alternate rapidly in time. PSR J1114-6100 has the fourth highest |RM| of any known pulsar despite its location far from the Galactic Centre. PSR J1825-1446 shows variations in both DM and RM likely due to its motion behind a foreground supernova remnant.

A critical question in astrobiology is whether exoEarth candidates (EECs) are Earth-like, in that they originate life that progressively oxygenates their atmospheres similarly to Earth. We propose answering this question statistically by searching for O2 and O3 on EECs with missions such as HabEx or LUVOIR. We explore the ability of these missions to constrain the fraction, fE, of EECs that are Earth-like in the event of a null detection of O2 or O3 on all observed EECs. We use the Planetary Spectrum Generator to simulate observations of EECs with O2 and O3 levels based on Earth's history. We consider four instrument designs: LUVOIR-A (15m), LUVOIR-B (8m), HabEx with a starshade (4m, "HabEx/SS"), HabEx without a starshade (4m, "HabEx/no-SS"); as well as three estimates of the occurrence rate of EECs (eta_earth): 24%, 5%, and 0.5%. In the case of a null-detection, we find that for eta_earth = 24%, LUVOIR-A, LUVOIR-B, and HabEx/SS would constrain fE to <= 0.094, <= 0.18, and <= 0.56, respectively. This also indicates that if fE is greater than these upper limits, we are likely to detect O3 on at least 1 EEC. Conversely, we find that HabEx/no-SS cannot constrain fE, due to the lack of an coronagraph ultraviolet channel. For eta_earth = 5%, only LUVOIR-A and LUVOIR-B would be able to constrain fE, to <= 0.45 and <= 0.85, respectively. For eta_earth = 0.5%, none of the missions would allow us to constrain fE, due to the low number of detectable EECs. We conclude that the ability to constrain fE is more robust to uncertainties in eta_earth for missions with larger aperture mirrors. However all missions are susceptible to an inconclusive null detection if eta_earth is sufficiently low.

We have observed the very low-mass Class 0 protostar IRAS 15398-3359 at scales ranging from 50 au to 1800 au, as part of the ALMA Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source along a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-like structure connected to the outer end of the linear feature and separated from the protostar, IRAS 15398-3359, by 1200 au. The arc-like structure is blue-shifted with respect to the systemic velocity. A velocity gradient of 1.2 km/s over 1200 au along the linear feature seen in the H2CO emission connects the protostar and the arc-like structure kinematically. SO, SiO, and CH3OH are known to trace shocks, and we interpret the arc-like structure as a relic shock region produced by an outflow previously launched by IRAS 15398-3359. The velocity gradient along the linear structure can be explained as relic outflow motion. The origins of the newly observed arc-like structure and extended linear feature are discussed in relation to turbulent motions within the protostellar core and episodic accretion events during the earliest stage of protostellar evolution.

In this work, we use the spectroscopy-based stellar color regression (SCR) method with ~ 0.7 million common stars between LAMOST DR7 and Gaia EDR3 to acquire color corrections in G - GRP and GBP - GRP. A sub-mmag precision is achieved. Our results demonstrate that improvements in the calibration process of the EDR3 have removed the color term in GBP - GRP and eliminated the discontinuity caused by the changes of instrument configurations to a great extent. However, modest systematic trends with G magnitude are still detected. The corresponding color correction terms as a function of G are provided for 9.5 < G < 17.5 and compared with other determinations. Possible applications as well as limitations of our results are discussed.

We present high-precision linear polarization observations of four bright hot Jupiter systems ($\tau$ Boo, HD 179949, HD 189733 and 51 Peg) and use the data to search for polarized reflected light from the planets. The data for 51 Peg are consistent with a reflected light polarization signal at about the level expected with 2.8$\sigma$ significance and a false alarm probability of 1.9 per cent. More data will be needed to confirm a detection of reflected light in this system. HD 189733 shows highly variable polarization that appears to be most likely the result of magnetic activity of the host star. This masks any polarization due to reflected light, but a polarization signal at the expected level of $\sim$20 ppm cannot be ruled out. $\tau$ Boo and HD 179949 show no evidence for polarization due to reflected light. The results are consistent with the idea that many hot Jupiters have low geometric albedos. Conclusive detection of polarized reflected light from hot Jupiters is likely to require further improvements in instrument sensitivity.

Major flares and coronal mass ejections (CMEs) tend to originate from the compact polarity inversion lines (PILs) in the solar active regions (ARs). Recently, a scenario named as "collisional shearing" is proposed by \citet{Chintzoglou_2019} to explain the phenomenon, which suggests that the collision between different emerging bipoles is able to form the compact PIL, driving the shearing and flux cancellation that are responsible to the subsequent large activities. In this work, through tracking the evolution of 19 emerging ARs from their birth until they produce the first major flares or CMEs, we investigated the source PILs of the activities, i.e., the active PILs, to explore the generality of "collisional shearing". We find that none of the active PILs is the self PIL (sPIL) of a single bipole. We further find that 11 eruptions originate from the collisional PILs (cPILs) formed due to the collision between different bipoles, 6 from the conjoined systems of sPIL and cPIL, and 2 from the conjoined systems of sPIL and ePIL (external PIL between the AR and the nearby preexisting polarities). Collision accompanied by shearing and flux cancellation is found developing at all PILs prior to the eruptions, with $84\%$ (16/19) cases having collisional length longer than 18~Mm. Moreover, we find that the magnitude of the flares is positively correlated with the collisional length of the active PILs, indicating that the intenser activities tend to originate from the PILs with severer collision. The results suggest that the "collisional shearing", i.e., bipole-bipole interaction during the flux emergence is a common process in driving the major activities in emerging ARs.

We report detection of PSR B0656$+$14 with the Gran Telescopio Canarias in narrow optical $F657$, $F754$, $F802$, and $F902$ and near-infrared $JHK_s$ bands. The pulsar detection in the $K_s$ band extends its spectrum to 2.2 $\mu$m and confirms its flux increase towards the infrared. We also present a thorough analysis of the optical spectrum obtained by us with the VLT. For a consistency check, we revised the pulsar near-infrared and narrow-band photometry obtained with the \textit{HST}. We find no narrow spectral lines in the optical spectrum. We compile available near-infrared-optical-UV and archival 0.3-20keV X-ray data and perform a self-consistent analysis of the rotation phase-integrated spectrum of the pulsar using unified spectral models. The spectrum is best fitted by the four-component model including two blackbodies, describing the thermal emission from the neutron star surface and its hot polar cap, the broken power-law, originating from the pulsar magnetosphere, and an absorption line near $\sim$0.5 keV detected previously. The fit provides better constraints on the model parameters than using only a single spectral domain. The derived surface temperature is $T_{NS}^{\infty}=7.9(3)\times10^5$K. The intrinsic radius (7.8-9.9 km) of the emitting region is smaller than a typical neutron star radius (13km) and suggests a nonuniform temperature distribution over the star surface. In contrast, the derived radius of the hot polar cap is about twice as large as the `canonical' one. The spectrum of the nonthermal emission steepens from the optical to X-rays and has a break near 0.1 keV. The X-ray data suggest the presence of another absorption line near 0.3keV.

We study galaxies undergoing ram pressure stripping in the Virgo cluster to examine whether we can identify any discernible trend in their star formation activity. We first use 48 galaxies undergoing different stages of stripping based on HI morphology, HI deficiency, and relative extent to the stellar disk, from the VIVA survey. We then employ a new scheme for galaxy classification which combines HI mass fractions and locations in projected phase space, resulting in a new sample of 365 galaxies. We utilize a variety of star formation tracers, which include g - r, WISE [3.4] - [12] colors, and starburstiness that are defined by stellar mass and star formation rates to compare the star formation activity of galaxies at different stripping stages. We find no clear evidence for enhancement in the integrated star formation activity of galaxies undergoing early to active stripping. We are instead able to capture the overall quenching of star formation activity with increasing degree of ram pressure stripping, in agreement with previous studies. Our results suggest that if there is any ram pressure stripping induced enhancement, it is at best locally modest, and galaxies undergoing enhancement make up a small fraction of the total sample. Our results also indicate that it is possible to trace galaxies at different stages of stripping with the combination of HI gas content and location in projected phase space, which can be extended to other galaxy clusters that lack high-resolution HI imaging.

One of the long-term goals of exoplanet science is the atmospheric characterization of dozens of small exoplanets in order to understand their diversity and search for habitable worlds and potential biosignatures. Achieving this goal requires a space mission of sufficient scale. We seek to quantify the exoplanet detection performance of a space-based mid-infrared nulling interferometer that measures the thermal emission of exoplanets. For this, we have developed an instrument simulator that considers all major astrophysical noise sources and coupled it with Monte Carlo simulations of a synthetic exoplanet population around main-sequence stars within 20 pc. This allows us to quantify the number (and types) of exoplanets that our mission concept could detect over a certain time period. Two different scenarios to distribute the observing time among the stellar targets are discussed and different apertures sizes and wavelength ranges are considered. Within a 2.5-year initial search phase, an interferometer consisting of four 2 m apertures covering a wavelength range between 4 and 18.5 $\mu$m could detect up to ~550 exoplanets with radii between 0.5 and 6 R$_\oplus$ with an integrated SNR$\ge$7. At least ~160 of the detected exoplanets have radii $\le$1.5 R$_\oplus$. Depending on the observing scenario, ~25-45 rocky exoplanets (objects with radii between 0.5 and 1.5 $_{\oplus}$) orbiting within the empirical habitable zone (eHZ) of their host stars are among the detections. With four times 3.5 m aperture size, the total number of detections can increase to up to ~770, including ~60-80 rocky, eHZ planets. With four times 1 m aperture size, the maximum detection yield is ~315 exoplanets, including $\le$20 rocky, eHZ planets. In terms of predicted detection yield, such a mission can compete with large single-aperture reflected light missions. (abridged)

The Large Area X-ray Proportional Counter (LAXPC) instrument on-board AstroSat has three nominally identical detectors for timing and spectral studies in the energy range of 3--80 keV. The performance of these detectors during the five years after the launch of AstroSat is described. Currently, only one of the detector is working nominally. The variation in pressure, energy resolution, gain and background with time are discussed. The capabilities and limitations of the instrument are described. A brief account of available analysis software is also provided.

We report the discovery of a bright, compact ultraviolet source at a projected separation of 1.1~kpc from the known active galactic nucleus (AGN) in Mrk~766 based on Astrosat/UVIT observations. We perform radial profile analysis and derive the UV flux almost free from the nearby contaminating sources. The new source is about 2.5 and 5.6 times fainter than the AGN in the far and near UV bands. The two sources appear as a pair of nuclei in Mrk~766. We investigate the nature of the new source based on the UV flux ratio, X-ray and optical emission. The new source is highly unlikely to be another accreting supermassive black hole in Mrk~766 as it lacks X-ray emission. We find that the UV/Optical flux of the new source measured at four different bands closely follow the shape of the template spectrum of starburst galaxies. This strongly suggests that the new source is a compact star-forming region.

We study the release of energy during the gradual phase of a flare, characterized by faint bursts of non-thermal hard X-ray (HXR) emission associated with decimetric radio spikes and type III radio bursts starting at high frequencies and extending to the heliosphere. We characterize the site of electron acceleration in the corona and study the radial evolution of radio source sizes in the high corona. Imaging and spectroscopy of the HXR emission with Fermi and RHESSI provide a diagnostic of the accelerated electrons in the corona as well as a lower limit on the height of the acceleration region. Radio observations in the decimetric range with the ORFEES spectrograph provide radio diagnostics close to the acceleration region. Radio spectro-imaging with LOFAR in the meter range provide the evolution of the radio source sizes with their distance from the Sun, in the high corona. Non-thermal HXR bursts and radio spikes are well correlated on short timescales. The spectral index of non-thermal HXR emitting electrons is -4 and their number is about $2\times 10^{33}$ electrons/s. The density of the acceleration region is constrained between $1-5 \times 10^9$ cm$^{-3}$. Electrons accelerated upward rapidly become unstable to Langmuir wave production, leading to high starting frequencies of the type III radio bursts, and the elongation of the radio beam at its source is between 0.5 and 11.4 Mm. The radio source sizes and their gradient observed with LOFAR are larger than the expected size and gradient of the size of the electron beam, assuming it follows the expansion of the magnetic flux tubes. These observations support the idea that the fragmentation of the radio emission into spikes is linked to the fragmentation of the acceleration process itself. The combination of HXR and radio diagnostics in the corona provides strong constrains on the site of electron acceleration.

A common assumption used in the study of accretion disks is that the magnetic energy density and the kinetic energy density should be in equipartition. This assumption relies on the faster growth rate of the magnetic field strength against the kinetic energy of the particles in the flow, for decreasing radius, combined with a dissipation mechanism that tends towards equipartition. In this paper, we examine this assumption by modeling the radio, mm and optical spectra of several black hole binaries in their quiescent state. We use a standard two-component disk model, consisting of an inner geometrically thick and optically thin disk, emitting thermal synchrotron radiation, along with an outer, thin disk, which radiates as a multicolor blackbody. We find that at the low accretion rates typical of the quiescent state, the spectral shape is qualitatively reproduced using magnetic fields that are between 0.1% and 1% of the equipartition value, considerably smaller than previously thought. We discuss our findings in view of (1) the launching of jets in these objects, which is commonly believed to rely on the presence of a strong magnetic field in the central region of the disk; and (2) the role of magnetic dissipation in the structure of the inflow.

BlueMUSE (Blue Multi Unit Spectroscopic Explorer) is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph proposed for the Very Large Telescope (VLT) and based on the MUSE concept. BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities in the 350 - 580 nm range: an optimised end-to-end transmission down to 350 nm, a larger FoV (up to $1.4 \times 1.4$ arcmin$^2$) sampled at 0.3 arcsec, and a higher spectral resolution ($\lambda/\Delta \lambda \sim 3500$) compared to MUSE. To our knowledge, achieving such capabilities with a comparable mechanical footprint and an identical detector format ($4\text{k} \times 4\text{k}$, 15 $\mathrm{\mu m}$ CCD) would not be possible with a conventional spectrograph design. In this paper, we present the optomechanical architecture and design of BlueMUSE at pre-phase A level, with a particular attention to some original aspects such as the use of curved detectors.

In recent years, the retrieval of entire asteroids has received significant attention, with many approaches leveraging the invariant manifolds of the Circular-Restricted Three-body Problem to capture an asteroid into a periodic orbit about the $L_1$ or $L_2$ points of the Sun-Earth system. Previous works defined an `Easily Retrievable Object' (ERO) as any Near-Earth Object (NEO) which is retrievable using these invariant manifolds with an impulsive $\Delta v$ of less than $500$ m/s. We extend the previous literature by analysing the Pareto fronts for the EROs discovered for the first time, using high-performance computing to lift optimisation constraints used in previous literature, and modifying the method used to filter unsuitable NEOs from the NEO catalogue. In doing so, we can demonstrate that EROs have approximately the same transfer cost for almost any possible transfer time, including single impulse transfers, which could offer significant flexibility to mission designers. We also identify $44$ EROs, of which $27$ are new, and improve on previously-known transfer solutions by up to $443$ m/s, including $17$ new capture trajectories with $\Delta v$ costs of less than $100$ m/s.

We investigate the possibility to induce double peaks of gravitational wave spectrum from primordial scalar perturbations in inflationary models with three inflection points. Where the inflection points can be generated from a polynomial potential or generated from Higgs like $\phi^4$ potential with the running of quartic coupling. In such models, the inflection point at large scales predicts the scalar spectral index and tensor-to-scalar ratio consistent with current CMB constraints, and the other two inflection points generate two large peaks in the scalar power spectrum at small scales, which can induce gravitational waves with double peaks energy spectrum. We calculate the induced gravitational wave energy spectrum and find that the double peaks signal can be detected in near future and can be distinguished with other single peak models.

Arrays of imaging atmospheric Cherenkov telescopes (IACT) are superb instruments to probe the very-high-energy gamma-ray sky. This type of telescope focuses the Cherenkov light emitted from air showers, initiated by very-high-energy gamma rays and cosmic rays, onto the camera plane. Then, a fast camera digitizes the longitudinal development of the air shower, recording its spatial, temporal, and calorimetric information. The properties of the primary very-high-energy particle initiating the air shower can then be inferred from those images: the primary particle can be classified as a gamma ray or a cosmic ray and its energy and incoming direction can be estimated. This so-called full-event reconstruction, crucial to the sensitivity of the array to gamma rays, can be assisted by machine learning techniques. We present a deep-learning driven, full-event reconstruction applied to simulated IACT events using CTLearn. CTLearn is a Python package that includes modules for loading and manipulating IACT data and for running deep learning models with TensorFlow, using pixel-wise camera data as input.

Water ice has a strong spectral feature at a wavelength of approximately $3~\mu$m, which plays a vital role in our understanding of the icy universe. In this study, we investigate the scattering polarization of this water-ice feature. The linear polarization degree of light scattered by $\mu$m-sized icy grains is known to be enhanced at the ice band; however, the dependence of this polarization enhancement on various grain properties is unclear. We find that the enhanced polarization at the ice band is sensitive to the presence of $\mu$m-sized grains as well as their ice abundance. We demonstrate that this enhancement is caused by the high absorbency of the water-ice feature, which attenuates internal scattering and renders the surface reflection dominant over internal scattering. Additionally, we compare our models with polarimetric observations of the low-mass protostar L1551 IRS 5. Our results show that scattering by a maximum grain radius of a few microns with a low water-ice abundance is consistent with observations. Thus, scattering polarization of the water-ice feature is a useful tool for characterizing ice properties in various astronomical environments.

We have performed Smoothed Particle Magneto-Hydrodynamics (SPMHD) calculations of colliding clouds to investigate the formation of massive stellar clusters, adopting a timestep criterion to prevent large divergence errors. We find that magnetic fields do not impede the formation of young massive clusters (YMCs), and the development of high star formation rates, although we do see a strong dependence of our results on the direction of the magnetic field. If the field is initially perpendicular to the collision, and sufficiently strong, we find that star formation is delayed, and the morphology of the resulting clusters is significantly altered. We relate this to the large amplification of the field with this initial orientation. We also see that filaments formed with this configuration are less dense. When the field is parallel to the collision, there is much less amplification of the field, dense filaments form, and the formation of clusters is similar to the purely hydrodynamical case. Our simulations reproduce the observed tendency for magnetic fields to be aligned perpendicularly to dense filaments, and parallel to low density filaments. Overall our results are in broad agreement with past work in this area using grid codes.

The European Open Science Cloud (EOSC) initiative faces the challenge of developing an agile, fit-for-purpose, and sustainable service-oriented platform that can address the evolving needs of scientific communities. The NEANIAS project plays an active role in the materialization of the EOSC ecosystem by actively contributing to the technological, procedural, strategic and business development of EOSC. We present the first outcomes of the NEANIAS activities relating to co-design and delivery of new innovative services for space research for data management and visualization (SPACE-VIS), map making and mosaicing (SPACE-MOS) and pattern and structure detection (SPACE-ML). We include a summary of collected user requirements driving our services and methodology for their delivery, together with service access details and pointers to future works.

VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high signal to noise, medium resolution spectra of galaxies at redshift between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We give a detailed description of sample selection, observations and data reduction procedures. The final catalogue reaches a target selection completeness of 40% at iAB = 25. The high Signal to Noise ratio of the spectra (above 7 in 80% of the spectra) and the dispersion of 2.5{\AA} allowed us to measure redshifts with high precision, the redshift measurement success rate reaching almost 100%. Together with the redshift catalogue and the reduced spectra, we also provide optical mid-IR photometry and physical parameters derived through SED fitting. The observed galaxy sample comprises both passive and star forming galaxies covering a stellar mass range 8.3< Log(M*/Msolar)<11.7. All catalogues and spectra are accessible through the survey database (this http URL) where all information can be queried interactively, and via the ESO Archive (https://www.eso.org/qi/).

Context. Prominence eruptions provide key observations to understand the launch of coronal mass ejections as their cold plasma traces a part of the unstable magnetic configuration. Aims. We select a well observed case to derive observational constraints for eruption models. Methods. We analyze the prominence eruption and loop expansion and contraction observed on 02 March 2015 associated with a GOES M3.7 class flare (SOL2015-03-02T15:27) using the data from Atmospheric Imaging Assembly (AIA) and the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). We study the prominence eruption and the evolution of loops using the time-distance techniques. Results. The source region is a decaying bipolar active region where magnetic flux cancellation is present for several days before the eruption. AIA observations locate the erupting prominence within a flux rope viewed along its local axis direction. We identify and quantify the motion of loops in contraction and expansion located on the side of the erupting flux rope. Finally, RHESSI hard X-ray observations identify the loop top and two foot-point sources. Conclusions. Both AIA and RHESSI observations support the standard model of eruptive flares. The contraction occurs 19 minutes after the start of the prominence eruption indicating that this contraction is not associated with the eruption driver. Rather, this prominence eruption is compatible with an unstable flux rope where the contraction and expansion of the lateral loop is the consequence of a side vortex developing after the flux rope is launched.

The differential rotation of the Sun is a crucial ingredient of the dynamo theory responsible for the generation of its magnetic field. Currently, the rotation profile of a star that hosts one or more transiting planets can be estimated. By detecting the same spot in a later transit, it is possible to infer the stellar rotation period at that latitude. In this work, we apply for the first time transit spot mapping to determine the differential rotation of Kepler-411, a K2V-type star with an average rotation period of 10.52 days, radius of 0.79 R$_\odot$ and mass of 0.83 M$_\odot$. Kepler-411 hosts at least four planets, the inner planet is a super-Earth with a radius of 1.88 R$_\oplus$ and an orbital period of 3.0051 days, whereas the two larger transiting planets are mini Neptunes with radii of 3.27 and 3.31 R$_\oplus$, and periods of 7.834435 and 58.0204 days, respectively. Their orbits are such that they transit the star at latitudes of -11$^{\circ}$, -21$^{\circ}$, and -49$^{\circ}$. Analysis of the transit light curves of the three planets resulted in the detection of a total of 198 spots. For each transit latitude, the rotation period of the star was estimated and the differential rotation pattern estimated independently. Then a solar like differential rotation profile was fit to the three rotation periods at the distinct latitudes, the result agreed extremely well with the previous ones, resulting in a differential shear of $0.0500\pm0.0006$ rd/d or a relative differential rotation of $8.4\pm0.1$\%.

We study the evolution of the binary black hole (BBH) mass distribution across cosmic time. The second gravitational-wave transient catalog (GWTC-2) from LIGO/Virgo contains BBH events out to redshifts $z \sim 1$, with component masses in the range $\sim5$--$80\,M_\odot$. In this catalog, the biggest black holes, with $m_1 \gtrsim 45\,M_\odot$, are only found at the highest redshifts, $z \gtrsim 0.4$. We ask whether the absence of high-mass BBH observations at low redshift indicates that the astrophysical BBH mass distribution evolves: the biggest BBHs only merge at high redshift, and cease merging at low redshift. Alternatively, this feature might be explained by gravitational-wave selection effects. Modeling the BBH primary mass spectrum as a power law with a sharp maximum mass cutoff (Truncated model), we find that the cutoff increases with redshift ($> 99.9\%$ credibility). An abrupt cutoff in the mass spectrum is expected from (pulsational) pair instability supernova simulations; however, GWTC-2 is only consistent with a Truncated mass model if the location of the cutoff increases from $45^{+13}_{-5}\,M_\odot$ at $z < 0.4$ to $80^{+16}_{-13}\,M_\odot$ at $z > 0.4$. Alternatively, if the primary mass spectrum has a break in the power law (Broken power law) at ${38^{+15}_{-8}\,M_\odot}$, rather than a sharp cutoff, the data are consistent with a non-evolving mass distribution. In this case, the overall rate of mergers, at all masses, increases with increasing redshift. Future observations will confidently distinguish between a sharp maximum mass cutoff that evolves with redshift and a non-evolving mass distribution with a gradual taper, such as a Broken power law. After $\sim 100$ BBH merger observations, a continued absence of high-mass, low-redshift events would provide a clear signature that the mass distribution evolves with redshift.

Clusters of galaxies can potentially produce cosmic rays (CRs) up to very-high energies via large-scale shocks and turbulent acceleration. Due to their unique magnetic-field configuration, CRs with energy $\leq 10^{17}$ eV can be trapped within these structures over cosmological time scales, and generate secondary particles, including neutrinos and gamma rays, through interactions with the background gas and photons. In this work, we compute the contribution from clusters of galaxies to the diffuse neutrino background. We employ three-dimensional cosmological magnetohydrodynamical simulations of structure formation to model the turbulent intergalactic medium. We use the distribution of clusters within this cosmological volume to extract the properties of this population, including mass, magnetic field, temperature, and density. We propagate CRs into this environment using multi-dimensional Monte Carlo simulations across different redshifts (from $z \sim 5$ to $z =0$), considering all relevant photohadronic, photonuclear, and hadronuclear interaction processes. For CRs injected with a power-law energy spectrum $\propto E^{-2}$, we find that clusters contribute to the diffuse flux observed by the IceCube Neutrino Observatory with up to $\sim 20 \;\%$. On the other hand, for a CR spectrum $\propto E^{-1.5}$, our results are compatible with the observed flux of IceCube, though most of the contribution comes from clusters with $M \gtrsim 10^{14} M_{\odot}$ and redshift $ z \lesssim 0.3$. These results suggest that clusters of galaxies can account for a sizeable fraction of the diffuse neutrino background.

The in-orbit results and lessons learned of the first Finnish satellite Aalto-1 are briefly presented in this paper. Aalto-1, a three-unit CubeSat which was launched in June 2017, performed AaSI (Aalto Spectral Imager), Radiation Monitor (RADMON), and Electrostatic Plasma Brake (EPB) missions. The satellite partly fulfilled its mission objectives and allowed to either perform or attempt the experiments. Although attitude control was partially functional, AaSI and RADMON were able to acquire valuable measurements. EPB was successfully commissioned but the tether deployment was not successful.

In this work we study a system of two galaxies, Astarte and Adonis, at z $\sim $2 when the Universe was undergoing its peak of star formation activity. Astarte is a dusty star-forming galaxy at the massive-end of the main sequence (MS) and Adonis is a less-massive, bright in ultraviolet (UV), companion galaxy with an optical spectroscopic redshift. We analyse the physical properties of this system, and probe the gas mass of Astarte with its ALMA CO emission, to investigate whether this ultra-massive galaxy is quenching or not. We use CIGALE - a spectral energy distribution modeling code - to derive the key physical properties of Astarte and Adonis, mainly their star formation rates (SFRs), stellar masses, and dust luminosities. We inspect the variation of the physical parameters depending on the assumed dust attenuation law. We also estimate the molecular gas mass of Astarte from its CO emission, using different $\alpha_{CO}$ and transition ratios ($r_{31}$) and discuss the implication of the various assumptions on the gas mass derivation. We find that Astarte exhibits a MS-like star formation activity, while Adonis is undergoing a strong starburst (SB) phase. The molecular gas mass of Astarte is far below the gas fraction of typical star-forming galaxies at z=2. This low gas content and high SFR, result in a depletion time of $0.22\pm0.07$ Gyrs, slightly shorter than what is expected for a MS galaxy at this redshift. The CO luminosity versus the total IR luminosity suggests a MS-like activity assuming a galactic conversion factor and a low transition ratio. The SFR of Astarte is of the same order using different attenuation laws, unlike its stellar mass that increases using shallow attenuation laws. We discuss these properties and suggest that Astarte might be experiencing a recent decrease of star formation activity and is quenching through the MS following a SB epoch.

Determining the black hole masses in active galactic nuclei (AGN) is of crucial importance to constrain the basic characteristics of their central engines and shed light on their growth and co-evolution with their host galaxies. While the black hole mass (MBH) can be robustly measured with dynamical methods in bright type 1 AGN, where the variable primary emission and the broad line region (BLR) are directly observed, a direct measurement is considerably more challenging if not impossible for the vast majority of heavily obscured type 2 AGN. In this work, we tested the validity of an X-ray-based scaling method to constrain the MBH in heavily absorbed AGN. To this end, we utilized a sample of type 2 AGN with good-quality hard X-ray data obtained by the nuSTAR satellite and with MBH dynamically constrained from megamaser measurements. Our results indicate that, when the X-ray broadband spectra are fitted with physically motivated self-consistent models that properly account for absorption, scattering, and emission line contributions from the putative torus and constrain the primary X-ray emission, then the X-ray scaling method yields MBH values that are consistent with those determined from megamaser measurements within their respective uncertainties. With this method we can therefore systematically determine the MBH in any type 2 AGN, provided that they possess good-quality X-ray data and accrete at a moderate to high rate.

A detailed study of comets active at large heliocentric distances (greater than 4 au) which enter the Solar System for the first time and are composed of matter in its elementary, unprocessed state, would help in our understanding of the history and evolution of the Solar System. In particular, contemporary giant planet formation models require the presence of accretion of volatile elements such as neon, argon, krypton, xenon and others, which initially could not survive at the distances where giant planets were formed. Nevertheless, the volatile components could be effectively delivered by the Kuiper-belt and Oort-cloud bodies, which were formed at temperatures below 30 K. This review is dedicated to the results of a multi-year comprehensive study of small bodies of the Solar System showing a comet-like activity at large heliocentric distances. The data were obtained from observations with the 6-meter telescope of SAO RAS equipped with multi-mode focal reducers SCORPIO and SCORPIO-2.

The Ingot wavefront sensor is a novel pupil-plane wavefront sensor, specifically designed to cope with the elongation typical of the extended nature of the Laser Guide Star (LGS). In the framework of the ELT, we propose an optical solution suitable for a Laser launch telescope, located outside the telescope pupil. In this paper, we present the current optical design, based on a reflective roof-shaped prism, which, at the level of the focal plane, splits the light from an LGS producing three beams. The three images of the telescope pupils can be then used for the retrieval of the first derivative of the wavefront. The 3D nature of such a device requires new alignment techniques to be determined theoretically and verified in the real world. A possible fully automated procedure, relying solely on the illumination observed at the three pupils, to align the prism to the image of the LGS is discussed. Careful attention needs to be put both on the telecentricity of the system and on the reference systems of the Ingot adjustments in the 3D space. This is crucial in order to disentangle all the possible misalignment effects. In this context, we devised a test-bench able to reproduce, in a scaled manner, the 3D illumination that the Ingot will face at the ELT, in order to validate the design and to perform preliminary tests of phase retrieval.

Bayesian inference is used to extract unknown parameters from gravitational wave signals. Detector noise is typically modelled as stationary, although data from the LIGO and Virgo detectors is not stationary. We demonstrate that the posterior of estimated waveform parameters is no longer valid under the assumption of stationarity. We show that while the posterior is unbiased, the errors will be under- or overestimated compared to the true posterior. A formalism was developed to measure the effect of the mismodelling, and found the effect of any form of non-stationarity has an effect on the results, but are not significant in certain circumstances. We demonstrate the effect of short-duration Gaussian noise bursts and persistent oscillatory modulation of the noise on binary-black-hole-like signals. In the case of short signals, non-stationarity in the data does not have a large effect on the parameter estimation, but the errors from non-stationary data containing signals lasting tens of seconds or longer will be several times worse than if the noise was stationary. Accounting for this limiting factor in parameter sensitivity could be very important for achieving accurate astronomical results, including an estimation of the Hubble parameter. This methodology for handling the non-stationarity will also be invaluable for analysis of waveforms that last minutes or longer, such as those we expect to see with the Einstein Telescope.

The ice-rich dwarf planet Ceres is the largest object in the main asteroid belt and is thought to have a brine or mud layer at a depth of tens of kilometers. Furthermore, recent surface deposits of brine-sourced material imply shallow feeder structures such as sills or dikes. Inductive sounding of Ceres can be performed using the solar wind as a source, as was done for the Moon during Apollo. However, the magnetotelluric method -- measuring both electric and magnetic fields at the surface -- is not sensitive to plasma effects that were experienced for Apollo, which used an orbit-to-surface magnetic transfer function. The highly conductive brine targets are readily separable from the resistive ice and rock interior, such that the depth to deep and shallow brines can be assessed simultaneously. The instrumentation will be tested on the Moon in 2023 and is ready for implementation on a Ceres landed mission.

There is a growing set of observational data demonstrating that cosmic rays exhibit small-scale anisotropies (5-30 deg) with amplitude deviations lying between 0.01-0.1 percent that of the average cosmic ray flux. A broad range of models have been proposed to explain these anisotropies ranging from finite-scale magnetic field structures to dark matter annihilation. The standard diffusion transport methods used in cosmic ray propagation do not capture the transport physics in a medium with finite-scale or coherent magnetic field structures. Here, we present a Monte Carlo transport method, applying it to a series of finite-scale magnetic field structures to determine the requirements of such fields in explaining the observed cosmic ray,small-scale anisotropies.

Ringed structures have been observed in a variety of protoplanetary discs. Among the processes that might be able to generate such features, the Secular Gravitational Instability (SGI) is a possible candidate. It has also been proposed that the SGI might lead to the formation of planetesimals during the non-linear phase of the instability. In this context, we employ two-fluid hydrodynamical simulations with self-gravity to study the non-axisymmetric, non-linear evolution of ringed perturbations that grow under the action of the SGI. We find that the non-linear evolution outcome of the SGI depends mainly on the initial linear growth rate. For SGI growth rates smaller than typically $\sigma \gtrsim 10^{-4}-10^{-5}\Omega$, dissipation resulting from dust feedback introduces a $m=1$ spiral wave in the gas, even for Toomre gas stability parameters $Q_g>2$ for which non-axisymmetric instabilities appear in a purely gaseous disc. This one-armed spiral subsequently traps dust particles until a dust-to-gas ratio $\epsilon \sim 1$ is achieved. For higher linear growth rates, the dust ring is found to undergo gravitational collapse until the bump in the surface density profile becomes strong enough to trigger the formation of dusty vortices through the Rossby Wave Instability (RWI). Enhancements in dust density resulting from this process are found to scale with the linear growth rate, and can be such that the dust density is higher than the Roche density, leading to the formation of bound clumps. Fragmentation of axisymmetric rings produced by the SGI might therefore appear as a possible process for the formation of planetesimals.

The H.E.S.S. Galactic Plane Survey (HGPS) revealed 78 TeV sources among which 47 are not clearly associated with a known object. We present a multiwavelength approach to constrain the origin of the emission from unidentified HGPS sources. We present a generic pipeline that explores a large database of multiwavelength archival data toward any region in the Galactic plane. Along with a visual inspection of the retrieved multiwavelength observations to search for faint and uncataloged counterparts, we derive a radio spectral index that helps disentangle thermal from nonthermal emission and a mean magnetic field through X-ray and TeV data in case of a leptonic scenario. We also search for a spectral connection between the GeV and the TeV regimes with the Fermi-LAT cataloged sources that may be associated with the unidentified HGPS source. We complete the association procedure with catalogs of known objects and with the source catalogs from instruments whose data are retrieved. The method is applied on two unidentified sources, namely HESS J1427$-$608 and HESS J1458$-$608, for which the multiwavelength constraints favor the pulsar wind nebula (PWN) scenario. We model their broadband nonthermal spectra in a leptonic scenario with a magnetic field $B \lesssim 10$ $\mu$G, which is consistent with that obtained from ancient PWNe. We place both sources within the context of the TeV PWN population to estimate the spin-down power and the characteristic age of the putative pulsar. We also shed light on two possibly significant $\gamma$-ray excesses in the HGPS: the first is located in the south of the unidentified source HESS J1632$-$478 and the second is spatially coincident with the synchrotron-emitting supernova remnant G28.6$-$0.1. The multiwavelength counterparts found toward both $\gamma$-ray excesses make these promising candidates for being new very-high energy $\gamma$-ray sources.

The indirect evidence for at least a dozen massive white dwarfs violating the Chandrasekhar mass-limit is considered to be one of the wonderful discoveries in astronomy for more than a decade. Researchers have already proposed a diverse amount of models to explain this astounding phenomenon. However, each of these models always carries some drawbacks. On the other hand, noncommutative geometry is one of the best replicas of quantum gravity, which is yet to be proved from observations. Madore introduced the idea of a fuzzy sphere to describe a formalism of noncommutative geometry. This article shows that the idea of a squashed fuzzy sphere can self-consistently explain the super-Chandrasekhar limiting mass white dwarfs. We further show that the length-scale beyond which the noncommutativity is prominent is an emergent phenomenon, and there is no prerequisite for an ad-hoc length-scale.

Black-hole spacetimes are known to possess closed light rings. We here present a remarkably compact theorem which reveals the physically intriguing fact that these unique null circular geodesics provide the {\it fastest} way, as measured by asymptotic observers, to circle around spinning Kerr black holes.

For scalar theories accommodating spherically symmetric Q-balls, there are also towers of quasi-stable composite Q-balls, called charge swapping Q-balls (CSQs). We investigate the properties, particularly the lifetimes, of these long-lived CSQs in 2+1D and 3+1D using numerical simulations with efficient second order absorbing boundary conditions. We find that the evolution of a CSQ typically consists of 4 distinct stages: initial relaxation, first plateau (CSQ stage), fast decay and second plateau (oscillon stage). We chart the lifetimes of CSQs for different parameters of the initial conditions and of the potential, and show the attractor behavior and other properties of the CSQs.

There are many uses for linear fitting; the context here is interpolation and denoising of data, as when you have calibration data and you want to fit a smooth, flexible function to those data. Or you want to fit a flexible function to de-trend a time series or normalize a spectrum. In these contexts, investigators often choose a polynomial basis, or a Fourier basis, or wavelets, or something equally general. They also choose an order, or number of basis functions to fit, and (often) some kind of regularization. We discuss how this basis-function fitting is done, with ordinary least squares and extensions thereof. We emphasize that it is often valuable to choose far more parameters than data points, despite folk rules to the contrary: Suitably regularized models with enormous numbers of parameters generalize well and make good predictions for held-out data; over-fitting is not (mainly) a problem of having too many parameters. It is even possible to take the limit of infinite parameters, at which, if the basis and regularization are chosen correctly, the least-squares fit becomes the mean of a Gaussian process. We recommend cross-validation as a good empirical method for model selection (for example, setting the number of parameters and the form of the regularization), and jackknife resampling as a good empirical method for estimating the uncertainties of the predictions made by the model. We also give advice for building stable computational implementations.

We derive the non-relativistic limit of a massive vector field. We show that the Cartesian spatial components of the vector behave as three identical, non-interacting scalar fields. We find classes of spherical, cylindrical, and planar self-gravitating vector solitons in the Newtonian limit. The gravitational properties of the lowest-energy vector solitons$\mathrm{-}$the gravitational potential and density field$\mathrm{-}$depend only on the net mass of the soliton and the vector particle mass. In particular, these self-gravitating, ground-state vector solitons are independent of the distribution of energy across the vector field components, and are indistinguishable from their scalar-field counterparts. Fuzzy Vector Dark Matter models can therefore give rise to halo cores with identical observational properties to the ones in scalar Fuzzy Dark Matter models. We also provide novel hedgehog vector soliton solutions, which cannot be observed in scalar-field theories. The gravitational binding of the lowest-energy hedgehog halo is about three times weaker than the ground-state vector soliton. Finally, we show that no spherically symmetric solitons exist with a divergence-free vector field.

We propose a combined mechanism to realize both winding inflation and de Sitter uplifts. We realize the necessary structure of competing terms in the scalar potential not via tuning the vacuum expectation values of the complex structure moduli, but by a hierarchy of the Gopakumar-Vafa invariants of the underlying Calabi-Yau threefold. To show that Calabi-Yau threefolds with the prescribed hierarchy actually exist, we explicitly create a database of all the genus $0$ Gopakumar-Vafa invariants up to total degree $10$ for all the complete intersection Calabi-Yau's up to Picard number $9$. As a side product, we also identify all the redundancies present in the CICY list, up to Picard number $13$. Both databases can be accessed at this link: https://www.desy.de/~westphal/GV_CICY_webpage/GVInvariants.html .

The Earth's magnetosphere represents a natural plasma laboratory that allows us to study the behavior of particle distribution functions in the absence of Coulomb collisions, typically described by the Kappa distributions. We have investigated the properties of these functions for ions and electrons in different magnetospheric regions, thereby making it possible to reveal the $\kappa$-parameters for a wide range of plasma beta ($\beta$) values (from $10^{-3}$ to $10^{2}$). This was done using simultaneous ion and electron measurements from the five Time History of Events and Macroscale Interactions during Substorms (THEMIS) spacecraft spanning the years 2008 to 2018. It was found that for a fixed plasma $\beta$, the $\kappa$-index and core energy ($E_c$) of the distribution can be modeled by the power-law $\kappa=AE_c^\gamma$ for both species, and the relation between $\beta$, $\kappa$, and $E_c$ is much more complex than earlier reported: both $A$ and $\gamma$ exhibit systematic dependencies with $\beta$. Our results indicate that $\beta \sim 0.1-0.3$ is a range where the plasma is more dynamic since it is influenced by both the magnetic field and temperature fluctuations, which suggests that the transition between magnetically dominated plasmas to kinetically dominated plasmas occurs at these values of $\beta$. For $\beta > 1 $, both $A$ and $\gamma$ take nearly constant values, a feature that is especially notable for the electrons and might be related to their demagnetization. The relation between $\beta$, $\kappa$, and $E_c$ that we present is an important result that can be used by theoretical models in the future.

The measured cosmological constant $\Lambda$ is usually interpreted as Dark Energy (DE) or modified gravity (MG). Here we propose instead that $\Lambda$ corresponds to a boundary term in the action of classical General Relativity. The action is zero for a perfect fluid solution and this fixes $\Lambda$ to the average density $\rho$ and pressure $p$ inside a primordial causal boundary set by inflation: $\Lambda = 4\pi G <\rho+3p>$. This explains both why the observed value of $\Lambda$ is related to the matter density today and also why other contributions to $\Lambda$, such as DE or MG, do not produce cosmic expansion. Cosmic acceleration results from the repulsive boundary force that occurs when the expansion reaches the causal horizon. This universe is similar to the $\Lambda$CDM universe, except on the largest observable scales, where we expect departures from homogeneity/isotropy, such as variations in cosmological parameters indicated by recent observations.

Image-to-image translation with Deep Learning neural networks, particularly with Generative Adversarial Networks (GANs), is one of the most powerful methods for simulating astronomical images. However, current work is limited to utilizing paired images with supervised translation, and there has been rare discussion on reconstructing noise background that encodes instrumental and observational effects. These limitations might be harmful for subsequent scientific applications in astrophysics. Therefore, we aim to develop methods for using unpaired images and preserving noise characteristics in image translation. In this work, we propose a two-way image translation model using GANs that exploits both paired and unpaired images in a semi-supervised manner, and introduce a noise emulating module that is able to learn and reconstruct noise characterized by high-frequency features. By experimenting on multi-band galaxy images from the Sloan Digital Sky Survey (SDSS) and the Canada France Hawaii Telescope Legacy Survey (CFHT), we show that our method recovers global and local properties effectively and outperforms benchmark image translation models. To our best knowledge, this work is the first attempt to apply semi-supervised methods and noise reconstruction techniques in astrophysical studies.

Constraints on the Covariant Canonical Gauge Gravity (CCGG) theory from low-redshift cosmology are studied. The formulation extends Einstein's theory of General Relativity (GR) by a quadratic Riemann-Cartan term in the Lagrangian, controlled by a "deformation" parameter. In the Friedman universe this leads to an additional geometrical stress energy and promotes, due to the necessary presence of torsion, the cosmological constant to a time-dependent function. The MCMC analysis of the combined data sets of Type Ia Supernovae, Cosmic Chronometers and Baryon Acoustic Oscillations yields a fit that is well comparable with the $\Lambda$CDM results. The modifications implied in the CCGG approach turn out to be subdominant in the low-redshift cosmology. However, a non-zero spatial curvature and deformation parameter are shown to be consistent with observations.

We revisit gravitational particle production from the Stokes phenomenon viewpoint, which is crucial to understand asymptotic behavior of mode functions with a time dependent frequency. One of our purposes of this work is to show how the analysis focusing on the Stokes phenomenon can be used to analytically estimate non-perturbative particle production rate. In particular, with several examples of time-dependent background, we examine some methods that make the analysis more practical. Specifically, we consider the particle production in simple expanding backgrounds, $R^2$ inflation, and a model with smoothly changing mass. Since some of our models show difficulties in analyzing Stokes phenomenon, we discuss simplification of the problem and the accuracy of analytic estimation within our approximations. We also propose an approximation to take into account the most important contribution among infinite number of turning points and poles, which greatly simplifies the problem, but still gives a good analytic estimation.

In the path integral formulation of the reduced phase space Loop Quantum Gravity (LQG), we propose a new approach to allow the spatial cubic lattice (graph) to change dynamically in the physical time evolution. The equations of motion of the path integral derive the effective dynamics of cosmology from the full LQG, when we focus on solutions with homogeneous and isotropic symmetry. The resulting cosmological effective dynamics with the dynamical lattice improves the effective dynamics obtained earlier from the path integral with fixed spatial lattice: The improved effective dynamics recovers the FLRW cosmology at low energy density and resolves the big-bang singularity with a bounce. The critical density $\rho_c$ at the bounce is Planckian $\rho_c\sim \Delta^{-1}$ where $\Delta$ is a Planckian area serving as certain UV cut-off of the effective theory. The effective dynamics gives the unsymmetric bounce and has the de-Sitter (dS) spacetime in the past of the bounce. The cosmological constant $\Lambda_{eff}$ of the dS spacetime is emergent from the quantum effect $\Lambda_{eff}\sim\Delta^{-1}$. These results are qualitatively similar to the properties of $\bar{\mu}$-scheme Loop Quantum Cosmology (LQC). Moreover, we generalize the earlier path integral formulation of the full LQG by taking into account the coupling with an additional real scalar field, which drives the slow-roll inflation of the effective cosmological dynamics. In addition, we discuss the cosmological perturbation theory on the dynamical lattice, and the relation to the Mukhanov-Sasaki equation.

Important efforts are currently made for understanding the so-called kinetic instabilities, driven by the anisotropy of different species of plasma particles present in the solar wind and terrestrial magnetosphere. These instabilities are fast enough to efficiently convert the free energy of plasma particles into enhanced (small-scale) fluctuations with multiple implications, regulating the anisotropy of plasma particles. In this paper we use both linear and quasilinear (QL) frameworks to describe complex unstable regimes, which realistically combine different temperature anisotropies of electrons and ions (protons). Thus parameterized are various instabilities, e.g., proton and electron firehose, electromagnetic ion cyclotron, and whistler instability, showing that their main linear properties are markedly altered by the interplay of anisotropic electrons and protons. Linear theory may predict a strong competition of two instabilities of different nature when their growth rates are comparable. In the QL phase wave fluctuations grow and saturate at different levels and temporal scales, by comparison to the individual excitation of the proton or electron instabilities. In addition, cumulative effects of the combined proton and electron induced fluctuations can markedly stimulate the relaxations of their temperature anisotropies. Only whistler fluctuations inhibit the efficiency of proton firehose fluctuations in the relaxation of anisotropic protons. These results offer valuable premises for further investigations in numerical simulations, to decode the full spectrum of kinetic instabilities resulting from the interplay of anisotropic electrons and protons in space plasmas.

Lower atmospheric global dust storms affect the small- and large-scale weather and variability of the whole Martian atmosphere. Analysis of the CO$_2$ density data from the Neutral Gas and Ion Mass Spectrometer instrument (NGIMS) on board NASA's Mars Atmosphere Volatile EvolutioN (MAVEN) spacecraft show a remarkable increase of GW-induced density fluctuations in the thermosphere during the 2018 major dust storm with distinct latitude and local time variability. The mean thermospheric GW activity increases by a factor of two during the storm event. The magnitude of relative density perturbations is around 20% on average and 40% locally. One and a half months later, the GW activity gradually decreases. Enhanced temperature disturbances in the Martian thermosphere can facilitate atmospheric escape. For the first time, we estimate that, for a 20% and 40% GW-induced disturbances, the net increase of Jeans escape flux of hydrogen is a factor of 1.3 and 2, respectively.

The flattening of galaxy rotation curves in the weak gravity regime has been the center of scientific debate for decades. The cold dark matter ($\rm CDM$) paradigm has been posited as the standard explanation, where the mass of the dark halo $m(r) = r/Gn$ roughly increases linearly with radial distance at large distances where $G$ the gravitational constant and $n$ a dimensionless parameter which depends on the amount of baryonic matter $M$ within the galaxy. Despite numerous advances in modeling galaxy formation and evolution within the dark energy-cold dark matter ($\Lambda{\rm CDM}$) model in cosmology, a scientific consensus on the origin of the observed dependence of the dimensionless parameter $n = (GMa_{0})^{-1/2}$ on the mass of baryonic matter $M$ within the galaxy (the Tully-Fisher relation), and the connection of the cosmological constant $\Lambda$ to the acceleration parameter $a_{0} \sim (\Lambda/3)^{1/2}$ remains elusive. Here, we propose relativistic equations of gravity ($\nabla_{\nu}\mathcal{K}^{\nu}_{\,\,\mu} = 8\pi GM\Psi^{*}\mathcal{D}_{\mu}\Psi$, where $\mathcal{K}_{\mu\nu}$ is a Hermitian tensor, $\mathcal{D}_{\mu}$ a covariant derivative and $\Psi$ is a complex-valued function) which we show to not only contain Einstein Field Equations but also satisfy the Tully-Fisher relation. In the weak field limit, the gravity equations reduce to a theory of $n$ bosons (Ginzburg Landau theory) where the order parameter $\Psi$ is normalized as $\int_{0}^{1/a_{0}} dr\,4\pi r^2\Psi^*\Psi = n$ and $1/a_{0} \sim (\Lambda/3)^{-1/2}$ is the cut-off radius comparable to the size of the de Sitter universe. Our investigations have significant implications for the dark matter versus MOdified Newtonian Dynamics (MOND) debate, since we provide a framework where the Tully-Fisher relation in galaxies is satisfied within the context of Einstein's general relativity.

Developing cyberinfrastructure for the growing needs of multi-messenger astrophysics requires expertise in both software development and domain science. However, due to the nature of scientific software development, many scientists neglect best practices for software engineering which results in software that is difficult to maintain. We present here a mitigation strategy where scientists adopt software development best practices by collaborating with professional software developers. Such a partnership brings inherent challenges. For the scientists, this can be a dependence on external resources and lack of control in the development process. For developers, this can be a reduction in effort available for core, non-scientific development goals. These issues can be alleviated by structuring the partnership using established software development practices, such as the Agile Scrum framework. This paper presents a case study wherein a scientist user group, the SuperNova Early Warning System (SNEWS), collaborated with a group of scientific software developers, the Scalable Cyberinfrastructure for Multi-Messenger Astrophysics (SCiMMA) project. The two organizations utilized an Agile Scrum framework to address the needs of each organization, mitigate the concerns of collaboration, and avoid pitfalls common to scientific software development. In the end, the scientists profited from a successful prototype and the software developers benefited from enhanced cyberinfrastructure and improved development skills. This suggests that structured collaborations could help address the prevailing difficulties in scientific computing.