More than 200 extrasolar planets have been discovered around relatively nearby stars, primarily through the Doppler line shifts owing to the reflex motions of their host stars, and more recently through transits of some planets across the face of the host stars. The detection of planets with the shortest known periods, 1.2 to 2.5 days, has mainly resulted from transit surveys which have generally targeted stars more massive than 0.75 M_sun. Here we report the results from a planetary transit search performed in a rich stellar field towards the Galactic bulge. We discovered 16 candidates with orbital periods between 0.4 and 4.2 days, five of which orbit stars of 0.44 to 0.75 M_sun. In two cases, radial-velocity measurements support the planetary nature of the companions. Five candidates have orbital periods below 1.0 day, constituting a new class of ultra-short-period planets (USPPs), which occur only around stars of less than 0.88 M_sun. This indicates that those orbiting very close to more luminous stars might be evaporatively destroyed, or that jovian planets around lower-mass stars might migrate to smaller radii.

The ocean and the atmosphere, and hence the climate, are governed at large scale by interactions between pressure gradient, Coriolis and buoyancy forces. This leads to a quasi-geostrophic balance in which, in a two-dimensional-like fashion, the energy injected by solar radiation, winds or tides goes to large scales in what is known as an inverse cascade. Yet, except for Ekman friction, energy dissipation and turbulent mixing occur at small scale implying the formation of such scales associated with breaking of geostrophic dynamics through wave-eddy interactions \cite{ledwell_00, vanneste_13} or frontogenesis \cite{hoskins_72, mcwilliams_10}, in opposition to the inverse cascade. Can it be both at the same time? We exemplify here this dual behavior of energy with the help of three-dimensional direct numerical simulations of rotating stratified Boussinesq turbulence. We show that efficient small-scale mixing and large-scale coherence develop simultaneously in such geophysical and astrophysical flows, both with constant flux as required by theoretical arguments, thereby clearly resolving the aforementioned contradiction.

Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombination times of ionised hydrogen are longer than the dynamical timescales, requiring a non-equilibrium (NE) treatment of hydrogen ionisation. The MURaM code is extended to include the physical process required for accurate simulation of the solar chromosphere, as implemented in the Bifrost code. This includes a time-dependent treatment of hydrogen ionisation, a scattering multi-group radiation transfer scheme and approximations for NLTE radiative cooling. The inclusion of NE and NLTE physics has a large impact on the structure of the chromosphere; the NE treatment of hydrogen ionisation leads to a higher ionisation fraction and enhanced populations in the first excited state throughout cold inter-shock regions of the chromosphere. Additionally this prevents hydrogen ioniation from buffering energy fluctuations, leading to hotter shocks and cooler inter-shock regions. The hydrogen populations in the ground and first excited state are enhanced by $10^2-10^3$ in the upper chromosphere and up to $10^9$ near the transition region. Including the necessary NLTE physics leads to significant differences in chromospheric structure and dynamics. The thermodynamics and hydrogen populations calculated using the extended version of the MURaM code are consistent with previous non-equilibrium simulations. The electron number and temperature calculated using the non-equilibrium treatment of the chromosphere are required to accurately synthesise chromospheric spectral lines.

We describe the design of the Sunrise Filter Imager (SuFI) and the Image Stabilization and Light Distribution (ISLiD) unit onboard the Sunrise balloon borne solar observatory. This contribution provides the necessary information which is relevant to understand the instruments working principles, the relevant technical data, and the necessary information about calibration issues directly related to the science data.

Inertial modes have been recently detected in the Sun via helioseismology, yet their origin, evolution, and role in the dynamics of the solar plasma and magnetic field remain poorly understood. In this study, we employ global numerical simulations to investigate the excitation mechanisms and dynamical consequences of inertial modes in the Sun and stellar interiors. We validate first our numerical setup by analyzing the evolution of sectoral and tesseral perturbations imposed on a rigidly rotating sphere. The results confirm that a perturbation of a given mode can excite neighboring modes with both smaller and larger wavenumbers along the dispersion relation of Rossby waves. Subsequently, we use a physically motivated forcing to impose differential rotation with varying shear amplitudes, and examine the spontaneous onset and nonlinear evolution of inertial modes. The simulations reveal that the growth of velocity perturbations is primarily driven by baroclinic instability. It gives rise to high-latitude inertial modes in the form of retrograde polar vortices whose properties depend on the imposed shear. Equatorial Rossby modes are also excited, albeit with lower intensity than their high-latitude counterpart. Perturbations with arbitrary azimuthal wavenumbers lead to the excitation of Rossby modes for all available wave numbers, sustained by both direct and inverse energy cascades. In simulations with stronger shear, the high latitude modes produce Reynolds stresses able to modify the imposed differential rotation and accelerate the rotation of the poles.

Non-linear force-free extrapolations are a common approach to estimate the 3D topology of coronal magnetic fields based on photospheric vector magnetograms. The force-free assumption is a valid approximation at coronal heights, but for the dense plasma conditions in the lower atmosphere, this assumption is not satisfied. In this study, we utilize multi-height magnetic field measurements in combination with physics-informed neural networks to advance solar magnetic field extrapolations. We include a flexible height-mapping, which allows us to account for the different formation heights of the observed magnetic field measurements. The comparison to analytical and simulated magnetic fields demonstrates that including chromospheric magnetic field measurements leads to a significant improvement of our magnetic field extrapolations. We also apply our method to chromospheric line-of-sight magnetograms, from the Vector Spectromagnetograph (VSM) on the Synoptic Optical Long-term Investigations of the Sun (SOLIS) observatory, in combination with photospheric vector magnetograms, from the Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamic Observatory (SDO). The comparison to observations in extreme ultraviolet wavelengths shows that the additional chromospheric information leads to a better agreement with the observed coronal structures. In addition, our method intrinsically provides an estimate of the corrugation of the observed magnetograms. With this new approach, we make efficient use of multi-height magnetic field measurements and advance the realism of coronal magnetic field simulations.

Coronagraphic observations enable direct monitoring of coronal mass ejections (CMEs) through scattered light from free electrons, but determining the 3D plasma distribution from 2D imaging data is challenging due to the optically-thin plasma and the complex image formation processes. We introduce SuNeRF-CME, a framework for 3D tomographic reconstructions of the heliosphere using multi-viewpoint coronagraphic observations. The method leverages Neural Radiance Fields (NeRFs) to estimate the electron density in the heliosphere through a ray-tracing approach, while accounting for the underlying Thomson scattering of image formation. The model is optimized by iteratively fitting the time-dependent observational data. In addition, we apply physical constraints in terms of continuity, propagation direction, and speed of the heliospheric plasma to overcome limitations imposed by the sparse number of viewpoints. We utilize synthetic observations of a CME simulation to fully quantify the model's performance for different viewpoint configurations. The results demonstrate that our method can reliably estimate the CME parameters from only two viewpoints, with a mean velocity error of $3.01\pm1.94\%$ and propagation direction errors of $3.39\pm1.94^\circ$ in latitude and $1.76\pm0.79^\circ$ in longitude. We further show that our approach can achieve a full 3D reconstruction of the simulated CME from two viewpoints, where we correctly model the three-part structure, deformed CME front, and internal plasma variations. Additional viewpoints can be seamlessly integrated, directly enhancing the reconstruction of the plasma distribution in the heliosphere. This study underscores the value of physics-informed methods for reconstructing the heliospheric plasma distribution, paving the way for unraveling the dynamic 3D structure of CMEs and enabling advanced space weather monitoring.

The activity of Sun-like stars is governed by the magnetic field, which is believed to be generated in a thin layer between convective and radiative envelopes. The dynamo layer, also called the tachocline, permits the existence of Rossby waves (r-modes) described by magnetohydrodynamic shallow water models, which may lead to short-term cycles in stellar activity. Convective cells penetrate into the layer creating an overshoot upper part, where they transport an additional energy for vigorous activity. The aim of this paper is to study the influence of overshooting convection on the dynamics of Rossby waves in the tachoclines of Sun-like stars. Here we write the magnetohydrodynamic shallow water equations with the effect of the penetrative convection and study the dynamics of wave modes in the layer. The formalism leads to the excitation of new oscillation modes connected with the dynamo coefficient, alpha, causing periodic modulations of all parameters in the tachocline. The modes are coupled with the Rossby waves resulting mutual exchange of convective and rotation energies. The timescales of Rossby-dynamo waves, for certain parameters, correspond to Schwabe (11 years) and Rieger (150-170 days) cycles as observed in solar activity. The waves provide a new paradigm for internal magnetism and may drive the dynamos of Sun-like stars. Theoretical properties of the waves and observations can be used for magneto-seismological sounding of stellar interiors.

Using the horizontal neutral wind observations from the MIGHTI instrument onboard NASA's ICON (Ionospheric Connection Explorer) spacecraft with continuous coverage, we determine the climatology of the mean zonal and meridional winds and the associated mean circulation at low- to middle latitudes ($10^\circ$S-40$^{\circ}$N) for Northern Hemisphere {summer} solstice conditions between 90 km and 200 km altitudes, specifically on 20 June 2020 solstice as well as for a one-month period from 8 June-7 July 2020 {and for Northern winter season from 16 December 2019-31 January 2020, which spans a 47-day period, providing full local time coverage}. The data are averaged within appropriate altitude, longitude, latitude, solar zenith angle, and local time bins to produce mean wind distributions. The geographical distributions and local time variations of the mean horizontal circulation are evaluated. The instantaneous horizontal winds exhibit a significant degree of spatiotemporal variability often exceeding $\pm 150$ m s$^{-1}$. The daily averaged zonal mean winds demonstrate day-to-day variability. Eastward zonal winds and northward (winter-to-summer) meridional winds are prevalent in the lower thermosphere, which provides indirect observational evidence of the eastward momentum deposition by small-scale gravity waves. The mean neutral winds and circulation exhibit smaller scale structures in the lower thermosphere (90-120 km), while they are more homogeneous in the upper thermosphere, indicating the increasingly dissipative nature of the thermosphere. The mean wind and circulation patterns inferred from ICON/MIGHTI measurements can be used to constrain and validate general circulation models, as well as input for numerical wave models.

In this work, we study the accuracy that can be achieved when inferring the atmospheric information from realistic numerical magneto-hydrodynamic simulations that reproduce the spatial resolution we will obtain with future observations made by the 4m class telescopes DKIST and EST. We first study multiple inversion configurations using the SIR code and the Fe I transitions at 630 nm until we obtain minor differences between the input and the inferred atmosphere in a wide range of heights. Also, we examine how the inversion accuracy depends on the noise level of the Stokes profiles. The results indicate that when the majority of the inverted pixels come from strongly magnetised areas, there are almost no restrictions in terms of the noise, obtaining good results for noise amplitudes up to 1$\times10^{-3}$ of $I_c$. At the same time, the situation is different for observations where the dominant magnetic structures are weak, and noise restraints are more demanding. Moreover, we find that the accuracy of the fits is almost the same as that obtained without noise when the noise levels are on the order of 1$\times10^{-4}$of $I_c$. We, therefore, advise aiming for noise values on the order of or lower than 5$\times10^{-4}$ of $I_c$ if observers seek reliable interpretations of the results for the magnetic field vector reliably. We expect those noise levels to be achievable by next-generation 4m class telescopes thanks to an optimised polarisation calibration and the large collecting area of the primary mirror.

We performed a detailed analysis of the magnetic structure and asymmetric eruption of a large (about 500 Mm) inverse S-shaped filament partially located in AR 13229 on February 24, 2023. We linked the filament's pre-eruptive magnetic configuration to its highly asymmetric eruption dynamics and the formation of a large-scale coronal dimming in a weak-field region (mean unsigned flux of about 5 G). To reconstruct the coronal magnetic field, we applied a physics-informed neural network (PINN)-based nonlinear force-free field (NLFFF) extrapolation method to a pre-eruption HMI vector magnetogram. The NLFFF extrapolation reveals a large-scale magnetic flux rope (MFR) of about 500 Mm in length, consistent with the filament. We identified an extended MFR footprint to the east that connects to the J-shaped flare ribbon, outlining where the coronal dimming began. Overlying strapping fields connect to the area into which the dimming and flare ribbon later expand. This configuration explains the formation of the dimming as a stationary flux rope and strapping flux dimming, with subsequent expansion driven by the growth of the MFR footprint through strapping-strapping reconnection. Conversely, the western filament leg shows multiple anchor points and strong overlying magnetic fields, which suppressed the dimming and partially confined the eruption on that side. The reconstructed pre-eruptive NLFFF configuration offers a clear physical explanation for the asymmetries seen in the eruption, flare geometry, and coronal dimming. This demonstrates that PINN-based NLFFF extrapolation can effectively model large-scale filaments extending into weak-field regions, enhancing our understanding of complex solar eruptions.

The NSF's Astronomy and Astrophysics Postdoctoral Fellowship (AAPF) is exceptional among the available postdoctoral awards in Astronomy and Astrophysics. The fellowship is one of the few that allows postdoctoral researchers to pursue an original research program, of their own design, at the U.S. institution of their choice. However, what makes this fellowship truly unique is the ability of Fellows to lead an equally challenging, original educational program simultaneously. The legacy of this singular fellowship has been to encourage and advance leaders in the field who are equally as passionate about their own research as they are about sharing that research and their passion for astronomy with students and the public. In this positional paper we address the importance of fellowships like the AAPF to the astronomical profession by identifying the science and educational contributions that Fellows have made to the community. Further, we recommend that fellowships that encourage leading postdoctoral researchers to also become leaders in Astronomy education be continued and expanded.

We report the distribution of planets as a function of planet radius (R_p), orbital period (P), and stellar effective temperature (Teff) for P < 50 day orbits around GK stars. These results are based on the 1,235 planets (formally "planet candidates") from the Kepler mission that include a nearly complete set of detected planets as small as 2 Earth radii (Re). For each of the 156,000 target stars we assess the detectability of planets as a function of R_p and P. We also correct for the geometric probability of transit, R*/a. We consider first stars within the "solar subset" having Teff = 4100-6100 K, logg = 4.0-4.9, and Kepler magnitude Kp < 15 mag. We include only those stars having noise low enough to permit detection of planets down to 2 Re. We count planets in small domains of R_p and P and divide by the included target stars to calculate planet occurrence in each domain. Occurrence of planets varies by more than three orders of magnitude and increases substantially down to the smallest radius (2 Re) and out to the longest orbital period (50 days, ~0.25 AU) in our study. For P < 50 days, the radius distribution is given by a power law, df/dlogR= k R^\alpha. This rapid increase in planet occurrence with decreasing planet size agrees with core-accretion, but disagrees with population synthesis models. We fit occurrence as a function of P to a power law model with an exponential cutoff below a critical period P_0. For smaller planets, P_0 has larger values, suggesting that the "parking distance" for migrating planets moves outward with decreasing planet size. We also measured planet occurrence over Teff = 3600-7100 K, spanning M0 to F2 dwarfs. The occurrence of 2-4 Re planets in the Kepler field increases with decreasing Teff, making these small planets seven times more abundant around cool stars than the hottest stars in our sample. [abridged]

In this paper we address the consumption of energy by astronomers while performing their professional duties. Although we find that astronomy uses a negligible fraction of the US energy budget, the rate at which energy is consumed by an average astronomer is similar to that of a typical high-flying businessperson. We review some of the ways in which astronomers are already acting to reduce their energy consumption. In the coming decades, all citizens will have to reduce their energy consumption to conserve fossil fuel reserves and to help avert a potentially catastrophic change in the Earth's climate. The challenges are the same for astronomers as they are for everyone: decreasing the distances we travel and investing in energy-efficient infrastructure. The high profile of astronomy in the media, and the great public interest in our field, can play a role in promoting energy-awareness to the wider population. Our specific recommendations are therefore to 1) reduce travel when possible, through efficient meeting organization, and by investing in high-bandwidth video conference facilities and virtual-world software, 2) create energy-efficient observatories, computing centers and workplaces, powered by sustainable energy resources, and 3) actively publicize these pursuits.

With the stagnation of processor core performance, further reductions in the time-to-solution for geophysical fluid problems are becoming increasingly difficult with standard time integrators. Parallel-in-time exposes and exploits additional parallelism in the time dimension which is inherently sequential in traditional methods. The rational approximation of exponential integrators (REXI) method allows taking arbitrarily long time steps based on a sum over a number of decoupled complex PDEs that can be solved independently massively parallel. Hence REXI is assumed to be well suited for modern massively parallel super computers which are currently trending. To date the study and development of the REXI approach has been limited to linearized problems on the periodic 2D plane. This work extends the REXI time stepping method to the linear shallow-water equations (SWE) on the rotating sphere, thus moving the method one step closer to solving fully nonlinear fluid problems of geophysical interest on the sphere. The rotating sphere poses particular challenges for finding an efficient solver due to the zonal dependence of the Coriolis term. Here we present an efficient REXI solver based on spherical harmonics, showing the results of: a geostrophic balance test, a comparison with alternative time stepping methods, an analysis of dispersion relations, indicating superior properties of REXI, and finally a performance comparison on Cheyenne supercomputer. Our results indicate that REXI is not only able to take larger time steps, but that REXI can also be used to gain higher accuracy and significantly reduced time-to-solution compared to currently existing time stepping methods.

We present a parametric space study of the decay of turbulence in rotating flows combining direct numerical simulations, large eddy simulations, and phenomenological theory. Several cases are considered: (1) the effect of varying the characteristic scale of the initial conditions when compared with the size of the box, to mimic "bounded" and "unbounded" flows; (2) the effect of helicity (correlation between the velocity and vorticity); (3) the effect of Rossby and Reynolds numbers; and (4) the effect of anisotropy in the initial conditions. Initial conditions include the Taylor-Green vortex, the Arn'old-Beltrami-Childress flow, and random flows with large-scale energy spectrum proportional to $k^4$. The decay laws obtained in the simulations for the energy, helicity, and enstrophy in each case can be explained with phenomenological arguments that separate the decay of two-dimensional from three-dimensional modes, and that take into account the role of helicity and rotation in slowing down the energy decay. The time evolution of the energy spectrum and development of anisotropies in the simulations are also discussed. Finally, the effect of rotation and helicity in the skewness and kurtosis of the flow is considered.

We analyze recent high resolution photospheric small-scale dynamo simulations that were computed with the MURaM radiative MHD code. We focus the analysis on newly forming downflow lanes in exploding granules since they show how weakly magnetized regions in the photosphere (center of granules) evolve into strongly magnetized regions (downflow lanes). We find that newly formed downflow lanes exhibit initially mostly a laminar converging flow that amplifies the vertical magnetic field embedded in the granule from a few 10 G to field strengths exceeding 800 G. This results in extended magnetic sheets that have a length comparable to granular scales. Field amplification by turbulent shear happens first a few 100 km beneath the visible layers of the photosphere. Shallow recirculation transports the resulting turbulent field into the photosphere within minutes, after which the newly formed downflow lane shows a mix of strong magnetic sheets and turbulent field components. We stress in particular the role of shallow and deep recirculation for the organization and strength of magnetic field in the photosphere and discuss the photospheric and sub-photospheric energy conversion associated with the small-scale dynamo process. While the energy conversion through the Lorentz force depends only weakly on the saturation field strength (and therefore deep or shallow recirculation), it is strongly dependent on the magnetic Prandtl number. We discuss the potential of these findings for further constraining small-scale dynamo models through high resolution observations.

Over the past decades and even centuries, the astronomical community has accumulated a signif-icant heritage of recorded observations of a great many astronomical objects. Those records con-tain irreplaceable information about long-term evolutionary and non-evolutionary changes in our Universe, and their preservation and digitization is vital. Unfortunately, most of those data risk becoming degraded and thence totally lost. We hereby call upon the astronomical community and US funding agencies to recognize the gravity of the situation, and to commit to an interna-tional preservation and digitization efforts through comprehensive long-term planning supported by adequate resources, prioritizing where the expected scientific gains, vulnerability of the origi-nals and availability of relevant infrastructure so dictates. The importance and urgency of this issue has been recognized recently by General Assembly XXX of the International Astronomical Union (IAU) in its Resolution B3: "on preservation, digitization and scientific exploration of his-torical astronomical data". We outline the rationale of this promotion, provide examples of new science through successful recovery efforts, and review the potential losses to science if nothing it done.

We report on the masses, sizes, and orbits of the planets orbiting 22 Kepler stars. There are 49 planet candidates around these stars, including 42 detected through transits and 7 revealed by precise Doppler measurements of the host stars. Based on an analysis of the Kepler brightness measurements, along with high-resolution imaging and spectroscopy, Doppler spectroscopy, and (for 11 stars) asteroseismology, we establish low false-positive probabilities for all of the transiting planets (41 of 42 have a false-positive probability under 1%), and we constrain their sizes and masses. Most of the transiting planets are smaller than 3X the size of Earth. For 16 planets, the Doppler signal was securely detected, providing a direct measurement of the planet's mass. For the other 26 planets we provide either marginal mass measurements or upper limits to their masses and densities; in many cases we can rule out a rocky composition. We identify 6 planets with densities above 5 g/cc, suggesting a mostly rocky interior for them. Indeed, the only planets that are compatible with a purely rocky composition are smaller than ~2 R_earth. Larger planets evidently contain a larger fraction of low-density material (H, He, and H2O).

The efficient release of magnetic energy in astrophysical plasmas, such as during solar flares, can in principle be achieved through magnetic diffusion, at a rate determined by the associated electric field. However, attempts at measuring electric fields in the solar atmosphere are scarce, and none exist for sites where the magnetic energy is presumably released. Here, we present observations of an energetic event using the National Science Foundation's Daniel K. Inouye Solar Telescope, where we detect the polarization signature of electric fields associated with magnetic diffusion. We measure the linear and circular polarization across the hydrogen H-epsilon Balmer line at 397 nm at the site of a brightening event in the solar chromosphere. Our spectro-polarimetric modeling demonstrates that the observed polarization signals can only be explained by the presence of electric fields, providing conclusive evidence of magnetic diffusion, and opening a new window for the quantitative study of this mechanism in space plasmas.