Source author record

J. B. Marston

J. B. Marston appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Recent Developments in Theories of Inhomogeneous and Anisotropic Turbulence

Understanding inhomogeneous and anisotropic fluid flows require mathematical and computational tools that are tailored to such flows and distinct from methods used to understand the canonical problem of homogeneous and isotropic turbulence. We review some recent developments in the theory of inhomogeneous and anisotropic turbulence, placing special emphasis on several kinds of quasilinear approximations and their corresponding statistical formulations. Aspects of quasilinear theory that have received insufficient attention in the literature are discussed, and open questions are framed.

preprint2021arXiv

Direct statistical simulation of the Lorenz63 system

We use direct statistical simulation (DSS) 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 systems, we solve the equations of motion for the statistics themselves 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. Statistics so obtained are compared to those found by the traditional approach.

preprint2021arXiv

Signatures of electronic correlations and spin-susceptibility anisotropy in nuclear magnetic resonance

We present a methodology for probing the details of electronic susceptibility through minimally-invasive nuclear magnetic resonance techniques. Specifically, we classify electron-mediated long-range interactions in an ensemble of nuclear spins by revealing their effect on simple spin echo experiments. We find that pulse strength and applied field orientation dependence of these spin echo measurements resolves the spatial extent and anisotropy of electronic spin susceptibility. This work provides an alternate explanation to NMR results in superconducting and magnetically-ordered systems. The methodology has direct applications for sensing and characterizing emergent electronic phases.

preprint2020arXiv

Dimensional Reduction of Direct Statistical Simulation

Direct Statistical Simulation (DSS) solves the equations of motion for the statistics of turbulent flows in place of the traditional route of accumulating statistics by Direct Numerical Simulation (DNS). That low-order statistics usually evolve slowly compared with instantaneous dynamics is one important advantage of DSS. Depending on the symmetry of the problem and the choice of averaging operation, however, DSS is usually more expensive computationally than DNS because even low order statistics typically have higher dimension than the underlying fields. Here we show that it is possible to go much further by using Proper Orthogonal Decomposition (POD) to address the "curse of dimensionality." We apply POD directly to DSS in the form of expansions in the equal-time cumulants to second order (CE2). We explore two averaging operations (zonal and ensemble) and test the approach on two idealized barotropic models on a rotating sphere (a jet that relaxes deterministically towards an unstable profile, and a stochastically-driven flow that spontaneously organizes into jets). Order-of-magnitude savings in computational cost are obtained in the reduced basis, potentially enabling access to parameter regimes beyond the reach of DNS.

preprint2020arXiv

Nontrivial topology in the continuous spectrum of a magnetized plasma

Classification of matter through topological phases and topological edge states between distinct materials has been a subject of great interest recently. While lattices have been the main setting for these studies, a relatively unexplored realm for this physics is that of continuum fluids. In the typical case of a fluid model with a point spectrum, nontrivial topology and associated edge modes have been observed previously. However, another possibility is that a continuous spectrum can coexist with the point spectrum. Here we demonstrate that a fluid plasma model can harbor nontrivial topology within its continuous spectrum, and that there are boundary modes at the interface between topologically distinct regions. We consider the ideal magnetohydrodynamics (MHD) model. In the presence of magnetic shear, we find nontrivial topology in the Alfvén continuum. For strong shear, the Chern number is $\pm 1$, depending on the sign of the shear. If the magnetic shear changes sign within the plasma, a topological phase transition occurs, and bulk-boundary correspondence then suggests a mode localized to the layer of zero magnetic shear. We confirm the existence of this mode numerically. Moreover, this reversed-shear Alfvén eigenmode (RSAE) is well known within magnetic fusion as it has been observed in several tokamaks. In examining how the MHD model might be regularized at small scales, we also consider the electron limit of Hall MHD. We show that the whistler band, which plays an important role in planetary ionospheres, has nontrivial topology. More broadly, this work raises the possibility that fusion devices could be carefully tailored to produce other topological states with potentially useful behavior.

preprint2020arXiv

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 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.

preprint2019arXiv

Reduced-Order Quasilinear Model of Ocean Boundary-Layer Turbulence

The combined effectiveness of model reduction and the quasilinear approximation for the reproduction of the low-order statistics of oceanic surface boundary-layer turbulence is investigated. Idealized horizontally homogeneous problems of surface-forced thermal convection and Langmuir turbulence are studied in detail. Model reduction is achieved with a Galerkin projection of the governing equations onto an subset of modes determined by proper orthogonal decomposition. For less than 0.2% of the modes retained, the reduced quasilinear model is able to reproduce vertical profiles of horizontal mean fields as well as certain energetically important second-order turbulent transport statistics and energies to within 30% error. For intermediate sizes of the basis truncation some statistics approach those found in the fully nonlinear simulations. Thus basis reduction can actually improve upon the accuracy of quasilinear dynamics. A test model with a small total number of modes demonstrates the non-monotonic convergence toward the correct statistics as the size of the basis is increased.

preprint2016arXiv

Cumulant expansions for atmospheric flows

The equations governing atmospheric flows are nonlinear. Consequently, the hierarchy of cumulant equations is not closed. But because atmospheric flows are inhomogeneous and anisotropic, the nonlinearity may manifest itself only weakly through interactions of mean fields with disturbances such as thermals or eddies. In such situations, truncations of the hierarchy of cumulant equations hold promise as a closure strategy. We review how truncations at second order can be used to model and elucidate the dynamics of atmospheric flows. Two examples are considered. First, we study the growth of a dry convective boundary layer, which is heated from below, leading to turbulent upward energy transport and growth of the boundary layer. We demonstrate that a quasilinear truncation of the equations of motion, in which interactions of disturbances among each other are neglected but interactions with mean fields are taken into account, can capture the growth of the convective boundary layer even if it does not capture important turbulent transport terms. Second, we study the evolution of two-dimensional large-scale waves representing waves in Earth's upper atmosphere. We demonstrate that a cumulant expansion truncated at second order (CE2) can capture the evolution of such waves and their nonlinear interaction with the mean flow in some circumstances, for example, when the wave amplitude is small enough or the planetary rotation rate is large enough. However, CE2 fails to capture the flow evolution when nonlinear eddy--eddy interactions in surf zones become important. Higher-order closures can capture these missing interactions. The results point to new ways in which the dynamics of turbulent boundary layers may be represented in climate models, and they illustrate different classes of nonlinear processes that can control wave dissipation and momentum fluxes in the troposphere.

preprint2016arXiv

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-63 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. 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.

preprint2016arXiv

The 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. In this paper 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 non-local 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 this approximation 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.

preprint2015arXiv

On the Late-Time Behaviour of a Bounded, Inviscid Two-Dimensional Flow

Using complementary numerical approaches at high resolution, we study the late-time behaviour of an inviscid, incompressible two-dimensional flow on the surface of a sphere. Starting from a random initial vorticity field comprised of a small set of intermediate wavenumber spherical harmonics, we find that -- contrary to the predictions of equilibrium statistical mechanics -- the flow does not evolve into a large-scale steady state. Instead, significant unsteadiness persists, characterised by a population of persistent small-scale vortices interacting with a large-scale oscillating quadrupolar vorticity field. Moreover, the vorticity develops a stepped, staircase distribution, consisting of nearly homogeneous regions separated by sharp gradients. The persistence of unsteadiness is explained by a simple point vortex model characterising the interactions between the four main vortices which emerge.

preprint2014arXiv

Energy and Charge Transfer for Na+ Ions Scattered from a Ag(001) Surface

We present energy- and charge-resolved measurements of low and hyperthermal energy Na+ ions scattered from a Ag(001) surface. With the primary ion beam oriented along the [110] crystal direction, distinct peaks in the energy distributions of the scattered beam flux that correspond to single or multiple collisions with target atoms are observed. A classical trajectory simulation reveals that these collisions can occur either at the surface or within the [110] channels, depending on incident beam energy. Within the simulation we probe the role of finite temperature and thermally displaced atoms on specific scattering events and show that contributions to the scattered distributions from single and double collisions dominate within the [110] channels. We also report velocity dependent measurements of the neutral/ion ratio of the scattered beam flux. A deviation between the data and simulated charge transfer results is observed for Na trajectories which penetrate the surface.

preprint2014arXiv

Hyperviscosity and statistical equilibria of Euler turbulence on the torus and the sphere

Coherent structures such as jets and vortices appear in two-dimensional (2D) turbulence. To gain insight into both numerical simulation and equilibrium statistical mechanical descriptions of 2D Euler flows, the Euler equation with added hyperviscosity is integrated forward in time on the square torus and on the sphere. Coherent structures that form are compared against a hierarchy of truncated Miller-Robert-Sommeria equilibria. The energy-circulation-enstrophy MRS-2 description produces a complete condensation of energy to the largest scales, and in the absence of rotation correctly predicts the number and polarity of coherent vortices. Perturbative imposition of the quartic Casimir constraint improves agreement with numerical simulation by sharpening the cores and transferring some energy to smaller-scale modes. MRS-2 cannot explain qualitative changes due to rotation, but descriptions that conserve higher Casimirs beyond enstrophy have the potential to do so. The result is in agreement with the somewhat paradoxical observation that hyperviscosity helps to remedy the non-conservation of the third and higher Casimirs in numerical simulation. For a rotating sphere, numerical simulation also demonstrates that coherent structures found at late times depend on initial conditions, limiting the usefulness of statistical mechanics.

preprint2012arXiv

Rectification by Doped Mott-Insulator Junctions

Junctions of doped Mott insulators offer a route to rectification at frequencies beyond the terahertz range. Mott insulators have strong electronic correlations and therefore short timescales for electron-electron scattering. It is this short time scale that allows for the possibility of rectification at frequencies higher than those of semiconductor devices that are limited by the slow diffusion of charge carriers. We model a junction by a one dimensional chain of electrons with p- and n-doping on the two halves of the chain. Two types of systems are investigated: spin polarized electrons with nearest-neighbor interaction, and spin-half electrons that interact via on-site repulsion (the Hubbard model). For short chains the many-body Schrodinger equation can be integrated numerically exactly, and when driven by an oscillating electromagnetic field such idealized junctions rectify, showing a preferred direction for charge transfer. Longer chains are studied by the time-dependent density-matrix renormalization-group method, and also shown to rectify.

preprint2011arXiv

Planetary Atmospheres as Non-Equilibrium Condensed Matter

Planetary atmospheres, and models of them, are discussed from the viewpoint of condensed matter physics. Atmospheres are a form of condensed matter, and many interesting phenomena of condensed matter systems are realized by them. The essential physics of the general circulation is illustrated with idealized 2-layer and 1-layer models of the atmosphere. Equilibrium and non-equilibrium statistical mechanics are used to directly ascertain the statistics of these models.

preprint2011arXiv

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.

preprint2010arXiv

Astrophysical Fluid Dynamics via Direct Statistical Simulation

In this paper we introduce the concept of Direct Statistical Simulation (DSS) for astrophysical flows. This technique may be appropriate for problems in astrophysical fluids where the instantaneous dynamics of the flows are of secondary importance to their statistical properties. We give examples of such problems including mixing and transport in planets, stars and disks. The method is described for a general set of evolution equations, before we consider the specific case of a spectral method optimised for problems on a spherical surface. The method is illustrated for the simplest non-trivial example of hydrodynamics and MHD on a rotating spherical surface. We then discuss possible extensions of the method both in terms of computational methods and the range of astrophysical problems that are of interest.

preprint2010arXiv

Statistics of the General Circulation from Cumulant Expansions

Large-scale atmospheric flows may not be so nonlinear as to preclude their statistical description by systematic expansions in cumulants. I extend previous work by examining a two-layer baroclinic model of the general circulation. The fixed point of the cumulant expansion describes the statistically steady state of the out-of-equilibrium model. Equal-time statistics so obtained agree well with those accumulated by direct numerical simulation.

preprint2007arXiv

Statistics of an Unstable Barotropic Jet from a Cumulant Expansion

Low-order equal-time statistics of a barotropic flow on a rotating sphere are investigated. The flow is driven by linear relaxation toward an unstable zonal jet. For relatively short relaxation times, the flow is dominated by critical-layer waves. For sufficiently long relaxation times, the flow is turbulent. Statistics obtained from a second-order cumulant expansion are compared to those accumulated in direct numerical simulations, revealing the strengths and limitations of the expansion for different relaxation times.

preprint2000arXiv

$κ-(BEDT-TTF)_2X$ organic crystals: superconducting versus antiferromagnetic instabilities in an anisotropic triangular lattice Hubbard model

A Hubbard model at half-filling on an anisotropic triangular lattice has been proposed as the minimal model to describe conducting layers of $κ-(BEDT-TTF)_2X$ organic materials. The model interpolates between the square lattice and decoupled chains. The $κ-(BEDT-TTF)_2X$ materials present many similarities with cuprates, such as the presence of unconventional metallic properties and the close proximity of superconducting and antiferromagnetic phases. As in the cuprates, spin fluctuations are expected to play a crucial role in the onset of superconductivity. We perform a weak-coupling renormalization-group analysis to show that a superconducting instability occurs. Frustration in the antiferromagnetic couplings, which arises from the underlying geometrical arrangement of the lattice, breaks the perfect nesting of the square lattice at half-filling. The spin-wave instability is suppressed and a superconducting instability predominates. For the isotropic triangular lattice, there are again signs of long-range magnetic order, in agreement with studies at strong-coupling.

preprint2000arXiv

Multidimensional Bosonization

Bosonization of degenerate fermions yields insight both into Landau Fermi liquids, and into non-Fermi liquids. We begin our review with a pedagogical introduction to bosonization, emphasizing its applicability in spatial dimensions greater than one. After a brief historical overview, we present the essentials of the method. Well known results of Landau theory are recovered, demonstrating that this new tool of many-body theory is robust. Limits of multidimensional bosonization are tested by considering several examples of non-Fermi liquids, in particular the composite fermion theory of the half-filled Landau level. Nested Fermi surfaces present a different challenge, and these may be relevant in the cuprate superconductors. We conclude by discussing the future of multidimensional bosonization.

preprint1995arXiv

Theory of Fermion Liquids

We develop a general theory of fermion liquids in spatial dimensions greater than one. The principal method, bosonization, is applied to the cases of short and long range longitudinal interactions, and to transverse gauge interactions. All the correlation functions of the system may be obtained with the use of a generating functional. Short-range and Coulomb interactions do not destroy the Landau Fermi fixed point. Novel fixed points are found, however, in the cases of a super-long range longitudinal interaction in two dimensions and transverse gauge interactions in two and three spatial dimensions. We consider in some detail the 2+1-dimensional problem of a Chern-Simons gauge action combined with a longitudinal two-body interaction $V({\bf q}) \propto |{\bf q}|^{y-1}$ which controls the density, and hence gauge, fluctuations. For $y < 0$ we find that the gauge interaction is irrelevant and the Landau fixed point is stable, while for $y > 0$ the interaction is relevant and the fixed point cannot be accessed by bosonization. Of special importance is the case $y = 0$ (Coulomb interaction) which describes the Halperin-Lee-Read theory of the half-filled Landau level. We obtain the full quasiparticle propagator which is of a marginal Fermi liquid form. Using Ward Identities, we show that neither the inclusion of nonlinear terms in the fermion dispersion, nor vertex corrections, alters our results: the fixed point is accessible by bosonization. As the two-point fermion Green's function is not gauge invariant, we also investigate the gauge-invariant density response function. Near momentum $Q = 2 k_F$, in addition to the Kohn anomaly we find singular behavior. In Appendices we present a numerical calculation of the spectral function for a Fermi liquid with Landau parameter $f_0 \neq 0$. We also show how Kohn's theorem is

preprint1993arXiv

Coulomb Interaction and the Fermi Liquid State: Solution by Bosonization

We investigate the effects of the Coulomb two-body interaction on Fermi liquids via bosonization. The Coulomb interaction is singular in the limit of low momentum transfer, and recent interest in the possibility that some singular interactions might destroy the Fermi liquid state motivate us to reexamine it. We calculate the exact boson correlation function to show that the Fermi liquid state is retained in the case of Coulomb interactions. Spin and charge degrees of freedom propagate together at the same velocity and collective charge excitations (plasmons) exhibit the expected energy gap in three dimensions. Non-Fermi liquid behavior occurs, however, for a super long range interaction studied recently by Bares and Wen.