skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Study of a mixed dispersal population dynamics model

Abstract

In this study, we consider a mixed dispersal model with periodic and Dirichlet boundary conditions and its corresponding linear eigenvalue problem. This model describes the time evolution of a population which disperses both locally and non-locally. We investigate how long time dynamics depend on the parameter values. Furthermore, we study the minimization of the principal eigenvalue under the constraints that the resource function is bounded from above and below, and with a fixed total integral. Biologically, this minimization problem is motivated by the question of determining the optimal spatial arrangement of favorable and unfavorable regions for the species to die out more slowly or survive more easily. Our numerical simulations indicate that the optimal favorable region tends to be a simply-connected domain. Numerous results are shown to demonstrate various scenarios of optimal favorable regions for periodic and Dirichlet boundary conditions.

Authors:
 [1];  [2];  [3];  [4];  [5];  [6]
  1. Claremont Graduate Univ., Claremont, CA (United States). Dept. of Mathematics
  2. Rochester Institute of Technology, Rochester, NY (United States). School of Mathematics Sciences
  3. Claremont McKenna College, Claremont, CA (United States). Dept. of Mathematical Sciences
  4. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Center for Applied Scientific Computing
  5. Univ. of Maryland Baltimore County (UMBC), Baltimore, MD (United States). Dept. of Mathematics and Statistics
  6. Brown Univ., Providence, RI (United States). Dept. of Applied Mathematics
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1305855
Report Number(s):
LLNL-JRNL-674461
Journal ID: ISSN 0940-6573
Grant/Contract Number:
AC52-07NA27344
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
The IMA Volumes in Mathematics and its Applications
Additional Journal Information:
Journal Volume: 160; Journal ID: ISSN 0940-6573
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE

Citation Formats

Chugunova, Marina, Jadamba, Baasansuren, Kao, Chiu -Yen, Klymko, Christine F., Thomas, Evelyn, and Zhao, Bingyu. Study of a mixed dispersal population dynamics model. United States: N. p., 2016. Web. doi:10.1007/978-1-4939-6399-7_3.
Chugunova, Marina, Jadamba, Baasansuren, Kao, Chiu -Yen, Klymko, Christine F., Thomas, Evelyn, & Zhao, Bingyu. Study of a mixed dispersal population dynamics model. United States. doi:10.1007/978-1-4939-6399-7_3.
Chugunova, Marina, Jadamba, Baasansuren, Kao, Chiu -Yen, Klymko, Christine F., Thomas, Evelyn, and Zhao, Bingyu. 2016. "Study of a mixed dispersal population dynamics model". United States. doi:10.1007/978-1-4939-6399-7_3. https://www.osti.gov/servlets/purl/1305855.
@article{osti_1305855,
title = {Study of a mixed dispersal population dynamics model},
author = {Chugunova, Marina and Jadamba, Baasansuren and Kao, Chiu -Yen and Klymko, Christine F. and Thomas, Evelyn and Zhao, Bingyu},
abstractNote = {In this study, we consider a mixed dispersal model with periodic and Dirichlet boundary conditions and its corresponding linear eigenvalue problem. This model describes the time evolution of a population which disperses both locally and non-locally. We investigate how long time dynamics depend on the parameter values. Furthermore, we study the minimization of the principal eigenvalue under the constraints that the resource function is bounded from above and below, and with a fixed total integral. Biologically, this minimization problem is motivated by the question of determining the optimal spatial arrangement of favorable and unfavorable regions for the species to die out more slowly or survive more easily. Our numerical simulations indicate that the optimal favorable region tends to be a simply-connected domain. Numerous results are shown to demonstrate various scenarios of optimal favorable regions for periodic and Dirichlet boundary conditions.},
doi = {10.1007/978-1-4939-6399-7_3},
journal = {The IMA Volumes in Mathematics and its Applications},
number = ,
volume = 160,
place = {United States},
year = 2016,
month = 8
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share:
  • We apply two approximate solutions of the quantum-classical Liouville equation (QCLE) in the mapping representation to the simulation of the laser-induced response of a quantum subsystem coupled to a classical environment. These solutions, known as the Poisson Bracket Mapping Equation (PBME) and the Forward-Backward (FB) trajectory solutions, involve simple algorithms in which the dynamics of both the quantum and classical degrees of freedom are described in terms of continuous variables, as opposed to standard surface-hopping solutions in which the classical degrees of freedom hop between potential energy surfaces dictated by the discrete adiabatic state of the quantum subsystem. The validitymore » of these QCLE-based solutions is tested on a non-trivial electron transfer model involving more than two quantum states, a time-dependent Hamiltonian, strong subsystem-bath coupling, and an initial energy shift between the donor and acceptor states that depends on the strength of the subsystem-bath coupling. In particular, we calculate the time-dependent population of the photoexcited donor state in response to an ultrafast, on-resonance pump pulse in a three-state model of an electron transfer complex that is coupled asymmetrically to a bath of harmonic oscillators through the optically dark acceptor state. Within this approach, the three-state electron transfer complex is treated quantum mechanically, while the bath oscillators are treated classically. When compared to the more accurate QCLE-based surface-hopping solution and to the numerically exact quantum results, we find that the PBME solution is not capable of qualitatively capturing the population dynamics, whereas the FB solution is. However, when the subsystem-bath coupling is decreased (which also decreases the initial energy shift between the donor and acceptor states) or the initial shift is removed altogether, both the PBME and FB results agree better with the QCLE-based surface-hopping results. These findings highlight the challenges posed by various conditions such as a time-dependent external field, the strength of the subsystem-bath coupling, and the degree of asymmetry on the accuracy of the PBME and FB algorithms.« less
  • In a previous study [F. A. Shakib and G. Hanna, J. Chem. Phys. 141, 044122 (2014)], we investigated a model proton-coupled electron transfer (PCET) reaction via the mixed quantum-classical Liouville (MQCL) approach and found that the trajectories spend the majority of their time on the mean of two coherently coupled adiabatic potential energy surfaces. This suggested a need for mean surface evolution to accurately simulate observables related to ultrafast PCET processes. In this study, we simulate the time-dependent populations of the three lowest adiabatic states in the ET-PT (i.e., electron transfer preceding proton transfer) version of the same PCET modelmore » via the MQCL approach and compare them to the exact quantum results and those obtained via the fewest switches surface hopping (FSSH) approach. We find that the MQCL population profiles are in good agreement with the exact quantum results and show a significant improvement over the FSSH results. All of the mean surfaces are shown to play a direct role in the dynamics of the state populations. Interestingly, our results indicate that the population transfer to the second-excited state can be mediated by dynamics on the mean of the ground and second-excited state surfaces, as part of a sequence of nonadiabatic transitions that bypasses the first-excited state surface altogether. This is made possible through nonadiabatic transitions between different mean surfaces, which is the manifestation of coherence transfer in MQCL dynamics. We also investigate the effect of the strength of the coupling between the proton/electron and the solvent coordinate on the state population dynamics. Drastic changes in the population dynamics are observed, which can be understood in terms of the changes in the potential energy surfaces and the nonadiabatic couplings. Finally, we investigate the state population dynamics in the PT-ET (i.e., proton transfer preceding electron transfer) and concerted versions of the model. The PT-ET results confirm the participation of all of the mean surfaces, albeit in different proportions compared to the ET-PT case, while the concerted results indicate that the mean of the ground- and first-excited state surfaces only plays a role, due to the large energy gaps between the ground- and second-excited state surfaces.« less
  • A plasma with a fairly high degree of ionization is shown to be formed at relatively low current densities that are typical for many experiments on measuring the atomic constants. In this case, the electron density appears to be sufficient for the rate of inelastic electron-atom collisions to exceed the rate of radiative transitions between highly excited states. A quasisteady population distribution is established in the manifold of states with n > 5, which makes it possible to use simple and evident relations for calculating the populations. 7 refs., 2 figs.
  • Using the Van Allen Probes data, we study the radiation belt seed population and it associated with the relativistic electron dynamics during 74 geomagnetic storm events. Based on the flux changes of 1 MeV electrons before and after the storm peak, these storm events are divided into two groups of “non-preconditioned” and “preconditioned”. The statistical study shows that the storm intensity is of significant importance for the distribution of the seed population (336 keV electrons) in the outer radiation belt. However, substorm intensity can also be important to the evolution of the seed population for some geomagnetic storm events. Formore » non-preconditioned storm events, the correlation between the peak fluxes and their L-shell locations of the seed population and relativistic electrons (592 keV, 1.0 MeV, 1.8 MeV, and 2.1 MeV) is consistent with the energy-dependent dynamic processes in the outer radiation belt. For preconditioned storm events, the correlation between the features of the seed population and relativistic electrons is not fully consistent with the energy-dependent processes. It is suggested that the good correlation between the radiation belt seed population and ≤1.0 MeV electrons contributes to the prediction of the evolution of ≤1.0 MeV electrons in the Earth’s outer radiation belt during periods of geomagnetic storms.« less