Search

When soft matter is driven out of equilibrium its constituents interact via effective interactions that escape Newton's action-reaction principle. Prominent examples include the hydrodynamic interactions between colloidal particles driven in viscous fluids, phoretic interactions between chemically active colloids, and quorum sensing interactions in bacterial colonies. Despite a recent surge of interest in non-reciprocal physics a fundamental question remains : do non-reciprocal interactions alter or strengthen the ordered phases of matter driven out of equilibrium? Here, through a combination of experiments and simulations, we show how nonreciprocal forces propel and fission dislocations formed in hydrodynamically driven Wigner crystals. We explain how dislocation motility results in the continuous reshaping of grain-boundary networks, and how their fission reaction melts driven crystals from their interfaces. Beyond the specifics of hydrodynamics, we argue theoretically that topological defects and nonreciprocal interactions should invariably conspire to deform and ultimately destroy crystals whose the elementary units defy Newton's third law, Comment: 11 pages, 5 figures

Here we report a custom design setup for the simultaneous measurement of magneto-resistance and MOKE imaging in longitudinal configuration for characterization of magnetic thin-films and sensors. The setup is designed so as to cope with a small signal to noise ratio of initial images and to allow sensitive magnetoresistance (MR) measurement at low field. An improved differential algorithm is used to get a good enough contrast of the magnetic images and get rid of beam illumination or small camera gain fluctuations. Home made power supply and pre-amplifier stage were designed so as to reduce the low frequency noise as well as the thermal electrical noise. A vector magnet is used to produce rotational magnetic field so as to study the magnetic anisotropy and calculate the anisotropic constants in magnetic thin films. Magnetoresistive sensors patterned on epitaxial La$_{0.67}$Sr$_{0.33}$MnO$_3$ (LSMO) thin film have been characterized with this setup. The Images of the magnetization reversal process as well as the local magnetization loops deduced from these images provided evidence of a magnetic uniaxial anisotropy induced by the vicinal substrate. The magnetic anistropic constant of the films was then inferred from the MR measurements., Comment: 9 Pages, 9 Figures

Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multi-dimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localised flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20M_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for timescales typically considered by two-dimensional numerical simulations., Comment: 18 pages, 13 figures; to be published in MNRAS

Terahertz (THz) generation via photomixing on photoconductive antenna using twin delayed chirped pulses provides a long THz pulse with a narrow bandwidth. To generate a long pulse with broad bandwidth, we propose a new method that combines two long optical pulses with opposite chirps. The pulses exhibit temporal distributions of their instantaneous frequencies with opposite slopes. As a result, interaction between the beat frequency evolving over time and a photoconductor produces a broad THz spectrum with temporal variations. In our experimental setup, we generate a 12 ps-long pulse with a 1 THz bandwidth spectrum, featuring a frequency ramp of 90 GHz/ps, resembling a chirped THz pulse. This approach signifies a major advancement toward integrating photomixer technology, particularly in THz ranging applications.

Theoretical descriptions of convective overshooting often rely on a one-dimensional parameterization of the flow called the filling factor for convection. Several definitions of the filling factor have been developed, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. We study pairs of identical two- and three-dimensional global simulations of stars produced with MUSIC, a fully compressible, time-implicit hydrodynamics code. We examine (1) a $3 M_\odot$ red giant star near the first dredge-up point, (2) a $1 M_\odot$ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a $20 M_\odot$ main-sequence star with a large convective core. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making theoretical interpretation difficult. We present two new alternatives to these standard definitions, which compare flows at two different radial points. The first is the penetration parameter of Anders et al. (2022). The second is a new statistic, the plume interaction parameter. We demonstrate that both of these parameters capture systematic differences between 2D and 3D simulations around the convective boundary., Comment: 15 pages, 16 figures

As a massive star evolves along the main sequence, its core contracts, leaving behind a stable stratification in helium. We simulate 2D convection in the core at three different stages of evolution of a $5M_{\odot}$ star, with three different stratifications in helium atop the core. We study the propagation of internal gravity waves in the stably-stratified envelope, along with the overshooting length of convective plumes above the convective boundary. We find that the stratification in helium in evolved stars hinders radial motions and effectively shields the radiative envelope against plume penetration. This prevents convective overshooting from being an efficient mixing process in the radiative envelope. In addition, internal gravity waves are less excited in evolved models compared to the zero-age-main-sequence model, and are also more damped in the stratified region above the core. As a result, the wave power is several orders of magnitude lower in mid- and terminal-main-sequence models compared to zero-age-main-sequence stars.

Dry-snow slab avalanches are considered to be the most difficult to predict, yet the deadliest avalanche types. The release of snow slab avalanches starts with a initial failure in a weak layer that may propagate across the slope until the slab fractures and slides. The evaluation of crack propagation area is a primary concern for avalanche forecasters. The purpose of this study is to test the hypothesis that the heterogeneity of snowpack properties is one of the primary factors that may potentially stop dynamic crack propagation. To test this assumption, we use a depth-averaged Material Point Method (DA-MPM) for efficient elasto-plastic modeling of snow slab avalanches. Our analysis includes scenarios involving i) pure-elastic slabs and ii) elasto-plastic slabs. In the first scenario, we report a significant decrease in slab tensile stress with increasing crack speed compared to quasi-static theory. In addition, we quantify the effect of weak layer heterogeneity and softening fracture energy on the crack stopping mechanism. In the second scenario, we analyse the interplay between weak layer heterogeneity and slab tensile fracture and quantify their combined effect on crack arrest. Results are interpreted through a scaling law relating the crack arrest distance to two dimensionless numbers related to weak layer strength variability and slab tensile fracture. Furthermore, the proposed model is applied to field campaigns in which spatial variations of weak layer shear strength were measured. Finally, DA-MPM simulations are performed on three-dimensional terrain with spatial variations revealing interesting release patterns. This research and the proposed methods can not only enhance our comprehension of the factors influencing avalanche release sizes,and possibly, the design of new mitigation measures for avalanche start zones.

Two-dimensional and van der Waals ferromagnets are ideal platform to study low dimensional magnetism and proximity effects in van der Waals heterostructures. Their ultimate two dimensional character offers also the opportunity to easily adjust their magnetic properties using strain or electric fields. Among 2D ferromagnets, the Cr$_{1+x}$Te$_2$ compounds with $x$=0-1 are very promising because their magnetic properties depend on the amount of self-intercalated Cr atoms between pure CrTe$_2$ layers and the Curie temperature (T$_C$) can reach room temperature for certain compositions. Here, we investigate the evolution of the composition, structural and magnetic properties of thin Cr$_{1.33}$Te$_2$ (Cr$_2$Te$_3$) films epitaxially grown on graphene upon annealing. We observe a transition above 450{\deg}C from the Cr$_{1.33}$Te$_2$ phase with perpendicular magnetic anisotropy and a T$_C$ of 180 K to a composition close to Cr$_{1.39}$Te$_2$ with in-plane magnetic anisotropy and a T$_C$ of 240-250 K. This phase remains stable up to 650{\deg}C above which a pure Cr film starts to form. This work demonstrates the complex interplay between intercalated Cr, lattice parameters and magnetic properties in Cr$_{1+x}$Te$_2$ compounds., Comment: 5 pages, 5 figures

We present a novel subtraction method to remove the soft and collinear divergences at next-to-leading order for processes involving an arbitrary number of fragmentation functions, where this method acts directly in the hadronic centre-of-mass frame. We provide the analytical formulae of the subtraction terms in the general case where all the final state partons fragment to hadrons and for the two special cases when one of the partons of the final state does not fragment, i.e. it is a photon or involved in a jet.

Recent observations in emission, extinction, and polarisation have at least partially invalidated most of the astronomical standard grain models for the diffuse ISM. Moreover, lab measurements on interstellar silicate analogues have shown differences with the optical properties used in these standard models. To address these issues, our objective is twofold: (i) to update the optical properties of silicates and (ii) to develop the THEMIS dust model to allow the calculation of polarised extinction and emission. Based on optical constants measured in the lab for amorphous silicates and on observational constraints in mid-IR extinction and X-ray scattering, we defined new optical constants for the THEMIS silicates. Absorption and scattering efficiencies for spheroidal grains were then derived with the discrete dipole approximation. These new optical properties make it possible to explain the dust emission and extinction, both total and polarised. The model is not yet pushed to its limits since it does not require the perfect alignment of all grains to explain the observations and it therefore has the potential to accommodate the highest polarisation levels inferred from extinction measures. Moreover, the dispersion of the optical properties of the different lab silicates naturally explain the variations in both the total and polarised emission and extinction observed in the diffuse ISM. A single, invariant model calibrated on one single set of observations is obsolete for explaining contemporary observations. We are proposing a completely flexible dust model based entirely on lab measurements that has the potential to make major advances in understanding the nature of ISM grains and how they evolve as a function of their environment. Even if challenging, this is also relevant for future missions that will aim to perform precise measurements of the CMB spectral distortions and polarisation., Comment: Accepted for publication by A&A ; abstract modified to match the arxiv standard ; typo corrected on 19/01/2024 in Tab. 3 ; units corrected on 07/02/2024 in Figs. 9 and C.2

We investigate the numerical performance of a Discontinuous Galerkin (DG) hydrodynamics implementation when applied to the problem of driven, isothermal supersonic turbulence. While the high-order element-based spectral approach of DG is known to efficiently produce accurate results for smooth problems (exponential convergence with expansion order), physical discontinuities in solutions, like shocks, prove challenging and may significantly diminish DG's applicability to practical astrophysical applications. We consider whether DG is able to retain its accuracy and stability for highly supersonic turbulence, characterized by a network of shocks. We find that our new implementation, which regularizes shocks at sub-cell resolution with artificial viscosity, still performs well compared to standard second-order schemes for moderately high Mach number turbulence, provided we also employ an additional projection of the primitive variables onto the polynomial basis to regularize the extrapolated values at cell interfaces. However, the accuracy advantage of DG diminishes significantly in the highly supersonic regime. Nevertheless, in turbulence simulations with a wide dynamic range that start with supersonic Mach numbers and can resolve the sonic point, the low numerical dissipation of DG schemes still proves advantageous in the subsonic regime. Our results thus support the practical applicability of DG schemes for demanding astrophysical problems that involve strong shocks and turbulence, such as star formation in the interstellar medium. We also discuss the substantial computational cost of DG when going to high order, which needs to be weighted against the resulting accuracy gain. For problems containing shocks, this favours the use of comparatively low DG order., Comment: 22 pages, 14 figures, accepted for publication in MNRAS

Bound states in the continuum (BIC) are of special interest in photonics due to their theoretically infinite radiative lifetime. Here, we design a structure composed of a GaN layer with guided exciton-polaritons and a TiO$_2$ 1D photonic crystal slab. The photonic BIC hosted by the photonic crystal slab couples with the excitons of GaN to form a polaritonic BIC with a negative mass. This allows condensation to be reached with a low threshold in a structure suitable for electrical injection, paving the way for room-temperature polariton microdevices. We study in detail how the repulsive interaction between exciton-polaritons affects the condensate distribution in reciprocal space and, consequently, the condensate's overlap with the BIC resonance and, therefore, the condensate lifetime. We study an intrinsic contribution related to the formation of a bright soliton and the extrinsic contribution related to the interaction with an excitonic reservoir induced by spatially focused non-resonant pumping. We then study the peculiar dynamics of the condensation process in a BIC state for interacting particles using Boltzmann equations and hybrid Boltzmann-Gross Pitaevskii equations. We find optimal conditions allowing one to benefit from the long lifetime of the BIC for polariton condensation in a real structure.

So far, exciton-polariton (polariton) lasers were mostly single mode lasers based on microcavities. Despite the large repulsive polariton-polariton interaction, pulsed mode-locked polariton laser was never reported. Here, we use a $60 \ \mu m-$long GaN-based waveguide surrounded by distributed Bragg reflectors forming a multi-mode horizontal cavity. We demonstrate experimentally and theoretically a polariton mode-locked micro-laser operating in the blue-UV, at room temperature, with a 300$\ GHz$ repetition rate and 100$\ fs$-long pulses. The mode locking is demonstrated by the compensation (linearization) of the mode dispersion by the self-phase modulation induced by the polariton-polariton interaction. It is also supported by the observation in experiment and theory of the typical envelope frequency profile of a bright soliton., Comment: 16 pages, 12 figures

We present the characterisation of an inner mini-Neptune in a 9.2292005$\pm$0.0000063 day orbit and an outer mono-transiting sub-Saturn planet in a 95.50$^{+0.36}_{-0.25}$ day orbit around the moderately active, bright (mv=8.9 mag) K5V star TOI-2134. Based on our analysis of five sectors of TESS data, we determine the radii of TOI-2134b and c to be 2.69$\pm$0.16 R$_{e}$ for the inner planet and 7.27$\pm$0.42 R$_{e}$ for the outer one. We acquired 111 radial-velocity spectra with HARPS-N and 108 radial-velocity spectra with SOPHIE. After careful periodogram analysis, we derive masses for both planets via Gaussian Process regression: 9.13$^{+0.78}_{-0.76}$ M$_{e}$ for TOI-2134b and 41.86$^{+7.69}_{-7.83}$ M$_{e}$ for TOI-2134c. We analysed the photometric and radial-velocity data first separately, then jointly. The inner planet is a mini-Neptune with density consistent with either a water-world or a rocky core planet with a low-mass H/He envelope. The outer planet has a bulk density similar to Saturn's. The outer planet is derived to have a significant eccentricity of 0.67$^{+0.05}_{-0.06}$ from a combination of photometry and RVs. We compute the irradiation of TOI-2134c as 1.45$\pm$0.10 times the bolometric flux received by Earth, positioning it for part of its orbit in the habitable sone of its system. We recommend further RV observations to fully constrain the orbit of TOI-2134c. With an expected Rossiter-McLaughlin (RM) effect amplitude of 7.2$\pm$1.3 m/s, we recommend TOI-2134c for follow-up RM analysis to study the spin-orbit architecture of the system. We calculate the Transmission Spectroscopy Metric, and both planets are suitable for bright-mode NIRCam atmospheric characterisation., Comment: 24 pages, 7 tables, 14 figures

GaN/AlxGa1-xN quantum wells were grown by molecular beam epitaxy on high quality bulk (0001) GaN substrates. The quantum well thickness was set in the 6-8 nm range to favor the photoluminescence emission of indirect excitons. Indeed, such excitons are known to be spatially indirect, due to the presence of the internal electric field which spatially separates the electron and hole wave functions. The growth conditions were optimized in view of minimizing the photoluminescence peak broadening. In particular, the impact of growth temperature (up to 900{\deg}C) on the surface morphology, structural and photoluminescence properties was studied. The diffusion of indirect excitons on the scale of tens of microns was measured with a micro-photoluminescence setup equipped with a spatially resolved detection. A dedicated model and its analysis allow us to extract from these measurements the exciton diffusion constant and to conclude on the optimum growth conditions for the GaN/AlxGa1-xN quantum well structures suited for studies of quantum collective effects in indirect exciton liquids., Comment: 22 pages, 10 figures, 2 tables

Magnetic fields are ubiquitous in the Universe. Recently, cosmological simulations of galaxies have successfully begun to incorporate magnetic fields and their evolution in galaxies and their haloes. However, so far they have mostly focused on Milky Way-like galaxies. Here we analyse a sample of high resolution cosmological zoom simulations of disc galaxies in haloes with mass $M_\mathrm{200c}$ from $10^{10}\,\mathrm{M}_\odot$ to $10^{13}\,\mathrm{M}_\odot$, simulated with the Auriga galaxy formation model. We show that with sufficient numerical resolution the magnetic field amplification and saturation is converged. The magnetic field strength reaches equipartition with turbulent energy density for galaxies in haloes with $M_\mathrm{200c}\gtrsim 10^{11.5}\,\mathrm{M_\odot}$. For galaxies in less massive haloes, the magnetic field strength saturates at a fraction of equipartition that decreases with decreasing halo mass. For our lowest mass haloes, the magnetic field saturates significantly below $10\%$ of equipartition. We quantify the resolution we need to obtain converged magnetic field strengths and discuss our resolution requirements also in the context of the IllustrisTNG cosmological box simulations. We show that, at $z=0$, rotation-dominated galaxies in our sample exhibit for the most part an ordered large scale magnetic field, with fewer field reversals in more massive galaxies. Finally, we compare the magnetic fields in our cosmological galaxies at $z=0$ with simulations of isolated galaxies in a collapsing halo setup. Our results pave the way for detailed studies of cosmic rays and other physical processes in similar cosmological galaxy simulations that crucially depend on the strength and structure of magnetic fields., Comment: 19 pages, 12 figures, accepted for publication by MNRAS, comments welcome

We present highly sampled photometry of the supernova (SN) 2023ixf, a Type II SN in M101, beginning 2 days before its first known detection. To gather these data, we enlisted the global Unistellar Network of citizen scientists. These 252 observations from 115 telescopes show the SN's rising brightness associated with shock emergence followed by gradual decay. We measure a peak $M_{V}$ = -18.18 $\pm$ 0.09 mag at 2023-05-25 21:37 UTC in agreement with previously published analyses., Comment: 4 pages, 1 figure

In this Letter, we use a model fluid mechanics experiment to elucidate the impact of curvature heterogeneities on two-dimensional fields deriving from harmonic potential functions. This result is directly relevant to explain the smooth stationary structures in physical systems as diverse as curved liquid crystal and magnetic films, heat and Ohmic transport in wrinkled two-dimensional materials and flows in confined channels. Combining microfluidic experiments and theory, we explain how curvature heterogeneities shape confined viscous flows. We show that isotropic bumps induce local distortions to Darcy's flows whereas anisotropic curvature heterogeneities disturb them algebraically over system-spanning scales. Thanks to an electrostatic analogy, we gain insight into this singular geometric perturbation, and quantitatively explain it using both conformal mapping and numerical simulations. Altogether our findings establish the robustness of our experimental observations and their broad relevance to all Laplacian problems beyond the specifics of our fluid mechanics experiment.

Excitons hosted by GaN/(Al,Ga)N quantum wells (QWs) are spatially indirect due to the giant built-in electric field that separates electrons and holes along the growth direction. This electric field, and thus exciton energy, can be reduced by depositing metallic layers on the sample surface. Using spatially-resolved micro-photoluminescence spectroscopy we compare the effects of two different materials, Nickel/Gold (NiAu) and few-layered graphene (FLG), on the potential landscape experienced by the excitons. We are able to (i) determine the potential barriers imposed on QW excitons by deposition of FLG and NiAu to be $14$ and $82$~meV, respectively, and (ii) to evidence their impact on the exciton transport at appropriate densities. Optical losses and inhomogeneous broadening induced by deposition of NiAu and FLG layers are similar, and their joined implementation constitute a promising tool for electrostatic modulation of IX densities even in the absence of any applied electric bias.

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91\% of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to $10^4$ times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure., Comment: 15 pages, 10 figures, accepted for publication in MNRAS

User equipment (UE) positioning accuracy is of paramount importance in current and future communications standard. However, traditional methods tend to perform poorly in non line of sight (NLoS) scenarios. As a result, deep learning is a candidate to enhance the UE positioning accuracy in NLoS environments. In this paper, we study the efficiency of deep learning on the 3GPP indoor factory (InF) statistical channel. More specifically, we analyse the impacts of several key elements on the positioning accuracy: the type of radio data used, the number of base stations (BS), the size of the training dataset, and the generalization ability of a trained model., Comment: Accepted for publication at the European Conference on Networks and Communications (EuCNC) 2023

On 26 September 2022 the Double Asteroid Redirection Test (DART) spacecraft impacted Dimorphos, a satellite of the asteroid 65803 Didymos. Because it is a binary system, it is possible to determine how much the orbit of the satellite changed, as part of a test of what is necessary to deflect an asteroid that might threaten Earth with an impact. In nominal cases, pre-impact predictions of the orbital period reduction ranged from ~8.8 - 17.2 minutes. Here we report optical observations of Dimorphos before, during and after the impact, from a network of citizen science telescopes across the world. We find a maximum brightening of 2.29 $\pm$ 0.14 mag upon impact. Didymos fades back to its pre-impact brightness over the course of 23.7 $\pm$ 0.7 days. We estimate lower limits on the mass contained in the ejecta, which was 0.3 - 0.5% Dimorphos' mass depending on the dust size. We also observe a reddening of the ejecta upon impact., Comment: Accepted by Nature

Achieving large-scale growth of two-dimensional (2D) ferromagnetic materials with high Curie temperature (TC) and perpendicular magnetic anisotropy (PMA) is highly desirable for the development of ultra-compact magnetic sensors and magnetic memories. In this context, van der Waals (vdW) Cr2Te3 appears as a promising candidate. Bulk Cr2Te3 exhibits strong PMA and a TC of 180 K. Moreover, both PMA and TC might be adjusted in ultrathin films by engineering composition, strain, or applying an electric field. In this work, we demonstrate the molecular beam epitaxy (MBE) growth of vdW heterostructures of five-monolayer quasi-freestanding Cr2Te3 on three classes of 2D materials: graphene (semimetal), WSe2 (semiconductor) and Bi2Te3 (topological insulator). By combining structural and chemical analysis down to the atomic level with ab initio calculations, we confirm the single crystalline character of Cr2Te3 films on the 2D materials with sharp vdW interfaces. They all exhibit PMA and TC close to the bulk Cr2Te3 value of 180 K. Ab initio calculations confirm this PMA and show how its strength depends on strain. Finally, Hall measurements reveal a strong anomalous Hall effect, which changes sign at a given temperature. We theoretically explain this effect by a sign change of the Berry phase close to the Fermi level. This transition temperature depends on the 2D material in proximity, notably as a consequence of charge transfer. MBE-grown Cr2Te3/2D material bilayers constitute model systems for the further development of spintronic devices combining PMA, large spin-orbit coupling and sharp vdW interface., Comment: 18 pages, 13 figures and Supplemental Material

Topological insulators (TIs) are a promising class of materials for manipulating the magnetization of an adjacent ferromagnet (FM) through the spin-orbit torque (SOT) mechanism. However, current studies combining TIs with conventional FMs present large device-to-device variations, resulting in a broad distribution of SOT magnitudes. It has been identified that the interfacial quality between the TI and the FM is of utmost importance in determining the nature and efficiency of the SOT. To optimize the SOT magnitude and enable ultra-low-power magnetization switching, an atomically smooth interface is necessary. To this end, we have developed the growth of a full van der Waals FM/TI heterostructure by molecular beam epitaxy. The compensated TI (Bi0.4Sb0.6)2Te3 and ferromagnetic Fe3GeTe2 (FGT) were chosen because of their exceptional crystalline quality, low carrier concentration in BST and relatively large Curie temperature and perpendicular magnetic anisotropy in FGT. We characterized the magnitude of the SOTs by using thorough harmonic magnetotransport measurements and showed that the magnetization of an ultrathin FGT film could be switched with a current density Jc < 10^10 A/m^2. In comparison to previous studies utilizing traditional FMs, our findings are highly reliable, displaying little to no variation between devices.

Shallow landslides pose a significant threat to people and infrastructure. While often modeled based on limit equilibrium analysis, finite or discrete elements, continuum particle-based approaches like the Material Point Method (MPM) have more recently been successful in modeling their full 3D elasto-plastic behavior. In this paper, we develop a depth-averaged Material Point Method (DAMPM) to efficiently simulate shallow landslides over complex topography based on both material properties and terrain characteristics. DAMPM is an adaptation of MPM with classical shallow water assumptions, thus enabling large-deformation elasto-plastic modeling of landslides in a computationally efficient manner. The model is here demonstrated on the release of snow slab avalanches, a specific type of shallow landslides which release due to crack propagation within a weak layer buried below a cohesive slab. Here, the weak layer is considered as an external shear force acting at the base of an elastic-brittle slab. We validate our model against previous analytical calculations and numerical simulations of the classical snow fracture experiment known as Propagation Saw Test (PST). Furthermore, large scale simulations are conducted to evaluate the shape and size of avalanche release zones over different topographies. Given the low computational cost compared to 3D MPM, we expect our work to have important operational applications in hazard assessment, in particular for the evaluation of release areas, a crucial input for geophysical mass flow models. Our approach can be easily adapted to simulate both the initiation and dynamics of various shallow landslides, debris and lava flows, glacier creep and calving., Comment: 30 pages, 12 figures

We perform two-dimensional numerical simulations of core convection for zero-age-main-sequence stars covering a mass range from 3 $M_\odot$ to 20 $M_\odot$. The simulations are performed with the fully compressible time-implicit code MUSIC. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mass and luminosity. We also study the impact of artificially increasing the stellar luminosity for 3 $M_\odot$ models. The simulations cover hundreds to thousands of convective turnover timescales. Applying the framework of extreme plume events previously developed for convective envelopes, we derive overshooting lengths as a function of stellar masses. We find that the overshooting distance ($d_{\rm ov}$) scales with the stellar luminosity ($L$) and the convective core radius ($r_{\rm conv}$). We derive a scaling law $d_{\rm ov} \propto L^{1/3} r_{\rm conv}^{1/2}$ which is implemented in a 1D stellar evolution code and the resulting stellar models are compared to observations. The scaling predicts values for the overshooting distance that significantly increase with stellar mass, in qualitative agreement with observations. Quantitatively, however, the predicted values are underestimated for masses $\gtrsim 10 M_\odot$. Our 2D simulations show the formation of a nearly-adiabatic layer just above the Schwarzschild boundary of the convective core, as exhibited in recent 3D simulations of convection. The most luminous models show a growth in size with time of the nearly-adiabatic layer. This growth seems to slow down as the upper edge of the nearly-adiabatic layer gets closer to the maximum overshooting length and as the simulation time exceeds the typical thermal diffusive timescale in the overshooting layer., Comment: 13 pages, accepted for publication in MNRAS

Dust grains influence many aspects of star formation, including planet formation, opacities for radiative transfer, chemistry, and the magnetic field via Ohmic, Hall, and ambipolar diffusion. The size distribution of the dust grains is the primary characteristic influencing all these aspects. Grain size increases by coagulation throughout the star formation process. We describe here numerical simulations of protostellar collapse using methods described in earlier papers of this series. We compute the evolution of the grain size distribution from coagulation and the non-ideal magnetohydrodynamics effects self-consistently and at low numerical cost. We find that the coagulation efficiency is mostly affected by the time spent in high-density regions. Starting from sub-micron radii, grain sizes reach more than 100 {\mu}m in an inner protoplanetary disk that is only 1000 years old. We also show that the growth of grains significantly affects the resistivities, and indirectly the dynamics and angular momentum of the disk., Comment: 11 pages, 9 figures, accepted in A&A

An important class of nonlinear weighted least-squares problems arises from the assimilation of observations in atmospheric and ocean models. In variational data assimilation, inverse error covariance matrices define the weighting matrices of the least-squares problem. For observation errors, a diagonal matrix (i.e., uncorrelated errors) is often assumed for simplicity even when observation errors are suspected to be correlated. While accounting for observationerror correlations should improve the quality of the solution, it also affects the convergence rate of the minimization algorithms used to iterate to the solution. If the minimization process is stopped before reaching full convergence, which is usually the case in operational applications, the solution may be degraded even if the observation-error correlations are correctly accounted for. In this article, we explore the influence of the observation-error correlation matrix (R) on the convergence rate of a preconditioned conjugate gradient (PCG) algorithm applied to a one-dimensional variational data assimilation (1D-Var) problem. We design the idealised 1D-Var system to include two key features used in more complex systems: we use the background error covariance matrix (B) as a preconditioner (B-PCG); and we use a diffusion operator to model spatial correlations in B and R. Analytical and numerical results with the 1D-Var system show a strong sensitivity of the convergence rate of B-PCG to the parameters of the diffusion-based correlation models. Depending on the parameter choices, correlated observation errors can either speed up or slow down the convergence. In practice, a compromise may be required in the parameter specifications of B and R between staying close to the best available estimates on the one hand and ensuring an adequate convergence rate of the minimization algorithm on the other.

We model the coagulation and fragmentation of dust grains during the protostellar collapse with our newly developed shark code. It solves the gas-dust hydrodynamics in a spherical geometry and the coagulation/fragmentation equation. It also computes the ionization state of the cloud and the Ohmic, ambipolar and Hall resistivities. We find that the dust size distribution evolves significantly during the collapse, large grain formation being controlled by the turbulent differential velocity. When turbulence is included, only ambipolar diffusion remains efficient at removing the small grains from the distribution, brownian motion is only efficient as a standalone process. The macroscopic gas-dust drift is negligible for grain growth and only dynamically significant near the first Larson core. At high density, we find that the coagulated distribution is unaffected by the initial choice of dust distribution. Strong magnetic fields are found to enhance the small grains depletion, causing an important increase of the ambipolar diffusion. This hints that the magnetic field strength could be regulated by the small grain population during the protostellar collapse. Fragmentation could be effective for bare silicates, but its modeling relies on the choice of ill-constrained parameters. It is also found to be negligible for icy grains. When fragmentation occurs, it strongly affects the magnetic resistivities profiles. Dust coagulation is a critical process that needs to be fully taken into account during the protostellar collapse. The onset and feedback of fragmentation remains uncertain and its modeling should be further investigated., Comment: Accepted for publication in MNRAS

More than 5,000 exoplanets have been confirmed and among them almost 4,000 were discovered by the transit method. However, few transiting exoplanets have an orbital period greater than 100 days. Here we report a transit detection of Kepler-167 e, a "Jupiter analog" exoplanet orbiting a K4 star with a period of 1,071 days, using the Unistellar ground-based telescope network. From 2021 November 18 to 20, citizen astronomers located in nine different countries gathered 43 observations, covering the 16 hour long transit. Using a nested sampling approach to combine and fit the observations, we detected the mid-transit time to be UTC 2021 November 19 17:20:51 with a 1$\sigma$ uncertainty of 9.8 minutes, making it the longest-period planet to ever have its transit detected from the ground. This is the fourth transit detection of Kepler-167 e, but the first made from the ground. This timing measurement refines the orbit and keeps the ephemeris up to date without requiring space telescopes. Observations like this demonstrate the capabilities of coordinated networks of small telescopes to identify and characterize planets with long orbital periods., Comment: 11 pages, 3 figures, 2 tables, accepted in ApJL