Angle of Repose and Angle of Marginal Stability: Molecular Dyanmics of Granular Particles

We present an implementation of realistic static friction in molecular dynamics (MD) simulations of granular particles. In our model, to break contacts between two particles, one has to apply a finite amount of force, determined by the Coulomb criterion. Using a two dimensional model, we show that piles generated by avalanches have a {\it finite} angle of repose $\theta_R$ (finite slopes). Furthermore, these piles are stable under tilting by an angle smaller than a non-zero tilting angle $\theta_T$, showing that $\theta_R$ is different from the angle of marginal stability $\theta_{MS}$, which is the maximum angle of stable piles. These measured angles are compared to a theoretical approximation. We also measure $\theta_{MS}$ by continuously adding particles on the top of a stable pile.Comment: 14 pages, Plain Te

Algorithm for numerical integration of the rigid-body equations of motion

A new algorithm for numerical integration of the rigid-body equations of motion is proposed. The algorithm uses the leapfrog scheme and the quantities involved are angular velocities and orientational variables which can be expressed in terms of either principal axes or quaternions. Due to specific features of the algorithm, orthonormality and unit norms of the orientational variables are integrals of motion, despite an approximate character of the produced trajectories. It is shown that the method presented appears to be the most efficient among all known algorithms of such a kind.Comment: 4 pages, 1 figur

Striation and convection in penumbral filaments

Observations with the 1-m Swedish Solar Telescope of the flows seen in penumbral filaments are presented. Time sequences of bright filaments show overturning motions strikingly similar to those seen along the walls of small isolated structures in the active regions. The filaments show outward propagating striations with inclination angles suggesting that they are aligned with the local magnetic field. We interpret it as the equivalent of the striations seen in the walls of small isolated magnetic structures. Their origin is then a corrugation of the boundary between an overturning convective flow inside the filament and the magnetic field wrapping around it. The outward propagation is a combination of a pattern motion due to the downflow observed along the sides of bright filaments, and the Evershed flow. The observed short wavelength of the striation argues against the existence of a dynamically significant horizontal field inside the bright filaments. Its intensity contrast is explained by the same physical effect that causes the dark cores of filaments, light bridges and `canals'. In this way striation represents an important clue to the physics of penumbral structure and its relation with other magnetic structures on the solar surface. We put this in perspective with results from the recent 3-D radiative hydrodynamic simulations.Comment: Accepted for publication in A&

Fine structure, magnetic field and heating of sunspot penumbrae

We interpret penumbral filaments as due to convection in field-free, radially aligned gaps just below the visible surface of the penumbra, intruding into a nearly potential field above. This solves the classical discrepancy between the large heat flux and the low vertical velocities observed in the penumbra. The presence of the gaps causes strong small-scale fluctuations in inclination, azimuth angle and field strength, but without strong forces acting on the gas. The field is nearly horizontal in a region around the cusp-shaped top of the gap, thereby providing an environment for Evershed flows. We identify this region with the recently discovered dark penumbral cores. Its darkness has the same cause as the dark lanes in umbral light-bridges, reproduced in numerical simulations by Nordlund and Stein (2005). We predict that the large vertical and horizontal gradients of the magnetic field inclination and azimuth in the potential field model will produce the net circular polarization seen in observations. The model also explains the significant elevation of bright filaments above their surroundings. It predicts that dark areas in the penumbra are of two different kinds: dark filament cores containing the most inclined (horizontal) fields, and regions between bright filaments, containing the least inclined field lines.Comment: submitted to A&

Yukawa potentials in systems with partial periodic boundary conditions I : Ewald sums for quasi-two dimensional systems

Yukawa potentials are often used as effective potentials for systems as colloids, plasmas, etc. When the Debye screening length is large, the Yukawa potential tends to the non-screened Coulomb potential ; in this small screening limit, or Coulomb limit, the potential is long ranged. As it is well known in computer simulation, a simple truncation of the long ranged potential and the minimum image convention are insufficient to obtain accurate numerical data on systems. The Ewald method for bulk systems, i.e. with periodic boundary conditions in all three directions of the space, has already been derived for Yukawa potential [cf. Y., Rosenfeld, {\it Mol. Phys.}, \bm{88}, 1357, (1996) and G., Salin and J.-M., Caillol, {\it J. Chem. Phys.}, \bm{113}, 10459, (2000)], but for systems with partial periodic boundary conditions, the Ewald sums have only recently been obtained [M., Mazars, {\it J. Chem. Phys.}, {\bf 126}, 056101 (2007)]. In this paper, we provide a closed derivation of the Ewald sums for Yukawa potentials in systems with periodic boundary conditions in only two directions and for any value of the Debye length. A special attention is paid to the Coulomb limit and its relation with the electroneutrality of systems.Comment: 40 pages, 5 figures and 4 table

Vulnerability of the British swine industry to classical swine fever

Classical swine fever (CSF) is a notifiable, highly contagious viral disease of swine which results in severe welfare and economic consequences in affected countries. To improve preparedness, it is critical to have some understanding of how CSF would spread should it be introduced. Based on the data recorded during the 2000 epidemic of CSF in Great Britain (GB), a spatially explicit, premises-based model was developed to explore the risk of CSF spread in GB. We found that large outbreaks of CSF would be rare and generated from a limited number of areas in GB. Despite the consistently low vulnerability of the British swine industry to large CSF outbreaks, we identified concerns with respect to the role played by the non-commercial sector of the industry. The model further revealed how various epidemiological features may influence the spread of CSF in GB, highlighting the importance of between-farm biosecurity in preventing widespread dissemination of the virus. Knowledge of factors affecting the risk of spread are key components for surveillance planning and resource allocation, and this work provides a valuable stepping stone in guiding policy on CSF surveillance and control in GB

Adaptive Management and the Value of Information: Learning Via Intervention in Epidemiology

Optimal intervention for disease outbreaks is often impeded by severe scientific uncertainty. Adaptive management (AM), long-used in natural resource management, is a structured decision-making approach to solving dynamic problems that accounts for the value of resolving uncertainty via real-time evaluation of alternative models. We propose an AM approach to design and evaluate intervention strategies in epidemiology, using real-time surveillance to resolve model uncertainty as management proceeds, with foot-and-mouth disease (FMD) culling and measles vaccination as case studies. We use simulations of alternative intervention strategies under competing models to quantify the effect of model uncertainty on decision making, in terms of the value of information, and quantify the benefit of adaptive versus static intervention strategies. Culling decisions during the 2001 UK FMD outbreak were contentious due to uncertainty about the spatial scale of transmission. The expected benefit of resolving this uncertainty prior to a new outbreak on a UK-like landscape would be £45–£60 million relative to the strategy that minimizes livestock losses averaged over alternate transmission models. AM during the outbreak would be expected to recover up to £20.1 million of this expected benefit. AM would also recommend a more conservative initial approach (culling of infected premises and dangerous contact farms) than would a fixed strategy (which would additionally require culling of contiguous premises). For optimal targeting of measles vaccination, based on an outbreak in Malawi in 2010, AM allows better distribution of resources across the affected region; its utility depends on uncertainty about both the at-risk population and logistical capacity. When daily vaccination rates are highly constrained, the optimal initial strategy is to conduct a small, quick campaign; a reduction in expected burden of approximately 10,000 cases could result if campaign targets can be updated on the basis of the true susceptible population. Formal incorporation of a policy to update future management actions in response to information gained in the course of an outbreak can change the optimal initial response and result in significant cost savings. AM provides a framework for using multiple models to facilitate public-health decision making and an objective basis for updating management actions in response to improved scientific understanding

Could changes in the agricultural landscape of northeastern China have influenced the long-distance transmission of highly pathogenic avian influenza H5Nx viruses?

In the last few years, several reassortant subtypes of highly pathogenic avian influenza viruses (HPAI H5Nx) have emerged in East Asia. These new viruses, mostly of subtype H5N1, H5N2, H5N6, and H5N8 belonging to clade 2.3.4.4, have been found in several Asian countries and have caused outbreaks in poultry in China, South Korea, and Vietnam. HPAI H5Nx also have spread over considerable distances with the introduction of viruses belonging to the same 2.3.4.4 clade in the U.S. (2014-2015) and in Europe (2014-2015 and 2016-2017). In this paper, we examine the emergence and spread of these new viruses in Asia in relation to published datasets on HPAI H5Nx distribution, movement of migratory waterfowl, avian influenza risk models, and land-use change analyses. More specifically, we show that between 2000 and 2015, vast areas of northeast China have been newly planted with rice paddy fields (3.21 million ha in Heilongjiang, Jilin, and Liaoning) in areas connected to other parts of Asia through migratory pathways of wild waterfowl. We hypothesize that recent land use changes in northeast China have affected the spatial distribution of wild waterfowl, their stopover areas, and the wild-domestic interface, thereby altering transmission dynamics of avian influenza viruses across flyways. Detailed studies of the habitat use by wild migratory birds, of the extent of the wild-domestic interface, and of the circulation of avian influenza viruses in those new planted areas may help to shed more light on this hypothesis, and on the possible impact of those changes on the long-distance patterns of avian influenza transmission

Molecular-dynamics simulations of the dynamical excitations in commensurate submonolayer films of nitrogen molecules on graphite

URL:http://link.aps.org/doi/10.1103/PhysRevB.54.14077 DOI:10.1103/PhysRevB.54.14077The dynamics of commensurate submonolayer solids of N2 molecules adsorbed on the basal planes of graphite have been studied using molecular-dynamics simulations. The calculations yielded the temperature dependence of the Brillouin-zone-center gap in the acoustic-phonon branches, for comparison with inelastic neutron-scattering experiments on the submonolayer solid. The calculated frequency gap was the same in submonolayer and monolayer films at low temperatures. At intermediate temperatures, the diffusive molecular motion associated with the presence of vacancies caused the gap mode to be less clearly defined in the coherent scattering function. Diffusion constants are calculated at submonolayer coverages, and temperatures up to 40 K for a population of molecules identified as mobile.This work was partially supported by the National Science Foundation under Grant No. DMR-9314235 (H.T.) and Nos. DMR-9120199 and DMR-9423307 (L.W.B.) and by The Danish Natural Science Foundation (F.Y.H.). L.W.B. thanks the Fysisk-Kemisk Institut and the Technical University of Denmark for hospitality during the period this work was completed

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.Comment: 73 pages, 19 figures, accepted by Solar Physic