Search

Recent observations of volume-limited samples of magnetic white dwarfs (WD) have revealed a higher incidence of magnetism in older WDs. Specifically, these studies indicate that magnetism is more prevalent in WDs with fully or partially crystallized cores compared to those with entirely liquid cores. This has led to the recognition of a crystallization-driven dynamo as an important mechanism for explaining magnetism in isolated WDs. However, recent simulations challenged the capability of this mechanism to match both the incidence of magnetism and the field strengths detected in WDs. In this letter, we explore an alternative hypothesis for the surface emergence of magnetic fields in isolated WDs. WDs with masses $\gtrsim 0.55 M_\odot$ are the descendants of main-sequence stars with convective cores capable of generating strong dynamo magnetic fields. This idea is supported by asteroseismic evidence of strong magnetic fields buried within the interiors of red giant branch stars. Assuming that these fields are disrupted by subsequent convective zones, we have estimated magnetic breakout times for WDs. Due to the significant uncertainties in breakout times stemming from the treatment of convective boundaries and mass loss rates, we cannot provide a precise prediction for the emergence time of the main-sequence dynamo field. However, we can predict that this emergence should occur during the WD phase for WDs with masses $\gtrsim 0.65 M_\odot$. We also find that the magnetic breakout is expected to occur earlier in more massive WDs, consistently with observations from volume-limited samples and the well-established fact that magnetic WDs tend to be more massive than non-magnetic ones. Moreover, within the uncertainties of stellar evolutionary models, we find that the emergence of main-sequence dynamo magnetic fields can account for a significant portion of the magnetic WDs., Comment: accepted for publication in A&A letters

We present the first measurement of the $^{47}$K($d,p\gamma$)$^{48}$K transfer reaction, performed in inverse kinematics using a reaccelerated beam of $^{47}$K. The level scheme of $^{48}$K has been greatly extended with nine new bound excited states identified and spectroscopic factors deduced. Detailed comparisons with SDPF-U and SDPF-MU shell-model calculations reveal a number of discrepancies with these results, and a preference for SDPF-MU is found. Intriguingly, an apparent systematic overestimation of spectroscopic factors and a poor reproduction of the energies for 1$^-$ states suggests that the mixing between the $\pi s^{\,\,\,1}_{1/2} d^{\,\,\,4}_{3/2}$ and $\pi s^{\,\,\,2}_{1/2} d^{\,\,\,3}_{3/2}$ proton configurations in $^{48}$K is not correctly described using current interactions, challenging our descriptions of light $N=28$ nuclei., Comment: 7 pages, 3 figures, 1 table

The Gemini South telescope is now equipped with a new high-resolution spectrograph called GHOST (the Gemini High-resolution Optical SpecTrograph). This instrument provides high-efficiency, high-resolution spectra covering 347-1060 nm in a single exposure of either one or two targets simultaneously, along with precision radial velocity spectroscopy utilizing an internal calibration source. It can operate at a spectral element resolving power of either 76000 or 56000, and can reach a SNR$\sim$5 in a 1hr exposure on a V$\sim$20.8 mag target in median site seeing, and dark skies (per resolution element). GHOST was installed on-site in June 2022, and we report performance after full integration to queue operations in November 2023, in addition to scientific results enabled by the integration observing runs. These results demonstrate the ability to observe a wide variety of bright and faint targets with high efficiency and precision. With GHOST, new avenues to explore high-resolution spectroscopy have opened up to the astronomical community. These are described, along with the planned and potential upgrades to the instrument., Comment: Accepted in the Astronomical Journal; 26 pages, 24 figures and 4 tables

Recent observational constraints on the internal structure of Jupiter and Saturn suggest that these planets have ``fuzzy" cores, i.e., gradients of the concentration of heavy elements that might span a large fraction of the planet's radius. These cores could be composed of a semi-convective staircase, i.e., multiple convective layers separated by diffusive interfaces arising from double-diffusive instabilities. However, to date, no study has demonstrated how such staircases can avoid layer mergers and persist over evolutionary time scales. In fact, previous work has found that these mergers occur rapidly, leading to only a single convective layer. Using 3D simulations, we demonstrate that rotation prolongs the lifetime of a convective staircase by increasing the timescale for both layer merger and erosion of the interface between the final two layers. We present an analytic model for the erosion phase, predicting that rotation increases the erosion time by a factor of approximately $\mathrm{Ro}^{-1/2}$, where $\mathrm{Ro}$ is the Rossby number of the convective flows (the ratio of the rotation period to the convective turnover time). For Jovian conditions at early times after formation (when convection is vigorous enough to mix a large fraction of the planet), we find the erosion time to be roughly $10^{9}~\mathrm{yrs}$ in the non-rotating case and $10^{11}~\mathrm{yrs}$ in the rotating case. If these timescales are confirmed with a larger suite of numerical simulations, the existence of convective staircases within the deep interiors of giant planets is a strong possibility, and rotation could be an important factor in the preservation of their fuzzy cores., Comment: Accepted for publication in ApJ Letters

The dynamics of a neutron star after a glitch involve the transfer of angular momentum from the crust (where the glitch is presumed to originate) to the liquid core, causing the core to spin up. The crust-core coupling, which determines how quickly this spin-up proceeds, can be achieved through various physical processes, including Ekman pumping, superfluid vortex-mediated mutual friction, and magnetic fields. Although the complex nature of these mechanisms has made it difficult to study their combined effects, analytical estimations for individual processes reveal that spin-up timescales vary according to the relative strength of Coriolis, viscous, and mutual friction forces, as well as the magnetic field. However, experimental and numerical validations of those analytical predictions are limited. In this paper, we focus on viscous effects and mutual friction. We conduct non-linear hydrodynamical simulations of the spin-up problem in a two-component fluid by solving the incompressible Hall-Vinen-Bekarevich-Khalatnikov (HVBK) equations in the full sphere (i.e., including $r=0$) for the first time. We find that the viscous (normal) component accelerates due to Ekman pumping, although the mutual friction coupling to the superfluid component alters the spin-up dynamics compared to the single-fluid scenario. Close to the sphere's surface, the response of the superfluid is accurately described by the mutual friction timescale irrespective of its coupling strength with the normal component. However, as we move deeper into the sphere, the superfluid accelerates on different timescales due to the slow viscous spin-up of the internal normal fluid layers. We discuss potential implications for neutron stars and requirements for future work to build more realistic models., Comment: Accepted for Publication in ApJ

The excited states of unstable $^{20}$O were investigated via $\gamma$-ray spectroscopy following the $^{19}$O$(d,p)^{20}$O reaction at 8 $A$MeV. By exploiting the Doppler Shift Attenuation Method, the lifetime of the 2$^+_2$ and 3$^+_1$ states were firmly established. From the $\gamma$-ray branching and E2/M1 mixing ratios for transitions deexciting the 2$^+_2$ and 3$^+_1$ states, the B(E2) and B(M1) were determined. Various chiral effective field theory Hamiltonians, describing the nuclear properties beyond ground states, along with a standard USDB interaction, were compared with the experimentally obtained data. Such a comparison for a large set of $\gamma$-ray transition probabilities with the valence space in medium similarity renormalization group ab-initio calculations was performed for the first time in a nucleus far from stability. It was shown that the ab-initio approaches using chiral EFT forces are challenged by detailed high-precision spectroscopic properties of nuclei. The reduced transition probabilities were found to be a very constraining test of the performance of the ab-initio models., Comment: Supplemental Material available

The origin of large magnetic fields ($\gtrsim 10^6~\mathrm{G}$) in isolated white dwarfs is not clear. One possible explanation is that crystallization of the star's core drives compositional convection, which when combined with the star's rotation, can drive a dynamo. However, whether convection is efficient enough to explain the large intensity of the observed magnetic fields is still under debate. Recent work has shown that convection in cooling white dwarfs spans two regimes: efficient convection at the onset of crystallization, and thermohaline convection during most of the star's cooling history. Here, we calculate the properties of crystallization-driven convection for cooling models of several white dwarfs of different masses. We combine mixing-length theory with scalings from magneto-rotational convection to estimate the typical magnitude of the convective velocity and induced magnetic field for both scenarios. In the thermohaline regime, we find velocities $\sim 10^{-6}$--$10^{-5}~\mathrm{cm~s^{-1}}$, with fields restricted to $\lesssim~100~\mathrm{G}$. However, when convection is efficient, the flow velocity can reach magnitudes of $\sim 10^2$--$10^3~\mathrm{cm~s^{-1}}$, with fields of $\sim 10^6$--$10^8~\mathrm{G}$, independent of the star's rotation rate. Thus, dynamos driven at the onset of crystallization could explain the large intensity magnetic fields measured for single white dwarfs., Comment: Accepted for publication in ApJ Letters

We investigate crystallization-driven convection in carbon-oxygen white dwarfs. We present a version of the mixing length theory (MLT) that self-consistently includes the effects of thermal diffusion and composition gradients, and provides solutions for the convective parameters based on the local heat and composition fluxes. Our formulation smoothly transitions between the regimes of fast adiabatic convection at large Peclet number and slow thermohaline convection at low Peclet number. It also allows for both thermally-driven and compositionally-driven convection, including correctly accounting for the direction of heat transport for compositionally-driven convection in a thermally-stable background. We use the MESA stellar evolution code to calculate the composition and heat fluxes during crystallization in different models of cooling white dwarfs, and determine the regime of convection and the convective velocity. We find that convection occurs in the regime of slow thermohaline convection during most of the cooling history of the star. However, at the onset of crystallization, the composition flux is large enough to drive fast overturning convection for a short time ($\sim 10\ \mathrm{Myr}$). We estimate the convective velocities in both of these phases and discuss the implications for explaining observed white dwarf magnetic fields with crystallization-driven dynamos., Comment: ApJ paper

Recent measurements of Jupiter's gravitational field (by Juno) and seismology of Saturn's rings (by Cassini) strongly suggest that both planets have a stably-stratified core that still possesses a primordial gradient in the concentration of heavy elements. The existence of such a "diffusely" stratified core has been a surprise as it was long expected that the Jovian planets should be fully convective and hence fully mixed. A vigorous zone of convection, driven by surface cooling, forms at the surface and deepens through entrainment of fluid from underneath. In fact, it was believed that this convection zone should grow so rapidly that the entire planet would be consumed in less than a million years. Here we suggest that two processes, acting in concert, present a solution to this puzzle. All of the giant planets are rapidly rotating and have a cooling rate that declines with time. Both of these effects reduce the rate of fluid entrainment into the convection zone. Through the use of an analytic prescription of entrainment in giant planets, we demonstrate that these two effects, rotation and dwindling surface cooling, result in a convection zone which initially grows but eventually stalls. The depth to which the convective interface asymptotes depends on the rotation rate and on the stratification of the stable interior. Conversely, in a nonrotating planet, or in a planet that maintains a higher level of cooling than current models suggest, the convection zone deepens forever, eventually spanning the entire planet., Comment: 7 pages, 2 figures, accepted for publication by Astrophysical Journal Letters

Meridional circulation regulates the Sun's interior dynamics and magnetism. While it is well accepted that meridional flows are poleward at the Sun's surface, helioseismic observations have yet to provide a definitive answer for the depth at which those flows return to the equator, or the number of circulation cells in depth. Here, we explore the observability of multiple circulation cells stacked in radius. Specifically, we examine the seismic signature of several meridional flow profiles by convolving time-distance averaging kernels with mean flows obtained from a suite of 3D hydrodynamic simulations. At mid and high latitudes, we find that weak flow structures in the deep convection zone can be obscured by signals from the much stronger surface flows. This contamination of 1--2 m s$^{-1}$ is caused by extended side lobes in the averaging kernels, which produce a spurious equatorward signal with flow speeds that are one order of magnitude stronger than the original flow speeds in the simulations. At low latitudes, the flows in the deep layers of the simulations are stronger ($> 2$ m s$^{-1}$) and multiple cells across the convection zone can produce a sufficiently strong signal to survive the convolution process. Now that meridional flows can be measured over two decades of data, the uncertainties arising from convective noise have fallen to a level where they are comparable in magnitude to the systematic biases caused by non-local features in the averaging kernels. Hence, these systematic errors are beginning to influence current helioseismic deductions and need broader consideration., Comment: New version including suggestions by the reviewers

Recent measurements of Jupiter's gravitational moments by the Juno spacecraft and seismology of Saturn's rings suggest that the primordial composition gradients in the deep interior of these planets have persisted since their formation. One possible explanation is the presence of a double-diffusive staircase below the planet's outer convection zone, which inhibits mixing across the deeper layers. However, hydrodynamic simulations have shown that these staircases are not long-lasting and can be disrupted by overshooting convection. In this paper we suggests that planetary rotation could be another factor for the longevity of primordial composition gradients. Using rotational mixing-length theory and 3D hydrodynamic simulations, we demonstrate that rotation significantly reduces both the convective velocity and the mixing of primordial composition gradients. In particular, for Jovian conditions at $t\sim10^{8}~\mathrm{yrs}$ after formation, rotation reduces the convective velocity by a factor of 6, and in turn, the kinetic energy flux available for mixing gets reduced by a factor of $6^3\sim 200$. This leads to an entrainment timescale that is more than two orders of magnitude longer than without rotation. We encourage future hydrodynamic models of Jupiter and other gas giants to include rapid rotation, because the decrease in the mixing efficiency could explain why Jupiter and Saturn are not fully mixed., Comment: Accepted for publication in the Astrophysical Journal Letters

The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements., Comment: White paper submitted to the Decadal Survey for Solar and Space Physics 2024-2033. 7 pages; 4 figures

In the forensic studies of painting masterpieces, the analysis of the support is of major importance. For plain weave fabrics, the densities of vertical and horizontal threads are used as main features, while angle deviations from the vertical and horizontal axis are also of help. These features can be studied locally through the canvas. In this work, deep learning is proposed as a tool to perform these local densities and angle studies. We trained the model with samples from 36 paintings by Vel\'azquez, Rubens or Ribera, among others. The data preparation and augmentation are dealt with at a first stage of the pipeline. We then focus on the supervised segmentation of crossing points between threads. The U-Net with inception and Dice loss are presented as good choices for this task. Densities and angles are then estimated based on the segmented crossing points. We report test results of the analysis of a few canvases and a comparison with methods in the frequency domain, widely used in this problem. We concluded that this new approach succeeds in some cases where the frequency analysis tools fail, while improving the results in others. Besides, our proposal does not need the labeling of part of the to-be-processed image. As case studies, we apply this novel algorithm to the analysis of two pairs of canvases by Vel\'azquez and Murillo, to conclude that the fabrics used came from the same roll., Comment: 21 pages, 20 figures, 1 table, manuscript submitted to journal

The low-lying structure of $^{15}$C has been investigated via the neutron-removal $^{16}$C$(d,t)$ reaction. Along with bound neutron sd-shell hole states, unbound p-shell hole states have been firmly confirmed. The excitation energies and the deduced spectroscopic factors of the cross-shell states are an important measure of the $[(p)^{-1}(sd)^{2}]$ neutron configurations in $^{15}$C. Our results show a very good agreement with shell-model calculations using the SFO-tls interaction for $^{15}$C. However, a modification of the $p$-$sd$ and $sd$-$sd$ monopole terms was applied in order to reproduce the $N=9$ isotone $^{17}$O. In addition, the excitation energies and spectroscopic factors have been compared to the first calculations of $^{15}$C with the $ab~ initio$ self-consistent Green's function method employing the NNLO$_{sat}$ interaction. The results show the sensitivity to the size of the $N=8$ shell gap and highlight the need of going beyond the current truncation scheme in the theory.

We investigate heat transport associated with compositionally-driven convection driven by crystallization at the ocean-crust interface in accreting neutron stars, or growth of the solid core in cooling white dwarfs. We study the effect of thermal diffusion and rapid rotation on the convective heat transport, using both mixing length theory and numerical simulations of Boussinesq convection. We determine the heat flux, composition gradient and P\'eclet number, $\mathrm{Pe}$ (the ratio of thermal diffusion time to convective turnover time) as a function of the composition flux. We find two regimes of convection with a rapid transition between them as the composition flux increases. At small Pe, the ratio between the heat flux and composition flux is independent of Pe,, because the loss of heat from convecting fluid elements due to thermal diffusion is offset by the smaller composition gradient needed to overcome the reduced thermal buoyancy. At large Pe, the temperature gradient approaches the adiabatic gradient, saturating the heat flux. We discuss the implications for neutron star and white dwarf cooling. Convection in neutron stars spans both regimes. We find rapid mixing of neutron star oceans, with a convective turnover time of order weeks to minutes depending on rotation. Except during the early stages of core crystallization, white dwarf convection is in the thermal-diffusion-dominated fingering regime. We find convective velocities much smaller than recent estimates for crystallization-driven dynamos. The small fraction of energy carried as kinetic energy calls into question the effectiveness of crystallization-driven dynamos as an explanation for observed white dwarf magnetic fields., Comment: Accepted for Publication in ApJ (Astrophysical Journal)

To support the expected increase in aquaculture production during the next years, a wider range of alternative ingredients to fishmeal is needed, towards contributing to an increase in production sustainability. This study aimed to test diets formulated with non-conventional feed ingredients on gilthead seabream (Sparus aurata) growth performance, feed utilization, apparent digestibility of nutrients and nutrient outputs to the environment. Four isonitrogenous and isoenergetic diets were formulated: a control diet (CTRL) similar to a commercial feed and three experimental diets containing, as main protein sources, plant by-products, glutens and concentrates (PLANT); processed animal proteins (PAP); or micro/macroalgae, insect meals and yeast (EMERG). Diets were tested in triplicate during 80 days. The (EMERG) treatment resulted in lower fish growth performance, higher FCR and lower nutrient and energy retentions than the other treatments. The lowest protein digestibility was found for the EMERG diet, which caused increased nitrogen losses. The PLANT and PAP treatments resulted in better fish growth performance, higher nutrient and energy retentions, and lower FCR than the CTRL treatment. The significant improvement in FCR found for fish fed PLANT and PAP diets and the high protein digestibility of these diets contribute towards minimizing the environmental impacts of seabream production

Radio images of protoplanetary disks demonstrate that dust grains tend to organize themselves into rings. These rings may be a consequence of dust trapping within gas pressure maxima wherein the local high dust-to-gas ratio is expected to trigger the formation of planetesimals and eventually planets. We revisit the behavior of dust near gas pressure perturbations enforced by a planet in two-dimensional, shearing box simulations. While dust grains collect into generally long-lived rings, particles with small Stokes parameter $\tau_s < 0.1$ tend to advect out of the ring within a few drift timescales. Scaled to the properties of ALMA disks, we find that rings composed of larger particles ($\tau_s \geq 0.1$) can nucleate a dust clump massive enough to trigger pebble accretion which proceeds to ingest the entire dust ring well within $\sim$1 Myr. To ensure the survival of the dust rings, we favor a non-planetary origin and typical grain size $\tau_s \lesssim 0.05$--0.1. Planet-driven rings may still be possible but if so we would expect the orbital distance of the dust rings to be larger for older systems., Comment: To appear in ApJ. Paper refocused to forcing by a planet

In 1D evolution models of gas giant planets, an outer convection zone advances into the interior as the surface cools, and multiple convective layers form beneath that convective front. To study layer formation below an outer convection zone in a similar scenario, we investigate the evolution of a stably-stratified fluid with a linear composition gradient that is constantly being cooled from above. We use the Boussinesq approximation in a series of 2D simulations at low and high Prandtl numbers ($\mathrm{Pr} = 0.5$ and 7), initialized with constant temperature everywhere, and cooled at different rates. We find that multiple convective layers form at $\mathrm{Pr} = 7$, {as the result of an instability in the} diffusive thermal boundary layer below the outer convection zone. At low Pr, layers do not form because the temperature gradient within the boundary layer is much smaller than at large Pr and, consequently, is not large enough to overcome the stabilizing effect of the composition gradient. For the stratification used in this study, on the long-term the composition gradient is an ineffective barrier against the propagation of the outer convection zone and the entire fluid becomes fully-mixed, whether layers form or not. Our results challenge 1D evolutionary models of gas giant planets, which predict that layers are long-lived and that the outer convective envelope stops advancing inwards. We discuss what is needed for future work to build more realistic models., Comment: Accepted for publication in Physical review fluids

Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions which are unstable according to the Ledoux criterion are convective. Initially, the radiative zone adjacent to the convection zone is Schwarzschild-unstable but Ledoux-stable due to a composition gradient. Over many convective overturn timescales the convection zone grows via entrainment. The convection zone saturates at the size originally predicted by the Schwarzschild criterion, although in this final state the Schwarzschild and Ledoux criteria agree. Therefore, the Schwarzschild criterion should be used to determine the size of stellar convection zones, except possibly during short-lived evolutionary stages in which entrainment persists., Comment: Manuscript accepted in ApJ Letters. Supplemental materials are in a Zenodo repository (https://zenodo.org/record/6344868). Code is in a Github repository (https://github.com/evanhanders/schwarzschild_or_ledoux)

We investigate the onset and evolution of zonal flows in a growing convective layer when a stably-stratified fluid with a composition gradient is cooled from above. This configuration allows the study of zonal flows for a wide range of values of the Rayleigh number, $Ra$, and aspect ratio of the convection zone within a given simulation. We perform a series of 2D simulations using the Boussinesq approximation, with aspect ratio of the computational domain between $1$ and $5$, and Prandtl number $Pr = 0.1$, 0.5, 1, and $7$. We find that for square domains zonal flows appear when the aspect ratio of the convective layer is smaller than two, and the evolution of the system depends on the Prandtl number. For $Pr\leq 1$, the fluid experiences bursts of convective transport with negligible convective transport between bursts. The magnitude and frequency of the bursts are smaller at low $Pr$, which suggests that the bursting regime is stronger in a narrow range around $Pr=1$, as observed in previous studies of thermal convection. For $Pr=7$, the structure of the flow consists of tilted convective plumes, and the convective transport is sustained at all times. In wider domains, the aspect ratio of the convective zone is always much larger than two and zonal flows do not appear. These results confirm and extend to fluids with stable composition gradients previous findings on thermal convection. The fact that zonal flows can be avoided by using computational domains with large aspect ratios opens up the possibility of 2D studies of convective overshoot, layer formation and transport properties across diffusive interfaces., Comment: Accepted for publication in Physical Review Fluids

We aim to demonstrate the scientific potential of the Gaia Early Data Release 3 (EDR3) for the study of the Milky Way structure and evolution. We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. We explore the disturbances of the current disc, the spatial and kinematical distributions of early accreted versus in-situ stars, the structures in the outer parts of the disc, and the orbits of open clusters Berkeley 29 and Saurer 1. We find that: i) the dynamics of the Galactic disc are very complex with vertical asymmetries, and new correlations, including a bimodality with disc stars with large angular momentum moving vertically upwards from below the plane, and disc stars with slightly lower angular momentum moving preferentially downwards; ii) we resolve the kinematic substructure (diagonal ridges) in the outer parts of the disc for the first time; iii) the red sequence that has been associated with the proto-Galactic disc that was present at the time of the merger with Gaia-Enceladus-Sausage is currently radially concentrated up to around 14 kpc, while the blue sequence that has been associated with debris of the satellite extends beyond that; iv) there are density structures in the outer disc, both above and below the plane, most probably related to Monoceros, the Anticentre Stream, and TriAnd, for which the Gaia data allow an exhaustive selection of candidate member stars and dynamical study; and v) the open clusters Berkeley~29 and Saurer~1, despite being located at large distances from the Galactic centre, are on nearly circular disc-like orbits. We demonstrate how, once again, the Gaia are crucial for our understanding of the different pieces of our Galaxy and their connection to its global structure and history., Comment: Gaia EDR3 performance verification paper, version 2 closer to published version in A&A, complete list of authors

We produce a clean and well-characterised catalogue of objects within 100\,pc of the Sun from the \G\ Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use. The selection of objects within 100\,pc from the full catalogue used selected training sets, machine-learning procedures, astrometric quantities, and solution quality indicators to determine a probability that the astrometric solution is reliable. The training set construction exploited the astrometric data, quality flags, and external photometry. For all candidates we calculated distance posterior probability densities using Bayesian procedures and mock catalogues to define priors. Any object with reliable astrometry and a non-zero probability of being within 100\,pc is included in the catalogue. We have produced a catalogue of \NFINAL\ objects that we estimate contains at least 92\% of stars of stellar type M9 within 100\,pc of the Sun. We estimate that 9\% of the stars in this catalogue probably lie outside 100\,pc, but when the distance probability function is used, a correct treatment of this contamination is possible. We produced luminosity functions with a high signal-to-noise ratio for the main-sequence stars, giants, and white dwarfs. We examined in detail the Hyades cluster, the white dwarf population, and wide-binary systems and produced candidate lists for all three samples. We detected local manifestations of several streams, superclusters, and halo objects, in which we identified 12 members of \G\ Enceladus. We present the first direct parallaxes of five objects in multiple systems within 10\,pc of the Sun., Comment: 45 Pages, 39 figures in main part and 18 in appendix, tables on CDS

Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions. Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar system barycentre with respect to the rest frame of the Universe. Apart from being an important scientific result by itself, the acceleration measured in this way is a good quality indicator of the Gaia astrometric solution. Methods. The effect of the acceleration is obtained as a part of the general expansion of the vector field of proper motions in Vector Spherical Harmonics (VSH). Various versions of the VSH fit and various subsets of the sources are tried and compared to get the most consistent result and a realistic estimate of its uncertainty. Additional tests with the Gaia astrometric solution are used to get a better idea on possible systematic errors in the estimate. Results. Our best estimate of the acceleration based on Gaia EDR3 is $(2.32 \pm 0.16) \times 10^{-10}$ m s${}^{-2}$ (or $7.33 \pm 0.51$ km s$^{-1}$ Myr${}^{-1}$) towards $\alpha = 269.1^\circ \pm 5.4^\circ$, $\delta = -31.6^\circ \pm 4.1^\circ$, corresponding to a proper motion amplitude of $5.05 \pm 0.35$ $\mu$as yr${}^{-1}$. This is in good agreement with the acceleration expected from current models of the Galactic gravitational potential. We expect that future Gaia data releases will provide estimates of the acceleration with uncertainties substantially below 0.1 $\mu$as yr${}^{-1}$., Comment: A&A, accepted

We compare the Gaia DR2 and Gaia EDR3 performances in the study of the Magellanic Clouds and show the clear improvements in precision and accuracy in the new release. We also show that the systematics still present in the data make the determination of the 3D geometry of the LMC a difficult endeavour; this is at the very limit of the usefulness of the Gaia EDR3 astrometry, but it may become feasible with the use of additional external data. We derive radial and tangential velocity maps and global profiles for the LMC for the several subsamples we defined. To our knowledge, this is the first time that the two planar components of the ordered and random motions are derived for multiple stellar evolutionary phases in a galactic disc outside the Milky Way, showing the differences between younger and older phases. We also analyse the spatial structure and motions in the central region, the bar, and the disc, providing new insights into features and kinematics. Finally, we show that the Gaia EDR3 data allows clearly resolving the Magellanic Bridge, and we trace the density and velocity flow of the stars from the SMC towards the LMC not only globally, but also separately for young and evolved populations. This allows us to confirm an evolved population in the Bridge that is slightly shift from the younger population. Additionally, we were able to study the outskirts of both Magellanic Clouds, in which we detected some well-known features and indications of new ones., Comment: This paper is part of the "demonstration papers" released with Gaia EDR3: https://www.cosmos.esa.int/web/gaia/earlydr3

We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motions, and the (G_BP-G_RP) colour are also available. The passbands for G, G_BP, and G_RP are provided as part of the release. For ease of use, the 7 million radial velocities from Gaia DR2 are included in this release, after the removal of a small number of spurious values. New radial velocities will appear as part of Gaia DR3. Finally, Gaia EDR3 represents an updated materialisation of the celestial reference frame (CRF) in the optical, the Gaia-CRF3, which is based solely on extragalactic sources. The creation of the source list for Gaia EDR3 includes enhancements that make it more robust with respect to high proper motion stars, and the disturbing effects of spurious and partially resolved sources. The source list is largely the same as that for Gaia DR2, but it does feature new sources and there are some notable changes. The source list will not change for Gaia DR3. Gaia EDR3 represents a significant advance over Gaia DR2, with parallax precisions increased by 30 percent, proper motion precisions increased by a factor of 2, and the systematic errors in the astrometry suppressed by 30--40 percent for the parallaxes and by a factor ~2.5 for the proper motions. The photometry also features increased precision, but above all much better homogeneity across colour, magnitude, and celestial position. A single passband for G, G_BP, and G_RP is valid over the entire magnitude and colour range, with no systematics above the 1 percent level., Comment: Accepted for A&A Special Issue on Gaia EDR3, 21 pages, 2 figures. This version includes the updates in the erratum (https://doi.org/10.1051/0004-6361/202039657e)

The first investigation of the single-particle structure of the bound states of 17C, via the d(16C, p) transfer reaction, has been undertaken. The measured angular distributions confirm the spin-parity assignments of 1/2+ and 5/2+ for the excited states located at 217 and 335 keV, respectively. The spectroscopic factors deduced for these states exhibit a marked single-particle character, in agreement with shell model and particle-core model calculations, and combined with their near degeneracy in energy provide clear evidence for the absence of the N = 14 sub-shell closure. The very small spectroscopic factor found for the 3/2+ ground state is consistent with theoretical predictions and indicates that the {\nu}1d3/2 strength is carried by unbound states. With a dominant l = 0 valence neutron configuration and a very low separation energy, the 1/2+ excited state is a one-neutron halo candidate., Comment: 7 pages, 5 figures, Accepted for publication in Physics Letters B

We study the formation and evolution of a convective layer when a stably-stratified fluid with a composition gradient is cooled from above. We perform a series of 2D simulations using the Bousinessq approximation with Prandtl number ranging from Pr = 0.1 to 7, extending previous work on salty water to low Pr. We show that the evolution of the convection zone is well-described by an entrainment prescription in which a fixed fraction of the kinetic energy of convective motions is used to mix fluid at the interface with the stable layer. We measure the entrainment efficiency and find that it grows with decreasing Prandtl number or increased applied heat flux. The kinetic energy flux that determines the entrainment rate is a small fraction of the total convective luminosity. In this time-dependent situation, the density ratio at the interface is driven to a narrow range that depends on the value of Pr, and with low enough values that advection dominates the interfacial transport. We characterize the interfacial flux ratio and how it depends on the interface stability. We present an analytic model that accounts for the growth of the convective layer with two parameters, the entrainment efficiency and the interfacial heat transport, both of which can be measure from the simulations., Comment: Published in Physical Review Fluids. The published version has some minor typos. We fixed them in this version

Full detector simulation was among the largest CPU consumer in all CERN experiment software stacks for the first two runs of the Large Hadron Collider (LHC). In the early 2010's, the projections were that simulation demands would scale linearly with luminosity increase, compensated only partially by an increase of computing resources. The extension of fast simulation approaches to more use cases, covering a larger fraction of the simulation budget, is only part of the solution due to intrinsic precision limitations. The remainder corresponds to speeding-up the simulation software by several factors, which is out of reach using simple optimizations on the current code base. In this context, the GeantV R&D project was launched, aiming to redesign the legacy particle transport codes in order to make them benefit from fine-grained parallelism features such as vectorization, but also from increased code and data locality. This paper presents extensively the results and achievements of this R&D, as well as the conclusions and lessons learnt from the beta prototype., Comment: 34 pages, 26 figures, 24 tables

Open clusters are key targets for both Galaxy structure and evolution and stellar physics studies. Since \textit{Gaia} DR2 publication, the discovery of undetected clusters has proven that our samples were not complete. Our aim is to exploit the Big Data capabilities of machine learning to detect new open clusters in \textit{Gaia} DR2, and to complete the open cluster sample to enable further studies on the Galactic disc. We use a machine learning based methodology to systematically search in the Galactic disc, looking for overdensities in the astrometric space and identifying them as open clusters using photometric information. First, we use an unsupervised clustering algorithm, DBSCAN, to blindly search for these overdensities in \textit{Gaia} DR2 $(l,b,\varpi,\mu_{\alpha^*},\mu_\delta)$. After that, we use a deep learning artificial neural network trained on colour-magnitude diagrams to identify isochrone patterns in these overdensities, and to confirm them as open clusters. We find $582$ new open clusters distributed along the Galactic disc, in the region $|b| < 20$. We can detect substructure in complex regions, and identify the tidal tails of a disrupting cluster UBC~$274$ of $\sim 3$ Gyr located at $\sim 2$ kpc. Adapting the methodology into a Big Data environment allows us to target the search driven by physical properties of the open clusters, instead of being driven by its computational requirements. This blind search for open clusters in the Galactic disc increases in a $45\%$ the number of known open clusters., Comment: 11 pages, 9 figures, accepted for publication by Astronomy and Astrophysics (A&A) the 18th January, 2020. Tables 1 and 2 available at the CDS

Glitches are rare spin-up events that punctuate the smooth slow-down of the rotation of pulsars. For the Vela pulsar and PSR J0537-6910, the glitch sizes and the times between consecutive events have clear preferred scales, contrary to the handful of other pulsars with enough glitches for such a study. Moreover, PSR J0537-6910 is the only pulsar showing a strong positive correlation between the size of each glitch and the waiting time until the following one. We attempt to understand this behavior through a detailed study of the distributions and correlations of glitch properties for the eight pulsars with at least ten detected glitches. We model the distributions of glitch sizes and times between consecutive glitches for this sample and perform Monte Carlo simulations to explore two hypotheses that could explain why the correlation is so much weaker in other pulsars than in PSR J0537-6910. We confirm the above results for the Vela pulsar and PSR J0537-6910, and verify that the latter is the only pulsar with a strong correlation between glitch size and waiting time to the following glitch. For the remaining six pulsars, the waiting time distributions are best fitted by exponentials, and the size distributions either by power laws, exponentials, or log-normal functions. Some pulsars in the sample yield significant Pearson and Spearman coefficients ($r_p$ and $r_s$) for the aforementioned correlation. Moreover, for all except the Crab, both coefficients are positive. Our simulations show that the weaker correlations in pulsars other than PSR J0537$-$6910 cannot be due to missing glitches too small to be detected. We also tested the hypothesis that each pulsar may have two kinds of glitches, namely large, correlated ones and small, uncorrelated ones. We found for all pulsars that their glitches are not consistent with a pure uncorrelated population either.

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$., Comment: 18 pages, 12 figures

We present a statistical study of the glitch population and the behaviour of the glitch activity across the known population of neutron stars. An unbiased glitch database was put together based on systematic searches of radio timing data of 898 rotation-powered pulsars obtained with the Jodrell Bank and Parkes observatories. Glitches identified in similar searches of 5 magnetars were also included. The database contains 384 glitches found in the rotation of 141 of these neutron stars. We confirm that the glitch size distribution is at least bimodal, with one sharp peak at approximately $20\, \rm{\mu\,Hz}$, which we call large glitches, and a broader distribution of smaller glitches. We also explored how the glitch activity $\dot{\nu}_{\rm{g}}$, defined as the mean frequency increment per unit of time due to glitches, correlates with the spin frequency $\nu$, spin-down rate $|\dot{\nu}|$, and various combinations of these, such as energy loss rate, magnetic field, and spin-down age. It is found that the activity is insensitive to the magnetic field and that it correlates strongly with the energy loss rate, though magnetars deviate from the trend defined by the rotation-powered pulsars. However, we find that a constant ratio $\dot\nu_{\rm{g}}/|\dot\nu| = 0.010 \pm 0.001$ is consistent with the behaviour of all rotation-powered pulsars and magnetars. This relation is dominated by large glitches, which occur at a rate directly proportional to $|\dot{\nu}|$. The only exception are the rotation-powered pulsars with the highest values of $|\dot{\nu}|$, such as the Crab pulsar and PSR B0540$-$69, which exhibit a much smaller glitch activity, intrinsically different from each other and from the rest of the population. The activity due to small glitches also shows an increasing trend with $|\dot\nu|$, but this relation is biased by selection effects., Comment: Accepted for publication in A&A