99 research outputs found

    Steered Transition Path Sampling

    Full text link
    We introduce a path sampling method for obtaining statistical properties of an arbitrary stochastic dynamics. The method works by decomposing a trajectory in time, estimating the probability of satisfying a progress constraint, modifying the dynamics based on that probability, and then reweighting to calculate averages. Because the progress constraint can be formulated in terms of occurrences of events within time intervals, the method is particularly well suited for controlling the sampling of currents of dynamic events. We demonstrate the method for calculating transition probabilities in barrier crossing problems and survival probabilities in strongly diffusive systems with absorbing states, which are difficult to treat by shooting. We discuss the relation of the algorithm to other methods.Comment: 11 pages, 8 figure

    Regional diversity in subsistence among early farmers in Southeast Europe revealed by archaeological organic residues

    Get PDF
    The spread of early farming across Europe from its origins in Southwest Asia was a culturally transformative process which took place over millennia. Within regions, the pace of the transition was probably related to the particular climatic and environmental conditions encountered, as well as the nature of localized hunter-gatherer and farmer interactions. The establishment of farming in the interior of the Balkans represents the first movement of Southwest Asian livestock beyond their natural climatic range, and widespread evidence now exists for early pottery being used extensively for dairying. However, pottery lipid residues from sites in the Iron Gates region of the Danube in the northern Balkans show that here, Neolithic pottery was being used predominantly for processing aquatic resources. This stands out not only within the surrounding region but also contrasts markedly with Neolithic pottery use across wider Europe. These findings provide evidence for the strategic diversity within the wider cultural and economic practices during the Neolithic, with this exceptional environmental and cultural setting offering alternative opportunities despite the dominance of farming in the wider region

    Practical recipes for the model order reduction, dynamical simulation, and compressive sampling of large-scale open quantum systems

    Full text link
    This article presents numerical recipes for simulating high-temperature and non-equilibrium quantum spin systems that are continuously measured and controlled. The notion of a spin system is broadly conceived, in order to encompass macroscopic test masses as the limiting case of large-j spins. The simulation technique has three stages: first the deliberate introduction of noise into the simulation, then the conversion of that noise into an equivalent continuous measurement and control process, and finally, projection of the trajectory onto a state-space manifold having reduced dimensionality and possessing a Kahler potential of multi-linear form. The resulting simulation formalism is used to construct a positive P-representation for the thermal density matrix. Single-spin detection by magnetic resonance force microscopy (MRFM) is simulated, and the data statistics are shown to be those of a random telegraph signal with additive white noise. Larger-scale spin-dust models are simulated, having no spatial symmetry and no spatial ordering; the high-fidelity projection of numerically computed quantum trajectories onto low-dimensionality Kahler state-space manifolds is demonstrated. The reconstruction of quantum trajectories from sparse random projections is demonstrated, the onset of Donoho-Stodden breakdown at the Candes-Tao sparsity limit is observed, a deterministic construction for sampling matrices is given, and methods for quantum state optimization by Dantzig selection are given.Comment: 104 pages, 13 figures, 2 table

    MIMO antenna design using characteristic mode concepts

    No full text
    In this thesis we show how characteristic mode (CM) theory can be used to in the design of the antennas with the low envelope correlation coefficients preferred for MIMO antennas. This is due to the fact that the radiated fields of CMs are orthogonal over the radiation sphere. We show how a handheld device structure can be excited at various points (which form the multi-antenna ports) to ensure that each port excites a set of CMs such that the mode indices in each set are different for each antenna. It is this important result which causes the natural orthogonality properties of the CMs to be transferred to the overall radiation patterns of the multiple antennas. The CM analysis we have performed includes the effects of structure losses; it appears to be the first time that the computation of such modes have been reported. A logical and satisfying methodology for MIMO antenna design is the result. The methodology is extended to include arbitrary MIMO scattering environments using polarization sensitive generalized characteristic modes, which is the first of its kind from both the theory and methodology perspective. This design approach for MIMO antennas is not restricted to handheld devices, and can be used in a variety of applications. Lastly we discuss the various port mode theories that yield far-field orthogonality for perfect scattering environments, the application of which has already been discussed in the literature. We show succinctly that physical port symmetry yields frequency insensitive modes, but very minor symmetry breaking yields very narrow low envelope correlation bandwidth, something that is not discussed in the literature, but discussed in detail in this thesis

    Petascale Parallelization of the Gyrokinetic Toroidal Code

    No full text
    The Gyrokinetic Toroidal Code (GTC) is a global, three-dimensional particle-in-cell application developed to study microturbulence in tokamak fusion devices. The global capability of GTC is unique, allowing researchers to systematically analyze important dynamics such as turbulence spreading. In this work we examine a new radial domain decomposition approach to allow scalability onto the latest generation of petascale systems. Extensive performance evaluation is conducted on three high performance computing systems: the IBM BG/P, the Cray XT4, and an Intel Xeon Cluster. Overall results show that the radial decomposition approach dramatically increases scalability, while reducing the memory footprint - allowing for fusion device simulations at an unprecedented scale. After a decade where high-end computing (HEC) was dominated by the rapid pace of improvements to processor frequencies, the performance of next-generation supercomputers is increasingly differentiated by varying interconnect designs and levels of integration. Understanding the tradeoffs of these system designs is a key step towards making effective petascale computing a reality. In this work, we examine a new parallelization scheme for the Gyrokinetic Toroidal Code (GTC) [?] micro-turbulence fusion application. Extensive scalability results and analysis are presented on three HEC systems: the IBM BlueGene/P (BG/P) at Argonne National Laboratory, the Cray XT4 at Lawrence Berkeley National Laboratory, and an Intel Xeon cluster at Lawrence Livermore National Laboratory. Overall results indicate that the new radial decomposition approach successfully attains unprecedented scalability to 131,072 BG/P cores by overcoming the memory limitations of the previous approach. The new version is well suited to utilize emerging petascale resources to access new regimes of physical phenomena

    Pioneer farming in southeast Europe during the early sixth millennium BC: Climate-related adaptations in the exploitation of plants and animals (vol 13, e0197225, 2018)

    No full text
    [This corrects the article DOI: 10.1371/journal.pone.0197225.].status: publishe
    corecore