Search

By combining ab initio downfolding with cluster dynamical mean-field theory, we study the degree of correlations in the low-temperature structures of the breathing-mode kagome van der Waals materials Nb$_3$(F,Cl,Br,I)$_8$. We show that the Coulomb correlation strength steadily increases from I to Br, Cl, and F, allowing us to identify Nb$_3$I$_8$ as a weakly correlated (obstructed atomic) insulator whose gap is only mildly affected by the local Coulomb interaction. Nb$_3$Br$_8$ and Nb$_3$Cl$_8$ are strongly correlated insulators, whose gaps are significantly influenced by Coulomb-induced vertex corrections. Nb$_3$F$_8$ is a prototypical bulk Mott-insulator whose gap is initially opened by strong correlation effects. Angle-resolved photoemission spectroscopy measurements comparing Nb$_3$Br$_8$ and Nb$_3$I$_8$ allow us to experimentally confirm these findings by revealing spectroscopic footprints of the degree of correlation. Our calculations further predict that the entire material family can be tuned into correlated charge-transfer or Mott-insulating phases upon electron or hole doping. Our magnetic property analysis additionally shows that inter-layer magnetic frustrations in the high-temperature phase drive the lattice phase transition to the low-temperature structures. The accompanying bilayer hybridization through inter-layer dimerization yields magnetic singlet ground states in the Cl, Br, and I compounds. Our findings establish Nb$_3$X$_8$ as a robust, versatile, and tunable class for van der Waals-based Coulomb and Mott engineering and allow us to speculate on the symmetry-breaking effects necessary for the recently observed Josephson diode effect in NbSe$_2$/Nb$_3$Br$_8$/NbSe$_2$ heterostructures., Comment: 23 pages, 18 figures

Two-dimensional (2D) transition metal dichalcogenides have emerged as a promising platform for next-generation optoelectronic and spintronic devices. Mechanical exfoliation using adhesive tape remains the dominant method for preparing 2D materials of highest quality, including transition metal dichalcogenides, but always results in small-sized flakes. This limitation poses a significant challenge for investigations and applications where large scale flakes are needed. To overcome these constraints, we explored the preparation of 2D WS2 and WSe2 using a recently developed kinetic in situ single-layer synthesis method (KISS). In particular, we focused on the influence of different substrates, Au and Ag, and chalcogen atoms, S and Se, on the yield and quality of the 2D films. The crystallinity and spatial morphology of the 2D films were characterized using optical microscopy and atomic force microscopy, providing a comprehensive assessment of exfoliation quality. Low-energy electron diffraction verified that there is no preferential orientation between the 2D film and the substrate, while optical microscopy revealed that WSe2 consistently outperformed WS2 in producing large monolayers, regardless of the substrate used. Finally, X-ray diffraction and X-ray photoelectron spectroscopy demonstrate that no covalent bonds are formed between the 2D material and the underlying substrate. These results identify KISS method as a non-destructive approach for a more scalable approach of high-quality 2D transition metal dichalcogenides., Comment: Article (13 pages, 5 figures) and supporting information (5 pages, 6 figures)

How turbulent convective fluctuations organise to form large-scale structures in planetary atmospheres remains a question that eludes quantitative answers. The assumption that this process is the result of an inverse cascade was suggested half a century ago in two-dimensional fluids, but its applicability to atmospheric and oceanic flows remains heavily debated, hampering our understanding of the energy balance in planetary systems. We show with direct numerical simulations of spatial resolutions of 122882 $\times$ 384 points that rotating and stratified flows can support a bidirectional cascade of energy, in three dimensions, with a ratio of Rossby to Froude numbers comparable to that of the Earth's atmosphere. Our results establish that in dry atmospheres spontaneous order can arise via an inverse cascade to the largest spatial scales.

Bilayer graphene twisted at the angle of about 1.1{\deg} better known as magic angle, exhibits ultra-flat moir\'e superlattice bands that are a source of highly-tunable, exotic quantum phenomena. Such phenomena, like superconductivity, correlated Mott-like insulating states and orbital ferromagnetism are driven by strong-correlation physics that defies classical theories. The inadequacy of such classical models and the lack of theoretical understanding of the recently observed exotic phenomena calls for revisiting the theory behind the material system and associating it with the observed behaviour. This article reviews the physics behind twisted bilayer graphene, focusing primarily on moir\'e physics and the importance of electronic (flat) band structure. In addition, this paper provides a brief overview of the emerging phenomena of correlated insulating states, superconductivity and orbital ferromagnetism. Finally, the most recent developments in controlling the interaction-driven states and tuning the electronic interactions are presented.

The transport of matter by turbulent flows plays an important role, in particular in a geophysical context. Here, we study the relative movement of groups of two (pairs) and four (tetrahedra) Lagrangian particles using direct numerical simulations of the stably-stratified Boussinsesq equations, with Brunt-V\"ais\"al\"a frequency $N$ and Coriolis parameter $f$. We cover regimes close to homogeneous isotropic turbulence, to flows dominated by stratification and rotation, keeping fixed the ratio $N/f = 5$. The flows studied are anisotropic, so the relative motion between two particles depends not only on the initial separation between the particles, but also on their orientation with respect to the vertical axis. In all cases considered, we demonstrate that the relative particle motion differs depending on whether dispersion is considered forward or backwards in time, although the asymmetry becomes less pronounced when stratification and rotation increase. On the other hand, the strong fluctuations in the dispersion between two particles become more extreme when $N$ and $f$ increase. We also find evidence for the formation of shear layers, which become more pronounced as $N$ and $f$ become larger. Finally, we show that the irreversibility on the dispersion of a set of particles forming initially a regular tetrahedron becomes weaker when the influence of stratification and rotation increase, a property that we relate to that of the perceived rate-of-strain tensor.

Turbulence in stratified and rotating turbulent flows is characterized by an interplay between waves and eddies, resulting in continuous exchanges between potential and kinetic energy. Here, we study how these processes affect the turbulent energy cascade from large to small scales, which manifests itself by an irreversible evolution of the relative kinetic energy between two tracer particles. We find that when $r_0$, the separation between particles, is below a characteristic length $\ell_t$, potential energy is on average transferred to kinetic energy, reducing time irreversibility, and conversely when $r_0 > \ell_t$. Our study reveals that the scale $\ell_t$ coincides with the buoyancy length scale $L_B$ over a broad range of configurations until a transitional wave-dominated regime is reached.

The diverse and intriguing phenomena observed in twisted bilayer systems, such as graphene and transition metal dichalcogenides, prompted new questions about the emergent effects that they may host. However, the practical challenge of realizing these structures on a scale large enough for spectroscopic investigation, remains a significant hurdle, resulting in a scarcity of direct measurements of the electronic band structure of twisted transition metal dichalcogenide bilayers. Here we present a systematic nanoscale angle-resolved photoemission spectroscopy investigation of bulk, single layer, and twisted bilayer WS2 with a small twist angle of 4.4{\deg}. The experimental results are compared with theoretical calculations based on density functional theory along the high-symmetry directions {\Gamma}-K and {\Gamma}-M. Surprisingly, the electronic band structure measurements suggest a structural relaxation occurring at 4.4{\deg} twist angle, and formation of large, untwisted bilayer regions., Comment: 23 pages, 11 figures

Simulating plasmas in the Hall-MagnetoHydroDynamics (Hall-MHD) regime represents a valuable {approach for the investigation of} complex non-linear dynamics developing in astrophysical {frameworks} and {fusion machines}. Taking into account the Hall electric field is {computationally very challenging as} it involves {the integration of} an additional term, proportional to $\bNabla \times ((\bNabla\times\mathbf{B})\times \mathbf{B})$ in the Faraday's induction {law}. {The latter feeds back on} the magnetic field $\mathbf{B}$ at small scales (between the ion and electron inertial scales), {requiring} very high resolution{s} in both space and time {in order to properly describe its dynamics.} The computational {advantage provided by the} kinetic Lattice Boltzmann (LB) approach is {exploited here to develop a new} code, the \textbf{\textsc{F}}ast \textbf{\textsc{L}}attice-Boltzmann \textbf{\textsc{A}}lgorithm for \textbf{\textsc{M}}hd \textbf{\textsc{E}}xperiments (\textsc{flame}). The \textsc{flame} code integrates the plasma dynamics in lattice units coupling two kinetic schemes, one for the fluid protons (including the Lorentz force), the other to solve the induction equation describing the evolution of the magnetic field. Here, the newly developed algorithm is tested against an analytical wave-solution of the dissipative Hall-MHD equations, pointing out its stability and second-order convergence, over a wide range of the control parameters. Spectral properties of the simulated plasma are finally compared with those obtained from numerical solutions from the well-established pseudo-spectral code \textsc{ghost}. Furthermore, the LB simulations we present, varying the Hall parameter, highlightthe transition from the MHD to the Hall-MHD regime, in excellent agreement with the magnetic field spectra measured in the solar wind.

The Ocean's turbulent energy cycle has a paradox; large-scale eddies under the control of Earth's rotation primarily transfer kinetic energy (KE) to larger scales via an inverse cascade, while a transfer to smaller scales is needed to accomplish dissipation. It has been argued, using numerical simulations, that fronts, waves and other turbulent structures can produce a forward cascade of KE toward dissipation scales. However, this forward cascade and its coexistence with known inverse cascade were not confirmed in observations. Here we present the first evidence of a dual KE cascade in the Ocean by analyzing velocity measurements from surface drifters released in the Gulf of Mexico. Our results show that KE is injected at two dominant scales and transferred to both large and small scales, with the downscale flux dominating at scales smaller than ~1-10km. The cascade rates are modulated seasonally, with stronger KE injection and forward transfer during winter., Comment: No comments

We observe the emergence of strong vertical drafts in direct numerical simulations of the Boussinesq equations in a range of parameters of geophysical interest. These structures, which appear intermittently in space and time, generate turbulence and enhance kinetic and potential energy dissipation, providing a possible explanation for the observed variability of the local energy dissipation in the bulk of oceanic flows, and the modulation of its probability distribution function. We show how, due to the extreme drafts, in runs with Froude numbers observable in geophysical scenarios, roughly 10% of the domain flow can account for up to 50% of the global volume dissipation, reminiscent of estimates based on oceanic models., Comment: 16 pages, 6 figures

Non-Gaussian statistics of large-scale fields are routinely observed in data from atmospheric and oceanic campaigns and global models. Recent direct numerical simulations (DNSs) showed that large-scale intermittency in stably stratified flows is due to the emergence of sporadic, extreme events in the form of bursts in the vertical velocity and the temperature. This phenomenon results from the interplay between waves and turbulent motions, affecting mixing. We provide evidence of the enhancement of the classical small-scale (or internal) intermittency due to the emergence of large-scale drafts, connecting large- and small-scale bursts. To this aim we analyze a large set of DNSs of the stably stratified Boussinesq equations over a wide range of values of the Froude number ($Fr\approx 0.01-1$). The variation of the buoyancy field kurtosis with $Fr$ is similar to (though with smaller values than) the kurtosis of the vertical velocity, both showing a non-monotonic trend. We present a mechanism for the generation of extreme vertical drafts and vorticity enhancements which follows from the exact equations for field gradients., Comment: 7 pages, 4 figures

Geophysical fluid flows are predominantly turbulent and often strongly affected by the Earth's rotation, as well as by stable density stratification. Using direct numerical simulations of forced Boussinesq equations, we study the influence of these effects on the motion of fluid particles, focusing on cases where the frequencies associated with rotation and stratification (RaS), $N$ and $f$ respectively, are held at a fixed ratio $N/f=5$. As the intensity of RaS increases, a sharp transition is observed between a regime dominated by eddies to a regime dominated by waves, which can also be seemingly described by simply comparing the time scale $1/N$ and $\tau_\eta$ (the Kolmogorov time scale). We perform a detailed study of Lagrangian statistics of acceleration, velocity and related quantities in the two regimes. The flow anisotropy induces a clear difference between particle motion in the horizontal and vertical directions. In the regime $N\tau_\eta <1$, acceleration statistics in both horizontal and vertical directions, exhibit well known characteristics of isotropic turbulence. In contrast for $N\tau_\eta >1$, they are directly influenced by imposed RaS. The Lagrangian velocity statistics exhibit visible anisotropy for all runs; nevertheless the degree of anisotropy becomes very strong in the regime $N\tau_\eta >1$. We find that in the regime $N\tau_\eta <1$, rotation enhances the mean displacement of particles in horizontal planes at short times, but inhibits them at longer times. This inhibition of horizontal displacement becomes stronger for $N\tau_\eta >1$, with no clear diffusive behavior. Displacements in the vertical direction are always inhibited. The inhibition becomes extremely strong when $N\tau_\eta >1$, with the particles almost being trapped horizontally., Comment: 24 pages, 10 figures

We investigate the large-scale intermittency of vertical velocity and temperature, and the mixing properties of stably stratified turbulent flows using both Lagrangian and Eulerian fields from direct numerical simulations, in a parameter space relevant for the atmosphere and the oceans. Over a range of Froude numbers of geophysical interest ($\approx 0.05-0.3$) we observe very large fluctuations of the vertical components of the velocity and the potential temperature, localized in space and time, with a sharp transition leading to non-Gaussian wings of the probability distribution functions. This behavior is captured by a simple model representing the competition between gravity waves on a fast time-scale and nonlinear steepening on a slower time-scale. The existence of a resonant regime characterized by enhanced large-scale intermittency, as understood within the framework of the proposed model, is then linked to the emergence of structures in the velocity and potential temperature fields, localized overturning and mixing. Finally, in the same regime we observe a linear scaling of the mixing efficiency with the Froude number and an increase of its value of roughly one order of magnitude.