Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

The generalized Langevin equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general nonequilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e. g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly nontrivial, especially if the memory kernels decay nonexponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this paper, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.

Predictors and dynamics of postpartum relapses in women with multiple sclerosis

Several studies have shown that pregnancy reduces multiple sclerosis (MS) relapses, which increase in the early postpartum period. Postpartum relapse risk has been predicted by pre-pregnancy disease activity in some studies.

Dynamics of the entanglement spectrum in spin chains

We study the dynamics of the entanglement spectrum, that is the time evolution of the eigenvalues of the reduced density matrices after a bipartition of a one-dimensional spin chain. Starting from the ground state of an initial Hamiltonian, the state of the system is evolved in time with a new Hamiltonian. We consider both instantaneous and quasi adiabatic quenches of the system Hamiltonian across a quantum phase transition. We analyse the Ising model that can be exactly solved and the XXZ for which we employ the time-dependent density matrix renormalisation group algorithm. Our results show once more a connection between the Schmidt gap, i.e. the difference of the two largest eigenvalues of the reduced density matrix and order parameters, in this case the spontaneous magnetisation.

The Enemy of My Enemy is My Friend... The Dynamics of Self-Defense Forces in Irregular War: The Case of the Sons of Iraq

This article assesses the effect that leveraging civilian defense force militias has on the dynamics of violence in civil war. We argue that the delegation of security and combat roles to local civilians shifts the primary targets of insurgent violence toward civilians, in an attempt to deter future defections, and re-establish control over the local population. This argument is assessed through an analysis of the Sunni Awakening and ancillary Sons of Iraq paramilitary program. The results suggest that at least in the Al-Anbar province of Iraq, the utilization of the civilian population in counterinsurgent roles had significant implications for the targets of insurgent violence.

Dynamics of interacting Dicke model in a coupled-cavity array

We consider the dynamics of an array of mutually interacting cavities, each containing an ensemble of N two-level atoms. By exploring the possibilities offered by ensembles of various dimensions and a range of atom-light and photon-hopping values, we investigate the generation of multisite entanglement, as well as the performance of excitation transfer across the array, resulting from the competition between on-site nonlinearities of the matter-light interaction and intersite photon hopping. In particular, for a three-cavity interacting system it is observed that the initial excitation in the first cavity completely transfers to the ensemble in the third cavity through the hopping of photons between the adjacent cavities. Probabilities of the transfer of excitation of the cavity modes and ensembles exhibit characteristics of fast and slow oscillations governed by coupling and hopping parameters, respectively. In the large-hopping case, by seeding an initial excitation in the cavity at the center of the array, a tripartite W state, as well as a bipartite maximally entangled state, is obtained, depending on the interaction time. Population of the ensemble in a cavity has a positive impact on the rate of excitation transfer between the ensembles and their local cavity modes. In particular, for ensembles of five to seven atoms, tripartite W states can be produced even when the hopping rate is comparable to the cavity-atom coupling rate. A similar behavior of the transfer of excitation is observed for a four-coupled-cavity system with two initial excitations.

Investigating the dynamics of laser induced sparks in atmospheric helium using Rayleigh and Thomson scattering

We have used optical Rayleigh and Thomson scattering to investigate the expansion dynamics of laser induced plasma in atmospheric helium and to map its electron parameters both in time and space. The plasma is created using 9 ns duration, 140 mJ pulses from a Nd:YAG laser operating at 1064 nm, focused with a 10 cm focal length lens, and probed with 7 ns, 80 mJ, and 532 nm Nd:YAG laser pulses. Between 0.4 μs and 22.5 μs after breakdown, the electron density decreases from 3.3 × 1017 cm−3 to 9 × 1013 cm−3, while the temperature drops from 3.2 eV to 0.1 eV. Spatially resolved Thomson scattering data recorded up to 17.5 μs reveal that during this time the laser induced plasma expands at a rate given by R ∼ t0.4 consistent with a non-radiative spherical blast wave. This data also indicate the development of a toroidal structure in the lateral profile of both electron temperature and density. Rayleigh scattering data show that the gas density decreases in the center of the expanding plasma with a central scattering peak reemerging after about 12 μs. We have utilized a zero dimensional kinetic global model to identify the dominant particle species versus delay time and this indicates that metastable helium and the He2 + molecular ion play an important role.

The Dynamics of Rapid Redshifted and Blueshifted Excursions in the Solar Hα Line

We analyze high temporal and spatial resolution time-series of spectralscans of the Hα line obtained with the CRisp Imaging SpectroPolarimeter instrument mounted on the Swedish Solar Telescope.The data reveal highly dynamic, dark, short-lived structures known asRapid Redshifted and Blueshifted Excursions (RREs, RBEs) that areon-disk absorption features observed in the red and blue wings ofspectral lines formed in the chromosphere. We study the dynamics of RREsand RBEs by tracking their evolution in space and time, measuring thespeed of the apparent motion, line of sight (LOS) Doppler velocity, andtransverse velocity of individual structures. A statistical study oftheir measured properties shows that RREs and RBEs have similaroccurrence rates, lifetimes, lengths, and widths. They also displaynon-periodic, nonlinear transverse motions perpendicular to their axesat speeds of 4-31 km s-1. Furthermore, both typesof structures either appear as high speed jets and blobs that aredirected outwardly from a magnetic bright point with speeds of50-150 km s-1, or emerge within a few seconds. Astudy of the different velocity components suggests that the transversemotions along the LOS of the chromospheric flux tubes are responsiblefor the formation and appearance of these redshifted/blueshiftedstructures. The short lifetime and fast disappearance of the RREs/RBEssuggests that, similar to type II spicules, they are rapidly heated totransition region or even coronal temperatures. We speculate that theKelvin-Helmholtz instability triggered by observed transversemotions of these structures may be a viable mechanism for their heating.

Ecological dynamics of extinct species in empty habitat networks. 1. The role of habitat pattern and quantity, stochasticity and dispersal

We examined a remnant host plant (Primula veris L.) habitat network that was last inhabited by the rare butterfly Hamearis lucina L. in north Wales in 1943, to assess the relative contribution of several spatial parameters to its regional extinction. We first examined relationships between P. veris characteristics and H. lucina eggs in surviving H. lucina populations, and used these to predict the suitability and potential carrying capacity of the habitat network in north Wales. This resulted in an estimate of roughly 4500 eggs (ca 227 adults). We developed a discrete space, discrete time metapopulation model to evaluate the relative contribution of dispersal distance, habitat and environmental stochasticity as possible causes of extinction. We simulated the potential persistence of the butterfly in the current network as well as in three artificial (historical and present) habitat networks that differed in the quantity (current and X3) and fragmentation of the habitat (current and aggregated). We identified that reduced habitat quantity and increased isolation would have increased the probability of regional extinction, in conjunction with environmental stochasticity and H. lucina's dispersal distance. This general trend did not change in a qualitative manner when we modified the ability of dispersing females to stay in, and find suitable habitats (by changing the size of the grid cells used in the model). Contrary to most metapopulation model predictions, system persistence declined with increasing migration rate, suggesting that the mortality of migrating individuals in fragmented landscapes may pose significant risks to system-wide persistence. Based on model predictions for the present landscape we argue that a major programme of habitat restoration would be required for a re-established metapopulation to persist for > 100 years.

Ecological dynamics of extinct species in empty habitat networks. 2. The role of host plant dynamics

This paper explores the relative effects of host plant dynamics and butterfly-related parameters on butterfly persistence. It considers an empty habitat network where a rare butterfly (Cupido minimus) became extinct in 1939 in part of its historical range in north Wales, UK. Surviving populations of the butterfly in southern Britain were visited to assess use of its host plant (Anthyllis vulneraria) in order to calibrate habitat suitability and carrying capacity in the empty network in north Wales. These data were used to deduce that only a portion ( similar to 19%) of the host plant network from north Wales was likely to be highly suitable for oviposition. Nonetheless, roughly 65,460 eggs (3273 adult equivalents) could be expected to be laid in north Wales, were the empty network to be populated at the same levels as observed on comparable plants in surviving populations elsewhere. Simulated metapopulations of C. minimus in the empty network revealed that time to extinction and patch occupancy were significantly influenced by carrying capacity, butterfly mean dispersal distance and environmental stochasticity, although for most reasonable parameter values, the model system persisted. Simulation outputs differed greatly when host plant dynamics was incorporated into the modelled butterfly dynamics. Cupido minimus usually went extinct when host plant were at low densities. In these simulations host plant dynamics appeared to be the most important determinant of the butterfly's regional extirpation. Modelling the outcome of a reintroduction programme to C. minimus variation at high quality locations, revealed that 65% of systems survived at least 100 years. Given the current amount of resources of the north Wales landscape, the persistence of C. minimus under a realistic reintroduction programme has a good chance of being successful, if carried out in conjunction with a host plant management programme.