Search

Chaos makes isolated systems of many interacting particles quickly thermalize and forget about their past. Here, we show that quantum mechanics hinders chaos in many-body systems: although the quantum eigenstates are thermal and strongly entangled, exponentially many of them are scarred, that is, have an enlarged weight along underlying classical unstable periodic orbits. Scarring makes the system more likely to be found on an orbit it was initialized on, retaining a memory of its past and thus weakly breaking ergodicity, even at long times and despite the system being fully thermal. We demonstrate the ubiquity of quantum scarring in many-body systems by considering a large family of spin models, including some of the most popular ones from condensed matter physics. Our findings, at hand for modern quantum simulators, prove structure in spite of chaos in many-body quantum systems., Comment: 5 pages, 3 figures

The accurate reconstruction of immiscible fluid-fluid interfaces from the volume fraction field is a critical component of geometric Volume of Fluid (VOF) methods. A common strategy is the Piecewise Linear Interface Calculation (PLIC), which fits a plane in each mixed-phase computational cell. However, recent work goes beyond PLIC by using two planes or even a paraboloid. To select such planes or paraboloids, complex optimization algorithms as well as carefully crafted heuristics are necessary. Yet, the potential exists for a well-trained machine learning model to efficiently provide broadly applicable solutions to the interface reconstruction problem at lower costs. In this work, the viability of a machine learning approach is demonstrated in the context of a single plane reconstruction. A feed-forward deep neural network is used to predict the normal vector of a PLIC plane given volume fraction and phasic barycenter data in a $3\times3\times3$ stencil. The PLIC plane is then translated in its cell to ensure exact volume conservation. Our proposed neural network PLIC reconstruction (PLIC-Net) is equivariant to reflections about the Cartesian planes. Training data is analytically generated with $\mathcal{O}(10^6)$ randomized paraboloid surfaces, which allows for sampling a wide range of interface shapes. PLIC-Net is tested in multiphase flow simulations where it is compared to standard (E)LVIRA reconstruction algorithms, and the impact of training data statistics on PLIC-Net's performance is also explored. It is found that PLIC-Net greatly limits the formation of spurious planes and generates cleaner numerical break-up of the interface. Additionally, the computational cost of PLIC-Net is lower than that of (E)LVIRA. These results establish that machine learning is a viable approach to VOF interface reconstruction and is superior to current reconstruction algorithms for some cases.

In the present paper, we model the velocity disturbance generated by a regularized forcing near a planar wall, which, along with the temporal nature of the forcing, provides an estimate of the unsteady velocity disturbance of the particle near a planar wall. We use the analytical solution for a singular in-time transient Stokeslet near a planar wall (Felderhof, 2009) and derive the corresponding time-persistent Stokeslets. The velocity disturbance due to a regularized forcing is then obtained numerically via a discrete convolution with the regularization kernel. The resulting Green's functions for parallel and perpendicular regularized forcing to the wall are stored as pre-computed temporal correction maps. By storing the time-dependent particle force on the fluid as fictitious particles, we estimate the unsteady velocity disturbance generated by the particle as a scalar product between the stored forces and the pre-computed Green's functions. Since the model depends on the analytical Green's function solution of the singular Stokeslet near a planar wall, the obtained velocity disturbance exactly satisfies the no-slip condition and does not require any fitted parameters to account for the rapid decay of the disturbance near the wall. The numerical evaluation of the convolution integral makes the present method suitable for arbitrary regularization kernels. Additionally, the generation of parallel and perpendicular correction maps enables to estimate the velocity disturbance due to particle motion in arbitrary directions relative to the local flow. The convergence of the method is studied on a fixed particle near a planar wall, and verification tests are performed on a settling particle parallel to a wall and a free-falling particle perpendicular to the wall., Comment: 40 pages

In volume-filtered Euler-Lagrange simulations of particle-laden flows, the fluid forces acting on a particle are estimated using reduced models, which rely on the knowledge of the local undisturbed flow for that particle. Since the two-way coupling between the particle and the fluid creates a local flow perturbation, the filtered fluid velocity interpolated to the particle location must be corrected prior to estimating the fluid forces, so as to subtract the contribution of this perturbation and recover the local undisturbed flow with good accuracy. In this manuscript, we present a new model for estimating a particle's self-induced flow disturbance that accounts for its transient development and for inertial effects related to finite particle Reynolds numbers. The model also does not require the direction of the momentum feedback to align with the direction of the particle's relative velocity, allowing force contributions other than the steady drag force to be considered. It is based upon the linearization of the volume-filtered equations governing the particle's self-induced flow disturbance, such that their solution can be expressed as a linear combination of regularized transient Stokeslet contributions. Tested on a range of numerical cases, the model is shown to consistently estimate the particle's self-induced flow disturbance with high accuracy both in steady and highly transient flow environments, as well as for finite particle Reynolds numbers.

We use time-resolved nonlinear pump--probe measurements to reveal features of semiconductor moir\'e materials not accessible to linear spectroscopy. With an intense, red-detuned pump pulse, we generate a high density of virtual excitons or exciton--polarons in various moir\'e minibands. A broadband probe pulse in turn measures the response of all optical resonances induced by the pump-generated excitations. We generically observe a coherent blue shift originating from contact-like exciton--exciton interactions. At charge neutrality, these measurements allow us to assess the spatial overlap between different optical excitations and to observe signatures of a bound biexciton state between two different moir\'e exciton modes. In contrast to electron doped monolayers, spatially confined moir\'e attractive polarons behave as an ensemble of non-interacting two-level emitters, exhibiting an electron-density-independent ac-Stark effect. Tuning the pump laser into resonance with the attractive polaron, we demonstrate the filling of the moir\'e lattice with localized polarons and thereby realize a nonequilibrium Bose--Fermi mixture in moir\'e flat bands.

The Local Volume Complete Cluster Survey (LoVoCCS) is an on-going program to observe nearly a hundred low-redshift X-ray-luminous galaxy clusters (redshifts $0.0310^{44}$ erg/s) with the Dark Energy Camera (DECam), capturing data in $u,g,r,i,z$ bands with a $5\sigma$ point source depth of approximately 25-26th AB magnitudes. Here, we map the aperture masses in 58 galaxy cluster fields using weak gravitational lensing. These clusters span a variety of dynamical states, from nearly relaxed to merging systems, and approximately half of them have not been subject to detailed weak lensing analysis before. In each cluster field, we analyze the alignment between the 2D mass distribution described by the aperture mass map, the 2D red-sequence (RS) galaxy distribution, and the brightest cluster galaxy (BCG). We find that the orientations of the BCG and the RS distribution are strongly aligned throughout the interiors of the clusters: the median misalignment angle is 19 deg within 2 Mpc. We also observe the alignment between the orientations of the RS distribution and the overall cluster mass distribution (by a median difference of 32 deg within 1 Mpc), although this is constrained by galaxy shape noise and the limitations of our cluster sample size. These types of alignment suggest long-term dynamical evolution within the clusters over cosmic timescales., Comment: 40 pages, 16 figures, 5 tables; revised and accepted for publication in ApJ

The volume-filtering of the Navier-Stokes equations allows to consider the effect that particles have on the fluid without further assumptions, but closures arise of which the implications are not fully understood. In the present paper, we carefully study every closure in the volume-filtered fluid momentum equation and investigate their impact on the momentum and energy transfer dependent on the filtering characteristics. We provide an analytical expression for the viscous closure that arises because filter and spatial derivative in the viscous term do not commute. An analytical expression for the regularization of the particle momentum source of a single sphere in the Stokes regime is derived. Furthermore, we propose a model for the subfilter stress tensor, which originates from filtering the advective term. The model for the subfilter stress tensor is shown to agree well with the subfilter stress tensor for small filter widths relative to the size of the particle. We show that the subfilter stress tensor requires modeling and should not be neglected. For small filter widths, we find that the commonly applied Gaussian regularization of the particle momentum source is a poor approximation of the spatial distribution of the particle momentum source, but for larger filter widths the spatial distribution approaches a Gaussian. Furthermore, we propose a modified advective term in the volume-filtered momentum equation that consistently circumvents the common stability issues observed at locally small fluid volume fractions and identify inconsistencies in previous studies of the phase-averaged kinetic energy of the volume-filtered fluid velocity. Finally, we propose a generally applicable form of the volume-filtered momentum equation and its closures based on clear and well-founded assumptions and propose guidelines for point-particle simulations based on the new findings.

Over the past decades, the volume-of-fluid (VOF) method has been the method of choice for simulating atomization processes, owing to its unique ability to discretely conserve mass. Current state-of-the-art VOF methods, however, rely on the piecewise-linear interface calculation (PLIC) to represent the interface used when calculating advection fluxes. This renders the estimated curvature of the transported interface zeroth-order accurate at best, adversely impacting the simulation of surface-tension-driven flows. In the past few years, there have been several attempts at using piecewise-parabolic interface approximations instead of piecewise-linear ones for computing advection fluxes, albeit all limited to two-dimensional cases or not inherently mass conservative. In this contribution, we present our most recent work on three-dimensional piecewise-parabolic interface reconstruction and apply it in the context of the VOF method. As a result of increasing the order of the interface representation, the reconstruction of the interface and the estimation of its curvature now become a single step instead of two separate ones. The performance of this new approach is assessed both in terms of accuracy and stability and compared to the classical PLIC-VOF approach on a range of canonical test-cases and cases of surface-tension-driven instabilities.

A precondition for the deployment of a Reinforcement Learning agent to a real-world system is to provide guarantees on the learning process. While a learning algorithm will eventually converge to a good policy, there are no guarantees on the performance of the exploratory policies. We study the problem of conservative exploration, where the learner must at least be able to guarantee its performance is at least as good as a baseline policy. We propose the first conservative provably efficient model-free algorithm for policy optimization in continuous finite-horizon problems. We leverage importance sampling techniques to counterfactually evaluate the conservative condition from the data self-generated by the algorithm. We derive a regret bound and show that (w.h.p.) the conservative constraint is never violated during learning. Finally, we leverage these insights to build a general schema for conservative exploration in DeepRL via off-policy policy evaluation techniques. We show empirically the effectiveness of our methods.

Using the photometric population prediction method {\bf Red Dragon}, we characterize the Red Sequence (RS) and Blue Cloud (BC) of DES galaxies in the COSMOS field. Red Dragon (RD) uses a redshift-evolving, error-corrected Gaussian mixture model to detail the distribution of photometric colors, smoothly parameterizing the two populations with relative weights, mean colors, intrinsic scatters, and inter-color correlations. This resulting fit of RS and BC yields RS membership probabilities $P_{\rm RS}$ for each galaxy. Even when training on only DES main bands $griz$, RD selects the quiescent population (defined here as galaxies with $\lg {\rm sSFR \cdot yr} < -11$) with $\gtrsim 90\%$ balanced accuracy out to $z=2$; augmenting with extended photometry from VIRCAM improves this accuracy to $\sim 95\%$ out to $z=3$. We measure redshift evolution of sSFR and galactic age in several stellar mass bins, finding that the BC is consistently more star-forming (by $\gtrsim 1~{\rm dex}$) and typically younger (by $\gtrsim 1~{\rm Gyr}$) than the RS (up to $z \sim 1.4$). This characterization of both RS and BC as functions of redshift and stellar mass improves our understanding of both populations and opens the door to more precise galaxy population characterization in future deep optical and IR systems., Comment: 16 pages, 19 figures; comments welcome

The locally noncentrosymmetric heavy fermion superconductor CeRh$_2$As$_2$ has attracted considerable interests due to its rich superconducting phases, accompanied by a quadrupole density wave and pronounced antiferromagnetic excitations. To understand the underlying physics, we here report measurements from high-resolution angle-resolved photoemission. Our results reveal fine splittings of the conduction bands related to the locally noncentrosymmetric structure, as well as a quasi-two-dimensional Fermi surface (FS) with strong $4f$ contributions. The FS exhibits nesting with an in-plane vector $({\pi}/a, {\pi}/a)$, which is facilitated by the van Hove singularity near $\bar X$ that arises from the characteristic conduction-$f$ hybridization. The FS nesting provides a natural explanation for the observed antiferromagnetic excitations at $({\pi}/a, {\pi}/a)$, which could be intimately connected to its unconventional superconductivity. Our experimental results are well supported by density functional theory plus dynamical mean field theory calculations, which can capture the strong correlation effects. Our study not only provides spectroscopic proof of the key factors underlying the field-induced superconducting transition, but also uncovers the critical role of FS nesting and lattice Kondo effect in the intertwined spin and charge fluctuations., Comment: v1 submitted on Sep 13th 2023

Massive halos hosting groups and clusters of galaxies imprint coherent, arcminute-scale features across the spectrophotometric sky, especially optical-IR clusters of galaxies, distortions in the sub-mm CMB, and extended sources of X-ray emission. Statistical modeling of such features often rely upon the evolving space-time density of dark matter halos -- the halo mass function (HMF) -- as a common theoretical ground for cosmological, astrophysical and fundamental physics studies. We propose a compact (eight parameter) representation of the HMF with readily interpretable parameters that stem from polynomial expansions, first in terms of log-mass, then expanding those coefficients similarly in redshift. We demonstrate good ($\sim \! 5\%$) agreement of this form, referred to as the dual-quadratic (DQ-HMF), with Mira-Titan N-body emulator estimates for halo masses above $10^{13.7} h^{-1} {\rm M}_\odot$ over the redshift range $0.1 < z < 1.5$, present best-fit parameters for a Planck 2018 cosmology, and present parameter variation in the $\sigma_8 - \Omega_{\rm m}$ plane. Convolving with a minimal mass-observable relation (MOR) yields closed-form expressions for counts, mean mass, and mass variance of cluster samples characterized by some observable property. Performing information-matrix forecasts of potential parameter constraints from existing and future surveys under different levels of systematic uncertainties, we demonstrate the potential for percent-level constraints on model parameters by an LSST-like optical cluster survey of 300,000 clusters and a richness-mass variance of $0.3^2$. Even better constraints could potentially be achieved by a survey with one-tenth the sample size but with a reduced selection property variance of $0.1^2$. Potential benefits and extensions to the basic MOR parameterization are discussed., Comment: 18 pages, 9 figures

A significant computational expense and source of numerical errors in front tracking is the remeshing of the triangulated front, required due to distortion and compaction of the front following the Lagrangian advection of its vertices. Additionally, in classic front tracking, the remeshing of the front mesh is required not only due to the deformation of the front shape, but also because the vertices of the front are translated in the direction tangential to the front, induced by the front advection. We present the normal-only advection (NOA) front-tracking method with the aim of preventing the tangential motion of the front vertices and the associated vertex clustering, in order to reduce the number of remeshing operations required to retain a high-quality triangulated interface. To this end, we reformulate the velocity used to advect the front at each discrete front-vertex position. The proposed method is validated and tested against the classic front-tracking method, comparing volume conservation, shape preservation, computational costs, and the overall need for front remeshing, as well as experimental results for canonical interfacial flows. The presented results demonstrate that the NOA front-tracking method leads to a typical reduction of remeshing operations by 80 % or more compared to the classic front-tracking method for well-resolved cases, and results in a smoother front mesh, which is essential for an accurate representation of the geometrical properties of the front. The volume conservation error is reduced by approximately one order of magnitude with the proposed method compared to the classic front-tracking method, at a similar computational cost.

Quantum many-body scars (QMBS) consist of a few low-entropy eigenstates in an otherwise chaotic many-body spectrum, and can weakly break ergodicity resulting in robust oscillatory dynamics. The notion of QMBS follows the original single-particle scars introduced within the context of quantum billiards, where scarring manifests in the form of a quantum eigenstate concentrating around an underlying classical unstable periodic orbit (UPO). A direct connection between these notions remains an outstanding problem. Here, we study a many-body spinor condensate that, owing to its collective interactions, is amenable to the diagnostics of scars. We characterize the system's rich dynamics, spectrum, and phase space, consisting of both regular and chaotic states. The former are low in entropy, violate the Eigenstate Thermalization Hypothesis (ETH), and can be traced back to integrable effective Hamiltonians, whereas most of the latter are scarred by the underlying semiclassical UPOs, while satisfying ETH. We outline an experimental proposal to probe our theory in trapped spin-1 Bose-Einstein condensates., Comment: 4 pages and 4 figures with SM (11 pages, 5 figures)

Laser induced shift of atomic states due to the AC-Stark effect has played a central role in cold-atom physics and facilitated their emergence as analog quantum simulators. Here, we explore this phenomena in an atomically thin layer of semiconductor MoSe$_2$, which we embedded in a heterostructure enabling charge tunability. Shining an intense pump laser with a small detuning from the material resonances, we generate a large population of virtual collective excitations, and achieve a regime where interactions with this background population is the leading contribution to the AC-Stark shift. Using this technique we study how itinerant charges modify -- and dramatically enhance -- the interactions between optical excitations. In particular, our experiments show that the interaction between attractive polarons could be more than an order of magnitude stronger than those between bare excitons.

We propose a novel model to obtain the subgrid-scale velocity in the context of large-eddy simulation (LES) of particle-laden turbulent flows, to recover accurate particle statistics. In the new wavelet enrichment model, the subgrid-scale velocity is discretized with a divergence-free wavelet vector basis, and the coefficients of the expansion are obtained by minimizing the squared error of the linearized subfilter Navier-Stokes equations (SFNSE). The compact support of the wavelet basis is exploited to achieve continuously varying subgrid-scale velocity statistics across the domain. The performance of the new wavelet enrichment model is evaluated in single-phase and particle-laden flow simulations, comparing the results with the results of direct numerical simulations (DNS). The simulations show that the model can generate inhomogeneous and anisotropic velocity statistics, accurate strain-rotation relations, and a good approximation of the kinetic energy spectrum of the corresponding DNS. Furthermore, the model significantly improves the prediction of the particle-pair dispersion, the clustering of the particles, and the turbulence modulation by particles in two-way coupled simulations. The proposed model recovers the most important interactions between fluid turbulence and the behavior of the particles, while maintaining the computational cost on the order of a LES.