2024-08-09 | | Total: 62

We present a near infrared (NIR) candidate star cluster catalog for the central kiloparsec of M82 based on new JWST NIRCam images. We identify star cluster candidates using the F250M filter, finding 1357 star cluster candidates with stellar masses $>10^4$ M$_\odot$. Compared to previous optical catalogs, nearly all (87%) of the candidates we identify are new. The star cluster candidates have a median intrinsic cluster radius of $\approx$1 pc and have stellar masses up to $10^6$ M$_\odot$. By comparing the color-color diagram to dust-free yggdrasil stellar population models, we estimate that the star cluster candidates have A$_{\rm V}\sim3-24$ mag, corresponding to A$_{\rm 2.5\mu m}\sim0.3-2.1$ mag. There is still appreciable dust extinction towards these clusters into the NIR. We measure the stellar masses of the star cluster candidates, assuming ages of 0 and 8 Myr. The slope of the resulting cluster mass function is $\beta=1.9\pm0.2$, in excellent agreement with studies of star clusters in other galaxies.

The Hubble tension is inherently multidimensional, and bears important implications for parameters beyond $H_0$. We discuss the key role of the matter density parameter $\Omega_m$ and the physical cold dark matter density $\omega_c$. We argue that once $\Omega_m$ and the physical baryon density $\omega_b$ are calibrated, through Baryon Acoustic Oscillations (BAO) and/or Type Ia Supernovae (SNeIa) for $\Omega_m$, and via Big Bang Nucleosynthesis for $\omega_b$, any model raising $H_0$ inevitably raises $\omega_c$ and, under minimal assumptions, the clustering parameter $S_8$. We explicitly verify that this behaviour holds when analyzing recent BAO and SNeIa data. We argue that a calibration of $\Omega_m$ which is as reliable and model-independent as possible should be a priority in the Hubble tension discussion, and an interesting possibility in this sense could be represented by galaxy cluster gas mass fraction measurements.

Over the last decade, much attention has been given to the study of modified gravity theories to find a more natural explanation for the late-time acceleration of the Universe. Particular attention has focused on the so-called $f(R)$ dark energy models. Instead of focusing on a particular f(R) model, we present a completely model-independent approach to study the background dynamics and the growth of matter density perturbations for those f(R) models that mimic the $\Lambda$CDM evolution at the background level. We do this by characterising the dynamics of the gravitational field using a set of dimensionless variables and using cosmography to determine the expansion history. We then illustrate the integrity of this method by fixing the cosmography to be the same as an exact $\Lambda$CDM model, allowing us to test the solution. We compare the exact evolution of the density contrast and growth index with what one obtains from various levels of the quasi-static approximation, without choosing the form of $f(R)$ dark energy.

We show how infinite derivative modifications of gravity impact on the stochastic background of Gravitational Waves from early Universe. The generic property of the ghost-free theory fixed on Minkowski space-time is the emergence of an infinite number of complex mass states when other classical backgrounds are considered. These additional states are shown to enhance the power spectrum of scalar perturbations generated during inflation. Current and future space-based and terrestrial interferometers offer indirect testing methods for the infinite derivative gravity action, enabling the exploration of new parameter spaces. In particular, we identify unconventional blue-tilted Gravitational Wave spectra, presenting a novel approach for testing infinite derivative quantum gravity in the future.

In this paper we do three things. First, we outline the conditions under which the interaction rate of inelastic processes that change the internal state of a system of $N$ targets scales as $N^2$. This is an effect distinct from coherent elastic scattering, but with the same scaling. Second, we compute example rates for such processes for various weakly interacting particles. Finally, we point to potential quantum observables for these processes that go beyond traditional energy exchange. Maximal coherence in inelastic processes is achieved when the targets are placed in an equal superposition of the ground and excited states. These coherent inelastic processes are analogous to Dicke superradiance, and we thus refer to them as superradiant interactions. We compute the superradiant interaction rates for the C$\nu$B, dark matter scattering and absorption, and late-universe particles, such as reactor neutrinos, when the two-level system is realized by nuclear or electron spins in a magnetic field. The rates we find can be quite sizable on macroscopic yet small targets. For example, the C$\nu$B interacts with a rate of $\mathcal{O}(\text{Hz})$ when scattering off a 10~cm liquid or solid-state density spin-polarized sphere. For QCD axion dark matter, similar rates can be achieved with much smaller samples, $N \sim \mathcal{O}(10^{15})\left(\frac{m}{2\times 10^{-8}~\text{eV}}\right)^{-1/2}$, where $m$ is the axion mass. Using the Lindblad formalism for open quantum systems, we show that these superradiant interactions can manifest as a source of noise on the system. This points to new observables, sensitive to the sum of the excitation and de-excitation rates, and can be viewed as introducing diffusion and decoherence to the system. The effects presented in this paper may point to a new class of ultra-low threshold detectors.

The particle contact model is important for powder simulations. Although several contact models have been proposed, their validity has not yet been well established. Therefore, we perform molecular dynamics (MD) simulations to clarify the particle interaction. We simulate head-on collisions of two particles with impact velocities less than a few percent of the sound velocity to investigate the dependence of the interparticle force and the coefficient of restitution (COR) on the impact velocity and particle radius. In this study, we treat particles with a radius of 10-100 nm and perform simulations. We find that the interparticle force exhibits hysteresis between the loading and unloading phases. Larger impact velocities result in strong hysteresis and plastic deformation. For all impact velocities and particle radii, the coefficient of restitution is smaller than that given by the Johnson-Kendall-Robert theory. An inelastic contact model cannot reproduce our MD simulations. In particular, the COR is significantly reduced when the impact velocity exceeds a certain value. This significant energy dissipation cannot be explained even by the contact models including plastic deformation. We also find that the COR increases with increasing particle radius. We also find that the previous contact models including plastic deformation cannot explain the strong energy dissipation obtained in our MD simulations, although they agree with the MD results for very low impact velocities. Accordingly, we have constructed a new dissipative contact model in which the dissipative force increases with the stress generated by collisions. The new stress dependent model successfully reproduces our MD results over a wider range of impact velocities than the conventional models do. In addition, we proposed another, simpler, dissipative contact model that can also reproduce the MD results.

We study the nonlinear theory of small-amplitude electron-acoustic solitons (EASs) in a relativistic astrophysical magnetoplasma consisting of two-temperature electrons: a sparse population of relativistic nondegenerate classical electrons and a group of fully degenerate dense relativistic electrons (main constituent) immersed in a static magnetic field with a neutralizing stationary ion background. By using the multiple-scale reductive perturbation technique with the Lorentz transformation, the Zakharov-Kuznetsov (ZK) and the modified Zakharov-Kuznetsov (mZK) equations are derived to describe the evolution of EASs in two different regimes of relativistic degeneracy: $r_{d0}<50$ and $r_{d0}\gtrsim50$. The characteristics of the plane soliton solutions of ZK and mZK equations and the soliton energy are studied. We show that the solitons moving at an angle $\alpha$ to the external magnetic field can be unstable under transverse long-wavelength perturbations. The growth rates of instabilities are obtained and analyzed with the effects of the relativity parameter $\beta_{\rm{cl}}=k_BT_{\rm{cl}}/m_ec^2$ and the degeneracy parameter $r_{d0}$, where $k_B$ is the Boltzmann constant and $T_{\rm{cl}}$ is the temperature of classical electrons. Interestingly, the ZK solitons, even if it is stable for the first-order perturbations, can be unstable in the second-order correction. Furthermore, while the first-order growth rates of perturbations for ZK solitons tend to vanish as $\alpha\rightarrow 38^\circ$, that for the mZK soliton goes to zero as $\alpha\rightarrow 90^\circ$. However, depending on the angle $\alpha$, the growth rates are found to be reduced either by increasing the values of $\beta_{\rm{cl}}$ or by decreasing the values of $r_{d0}$. The applications of our results to astrophysical plasmas, such as those in the environments of white dwarfs are discussed.

Loop quantum cosmology has achieved great successes, in which the polymerization plays a crucial role. In particular, the phase-space-variable dependent polymerization turns out to be the unique one that leads to consistent quantization of the homogeneous and isotropic universe. However, when applying the same scheme to the quantization of black holes, it meets resistances. In this paper, we continue to study the quantum effects of the polymerization near the location that a classical black hole horizon used to exist, from the point of view of effective loop quantum gravity. In particular, we find a phase-space-variable dependent polymerization scheme that leads to negligible quantum effects near the location of the classical black hole horizon, but significantly alters the spacetime structure near the origin, so that the classical singularity is finally replaced by a finite and regular transition surface. The final geodesically-complete spacetime consists of a regular transition surface that connects a black hole in one side and an anti-trapped region in the other side. In the anti-trapped region, no white hole horizons are found and the spacetime is extended to infinity, at which the geometric radius of the two-spheres becomes infinitely large.

Following the work of H. Ellis, we study warp drives in the gravitational field of a Schwarzschild black hole. We find that as long as the warp drive crosses the black hole horizon at a subluminal speed, the horizon would be effectively absent inside the warp bubble. Moreover, we discover that the black hole's gravitational field can alleviate the violations of the weak energy condition (WEC) and the null energy condition (NEC) and therefore decrease the amount of negative energy required to sustain a warp drive, which may be instrumental for creating microscopic warp drives in lab experiments. We also consider the thermodynamics of a warp bubble interacting with a black hole and point out some paradoxes that may indicate a gap in our understanding of them from the thermodynamic point of view.

In the strong magnetic field of a neutron star's magnetosphere, axions coupled to electromagnetism develop a nonzero probability to convert into photons. Past studies have revealed that the axion-photon conversion can be resonantly enhanced. We recognize that the axion-photon resonance admits two parametrically distinct resonant solutions, which we call the mass-matched resonance and the Euler-Heisenberg assisted resonance. The mass-matched resonance occurs at a point in the magnetosphere where the radially-varying plasma frequency crosses the axion mass $\omega_\mathrm{pl} \approx m_a$. The Euler-Heisenberg assisted resonance occurs where the axion energy satisfies $\omega \approx (2 \omega_\mathrm{pl}^2 / 7 g_{\gamma\gamma\gamma\gamma} \bar{B}^2 )^{1/2}$. This second resonance is made possible though the strong background magnetic field $\bar{B}$ as well as the nonzero Euler-Heisenberg four-photon self interaction, which has the coupling $g_{\gamma\gamma\gamma\gamma} = 8 \alpha^2 / 45 m_e^4$. We study the resonant conversion of relativistic axion dark radiation into photons via the Euler-Heisenberg assisted resonance, and we calculate the expected electromagnetic radiation assuming different values for the axion-photon coupling $g_{a\gamma\gamma}$ and different amplitudes for the axion flux onto the neutron star $\Phi_a$. We briefly discuss several possible sources of axion dark radiation. Achieving a sufficiently strong axion flux to induce a detectable electromagnetic signal seems unlikely.

We study the gravitational wave (GW) phase shift arising from R{\o}mer delay in binary black hole (BBH) mergers formed dynamically in three-body systems, where both the inner orbit of the merging binary and the outer orbit are eccentric. We provide a semi-analytical model and several analytical approximations that allow for fast evaluation of both the temporal evolution and the maximum value of the phase shift. The highest phase shifts occur when the binary merges close to the pericentre of the outer orbit, and can in this case be orders-of-magnitude larger compared to the circular limit. At high outer orbit eccentricities, the orbital curvature leaves distinct imprints onto the phase shift if the binary passes the outer pericentre during its inspiral. By comparing with phase-shifts measured in numerical chaotic 3-body scatterings, we show that our model accurately describes the observed phase of dynamically assembled binary systems in realistic astrophysical scenarios, providing a way to directly determine their formation channel via single GW observations.

Coronal mass ejections (CMEs) are magnetized plasma systems with highly complex magnetic topology and evolution. Methods developed to assess their magnetic configuration have primarily focused on reconstructing three-dimensional representations from one-dimensional time series measurements taken in situ using techniques based on the "highly twisted magnetic flux rope" approximations. However, the magnetic fields of CMEs is know to have more complicated geometries. Their structure can be quantified using measures of field line topology, which have been primarily used for solar physics research. In this work, we introduce a novel technique of directly quantifying the various form of magnetic helicity within a CME in the interplanetary space using synthetic in situ measurements. We use a relatively simple three-dimensional simulation of a CME initiated with a highly-twisted flux rope. We find that a significant portion of the magnetic helicity near 1~au is contained in writhe and mutual helicity rather than just in twist. We discuss the implications of this finding for fitting and reconstruction techniques.

The Peccei-Quinn (PQ) symmetry does not need to be exact, and even a tiny source of PQ breaking not aligned with the QCD anomaly might have significant phenomenological implications. In this study, we examine the effects of a general class of PQ-breaking operators on the axion cosmological production via misalignment, focussing on both temperature-dependent and independent PQ-breaking potentials. In particular, we show that a variant of the trapped misalignment mechanism can delay the onset of axion oscillation, leading to an axion dark matter window with $m_a \gg 10^{-5}$ eV. This scenario is testable through various experimental approaches, including standard axion haloscopes and helioscopes, as well as searches for electric dipole moments and axion-mediated forces.

We present NIRSpec/MSA observations from the JWST large-area survey WIDE, targeting the rest-frame UV-optical spectrum of Ulema, a radio-AGN host at redshift z=4.6348. The low-resolution prism spectrum displays high equivalent width nebular emission, with remarkably high ratios of low-ionisation species of oxygen, nitrogen and sulphur, relative to hydrogen; auroral O$^+$ emission is clearly detected, possibly also C$^+$. From the high-resolution grating spectrum, we measure a gas velocity dispersion $\sigma$~400 km s$^{-1}$, broad enough to rule out star-forming gas in equilibrium in the gravitational potential of the galaxy. Emission-line ratio diagnostics suggest that the nebular emission is due to a shock which ran out of pre-shock gas. To infer the physical properties of the system, we model simultaneously the galaxy spectral energy distribution (SED) and shock-driven line emission under a Bayesian framework. We find a relatively low-mass, star-forming system (M* = 1.4$\times$10^{10} M$_\odot$, SFR = 70 M$_\odot$ yr$^{-1}$), where shock-driven emission contributes 50 per cent to the total H$\beta$ luminosity. The nebular metallicity is near solar - three times higher than that predicted by the mass-metallicity relation at z=4.6, possibly related to fast-paced chemical evolution near the galaxy nucleus. We find no evidence for a recent decline in the SFR of the galaxy, meaning that, already at this early epoch, fast radio-mode AGN feedback was poorly coupled with the bulk of the star-forming gas; therefore, most of the feedback energy must end up in the galaxy halo, setting the stage for future quenching.

It has been suggested that giant planet occurrence peaks for stars with $M_{\ast}~\approx~3~M_{\odot}$ at a value a factor of four higher than observed for solar-mass stars. This population of giant planets predicted to frequently orbit main sequence B stars at $a~\approx~10$ AU is difficult to characterize observationally while fusion persists in their host stars. Fortunately, stars with $M_{\ast}~\gtrsim~3~M_{\odot}$ sustain fusion for only a few hundred million years. By the time those stars become massive, young white dwarfs, any giant planets present would still be luminous as a consequence of their recent formation. From an initial sample of 3268 Gaia-identified massive, young white dwarfs, we use homogeneous Spitzer Infrared Array Camera (IRAC) photometry to search for evidence of unresolved giant planets. For 30 systems, these IRAC data provide sensitivity to objects with $M~\lesssim~10~M_{\text{Jup}}$, and we identify one candidate with $M~\approx~4~M_{\text{Jup}}$ orbiting the white dwarf GALEX J071816.4+373139. Correcting for the possibility that some of the white dwarfs in our sample result from mergers, we find a giant planet occurrence $\eta_{\text{GP}}~=~0.11_{-0.07}^{+0.13}$ for stars with initial masses $M_{\ast}~\gtrsim~3~M_{\odot}$. Our occurrence inference is consistent with both the Doppler-inferred occurrence of giant planets orbiting $M_{\ast}~\approx~2~M_{\odot}$ giant stars and the theoretically predicted factor-of-four enhancement in the occurrence of giant planets orbiting $M_{\ast}~\approx~3~M_{\odot}$ stars relative to solar-mass stars. Future James Webb Space Telescope Near Infrared Camera observations of our sample would provide sensitivity to Saturn-mass planets and thereby a definitive estimate of the occurrence of giant planets orbiting stars with $M_{\ast}~\gtrsim~3~M_{\odot}$.

Galaxy clusters show large-scale azimuthal X-ray surface brightness fluctuations known as cold fronts. These are overdense (average density jumps $\sim 30\%$ or post-jump density $\sim 130\%$) and have milder discontinuity in pressure. Cold fronts are argued to originate due to sloshing driven by sub-halo passage at close proximity to the cluster center. While this is a viable source of large-scale perturbations, the physical mechanisms that can sustain such density structures (of specific geometry) are not clear. In this work, we explore whether long wavelength thermal instability is an explanation for cold front formation in a cluster core which is perturbed by sub-halos or AGN activity. Using global linear perturbation analysis, we show that internal gravity waves (thermally unstable) can form large-scale three-dimensional spiral structures, akin to observed cold fronts. We explore if the presence of magnetic field (along spherical $\hat{\phi}$) may support such structures (by suppressing small scale Kelvin-Helmholtz modes) or disrupt them (by promoting additional thermal instability). We find that latter happens at shorter wavelengths and only at frequencies above the characteristic buoyancy or Brunt V\"ais\"al\"a frequency ($>N_{\rm BV}$). Our work implies, firstly, that large-scale spirals may be formed and sustained over a long timescale ($>N^{-1}_{\rm BV}$) even in presence of aligned magnetic fields that is otherwise supportive against mixing at the interface. Secondly, short-wavelength (but relatively longer along the field) unstable compressive modes may form within or in the vicinity of such spirals. The instability is an overstable slow wave, and grows in 2D at timescales $\gtrsim 2-3$ times longer than the spiral growth timescale (via thermal instability). Thus we claim that this instability cannot destroy the large scale coherence.

The two main competing theories proposed to explain the formation of massive ($>10$M$_\odot$) stars -- competitive accretion and monolithic core collapse -- make different observable predictions for the environment of the massive stars during, and immediately after, their formation. Proponents of competitive accretion have long predicted that the most massive stars should have a different spatial distribution to lower-mass stars, either through the stars being mass segregated, or being in areas of higher relative densities, or sitting deeper in gravitational potential wells. We test these predictions by analysing a suite of SPH simulations where star clusters form massive stars via competitive accretion with and without feedback. We find that the most massive stars have higher relative densities, and sit in deeper potential wells, only in simulations in which feedback is not present. When feedback is included, only half of the simulations have the massive stars residing in deeper potential wells, and there are no other distinguishing signals in their spatial distributions. Intriguingly, in our simple models for monolithic core collapse, the massive stars may also end up in deeper potential wells, because if massive cores fragment the stars are still massive, and dominate their local environs. We find no robust diagnostic test in the spatial distributions of massive stars that can distinguish their formation mechanisms, and so other predictions for distinguishing between competitive accretion and monolithic collapse are required.

SN 2023ixf was discovered in M101 within a day of explosion and rapidly classified as a Type II supernova with flash features. Here we present ultraviolet (UV) spectra obtained with the Hubble Space Telescope 14, 19, 24, and 66 days after explosion. Interaction between the supernova ejecta and circumstellar material (CSM) is seen in the UV throughout our observations in the flux of the first three epochs and asymmetric MgII emission at day 66. We compare our observations to CMFGEN supernova models which include CSM interaction ($\dot{M}<10^{-3}$ $M_{\odot}$/yr) and find that the power from CSM interaction is decreasing with time, from $L_{\rm sh}\approx5\times10^{42}$ erg/s at day 14 to $L_{sh}\approx1\times10^{40}$ erg/s at day 66. We examine the contribution of individual atomic species to the spectra at day 14 and 19, showing that the majority of the features are dominated by iron, nickel, magnesium, and chromium absorption in the ejecta. The UV spectral energy distribution of SN 2023ixf sits between that of supernovae which show no definitive signs of CSM interaction and those with persistent signatures assuming the same progenitor radius and metallicity. Finally, we show that the evolution and asymmetric shape of the MgII emission are not unique to SN 2023ixf. These observations add to the early measurements of dense, confined CSM interaction, tracing the mass-loss history of SN 2023ixf to $\sim33$ yr prior to explosion and the density profile to a radius of $\sim5.7\times10^{15}$ cm. They show the relatively short evolution from quiescent red supergiant wind to high mass loss.

Recent searches for parity breaking in the galaxy four-point correlation function, as well as the prospects for greatly improved sensitivity to parity breaking in forthcoming surveys, motivate the search for physical mechanisms that could produce such a signal. Here we show that a parity-violating galaxy four-point correlation function may be induced by lensing by a chiral gravitational-wave background. We estimate the amplitude of a signal that would be detectable with a current galaxy survey, taking into account constraints to the primordial gravitational-wave-background amplitude. We find that this mechanism is unlikely to produce a signal large enough to be seen with a galaxy survey but note that it may come within reach with future 21cm observations.

Since the universe is not transparent to gamma rays with energies above around one hundred GeV, it is necessary to account for the interaction of high-energy photons with intergalactic radiation fields in order to model gamma-ray propagation. Here, we present a public numerical software for the modeling of gamma-ray observables. This code computes the effects on gamma-ray spectra from the development of electromagnetic cascades and cosmological redshifting. The code introduced here is based on the original $\gamma$-Cascade, and builds on it by improving its performance at high redshifts, introducing new propagation modules, and adding many more extragalactic radiation field models, which enables the ability to estimate the uncertainties inherent to EBL modeling. We compare the results of this new code to existing electromagnetic transport models.

Galaxy scaling relations provide insights into the processes that drive galaxy evolution. The extension of these scaling relations into the dwarf galaxy regime is of particular interest. This is because dwarf galaxies represent a crucial stage in galaxy evolution, and understanding them could also shed light on their role in reionising the early Universe. There is currently no consensus on the processes that dominate the evolution of dwarfs. In this work we constrain the atomic gas sequence (stellar mass vs. atomic gas fraction) and mass-metallicity relation (stellar mass vs. gas phase metallicity) from dwarf ($10^{6.5}$ $\textrm{M}_{\odot}$) to massive ($10^{11.5}$ $\textrm{M}_{\odot}$) galaxies in the local Universe. The combined optical and 21-cm spectroscopic observations of the DESI and ALFALFA surveys allow us to simultaneously constrain both scaling relations. We find a slope change of the atomic gas sequence at a stellar mass of $\sim 10^{9} ~\textrm{M}_{\odot}$. We also find that the shape and scatter of the atomic gas sequence and mass-metallicity relation are strongly linked for both dwarfs and more massive galaxies. Consequently, the low mass slope change of the atomic gas sequence is imprinted onto the mass-metallicity relation of dwarf galaxies. The mass scale of the measured slope change is consistent with a predicted escape velocity threshold below which low mass galaxies experience significant supernova-driven gas loss, as well as with a reduction in cold gas accretion onto more massive galaxies.

Carbonyl sulfide (OCS) is widely observed in the gas phase towards star-forming regions and was the first of the only two sulfur-bearing species detected in interstellar ices so far. However, the chemical network governing its formation is still not fully understood. While the sulfurization of CO and the oxidation of CS are often invoked to form OCS, other mechanisms could have a significant contribution. In particular, the multistep reaction involving CO and SH is a good candidate to forming OCS in dense cloud environments. We aim to constrain the viability of the CO + SH route to forming solid OCS in the interstellar medium, in a similar manner as CO + OH is known to produce CO2 ice. This is achieved by conducting a systematic laboratory investigation of the targeted reactions on interstellar ice analogues under dense cloud conditions. An ultrahigh vacuum chamber is utilized to simultaneously deposit CO, H2S, and atomic H at 10 K. SH radicals produced in situ via hydrogen abstraction from H2S react with CO to form OCS. OCS is efficiently formed through surface reactions involving CO, H2S, and H atoms. The suggested underlying mechanism behind OCS formation is CO + SH -> HSCO followed by HSCO + H -> OCS + H2. The OCS yield reduces slowly, but remains significant with increasing CO:H2S mixing ratios (CO:H2S = 1:1, 5:1, 10:1, and 20:1). Our experiments provide unambiguous evidence that OCS can be formed from CO + SH in the presence of H atoms. This route remains efficient for large H2S dilutions (5% w.r.t CO), suggesting that it is a viable mechanism in interstellar ices. Given that SH radicals can be created in clouds throughout a wide evolutionary timescale, this mechanism could have a non-negligible contribution to forming interstellar OCS ice.

The very long-term evolution of the hierarchical restricted three-body problem with a massive perturber is analyzed analytically in the high eccentricity regime. Perturbations on the time scale of the outer orbit can accumulate over long timescales and be comparable to the effect of the octupole term. These perturbations are described by Brown's Hamiltonian - having different forms in the literature. We show that at the high eccentricity regime - the effect of Brown's Hamiltonian is an azimuthal precesssion of the eccentricity vector and can be solved analytically. In fact, the dynamics are equivalent to a simple pendulum model allowing an explicit flip criterion.

Crush curves are of fundamental importance to numerical modeling of small and porous astrophysical bodies. The empirical literature often measures them for silica grains, and different studies have used various methods, sizes, textures, and pressure conditions. Here we review past studies and supplement further experiments in order to develop a full and overarching understanding of the silica crush curve behavior. We suggest a new power-law function that can be used in impact simulations of analog materials similar to micro-granular silica. We perform a benchmarking study to compare this new crush curve to the parametric quadratic crush curve often used in other studies, based on the study case of the DART impact onto the asteroid Dimorphos. We find that the typical quadratic crush curve parameters do not closely follow the silica crushing experiments, and as a consequence they under (over) estimate compression close (far) from the impact site. The new crush curve presented here, applicable to pressures between a few hundred Pa and up to 1.1 GPa, might therefore be more precise. Additionally, it is not calibrated by case-specific parameters, and can be used universally for comet- or asteroid-like bodies, given an assumed composition similar to micro-granular silica.

The R-value is a measure of the strength of photospheric magnetic Polarity Inversion Lines (PILs) in Active Regions (ARs). This work investigates the possibility of a relation between R-value variations and the occurrence of X-class flares in ARs, not in the solar photosphere, as usual, but above it in regions, closer to where flares occur. The modus operandi is to extrapolate the Solar Dynamic Observatory's (SDO) Helioseismic and Magnetic Imager (HMI) magnetogram data up to a height of 3.24 Mm above the photosphere and then compute the R-value based on the extrapolated magnetic field. Recent studies have shown that certain flare-predictive parameters such as the horizontal gradient of the vertical magnetic field and magnetic helicity may improve flare prediction lead times significantly if studied at a specific height range above the photosphere, called the Optimal Height Range (OHR). Here we define the OHR as a collection of heights where a sudden but sustained increase in R-value is found. For the eight case studies discussed in this paper, our results indicate that it is possible for OHRs to exist in the low solar atmosphere (between 0.36 - 3.24 Mm), where R-value spikes occur 48-68 hrs before the first X-class flare of an emerging AR. The temporal evolution of R-value before the first X-class flare for an emerging AR is also found to be distinct from that of non-flaring ARs. For X-class flares associated with non-emerging ARs, an OHR could not be found.