GREX-PLUS (Galaxy Reionization EXplorer and PLanetary Universe Spectrometer) is a mission candidate for a JAXA's strategic L-class mission to be launched in the 2030s. Its primary sciences are two-fold: galaxy formation and evolution and planetary system formation and evolution. The GREX-PLUS spacecraft will carry a 1.2 m primary mirror aperture telescope cooled down to 50 K. The two science instruments will be onboard: a wide-field camera in the 2-8 $\mu$m wavelength band and a high resolution spectrometer with a wavelength resolution of 30,000 in the 10-18 $\mu$m band. The GREX-PLUS wide-field camera aims to detect the first generation of galaxies at redshift $z>15$. The GREX-PLUS high resolution spectrometer aims to identify the location of the water ``snow line'' in proto-planetary disks. Both instruments will provide unique data sets for a broad range of scientific topics including galaxy mass assembly, origin of supermassive blackholes, infrared background radiation, molecular spectroscopy in the interstellar medium, transit spectroscopy for exoplanet atmosphere, planetary atmosphere in the Solar system, and so on.

In this paper, we investigate the impact of higher-order distortions on the precise measurement of weak gravitational lensing shear and flexion. We begin by defining generalized higher-order distortions and outlining methods for measuring them. Then, using several lens models, we examine how these distortions affect shear and flexion measurements. Our results show that neglecting higher-order distortions can introduce systematic errors of a few percent in both shear and flexion measurements, indicating that these effects cannot be ignored. Although the strength of these errors depends on factors such as lensing strength and the size of background sources, we demonstrate that simultaneous measurement of higher-order distortions can reduce the systematic errors to below 1% in most cases.

Euclid will provide deep NIR imaging to $\sim$26.5 AB magnitude over $\sim$59 deg$^2$ in its deep and auxiliary fields. The Cosmic DAWN survey complements the deep Euclid data with matched depth multiwavelength imaging and spectroscopy in the UV--IR to provide consistently processed Euclid selected photometric catalogs, accurate photometric redshifts, and measurements of galaxy properties to a redshift of $z\sim 10$. In this paper, we present an overview of the survey, including the footprints of the survey fields, the existing and planned observations, and the primary science goals for the combined data set.

The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed: the main science targets, Sgr A* and M87 along with various calibrators. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Königl (BK) jet model. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2-86 GHz range. Then, we studied the dependences of the VLBI core flux density, size, and brightness temperature on the frequency measured in the AGN host frame. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Our results indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations.

We explore the properties of interferometric data from high-redshift 21~cm measurements using the Murchison Widefield Array. These data contain redshifted 21~cm signal, contamination from continuum foreground sources, and radiometric noise. The 21~cm signal from the Epoch of Reionization is expected to be highly-Gaussian, which motivates the use of the power spectrum as an effective statistical tool for extracting astrophysical information. We find that foreground contamination introduces non-Gaussianity into the distribution of measurements, and then use this information to separate Gaussian from non-Gaussian signal. We present improved upper limits on the 21cm EoR power spectrum from the MWA using a Gaussian component of the data, based on the existing analysis from Nunhokee et al (2025). This is extracted as the best-fitting Gaussian to the measured data. Our best 2 sigma (thermal+sample variance) limit for 268 hours of data improves from (30.2~mK)^2 to (23.0~mK)^2 at z=6.5 for the EW polarisation, and from (39.2~mK)^2 to (21.7~mK)^2 = 470~mK^2 in NS. The best limits at z=6.8 (z=7.0) improve to P < (25.9~mK)^2 (P < (32.0~mK)^2), and k = 0.18h/Mpc (k = 0.21h/Mpc). Results are compared with realistic simulations, which indicate that leakage from foreground contamination is a source of the non-Gaussian behaviour.

Little Red Dots (LRDs) are compact, point-like sources characterized by their red color and broad Balmer lines, which have been debated to be either dominated by active galactic nuclei (AGN) or dusty star-forming galaxies (DSFGs). Here we report two LRDs (ID9094 and ID2756) at z$_{\rm spec}$>7, recently discovered in the JWST FRESCO GOODS-North field. Both satisfy the "v-shape" colors and compactness criteria for LRDs and are identified as Type-I AGN candidates based on their broad H$\beta$ emission lines (full width at half maximum: 2280$\pm$490 km/s for ID9094 and 1070$\pm$240 km/s for ID2756) and narrow [OI] lines ($\sim$ 300-400 km/s). To investigate their nature, we conduct deep NOEMA follow-up observations targeting the [CII] 158${\rm \mu m}$ emission line and the 1.3 mm dust continuum. We do not detect [CII] or 1.3 mm continuum emission for either source. Notably, in the scenario that the two LRDs were DSFGs, we would expect significant detections: $>16\sigma$ for [CII] and $>3\sigma$ for the 1.3 mm continuum of ID9094, and $>5\sigma$ for [CII] of ID2756. Using the 3$\sigma$ upper limits of [CII] and 1.3 mm, we perform two analyses: (1) UV-to-FIR spectral energy distribution (SED) fitting with and without AGN components, and (2) comparison of their properties with the L$_{[CII]}$-SFR$_{tot}$ empirical relation. Both analyses are consistent with a scenario where AGN activity may contribute to the observed properties, though a dusty star-forming origin cannot be fully ruled out. Our results highlight the importance of far-infrared observations for studying LRDs, a regime that remains largely unexplored.

The Indian Pulsar Timing Array (InPTA) employs unique features of the upgraded Giant Metrewave Radio Telescope (uGMRT) to monitor dozens of the International Pulsar Timing Array (IPTA) millisecond pulsars (MSPs), simultaneously in the 300-500 MHz and the 1260-1460 MHz bands. This dual-band approach ensures that any frequency-dependent delays are accurately characterized, significantly improving the timing precision for pulsar observations, which is crucial for pulsar timing arrays. We present details of InPTA's second data release that involves 7 yrs of data on 27 IPTA MSPs. This includes sub-banded Times of Arrival (ToAs), Dispersion Measures (DM), and initial timing ephemerides for our MSPs. A part of this dataset, originally released in InPTA's first data release, is being incorporated into IPTA's third data release which is expected to detect and characterize nanohertz gravitational waves in the coming years. The entire dataset is reprocessed in this second data release providing some of the highest precision DM estimates so far and interesting solar wind related DM variations in some pulsars. This is likely to characterize the noise introduced by the dynamic inter-stellar ionised medium much better than the previous release thereby increasing sensitivity to any future gravitational wave search.

We present new high resolution imaging of a light-scattering dust ring and halo around the young star HD 35841. Using spectroscopic and polarimetric data from the Gemini Planet Imager in H-band (1.6 microns), we detect the highly inclined (i=85 deg) ring of debris down to a projected separation of ~12 au (~0.12") for the first time. Optical imaging from HST/STIS shows a smooth dust halo extending outward from the ring to >140 au (>1.4"). We measure the ring's scattering phase function and polarization fraction over scattering angles of 22-125 deg, showing a preference for forward scattering and a polarization fraction that peaks at ~30% near the ansae. Modeling of the scattered-light disk indicates that the ring spans radii of ~60-220 au, has a vertical thickness similar to that of other resolved dust rings, and contains grains as small as 1.5 microns in diameter. These models also suggest the grains have a low porosity, are more likely to consist of carbon than astrosilicates, and contain significant water ice. The halo has a surface brightness profile consistent with that expected from grains pushed by radiation pressure from the main ring onto highly eccentric but still bound orbits. We also briefly investigate arrangements of a possible inner disk component implied by our spectral energy distribution models, and speculate about the limitations of Mie theory for doing detailed analyses of debris disk dust populations.

Solar flares are often accompanied by filament/prominence eruptions ($\sim10^{4}$ K and $\sim 10^{10-11}$ cm$^{-3}$), sometimes leading to coronal mass ejections (CMEs) that directly affect the Earth's environment. `Superflares' are found on some active solar-type (G-type main-sequence) stars, but the association of filament eruptions/CMEs has not been established. Here we show that our optical spectroscopic observation of the young solar-type star EK Draconis reveals the evidence for a stellar filament eruption associated with a superflare. This superflare emitted a radiated energy of $2.0\times10^{33}$ erg, and blue-shifted hydrogen absorption component with a large velocity of $-510$ km s$^{-1}$ was observed shortly after. The temporal changes in the spectra greatly resemble those of solar filament eruptions. Comparing this eruption with solar filament eruptions in terms of the length scale and velocity strongly suggests that a stellar CME occurred. The erupted filament mass of $1.1\times10^{18}$ g is 10 times larger than those of the largest solar CMEs. The massive filament eruption and an associated CME provide the opportunity to evaluate how they affect the environment of young exoplanets/young Earth and stellar mass/angular-momentum evolution.

The Laser Interferometer Gravitational-wave Observatory Scientific Collaboration and Virgo Collaboration (LVC) sent out 56 gravitational-wave (GW) notices during the third observing run (O3). Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM) performed optical and near-infrared observations to identify and observe an electromagnetic (EM) counterpart. We constructed web-based system which enabled us to obtain and share information of candidate host galaxies for the counterpart, and status of our observations. Candidate host galaxies were selected from the GLADE catalog with a weight based on the three-dimensional GW localization map provided by LVC. We conducted galaxy-targeted and wide-field blind surveys, real-time data analysis, and visual inspection of observed galaxies. We performed galaxy-targeted follow-ups to 23 GW events during O3, and the maximum probability covered by our observations reached to 9.8%. Among them, we successfully started observations for 10 GW events within 0.5 days after the detection. This result demonstrates that our follow-up observation has a potential to constrain EM radiation models for a merger of binary neutron stars at a distance of up to $\sim$100~Mpc with a probability area of $\leq$ 500~deg$^2$.

The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets' orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{\rm{Jup}}$ and $7~M_{\rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 \pm 5$~Myr, we found the mass of planet b to be $5.8 \pm 0.5~M_{\rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{\rm{Jup}}$ each.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

CASA, the Common Astronomy Software Applications package, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and NSF's Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can process data from both single-dish and aperture-synthesis telescopes, and one of its core functionalities is to support the data reduction and imaging pipelines for ALMA, VLA and the VLA Sky Survey (VLASS). CASA has recently undergone several exciting new developments, including an increased flexibility in Python (CASA 6), support of Very Long Baseline Interferometry (VLBI), performance gains through parallel imaging, data visualization with the new Cube Analysis Rendering Tool for Astronomy (CARTA), enhanced reliability and testing, and modernized documentation. These proceedings of the 2019 Astronomical Data Analysis Software & Systems (ADASS) conference give an update of the CASA project, and detail how these new developments will enhance user experience of CASA.

We present the results from a full polarization study carried out with ALMA during the first VLBI campaign, which was conducted in Apr 2017 in the $\lambda$3mm and $\lambda$1.3mm bands, in concert with the Global mm-VLBI Array (GMVA) and the Event Horizon Telescope (EHT), respectively. We determine the polarization and Faraday properties of all VLBI targets, including Sgr A*, M87, and a dozen radio-loud AGN. We detect high linear polarization fractions (2-15%) and large rotation measures (RM $>10^{3.3}-10^{5.5}$ rad m$^{-2}$). For Sgr A* we report a mean RM of $(-4.2\pm0.3) \times10^5$ rad m$^{-2}$ at 1.3 mm, consistent with measurements over the past decade, and, for the first time, an RM of $(-2.1\pm0.1) \times10^5$ rad m$^{-2}$ at 3 mm, suggesting that about half of the Faraday rotation at 1.3 mm may occur between the 3 mm photosphere and the 1.3 mm source. We also report the first unambiguous measurement of RM toward the M87 nucleus at mm wavelengths, which undergoes significant changes in magnitude and sign reversals on a one year time-scale, spanning the range from -1.2 to 0.3 $\times\,10^5$ rad m$^{-2}$ at 3 mm and -4.1 to 1.5 $\times\,10^5$ rad m$^{-2}$ at 1.3 mm. Given this time variability, we argue that, unlike the case of Sgr A*, the RM in M87 does not provide an accurate estimate of the mass accretion rate onto the black hole. We put forward a two-component model, comprised of a variable compact region and a static extended region, that can simultaneously explain the polarimetric properties observed by both the EHT and ALMA. These measurements provide critical constraints for the calibration, analysis, and interpretation of simultaneously obtained VLBI data with the EHT and GMVA.

Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction independent calibration error due to ionospheric offsets. Our calibration strategy is optimized for the lowest frequency bands by reducing the number of direction dependent calibrators and taking into account radio sources within a wider field of view. We remove data polluted by systematics based on the RFI occupancy and ionospheric conditions, finally selecting 5.5 hours of the cleanest data. Using these data, we obtain two sigma upper limits on the 21 cm power spectrum in the range of $0.1\lessapprox k \lessapprox 1 ~\rm ~h~Mpc^{-1}$ and at $z$=14.2, 15.2 and 16.5, with the lowest limit being $6.3\times 10^6 ~\rm mK^2$ at $\rm k=0.14 \rm ~h~Mpc^{-1}$ and at $z=15.2$ with a possibility of a few \% of signal loss due to direction independent calibration.

In 2017, the next generation Very Large Array (ngVLA) Science Advisory Council, together with the international astronomy community, developed a set of five Key Science Goals (KSGs) to inform, prioritize and refine the technical capabilities of a future radio telescope array for high angular resolution operation from 1.2 - 116 GHz with 10 times the sensitivity of the Jansky VLA and ALMA. The resulting KSGs, which require observations at centimeter and millimeter wavelengths that cannot be achieved by any other facility, represent a small subset of the broad range of astrophysical problems that the ngVLA will be able address. This document presents an update to the original ngVLA KSGs, taking account of new results and progress in the 7+ years since their initial presentation, again drawing on the expertise of the ngVLA Science Advisory Council and the broader community in the ngVLA Science Working Groups. As the design of the ngVLA has also matured substantially in this period, this document also briefly addresses initial expectations for ngVLA data products and processing that will be needed to achieve the KSGs. The original ngVLA KSGs endure as outstanding problems of high priority. In brief, they are: (1) Unveiling the Formation of Solar System Analogues; (2) Probing the Initial Conditions for Planetary Systems and Life with Astrochemistry; (3) Charting the Assembly, Structure, and Evolution of Galaxies from the First Billion Years to the Present; (4) Science at the Extremes: Pulsars as Laboratories for Fundamental Physics; (5) Understanding the Formation and Evolution of Stellar and Supermassive Black Holes in the Era of Multi-Messenger Astronomy.

The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., prompt-magnetic-jet and delayed-magnetic-jet, for which the magnetic fields of proto-neutron stars (PNSs) during collapse and the core-bounce are strong and comparatively moderate, respectively. Following the hydrodynamical trajectories of each explosion model, we confirmed that r-processes successfully occur in the prompt-magnetic-jets, which produce heavy nuclei including actinides. On the other hand, the r-process in the delayed-magnetic-jet is suppressed, which synthesizes only nuclei up to the second peak ($A \sim 130$). Thus, the r-process in the delayed-magnetic-jets could explain only "weak r-process" patterns observed in metal-poor stars rather than the "main r-process", represented by the solar abundances. Our results imply that core-collapse supernovae are possible astronomical sources of heavy r-process elements if their magnetic fields are strong enough, while weaker magnetic explosions may produce "weak r-process" patterns ($A \lesssim 130$). We show the potential importance and necessity of magneto-rotational supernovae for explaining the galactic chemical evolution, as well as abundances of r-process enhanced metal-poor stars. We also examine the effects of the remaining uncertainties in the nature of PNSs due to weak interactions that determine the final neutron-richness of ejecta. Additionally, we briefly discuss radioactive isotope yields in primary jets (e.g., $^{56}$Ni), with relation to several optical observation of SNe and relevant high-energy astronomical phenomena.

Red giants undergo dramatic and complex structural transformations as they evolve. Angular momentum is transported between the core and envelope during this epoch, a poorly understood process. Here, we infer envelope and core rotation rates from Kepler observations of $\sim$1517 red giants. While many measurements are consistent with the existing studies, our investigation reveals systematic changes in the envelope-to-core rotation ratio and we report the discovery of anomalies such as clump stars with rapidly rotating cores, and red giants with envelopes rotating faster than their cores. We propose binary interactions as a possible mechanism by which some of these cores and envelopes are spun up. These results pose challenges to current theoretical expectations and can have major implications for compact remnants born from stellar cores.