Binding energy (BE) is a critical parameter in astrochemical modeling, governing the retention of species on interstellar dust grains and their subsequent chemical evolution. However, conventional models often rely on single-valued BEs, overlooking the intrinsic distribution arising from diverse adsorption sites. In this study, we present BEs for monohydric alcohols, thiols, and their plausible precursors, including aldehydes and thioaldehydes. We incorporate a distribution of BEs to capture the realistic variation in adsorption strengths. The quantum chemical calculations provide a range of BE values rather than a single estimate, ensuring a more precise description of molecular diffusion and surface chemistry. The BE trend of analogous species provides qualitative insight into the dominant reaction pathways and key precursors that drive the formation of larger molecules under interstellar conditions. Oxygen-bearing species generally exhibit higher BEs than their sulfur analogues, primarily due to stronger interactions, further influencing molecular adsorption and reactivity. We implemented BE distributions in astrochemical models, revealing significant effects on predicted abundances and establishing a more accurate framework for future astrochemical modeling.

Solar flares are complex phenomena driven by the release of magnetic energy, but a large energy reservoir is not sufficient to determine their eruptive potential; the magnetic topology and plasma dynamics play a key role. We investigate the thermodynamic and magnetic properties of the solar atmosphere during the rise, peak, and decay phases of a C5.1-class flare and filament eruption in active region NOAA 12561 on 2016 July 7, to understand the origin and atmospheric response of this event. High spatial and spectral resolution spectropolarimetric observations of the chromospheric Ca II 8542A line and nearby photospheric lines were obtained with the TRIPPEL-SP spectropolarimeter at the Swedish 1-m Solar Telescope. Using non-local thermodynamic equilibrium (NLTE) inversions and non-force-free field (NFFF) magnetic extrapolations, we followed the event's evolution from its precursor to its decay. Before the flare, our analysis reveals a complex, sheared magnetic topology with a high free energy content ($\sim2\times10^{30}$ erg). In this precursor phase, we detected persistent, localized heating (temperature increase of $\sim$2000 K) with strong downflows ($\sim$10-20 km/s) deep in the atmosphere. This heating was co-spatial with a bald-patch region, suggesting that low-altitude magnetic reconnection could destabilize the filament of the region. The flare's rise phase was marked by the filament's eruption, with a total speed larger than $\sim$70 km/s, when combining inversions and plane-of-sky motions. Following the eruption, the free energy decreased by $\sim$30$\%$ as post-flare loops formed, connecting the flare ribbons and channeling the released energy into the lower atmosphere. The flare ribbons exhibited significant heating to $\sim$8500 K and downflows up to $\sim$10 km/s, consistent with energy deposition along reconnected loops.

Understanding the origin and evolution of carbon-enhanced metal-poor (CEMP) stars is key to tracing the Galaxy's early chemical enrichment. We investigate how realistic 3D radiation-hydrodynamics (RHD) model atmospheres affect carbon abundances in CEMP stars and implications for their classification and Galactic chemical evolution (GCE). We focus on biases from traditional 1D hydrostatic models. We use the M3DIS code to compute 3D RHD model atmospheres for main-sequence and sub-giant stars over a wide range of metallicities and carbon enhancements. Synthetic spectra of the CH G-band are calculated with 3D radiative transfer and compared to spectra from classical 1D MARCS models. We derive abundance corrections and apply them to a large SAGA database sample to quantify effects on the carbon abundance distribution and CEMP classification. Our new 3D CEMP models predict cooler upper atmospheric layers than in 1D models, resulting in stronger CH absorption and lower inferred carbon abundances by up to -0.9 dex at the lowest metallicities. Carbon enhancement in the atmosphere itself increases molecular opacities and leads to radiative re-heating, partly offsetting adiabatic cooling in 3D models and reducing 3D-1D abundance corrections. Applying these corrections lowers the CEMP fraction by up to 20% below [Fe/H]=-3 and alters the relative contribution of CEMP sub-classes. The fraction of CEMP-no stars increases while the number of CEMP-r/s stars decreases, due to the downward revision of absolute carbon abundances. These changes bring the Galactic carbon distribution into better agreement with GCE models assuming a 20% contribution from faint supernovae. Realistic model atmospheres are essential to reliably reconstruct the Galaxy's early chemical enrichment history.

The nature and detailed properties of the heating of the million-degree solar corona are important issues that are still largely unresolved. Nanoflare heating might be dominant in active regions and quiet Sun, although direct signatures of such small-scale events are difficult to observe in the highly conducting, faint corona. The aim of this work is to test the theory of coronal heating by nanoflares in braided magnetic field structures. We analyze a 3D MHD model of a multistrand flux tube in a stratified solar atmosphere, driven by twisting motions at the boundaries. We show how the magnetic structure is maintained at high temperature and for an indefinite time, by intermittent episodes of local magnetic energy release due to reconnection. We forward-modelled optically thin emission with SDO/AIA and MUSE and compared the synthetic observations with the intrinsic coronal plasma properties, focusing on the response to impulsive coronal heating. Currents build up and their impulsive dissipation into heat are also investigated through different runs. In this first paper, we describe the proliferation of heating from the dissipation of narrow current sheets in realistic simulations of braided coronal flux tubes at unprecedented high spatial resolutions.

Scale invariance is a hallmark of many natural systems, including solar flares, where energy release spans a vast range of scales. Recent computational advances, at the level of both algorithmics and hardware, have enabled high-resolution magnetohydrodynamical (MHD) simulations to span multiple scales, offering new insights into magnetic energy dissipation processes. Here, we study scale invariance of magnetic energy dissipation in two distinct MHD simulations. Current sheets are identified and analyzed over time. Results demonstrate that dissipative events exhibit scale invariance, with power-law distributions characterizing their energy dissipation and lifetimes. Remarkably, these distributions are consistent across the two simulations, despite differing numerical and physical setups, suggesting universality in the process of magnetic energy dissipation. Comparisons between the evolution of dissipation regions reveals distinct growth behaviors in high plasma-beta regions (convective zone) and low plasma-beta regions (atmosphere). The latter display spatiotemporal dynamics similar to those of avalanche models, suggesting self-organized criticality and a common universality class.

Advancements in instrumentation have revealed a multitude of small-scale EUV events in the solar atmosphere. Our aim is to employ high-resolution magnetograms to gain a detailed understanding of the magnetic origin of such phenomena. We have used coordinated observations from SST, IRIS, and SDO to analyze an ephemeral magnetic flux emergence episode and the following chain of small-scale energetic events. These unique observations clearly link these phenomena together. The high-resolution (0."057/pixel) magnetograms obtained with SST/CRISP allows us to reliably measure the magnetic field at the photosphere and detect the emerging bipole that causes the subsequent eruptive atmospheric events. Notably, this small-scale emergence episode remains indiscernible in the lower resolution SDO/HMI magnetograms (0."5/pixel). We report the appearance of a dark bubble in Ca II K related to the emerging bipole, a sign of the canonical expanding magnetic dome predicted in flux emergence simulations. Evidences of reconnection are also found: first through an Ellerman bomb, and later by the launch of a surge next to a UV burst. The UV burst exhibits a weak EUV counterpart in the coronal SDO/AIA channels. By calculating DEM, its plasma is shown to reach a temperature beyond 1 MK and have densities between the upper chromosphere and transition region. Our study showcases the importance of high-resolution magnetograms to unveil the mechanisms triggering phenomena such as EBs, UV bursts, and surges. This could hold implications for small-scale events akin to those recently reported in EUV using Solar Orbiter. The finding of temperatures beyond 1 MK in the UV burst plasma strongly suggests that we are examining analogous features. Therefore, we signal caution regarding drawing conclusions from full-disk magnetograms that lack the necessary resolution to reveal their true magnetic origin.

Millimter (mm) frequencies are primarily sensitive to thermal emission from layers across the stellar chromosphere up to the transition region, while metrewave (radio) frequencies probe the coronal heights. Together the mm and radio band spectroscopic snapshot imaging enables the tomographic exploration of the active atmospheric layers of the cool main-sequence stars (spectral type: FGKM), including our Sun. Sensitive modern mm and radio interferometers let us explore solar/stellar activity covering a range of energy scales at sub-second and sub-MHz resolution over wide operational bandwidths. The superior uv-coverage of these instruments facilitate high dynamic range imaging, letting us explore the morphological evolution of even energetically weak events on the Sun at fine spectro-temporal cadence. This article will introduce the current advancements, the data analysis challenges and available tools. The impact of these tools and novel data in field of solar/stellar research will be summarised with future prospects.

This paper proposes that the ion temperature is several times the local electron temperature in the hot onset phase and at the above-the-loop region of solar flares. The paper considers: the evidence of spectral line Doppler widths ("non-thermal" broadening); evidence for "universal" ion and electron temperature increase scaling relations for magnetic reconnection in the solar wind, Earth's magnetopause, Earth's magnetotail and numerical simulations; and thermal equilibration times for onset and above-the-loop densities, which are much longer than previous estimates based on soft X-ray flare loops. We conclude that the ion temperature is likely to reach 60 MK or greater and that it may represent a substantial part of spectral line widths, significantly contributing to solving the long-standing issue of the excess nonthermal broadening in flare lines.

Large amplitude oscillations commonly occur in solar prominences, triggered by energetic phenomena such as jets and flares. On March 14-15, 2015, a filament partially erupted in two stages, leading to oscillations in different parts. This study explores longitudinal oscillations from the eruption, focusing on the mechanisms behind their initiation, with special attention to the large oscillation on March 15. The oscillations and jets are analyzed using the time-distance technique. For flares and their interaction with the filament, we analyze AIA channels and use the DEM technique. Initially, a jet fragments the filament, splitting it into two segments. One remains in place, while the other detaches and moves. This causes oscillations in both segments: (a) the position change causes the detached segment to oscillate with a period of $69 \pm 3$ minutes; (b) the jet flows cause the remaining filament to oscillate with a period of $62 \pm 2$ minutes. In the second phase, on March 15, another jet seemingly activates the detached filament eruption, followed by a flare. A large longitudinal oscillation occurs in the remnant segment with a period of $72 \pm 2$ minutes and velocity amplitude $73 \pm 1 \, \mathrm{km s^{-1}}$. During the oscillation trigger, bright field lines connect the flare with the filament, appearing only in the AIA 131$Å$ and 94$Å$ channels, indicating the presence of hot plasma. DEM analysis confirms this, showing plasma around 10 MK pushing the prominence from its southeastern side, displacing it along the field lines and starting the oscillation. From this, the flare -- not the preceding jet-triggers the oscillation. The hot plasma flows into the filament channel. We explain how flares trigger large oscillations in filaments by proposing that post-flare loops reconnect with the filament channel's magnetic field.

Total solar eclipses (TSEs) provide a unique opportunity to observe the large-scale solar corona. The solar wind plays an important role in forming the large-scale coronal structure and magnetohydrodynamic (MHD) simulations are used to reproduce it for further studying coronal mass ejections (CMEs). We conduct a data-constrained MHD simulation of the global solar corona including solar wind effects of the 2024 April 8 TSE with observed magnetograms using the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC) within 2.5 $R_\odot$. This TSE happened within the solar maximum, hence the global corona was highly structured. Our MHD simulation includes the energy equation with a reduced polytropic index $\gamma=1.05$. We compare the global magnetic field for multiple magnetograms and use synchronic frames from the Solar Dynamics Observatory/Helioseismic and Magnetic Imager to initialize the magnetic field configuration from a magneto-frictionally equilibrium solution, called the Outflow field. We detail the initial and boundary conditions employed to time-advance the full set of ideal MHD equations such that the global corona is relaxed to a steady state. The magnetic field, the velocity field, and distributions of the density and thermal pressure are successfully reproduced. We demonstrate direct comparisons with TSE images in white-light and Fe XIV emission augmented with quasi-separatrix layers, the integrated current density, and the synthetic white-light radiation, and find a good agreement between simulations and observations. This provides a fundamental background for future simulations to study the triggering and acceleration mechanisms of CMEs under solar wind effects.

Coronal jets are ubiquitous, collimated million-degree ejections that contribute to the energy and mass supply of the upper solar atmosphere and the solar wind. Solar Orbiter provides an unprecedented opportunity to observe fine-scale jets from a unique vantage point close to the Sun. We aim to uncover thin jets originating from Coronal Bright Points (CBPs) and investigate observable features of plasmoid-mediated reconnection. We analyze eleven datasets from the High Resolution Imager 174 Å of the Extreme Ultraviolet Imager (HRIEUV) onboard Solar Orbiter, focusing on narrow jets from CBPs and signatures of magnetic reconnection within current sheets and outflow regions. To support the observations, we compare with CBP simulations performed with the Bifrost code. We have identified thin coronal jets originating from CBPs with widths ranging from 253 km to 706 km: scales that could not be resolved with previous EUV imaging instruments. Remarkably, these jets are 30-85% brighter than their surroundings and can extend up to 22 Mm while maintaining their narrow form. In one of the datasets, we directly identify plasmoid-mediated reconnection through the development within the current sheet of a small-scale plasmoid that reaches a size of 332 km and propagates at 40 km/s. In another dataset, we infer plasmoid signatures through the intermittent boomerang-like pattern that appears in the outflow region. Both direct and indirect plasmoid-mediated reconnection signatures are supported by comparisons with the synthetic HRIEUV emission from the simulations.

On April 17, 2024, the third successful Hi-C sounding rocket flight, Hi-C Flare, recorded coronal images in Fe XXI 129 A emission from 11 MK plasma during the post-maximum phase of an M1.6-class solar flare, achieving unprecedented spatial (~300 km) and temporal (1.3 s) resolutions. The flare started at 21:55 UT, peaked at 22:08 UT, and lasted ~40 minutes. Hi-C observed for over five minutes (22:15:45 to 22:21:25), starting roughly eight minutes after flare maximum. A sudden compact bright burst - 875 +/- 25 km wide, lasting 90 +/- 1.3 s, displaying a plane-of-sky motion of ~50 km/s toward the loop apex, and splitting into two toward the end - occurs near the foot of some post-flare loops. Its size and brightness are reminiscent of flare-ribbon kernels during a flare's rapid rise phase, kernels marking sites of sudden heating and hot plasma upflow, making its occurrence during the late phase surprising. Such isolated brightenings in a flare's post-maximum phase are rare, and have not been previously reported. The kernel was detected in all SDO/AIA channels. Its 1600 A light curve peaked ~50 s earlier than its 131 A light curve, similar to that of flare-ribbon kernels, albeit with a smaller delay of ~25 s, during the impulsive phase of the flare. In SDO/HMI magnetograms, the kernel sits in unipolar positive magnetic flux near an embedded clump of negative flux. Although localized magnetic reconnection within the kernel (a microflare) cannot be ruled out for its cause, the observations favor the localized brightening being an isolated, exceptionally late flare-ribbon kernel, resulting from an exceptionally late burst of the flare's coronal reconnection.

Context. EUV late phase (ELP) flares exhibit a second peak in warm coronal emissions minutes to hours after the main peak of the flare. This phase is all but negligible, yet it is still poorly understood what role it plays across the solar cycle and what governs it. Aims. We present a statistical analysis of ELP flares over four years between May 2010 and May 2014 based on properties like eruptivity, magnetic configuration, and late-phase duration, delay and strength in order to understand what influences the likelihood of this class of flares and their behavior on a general scale. Methods. We primarily make use of data from the Solar Dynamics Observatory's (SDO) Extreme-ultraviolet Variability Experiment (EVE), as well as complimentary spatial information provided by the Atmospheric Imaging Assembly (AIA), to assess relationships between the various parameters and see if ELP flares differ from the general flare population. We quantify the criteria for ELP flare definition and determine its characteristics. Results. Our analysis shows that about 10% of all flares with a GOES class greater than or equal to C3.0 experience an EUV late phase (179 out of 1803). This percentage decreases from solar minimum to solar maximum. C-class flares are considerably less likely to be identified as ELP flares than their higher-energetic counterparts, which is in line with previous investigations. The majority of this type of flares is confined (67%), more so than in the general flare population (greater than or equal to C5.0). There appears to be a (linear) relationship between the late-phase delay and its duration. The ratio of the emission peak of the late and main flare phase lies between 0.3 and 5.9, and exceeds 1 in 71.5% of cases, which is considerably higher than previously reported.

High-mass stars, born in massive dense cores (MDCs), profoundly impact the cosmic ecosystem through feedback processes and metal enrichment, yet little is known about how MDCs assemble and transfer mass across scales to form high-mass young stellar objects (HMYSOs). Using multi-scale (40-2500 au) observations of an MDC hosting an HMYSO, we identify a coherent dynamical structure analogous to barred spiral galaxies: three 20,000 au spiral arms feed a 7,500 au central bar, which channels gas to a 2,000 au pseudodisk. Further accretion proceeds through the inner structures, including a Keplerian disk and an inner disk (100 au), which are thought to be driving a collimated bipolar outflow. This is the first time that these multi-scale structures (spiral arms, bar, streamers, envelope, disk, and outflow) have been simultaneously observed as a physically coherent structure within an MDC. Our discovery suggests that well-organized hierarchical structures play a crucial role during the gas accretion and angular momentum build-up of a massive disk.

Active region recurrent jets are manifestations of episodic magnetic energy release processes driven by complex interactions in the lower solar atmosphere. While magnetic flux emergence and cancellation are widely recognized as key contributors to jet formation, the mechanisms behind repeated magnetic reconnection remain poorly understood. In this letter, we report a sequence of nine recurrent jets originating from active region AR 12715 during its decay phase, where the jet activity was associated with a complex distribution of fragmented magnetic flux. Non-linear force-free field (NLFFF) extrapolations reveal the presence of low-lying, current-carrying loops beneath overarching open magnetic fields near the jet footpoints. These magnetic structures were perturbed by (i) emerging flux elements and (ii) interactions between oppositely polarized moving magnetic features (MMFs). To interpret these observations, we compare them with 3D radiative MHD simulation from the Bifrost model, which reproduce jet formation driven by interacting bipolar MMFs, leading to subsequent flux cancellation in the photosphere. Our results emphasize the critical role of MMF-driven flux interactions in initiating and sustaining recurrent jet activity in active regions.

Explosive transient events occur throughout the solar atmosphere. The differing manifestations range from coronal mass ejections to Ellerman bombs. The former may have negligible signatures in the lower atmosphere, and the latter may have negligible nonthermal emissions such as hard X-radiation. A solar flare generally involves a broad range of emission signatures. Using a suite of four space-borne telescopes, we report a solar event that combines aspects of simple UV bursts and hard X-ray emitting flares at the same time. The event is a compact C-class flare in active region AR11861, SOL2013-10-12T00:30. By fitting a combined isothermal and nonthermal model to the hard X-ray spectrum, we inferred plasma temperatures in excess of 15\,MK and a nonthermal power of about $3\times10^{27}$\,erg\,s$^{-1}$ in this event. Despite these high temperatures and evidence for nonthermal particles, the flare was mostly confined to the chromosphere. However, the event lacked clear signatures of UV spectral lines, such as the Fe\,{\sc xii} 1349\,\AA\ and Fe\,{\sc xxi} 1354\,\AA\ emission lines, which are characteristic of emission from hotter plasma with a temperature over 1\,MK. Moreover, the event exhibited very limited signatures in the extreme-UV wavelengths. Our study indicates that a UV burst -- hard X-ray flare hybrid phenomenon exists in the low solar atmosphere. Plasma that heats to high temperatures coupled with particle acceleration by magnetic energy that is released directly in the lower atmosphere sheds light on the nature of active region core heating and on inferences of flare signatures.

Full-Stokes polarimetric datasets, originating from slit-spectrograph or narrow-band filtergrams, are routinely acquired nowadays. The data rate is increasing with the advent of bi-dimensional spectropolarimeters and observing techniques that allow long-time sequences of high-quality observations. There is a clear need to go beyond the traditional pixel-by-pixel strategy in spectropolarimetric inversions by exploiting the spatiotemporal coherence of the inferred physical quantities. We explore the potential of neural networks as a continuous representation of the physical quantities over time and space (also known as neural fields), for spectropolarimetric inversions. We have implemented and tested a neural field to perform the inference of the magnetic field vector (approach also known as physics-informed neural networks) under the weak-field approximation (WFA). By using a neural field to describe the magnetic field vector, we can regularize the solution in the spatial and temporal domain by assuming that the physical quantities are continuous functions of the coordinates. We investigated the results in synthetic and real observations of the Ca II 8542 A line. We also explored the impact of other explicit regularizations, such as using the information of an extrapolated magnetic field, or the orientation of the chromospheric fibrils. Compared to the traditional pixel-by-pixel inversion, the neural field approach improves the fidelity of the reconstruction of the magnetic field vector, especially the transverse component. This implicit regularization is a way of increasing the effective signal-to-noise of the observations. Although it is slower than the pixel-wise WFA estimation, this approach shows a promising potential for depth-stratified inversions, by reducing the number of free parameters and inducing spatio-temporal constraints in the solution.

The solar atmosphere is a complex environment with diverse species and varying ionization states, especially in the chromosphere, where significant ionization variations occur. This region transitions from highly collisional to weakly collisional states, leading to complex plasma state transitions influenced by magnetic strengths and collisional properties. These processes introduce numerical stiffness in multi-fluid models, imposing severe timestep restrictions on standard time integration methods. New numerical methods are essential to address these computational challenges, effectively managing the diverse timescales in multi-fluid and multi-physics models. The widely used time operator splitting technique offers a straightforward approach but requires careful timestep management to avoid stability issues and errors. Despite some studies on splitting errors, their impact on solar and stellar astrophysics is often overlooked. We focus on a Multi-Fluid Multi-Species (MFMS) model, which presents significant challenges for time integration. We propose a second-order Partitioned Implicit-Explicit Runge-Kutta (PIROCK) method that combines efficient explicit and implicit integration techniques with variable time-stepping and error control. Compared to a standard third-order explicit method and a first-order Lie splitting approach, the PIROCK method shows robust advantages in accuracy, stability, and computational efficiency. Our results reveal PIROCK's capability to solve multi-fluid problems with unprecedented efficiency. Preliminary results on chemical fractionation represent a significant step toward understanding the First-Ionization-Potential (FIP) effect in the solar atmosphere.

Context: Flux emergence in the solar atmosphere is a complex process that causes a release of magnetic energy as heat and acceleration of solar plasma at a variety of spatial scales. Methods: We analysed imaging spectropolarimetric data taken in the He 1083 nm line. This line is sensitive to temperatures larger than 15 kK, unlike diagnostics such as \MgIIhk, \CaIIHK, and \Halpha, which lose sensitivity already at 15 kK. The He I data is complemented by imaging spectropolarimetry in the \CaIIK, \Feline, and \Caline\ lines. We employed inversions to determine the magnetic field and vertical velocity in the solar atmosphere. We computed He 1083 nm profiles from a radiation-MHD simulation of the solar atmosphere to help interpretation of the observations. Results: We find fast-evolving blob-like emission features in the He 1083 nm triplet at locations where the magnetic field is rapidly changing direction, and these are likely sites of magnetic reconnection. We fit the line with a model consisting of an emitting layer located below a cold layer representing the fibril canopy. The modelling provides evidence that this model, while simple, catches the essential characteristics of the line formation. The morphology of the emission in the He 1083 nm line is localized and blob-like, unlike the emission in the \CaIIK\ line, which is more filamentary. Conclusions: The modelling shows that the \Heline\ emission features and their Doppler shifts can be caused by opposite-polarity reconnection and/or horizontal current sheets below the canopy layer in the chromosphere. Based on the high observed Doppler width and the blob-like appearance of the emission features, we conjecture that at least a fraction of them are produced by plasmoids. We conclude that transition-region-like temperatures in the deeper layers of the active region chromosphere are more common than previously thought.

We present centimeter continuum observations of seven high luminosity massive protostars and their surroundings as part of the SOFIA Massive (SOMA) Star Formation Survey. With data from the Very Large Array and the Australia Telescope Compact Array, we analyze the spectral index, morphology and multiplicity of the detected radio sources. The high sensitivity and high resolution observations allow us to resolve many sources. We report thirteen new detections and two previously known detections that we observed for the first time in radio frequencies. We use the observations to build radio spectral energy distributions (SEDs) to calculate spectral indices. With radio morphologies and the spectral indices, we give assessments on the nature of the sources, highlighting five sources that display a radio jet-like morphology and a spectral index consistent with ionized jets. Combining with the SOMA Radio I sample, we present the radio - bolometric luminosity relation, especially probing the regime from $L_{\rm bol}\sim 10^4$ to $10^6\:L_\odot$. Here we find a steep rise in radio luminosity, which is expected by models that transition from shock ionization to photoionization.