Variational Formulation of Macro-Particle Models for Electromagnetic Plasma Simulations

A variational method is used to derive a self-consistent macro-particle model for relativistic electromagnetic kinetic plasma simulations. Extending earlier work [E. G. Evstatiev and B. A. Shadwick, J. Comput. Phys., vol. 245, pp. 376-398, 2013], the discretization of the electromagnetic Low Lagrangian is performed via a reduction of the phase-space distribution function onto a collection of finite-sized macro-particles of arbitrary shape and discretization of field quantities onto a spatial grid. This approach may be used with both lab frame coordinates or moving window coordinates; the latter can greatly improve computational efficiency for studying some types of laser-plasma interactions. The primary advantage of the variational approach is the preservation of Lagrangian symmetries, which in our case leads to energy conservation and thus avoids difficulties with grid heating. Additionally, this approach decouples particle size from grid spacing and relaxes restrictions on particle shape, leading to low numerical noise. The variational approach also guarantees consistent approximations in the equations of motion and is amenable to higher order methods in both space and time. We restrict our attention to the 1-1/2 dimensional case (one coordinate and two momenta). Simulations are performed with the new models and demonstrate energy conservation and low noise.Comment: IEEE Transaction on Plasma Science (TPS) Special Issue: Plenary and Invited Papers of the Pulsed Power and Plasma Science Conference (PPPS 2013

Late-summer biogeochemistry in the Mertz Polynya: East Antarctica

A marked reconfiguration of the Mertz Polynya following the 2010 calving of the Mertz Glacier Tongue has been associated with a decrease in the size and activity of the polynya. We report observations of the oceanic carbonate (CO2) system in late-summer 2013, the third post-calving summer season. Estimates of seasonal net community production (NCP) based on inorganic carbon deficits and the oxygen-argon ratio indicate that the waters on the shelf to the east of Commonwealth Bay (adjacent to the Mertz Glacier) remain productive compared to pre-calving conditions. The input of residual or excess alkalinity from melting sea ice is found to contribute to the seasonal enhancement of carbonate saturation state and pH in shelf waters. Mean rates of NCP in 2012-2013 are more than twice as large as those observed in the pre-calving summers of 2001 and 2008 and suggest that the new (post-calving) configuration of the polynya favors enhanced net community production and a stronger surface ocean sink for atmospheric CO2 due at least in part to the redistribution of sea ice and associated changes in summer surface stratification

Seasonality of biological and physical controls on surface ocean CO2 from hourly observations at the Southern Ocean Time Series site south of Australia

The Subantarctic Zone (SAZ), which covers the northern half of the Southern Ocean between the Subtropical and Subantarctic Fronts, is important for air-sea CO2 exchange, ventilation of the lower thermocline, and nutrient supply for global ocean productivity. Here we present the first high-resolution autonomous observations of mixed layer CO2 partial pressure (pCO(2)) and hydrographic properties covering a full annual cycle in the SAZ. The amplitude of the seasonal cycle in pCO(2) (similar to 60 mu atm), from near-atmospheric equilibrium in late winter to similar to 330 mu atm in midsummer, results from opposing physical and biological drivers. Decomposing these contributions demonstrates that the biological control on pCO(2) (up to 100 mu atm), is 4 times larger than the thermal component and driven by annual net community production of 2.45 +/- 1.47 mol C m(-2) yr(-1). After the summer biological pCO(2) depletion, the return to near-atmospheric equilibrium proceeds slowly, driven in part by autumn entrainment into a deepening mixed layer and achieving full equilibration in late winter and early spring as respiration and advection complete the annual cycle. The shutdown of winter convection and associated mixed layer shoaling proceeds intermittently, appearing to frustrate the initiation of production. Horizontal processes, identified from salinity anomalies, are associated with biological pCO(2) signatures but with differing impacts in winter (when they reflect far-field variations in dissolved inorganic carbon and/or biomass) and summer (when they suggest promotion of local production by the relief of silicic acid or iron limitation). These results provide clarity on SAZ seasonal carbon cycling and demonstrate that the magnitude of the seasonal pCO(2) cycle is twice as large as that in the subarctic high-nutrient, low-chlorophyll waters, which can inform the selection of optimal global models in this region

Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime

Electron self-injection and acceleration until dephasing in the blowout regime is studied for a set of initial conditions typical of recent experiments with 100 terawatt-class lasers. Two different approaches to computationally efficient, fully explicit, three-dimensional particle-in-cell modelling are examined. First, the Cartesian code VORPAL using a perfect-dispersion electromagnetic solver precisely describes the laser pulse and bubble dynamics, taking advantage of coarser resolution in the propagation direction, with a proportionally larger time step. Using third-order splines for macroparticles helps suppress the sampling noise while keeping the usage of computational resources modest. The second way to reduce the simulation load is using reduced-geometry codes. In our case, the quasi-cylindrical code CALDER-CIRC uses decomposition of fields and currents into a set of poloidal modes, while the macroparticles move in the Cartesian 3D space. Cylindrical symmetry of the interaction allows using just two modes, reducing the computational load to roughly that of a planar Cartesian simulation while preserving the 3D nature of the interaction. This significant economy of resources allows using fine resolution in the direction of propagation and a small time step, making numerical dispersion vanishingly small, together with a large number of particles per cell, enabling good particle statistics. Quantitative agreement of the two simulations indicates that they are free of numerical artefacts. Both approaches thus retrieve physically correct evolution of the plasma bubble, recovering the intrinsic connection of electron self-injection to the nonlinear optical evolution of the driver

An Exactly Conservative Integrator for the n-Body Problem

The two-dimensional n-body problem of classical mechanics is a non-integrable Hamiltonian system for n > 2. Traditional numerical integration algorithms, which are polynomials in the time step, typically lead to systematic drifts in the computed value of the total energy and angular momentum. Even symplectic integration schemes exactly conserve only an approximate Hamiltonian. We present an algorithm that conserves the true Hamiltonian and the total angular momentum to machine precision. It is derived by applying conventional discretizations in a new space obtained by transformation of the dependent variables. We develop the method first for the restricted circular three-body problem, then for the general two-dimensional three-body problem, and finally for the planar n-body problem. Jacobi coordinates are used to reduce the two-dimensional n-body problem to an (n-1)-body problem that incorporates the constant linear momentum and center of mass constraints. For a four-body choreography, we find that a larger time step can be used with our conservative algorithm than with symplectic and conventional integrators.Comment: 17 pages, 3 figures; to appear in J. Phys. A.: Math. Ge

Mechanisms Driving Decadal Changes in the Carbonate System of a Coastal Plain Estuary

Understanding decadal changes in the coastal carbonate system is essential for predicting how the health of these waters responds to anthropogenic drivers, such as changing atmospheric conditions and riverine inputs. However, studies that quantify the relative impacts of these drivers are lacking. In this study, the primary drivers of decadal trends in the surface carbonate system, and the spatiotemporal variability in these trends, are identified for a large coastal plain estuary: the Chesapeake Bay. Experiments using a coupled three-dimensional hydrodynamic-biogeochemical model highlight that, over the past three decades, the changes in the surface carbonate system of Chesapeake Bay have strong seasonal and spatial variability. The greatest surface pH and aragonite saturation state (ΩAR) reductions have occurred in the summer in the middle (mesohaline) Bay: −0.24 and −0.9 per 30 years, respectively, with increases in atmospheric CO2 and reductions in nitrate loading both being primary drivers. Reductions in nitrate loading have a strong seasonal influence on the carbonate system, with the most pronounced decadal decreases in pH and ΩAR occurring during the summer when primary production is strongly dependent on nutrient availability. Increases in riverine total alkalinity and dissolved inorganic carbon have raised surface pH in the upper oligohaline Bay, while other drivers such as atmospheric warming and input of acidified ocean water through the Bay mouth have had comparatively minor impacts on the estuarine carbonate system. This work has significant implications for estuarine ecosystem services, which are typically most sensitive to surface acidification in the spring and summer seasons

Effective Hamiltonian approach to adiabatic approximation in open systems

The adiabatic approximation in open systems is formulated through the effective Hamiltonian approach. By introducing an ancilla, we embed the open system dynamics into a non-Hermitian quantum dynamics of a composite system, the adiabatic evolution of the open system is then defined as the adiabatic dynamics of the composite system. Validity and invalidity conditions for this approximation are established and discussed. A High-order adiabatic approximation for open systems is introduced. As an example, the adiabatic condition for an open spin-$\frac 1 2$ particle in time-dependent magnetic fields is analyzed.Comment: 6 pages, 2 figure

Efficient electron injection into plasma waves using higher-orderlaser modes

Using higher-order transverse laser modes as drivers forplasma wave excitation, and, in particular, using a ring laser beam withmaximum intensity off-axis, results in reversal of the focusinganddefocusing phase regions in a laser wakefield accelerator. Thisresults in improved performance of self-trapping and laser injectionschemes. Specifically, the trapping threshold required foropticalinjection is decreased and the maximum energy gain of the trappedelectrons is increased. This scheme could also be of interest for thegeneration of ring electron beams or for beam conditioning

Measles among migrants in the European Union and the European Economic Area

Aims: Progress towards meeting the goal of measles elimination in the EU and the European Economic Area (EEA) by 2015 is being obstructed, as some children are either not immunized on time or never immunized. One group thought to be at increased risk of measles is migrants; however, the extent to which this is the case is poorly understood, due to a lack of data. This paper addresses this evidence gap by providing an overview of the burden of measles in migrant populations in the EU/EEA. Methods: Data were collected through a comprehensive literature review, a country survey of EU/EEA member states and information from measles experts gathered at an infectious disease workshop. Results: Our results showed incomplete data on measles in migrant populations, as national surveillance systems do not systematically record migration-specific information; however, evidence from the literature review and country survey suggested that some measles outbreaks in the EU/EEA were due to sub-optimal vaccination coverage in migrant populations. Conclusions: We conclude that it is essential that routine surveillance of measles cases and measles, mumps and rubella (MMR) vaccination coverage become strengthened, to capture migrant-specific data. These data can help to inform the provision of preventive services, which may need to reach out to vulnerable migrant populations that currently face barriers in accessing routine immunization and health services