Search

Identified as parsec-size, gas clumps at the junction of multiple filaments, hub-filament systems (HFS) play a crucial role during the formation of young clusters and high-mass stars. These HFS appear nevertheless to be detached from most galactic filaments when compared in the mass-length (M-L) phase-space. We aim to characterize the early evolution of HFS as part of the filamentary description of the interstellar medium. Combining previous scaling relations with new analytic calculations, we created a toy model to explore the different physical regimes described by the M-L diagram. Despite its simplicity, our model accurately reproduces several observational properties reported for filaments and HFS such as their expected typical aspect ratio ($A$), mean surface density ($\Sigma$), and gas accretion rate ($\dot{m}$). Moreover, this model naturally explains the different mass and length regimes populated by filaments and HFS, respectively. Our model predicts a dichotomy between filamentary ($A\geq 3$) and spheroidal ($A<3$) structures connected to the relative importance of their fragmentation, accretion, and collapse timescales. Individual filaments with low accretion rates are dominated by an efficient internal fragmentation. In contrast, the formation of compact HFS at the intersection of filaments triggers a geometric phase-transition leading to the gravitational collapse of these structures at parsec-scales in $\sim$1Myr also inducing higher accretion rates., Comment: Accepted by A&A, 17 pages, 11 figures

The reported discovery of a cold ($\sim$10$^4~\text{K}$) disc-like structure around the super-massive black hole at the centre of the Milk Way, Sagittarius A* (Sgr A*), has challenged our understanding of the gas dynamics and thermodynamic state of the plasma in its immediate vicinity. State-of-the-art simulations do not agree on whether or not such a disc can indeed be a product of the multiple stellar wind interactions of the mass-losing stars in the region. This study aims to constrain the conditions for the formation of a cold disc as a natural outcome of the system of the mass-losing stars orbiting around Sgr A*, to investigate if the disc is a transient or long-lasting structure, and to assess the validity of the model through direct comparisons with observations. We conduct a set of hydrodynamic simulations of the observed Wolf-Rayet (WR) stars feeding Sgr A* using the finite-volume adaptive mesh-refinement code Ramses. We focus, for the first time, on the impact of the chemical composition of the plasma emanating from the WR stars. The simulations show that the chemical composition of the plasma affects the radiative cooling enough to impact the properties of the medium such as density and temperature and, as a consequence, the rate at which the material inflows onto Sgr A*. We demonstrated that the formation of a cold disc from the stellar winds is possible for certain chemical compositions that are consistent with the current observational constraints. However, even in such a case, it is not possible to reproduce the reported properties of the observed disc-like structure, namely its inclination and hydrogen recombination line fluxes. We conclude that the stellar winds on their own cannot form the cold disc around Sgr A* inferred from the observations. Either relevant ingredients are still missing in the model, or the interpretation of the observed data needs to be revised., Comment: Submitted to A&A, 17 pages, 11 figures (+2 pages, +2 figures in Appendix)

The spin structure functions of the proton and the deuteron were measured during the EG4 experiment at Jefferson Lab in 2006. Data were collected for longitudinally polarized electron scattering off longitudinally polarized NH$_3$ and ND$_3$ targets, for $Q^2$ values as small as 0.012 and 0.02 GeV$^2$, respectively, using the CEBAF Large Acceptance Spectrometer (CLAS). This is the archival paper of the EG4 experiment that summaries the previously reported results of the polarized structure functions $g_1$, $A_1F_1$, and their moments $\overline \Gamma_1$, $\overline \gamma_0$, and $\overline I_{TT}$, for both the proton and the deuteron. In addition, we report on new results on the neutron $g_1$ extracted by combining proton and deuteron data and correcting for Fermi smearing, and on the neutron moments $\overline \Gamma_1$, $\overline \gamma_0$, and $\overline I_{TT}$ formed directly from those of the proton and the deuteron. Our data are in good agreement with the Gerasimov-Drell-Hearn sum rule for the proton, deuteron, and neutron. Furthermore, the isovector combination was formed for $g_1$ and the Bjorken integral $\overline \Gamma_1^{p-n}$, and compared to available theoretical predictions. All of our results provide for the first time extensive tests of spin observable predictions from chiral effective field theory ($\chi$EFT) in a $Q^2$ range commensurate with the pion mass. They motivate further improvement in $\chi$EFT calculations from other approaches such as the lattice gauge method., Comment: 33 pages. 26 figures. Data table provided in supplementary material (30 pages)

If primordial black holes (PBHs) of asteroidal mass make up the entire dark matter they could be detectable through their gravitational influence in the solar system. In this work, we study the perturbations that PBHs induce on the orbits of planets. Detailed numerical simulations of the solar system, embedded in a halo of PBHs, are performed. Using the Earth-Mars distance as an observational probe, we show that the perturbations are below the current detection limits and thus PBHs are not directly constrained by solar system ephemerides. We estimate that an improvement in the measurement accuracy by more than an order of magnitude or the extraction of signals well below the noise level are required to detect the gravitational influence of PBHs in the solar system in the foreseeable future., Comment: 12 pages, 6 figures, submitted to APJ

We provide a detailed comparison between the ``gravoturbulent'' (GT) and ``global hierarchical collapse'' (GHC) models for molecular clouds and star formation, their respective interpretations of the observational data, the features they share, and suggested tests and observations to discern between them. Also, we clarify common misconceptions in recent literature about the global and hierarchical nature of the GHC scenario, and briefly discuss the evolution of some aspects of both models toward convergence. GT assumes that molecular clouds and their substructures are in approximate virial equilibrium and are in a near-stationary state, interprets the linewidth-derived nonthermal motions exclusively as turbulence, which provides additional pressure against self-gravity. Conversely, GHC assumes that most star-forming molecular clouds and their substructures are part of a continuous gravitationally-driven flow, each accreting from their parent structure. Thus, the clouds and their star formation rate evolve in time. GHC interprets nonthermal motions as a mixture of infall and turbulent components, with the relative importance of the former increasing as the objects become denser and/or more massive. Tests that may provide clues to distinguishing between GT and GHC must take into account that the innermost parts of globally gravitationally bound structures may not locally appear bound, and thus the binding may have to be searched for at the largest scale of the structure., Comment: Submitted to MNRAS. Comments welcome. Animation for Fig. 3 available at https://www.irya.unam.mx/gente/g.gomez/publica/24-GT_vs_GHC-f3.mpg

Understanding star formation in galaxies requires resolving the physical scale on which star formation often occurs: the scale of star clusters. We present a multiwavelength, eight-parsec resolution study of star formation in the circumnuclear star cluster and molecular gas rings of the early-type spiral NGC 1386. The cluster ring formed simultaneously ~ 4 Myr ago. The clusters have similar properties in terms of mass and star formation rate, resembling those of H II regions in the Milky Way disc. The molecular CO gas resolves into long filaments, which define a secondary ring detached from the cluster ring. Most clusters are in CO voids. Their separation with respect the CO filaments is reminiscent of that seen in galaxy spiral arms. By analogy, we propose that a density wave through the disc of this galaxy may have produced this gap in the central kpc. The CO filaments fragment into strings of dense, unresolved clouds with no evidence of a stellar counterpart. These clouds may be the sites of a future population of clusters in the ring. The free-fall time of these clouds, ~ 10 Myr, is close to the orbital time of the CO ring. This coincidence could lead to a synchronous bursting ring, as is the case for the current ring. The inward spiralling morphology of the CO filaments and co-spatiality with equivalent kpc-scale dust filaments are suggestive of their role as matter carriers from the galaxy outskirts to feed the molecular ring and a moderate active nucleus., Comment: 22 pages, 9 figures. Accepted for publication in MNRAS

Measuring Deeply Virtual Compton Scattering on the neutron is one of the necessary steps to understand the structure of the nucleon in terms of Generalized Parton Distributions (GPDs). Neutron targets play a complementary role to transversely polarized proton targets in the determination of the GPD $E$. This poorly known and poorly constrained GPD is essential to obtain the contribution of the quarks' angular momentum to the spin of the nucleon. DVCS on the neutron was measured for the first time selecting the exclusive final state by detecting the neutron, using the Jefferson Lab longitudinally polarized electron beam, with energies up to 10.6 GeV, and the CLAS12 detector. The extracted beam-spin asymmetries, combined with DVCS observables measured on the proton, allow a clean quark-flavor separation of the imaginary parts of the GPDs $H$ and $E$., Comment: 7 pages, 6 figures

Recent hydrodynamical simulations of the late stages of supernova remnant (SNR) evolution have revealed that as they merge with the ambient medium, SNRs implode, leading to the formation of dense clouds in their center. While being highly chemically enriched by their host SNR, these clouds appear to have similar properties as giant molecular clouds, which are believed to be the main site of star formation. Here, we develop a simple model, in order to estimate the efficiency of the star formation that might be triggered by the implosion of SNRs. We separately consider two cases, cyclic star formation, maintained by the episodic driving of feedback from new generations of stars; and a single burst of star formation, triggered by a single explosion. We find that in the cyclic case, star formation is inefficient, with implosion-triggered star-formation contributing a few percent of the observed star-formation efficiency per free-fall timescale. In the single-burst case, higher star-formation efficiencies can be obtained. However, while the implosion-triggered process might not contribute much to the overall star-formation, due to the high chemical enrichment of the birth clouds, it can explain the formation of a significant fraction of metal-rich stars., Comment: 9 pages, 3 figures, accepted for publication in ApJL

(Abridged) Recent molecular surveys have revealed a rich gas organization of sonic-like fibers in all kind of environments prior to the formation of low- and high-mass stars. This paper introduces the EMERGE project aiming to investigate whether complex fiber arrangements could explain the origin of high-mass stars and clusters. We analyzed the EMERGE Early ALMA Survey including 7 star-forming regions in Orion (OMC-1/2/3/4 South, L1641N, NGC2023, and Flame Nebula) homogeneously surveyed in both molecular lines (N$_2$H$^+$ J=1-0, HNC J=1-0, plus HC3N J=10-9) and 3mm-continuum using a combination of interferometric ALMA mosaics and IRAM-30m single-dish (SD) maps. Based on our low-resolution (SD) observations, we describe the global properties of our sample covering a wide range of physical conditions including low-, intermediate, and high-mass star-forming regions in different evolutionary stages. Their comparison with ancillary YSO catalogs denotes N$_2$H$^+$ as the best proxy for the dense, star-forming gas in our targets showing a constant star formation efficiency and a fast time evolution of <1 Myr. While apparently clumpy and filamentary in our SD data, all targets show a much more complex fibrous substructure at the enhanced resolution of our ALMA+IRAM-30m maps. A large number of filamentary features at sub-parsec scales are clearly recognized in the high-density gas traced by N$_2$H$^+$ directly connected to the formation of individual protostars. This complex gas organization appears to extend further into the more diffuse gas traced by HNC. This paper presents the EMERGE Early ALMA survey including a first data release of continuum maps and spectral products for this project to be analysed in future papers of this series. A first look at these results illustrates the need of advanced data combination techniques to investigate the intrinsic multi-scale, gas structure of the ISM., Comment: 37 pages, 28 fgures, Accepted by A&A

The energy injection through Hawking evaporation has been used to put strong constraints on primordial black holes as a dark matter candidate at masses below $10^{17}\,\rm{g}$. However, Hawking's semiclassical approximation breaks down at latest after half-decay. Beyond this point, the evaporation could be significantly suppressed, as was shown in recent work. In this study, we review existing cosmological and astrophysical bounds on primordial black holes, taking this effect into account. We show that the constraints disappear completely for a reasonable range of parameters, which opens a new window below $10^{10}\,\rm{g}$ for light primordial black holes as a dark matter candidate., Comment: 9 pages, 5 figures, published in MNRAS

We investigate a possible close encounter between M33 and M31 in the past to understand the role of galaxy-galaxy interactions in shaping the matter distribution in galaxy outskirts. We recovered possible orbital trajectories of M33, M31 and the Milky Way in the past, which are compatible with the Early Third Data Release of the Gaia mission and with mass estimates of Local Group spirals, after tuning mass losses and the dynamical friction term with the help of N-body numerical simulations. A close encounter of M33 and M31 in the past has a low but non-negligible probability. If the two galaxies had been closer in the past, their minimum distance would be of the order of 100 kpc or larger, and this happened earlier than 3 Gyr ago. During this encounter, 35-40% of the dark matter mass of M33 might have been removed from the halo due to tidal stripping. A detailed comparison of the results of test-particle simulations with the observed disk warp or with the spatial distribution of candidate dark satellites of M33 suggests that a closer passage of M33 around M31 cannot, however, be responsible for the observed morphological features. We suggest that more recent gas accretion events, possibly from a cosmic filament, might cause the misalignment of the outer disk of M33 after the rapid inner disk formation., Comment: 17 pages, 10 figures. Accepted for publication in A&A

Our Sun lies within 300 pc of the 2.7-kpc-long sinusoidal chain of dense gas clouds known as the Radcliffe Wave. The structure's wave-like shape was discovered using 3D dust mapping, but initial kinematic searches for oscillatory motion were inconclusive. Here we present evidence that the Radcliffe Wave is oscillating through the Galactic plane while also drifting radially away from the Galactic Center. We use measurements of line-of-sight velocity for 12CO and 3D velocities of young stellar clusters to show that the most massive star-forming regions spatially associated with the Radcliffe Wave (including Orion, Cepheus, North America, and Cygnus X) move as if they are part of an oscillating wave driven by the gravitational acceleration of the Galactic potential. By treating the Radcliffe Wave as a coherently oscillating structure, we can derive its motion independently of the local Galactic mass distribution, and directly measure local properties of the Galactic potential as well as the Sun's vertical oscillation period. In addition, the measured drift of the Radcliffe Wave radially outward from the Galactic Center suggests that the cluster whose supernovae ultimately created today's expanding Local Bubble may have been born in the Radcliffe Wave., Comment: Published in Nature on 20 February 2024. For data, interactive visualizations, and more information see https://sites.google.com/cfa.harvard.edu/radcliffewave

The deposition of energy and momentum by supernova explosions has been subject to numerous studies in the past few decades. However, while there has been some work that focused on the transition from the adiabatic to the radiative stage of a supernova remnant (SNR), the late radiative stage and merging with the interstellar medium (ISM) have received little attention. Here, we use three-dimensional, hydrodynamic simulations, focusing on the evolution of SNRs during the radiative phase, considering a wide range of physical explosion parameters ($n_{\text{H, ISM}} \in \left[0.1, 100\right] \text{cm}^{-3}$ and $E_{\text{SN}} \in \left[1, 14\right]\times 10^{51} \text{erg}$). We find that the radiative phase can be subdivided in four stages: A pressure driven snowplow phase during which the hot overpressurized bubble gas is evacuated and pushed into the cold shell, a momentum conserving snowplow phase which is accompanied by a broadening of the shell, an implosion phase where cold material from the back of the shell is flooding the central vacuum and a final cloud phase, during which the imploding gas is settling as a central, compact overdensity. The launching timescale for the implosion ranges from a few 100 kyr to a few Myr, while the cloud formation timescale ranges from a few to about 10 Myr. The highly chemically enriched clouds can become massive ($M_{\text{cl}} \sim 10^3 - 10^4 \, \text{M}_{\odot}$) and self-gravitating within a few Myr after their formation, providing an attractive, novel pathway for supernova induced star and planet formation in the ISM., Comment: 19 pages, 13 figures, resubmitted to ApJ after first revision, comments are welcome

Recent observations of protostellar cores suggest that most of the material in the protostellar phase is accreted along streamers. Streamers in this context are defined as velocity coherent funnels of denser material potentially connecting the large scale environment to the small scales of the forming accretion disk. Using simulations which simultaneously resolve the driving of turbulence on the filament scale as well as the collapse of the core down to protostellar disk scales, we aim to understand the effect of the turbulent velocity field on the formation of overdensities in the accretion flow. We perform a three-dimensional numerical study on a core collapse within a turbulent filament using the RAMSES code and analyse the properties of overdensities in the accretion flow. We find that overdensities are formed naturally by the initial turbulent velocity field inherited from the filament and subsequent gravitational collimation. This leads to streams which are not really filamentary but show a sheet-like morphology. Moreover, they have the same radial infall velocities as the low density material. As a main consequence of the turbulent initial condition, the mass accretion onto the disk does not follow the predictions for solid body rotation. Instead, most of the mass is funneled by the overdensities to intermediate disk radii.

Age is one of the most fundamental parameters of stars, yet it is one of the hardest to determine as it requires modelling various aspects of stellar formation and evolution. When we compare the ages derived from isochronal and dynamical traceback methods for six young stellar associations, we find a systematic discrepancy. Specifically, dynamical traceback ages are consistently younger by an average of $\langle\Delta_{\rm Age}\rangle = 5.5 \pm 1.1$ Myr. We rule out measurement errors as the cause of the age mismatch and propose that $\Delta_{\rm Age}$ indicates the time a young star remains bound to its parental cloud before moving away from its siblings. In this framework, the dynamical traceback "clock" starts when a stellar cluster or association begins to expand after expelling most of the gas, while the isochronal "clock" starts earlier when most stars form. The difference between these two age-dating techniques is a powerful tool to constraint evolutionary models, as isochronal ages cannot be younger than dynamical traceback ages. Measuring the $\Delta_{\rm Age}$ accurately and understanding its variations across different environments will provide further information on the impact of local conditions and stellar feedback on the formation and dispersal of stellar clusters.

We analyze TYPHOON long slit absorption line spectra of the starburst barred spiral galaxy NGC 1365 obtained with the Progressive Integral Step Method covering an area of 15 square kpc. Applying a population synthesis technique, we determine the spatial distribution of ages and metallicity of the young and old stellar population together with star formation rates, reddening, extinction and the ratio R$_V$ of extinction to reddening. We detect a clear indication of inside-out growth of the stellar disk beyond 3 kpc characterized by an outward increasing luminosity fraction of the young stellar population, a decreasing average age and a history of mass growth, which was finished 2 Gyrs later in the outermost disk. The metallicity of the young stellar population is clearly super solar but decreases towards larger galactocentric radii with a gradient of -0.02 dex/kpc. On the other hand, the metal content of the old population does not show a gradient and stays constant at a level roughly 0.4 dex lower than that of the young population. In the center of NGC 1365 we find a confined region where the metallicity of the young population drops dramatically and becomes lower than that of the old population. We attribute this to infall of metal poor gas and, additionally, to interrupted chemical evolution where star formation is stopped by AGN and supernova feedback and then after several Gyrs resumes with gas ejected by stellar winds from earlier generations of stars. We provide a simple model calculation as support for the latter.

Recent observations with JWST and ALMA have revealed extremely massive quiescent galaxies at redshifts of z=3 and higher, indicating both rapid onset and quenching of star formation. Using the cosmological simulation suite Magneticum Pathfinder we reproduce the observed number densities and stellar masses, with 36 quenched galaxies of stellar mass larger than 3e10Msun at z=3.42. We find that these galaxies are quenched through a rapid burst of star-formation and subsequent AGN feedback caused by a particularly isotropic collapse of surrounding gas, occurring on timescales of around 200Myr or shorter. The resulting quenched galaxies host stellar components which are kinematically fast rotating and alpha-enhanced, while exhibiting a steeper metallicity and flatter age gradient compared to galaxies of similar stellar mass. The gas of the galaxies has been metal enriched and ejected. We find that quenched galaxies do not inhabit the densest nodes, but rather sit in local underdensities. We analyze observable metrics to predict future quenching at high redshifts, finding that on shorter timescales <500Myr the ratio M_bh/M_* is the best predictor, followed by the burstiness of the preceding star-formation, t50-t90 (time to go from 50% to 90% stellar mass). On longer timescales, >1Gyr, the environment becomes the strongest predictor, followed by t50-t90, indicating that at high redshifts the consumption of old and lack of new gas are more relevant for long-term prevention of star-formation than the presence of a massive AGN. We predict that relics of such high-z quenched galaxies should best be characterized by a strong alpha enhancement., Comment: 22 pages, 13 figures, Submitted to ApJ, Comments welcome

We study the kinematics of the recently discovered Corona Australis (CrA) chain of clusters by examining the 3D space motion of its young stars using Gaia DR3 and APOGEE-2 data. While we observe linear expansion between the clusters in the Cartesian XY directions, the expansion along Z exhibits a curved pattern. To our knowledge, this is the first time such a nonlinear velocity-position relation has been observed for stellar clusters. We propose a scenario to explain our findings, in which the observed gradient is caused by stellar feedback, accelerating the gas away from the Galactic plane. A traceback analysis confirms that the CrA star formation complex was located near the central clusters of the Scorpius Centaurus (Sco-Cen) OB association 10-15 Myr ago. It contains massive stars and thus offers a natural source of feedback. Based on the velocity of the youngest unbound CrA cluster, we estimate that a median number of about two supernovae would have been sufficient to inject the present-day kinetic energy of the CrA molecular cloud. This number agrees with that of recent studies. The head-tail morphology of the CrA molecular cloud further supports the proposed feedback scenario, in which a feedback force pushed the primordial cloud from the Galactic north, leading to the current separation of 100 pc from the center of Sco-Cen. The formation of spatially and temporally well-defined star formation patterns, such as the CrA chain of clusters, is likely a common process in massive star-forming regions., Comment: Accepted for publication as a Letter in Astronomy and Astrophysics

We present the first determination of the proton mechanical radius. The result was obtained by employing a novel theoretical approach that connects experimental data of deeply virtual Compton scattering with the spin = 2 interaction that is characteristic of gravity coupling with matter. We find that the proton mechanical radius is significantly smaller than its charge radius, consistent with the latest Lattice QCD computation., Comment: 4 pages, 3 figures, 1 table

The parameterization of the nucleon structure through Generalized Parton Distributions (GPDs) shed a new light on the nucleon internal dynamics. For its direct interpretation, Deeply Virtual Compton Scattering (DVCS) is the golden channel for GPDs investigation. The DVCS process interferes with the Bethe-Heitler (BH) mechanism to constitute the leading order amplitude of the $eN \to eN\gamma$ process. The study of the $ep\gamma$ reaction with polarized positron and electron beams gives a complete set of unique observables to unravel the different contributions to the $ep \gamma$ cross section. This separates the different reaction amplitudes, providing a direct access to their real and imaginary parts which procures crucial constraints on the model dependences and associated systematic uncertainties on GPDs extraction. The real part of the BH-DVCS interference amplitude is particularly sensitive to the $D$-term which parameterizes the Gravitational Form Factors of the nucleon. The separation of the imaginary parts of the interference and DVCS amplitudes provides insights on possible higher-twist effects. We propose to measure the unpolarized and polarized Beam Charge Asymmetries (BCAs) of the $\vec{e}^{\pm}p \to e^{\pm}p \gamma$ process on an unpolarized hydrogen target with {\tt CLAS12}, using polarized positron and electron beams at 10.6 GeV. The azimuthal and $t$-dependences of the unpolarized and polarized BCAs will be measured over a large $(x_B,Q^2)$ phase space using a 100 day run with a luminosity of 0.66$\times 10^{35}$cm$^{-2}\cdot$s$^{-1}$., Comment: Proposal to the Jefferson Lab Program Advisory Committee (PAC51)

One of the most surprising discoveries made at Jefferson Lab has been the discrepancy in the determinations of the proton's form factor ratio $\mu_p G_E^p/G_M^p$ between unpolarized cross section measurements and the polarization transfer technique. Over two decades later, the discrepancy not only persists but has been confirmed at higher momentum transfers now accessible in the 12-GeV era. The leading hypothesis for the cause of this discrepancy, a non-negligible contribution from hard two-photon exchange, has neither been conclusively proven or disproven. This state of uncertainty not only clouds our knowledge of one-dimensional nucleon structure but also poses a major concern for our field's efforts to map out the three-dimensional nuclear structure. A better understanding of multi-photon exchange over a wide phase space is needed. We propose making comprehensive measurements of two-photon exchange over a wide range in momentum transfer and scattering angle using the CLAS12 detector. Specifically, we will measure the ratio of positron-proton to electron-proton elastic scattering cross sections, using the proposed positron beam upgrade for CEBAF. The experiment will use 2.2, 4.4, and 6.6 GeV lepton beams incident on the standard CLAS12 unpolarized hydrogen target. Data will be collected by the CLAS12 detector in its standard configuration, except for a modified trigger to allow the recording of events with beam leptons scattered into the CLAS12 central detector. The sign of the beam charge, as well as the polarity of the CLAS12 solenoid and toroid, will be reversed several times in order to suppress systematics associated with local detector efficiency and time-dependent detector performance. The proposed high-precision determination of two-photon effects will be..., Comment: Experimental Proposal E12+23-008 submitted to Jefferson Lab PAC 51, 34 pages, 18 figures

The inner kiloparsec regions surrounding sub-Eddington (luminosity less than 10$^{-3}$ in Eddington units, L$_{Edd}$) supermassive black holes (BHs) often show a characteristic network of dust filaments that terminate in a nuclear spiral in the central parsecs. Here we study the role and fate of these filaments in one of the least accreting BHs known, M31 (10$^{-7}$ L$_{Edd}$) using hydrodynamical simulations. The evolution of a streamer of gas particles moving under the barred potential of M31 is followed from kiloparsec distance to the central parsecs. After an exploratory study of initial conditions, a compelling fit to the observed dust/ionized gas morphologies and line-of-sight velocities in the inner hundreds of parsecs is produced. After several million years of streamer evolution, during which friction, thermal dissipation, and self-collisions have taken place, the gas settles into a disk tens of parsecs wide. This is fed by numerous filaments that arise from an outer circumnuclear ring and spiral toward the center. The final configuration is tightly constrained by a critical input mass in the streamer of several 10$^3$ M$_{\odot}$ (at an injection rate of 10$^{-4}$ M$_{\odot}$ yr$^{-1}$); values above or below this lead to filament fragmentation or dispersion respectively, which are not observed. The creation of a hot gas atmosphere in the region of $\sim$10$^6$ K is key to the development of a nuclear spiral during the simulation. The final inflow rate at 1pc from the center is $\sim$1.7 $\times$ 10$^{-7}$ M$_{\odot}$ yr$^{-1}$, consistent with the quiescent state of the M31 BH., Comment: 13 pages, 11 figures

Gravity is the driving force of star formation. Although gravity is caused by the presence of matter, its role in complex regions is still unsettled. One effective way to study the pattern of gravity is to compute the accretion it exerts on the gas by providing gravitational acceleration maps. A practical way to study acceleration is by computing it using 2D surface density maps, yet whether these maps are accurate remains uncertain. Using numerical simulations, we confirm that the accuracy of the acceleration maps $\mathbf a_{\rm 2D}(x,y)$ computed from 2D surface density are good representations for the mean acceleration weighted by mass. Due to the under-estimations of the distances from projected maps, the magnitudes of accelerations will be over-estimated $|\mathbf a_{\rm 2D}(x,y)| \approx 2.3 \pm 1.8 \; |\mathbf a_{\rm 3D}^{\rm proj}(x,y)|$, where $\mathbf a_{\rm 3D}^{\rm proj}(x,y)$ is mass-weighted projected gravitational acceleration, yet $\mathbf a_{\rm 2D}(x,y)$ and $ \mathbf a_{\rm 3D}^{\rm proj}(x,y)$ stay aligned within 20$^{\circ}$. Significant deviations only occur in regions where multiple structures are present along the line of sight. The acceleration maps estimated from surface density provide good descriptions of the projection of 3D acceleration fields. We expect this technique useful in establishing the link between cloud morphology and star formation, and in understanding the link between gravity and other processes such as the magnetic field. A version of the code for calculating surface density gravitational potential is available at \url{https://github.com/zhenzhen-research/phi_2d}., Comment: Accepted by MNRAS

The onset of star formation is set by the collapse of filaments in the interstellar medium. From a theoretical point of view, an isolated cylindrical filament forms cores via the edge effect. Due to the self-gravity of a filament, the strong increase in acceleration at both ends leads to a pile-up of matter which collapses into cores. However, this effect is rarely observed. Most theoretical models consider a sharp density cut-off at the edge of the filament, whereas a smoother transition is more realistic and would also decrease the acceleration at the ends of the filament. We show that the edge effect can be significantly slowed down by a density gradient, although not completely avoided. However, this allows perturbations inside the filament to grow faster than the edge. We determine the critical density gradient for which the timescales are equal and find it to be of the order of several times the filament radius. Hence, the density gradient at the ends of a filament is an essential parameter for fragmentation and the low rate of observed cases of the edge effect could be naturally explained by shallow gradients., Comment: 10 pages, 9 figures, published in MNRAS

The new experimental measurements of beam spin asymmetry were performed for the deeply virtual exclusive $\pi^0$ production in a wide kinematic region with the photon virtualities $Q^2$ up to 8 GeV$^2$ and the Bjorken scaling variable $x_B$ in the valence regime. The data were collected by the CEBAF Large Acceptance Spectrometer (CLAS12) at Jefferson Lab with longitudinally polarized 10.6 GeV electrons scattered on an unpolarized liquid-hydrogen target. Sizable asymmetry values indicate a substantial contribution from transverse virtual photon amplitudes to the polarized structure functions.The interpretation of these measurements in terms of the Generalized Parton Distributions (GPDs) demonstrates their sensitivity to the chiral-odd GPD $\bar E_T$, which contains information on quark transverse spin densities in unpolarized and polarized nucleons and provides access to the proton's transverse anomalous magnetic moment. Additionally, the data were compared to a theoretical model based on a Regge formalism that was extended to the high photon virtualities., Comment: arXiv admin note: substantial text overlap with arXiv:2210.14557

The electroexcitation amplitudes or $\gamma_vpN^*$ electrocouplings of the $N(1440)1/2^+$, $N(1520)3/2^-$, and $\Delta(1600)3/2^+$ resonances were obtained for the first time from the $ep \to e'\pi^+\pi^-p'$ differential cross sections measured with the CLAS detector at Jefferson Lab within the range of invariant mass $W$ of the final state hadrons from 1.4--1.7 GeV for photon virtualities $Q^2$ from 2.0--5.0 GeV$^2$. The electrocouplings were determined in independent fits of the $\pi^+\pi^-p$ cross sections within three overlapping $W$ intervals with a substantial contribution from each of the three resonances listed above. Consistent results on the electrocouplings extracted from the data in these $W$ intervals provide evidence for their reliable extraction. These studies extend information on the electrocouplings of the $N(1440)1/2^+$ and $N(1520)3/2^-$ available from this channel over a broader range of $Q^2$. The electrocouplings of the $\Delta(1600)3/2^+$, which decays preferentially into $\pi\pi N$ final states, have been determined for the first time. Our results provide further evidence for the structure of these resonances in terms of an interplay between the inner core of three dressed quarks and an external meson-baryon cloud., Comment: 27 pages, 18 figures, 9 tables

This document presents the initial scientific case for upgrading the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab (JLab) to 22 GeV. It is the result of a community effort, incorporating insights from a series of workshops conducted between March 2022 and April 2023. With a track record of over 25 years in delivering the world's most intense and precise multi-GeV electron beams, CEBAF's potential for a higher energy upgrade presents a unique opportunity for an innovative nuclear physics program, which seamlessly integrates a rich historical background with a promising future. The proposed physics program encompass a diverse range of investigations centered around the nonperturbative dynamics inherent in hadron structure and the exploration of strongly interacting systems. It builds upon the exceptional capabilities of CEBAF in high-luminosity operations, the availability of existing or planned Hall equipment, and recent advancements in accelerator technology. The proposed program cover various scientific topics, including Hadron Spectroscopy, Partonic Structure and Spin, Hadronization and Transverse Momentum, Spatial Structure, Mechanical Properties, Form Factors and Emergent Hadron Mass, Hadron-Quark Transition, and Nuclear Dynamics at Extreme Conditions, as well as QCD Confinement and Fundamental Symmetries. Each topic highlights the key measurements achievable at a 22 GeV CEBAF accelerator. Furthermore, this document outlines the significant physics outcomes and unique aspects of these programs that distinguish them from other existing or planned facilities. In summary, this document provides an exciting rationale for the energy upgrade of CEBAF to 22 GeV, outlining the transformative scientific potential that lies within reach, and the remarkable opportunities it offers for advancing our understanding of hadron physics and related fundamental phenomena., Comment: Updates to the list of authors; Preprint number changed from theory to experiment; Updates to sections 4 and 6, including additional figures

The orbital distribution of the S-star cluster surrounding the supermassive black hole in the center of the Milky Way is analyzed. A tight, roughly exponential dependence of the pericenter distance r$_{p}$ on orbital eccentricity e$_{\star}$ is found, $\log ($r$_p)\sim$(1-e$_{\star}$), which cannot be explained simply by a random distribution of semi-major axes and eccentricities. No stars are found in the region with high e$_{\star}$ and large log r$_{p}$ or in the region with low e$_{\star}$ and small log r$_{p}$. G-clouds follow the same correlation. The likelihood P(log r$_p$,(1-e$_{\star}$)) to determine the orbital parameters of S-stars is determined. P is very small for stars with large e$_{\star}$ and large log r$_{p}$. S-stars might exist in this region. To determine their orbital parameters, one however needs observations over a longer time period. On the other hand, if stars would exist in the region of low log r$_{p}$ and small e$_{\star}$, their orbital parameters should by now have been determined. That this region is unpopulated therefore indicates that no S-stars exist with these orbital characteristics, providing constraints for their formation. We call this region, defined by $\log$ (r$_p$/AU) $<$ 1.57+2.6(1-e$_{\star})$, the zone of avoidance. Finally, it is shown that the observed frequency of eccentricities and pericenter distances is consistent with a random sampling of log r$_{p}$ and e$_{\star}$. However, only if one takes into account that no stars exist in the zone of avoidance and that orbital parameters cannot yet be determined for stars with large r$_{p}$ and large e$_{\star}$., Comment: 15 pages, 4 figures, submitted to ApJ, comments very welcome

We report high-quality H${\alpha}$/CO, imaging spectroscopy of nine massive (log median stellar mass = 10.65 $M_{\odot}$), disk galaxies on the star-forming, main sequence (henceforth `SFGs'), near the peak of cosmic galaxy evolution ($z\sim$1.1-2.5), taken with the ESO-Very Large Telescope, IRAM-NOEMA and Atacama Large Millimeter/submillimeter Array. We fit the major axis position-velocity cuts with beam-convolved, forward models with a bulge, a turbulent rotating disk, and a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from our previous rotation curve analyses of these galaxies. We then subtract the inferred 2D model-galaxy velocity and velocity dispersion maps from those of the observed galaxies. We investigate whether the residual velocity and velocity dispersion maps show indications for radial flows. We also carry out kinemetry, a model-independent tool for detecting radial flows. We find that all nine galaxies exhibit significant non-tangential flows. In six SFGs, the inflow velocities ($v_r\sim$30-90 km s$^{-1}$, 10%-30% of the rotational component) are along the minor axis of these galaxies. In two cases the inflow appears to be off the minor axis. The magnitudes of the radial motions are in broad agreement with the expectations from analytic models of gravitationally unstable, gas-rich disks. Gravitational torques due to clump and bar formation, or spiral arms, drive gas rapidly inward and result in the formation of central disks and large bulges. If this interpretation is correct, our observations imply that gas is transported into the central regions on ~10 dynamical time scales., Comment: 31 pages, 27 figures, 3 tables. The Astrophysical Journal in press

The polarized cross section ratio $\sigma_{LT'}/\sigma_{0}$ from hard exclusive $\pi^{-} \Delta^{++}$ electroproduction off an unpolarized hydrogen target has been extracted based on beam-spin asymmetry measurements using a 10.2 GeV / 10.6 GeV incident electron beam and the CLAS12 spectrometer at Jefferson Lab. The study, which provides the first observation of this channel in the deep-inelastic regime, focuses on very forward-pion kinematics in the valence regime, and photon virtualities ranging from 1.5 GeV$^{2}$ up to 7 GeV$^{2}$. The reaction provides a novel access to the $d$-quark content of the nucleon and to $p \rightarrow \Delta^{++}$ transition generalized parton distributions. A comparison to existing results for hard exclusive $\pi^{+} n$ and $\pi^{0} p$ electroproduction is provided, which shows a clear impact of the excitation mechanism, encoded in transition generalized parton distributions, on the asymmetry.

The physics of the gravitational form factors of the proton, and their understanding within quantum chromodynamics, has advanced significantly in the past two decades through both theory and experiment. This Colloquium provides an overview of this progress, highlights the physical insights unveiled by studies of gravitational form factors, and reviews their interpretation in terms of the mechanical properties of the proton., Comment: 22 pages, 15 figures

This White Paper presents the community inputs and scientific conclusions from the Hot and Cold QCD Town Meeting that took place September 23-25, 2022 at MIT, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 424 physicists registered for the meeting. The meeting highlighted progress in Quantum Chromodynamics (QCD) nuclear physics since the 2015 LRP (LRP15) and identified key questions and plausible paths to obtaining answers to those questions, defining priorities for our research over the coming decade. In defining the priority of outstanding physics opportunities for the future, both prospects for the short (~ 5 years) and longer term (5-10 years and beyond) are identified together with the facilities, personnel and other resources needed to maximize the discovery potential and maintain United States leadership in QCD physics worldwide. This White Paper is organized as follows: In the Executive Summary, we detail the Recommendations and Initiatives that were presented and discussed at the Town Meeting, and their supporting rationales. Section 2 highlights major progress and accomplishments of the past seven years. It is followed, in Section 3, by an overview of the physics opportunities for the immediate future, and in relation with the next QCD frontier: the EIC. Section 4 provides an overview of the physics motivations and goals associated with the EIC. Section 5 is devoted to the workforce development and support of diversity, equity and inclusion. This is followed by a dedicated section on computing in Section 6. Section 7 describes the national need for nuclear data science and the relevance to QCD research., Comment: QCD Town Meeting White Paper, as submitted to 2023 NSAC LRP committee on Feb. 28, 2023

We study the relation between the metallicities of ionised and atomic gas in star-forming galaxies at z=0-3 using the EAGLE cosmological, hydrodynamical simulations. This is done by constructing a dense grid of sightlines through the simulated galaxies and obtaining the star formation rate- and HI column density-weighted metallicities, Z_{SFR} and Z_{HI}, for each sightline as proxies for the metallicities of ionised and atomic gas, respectively. We find Z_{SFR} > Z_{HI} for almost all sightlines, with their difference generally increasing with decreasing metallicity. The stellar masses of galaxies do not have a significant effect on this trend, but the positions of the sightlines with respect to the galaxy centres play an important role: the difference between the two metallicities decreases when moving towards the galaxy centres, and saturates to a minimum value in the central regions of galaxies, irrespective of redshift and stellar mass. This implies that the mixing of the two gas phases is most efficient in the central regions of galaxies where sightlines generally have high column densities of HI. However, a high HI column density alone doesn't guarantee a small difference between the two metallicities. In galaxy outskirts, the inefficiency of the mixing of star-forming gas with HI seems to dominate over the dilution of heavy elements in HI through mixing with the pristine gas. We find good agreement between the available observational data and the Z_{SFR}-Z_{HI} relation predicted by the EAGLE simulations. Though, observed regions with a nuclear-starburst mode of star formation appear not to follow the same relation., Comment: Accepted for publication in ApJ

This paper presents a comprehensive review of both the theory and experimental successes of Quantum Chromodynamics, starting with its emergence as a well defined theory in 1972-73 and following developments and results up to the present day. Topics include a review of the earliest theoretical and experimental foundations; the fundamental constants of QCD; an introductory discussion of lattice QCD, the only known method for obtaining exact predictions from QCD; methods for approximating QCD, with special focus on effective field theories; QCD under extreme conditions; measurements and predictions of meson and baryon states; a special discussion of the structure of the nucleon; techniques for study of QCD at high energy, including treatment of jets and showers; measurements at colliders; weak decays and quark mixing; and a section on the future, which discusses new experimental facilities or upgrades currently funded. The paper is intended to provide a broad background for Ph.D. students and postdocs starting their career. Some contributions include personal accounts of how the ideas or experiments were developed., Comment: Invited volume for the EJPC; 567 pages if text/figures and 4783 references occupying about 160 additional pages. arXiv abstract abridged, for the complete abstract please see the full text

The Heavy Photon Search experiment (HPS) at the Thomas Jefferson National Accelerator Facility searches for electro-produced dark photons. We report results from the 2016 Engineering Run consisting of 10608/nb of data for both the prompt and displaced vertex searches. A search for a prompt resonance in the $e^+e^-$ invariant mass distribution between 39 and 179 MeV showed no evidence of dark photons above the large QED background, limiting the coupling of {\epsilon}^2 {\geq} 10^-5, in agreement with previous searches. The search for displaced vertices showed no evidence of excess signal over background in the masses between 60 and 150 MeV, but had insufficient luminosity to limit canonical heavy photon production. This is the first displaced vertex search result published by HPS. HPS has taken high-luminosity data runs in 2019 and 2021 that will explore new dark photon phase space., Comment: 28 pages, 46 figures

This is a contribution to the review 50 Years of Quantum Chromodynamics edited by F. Gross and E. Klempt, to be published in Journal EPJC. This contribution reviews the nucleon resonance transition form factors determined from meson electro-production experiments at electron accelerator facilities, i.e. this contribution focuses on space-like transition form factors and amplitudes. Comparisons are made when available to LQCD and to approaches with traceable links to strong QCD and to advanced quark model calculations, Comment: 14 pages, 12 figures