Nonlinear dynamics of multidimensional electrostatic excitations in nonthermal plasmas

Electrostatic solitary waves in plasmas are the focus of many current studies of localized electrostatic disturbances in both laboratory and astrophysical plasmas. Motivated by recent experimental observations, in which electrostatic solitary structures were detected in laser-plasma experiments, we have undertaken an investigation of the nonlinear dynamics of plasma evolving in two dimensions, in the presence of excess superthermal background electrons. We investigate the effect of a magnetic field on weakly nonlinear ion-acoustic waves. Deviation from the Maxwellian distribution is effectively modelled by the kappa model. A linear dispersion relation is derived, and a decrease in frequency and phase speed in both parallel and perpendicular modes can be seen, which is due to excess superthermal electrons, and which is stronger in the upper mode, and hardly noticeable in the lower (acoustic) mode. We show that ion-acoustic solitary waves can be generated during the nonlinear evolution of a plasma fluid, and their nonlinear propagation is governed by a Zakharov-Kuznetsov (ZK) type equation. A multiple scales perturbation technique is used to derive the ZK equation. Shock excitations can be produced if we allow for dissipation in the model, resulting in a Zakharov-Kuznetsov Burgers type equation. Different types of shock solutions and solitary waves are obtained, depending on the relation between the system parameters, and the effect of these on electrostatic shock structures is investigated numerically. A parametric investigation is conducted into the role of plasma nonthermality and magnetic field strength. © 2013 IOP Publishing Ltd.

Micromaser spectrum and phase diffusion dynamics

We discuss a scheme to relate the phase diffusion dynamics of the micromaser field to the measured atomic population statistics. This can allow us to measure the linewidth of the micromaser spectrum and to solve a relevant decoherence problem in cavity quantum electrodynamics. The main steps are (i) a suitable preparation of the cavity field state to generate coherences, (ii) the transfer of information on the dynamics of field coherences to probe atoms by the action of an external resonant coherent field and (iii) the derivation of the phase diffusion rate, hence the micromaser linewidth, from the measured population statistics of the probe atoms. The method can be applied even in the presence of trapping states, where peculiar linewidth oscillations are expected for increasing pump rate, due to the quantum nature of the micromaser field.

How to measure the phase diffusion dynamics in the micromaser

We propose a realistic scheme for measuring the micromaser linewidth by monitoring the phase diffusion dynamics of the cavity field. Our strategy consists of exciting an initial coherent state with the same photon number distribution as the micromaser steady-state field, singling out a purely diffusive process in the system dynamics. After the injection of a counterfield, measurements of the population statistics of a probe atom allow us to derive the micromaser linewidth in all ranges of the relevant parameters, establishing experimentally the distinctive features of the micromaser spectrum due to the discreteness of the electromagnetic field.

Assessing the Nonequilibrium Thermodynamics in a Quenched Quantum Many-Body System via Single Projective Measurements

We analyze the nature of the statistics of the work done on or by a quantum many-body system brought out of equilibrium. We show that, for the sudden quench and for an initial state that commutes with the initial Hamiltonian, it is possible to retrieve the whole nonequilibrium thermodynamics via single projective measurements of observables. We highlight, in a physically clear way, the qualitative implications for the statistics of work coming from considering processes described by operators that either commute or do not commute with the unperturbed Hamiltonian of a given system. We consider a quantum many-body system and derive an expression that allows us to give a physical interpretation, for a thermal initial state, to all of the cumulants of the work in the case of quenched operators commuting with the unperturbed Hamiltonian. In the commuting case, the observables that we need to measure have an intuitive physical meaning. Conversely, in the noncommuting case, we show that, although it is possible to operate fully within the single-measurement framework irrespectively of the size of the quench, some difficulties are faced in providing a clear-cut physical interpretation to the cumulants. This circumstance makes the study of the physics of the system nontrivial and highlights the nonintuitive phenomenology of the emergence of thermodynamics from the fully quantum microscopic description. We illustrate our ideas with the example of the Ising model in a transverse field showing the interesting behavior of the high-order statistical moments of the work distribution for a generic thermal state and linking them to the critical nature of the model itself.

Robust nonadiabatic molecular dynamics for metals and insulators

We present a new formulation of the correlated electron-ion dynamics (CEID) scheme, which systematically improves Ehrenfest dynamics by including quantum fluctuations around the mean-field atomic trajectories. We show that the method can simulate models of nonadiabatic electronic transitions and test it against exact integration of the time-dependent Schrodinger equation. Unlike previous formulations of CEID, the accuracy of this scheme depends on a single tunable parameter which sets the level of atomic fluctuations included. The convergence to the exact dynamics by increasing the tunable parameter is demonstrated for a model two level system. This algorithm provides a smooth description of the nonadiabatic electronic transitions which satisfies the kinematic constraints (energy and momentum conservation) and preserves quantum coherence. The applicability of this algorithm to more complex atomic systems is discussed.

Later prehistoric vegetation dynamics and Bronze Age agriculture at Hobbister, Orkney, Scotland

The Bronze Age in Britain was a time of major social and cultural changes, reflected in the division of the landscape into field systems and the establishment of new belief systems and ritual practices. Several hypotheses have been advanced to explain these changes, and assessment of many of them is dependent on the availability of detailed palaeoenvironmental data from the sites concerned. This paper explores the development of a later prehistoric landscape in Orkney, where a Bronze Age field system and an apparently ritually-deposited late Bronze Age axe head are located in an area of deep blanket peat from which high-resolution palaeoenvironmental sequences have been recovered. There is no indication that the field system was constructed to facilitate agricultural intensification, and it more likely reflects a cultural response to social fragmentation associated with a more dispersed settlement pattern. There is evidence for wetter conditions during the later Bronze Age, and the apparent votive deposit may reflect the efforts of the local population to maintain community integrity during a time of perceptible environmental change leading to loss of farmland. The study emphasises the advantages of close integration of palaeoenvironmental and archaeological data for interpretation of prehistoric human activity. The palaeoenvironmental data also provide further evidence for the complexity of prehistoric woodland communities in Orkney, hinting at greater diversity than is often assumed. Additionally, differing dates for woodland decline in the two sequences highlight the dangers of over-extrapolation from trends observed in a single pollen profile, even at a very local scale.

Nonequilibrium Quantum Landauer Principle

Using the operational framework of completely positive, trace preserving operations and thermodynamic fluctuation relations, we derive a lower bound for the heat exchange in a Landauer erasure process on a quantum system. Our bound comes from a nonphenomenological derivation of the Landauer principle which holds for generic nonequilibrium dynamics. Furthermore, the bound depends on the nonunitality of dynamics, giving it a physical significance that differs from other derivations. We apply our framework to the model of a spin-1/2 system coupled to an interacting spin chain at finite temperature.

SIMLA: Simulating particle dynamics in intense laser and other electromagnetic fields via classical and quantum electrodynamics

We present the Fortran program SIMLA, which is designed for the study of charged particle dynamics in laser and other background fields. The dynamics can be determined classically via the Lorentz force and Landau–Lifshitz equations or, alternatively, via the simulation of photon emission events determined by strong-field quantum-electrodynamics amplitudes and implemented using Monte-Carlo routines. Multiple background fields can be included in the simulation and, where applicable, the propagation direction, field type (plane wave, focussed paraxial, constant crossed, or constant magnetic), and time envelope of each can be independently specified.

Pathways Project Final Report Volume 1: Field Investigation and Catchment Conceptual Models

viii
Executive Summary
The Pathways Project field studies were targeted at improving the understanding of contaminant transport along different hydrological pathways in Irish catchments, including their associated impacts on water quality and river ecology. The contaminants of interest were phosphorus, nitrogen and sediment. The working Pathways conceptual model included overland flow, interflow, shallow groundwater flow, and deep groundwater flow. This research informed the development of a set of Catchment Management Support Tools (CMSTs) comprising an Exploratory Tool, Catchment Characterization Tool (CCT) and Catchment Modelling Tool (CMT) as outlined in Pathways Project Final Reports Volumes 3 and 4.
In order to inform the CMST, four suitable study catchments were selected following an extensive selection process, namely the Mattock catchment, Co. Louth/Meath; Gortinlieve catchment, Co. Donegal; Nuenna catchment, Co. Kilkenny and the Glen Burn catchment, Co. Down. The Nuenna catchment is well drained as it is underlain by a regionally important karstified limestone aquifer with permeable limestone tills and gravels, while the other three catchments are underlain by poorly productive aquifers and low permeability clayey tills, and are poorly drained.
All catchments were instrumented, and groundwater, surface and near-surface water and aquatic ecology were monitored for a period of two years. Intensive water quality sampling during rainfall events was used to investigate the pathways delivering nutrients. The proportion of flow along each pathway was determined using chemical and physical hydrograph separation techniques, supported by numerical modelling.
The outcome of the field studies broadly supported the use of the initial four-pathway conceptual model used in the Pathways CMT (time-variant model). The artificial drainage network was found to be a significant contributing pathway in the poorly drained catchments, at low flows and during peak flows in wet antecedent conditions. The transition zone (TZ), i.e. the broken up weathered zone at the top of the bedrock, was also found to be an important pathway. It was observed to operate in two contrasting hydrogeological scenarios: in groundwater discharge zones the TZ can be regarded as being part of the shallow groundwater pathway, whereas in groundwater recharge zones it behaves more like interflow.
In the catchments overlying poorly productive aquifers, only a few fractures or fracture zones were found to be hydraulically active and the TZ, where present, was the main groundwater pathway. In the karstified Nuenna catchment, the springs, which are linked to conduits as well as to a diffuse fracture network, delivered the majority of the flow. These findings confirm the two-component groundwater contribution from bedrock but suggest that the size and nature of the hydraulically active fractures and the nature of the TZ are the dominant factors at the scale of a stream flow event.
Diffuse sources of nitrate were found to be typically delivered via the subsurface pathways, especially in the TZ and land drains in the poorly productive aquifer catchments, and via the bedrock groundwater in the Nuenna. Phosphorus was primarily transported via overland flow in both particulate and soluble forms. Where preferential flow paths existed in the soil and subsoil, soluble P, and to a lesser extent particulate P, were also transported via the TZ and in drains and ditches. Arable land was found to be the most important land use for
ix
the delivery of sediment, although channel bank and in-stream sources were the most significant in the Glen Burn catchment. Overland flow was found to be the predominant transport sediment pathway in the poorly productive catchments. These findings informed the development of the transport and attenuation equations used in the CCT and CMT. From an assessment of the relationship between physico-chemical and biological conditions, it is suggested that in the Nuenna, Glen Burn and Gortinlieve catchments, a relationship may exist between biological water quality and nitrogen concentrations, as well as with P. In the Nuenna, there was also a relationship between macroinvertebrate status and alkalinity.
Further research is recommended on the transport and delivery of phosphorus in groundwater, the transport and attenuation dynamics in the TZ in different hydrogeological settings and the relationship between macroinvertebrates and co-limiting factors. High resolution temporal and spatial sampling was found to be important for constraining the conceptual understanding of nutrient and sediment dynamics which should also be considered in future studies.

An evaluation of strain softening of glacial tills as a result of pore pressure dynamics.

Deep-seated progressive failures of cuttings in heavily overconsolidated clays have been observed in the field and are well documented, especially for London Clays (Potts, Kovacevic, & Vaughan, 1997; Smethurst, Powrie, & Clarke, 2006; Take, 2003), however, the process of softening and the development of a rupture surface in other clays, including the clay fraction of glacial tills, is still to be established. Recent decades have witnessed extreme weather conditions in Northern Ireland with dry summers and wet winters. The dynamics of this pore pressure variation can trigger strength reduction and progressive plastic straining, both of which will lead to slope failure. The aim of this research is to evaluate the effect of pore pressure variations on the deformation and long-term stability of large cuttings in glacial tills in Northern Ireland. This paper outlines the overall research program and presents initial laboratory findings (Carse, 2013).

Nonconservative current-driven dynamics: beyond the nanoscale

Long metallic nanowires combine crucial factors for nonconservative current-driven atomic motion. These systems have degenerate vibrational frequencies, clustered about a Kohn anomaly in the dispersion relation, that can couple under current to form nonequilibrium modes of motion growing exponentially in time. Such motion is made possible by nonconservative current-induced forces on atoms, and we refer to it generically as the waterwheel effect. Here the connection between the waterwheel effect and the stimulated directional emission of phonons propagating along the electron flow is discussed in an intuitive manner. Nonadiabatic molecular dynamics show that waterwheel modes self-regulate by reducing the current and by populating modes in nearby frequency, leading to a dynamical steady state in which nonconservative forces are counter-balanced by the electronic friction. The waterwheel effect can be described by an appropriate effective nonequilibrium dynamical response matrix. We show that the current-induced parts of this matrix in metallic systems are long-ranged, especially at low bias. This nonlocality is essential for the characterisation of nonconservative atomic dynamics under current beyond the nanoscale.

Dynamics and structure of self-generated magnetics fields on solids following high contrast, high intensity laser irradiation

The dynamics of self-generated magnetic B-fields produced following the interaction of a high contrast, high intensity (I > 10¹⁹W cm^-2) laser beam with thin (3 μm thick) solid (Al or Au) targets is investigated experimentally and numerically. Two main sources drive the growth of B-fields on the target surfaces. B-fields are first driven by laser-generated hot electron currents that relax over ∼10-20 ps. Over longer timescales, the hydrodynamic expansion of the bulk of the target into vacuum also generates B-field induced by non-collinear gradients of density and temperature. The laser irradiation of the target front side strongly localizes the energy deposition at the target front, in contrast to the target rear side, which is heated by fast electrons over a much larger area. This induces an asymmetry in the hydrodynamic expansion between the front and rear target surfaces, and consequently the associated B-fields are found strongly asymmetric. The sole long-lasting (>30 ps) B-fields are the ones growing on the target front surface, where they remain of extremely high strength (∼8-10 MG). These B-fields have been recently put by us in practical use for focusing laser-accelerated protons [B. Albertazzi et al., Rev. Sci. Instrum. 86, 043502 (2015)]; here we analyze in detail their dynamics and structure.

Electron-phonon thermalization in a scalable method for real-time quantum dynamics

We present a quantum simulation method that follows the dynamics of out-of-equilibrium many-body systems of electrons and oscillators in real time. Its cost is linear in the number of oscillators and it can probe time scales from attoseconds to hundreds of picoseconds. Contrary to Ehrenfest dynamics, it can thermalize starting from a variety of initial conditions, including electronic population inversion. While an electronic temperature can be defined in terms of a nonequilibrium entropy, a Fermi-Dirac distribution in general emerges only after thermalization. These results can be used to construct a kinetic model of electron-phonon equilibration based on the explicit quantum dynamics.

Hα Line Profile Asymmetries and the Chromospheric Flare Velocity Field

The asymmetries observed in the line profiles of solar flares can provide important diagnostics of the properties and dynamics of the flaring atmosphere. In this paper the evolution of the Hα and Ca ii λ8542 lines are studied using high spatial, temporal, and spectral resolution ground-based observations of an M1.1 flare obtained with the Swedish 1 m Solar Telescope. The temporal evolution of the Hα line profiles from the flare kernel shows excess emission in the red wing (red asymmetry) before flare maximum and excess in the blue wing (blue asymmetry) after maximum. However, the Ca ii λ8542 line does not follow the same pattern, showing only a weak red asymmetry during the flare. RADYN simulations are used to synthesize spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. We show that the red asymmetry observed in Hα is not necessarily associated with plasma downflows, and the blue asymmetry may not be related to plasma upflows. Indeed, we conclude that the steep velocity gradients in the flaring chromosphere modify the wavelength of the central reversal in the Hα line profile. The shift in the wavelength of maximum opacity to shorter and longer wavelengths generates the red and blue asymmetries, respectively.

Dielectrics in a time-dependent electric field: a real-time approach based on density-polarization functional theory

In the presence of a (time-dependent) macroscopic electric field the electron dynamics of dielectrics cannot be described by the time-dependent density only. We present a real-time formalism that has the density and the macroscopic polarization P as key quantities. We show that a simple local function of P already captures long-range correlation in linear and nonlinear optical response functions. Specifically, after detailing the numerical implementation, we examine the optical absorption, the second- and third-harmonic generation of bulk Si, GaAs, AlAs and CdTe at different level of approximation. We highlight links with ultranonlocal exchange-correlation functional approximations proposed within linear response time-dependent density functional theory framework.