Source author record

Stephen Pankavich

Stephen Pankavich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
20topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

Global Sensitivity Analysis of Plasma Instabilities via Active Subspaces

Active subspace analysis is a useful computational tool to identify and exploit the most important linear combinations in the space of a model's input parameters. These directions depend inherently on a quantity of interest, which can be represented as a function from input parameters to model outputs. As the dynamics of many plasma models are driven by potentially uncertain parameter values, the utilization of active subspaces to perform global sensitivity analysis represents an important tool to understand how certain physical phenomena depend upon fluctuations in the values of these parameters. In the current paper, we construct and implement new computational methods to quantify the induced uncertainty within the growth rate generated by perturbations in a collisionless plasma modeled by the one-dimensional Vlasov-Poisson system near an unstable, spatially-homogeneous equilibrium in the linear regime.

preprint2021arXiv

Asymptotic Dynamics of Dispersive, Collisionless Plasmas

A multispecies, collisionless plasma is modeled by the Vlasov-Poisson system. Assuming that the electric field decays with sufficient rapidity as $t \to\infty$, we show that the velocity characteristics and spatial averages of the particle distributions converge as time grows large. Using these limits we establish the precise asymptotic profile of the electric field and its derivatives, as well as, the charge and current densities. Modified spatial characteristics are then shown to converge using the limiting electric field. Finally, we establish a modified $L^\infty$ scattering result for each particle distribution function, namely we show that they converge as $t \to \infty$ along the modified spatial characteristics. When the plasma is non-neutral, the estimates of these quantities are sharp, while in the neutral case they may imply faster rates of decay.

preprint2016arXiv

A Modified SEIR Model for the Spread of Ebola in Western Africa and Metrics for Resource Allocation

A modified, deterministic SEIR model is developed for the 2014 Ebola epidemic occurring in the West African nations of Guinea, Liberia, and Sierra Leone. The model describes the dynamical interaction of susceptible and infected populations, while accounting for the effects of hospitalization and the spread of disease through interactions with deceased, but infectious, individuals. Using data from the World Health Organization (WHO), parameters within the model are fit to recent estimates of infected and deceased cases from each nation. The model is then analyzed using these parameter values. Finally, several metrics are proposed to determine which of these nations is in greatest need of additional resources to combat the spread of infection. These include local and global sensitivity metrics of both the infected population and the basic reproduction number with respect to rates of hospitalization and proper burial.

preprint2016arXiv

Mathematical Analysis and Dynamic Active Subspaces for a Long term model of HIV

Recently, a long-term model of HIV infection dynamics was developed to describe the entire time course of the disease. It consists of a large system of ODEs with many parameters, and is expensive to simulate. In the current paper, this model is analyzed by determining all infection-free steady states and studying the local stability properties of the unique biologically-relevant equilibrium. Active subspace methods are then used to perform a global sensitivity analysis and study the dependence of an infected individual's T-cell count on the parameter space. Building on these results, a global-in-time approximation of the T-cell count is created by constructing dynamic active subspaces and reduced order models are generated, thereby allowing for inexpensive computation.

preprint2016arXiv

Principal Component Analysis: Resources for an Essential Application of Linear Algebra

Principal Component Analysis (PCA) is a highly useful topic within an introductory Linear Algebra course, especially since it can be used to incorporate a number of applied projects. This method represents an essential application and extension of the Spectral Theorem and is commonly used within a variety of fields, including statistics, neuroscience, and image compression. We present a synopsis of PCA and include a number of examples that can be used within upper-level mathematics courses to engage undergraduate students while introducing them to one of the most widely-used applications of linear algebra.

preprint2016arXiv

Separated Characteristics and Global Solvability for the one and one-half dimensional Vlasov Maxwell System

The motion of a collisionless plasma - a high-temperature, low-density, ionized gas - is described by the Vlasov-Maxwell (VM) system. These equations are considered in one space dimension and two momentum dimensions without the assumption of relativistic velocity corrections. The main results are bounds on the spatial and velocity supports of the particle distribution function and uniform estimates on derivatives of this function away from the critical velocity $| v_1 | = 1$. Additionally, for initial particle distributions that are even in the second velocity argument $v_2$, the global-in-time existence of solutions is shown.

preprint2015arXiv

A Short Proof of Increased Parabolic Regularity

We present a new, short proof of the increased regularity obtained by solutions to uniformly parabolic partial differential equations. Though this setting is fairly introductory, our new method of proof, which uses a priori estimates, can be extended to prove analogous results for problems with time-dependent coefficients, transport equations, and nonlinear equations even when other tools, such as semigroup methods or the use of explicit fundamental solutions, are unavailable.

preprint2015arXiv

An in-host model of HIV incorporating latent infection and viral mutation

We construct a seven-component model of the in-host dynamics of the Human Immunodeficiency Virus Type-1 (i.e, HIV) that accounts for latent infection and the propensity of viral mutation. A dynamical analysis is conducted and a theorem is presented which characterizes the long time behavior of the model. Finally, we study the effects of an antiretroviral drug and treatment implications.

preprint2015arXiv

Instability conditions for some periodic BGK waves in the Vlasov-Poisson system

A one-dimensional, collisionless plasma given by the Vlasov-Poisson system is considered and the stability properties of periodic steady state solutions known as Bernstein-Greene-Kruskal (BGK) waves are investigated. Sufficient conditions are determined under which BGK waves are linearly unstable under perturbations that share the same period as the equilibria. It is also shown that such solutions cannot support a monotonically decreasing particle distribution function.

preprint2015arXiv

Mathematical analysis of an in-host model of viral dynamics with spatial heterogeneity

We consider a spatially-heterogeneous generalization of a well-established model for the dynamics of the Human Immunodeficiency Virus-type 1 (HIV) within a susceptible host. The model consists of a nonlinear system of three coupled reaction-diffusion equations with parameters that may vary spatially. Upon formulating the model, we prove that it preserves the positivity of initial data and construct global-in-time solutions that are both bounded and smooth. Finally, additional results concerning the local and global asymptotic behavior of these solutions are also provided.

preprint2014arXiv

Nanosystem Self-Assembly Pathways Discovered via All-Atom Multiscale Analysis

We consider the self-assembly of composite structures from a group of nanocomponents, each consisting of particles within an $N$-atom system. Self-assembly pathways and rates for nanocomposites are derived via a multiscale analysis of the classical Liouville equation. From a reduced statistical framework, rigorous stochastic equations for population levels of beginning, intermediate, and final aggregates are also derived. It is shown that the definition of an assembly type is a self-consistency criterion that must strike a balance between precision and the need for population levels to be slowly varying relative to the time scale of atomic motion. The deductive multiscale approach is complemented by a qualitative notion of multicomponent association and the ensemble of exact atomic-level configurations consistent with them. In processes such as viral self-assembly from proteins and RNA or DNA, there are many possible intermediates, so that it is usually difficult to predict the most efficient assembly pathway. However, in the current study, rates of assembly of each possible intermediate can be predicted. This avoids the need, as in a phenomenological approach, for recalibration with each new application. The method accounts for the feedback across scales in space and time that is fundamental to nanosystem self-assembly. The theory has applications to bionanostructures, geomaterials, engineered composites, and nanocapsule therapeutic delivery systems.

preprint2014arXiv

Spatially homogeneous solutions of the Vlasov-Nordström-Fokker-Planck system

The Vlasov-Nordström-Fokker-Planck system describes the evolution of self-gravitating matter experiencing collisions with a fixed background of particles in the framework of a relativistic scalar theory of gravitation. We study the spatially-homogeneous system and prove global existence and uniqueness of solutions for the corresponding initial value problem in three momentum dimensions. Additionally, we study the long time asymptotic behavior of the system and prove that even in the absence of friction, solutions possess a non-trivial asymptotic profile. An exact formula for the long time limit of the particle density is derived in the ultra-relativistic case.

preprint2013arXiv

A One-dimensional kinetic model of plasma dynamics with a transport field

Motivated by the fundamental model of a collisionless plasma, the Vlasov-Maxwell (VM) system, we consider a related, nonlinear system of partial differential equations in one space and one momentum dimension. As little is known regarding the regularity properties of solutions to the non-relativistic version of the (VM) equations, we study a simplified system which also lacks relativistic velocity corrections and prove local-in-time existence and uniqueness of classical solutions to the Cauchy problem. For special choices of initial data, global-in-time existence of these solutions is also shown. Finally, we provide an estimate which, independent of the initial data, yields additional global-in-time regularity of the associated field.

preprint2013arXiv

Analysis and Simulation of the Three-Component Model of HIV Dynamics

Mathematical modeling of biological systems is crucial to effectively and efficiently developing treatments for medical conditions that plague humanity. Often, systems of ordinary differential equations are a traditional tool used to describe the spread of disease within the body. We consider the dynamics of the Human Immunodeficiency Virus (HIV) in vivo during the initial stage of infection. In particular, we examine the well-known three-component model and prove the existence, uniqueness, and boundedness of solutions. Furthermore, we prove that solutions remain biologically meaningful, i.e., are positivity preserving, and perform a thorough, local stability analysis for the equilibrium states of the system. Finally, we incorporate random coefficients within the model and obtain numerical results to predict the probability of infection given the transmission of the virus to a new individual.

preprint2013arXiv

Global Existence for the "One and one-half" dimensional relativistic Vlasov-Maxwell-Fokker-Planck system

In a recent paper Calogero and Alcantara derived a Lorentz-invariant Fokker-Planck equation, which corresponds to the evolution of a particle distribution associated with relativistic Brownian Motion. We study the "one and one-half" dimensional version of this problem with nonlinear electromagnetic interactions - the relativistic Vlasov-Maxwell-Fokker-Planck system - and obtain the first results concerning well-posedness of solutions. Specifically, we prove the global-in-time existence and uniqueness of classical solutions to the Cauchy problem and a gain in regularity of the distribution function in its momentum argument.

preprint2013arXiv

Hierarchical Multiscale Modeling of Macromolecules and their Assemblies

Soft materials (e.g., enveloped viruses, liposomes, membranes and supercooled liquids) simultaneously deform or display collective behaviors, while undergoing atomic scale vibrations and collisions. While the multiple space-time character of such systems often makes traditional molecular dynamics simulation impractical, a multiscale approach has been presented that allows for long-time simulation with atomic detail based on the co-evolution of slowly-varying order parameters (OPs) with the quasi-equilibrium probability density of atomic configurations. However, this approach breaks down when the structural change is extreme, or when nearest-neighbor connectivity of atoms is not maintained. In the current study, a self-consistent approach is presented wherein OPs and a reference structure co-evolve slowly to yield long-time simulation for dynamical soft-matter phenomena such as structural transitions and self assembly. The development begins with the Liouville equation for N classical atoms and an ansatz on the form of the associated N-atom probability density. Multiscale techniques are used to derive Langevin equations for the coupled OP configurational dynamics. The net result is a set of equations for the coupled stochastic dynamics of the OPs and centers of mass of the subsystems that constitute a soft material body. The theory is based on an all-atom methodology and an interatomic force field, and therefore enables calibration-free simulations of soft matter, such as macromolecular assemblies.

preprint2013arXiv

Nonlinear instability of solutions in parabolic and hyperbolic diffusion

We consider semilinear evolution equations of the form $a(t)\partial_{tt}u + b(t) \partial_t u + Lu = f(x,u)$ and $b(t) \partial_t u + Lu = f(x,u),$ with possibly unbounded $a(t)$ and possibly sign-changing damping coefficient $b(t)$, and determine precise conditions for which linear instability of the steady state solutions implies nonlinear instability. More specifically, we prove that linear instability with an eigenfunction of fixed sign gives rise to nonlinear instability by either exponential growth or finite-time blow-up. We then discuss a few examples to which our main theorem is immediately applicable, including evolution equations with supercritical and exponential nonlinearities.

preprint2013arXiv

The Effects of Latent Infection on the Dynamics of HIV

One way in which the human immunodeficiency virus (HIV-1) replicates within a host is by infecting activated CD4+ T-cells, which then produce additional copies of the virus. Even with the introduction of antiretroviral drug therapy, which has been very successful over the past decade, a large obstacle to the complete eradication of the virus is the presence of viral reservoirs in the form of latently infected CD4+ T-cells. We consider a model of HIV infection that describes T-cell and viral interactions, as well as, the production and activation of latently infected T-cells. Upon determining equilibrium states of the latent cell model, the local and global asymptotic behavior of solutions is examined, and the basic reproduction number of the system is computed to be strictly less than that of the corresponding three-component model, which omits the effects of latent infection. In particular, this implies that a wider variety of parameter values will lead to viral eradication as $t \to \infty$ due to the appearance of latent CD4+ T-cells. With this realization we discuss possible alternative notions for eradication and persistence of infection other than traditional tools. These results are further illustrated by a number of numerical simulations.

preprint2012arXiv

Multiscale Theory of Finite Size Bose Systems: Implications for Collective and Single-Particle Excitations

Boson droplets (i.e., dense assemblies of bosons at low temperature) are shown to mask a significant amount of single-particle behavior and to manifest collective, droplet-wide excitations. To investigate the balance between single-particle and collective behavior, solutions to the wave equation for a finite size Bose system are constructed in the limit where the ratio \varepsilon of the average nearest-neighbor boson distance to the size of the droplet or the wavelength of density disturbances is small. In this limit, the lowest order wave function varies smoothly across the system, i.e., is devoid of structure on the scale of the average nearest-neighbor distance. The amplitude of short range structure in the wave function is shown to vanish as a power of \varepsilon when the interatomic forces are relatively weak. However, there is residual short range structure that increases with the strength of interatomic forces. While the multiscale approach is applied to boson droplets, the methodology is applicable to any finite size bose system and is shown to be more direct than field theoretic methods. Conclusions for Helium-4 nanodroplets are drawn.

preprint2011arXiv

Multiscaling for Systems with a Broad Continuum of Characteristic Lengths and Times: Structural Transitions in Nanocomposites

The multiscale approach to N-body systems is generalized to address the broad continuum of long time and length scales associated with collective behaviors. A technique is developed based on the concept of an uncountable set of time variables and of order parameters (OPs) specifying major features of the system. We adopt this perspective as a natural extension of the commonly used discrete set of timescales and OPs which is practical when only a few, widely-separated scales exist. The existence of a gap in the spectrum of timescales for such a system (under quasiequilibrium conditions) is used to introduce a continuous scaling and perform a multiscale analysis of the Liouville equation. A functional-differential Smoluchowski equation is derived for the stochastic dynamics of the continuum of Fourier component order parameters. A continuum of spatially non-local Langevin equations for the OPs is also derived. The theory is demonstrated via the analysis of structural transitions in a composite material, as occurs for viral capsids and molecular circuits.

preprint2011arXiv

Time Decay for solutions to One-Dimensional Two-Component Plasma Equations

We represent three generations of students: Bob Glassey, Walter's student finishing at Brown in 1972, Jack Schaeffer, Bob's student finishing at Indiana University in 1983, and Steve Pankavich, Jack's student finishing at Carnegie Mellon in 2005. We have all thrived professionally from our association with Walter and are delighted to dedicate this note to him on the occasion of his 70th birthday. The problem we study concerns the asymptotic behavior of solutions to Vlasov equations, an area to which Walter has contributed greatly.

preprint2010arXiv

A Particle Method for a Collisionless Plasma with Infinite Mass

The one-dimensional Vlasov-Poisson system is considered and a particle method is developed to approximate solutions without compact support which tend to a fixed background of charge as $| x | \to \infty$. Such a system of equations can be used to model kinetic phenomena occurring in plasma physics, such as the solar wind. The particle method is constructed, implemented, and used to determine information regarding the time asymptotics of the electrostatic field.

preprint2010arXiv

Decay in Time for a One-Dimensional Two-Component Plasma

The motion of a collisionless plasma is described by the Vlasov-Poisson system, or in the presence of large velocities, the relativistic Vlasov-Poisson system. Both systems are considered in one space and one momentum dimension, with two species of oppositely charged particles. A new identity is derived for both systems and is used to study the behavior of solutions for large times.

preprint2010arXiv

Explicit Solutions of the One-dimensional Vlasov-Poisson System with Infinite Mass and Energy

A collisionless plasma is modeled by the Vlasov-Poisson system in one-dimension. A fixed background of positive charge, dependent only upon velocity, is assumed and the situation in which the mobile negative ions balance the positive charge as x tends to positive or negative infinity. Thus, the total positive charge and the total negative charge are infinite. In this paper, the charge density of the system is shown to be compactly supported. More importantly, both the electric field and the number density are determined explicitly for large values of x.

preprint2010arXiv

Global Existence and Increased Spatial Decay for the Radial Vlasov-Poisson System with Steady Spatial Asymptotics

A collisionless plasma is modeled by the Vlasov-Poisson system in three space dimensions. A fixed background of positive charge, which is independent of time and space, is assumed. The situation in which mobile negative ions balance the positive charge as x tends to infinity is considered. Hence the total positive charge, total negative charge, and total energy are all infinite. Smooth solutions with appropriate asymptotic behavior for large x, which were previously shown to exist locally in time, are continued globally for spherically symmetric data. This is done by showing that the charge density decays at least as fast as x^{-4}. Finally, an increased decay rate of x^{-6} is shown in the general case without the assumption of spherical symmetry.

preprint2010arXiv

Global Existence for the Vlasov-Poisson System with Steady Spatial Asymptotics

A collisionless plasma is modeled by the Vlasov-Poisson system in three space dimensions. A fixed background of positive charge - dependant upon only velocity - is assumed. The situation in which mobile negative ions balance the positive charge as x tends to infinity is considered. Thus, the total positive charge and the total negative charge are both infinite. Smooth solutions with appropriate asymptotic behavior for large x, which were previously shown to exist locally in time, are continued globally. This is done by showing that the charge density decays at least as fast as x^{-6}. This article also establishes spatial decay estimates for the electrostatic field and its derivatives.

preprint2010arXiv

Liquid-Crystal Transitions: A First Principles Multiscale Approach

A rigorous theory of liquid-crystal transitions is developed starting from the Liouville equation. The starting point is an all-atom description and a set of order parameter field variables that are shown to evolve slowly via Newton's equations. The separation of timescales between that of atomic collisions and the order parameter fields enables the derivation of rigorous equations for stochastic order parameter field dynamics. When the fields provide a measure of the spatial profile of the probability of molecular position, orientation, and internal structure, a theory of liquid-crystal transitions emerges. The theory uses the all-atom/continuum approach developed earlier to obtain a functional generalization of the Smoluchowski equation wherein key atomic details are embedded. The equivalent non-local Langevin equations are derived and computational aspects are discussed. The theory enables simulations that are much less computationally intensive than molecular dynamics and thus does not require oversimplification of the system's constituent components. The equations obtained do not include factors that require calibration and can thus be applicable to various phase transitions which overcomes the limitations of phenomenological field models. The relation of the theory to phenomenological descriptions of Nematic and Smectic phase transitions, and the possible existence of other types of transitions involving intermolecular structural parameters are discussed.

preprint2010arXiv

Local Existence for the One-dimensional Vlasov-Poisson System with Infinite Mass

A collisionless plasma is modeled by the Vlasov-Poisson system in one dimension. We consider the situation in which mobile negative ions balance a fixed background of positive charge, which is independent of space and time, as x tends to positive or negative infinity. Thus, the total positive charge and the total negative charge are both infinite. Smooth solutions with appropriate asymptotic behavior are shown to exist locally in time, and criteria for the continuation of these solutions are established.

preprint2010arXiv

Multiscaling for Classical Nanosystems: Derivation of Smoluchowski and Fokker-Planck Equations

Using multiscale analysis and methods of statistical physics, we show that a solution to the N-atom Liouville Equation can be decomposed via an expansion in terms of a smallness parameter epsilon, wherein the long scale time behavior depends upon a reduced probability density that is a function of slow-evolving order parameters. This reduced probability density is shown to satisfy the Smoluchowski equation up to order epsilon squared for a given range of initial conditions. Furthermore, under the additional assumption that the nanoparticle momentum evolves on a slow time scale, we show that this reduced probability density satisfies a Fokker-Planck equation up to the same order in epsilon. This approach applies to a broad range of problems in the nanosciences.

preprint2010arXiv

Self-Assembly of Nanocomponents into Composite Structures: Derivation and Simulation of Langevin Equations

The kinetics of the self-assembly of nanocomponents into a virus, nanocapsule, or other composite structure is analyzed via a multiscale approach. The objective is to achieve predictability and to preserve key atomic-scale features that underlie the formation and stability of the composite structures. We start with an all-atom description, the Liouville equation, and the order parameters characterizing nanoscale features of the system. An equation of Smoluchowski type for the stochastic dynamics of the order parameters is derived from the Liouville equation via a multiscale perturbation technique. The self-assembly of composite structures from nanocomponents with internal atomic structure is analyzed and growth rates are derived. Applications include the assembly of a viral capsid from capsomers, a ribosome from its major subunits, and composite materials from fibers and nanoparticles. Our approach overcomes errors in other coarse-graining methods which neglect the influence of the nanoscale configuration on the atomistic fluctuations. We account for the effect of order parameters on the statistics of the atomistic fluctuations which contribute to the entropic and average forces driving order parameter evolution. This approach enables an efficient algorithm for computer simulation of self-assembly, whereas other methods severely limit the timestep due to the separation of diffusional and complexing characteristic times. Given that our approach does not require recalibration with each new application, it provides a way to estimate assembly rates and thereby facilitate the discovery of self-assembly pathways and kinetic dead-end structures.

preprint2010arXiv

Stochastic Dynamics of Bionanosystems: Multiscale Analysis and Specialized Ensembles

An approach for simulating bionanosystems, such as viruses and ribosomes, is presented. This calibration-free approach is based on an all-atom description for bionanosystems, a universal interatomic force field, and a multiscale perspective. The supramillion-atom nature of these bionanosystems prohibits the use of a direct molecular dynamics approach for phenomena like viral structural transitions or self-assembly that develop over milliseconds or longer. A key element of these multiscale systems is the cross-talk between, and consequent strong coupling of, processes over many scales in space and time. We elucidate the role of interscale cross-talk and overcome bionanosystem simulation difficulties with automated construction of order parameters (OPs) describing supra-nanometer scale structural features, construction of OP dependent ensembles describing the statistical properties of atomistic variables that ultimately contribute to the entropies driving the dynamics of the OPs, and the derivation of a rigorous equation for the stochastic dynamics of the OPs. Since the atomic scale features of the system are treated statistically, several ensembles are constructed that reflect various experimental conditions. The theory provides a basis for a practical, quantitative bionanosystem modeling approach that preserves the cross-talk between the atomic and nanoscale features. A method for integrating information from nanotechnical experimental data in the derivation of equations of stochastic OP dynamics is also introduced.

preprint2009arXiv

Large Time Behavior of the Relativistic Vlasov Maxwell System in Low Space Dimension

When particle speeds are large the motion of a collisionless plasma is modeled by the relativistic Vlasov Maxwell system. Large time behavior of solutions which depend on one position variable and two momentum variables is considered. In the case of a single species of charge it is shown that there are solutions for which the charge density does not decay in time. This is in marked contrast to results for the non-relativistic Vlasov Poisson system in one space dimension. The case when two oppositely charged species are present and the net total charge is zero is also considered. In this case, it is shown that the support in the first component of momentum can grow at most like t to the three-fourths power.