高能暂现源 · 每日文献综述

M-dwarfs produce frequent flares, and their associated mass ejections are expected to significantly affect the habitability of close-in exoplanets. Recent spectroscopic observations have revealed several prominence eruptions-indicative of stellar mass ejections-on M-dwarfs through Doppler shifts of the H$α$ line. However, systematic and statistical studies, particularly regarding their association with white-light flares, have been limited due to the lack of intensive and continuous simultaneous photometric and spectroscopic monitoring of the same target star. We conducted one month of continuous spectroscopic observations of the active M-dwarf YZ CMi using the 3.8-m Seimei telescope with an unprecedentedly high time cadence of $\sim$1 min, simultaneously with TESS. We detected four prominence eruptions, among which two events showed rapid, short-duration eruptions with velocities of $\sim$300-500 km s$^{-1}$ and durations of $\sim$5 min. Such short-duration events may have been missed in previous observations due to insufficient time cadence. We further performed a systematic analysis using 35 H$α$ flares on YZ CMi observed simultaneously with TESS. Notably, most prominence eruptions (6 out of 7) were not associated with detectable white-light flares. This result suggests that most observed prominence eruptions on M-dwarfs may have occurred near the stellar limb, where white-light flares are difficult to detect, which may imply a potential observational bias in their detectability due to low contrast with the background emission. These first statistical constraints, together with the discovery of rapid, short-duration prominence eruptions, indicate that previous observations may have underestimated both the frequency and velocities of mass ejections on M-dwarfs due to observational biases and highlight the necessity of reassessing their impact on close-in exoplanets.

Using pulse profile modeling, we obtain the mass-radius measurement of a millisecond pulsar (MSP) with data from the Neutron Star Interior Composition ExploreR, XMM-Newton and the Chandra X-ray Observatory. We report here the radius of PSR J1614$-$2230, the second most massive MSP confirmed by radio timing. All of the data sets are well described by a simple model composed of two circular hot spots. The final result yields an equatorial radius of $R_{\rm eq}=10.06 ^{+1.25}_{-0.87}\,$km, and a gravitational mass of $M=1.937^{+0.012}_{-0.013}\,M_{\odot}$ (equally tailed 68% credible intervals). Although a non-thermal component was previously reported at higher energies, we find no sign of it in either our phase-averaged or phase-resolved spectral analyses. Using new relations linking the compactness to oblateness or surface gravity, and tailored to the spin frequency of PSR J1614$-$2230, we infer a configuration with one hot spot near the pole, and another near the equator. The tight mass posterior is essentially informed by radio timing, while the radius constraint is not as tight due to the low source signal (8.5$σ$ X-ray pulse significance). However, over all geometries and atmosphere models tested, the radius posterior tends toward low values ($\lesssim12.08\,$km, 90th percentile in all cases).

We numerically study the collapse of rotating, very-massive stellar cores with masses of $\approx 200$, $300$, $500$, and $1100M_\odot$ into black holes using both axisymmetric and three-dimensional (3D) numerical relativity. Our results indicate that when the dimensionless spin of the resulting black hole exceeds 0.8, a massive disk consistently forms around it. These massive, compact disks, carrying more than about 10\% of the black hole's mass, are prone to non-axisymmetric deformations that trigger gravitational-wave bursts with frequencies around 10-50 Hz. Such waves could be detected by the Einstein Telescope and Cosmic Explorer, even from sources a few Gpc away. We also summarize the gravitational-wave signals from axisymmetric collapse, which tend to have lower amplitudes and higher frequencies than those caused by non-axisymmetric instabilities.

With the projected increase in the detection rate of compact-binary coalescences in the coming decade, there is critical need to develop fast, robust, and scalable alternatives to matched filtering for gravitational-wave searches. Transformer models have revolutionized natural language and audio processing but their application to gravitational-wave astronomy is still largely unexplored. In this work, we introduce \castor, a transformer-based coincident search pipeline for detecting binary black hole gravitational-wave signals from Advanced LIGO detectors. One of the major features of our model is that it allows the false-alarm rate to be estimated via time slides cheaply without requiring repeated evaluations of the neural network. We evaluate \castor\ on datasets from the Machine-Learning Gravitational-Wave Search Challenge (MLGWSC-1) and on approximately five months of real O3b observing strain. When tested on benchmark datasets, \castor\ ranks among the most sensitive machine-learning pipelines and successfully recovers the majority of confident events from the GWTC-3 catalog that lie within its training range. We also benchmark \castor\ against another transformer architecture, GW-Whisper, a domain-adaptation of OpenAI's audio foundation model. We find that \castor\ substantially outperforms the repurposed audio model in sensitivity and also reduces the computational cost of background estimation by a factor of 20. Our results demonstrate a highly practical, scalable approach for deep-learning gravitational wave searches and empirical background estimation for future observing runs.

We present JWST/NIRSpec IFU observations of SPIRITS 19q, the highly dust-producing carbon-rich (WC) binary candidate located in a subsolar-metallicity region of the nearby spiral galaxy NGC 2403. The observations, taken in April of 2024, confirm the association of a dusty outburst observed in 2019 by the Spitzer Space Telescope with an early-type WC star. Using models from the Potsdam Wolf-Rayet (PoWR) LMC model grid we find that the WC star of SPIRITS 19q likely has an especially high mass-loss rate ($\gtrsim$ 10$^{-4}$ $M_{\odot}$ yr$^{-1}$). From the flux peak of the IR transient as measured by Spitzer/IRAC as well as constraints on dust composition and size from the JWST spectrum, we estimate a total dust mass formed in the outburst of 6.6 $\pm$ 0.4 $\times$ 10$^{-6}$ $M_{\odot}$. Assuming a minimum orbital period of 12 years, this corresponds to a period-averaged dust production rate of $\lesssim$ 5.5 $\times$ 10$^{-7}$ $M_{\odot}$ yr$^{-1}$. These observations suggest that even a single WC system can contribute to the dust budget at subsolar metallicities, and that such systems are an important source of carbonaceous dust grains in the early universe.

It has been proposed that the radio-weak Seyfert galaxy GRS 1734-292 is the counterpart of 3EG J1736-2908 and of the nearby Fermi/LAT source. I examine this identification, considering that the radio weakness could result from free-free absorption of a relativistic jet. An association with the EGRET detections is plausible, whereas the latest Fermi/LAT localization no longer favors GRS 1734-292. New observations, especially at radio frequencies of tens of GHz, would help to settle the issue.

The coefficients of the transport equations, describing the propagation features, are particularly important in astrophysics, interplanetary physics, and experimental plasma physics. In this paper, the variable-dependent diffusion coefficients are investigated. For the momentum transport equations, we find that the coefficients are related to statistical quantities, some of which are new. In addition, these coefficients take logarithmic forms, which are different from the results derived in previous papers. For isotropic pitch-angle scattering, we also obtain a formula that takes a logarithmic form. For fractional transport equations, the coefficients expressed in terms of the moments are obtained from the governing equations.

Neutrinos produced at the highest energies can interact with cosmic neutrino and radiation backgrounds during their propagation to Earth. The many available $νν$, $ν\barν$, and $νγ$ channels can lead to their absorption or energy redistribution, whilst the leptonic and hadronic final states may feed secondary fluxes of neutrinos, protons, and electromagnetic particles through the decay or hadronisation of the heavy leptons, bosons, and quarks produced. We present a framework to characterise these propagation effects in detail, \texttt{$ν$propa}, an extension of the CRPropa Monte Carlo code that interfaces with event generators and to dedicated computations of the relevant cross sections. It also treats flavour oscillations in vacuum. Using this code, we investigate sources at high redshifts ($z = 10$), and find a strong absorption of the prompt flux beyond~$\sim 10^{21} \; \text{eV}$, although the copious secondary neutrinos partially compensate this depletion, also contributing to the spectrum at lower energies. The framework is designed to study scenarios of cosmological neutrino production beyond~EeV energies such as superheavy dark matter, cosmic strings, and primordial black holes, and to yield reliable predictions for the forthcoming neutrino observatories.

Recently, the Large High Altitude Air Shower Observatory (LHAASO) reported the first TeV gamma-ray detection from a low-luminosity active galactic nucleus (LLAGN) lacking a strong jet. This source, NGC 4278, is one of the nearest LLAGNs that can be studied in detail as a template for more distant objects. Most previous studies have focused on the month-scale-variable gamma-ray activities, seeking to explain them with weak jets. We focus instead on NGC 4278's quiet-period emission, proposing a different mechanism based on a wind launched from a radiatively inefficient accretion flow (RIAF). In particular, we model particle acceleration at a wind-driven shock and calculate the resulting multimessenger emission. To match the TeV gamma-ray fluxes reported by LHAASO in the quiet period, our model requires a wind power of $\lesssim2\times10^{43}~\mathrm{erg~s^{-1}}$ with a spectral index of injected cosmic rays of 2.0. We show that this is achievable through Bondi accretion of the ionized gas observed on the 100-pc scale, for a shallow accretion profile. We discuss the testability of our wind-driven shock scenario with multimessenger observations. As another important application, we extend our model to Sgr A*, the nearest LLAGN, and show that its TeV gamma-ray emission is reproduced with reasonable parameters. Our results suggest that RIAF-driven winds may be a generic emission mechanism in LLAGNs, making this class a potentially important population of both gamma-ray emitters and cosmic-ray accelerators.

Three-minute ($\sim$5 mHz) oscillations are a prominent signature of wave propagation in the solar chromosphere, but their measured power can vary substantially between full-disc observations and small, high-resolution fields of view. We investigate whether chromospheric magnetic canopy topology, particularly its lateral continuity and areal filling, modulates field-of-view-averaged three-minute power. We analyse H$α$ and Ca II 8542 Å spectral imaging from the Swedish 1-m Solar Telescope, comparing canopy-poor quiet-Sun scenes with canopy-dominated active-region environments, complemented by radiative-MHD Bifrost simulations of canopy-poor and canopy-rich atmospheres. We compute field-of-view-averaged power spectra from fixed-wavelength intensities and bisector-derived velocities in the observations, and from synthetic intensities, velocities, and temperature proxies in the simulations. We also quantify the relation between integrated 3-5 mHz excess power and canopy filling factor, using HMI-based magnetic-field extrapolations for the observations. In both observations and simulations, canopy-poor scenes show clear 3-5 mHz power enhancements, whereas canopy-dominated scenes show strongly reduced measured power. The trend is present in H$α$ and Ca II 8542 Å intensities and is supported by complementary velocity and temperature proxies. Integrated 3-5 mHz excess power decreases with increasing canopy filling factor. We conclude that chromospheric field-of-view-averaged three-minute power is a topology-weighted observable: dense, laterally continuous magnetic canopies can strongly reduce the measured large-scale three-minute signature without requiring the absence of local wave activity. This provides a practical criterion for interpreting chromospheric oscillation power as a proxy for atmospheric wave energy in the Sun and in magnetised cool-star atmospheres.

The design of the IceCube-Gen2 observatory, proposed as a next-generation extension of IceCube, includes a surface array consisting of scintillators and radio antennas. In addition to several such detectors already deployed at IceCube's surface array at the South Pole, a complete prototype station including three SKALA antennas has been operating at the Pierre Auger Observatory for several years. This setup has been used to successfully observe radio signals from air showers, demonstrated through coincident observations with the Auger Surface Detector. In this contribution, we present an updated analysis of these radio signals, and compare them to CoREAS simulations using the reconstruction from the Auger Surface Detector as input.

Layering is a ubiquitous feature of astrophysical objects. Motivated by the fact that physical black holes retain the layered structure of their progenitor stars when viewed in the time concepts synchronizable with the clock of an outside fixed-position probe, we investigate linear perturbations and gravitational-wave (GW) echoes from a compact object composed of two concentric thin shells. Compared with the single-shell case, the double-shell structure introduces an extra barrier in the effective potential and partitions the wave propagation space into four coupled effective cavities. As the mass ratio of the inner shell increases, new spectral peaks enter from the high-frequency side of the echo spectrum; the second and later peaks shift toward higher frequencies; and the lowest-frequency peak first shifts toward lower frequencies and then returns to its $q=0$ position. We call this variation pattern spectral-peak queueing (SQ). Its existence suggests that GW echoes can be used as probes for the internal structure of compact objects under consideration.

In this work, we investigate the resonant interaction between stochastic gravitational waves (GWs) and chiral Alfvén waves in a magnetized chiral plasma in the early Universe. Starting from covariant chiral magnetohydrodynamics coupled to linearized gravity, we derive a closed system of equations for the coupled chiral-Alfvén velocity and magnetic-field perturbations. Both analytics and numerics show a parametric resonance at the sum frequency of the two CVE-split branches, with an instability band that widens as the GW strain grows. This CVE-to-Alfvén ratio is not a free parameter. In the early Universe, it is fixed by the Standard Model's relativistic degrees of freedom. Including the plasma's backreaction on the GW makes the energy exchange nonlinear, breaking the usual single-frequency Manley-Rowe picture. For a self-consistent choice of parameters, where backreaction is only a small correction, the unsuppressed resonance still drives the system to a finite-time blow-up within a few Hubble times. We then work out the resonant frequencies and growth rates involved, mapping out what current and future gravitational-wave detectors could observe.

The Askaryan Radio Array (ARA) is an ultrahigh energy (UHE) neutrino detector at the South Pole that searches for impulsive broadband radio signals from neutrino-induced particle showers in glacial ice. ARA consists of five autonomous stations with receiving antennas deployed up to 200 m deep and has accumulated the largest livetime of any in-ice radio array. We present the first array-wide diffuse UHE neutrino search using data collected from January 2013 to December 2023, incorporating improved detector characterization and simulation, including data-driven noise and electronics models, revised antenna responses, and updated neutrino interaction and lepton propagation modeling, within a unified framework for event processing, background rejection, and cut optimization across independently operating stations. This search is expected either to identify the first UHE neutrino candidates observed by an in-ice radio detector or to set the most stringent UHE neutrino diffuse flux limits above a few EeV, while informing analysis strategies for under-construction and future radio arrays, such as RNO-G and IceCube-Gen2 Radio.

Black holes link gravity, thermodynamics and information via limits on the relaxation rate of a perturbed system. The Bekenstein-Hod bound imposes a minimum relaxation time at fixed temperature, but its observational test demands both black hole thermodynamic characterization and decay time measurement. Here we test this bound with GW250114, the loudest gravitational-wave signal yet observed from a binary black-hole merger. We infer the remnant temperature from pre-merger data truncated at least $10\,M$ before the peak and its longest-lived decay time from post-merger data, thereby avoiding direct reuse of the same strain samples. The ringdown frequencies and damping times are allowed to vary independently rather than being fixed to the Kerr spectrum. For the primary $t_{<}=-10\,M$ analysis, the bound is verified at $3.3-3.6σ$ across the focal ringdown start times, representing a substantial improvement over the $91\%$ confidence level set by GW150914. The conclusion remains robust under varied pre-merger cutoffs and explicit inclusion of the short-lived first overtone in waveform modelling. This separated-data measurement converts an information-theoretic relaxation bound into a precision test of a single astrophysical black hole.

For eccentric systems with aligned spins, we compute all spin contributions to the gauge-invariant conservative map between the constants of motion (energy and angular momentum) and the fundamental (radial and azimuthal) frequencies at the fourth post-Newtonian order. This extends the nonspinning map derived in [arXiv:2511.10735] to include terms that are linear, quadratic, cubic, and quartic in spin, as well as the associated spin-deformability parameters. We also obtain derived quantities such as the redshift and gyroscopic invariants, as well as circular links. We notice that the Blanchet-Iyer-Favata post-Newtonian stability criterion that determines the dimensionless frequency of the innermost stable circular orbit (ISCO) is in fact intricately related to the periastron advance for circular orbits. Finally, in the case of unbound orbits, we complete the 4PN scattering angle with all spin contributions.

Astrophysical observations of our universe have been key to our understanding of how the universe works. Shortly after the turn of the millennium, the National Research Council delivered \textit{Connecting Quarks with the Cosmos: Eleven Science Questions for the New Century}. In the subsequent quarter-century, we have made substantial progress in answering each question. These advancements have, in part, arisen because of the success of major US facilities across several domains of physics, guided by long-term planning documents which still largely focus on these questions. This report seeks to provide a status update on each question, and to outline what space-based facilities are crucial for future progress, intended to guide NASA's preparatory work for the Astro2030 Decadal.

Nambu-Goto cosmic string loops generically develop cusps, points where the string momentarily reaches the speed of light. These cusps produce strong gravitational wave bursts with a characteristic strain spectrum $\tilde{h}(ω)\propto (Gμ)\,ω^{-4/3}$, making them prime targets for gravitational wave searches, where $μ$ is the string mass per unit length. However, this picture is modified when one accounts for gravitational backreaction. Using a convenient gauge, we reformulate the Nambu-Goto equations of motion for a loop moving in its own dynamically sourced gravitational field, enabling the first continuous numerical evolution of loops under this backreaction. A key finding is that locally cusps survive backreaction. Nevertheless, the gravitational waveform calculated in this formalism is significantly modified, weakening the cusp burst and introducing a high-frequency cutoff at $f_c\propto(Gμ)^{-3/2}$. This suppression above $f_c$ reduces the expected signal-to-noise ratio in current and near-future detectors relative to unperturbed waveform predictions.

Supermassive black holes/active galactic nuclei (SMBHs/AGN), forming only a few hundred million years after the Big Bang as observed with the James Webb Space Telescope (JWST), challenge theoretical understanding. How could they grow so massive $(M_{\rm BH} > 10^6 {\rm \,M}_\odot)$ so quickly after initial seeding? Is this rapid growth related to the numerous and enigmatic Little Red Dots (LRDs), compact sources with AGN-like characteristics, discovered by JWST? To address these mysteries, we consider the first-order constraint on SMBH growth: enough baryonic material has to reach the vicinity of the SMBH seed, located near the bottom of the gravitational potential well of the host dark matter halo. We specifically examine cold-mode accretion, where gas from the cosmic environment flows into the virialized halo in cold streams without being shock-heated, efficiently reaching the center on a free-fall timescale. We find that cold mode accretion is necessary to supply material for the SMBHs to reach the observed masses, whereas for shock-heated gas inflow the required amount could only be supplied by implausibly rare halos. Moreover, cold-mode inflow in rare $(\sim1$ Gpc$^{-3}$) halos matches the mass and number of the massive quasars, and halos able to support super-Eddington accretion for massive SMBHs ($\sim10^7$ M$_\odot$) match LRD number densities. The decreasing LRD abundance at lower redshifts may then reflect the termination of cold-mode accretion in the growing host halos. The populations of massive SMBHs and LRDs at early times may thus arise naturally from cosmological structure formation, based on the abundance of halos capable of supplying sufficient material through cold accretion.

Observations of magnetic polluted white dwarfs indicate that most of them have higher concentrations of metals near the poles compared to lower latitudes. Maintaining such abundance gradients requires gravitational sinking times $t_\downarrow$ below the convection zone that are shorter than the horizontal spreading time $t_\leftrightarrow$ across the surface by convective eddies. We show analytically that $t_\downarrow/t_\leftrightarrow\propto T^{8/3}$, where $T$ is the temperature at the base of the convection zone, which rises by more than an order of magnitude as convection penetrates deeper into the atmosphere. Correspondingly, $t_\downarrow/t_\leftrightarrow$ jumps by several orders of magnitude, clearly delineating between hot white dwarfs with abundance variations and cold ones with homogenous surfaces. We incorporate magnetic fields self-consistently into the stellar structure and compute for the first time the sinking to spreading time-scale ratio as a function of $B$. Magnetic white dwarfs have shallower convection zones at a given $T_{\rm eff}$, and the inward penetration of convection shifts to lower $T_{\rm eff}$. Quantitatively, $B\sim 10^5\textrm{ G}$ extends the range for enhanced abundance patches from $T_{\rm eff}\gtrsim 30\,000 \textrm{ K}$ ($13\,000\textrm{ K}$) to $T_{\rm eff}\gtrsim 15\,000 \textrm{ K}$ ($6000\textrm{ K}$) for helium (hydrogen) dominated atmospheres, which is insufficient to explain the observed patchy white dwarfs, which are even colder. A possible solution is variability in the accretion rate on a time-scale of a decade, as recently detected in another white dwarf.

The Sun has long served as a natural dark matter detector, capturing halo particles that scatter with solar nuclei and electrons and subsequently produce indirect signals from dark matter annihilation that can be probed by either neutrino or $γ$-ray instruments. Previous efforts have computed cross-section constraints for dark matter scattering with electrons, as well as spin-independent and spin-dependent dark matter-proton scattering. However, the spin-dependent dark matter-neutron scattering scenario has been neglected, despite its importance in direct detection. For the first time, we compute the capture and evaporation rates for spin-dependent dark matter-neutron scattering in the Sun, including odd-neutron isotope targets in the current standard solar model. Comparing the predicted annihilation signals against current neutrino and $γ$-ray observations, and projecting the reach of upcoming detectors, we find that solar constraints exceed direct detection limits for several annihilation channels. For models in which the dark matter annihilates into long-lived mediators, these constraints can extend below the neutrino fog.

Axion fields can form exponentially growing gravitational clouds around compact objects through self-interaction-driven relaxation of ambient axion waves. As the field amplitude approaches the axion decay constant, nonlinear effects become important. We identify two distinct regimes of late-time evolution, determined by the gravitational fine-structure constant and the cloud growth rate: a Bosenova regime, characterized by collapse accompanied by explosive axion bursts, and a saturation regime, in which self-interaction-induced axion emission balances accretion. In the latter regime, the emitted axion radiation exhibits stable discrete spectral lines at odd multiples of the bound-state energy, directly probing the global structure of the axion potential beyond its quadratic minimum. We show that single-cosine potentials and QCD axion-like potentials predict distinct emission spectra, enabling probes of the underlying axion self-interaction structure and its ultraviolet completion through terrestrial detection of relativistic axion fluxes from compact objects.

We present results from a deep MeerKAT 1.28 GHz observation of Abell 520, a galaxy cluster in the middle of a merger and a treasure trove of merger-generated phenomena. MeerKAT's combination of high angular resolution and sensitivity to extended radio emission lets us uncover unprecedented detail of the cluster's giant halo and puzzling features beyond the halo. We find that the main radio halo is sharply bounded along the direction of the NE-SW merger axis by the two X-ray shock fronts. In the perpendicular direction, the halo wings smoothly span over 2 Mpc, with the spectrum steepening in the outskirts. We detect faint radio emission upstream of the main (SW) bow shock; the emissivity jump across this shock can be used to constrain the physical mechanisms at work there. We discover a tail of radio emission extending past the NE shock to the gas-poor, high-entropy ``dark subcluster'' in the NE outskirts. The tail and subcluster have a higher radio/X-ray brightness ratio than the main halo's, suggesting efficient reacceleration of cosmic rays. We also discover two peripheral radio relics beyond the cluster virial radius, and a faint, $\sim2$ Mpc long, broad radio bridge connecting the dark subcluster to one of the relics along the direction perpendicular to the merger axis. The spectral slopes of the dark subcluster and the bridge are consistent with that of the main halo. This bridge and both relics may be caused by another giant merger-generated shock front at the virial radius.

Using first-principles field-theoretic methods, we investigate neutrino emission from strongly magnetized dense quark matter under conditions relevant to compact stars. We account for Landau-level quantization of both electron and quark states and show that it strongly modifies the kinematics of Urca processes. In particular, quark quantization restricts the available phase space in which both quark and electron energies can simultaneously lie near their respective Fermi surfaces. At moderately strong magnetic fields, before pronounced quark quantization sets in, the emission rate tends to increase on average with increasing field strength. In the regime of very strong fields, however, the increasingly restricted phase space first gives rise to Shubnikov--de Haas-type oscillations and then to resonance-like spikes near a discrete sequence of Urca-resonant magnetic field values, separated by regions of strong suppression. Finally, the emission rate becomes nearly completely quenched once $|eB| \gtrsim 6μ_{e}μ_{u}$, corresponding to approximately $B\gtrsim 1.5\times 10^{19}~\mbox{G}$ for the representative set of model parameters considered. We also find significant anisotropy in the longitudinal momentum emission near the Urca-resonant magnetic field values.

We study the impact of new mass-radius information from PSR J1614-2230, PSR J2124-3358, and 47 Tuc X7 in a combined equation of state inference based on chiral effective field theory constraints at nuclear densities and using different high-density extensions, including perturbative QCD constraints. The largest impact stems from the heavy-mass PSR J1614-2230 star, which shifts heavy neutron stars to smaller radii by around 0.4 km. Moreover, the combined astrophysical NICER, LIGO/Virgo, and X-ray information drives the radius posterior to a more data-driven distribution, which is less sensitive to the high-density extension. For the equation of state, the new mass-radius information significantly tightens the pressure and speed-of-sound posterior distributions, especially around three times saturation density. Finally, we make predictions for the poorly constrained masses of PSR J2124-3358 and 47 Tuc X7 based on the combined equation of state analysis and the other astrophysical sources.

We present the first measurement of the Hubble constant ($H_0$) from the quadruply-lensed quasar J1537$-$3010, which has optically-measured time delays at $\sim2 \%$ precision. We combine these delays with multi-band imaging data from the Hubble Space Telescope (HST) and model the system with two independent software and teams. We adopt a mass profile that is maximally degenerate with $H_0$ to fully incorporate the mass-sheet degeneracy in the error budget, with nuisance parameters constrained by spatially resolved stellar kinematics from the Multi Unit Spectroscopic Explorer (MUSE) and a line-of-sight (LoS) analysis using the Euclid Flagship simulation. The entire analysis is performed blindly to $H_0$, distances and mass density slope of the main deflector. After unblinding, we measure $H_0 = 75.5^{+9.3}_{-5.8}\ {\rm km\,s^{-1}\,Mpc^{-1}}$, corresponding to a $10\%$ precision measurement from a single system. This precision is driven by conservative lens modeling assumptions including differences between lens modeling methods and the limited stellar kinematics constraints, from which we infer a total mass-sheet parameter $λ\equiv(1-κ_{\rm ext})λ_{\rm int} = 0.89^{+0.11}_{-0.06}$ ($λ_{\rm int}=0.89^{+0.09}_{-0.07}$) that is consistent with current results for elliptical galaxies ($λ\approx 1$). Our lensing constraints, stellar kinematic measurements and LoS characterization will be included in subsequent population-level measurements of $H_0$. Moreover, our lens models combined with future near-infrared spectroscopy and imaging from the James Webb Space Telescope will further reduce the uncertainties on $H_0$ from J1537$-$3010 alone, bringing it closer to the few percents precision of the time delays.

Using the elemental abundances in 68 metal-poor (MP) stars, we present a data-driven model for $r$-process production patterns covering Sr to U. We show that essentially all the $r$-process patterns in those and other test stars can be adequately explained as mixtures of Patterns 1 and 2, which provides theoretical insights into the empirical categories of limited-$r$, $r$-I, and $r$-II stars. We carry out an extensive survey of the yield templates produced by parametric $r$-process calculations. We propose that Pattern 2 may be dominated by a single template and points to regularity in $r$-process production by a subset of neutron star mergers (NSMs), while Pattern 1 is the average superposition of multiple templates and reflects production by other NSMs and perhaps also some magneto-rotational supernovae. We raise possible systematic issues with abundance ratios for elements measured in different ionization states, and highlight the need for examining the Os, Ir, and Pt measurements for HD~122563, which is dominated by Pattern 1 with a prominent Pt peak in our model. If this result is confirmed, the meaning of the limited-$r$ category requires drastic revision. Further measurements of a wider range of $r$-process elements in a larger sample of MP stars are critical to test and improve our model.

The early evolution of fast polar coronal hole (PCH) solar wind remains largely unconstrained by in situ measurements. In March 2025, Parker Solar Probe (Parker) at its closest approach of 9.86 Solar Radii ($R_\odot$) measured outflow from a large equatorial coronal hole (ECH) which was also measured at 1\,au and at intermediate distances by Solar Orbiter (also near its perihelion). At 1\,au the stream properties are consistent with PCH properties established by Ulysses. The stream was measured by Parker substantially below the Alfvén surface, with proton temperatures in excess of 2\,MK and a speed at $\sim$10\,$R_\odot$ which was only $\sim$60\% of its asymptotic value. The Solar Orbiter data indicates that the acceleration is largely complete by 60~$R_{\odot}$. Spherically-polarized fluctuations in the stream are observed to develop from near-transverse and small-angle at Parker to full reversal ``switchbacks'' at Solar Orbiter. Comparison of the implied acceleration profile to historical doppler-dimming measurements suggests that the stream's low coronal acceleration is similar to that of PCH flows. Consistent with previous work, this acceleration requires significantly more energy than can be provided by the observed thermal pressure gradients, with a significant contribution likely from the abundant Alfvénic fluctuation energy observed at Parker. These observations provide unique constraints on models of the radial evolution of the fastest solar wind, and indicate that these wind streams experience gradual, steady acceleration over their first few tens of solar radii of evolution.

We present a Geant4 Monte Carlo simulation of cosmic-ray (CR) hadronic interactions in a spherical molecular cloud (<inline-formula><mml:math><mml:mo>~</mml:mo><mml:mn>5400</mml:mn><mml:msub><mml:mi>M</mml:mi><mml:mo>⊙</mml:mo></mml:msub></mml:math></inline-formula>, 13 concentric shells, <inline-formula><mml:math><mml:msup><mml:mi>r</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn>1.2</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> density profile) using six primary CR species, protons, alpha particles, carbon, oxygen, silicon, and iron, drawn from GALPROP-HELMOD local interstellar spectra over <inline-formula><mml:math><mml:mrow><mml:mn>1</mml:mn><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>GeV</mml:mi><mml:mtext> </mml:mtext><mml:msup><mml:mrow><mml:mi>nuc</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> to <inline-formula><mml:math><mml:mrow><mml:mn>1</mml:mn><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>PeV</mml:mi><mml:mtext> </mml:mtext><mml:msup><mml:mrow><mml:mi>nuc</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, the hadronic interaction regime in which pion production and nuclear spallation dominate the secondary particle output. We compute, for the first time within this framework, radial energy deposition profiles, per-layer secondary production rates, and the differential physical flux of all emerging stable secondaries for the full nuclear CR composition. Energy deposition scales as <inline-formula><mml:math><mml:msup><mml:mi>Z</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:math></inline-formula>, with the secondary cascade dominating <inline-formula><mml:math><mml:mo>&gt;=</mml:mo><mml:mn>92</mml:mn><mml:mo>%</mml:mo></mml:math></inline-formula> of the deposited energy in the dense core for all species. Alpha particles produce the highest integrated gamma-ray flux (<inline-formula><mml:math><mml:mrow><mml:mn>0.125</mml:mn><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:msup><mml:mrow><mml:mi>cm</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mtext> </mml:mtext><mml:msup><mml:mrow><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mtext> </mml:mtext><mml:msup><mml:mrow><mml:mi>sr</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> above 100 MeV), followed by protons (<inline-formula><mml:math><mml:mrow><mml:mn>6.71</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:msup><mml:mrow><mml:mi>cm</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mtext> </mml:mtext><mml:msup><mml:mrow><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mtext> </mml:mtext><mml:msup><mml:mrow><mml:mi>sr</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>); together, they account for <inline-formula><mml:math><mml:mo>≈</mml:mo><mml:mn>96</mml:mn><mml:mo>%</mml:mo></mml:math></inline-formula> of the total emerging gamma-ray and neutrino fluxes. A robust composition-diagnostic signature emerges in the neutrino sector: the <inline-formula><mml:math><mml:msub><mml:mi>ν</mml:mi><mml:mi>e</mml:mi></mml:msub><mml:mo>/</mml:mo><mml:msub><mml:mover><mml:mi>ν</mml:mi><mml:mo>̄</mml:mo></mml:mover><mml:mi>e</mml:mi></mml:msub></mml:math></inline-formula> ratio above 100 MeV decreases from <inline-formula><mml:math><mml:mo>≈</mml:mo><mml:mn>6.4</mml:mn></mml:math></inline-formula> for protons, driven by the <inline-formula><mml:math><mml:mrow><mml:msup><mml:mrow><mml:mi>π</mml:mi></mml:mrow><mml:mrow><mml:mo>+</mml:mo></mml:mrow></mml:msup><mml:mo>/</mml:mo><mml:msup><mml:mrow><mml:mi>π</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo></mml:mrow></mml:msup><mml:mo>≈</mml:mo><mml:mn>8</mml:mn><mml:mi>--</mml:mi><mml:mn>11</mml:mn></mml:mrow></mml:math></inline-formula> charge asymmetry, to <inline-formula><mml:math><mml:mo>≈</mml:mo><mml:mn>2.8</mml:mn></mml:math></inline-formula> for <inline-formula><mml:math><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mi>Z</mml:mi></mml:math></inline-formula> nuclei (C, O, Si), directly encoding the isospin content of the primary CR beam in an observable flux ratio. Secondary nuclear clusters (<inline-formula><mml:math><mml:mi>d</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math><mml:mi>t</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math><mml:mrow><mml:mmultiscripts><mml:mrow><mml:mi>He</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:mmultiscripts></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math><mml:mrow><mml:mi>α</mml:mi></mml:mrow></mml:math></inline-formula>, and heavier ions up to <inline-formula><mml:math><mml:mi>Z</mml:mi><mml:mo>=</mml:mo><mml:mn>26</mml:mn></mml:math></inline-formula>) emerge exclusively for nuclei heavier than protons, contributing at the <inline-formula><mml:math><mml:mo>≈</mml:mo><mml:mn>32</mml:mn><mml:mo>%</mml:mo></mml:math></inline-formula> level of the secondary proton flux for light clusters, while heavy-ion fragments remain 1 to 2 orders of magnitude below the secondary proton channel.

Within the entropic gravity paradigm introduced by Verlinde, one assumes that the microscopic degrees of freedom residing on the holographic screen obey the equipartition law of energy. Nevertheless, implications of statistical mechanics suggest that this energy sharing can acquire corrections that depend on temperature. Taking such modifications into account leads to altered gravitational field equations when derived from thermodynamic considerations. We solve the resulting modified Einstein equations in the case of a static, spherically symmetric spacetime and determine the general structure of the metric components. Our findings indicate that if the temperature correction function behaves as $f(T)\propto T^{2}$, the corresponding spacetime geometry reduces to the flat Minkowski metric. Therefore, to obtain small deviations from flatness, it is necessary to consider slight departures from this purely quadratic temperature dependence. Such deviations can be interpreted as encoding additional gravitational effects, potentially associated with matter contributions that are not entirely described by the holographic screen. As a result of these effects, the effective gravitational potential acquires a logarithmic correction, which can give rise to deviations from Newtonian gravity. This aspect is especially significant in the regime of very weak gravitational fields, where these corrections may affect particle motion and could carry implications for both astrophysical and cosmological contexts.

First-order phase transitions in the early universe can generate a stochastic background of gravitational waves, offering a unique probe of high-energy physics. In this work, we investigate aspects of bubble nucleation and transition timescales, which play a central role in shaping the resulting gravitational wave spectrum. Many common approaches characterise the transition rate via a Taylor expansion of the false vacuum decay rate. We argue that a more fundamental description is instead given by the distribution of bubble lifetimes, and define a new timescale, $β_ν$, as the first moment of this distribution. We show that $β_ν$ reproduces the behaviour of previous timescale definitions in the appropriate limits, while avoiding their pathologies, and offers a more natural description of the ensemble of nucleated bubbles. We then quantify the impact of this improved timescale on the predicted gravitational wave spectrum, finding that it shifts the peak amplitude by up to an order of magnitude relative to previous definitions. As next-generation gravitational wave detectors come online, robust theoretical predictions will be essential; we hope this work represents a step in that direction.

Within the framework of quasi-topological gravity, we systematically investigate the correspondences among three classes of observables---the shadow, the quasinormal modes, and the graybody factors---of regular black holes in $D = 5$ dimensions. Using the WKB approximation and GrayHawk numerical integration, we test the applicability of the GBF--QNM correspondence and the shadow--GBF correspondence in the low-mode ($l=2$) and higher-dimensions ($D=5$): the difference between the GBF--QNM correspondence and the numerical results remains small in magnitude; we apply the known Langer correction (replacing $l$ by the effective angular momentum $κ= \sqrt{l(l+D-3)} \simeq l + (D-3)/2$), which reduces the error of the shadow--GBF correspondence at low modes to a level comparable with that of the GBF--QNM correspondence. Our work extends the shadow--GBF correspondence to higher-dimensional regular black holes, providing a viable route to predict the black hole spectrum and the Hawking radiation profile from a single observable, and thereby to realize multi-messenger tests of gravity.

A large sample of OB stars plays an important role in studying the stellar parameters of massive stars, as well as the formation and evolution of the Milky Way. With the help of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Data Release 10 (DR10), we are able to construct a large sample of OB stars with spectroscopic data. In this study, we identify 48,463 spectra of 34,550 OB stars from LAMOST DR10, based on the Hertzsprung-Russell (H-R) diagram constructed with Gaia DR3 data and spectral line indices measured from LAMOST DR10 low-resolution spectra. Among these, 6907 OB stars are newly identified. We use the MKCLASS tool to derive the spectral subtypes of the OB sample. The spatial distribution of 25,287 OB stars and the Toomre diagram of 20,397 OB stars indicate that the majority of these stars are located in the Galactic disk. Based on their peculiar velocities, we identify 1960 runaway star candidates.

The distribution of bio-essential elements is key to assessing the chemical fertility of galactic environments. We analyse phosphorus abundances for 233 nearby FGK stars from the Hypatia Catalog, separating thin- and thick-disk populations. We study [P/Fe] as a function of [Fe/H] and other elemental abundances using non-parametric tests and bootstrap resampling, and we examine residual correlations at fixed [Fe/H] together with molar ratios involving P. We find a [P/Fe]-[Fe/H] trend broadly similar to that of alpha elements, supporting a dominant origin in core-collapse supernovae, while deviations point to additional nucleosynthetic channels. Phosphorus shows strong positive correlations with several CHNOPS, alpha, and Fe-peak elements, especially in the thin disk. These co-enrichment patterns remain significant at fixed [Fe/H], indicating that they are not a trivial consequence of global metallicity. Thin-disk stars are also more enriched than thick-disk stars in many elements. The solar P/Mg and P/Si ratios are representative of the local stellar population, whereas solar P/O is relatively elevated. Overall, phosphorus emerges as a tracer of genuine multi-element co-enrichment and of the chemical fertility of the solar neighbourhood.

Context. The Pennsylvania-Toruń Planet Search, which operates on a sample of $\sim$1000 northern-hemisphere stars, has been underway since 2004. Stars with radial velocity (RV) amplitudes exceeding 2 km s$^{-1}$ represent a backup subprogram sample dedicated to monitoring binary stars with multiple instruments. Aims. We used almost 21 years of combined RV measurements to search for Doppler signals consistent with stellar or brown dwarf companions and to produce a catalog of both known and previously unpublished binary stars in our planet-search sample. Methods. We analyzed the combined RV measurements, searching for stellar companions and obtaining orbital solutions for both known and new binary systems. We also searched for periods in available long-term photometric monitoring data from All Sky Automated Survey (ASAS), and applied Gaia astrometric data whenever available. Results. We report the results of long-term RV monitoring of 16 systems: new detections of low-mass companions to 11 stars and updated orbital elements based on combined sets of our new and literature data for five systems. For two objects (TYC 3318-00789-1 and TYC 3318-01538-1), we find evidence of activity or the influence of another unseen companion, possibly due to line-profile variations or unresolved spectral contamination. We present true masses for three systems in our sample: TYC 1931-1040-1, TYC 3451- 1449-1, and TYC 3667-1636-1. We obtained these masses using Gaia DR3 non-single star data.

The estimation of stellar angular diameters can be performed from the surface brightness - color relation (SBCR) and photometry. The SBCR have been considered by the PLATO space mission as an independent empirical alternative for estimating stellar radii of FGK stars. In this context, we have implemented an homogeneous approach not only for calibrating the SBCR for FGK-IV/V stars but, also for determining their fundamental parameters in order to place the stars reliably on the HR diagram and to study their impact on the SBCR. We have performed interferometric observations of 18 quiescent FGK-IV/V stars in the Gaia color range of $3.088 \leq G \leq 5.498$. For the first time, we used 3 different interferometers operating in the $R$, $H$ , and $K$ bands to measure polychromatic limb-darkened angular diameters (LDAD). In parallel, by using public domain spectra we have estimated the stellar parameters ($T_{eff}$, $log g$ and Z) by using the open python tool iSpec. We achieved an average accuracy of 2.3% for the LDAD based on polychromatic observations. However, we have observed that our SBCR does not follow the calibration of the relation between surface brightness and color in Gaia found in the literature. Furthermore, we found that $log g$ and $Z$ have no impact on the SBC relation; given the characteristics of our sample of quiescent stars, it constitutes an ideal set of targets for conducting a new SBCR calibration within the framework of the PLATO space mission. In this context, we reported a SBCR calibration with $σ_{RMS} = $ 0.012, 0.009, 0.009 in the $G, G_{BP}$ and $G_{RP}$ Gaia bands, respectively. This article is part of a series of papers reporting the first results obtained using a polychromatic approach to measure LDAD, employing a fully homogeneous methodology for both determining the fundamental parameters of stars and measuring the LDAD, $θ_{LD}$.

Strong magnetic fields affect particle dynamics through two distinct channels: gravitational backreaction deforms the spacetime, while direct coupling to an intrinsic magnetic moment depends on the relative orientation of the field and the dipole. We disentangle these effects by studying equatorial orbits and near-horizon collisions of electrically neutral magnetized particles in the exact Kerr--Bertotti--Robinson spacetime. The field-induced geometric deformation shifts turning points and circular-orbit domains and can eliminate a finite effective-potential well together with its bound orbits. Even without direct dipole coupling, it can also offset the Kerr periapsis advance and produce a finite-radius zero-precession orbit. Direct dipole coupling breaks the symmetry under magnetic-field reversal and shifts the radius, energy, and angular momentum of the innermost stable circular orbit in an orientation-dependent manner. The formal ultrarelativistic endpoints of these orbit branches, however, remain fixed by the background geometry and approach circular null orbits. In the Bañados--Silk--West mechanism, both magnetic effects modify finite-radius potential barriers and hence the ability of a critical particle to reach the near-horizon collision region. In the representative nonzero-field cases examined here, an exactly critical particle released from infinity is blocked before reaching an extremal horizon, although a locally admissible collision between critical and usual particles can still produce unbounded center-of-mass energy. Finite dipole coupling shifts the barriers but does not change the leading near-horizon divergence. Near a nonextremal horizon, an exactly critical particle is excluded and the collision energy remains finite.

The Epoch of Reionisation imprints its history and morphology on the Cosmic Microwave Background temperature and polarisation anisotropies through two effects: The kinetic Sunyaev Zel'dovich (kSZ) effect and Thomson scattering. To perform joint analyses of these observables and put tight constraints on the reionisation history and morphology, a common reionisation model is necessary. We present a quick, physically-accurate analytical approach to derive consistently the angular power spectra of the three main resulting imprints that are the spatial fluctuations of the Thomson optical depth, the patchy kSZ effect, and the scattering and screening B-modes. The approach differs from existing analytical models, as it is calibrated on high-resolution hydrodynamical simulations and depends on physical parameters of reionisation. Our predicted power spectra are consistent with the ones derived from cosmological simulations, whilst including more physical information -- such as the typical size of ionised bubbles and the duration of reionisation. We provide an updated detectability assessment and show that both $C_\ell^{ττ}$ and $C_\ell^{BB}$ could be measured by experiments similar to CMB Stage 4 and CMB-HD, respectively. We show the complementarity of the three observables in constraining reionisation history and morphology, and illustrate the potential of their joint analysis to get a global picture of reionisation. This analytical model will be a powerful tool in the analysis of upcoming CMB data, either alone or in combination with independent datasets such as measurements of the high-redshift 21cm signal. The code to derive these spectra is publicly available online as the preion Python package.

We studied the evolution of a butterfly-shaped coronal hole (CH) observed from 2025 September 8 to 14 using SDO/AIA 193~Å images and SDO/HMI line-of-sight magnetograms. The CH boundary was identified using a fixed 100~DN intensity threshold, and its area, mean EUV intensity, and total unsigned magnetic flux were tracked across the solar disk. As the CH approached the central meridian, its area increased from approximately $7.70\times10^{16}$ to $1.63\times10^{17}$~m$^{2}$ and its unsigned magnetic flux from $8.68\times10^{21}$ to $1.91\times10^{22}$~Mx, while the mean 193~Å intensity decreased from about 61 to 45~DN. The area and unsigned flux showed similar temporal trends, whereas the intensity varied oppositely. Near central-meridian passage, these quantities remained relatively stable for several days. A high-speed solar-wind stream was detected at 1~AU about 2--4 days later, reaching 600--750~km~s$^{-1}$, consistent with the CH being its likely solar source. These results link the morphological, radiative, and magnetic evolution of the CH with the subsequent high-speed solar-wind stream.

We consider the difference between ensemble and volume average in cosmology. It is known that for sufficiently weak long-range correlations the root mean square of the difference, which we call ergodicity bias, decays like $R^{-3/2}$ in the limit of large volume $R^3$. We calculate the condition this imposes on the power spectrum of a Gaussian random field. We quantify the bias for finite $R$, and show that the $R\to\infty$ limit is of little relevance for cosmological observations when the measured scales and correlations extend to the size of the observable universe. We consider curvature, density, and velocity perturbations. On large scales the bias is important in all three cases. For the density perturbations, which are observationally the most relevant, the relative bias first exceeds 100% at the separation $r=177$ Mpc, and is larger than 100% for all $r>560$ Mpc. It should be taken into account when comparing ensemble and volume averages for large-scale structure. The bias is also large for the cosmic microwave background temperature perturbations on large angular scales, but this is not relevant for observations, as their analysis does not involve volume averaging.

We investigate the Saha ionization equilibrium in the exterior Schwarzschild spacetime. Starting from the mass-shell condition for a massive particle, we derive the conserved energy associated with the timelike Killing symmetry and obtain its nonrelativistic limit in terms of the Schwarzschild lapse function. This energy is then incorporated into the Maxwell--Boltzmann distribution to construct the ionization balance as seen by an asymptotic observer. We show that the gravitational field modifies the Saha relation through the redshift of both the particle energies and the thermodynamic variables. By imposing the Tolman--Ehrenfest and Klein equilibrium conditions, the formulation expressed entirely in terms of locally measured quantities reduces to the standard flat-spacetime Saha equation. We further analyze the gravitational redshift of the ionization energy, the relative departure from the flat-spacetime ionization ratio, and the corresponding shift of the ionization fraction as a function of the asymptotic temperature. We also present the relativistic classical counterpart based on the Maxwell--Jüttner distribution and show explicitly that its leading correction is negligible in the temperature range where neutral hydrogen is appreciably abundant. These results provide a consistent description of ionization equilibrium in a static gravitational field and clarify the distinction between gravitational redshift effects and intrinsic modifications of local atomic physics.

The Auger Radio Infill SKALA Extension (ARISE) at the Pierre Auger Observatory in Argentina was deployed in 2025 and measures cosmic-ray air showers in the energy region of the Galactic-to-extragalactic transition. ARISE is comprised of 18 SKALA-2 antennas featuring two polarization channels each, deployed within $100\,$m of a surface detector station in the enhancement area of the Pierre Auger Observatory. This area of the surface array features a denser spacing of $433\,$m between surface stations, each equipped with underground muon detectors. One of these surface detector stations provides a trigger for simultaneous readout of all ARISE antenna channels. The wide frequency range of ARISE, from $50$ to $350\,$MHz, includes the sub-band of optimum signal-to-noise ratio for air-shower radio emission against the Galactic radio background. In combination with the dense antenna spacing, this enables a relatively low detection threshold, and ARISE aims at demonstrating full detection efficiency for near-vertical air showers above $100\,$ PeV. As an advantage over the current radio detectors at Auger, which are more efficient for inclined air showers, this would enable low systematic uncertainties for physics analysis combining ARISE radio measurements with coincident measurements of the underground muon detectors in the same area. In this presentation, we will provide an overview over the ARISE instrumentation operating at the Pierre Auger Observatory and will outline the science goals.

Polluted white dwarfs provide a unique window into the composition of exoplanetary material, but interpreting their surface abundances requires a quantitative understanding of the mixing processes in their envelopes. Accretion of metal-rich debris onto hydrogen-rich white dwarfs drives thermohaline (fingering) convection in the underlying radiative zone. Previous studies have modeled this as enhanced diffusive transport, often inferring mass accretion rates, $\dot{M}_{\rm acc}$, that exceed constraints from X-ray observations. Here we develop a time-dependent model that treats thermohaline mixing as a propagating front whose evolution is coupled self-consistently to the surface abundance of heavy elements. We perform hydrodynamical simulations which demonstrate that thermohaline mixing drives the composition gradient towards marginal stability. Using this condition, we derive analytic solutions in shallow and deep stratification regimes, finding that the mixed layer depth, $h$, and surface abundance, $X_{\rm surf}$, both scale as $t^{1/2}$ at early times, transitioning to $h \propto t^{4/27}$ and $X_{\rm surf} \propto t^{10/27}$ at later times. The mixed layer reaches depths of only $\sim 100$ km over $10^5$ yr and, without gravitational settling, never reaches a steady state. Including settling leads to an equilibrium surface abundance scaling as $X_{\rm surf, eq} \propto \dot{M}_{\rm acc}^{30/43}$. Applied to G 29--38, the model matches the observed heavy-element abundance for $\dot{M}_{\rm acc} \sim 2\times 10^9$--$4\times 10^{9}~\mathrm{g~s^{-1}}$, in better agreement with X-ray constraints than prior estimates incorporating thermohaline mixing. We also present a general scaling relation for inferring accretion rates directly from observed surface abundances and stellar properties, giving a practical prescription for other polluted white dwarfs.

We present a map of the cosmic microwave background (CMB) lensing potential reconstructed from observations taken during the 2019 and 2020 seasons with the third-generation camera on the South Pole Telescope (SPT), covering the $1500\,{\rm deg}^{2}$ SPT-3G Main field, referred to as the SPT-3G D1 dataset. From the multi-frequency temperature and polarization data, we reconstruct the CMB lensing field using a quadratic estimator that jointly accounts for the $T$, $E$, and $B$ fields and their covariance. The resulting lensing map is dominated by polarization information for $L \lesssim 600$ and provides the highest signal-to-noise measurement per mode reported to date. With nuisance parameters fixed to their best-fit values, we measure a lensing amplitude consistent with unity at $2\%$ precision relative to the $Λ$CDM model that best fits the combined Planck, ACT DR6, and SPT-3G D1 $TT/TE/EE$ likelihoods (${\rm CMB}_{\rm SPA}$). We further measure the structure-growth parameter $σ_{8}Ω_{\rm m}^{0.25}$ to be $0.6046\pm0.0096$ from the SPT-3G D1 lensing spectrum alone and $0.6020\pm0.0084$ when combined with ACT DR6 and Planck PR4 CMB lensing. By further combining this with ${\rm CMB}_{\rm SPA}$ and the latest DESI DR2 BAO data, we obtain $\sum m_ν < 0.072\,\mathrm{eV}$ (95% C.L.) when allowing the neutrino mass to vary within $Λ$CDM. Compared with previous work, the better agreement of our measurement with DESI DR2 BAO yields both this relaxed upper bound and reduced ($\mathord{\sim}2σ$) preferences for nonzero spatial curvature and for deviations of $(w_0,w_a)$ from the $Λ$CDM expectation. When we combine CMB lensing with the DES Y3 3$\times$2pt analysis, we obtain $S_{8}=0.811\pm0.011$, corresponding to a $1.4\%$ constraint on the late-time clustering amplitude. This precision is competitive with that obtained from the primary CMB within $Λ$CDM.

It was suggested by 't Hooft that to avoid the divergence of the free energy or entropy of the probe fields, one has to introduce an infrared (IR) cutoff. This introduction of the IR cutoff, namely, the brick wall model, has been studied for the Schwarzschild black hole and has been shown by 't Hooft that to have the entropy of the order of Bekenstein-Hawking entropy, the cutoff must be of the order of Planck length. In this study, we investigated the brick wall model in the context of a Reissner-Nordstróm, quantum corrected Schwarzschild and quantum corrected Reissner-Nordström black holes. In particular, we have shown that the brick wall model can give logarithmic correction to the black hole entropy in the subleading order. The important aspect of our analysis is that we have not considered the near horizon approximation of the lapse function and get more accurate correction to the entropy. This corrected entropy has an important feature that it is larger than the usual Bekenstein-Hawking entropy which aligns with the result of Barrow entropy.

Gyrochronology ties stellar rotation periods to ages. It is well studied in the young, open cluster age domain, but there are few constraints on gyrochronology in the old field star regime. In this work we use gyro-kinematic ages to explore the spin down of stars in this formerly inaccessible domain. Using forward modeling techniques which relax strict Rossby scaled assumptions, we find evidence for a departure from a standard spin down models. This departure can be explained with a mass-dependent term either in the global strength of the stellar wind or in the relationship between angular velocity and wind strength. Models with this additional mass- dependence help explain prior difficulty in fitting open cluster rotation distributions across the full mass range. Additionally, we use rotation models to identify an mass-dependent age domain over which wind driven stellar spin down can be isolated from other physical effects. For lower mass stars, this age domain is most affected by the core-envelope coupling timescale, where as higher mass stars are more subject to inertial effects late in their main sequence lifetimes. Using simple, analytic models of stellar spin down in this domain is unable to determine whether the stellar wind has a mass-dependence or if the stellar wind is strictly Skumanich-like in nature. In future studies of rotational and dynamo evolution we advocate for a forward modeling approach to gyrochronology over purely empirical approaches due to its superior ability to trace physical effects governing stellar spin down.

We show that the subleading collinear sector of Yang-Mills (YM) amplitudes is controlled by ordinary Einstein-Yang-Mills (EYM) amplitudes and their higher-derivative corrections. The complete set of strict subleading collinear limits associated with an equal-helicity collinear pair can be extracted from open-closed string disk amplitudes. The latter generate $\frac12 (n-3)!$ BCJ-like relations with non-linear kinematic coefficients, reducing the $(n-3)!$ collinear data to a basis of dimension $\frac12 (n-3)!$. This reduced basis can be represented by gravitational amplitudes in EYM theories and by their higher-derivative open-string corrections. At multiplicity $n$ this gravitational basis consists of the $(n-4)!$ independent EYM subamplitudes together with $\frac12 (n-5)(n-4)!$ higher-order corrections involving one graviton and $n-2$ gluons. The resulting decomposition is governed by unsigned Stirling numbers of the first kind: the ordinary EYM amplitudes correspond to the sector $\left[{n-3\atop 1}\right]$, the higher-order BCJ-like relations to even sectors $\left[{n-3\atop 2j}\right]$, and the higher-derivative EYM corrections to odd sectors $\left[{n-3\atop 2j+1}\right]$ with $j\geq 1$.

Surface Flux Transport (SFT) models are routinely used to model the Sun's photospheric magnetic field and provide inner boundary conditions for coronal and heliospheric models, yet they remain fundamentally limited by the lack of direct information about flux emergence on the far side of the Sun. To address this, we develop a framework for incorporating helioseismic images of far-side active regions (HIFAR), into the Advective Flux Transport (AFT) model. Using near-simultaneous magnetic flux maps (HIFARM) inferred from HIFAR, and STEREO/EUVI 304 A observations from 2010 May 13 to 2014 August 18, we develop a logistic regression model to estimate the probability that a helioseismic detection corresponds to an AR based on descriptors of magnetic flux, field strength, and location. The model achieves a precision of 0.93 for active region (AR) detections at a probability threshold of 0.73, selected to reject 90% of "ghost" ARs. We also derive an empirical scaling relationship between HIFARM and AFT flux to place HIFARM on the AFT flux scale and demonstrate that the framework reproduces the flux evolution of ARs consistent with near-side observations over multiple solar rotations. At the global scale, we find that conventional near-side-driven AFT simulations underestimate the total unsigned magnetic flux by ~10-20% of the Sun's total magnetic flux budget. The probabilistic AR detection model and the complete framework developed in this work provide a practical pathway for incorporating helioseismic far-side ARs into SFT models, a step toward realistic modeling of full-Sun magnetic field.

We present a study of stellar limb-darkening using JWST transit observations of seven exoplanets orbiting FGK host stars, spanning 4200-6800 K. The wide wavelength coverage and high S/N of NIRISS/SOSS and NIRSpec/PRISM enable precise constraints on the wavelength-dependent limb-darkening. Using Bayesian model selection, we find that the quadratic limb-darkening law is statistically preferred over higher-order laws for most FGK stars, recovers consistent intensity profiles for $>$95% of the stellar disk area, and introduces a minimal bias in the derived transit depths of only $\sim$14 ppm (1$σ$) for 6/7 targets. This contrasts with previous studies relying on stellar models. We compare the empirically derived quadratic coefficients to predictions from the PHOENIX, MPS-ATLAS, MURaM, and Stagger stellar atmosphere grids. We introduce a quadratic limb-darkening parameterization in terms of limb intensity ($\ell$) and curvature at mid-$μ$ ($δ$), finding that empirical FGK limb darkening is generally more linear than models predict ($δ\lesssim 0.1$). We identify wavelength-independent offsets between data and model quadratic coefficients, minimized by adopting $μ_{\rm min} = 0.2$ in intensity calculations; we attribute this in part to models overpredicting limb-darkening near the limb where the plane-parallel approximation breaks down. For spherical PHOENIX models, we derive a $μ$ rescaling method using a $τ= 1$ photospheric radius. With these corrections, residual offsets are minimized and all stellar models achieve statistically acceptable fits. We provide recommended limb-darkening offset priors for use in JWST transit analyses, enabling more accurate constraints on exoplanet transmission spectra while accounting for residual stellar model uncertainties.

All observed planetary systems orbiting single white dwarfs have lived through the hot stellar transition from an asymptotic giant branch star. In this post-nebular transition period, the initial conditions for planetary system evolution throughout white dwarf cooling are established. The hottest ($\gtrsim$ 20,000 K) and youngest ($\lesssim$ 20 Myr-old) white dwarf planetary system host stars differ significantly from their canonical older and colder counterparts by failing to support solid body accumulation in the immediate vicinity (a few $R_{\odot}$) of the white dwarf. Here, we analyse the likely locations of both solid-body survival and the sublimated gaseous content during this pivotal epoch, and the consequences. We find that (i) reservoirs of iron-rich, rocky and water-rich asteroids of radius $R$ that later observably enrich, or pollute, the white dwarf need to remain parked for the first tens of Myr of white dwarf cooling beyond critical distances of (16 au)$\times(1 {\rm km}/R)^{1/2}$ (for iron), (30 au)$\times(1 {\rm km}/R)^{1/2}$ (for rock) and (130 au)$\times(1 {\rm km}/R)^{1/2}$ (for snow), (ii) sublimation acts much more quickly than radiatively-driven orbital drifts from Poynting-Robertson drag or the Yarkovsky effect, and (iii) although large asteroids ($R \approx$ 10-1000 km) that are kicked on highly eccentric orbits around newly born white dwarfs could survive sublimation, they may fragment into debris which will sublime before the white dwarf cools. These results support, but do not necessitate, dynamical origin scenarios of polluted white dwarfs that feature delayed gravitational instability subsequent to the host star's asymptotic giant branch phase at Kuiper Belt-like distances, and beyond.

The majority of very close binaries (orbital periods P < 5 days) are in triples. Previous studies surmised that the tertiary hardened the inner binary via Kozai-Lidov cycles coupled with tidal friction. Alternatively, the trend could be due to a common origin, e.g., more massive molecular cores and protostellar disks are more likely to fragment twice into a triple and harden the inner binary to very short periods via disk migration. Disk migration occurs during the embedded Class I protostellar phase while secular evolution with a tertiary occurs on much longer timescales. To determine the dominant migration mechanism, I accurately measure the very close binary fractions of both late-type and early-type primaries in star-forming environments. After combining five surveys, the late-type Class II/III T Tauri very close binary fraction is F_VCB = 2.1_{-0.5}^{+1.0}%, which is consistent with the measurement for late-type field main-sequence stars (1.7% +/- 0.1%). More than 85% of very close late-type binaries hardened via disk migration during the Class I protostellar phase. I next measure the OGLE-III eclipsing binary (EB) fractions of young massive stars in the LMC. The very close EB fractions of massive young stellar objects (YSOs; F_VCEB = 8.3_{-3.2}^{+4.9}%), OB stars near massive YSOs (4.7% +/- 0.8%), and OB stars in small, young H II regions (6.4% +/- 1.5%) are consistent with each other and the measurement for OB stars in slightly older clusters (4.7% +/- 0.3%). More than 90% of very close massive binaries hardened via disk migration during the embedded protostellar phase.

The origin of supermassive black holes (SMBHs) in the early Universe remains uncertain. While the classical direct-collapse scenario requires pristine gas exposed to intense Lyman-Werner (LW) radiation, recent JWST observations suggest that early growing BHs may reside in compact, metal-enriched environments. We investigate early BH formation using a cosmological radiation-hydrodynamic simulation that self-consistently follows the formation and evolution of rapidly accreting supermassive stars (SMSs). We find that intense LW radiation develops in clustered star-forming regions; by the time the first heavy seed forms, the mass-weighted internal intensity has reached $J_{21}\sim1000$. Most seeds form in weakly metal-enriched gas with [Z/H]$\sim-3$ to $-2$, because the strongest LW fields occur in regions already enriched by previous stellar feedback. All halos above $10^9M_\odot$ in our zoom-in region host a BH more massive than $10^5M_\odot$. We construct stellar spectra using MIST/PHOENIX libraries for normal stars and $6000~\mathrm{K}$ blackbody spectra for bloated SMSs. Five of eight massive clusters with stellar masses larger than $10^6~M_\odot$ exhibit intrinsically V-shaped spectra similar to those of Little Red Dots (LRDs). In the most prominent cases, luminous SMSs dominate the red optical emission, although Balmer breaks in normal stellar populations also contribute. We also identify off-center LRD-like systems in our simulation resembling those recently reported in strongly lensed observations. Our results suggest that heavy-seed formation, compact-cluster formation, and LRD-like continuum emission are linked outcomes of clustered, weakly metal-enriched, strongly irradiated environments.

Quantum mutual information is a fundamental quantity for characterizing correlations among quantum systems. In this work, we investigate the dynamic evolution of the quantum mutual information based on Rényi-2 entropy for three-mode Gaussian states in the background of a Schwarzschild black hole. We find that the physically inaccessible mutual information increases with Hawking temperature and eventually saturates, while the physically accessible mutual information exhibits nonmonotonic temperature dependence: it first rises, reaches a peak, then declines toward a finite asymptotic value. This nonmonotonic behavior differs from the monotonic degradation typically observed for Gaussian entanglement and steering in curved spacetime. Furthermore, we establish several constraint relations governing the distribution of mutual information among subsystems. These findings contribute to understanding the quantumness of quantum correlations for continuous variable in Schwarzschild spacetime.

CO isotopologue emission is widely used to infer gas masses and volatile carbon abundances in protoplanetary disks, but converting line emission into a CO column depends on optical depth, temperature structure, linewidth, and isotope ratios. Inferring CO/H$_2$ additionally requires an independent constraint on the local hydrogen column. We use high-resolution $^{13}$C$^{18}$O $2$-$1$ observations of HD 163296 to derive a spatially localized empirical constraint on the CO column at the resolved CO snowline edge. We focus on the 70-75 au annulus, on the inner, high-column side of the observed profile steepening near 75 au. Using RADEX slab calculations conditioned on a two-dimensional temperature structure, we infer an effective beam-averaged CO column from the absolute integrated intensity. Across representative temperatures, isotope-ratio pairs, and effective local linewidths of 0.30 and 0.50 km s$^{-1}$, we find $N_{\rm CO}^{\rm beam}=(1.6-2.4)\times10^{20}$ cm$^{-2}$ and $^{13}$C$^{18}$O line-center optical depths $τ=0.39-0.97$. Thus, even this rare isotopologue is not safely optically thin at the snowline edge. Adopting $N_{\rm H2}=2.5\times10^{24}$ cm$^{-2}$ from a published parametric gas surface-density profile gives the conditional abundance ${\rm CO/H_2}=(6.4-9.5)\times10^{-5}$. A beam-forward radial-profile analysis gives consistent columns, and a physical disk model with a near-canonical warm-layer CO abundance supplies a comparable CO column. The measurement is consistent with the higher C$^{17}$O-based MAPS estimate. For the adopted hydrogen column, the inferred CO/H$_2$ ratio is consistent with a near-canonical abundance on the warm side of the snowline.

Helicity in gravity is a multifaceted concept, and several inequivalent definitions exist in the literature. In this work, we reduce the topological Pontryagin and Nieh-Yan terms to a constant time spacelike hypersurface. Remarkably, upon applying the SVT decomposition at the linearized level, we obtain four distinct helicity densities. The helicity density arising from the Nieh-Yan term is termed the spin-1 and spin-2 gravitomagnetic helicity density, depending on whether it originates from vector or tensor modes. Meanwhile, the helicity density arising from the Pontryagin term is termed the spin-1 and spin-2 gravito-current helicity density, according to the corresponding mode contributions. We study the mode and multipole expansions of these helicity functionals and apply them to leading order Newtonian two-body systems, weak field boosted Kerr, a locally defined slowly varying gyratonic pp-wave model, and a linearized gravitational Hopfion.

We study the algebraic and differential structure of the Riemann curvature tensor of the higher-dimensional Kerr-NUT-(A)dS spaces, motivated by the eventual goal of finding an IDEAL (Intrinsic, Deductive, Explicit and ALgorithmic) characterization of this family of metrics. The special geometry of these spaces is governed by the principal tensor $h_{ab}$, a non-degenerate closed conformal Killing-Yano $2$-form, whose existence singles out the Kerr-NUT-(A)dS family; we therefore focus on the algebraic relationship between the Riemann curvature $R_{abcd}$ and $h_{ab}$. In our investigations we encountered an obstruction, which results in a no-go theorem: no non-trivial $2$-form, including $h_{ab}$ itself, can be constructed covariantly from the undifferentiated Riemann tensor alone. Despite that, we do characterize the family of Riemann-symmetric tensors that could be the curvature of a Kerr-NUT-(A)dS metric by an algebraic and a differential condition, stemming from the integrability of the defining equation of $h_{ab}$ and the second Bianchi identity. Our calculations apply in all higher dimensions, which becomes feasible by a judicious application of representation theoretic techniques related to a semi-direct product group $U(1)^n\rtimes S_n$ that stabilizes $h_{ab}$.

We investigate the equilibrium abundance ratio of the molecular hydrogen ion \(\mathrm{H}_2^+\) and the negative hydrogen ion \(\mathrm{H}^-\) in stellar atmospheres using compact Saha-type equilibrium expressions and layer-by-layer partial pressures from MARCS model atmospheres. The analysis distinguishes the temperature-dependent coefficient ratio \(K_{\mathrm{H}_2^+}(T)/K_{\mathrm{H}^-}(T)\) from the full abundance ratio, which also depends on the local ionisation factor \(N_{\mathrm{H}^+}/N_e\). The internal partition function \(Z_{\mathrm{H}_2^+}(T)\) is included in the \(\mathrm{H}_2^+\) equilibrium coefficient. We processed \(51{,}994\) MARCS models, corresponding to \(2{,}911{,}664\) atmospheric layers, and computed \(N(\mathrm{H}_2^+)/N(\mathrm{H}^-)\) directly from the tabulated equilibrium partial pressures, constructing the first systematic map of this ratio across the MARCS grid. \(\mathrm{H}_2^+\) exceeds \(\mathrm{H}^-\) in \(827{,}350\) layers, corresponding to \(F=0.284\), or \(28.4\%\) of all analysed layers. Although \(\mathrm{H}^-\) remains dominant in the global logarithmic mean, with \(\langle \log_{10}[N(\mathrm{H}_2^+)/N(\mathrm{H}^-)] \rangle=-4.229\), \(\mathrm{H}_2^+\) becomes comparable to or more abundant than \(\mathrm{H}^-\) in a substantial subset of layers. Its relative significance increases toward higher \(T_{\mathrm{eff}}\), lower \([\mathrm{Fe}/\mathrm{H}]\), and deeper layers. The balance is therefore controlled by the combined local thermodynamic and ionisation structure, not by temperature alone. These ratios are equilibrium abundance diagnostics and do not imply opacity dominance; assessing opacity requires wavelength-dependent cross-sections and radiative-transfer calculations.

We derive a $2d$ effective action for semi-classical AdS$_3$ gravity with an EOW brane of arbitrary tension by solving the radial Hamiltonian flow in a derivative expansion. Interpreting the radial flow as an RG flow, the effective theory is identified with a $T\bar T$-like deformed CFT. This leads us to propose that the holographic dual for semi-classical AdS$_3$ gravity with Neumann or conformal boundary conditions is obtained by setting the $T\bar T$-like deformed CFT free. We verify our proposal explicitly for saddles corresponding to $(A)dS_2$ branes in empty AdS$_3$. We discuss three useful ways to decompose the effective brane theory into a matter sector ($T\bar T$-deformed, $T\bar T$-like deformed, and undeformed CFT, respectively) and a corresponding braneworld gravity sector. In particular, the split into $T\bar T$-deformed CFT coupled to $T\bar T$-like deformed Liouville gravity directly extends the Liouville description of the integrated Weyl anomaly into the bulk. The effective action is then used to perform an island-rule calculation for a non-asymptotic brane in an AdS/bCFT set-up, finding exact agreement with the Ryu-Takayanagi result. This extends previous holographic checks on the island rule to non-asymptotic regimes of application, as well as provides a Wald entropy interpretation for $T\bar T$ entanglement. Finally, we relate our effective theory and its higher derivative corrections to AdS$_3$ gravity with a fluctuating radial cutoff by explicit calculation of the on-shell action.

(Abbreviated) The interpretation of high-precision astrometry for intrinsically variable stars remains challenging, particularly for unresolved binary systems containing asymptotic giant branch (AGB) stars. In such systems, large-amplitude pulsations and evolving convective surface structures induce time-dependent photocentre displacements that perturb the observed orbital motion. At the same time, Gaia DR4 epoch astrometry offers the prospect of deriving accurate orbital solutions and parallaxes for nearby unresolved AGB binaries. Building upon the variability-induced mover (VIM) framework, we extend the astrometric model for unresolved evolved binaries by combining Keplerian photocentre motion, pulsation-induced flux-dependent photocentre shifts, and convection-induced photocentre motion. In the forward simulations, convection is represented either by stochastic photocentre displacements drawn from an exponential correlation function or by photocentre time series extracted from 3D radiation-hydrodynamic simulations. We then test a retrieval framework that fits the orbital and VIM signal while incorporating a red-noise covariance matrix to account for correlated astrometric residuals produced by convection. We show that Gaia epoch astrometry can recover reliable orbital and astrometric parameters for unresolved AGB binaries despite strong pulsation- and convection-induced photocentre variability. Pulsation-induced variability produces a coherent VIM signal that can be modeled jointly with the Keplerian photocentre orbit, while convection-induced photocentre motion behaves primarily as temporally correlated astrometric noise. Retrievals that ignore this correlated component lead to biased proper motions and parallaxes, whereas the inclusion of a physically motivated red-noise covariance model enables accurate recovery of the underlying orbital solution and astrometric parameters.

Interstellar grains play key roles in star and planet formation, including the coupling of the gas to the magnetic field during the protostellar collapse. These roles depend on the local grain size distribution, which requires a multifluid treatment of charged dust. We aim to understand the fundamental physics of the dynamics of a dust distribution in interaction with the gas and the magnetic field. In particular, the purpose is to characterize the (de)coupling conditions of these different components. We provide a multifluid model of charged dust which accounts for the inertia of the grains. A chemical network is used to simulate the charge equilibrium in collapsing protostellar cores. We compute the Alfven modes and the magnetosonic modes to understand the coupling regimes between the gas, the dust fluids, and the magnetic field. We also analyze and compare to the predictions of existing models, that are the neutral dust multifluid and the standard non-ideal magnetohydrodynamics. The charged multifluid model agrees with non-ideal magnetohydrodynamics on the larger scales of a collapsing dense core and the forming disk, while we successfully extend to new regimes where the inertia of dust grains matters. We found that high charge-to-mass dust grains carry the propagation of magnetohydrodynamical waves in protostellar envelopes. We provide analytical expressions of the speed of these waves depending on the dust distribution. The magnetocompressive perturbations lead to local dust-to-gas ratio variations at au scales. A theoretical understanding of the dynamics of a charged dust distribution is provided in the linear regime. The closed set of magnetohydrodynamics equations can be implemented in numerical codes to explore nonlinear effects during the protostellar collapse such as turbulence, angular momentum transport and magnetic dust clumping.

Our work clarifies the origin of enhanced symmetry of near horizon geometries. We relate this property to an existence of specific foliation by non-expanding horizons. The existence of this foliation can be shown in any spacetime satisfying null convergence condition. Our proof is geometric and based on analyzing null hypersurfaces in near horizon geometry spacetimes.

Large-scale, high-resolution spectroscopic surveys map detailed chemical abundance patterns for millions of stars. However, libraries of stellar evolution models and isochrones still rely on a simplified description of chemical abundances, generally assuming either solar-scaled compositions or a constant enrichment for all $α$-elements. Moreover, the definition of $α$-enrichment is survey-dependent and is not necessarily the most relevant quantity for stellar structure and evolution. In the era of precision stellar astrophysics, stellar models need to better account for the wealth of chemical information provided by spectroscopic surveys. We develop a perturbative method to compute stellar evolution tracks for stars with arbitrary patterns of $α$-elements. We quantify the response of stellar evolutionary tracks to variations of individual $α$-elements through polynomial expansions, which can be linearly combined to produce synthetic tracks for a given abundance pattern. We also introduce the evolutionary $α$-enhancement, [$α$/Fe]$_{\rm ev}$, a constant enhancement for all $α$-elements that reproduces the impact of an observed abundance pattern on stellar models. We validate both methods using stars with detailed APOGEE abundances and quantify the differences between spectroscopic [$α$/Fe] and [$α$/Fe]$_{\rm ev}$ for different surveys. We identify a minimal set of chemical elements required for accurate stellar models and present ALCHEMY, a simple and fast tool to compute [$α$/Fe]$_{\rm ev}$ for general $α$-element abundance patterns.

We study the existence of gravitational radiation at infinity for pure-radiation Robinson-Trautman metrics with one spacelike Killing vector and a cosmological constant $Λ$ of arbitrary sign. These metrics describe "photon rockets" with point, string, and sheet sources moving through spacetime while emitting null radiation. For $Λ\neq 0$, we use criteria for the existence of gravitational radiation based on the asymptotic super-Poynting vector. For point sources, we extend to $Λ\neq 0$ the known $Λ=0$ result that the Kinnersley rocket is the only rocket without gravitational radiation. However, for string and sheet sources, new possibilities without gravitational radiation arise. For $Λ>0$, we show that the vanishing of the canonical asymptotic super-Poynting vector---computed with respect to the unit normal to scri---exactly characterizes the absence of radiation. On the other hand, for $Λ<0$, we show that the appropriate criterion is the vanishing of the components normal to scri of the asymptotic super-Poynting vectors associated with any unit timelike vector tangent to scri or, equivalently, the proportionality between the Cotton-York tensor and the holographic stress tensor at scri. In particular, we provide explicit examples of metrics for which the Cotton-York and holographic stress tensors commute, yet gravitational radiation is present because the two tensors are not proportional. We also determine the principal null directions (PNDs) of the rescaled Weyl tensor at scri for both signs of $Λ$ and relate their geometry relative to the normal to scri to the tensorial criteria for gravitational radiation.

We introduce a family of effective theories and use them to define a specific type of coarse-grained models for canonical loop quantum gravity. Each effective theory is characterized by two parameters and it consists mainly of a Hilbert space of coarse states and an effective Hamiltonian operator. The construction of the coarse Hilbert spaces relies on a systematic coarse-graining procedure for spin network states. The effective Hamiltonians are then obtained from an analysis of the interplay between this coarse-graining procedure and the action of various Hamiltonian operators defined in loop quantum gravity. Finally, we provide a prescription for implementing the Hamiltonian non-perturbative renormalization framework for the coarse-grained models, and we explicitly write down the renormalization flow equations. The aim of this work is to build an arena for studying the emergence of the continuum limit and the derivation of phenomenological models in the context of canonical loop quantum gravity.

Complex periodic variables (CPVs) are young low-mass stars whose light curves show periodic dips indicative of transiting corotating material. The origin and composition of this material are unclear. Here we present new optical and near-infrared spectroscopy and photometry of four CPVs from Magellan, Keck, Hale, MuSCAT1, MuSCAT2, Tierras, KeplerCam, and TESS. The spectra imply that CPVs host magnetically bound circumstellar plasma clumps, on the basis of sinusoidal-in-time Balmer emission out of transit, and Balmer dimming during transit. Yet large night-to-night changes in circumstellar hydrogen emissivity occur without clear changes in light curve morphology, suggesting that the sharp flux dips are caused not by circumstellar plasma but by dust. Optical chromaticities (depth proportional to $λ^{-β}$, with $β$ = 0.79 $\pm$ 0.16) support this, but the power law breaks in the near-infrared, where a single power law under-predicts the depths observed at 2.1 microns. We therefore favor dips caused by dusty plasma clumps with opaque cores and optically thin halos, though we cannot rule out models in which the dust properties vary per-star or per-epoch. Observations at wavelengths greater than 2 microns and less than 0.4 microns would test this interpretation and clarify the dust's origin.

Observations with a high angular resolution by ALMA have revealed that substructures are ubiquitous in protoplanetary disks. Axisymmetric dust rings are the most common morphology. The profiles of some observed disks, including young disks, are smooth overall, with a localized dip or bump near the outer edge of the continuum disk that manifests as an edge ring. While embedded planets might contribute to these structures, their physical origin remains unclear. We investigated the possibility that these edge rings arise purely from radiative transfer effects at the outer edge of a protoplanetary disk. A steep surface density dust gradient at the outer edge of the disk allows stellar irradiation to penetrate more efficiently beyond the disk edge, producing a non-monotonic temperature profile characterized by a dip that is followed by a bump. We tested whether this non-monotonic temperature structure can generate and maintain a localized continuum enhancement. We coupled radiative transfer and dust evolution by iterating between the Monte Carlo radiative transfer code RADMC-3D and the dust evolution code DustPy. This framework self-consistently follows the coupled evolution of temperature, grain growth, and dust dynamics. Thermodynamic feedback at the disk edge can naturally generate and maintain localized dust enhancements resembling the edge rings that are observed in some extremely young disks. Without invoking planets or additional dynamical perturbations, this mechanism offers a plausible explanation for the first-generation ring formation. Our results highlight the importance of coupling thermodynamics and dust evolution when modeling protoplanetary disks, suggesting that thermodynamic feedback probably plays a role in shaping disk substructures.

We revisit the quantum Hall analog of Hawking radiation in which a Fermi sea of electrons occupying half of a plane, subjected to a quadrupolar electric potential, gives rise to analog Hawking radiation of chiral edge modes. We show that the phenomenon is accurately captured by quantization of the classical guiding center theory, where the gyroscopic motion of the electrons is coarse grained and only the drift motion is resolved. The quantum dynamics takes place on a non-commutative plane, which requires an electron localized in the half-plane to be supported everywhere along the edge direction. This kinematical constraint on the quantum state leads directly to the radiation, which propagates in both directions away from the origin along the edge. The radiation is thermal with respect to laboratory time, which is equal (up to a constant factor) to the boost angle in the analog Minkowski spacetime in which the chiral edge modes propagate, so in fact it corresponds to analog Unruh radiation.

We investigate cosmological diffusion models in unimodular gravity within a dynamical systems reconstruction framework. By treating the logarithmic slope of the diffusion sector as an invertible dynamical variable, the diffusion function can be systematically reconstructed from the phase-space structure of the cosmological evolution. Under these conditions, we determine the physically admissible fixed points of the system, identifying novel matter--diffusion scaling solutions associated with power-law diffusion sectors, as well as purely diffusion-dominated configurations capable of driving late-time accelerated expansion without requiring an explicit cosmological constant term. The local behavior around the fixed points is then extended to the full cosmological evolution, providing a framework to explore the global implications of diffusion cosmologies. Beyond the asymptotic fixed-point structure, we further develop a reconstruction formalism based on the dynamical evolution of the diffusion slope, allowing for trajectories interpolating between different diffusion regimes during the cosmic history. Our results establish a systematic framework for constructing and classifying viable diffusion cosmologies in unimodular gravity directly from the phase-space dynamics.

We consider the effects of radiative corrections in Starobinsky-like models of inflation, concentrating on models of inflation formulated in N=1 no-scale supergravity. Inflaton couplings to matter fields are necessary for reheating and these have an impact on loop corrections to the inflaton potential. Whilst corrections due to the supergravity couplings of the inflaton to Standard Model (MSSM) fields are negligible, we use {\it Planck} data to obtain interesting constraints on GUT bilinear couplings, vevs, and gauge boson masses that could be sharpened by future CMB measurements.

In this article, we establish a new dimensional reduction result in (Galilean or Lorentzian) Robertson Walker spacetimes with Euclidean fiber. In this context, we prove that every unchanged direction observer can be embedded into a three dimensional totally geodesic timelike submanifold.

Higher-gradient modifications of Newtonian gravity remove the point-mass singularity, but their Lagrangian dynamics is haunted by the Ostrogradsky instability. We show that the same field equations follow from the Second Law of thermodynamics applied to a self-gravitating fluid with a weakly nonlocal state space, and that on this route the instability is absent by construction: the potential obeys a first-order relaxation equation that is a gradient flow, whose Lyapunov function is the entropy density, and whose relaxation spectrum is negative definite at every wavenumber precisely when the entropy is concave. The derivation also yields the pressure tensor and the energy current, an exact identity that reproduces Newton's force law without further input, and a dissipation-free cross-coupling to the bulk viscous pressure that no variational principle can produce. Regularity constrains the Yukawa-type corrections by an amplitude sum rule, $\sum_iα_i=-1$, so that laboratory tests of the inverse-square law bound the largest internal length below $4\times10^{-5}$~m, and identify the complex-root regime as a classically consistent Lee--Wick-like sector.

The energy of various asymptotically AdS and asymptotically flat solutions to the Einstein equations with time-translational invariance and a timelike singularity is computed, spaces with Kasner timelike singularities and black holes with negative masses among them. The definition used for gravitational energy is the on-shell Hamiltonian of the physical classical action proposed by Hawking and Horowitz thirty years ago and which agrees with the Abbott-Deser and ADM expressions. To deal with the singularity, a timelike regulating cut-off boundary surrounding the singularity compatible with a well-defined foliation in constant-time slices is introduced, so calculations are performed at a finite value of the regulating parameter, which is set to zero in the final expression for the energy. In all instances considered, the energy contribution from the singularity is finite and positive, its value depends on the geometry around the singularity and is such that the total gravitational energy is positive. Of the solutions reviewed, only the AdS soliton (which has no singularity) has negative gravitational energy relative to AdS space. The AdS and the Schwarzschild black holes with negative masses have respectively the same energies as AdS and Minkowski spaces.

SV Per, RW Cam, and KN Cen are Galactic classical Cepheids found in binary systems along with B-type companions. Each system exhibits an ``apparent age discrepancy", where the Cepheid appears significantly younger than its companion. We compute stellar evolution tracks using the MESA code and compare them with these stars on the Hertzsprung-Russell diagram to estimate the masses and the corresponding ages. We find that, for two of the systems, the companions appear to be between 1.5 times to double the age of the Cepheids, suggesting that the Cepheids may be rejuvenated merger products of two main-sequence stars that have since evolved. We find that the third system KN Cen is consistent with single star evolution when we reevaluate ultraviolet observations and find that the companion is hotter than previously measured. We test this idea using the MESA code for each system by evolving a single star model and rapidly accreting mass on the main sequence to emulate a merger, and then continuing the evolution until the model is consistent with the parameters of the Cepheids. The results are consistent with the hypothesis that the Cepheids were originally merger products. We conclude with a discussion of potential observational tests of the merger hypothesis.

The use of agentic large language models obviates human interpretation of scientific results, and will lead to substantial distrust in the literature.