Direct statistical simulation of the Lorenz63 system.

We use direct statistical simulation to find the low-order statistics of the well-known dynamical system, the Lorenz63 model. Instead of accumulating statistics from numerical simulation of the dynamical system or solving the Fokker-Planck equation for the full probability distribution of the dynamical system, we directly solve the equations of motion for the low-order statistics after closing them by making several different choices for the truncation. Fixed points of the statistics are obtained either by time evolving or by iterative methods. The stability and statistical realizability of the fixed points of the statistics are analyzed, and the statistics so obtained are compared to those found by the traditional approach. Low-order statistics of the chaotic Lorenz63 system can be obtained from cumulant expansions more efficiently than by accumulation via direct numerical simulation or by solution of the Fokker-Planck equation.

Topological Gaseous Plasmon Polariton in Realistic Plasma.

Nontrivial topology in bulk matter has been linked with the existence of topologically protected interfacial states. We show that a gaseous plasmon polariton (GPP), an electromagnetic surface wave existing at the boundary of magnetized plasma and vacuum, has a topological origin that arises from the nontrivial topology of magnetized plasma. Because a gaseous plasma cannot sustain a sharp interface with discontinuous density, one must consider a gradual density falloff with scale length comparable to or longer than the wavelength of the wave. We show that the GPP may be found within a gapped spectrum in present-day laboratory devices, suggesting that platforms are currently available for experimental investigation of topological wave physics in plasmas.

Topological origin of equatorial waves.

Topology sheds new light on the emergence of unidirectional edge waves in a variety of physical systems, from condensed matter to artificial lattices. Waves observed in geophysical flows are also robust to perturbations, which suggests a role for topology. We show a topological origin for two well-known equatorially trapped waves, the Kelvin and Yanai modes, owing to the breaking of time-reversal symmetry by Earth's rotation. The nontrivial structure of the bulk Poincaré wave modes encoded through the first Chern number of value 2 guarantees the existence of these waves. This invariant demonstrates that ocean and atmospheric waves share fundamental properties with topological insulators and that topology plays an unexpected role in Earth's climate system.

Statistics of the stochastically forced Lorenz attractor by the Fokker-Planck equation and cumulant expansions.

We investigate the Fokker-Planck description of the equal-time statistics of the three-dimensional Lorenz attractor with additive white noise. The invariant measure is found by computing the zero (or null) mode of the linear Fokker-Planck operator as a problem of sparse linear algebra. Two variants are studied: a self-adjoint construction of the linear operator and the replacement of diffusion with hyperdiffusion. We also access the low-order statistics of the system by a perturbative expansion in equal-time cumulants. A comparison is made to statistics obtained by the standard approach of accumulation via direct numerical simulation. Theoretical and computational aspects of the Fokker-Planck and cumulant expansion methods are discussed.

Generalized Quasilinear Approximation: Application to Zonal Jets.

Quasilinear theory is often utilized to approximate the dynamics of fluids exhibiting significant interactions between mean flows and eddies. We present a generalization of quasilinear theory to include dynamic mode interactions on the large scales. This generalized quasilinear (GQL) approximation is achieved by separating the state variables into large and small zonal scales via a spectral filter rather than by a decomposition into a formal mean and fluctuations. Nonlinear interactions involving only small zonal scales are then removed. The approximation is conservative and allows for scattering of energy between small-scale modes via the large scale (through nonlocal spectral interactions). We evaluate GQL for the paradigmatic problems of the driving of large-scale jets on a spherical surface and on the beta plane and show that it is accurate even for a small number of large-scale modes. As GQL is formally linear in the small zonal scales, it allows for the closure of the system and can be utilized in direct statistical simulation schemes that have proved an attractive alternative to direct numerical simulation for many geophysical and astrophysical problems.

Feedbacks and landscape-level vegetation dynamics.

Alternative stable-state theory (ASS) is widely accepted as explaining landscape-level vegetation dynamics, such as switches between forest and grassland. This theory argues that webs of feedbacks stabilise vegetation composition and structure, and that abrupt state shifts can occur if stabilising feedbacks are weakened. However, it is difficult to identify stabilising feedback loops and the disturbance thresholds beyond which state changes occur. Here, we argue that doing this requires a synthetic approach blending observation, experimentation, simulation, conceptual models, and narratives. Using forest boundaries and large mammal extinctions, we illustrate how a multifaceted research program can advance understanding of feedback-driven ecosystem change. Our integrative approach has applicability to other complex macroecological systems controlled by numerous feedbacks where controlled experimentation is impossible.

Direct statistical simulation of out-of-equilibrium jets.

We present direct statistical simulation of jet formation on a ß plane, solving for the statistics of a fluid flow via an expansion in cumulants. Here we compare an expansion truncated at second order (CE2) to statistics accumulated by direct numerical simulations. We show that, for jets near equilibrium, CE2 is capable of reproducing the jet structure (although some differences remain in the second cumulant). However, as the degree of departure from equilibrium is increased (as measured by the zonostrophy parameter), the jets meander more and CE2 becomes less accurate. We discuss a possible remedy by inclusion of higher cumulants.

Strong correlations in actinide redox reactions.

Reduction-oxidation (redox) reactions of the redox couples An(VI)/An(V), An(V)/An(IV), and An(IV)/An(III), where An is an element in the family of early actinides (U, Np, and Pu), as well as Am(VI)/Am(V) and Am(V)/Am(III), are modeled by combining density functional theory with a generalized Anderson impurity model that accounts for the strong correlations between the 5f electrons. Diagonalization of the Anderson impurity model yields improved estimates for the redox potentials and the propensity of the actinide complexes to disproportionate.

Weak ferromagnetic exchange and anomalous specific heat in ZnCu3(OH)6Cl2.

Experimental evidence for a plethora of low energy spin excitations in the spin-1/2 kagome antiferromagnet ZnCu3(OH)6Cl2 may be understandable in terms of an extended Fermi surface of spinons coupled to a U(1) gauge field. We carry out variational calculations to examine the possibility that such a state may be energetically viable. A Gutzwiller-projected wave function reproduces the dimerization of a kagome strip found previously by the density matrix renormalization group. Application to the full kagome lattice shows that the inclusion of a small ferromagnetic next-nearest-neighbor interaction favors a ground state with a spinon Fermi surface.

Broken time-reversal symmetry in strongly correlated ladder structures.

We provide, for the first time, in a doped strongly correlated system (two-leg ladder), a controlled theoretical demonstration of the existence of a state in which long-range ordered orbital currents are arranged in a staggered pattern, coexisting with a charge density wave. The method used is the highly accurate density-matrix renormalization group technique. This brings us closer to recent proposals that this order is realized in the enigmatic pseudogap phase of the cuprate high temperature superconductors.

Thermally enhanced neutralization in hyperthermal energy ion scattering.

Neutralization probabilities are presented for hyperthermal energy Na+ ions scattered from a Cu(001) crystal as a function of surface temperature and scattered velocity. A large enhancement in neutralization is observed as the temperature is increased. Velocity-dependent charge transfer regimes are probed by varying the incident energy, with the most prominent surface temperature effects occurring at the lowest energies. The data agree well with results obtained from a model based on the Newns-Anderson Hamiltonian, where the effects of both temperature and velocity are incorporated.

Staggered flux phase in a model of strongly correlated electrons.

We present numerical evidence for the existence of a staggered flux (SF) phase in the half-filled two-leg t-U-V-J ladder, with true long-range order in the countercirculating currents. The density-matrix renormalization-group finite-size scaling approach, generalized to describe complex-valued Hamil-tonians and wave functions, is employed. The SF phase exhibits robust currents at intermediate values of the interaction strength.

Time-dependent density-matrix renormalization group: a systematic method for the study of quantum many-body out-of-equilibrium systems.

The density-matrix renormalization-group algorithm is extended to treat time-dependent problems. The method provides a systematic and robust tool to explore out-of-equilibrium phenomena in quantum many-body systems. We illustrate the method by showing that attractive interactions enhance the tunneling current between two Luttinger liquids, whereas repulsive interactions suppress it, in qualitative agreement with analytical predictions. Enhancement of the transport current through a quantum dot in the Kondo regime is also exhibited.