Kinetic and spectral descriptions of autoionization phenomena associated with atomic processes in plasmas

NASA Astrophysics Data System (ADS)

Jacobs, Verne L.

2017-06-01

This investigation has been devoted to the theoretical description and computer modeling of atomic processes giving rise to radiative emission in energetic electron and ion beam interactions and in laboratory plasmas. We are also interested in the effects of directed electron and ion collisions and of anisotropic electric and magnetic fields. In the kinetic-theory description, we treat excitation, de-excitation, ionization, and recombination in electron and ion encounters with partially ionized atomic systems, including the indirect contributions from processes involving autoionizing resonances. These fundamental collisional and electromagnetic interactions also provide particle and photon transport mechanisms. From the spectral perspective, the analysis of atomic radiative emission can reveal detailed information on the physical properties in the plasma environment, such as non-equilibrium electron and charge-state distributions as well as electric and magnetic field distributions. In this investigation, a reduced-density-matrix formulation is developed for the microscopic description of atomic electromagnetic interactions in the presence of environmental (collisional and radiative) relaxation and decoherence processes. Our central objective is a fundamental microscopic description of atomic electromagnetic processes, in which both bound-state and autoionization-resonance phenomena can be treated in a unified and self-consistent manner. The time-domain (equation-of-motion) and frequency-domain (resolvent-operator) formulations of the reduced-density-matrix approach are developed in a unified and self-consistent manner. This is necessary for our ultimate goal of a systematic and self-consistent treatment of non-equilibrium (possibly coherent) atomic-state kinetics and high-resolution (possibly overlapping) spectral-line shapes. We thereby propose the introduction of a generalized collisional-radiative atomic-state kinetics model based on a reduced-density-matrix formulation. It will become apparent that the full atomic data needs for the precise modeling of extreme non-equilibrium plasma environments extend beyond the conventional radiative-transition-probability and collisional-cross-section data sets.

Self-consistent electrostatic potential due to trapped plasma in the magnetosphere

NASA Technical Reports Server (NTRS)

Miller, Ronald H.; Khazanov, George V.

1993-01-01

A steady state solution for the self-consistent electrostatic potential due to a plasma confined in a magnetic flux tube is considered. A steady state distribution function is constructed for the trapped particles from the constants of the motion, in the absence of waves and collisions. Using Liouville's theorem, the particle density along the geomagnetic field is determined and found to depend on the local magnetic field, self-consistent electric potential, and the equatorial plasma distribution function. A hot anisotropic magnetospheric plasma in steady state is modeled by a bi-Maxwellian at the equator. The self-consistent electric potential along the magnetic field is calculated assuming quasineutrality, and the potential drop is found to be approximately equal to the average kinetic energy of the equatorially trapped plasma. The potential is compared with that obtained by Alfven and Faelthammar (1963).

Effects of electric field methods on modeling the midlatitude ionospheric electrodynamics and inner magnetosphere dynamics

DOE PAGES

Yu, Yiqun; Jordanova, Vania Koleva; Ridley, Aaron J.; ...

2017-05-10

Here, we report a self-consistent electric field coupling between the midlatitude ionospheric electrodynamics and inner magnetosphere dynamics represented in a kinetic ring current model. This implementation in the model features another self-consistency in addition to its already existing self-consistent magnetic field coupling with plasma. The model is therefore named as Ring current-Atmosphere interaction Model with Self-Consistent magnetic (B) and electric (E) fields, or RAM-SCB-E. With this new model, we explore, by comparing with previously employed empirical Weimer potential, the impact of using self-consistent electric fields on the modeling of storm time global electric potential distribution, plasma sheet particle injection, andmore » the subauroral polarization streams (SAPS) which heavily rely on the coupled interplay between the inner magnetosphere and midlatitude ionosphere. We find the following phenomena in the self-consistent model: (1) The spatially localized enhancement of electric field is produced within 2.5 < L < 4 during geomagnetic active time in the dusk-premidnight sector, with a similar dynamic penetration as found in statistical observations. (2) The electric potential contours show more substantial skewing toward the postmidnight than the Weimer potential, suggesting the resistance on the particles from directly injecting toward the low-L region. (3) The proton flux indeed indicates that the plasma sheet inner boundary at the dusk-premidnight sector is located further away from the Earth than in the Weimer potential, and a “tongue” of low-energy protons extends eastward toward the dawn, leading to the Harang reversal. (4) SAPS are reproduced in the subauroral region, and their magnitude and latitudinal width are in reasonable agreement with data.« less

Effects of electric field methods on modeling the midlatitude ionospheric electrodynamics and inner magnetosphere dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yu, Yiqun; Jordanova, Vania Koleva; Ridley, Aaron J.

Here, we report a self-consistent electric field coupling between the midlatitude ionospheric electrodynamics and inner magnetosphere dynamics represented in a kinetic ring current model. This implementation in the model features another self-consistency in addition to its already existing self-consistent magnetic field coupling with plasma. The model is therefore named as Ring current-Atmosphere interaction Model with Self-Consistent magnetic (B) and electric (E) fields, or RAM-SCB-E. With this new model, we explore, by comparing with previously employed empirical Weimer potential, the impact of using self-consistent electric fields on the modeling of storm time global electric potential distribution, plasma sheet particle injection, andmore » the subauroral polarization streams (SAPS) which heavily rely on the coupled interplay between the inner magnetosphere and midlatitude ionosphere. We find the following phenomena in the self-consistent model: (1) The spatially localized enhancement of electric field is produced within 2.5 < L < 4 during geomagnetic active time in the dusk-premidnight sector, with a similar dynamic penetration as found in statistical observations. (2) The electric potential contours show more substantial skewing toward the postmidnight than the Weimer potential, suggesting the resistance on the particles from directly injecting toward the low-L region. (3) The proton flux indeed indicates that the plasma sheet inner boundary at the dusk-premidnight sector is located further away from the Earth than in the Weimer potential, and a “tongue” of low-energy protons extends eastward toward the dawn, leading to the Harang reversal. (4) SAPS are reproduced in the subauroral region, and their magnitude and latitudinal width are in reasonable agreement with data.« less

Effects of electric field methods on modeling the midlatitude ionospheric electrodynamics and inner magnetosphere dynamics

NASA Astrophysics Data System (ADS)

Yu, Yiqun; Jordanova, Vania K.; Ridley, Aaron J.; Toth, Gabor; Heelis, Roderick

2017-05-01

We report a self-consistent electric field coupling between the midlatitude ionospheric electrodynamics and inner magnetosphere dynamics represented in a kinetic ring current model. This implementation in the model features another self-consistency in addition to its already existing self-consistent magnetic field coupling with plasma. The model is therefore named as Ring current-Atmosphere interaction Model with Self-Consistent magnetic (B) and electric (E) fields, or RAM-SCB-E. With this new model, we explore, by comparing with previously employed empirical Weimer potential, the impact of using self-consistent electric fields on the modeling of storm time global electric potential distribution, plasma sheet particle injection, and the subauroral polarization streams (SAPS) which heavily rely on the coupled interplay between the inner magnetosphere and midlatitude ionosphere. We find the following phenomena in the self-consistent model: (1) The spatially localized enhancement of electric field is produced within 2.5 < L < 4 during geomagnetic active time in the dusk-premidnight sector, with a similar dynamic penetration as found in statistical observations. (2) The electric potential contours show more substantial skewing toward the postmidnight than the Weimer potential, suggesting the resistance on the particles from directly injecting toward the low-L region. (3) The proton flux indeed indicates that the plasma sheet inner boundary at the dusk-premidnight sector is located further away from the Earth than in the Weimer potential, and a "tongue" of low-energy protons extends eastward toward the dawn, leading to the Harang reversal. (4) SAPS are reproduced in the subauroral region, and their magnitude and latitudinal width are in reasonable agreement with data.

Development of high performance particle in cell code for the exascale age

NASA Astrophysics Data System (ADS)

Lapenta, Giovanni; Amaya, Jorge; Gonzalez, Diego; Deep-Est H2020 Consortium Collaboration

2017-10-01

Magnetized plasmas are most effectively described by magneto-hydrodynamics, MHD, a fluid theory based on describing some fields defined in space: electromagnetic fields, density, velocity and temperature of the plasma. However, microphysics processes need kinetic theory, where statistical distributions of particles are governed by the Boltzmann equation. While fluid models are based on the ordinary space and time, kinetic models require a six dimensional space, called phase space, besides time. The two methods are not separated but rather interact to determine the system evolution. Arriving at a single self-consistent model is the goal of our research. We present a new approach developed with the goal of extending the reach of kinetic models to the fluid scales. Kinetic models are a higher order description and all fluid effects are included in them. However, the cost in terms of computing power is much higher and it has been so far prohibitively expensive to treat space weather events fully kinetically. We have now designed a new method capable of reducing that cost by several orders of magnitude making it possible for kinetic models to study macroscopic systems. H2020 Deep-EST consortium (European Commission).

Localized mRNA translation and protein association

NASA Astrophysics Data System (ADS)

Zhdanov, Vladimir P.

2014-08-01

Recent direct observations of localization of mRNAs and proteins both in prokaryotic and eukaryotic cells can be related to slowdown of diffusion of these species due to macromolecular crowding and their ability to aggregate and form immobile or slowly mobile complexes. Here, a generic kinetic model describing both these factors is presented and comprehensively analyzed. Although the model is non-linear, an accurate self-consistent analytical solution of the corresponding reaction-diffusion equation has been constructed, the types of localized protein distributions have been explicitly shown, and the predicted kinetic regimes of gene expression have been classified.

Role of filament annealing in the kinetics and thermodynamics of nucleated polymerization.

PubMed

Michaels, Thomas C T; Knowles, Tuomas P J

2014-06-07

The formation of nanoscale protein filaments from soluble precursor molecules through nucleated polymerization is a common form of supra-molecular assembly phenomenon. This process underlies the generation of a range of both functional and pathological structures in nature. Filament breakage has emerged as a key process controlling the kinetics of the growth reaction since it increases the number of filament ends in the system that can act as growth sites. In order to ensure microscopic reversibility, however, the inverse process of fragmentation, end-to-end annealing of filaments, is a necessary component of a consistent description of such systems. Here, we combine Smoluchowski kinetics with nucleated polymerization models to generate a master equation description of protein fibrillization, where filamentous structures can undergo end-to-end association, in addition to elongation, fragmentation, and nucleation processes. We obtain self-consistent closed-form expressions for the growth kinetics and discuss the key physics that emerges from considering filament fusion relative to current fragmentation only models. Furthermore, we study the key time scales that describe relaxation to equilibrium.

UCLA IGPP Space Plasma Simulation Group

NASA Technical Reports Server (NTRS)

1998-01-01

During the past 10 years the UCLA IGPP Space Plasma Simulation Group has pursued its theoretical effort to develop a Mission Oriented Theory (MOT) for the International Solar Terrestrial Physics (ISTP) program. This effort has been based on a combination of approaches: analytical theory, large scale kinetic (LSK) calculations, global magnetohydrodynamic (MHD) simulations and self-consistent plasma kinetic (SCK) simulations. These models have been used to formulate a global interpretation of local measurements made by the ISTP spacecraft. The regions of applications of the MOT cover most of the magnetosphere: the solar wind, the low- and high-latitude magnetospheric boundary, the near-Earth and distant magnetotail, and the auroral region. Most recent investigations include: plasma processes in the electron foreshock, response of the magnetospheric cusp, particle entry in the magnetosphere, sources of observed distribution functions in the magnetotail, transport of oxygen ions, self-consistent evolution of the magnetotail, substorm studies, effects of explosive reconnection, and auroral acceleration simulations.

Object-oriented code SUR for plasma kinetic simulation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Levchenko, V.D.; Sigov, Y.S.

1995-12-31

We have developed a self-consistent simulation code based on object-oriented model of plasma (OOMP) for solving the Vlasov/Poisson (V/P), Vlasov/Maxwell (V/M), Bhatnagar-Gross-Krook (BGK) as well as Fokker-Planck (FP) kinetic equations. The application of an object-oriented approach (OOA) to simulation of plasmas and plasma-like media by means of splitting methods permits to uniformly describe and solve the wide circle of plasma kinetics problems, including those being very complicated: many-dimensional, relativistic, with regard for collisions, specific boundary conditions etc. This paper gives the brief description of possibilities of the SUR code, as a concrete realization of OOMP.

Three dimensional fluid-kinetic model of a magnetically guided plasma jet

NASA Astrophysics Data System (ADS)

Ramos, Jesús J.; Merino, Mario; Ahedo, Eduardo

2018-06-01

A fluid-kinetic model of the collisionless plasma flow in a convergent-divergent magnetic nozzle is presented. The model combines the leading-order Vlasov equation and the fluid continuity and perpendicular momentum equation for magnetized electrons, and the fluid equations for cold ions, which must be solved iteratively to determine the self-consistent plasma response in a three-dimensional magnetic field. The kinetic electron solution identifies three electron populations and provides the plasma density and pressure tensor. The far downstream asymptotic behavior shows the anisotropic cooling of the electron populations. The fluid equations determine the electric potential and the fluid velocities. In the small ion-sound gyroradius case, the solution is constructed one magnetic line at a time. In the large ion-sound gyroradius case, ion detachment from magnetic lines makes the problem fully three-dimensional.

Physical Kinetics of Electrons in a High-Voltage Pulsed High-Pressure Discharge with Cylindrical Geometry

NASA Astrophysics Data System (ADS)

Kozhevnikov, V. Yu.; Kozyrev, A. V.; Semeniuk, N. S.

2017-12-01

Results of theoretical modeling of the phenomenon of a high-voltage discharge in nitrogen at atmospheric pressure are presented, based on a consistent kinetic theory of the electrons. A mathematical model of a nonstationary high-pressure discharge has been constructed for the first time, based on a description of the electron component from first principles. The physical kinetics of the electrons are described with the help of the Boltzmann kinematic equation for the electron distribution function over momenta with only ionization and elastic collisions taken into account. A detailed spatiotemporal picture of a nonstationary discharge with runaway electrons under conditions of coaxial geometry of the gas diode is presented. The model describes in a self-consistent way both the process of formation of the runaway electron flux in the discharge and the influence of this flux on the rate of ionization processes in the gas. Total energy spectra of the electron flux incident on the anode are calculated. The obtained parameters of the current pulse of the beam of fast electrons correlate well with the known experimental data.

Berkeley Lab - Materials Sciences Division

Science.gov Websites

. Orlita, L. Z. Tan, M. Potemski, M. Sprinkle, C. Berger, W. A. de Heer, S. G. Louie and G. Martinez . DePaolo and J. J. De Yoreo. Self-consistent ion-by-ion growth model for kinetic isotopic fractionation Biosynthesis Restricts Mycobacterium tuberculosis Growth in Human Macrophages. ACS Chemical Biology 7, 863

Interfacial mixing in high energy-density matter with a multiphysics kinetic model

NASA Astrophysics Data System (ADS)

Haack, Jeff; Hauck, Cory; Murillo, Michael

2017-10-01

We have extended a recently-developed multispecies, multitemperature BGK model to include multiphysics capability that allows modeling of a wider range of plasma conditions. In particular, we have extended the model to describe one spatial dimension, and included a multispecies atomic ionization model, accurate collision physics across coupling regimes, self-consistent electric fields, and degeneracy in the electronic screening. We apply the new model to a warm dense matter scenario in which the ablator-fuel interface of an inertial confinement fusion target is heated, similar to a recent molecular dynamics study, but for larger length and time scales and for much higher temperatures. From our numerical results we are able to explore a variety of phenomena, including hydrogen jetting, kinetic effects (non-Maxwellian and anisotropic distributions), plasma physics (size, persistence and role of electric fields) and transport (relative sizes of Fickean diffision, electrodiffusion and barodiffusion). As compared with the recent molecular dynamics work the kinetic model greatly extends the accessible physical domains we are able to model.

Study of electron transport in a Hall thruster by axial–radial fully kinetic particle simulation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cho, Shinatora, E-mail: choh.shinatora@jaxa.jp; Kubota, Kenichi; Funaki, Ikkoh

2015-10-15

Electron transport across a magnetic field in a magnetic-layer-type Hall thruster was numerically investigated for the future predictive modeling of Hall thrusters. The discharge of a 1-kW-class magnetic-layer-type Hall thruster designed for high-specific-impulse operation was modeled using an r-z two-dimensional fully kinetic particle code with and without artificial electron-diffusion models. The thruster performance results showed that both electron transport models captured the experimental result within discrepancies less than 20% in thrust and discharge current for all the simulated operation conditions. The electron cross-field transport mechanism of the so-called anomalous diffusion was self-consistently observed in the simulation without artificial diffusion models;more » the effective electron mobility was two orders of magnitude higher than the value obtained using the classical diffusion theory. To account for the self-consistently observed anomalous transport, the oscillation of plasma properties was speculated. It was suggested that the enhanced random-walk diffusion due to the velocity oscillation of low-frequency electron flow could explain the observed anomalous diffusion within an order of magnitude. The dominant oscillation mode of the electron flow velocity was found to be 20 kHz, which was coupled to electrostatic oscillation excited by global ionization instability.« less

Global Particle-in-Cell Simulations of Mercury's Magnetosphere

NASA Astrophysics Data System (ADS)

Schriver, D.; Travnicek, P. M.; Lapenta, G.; Amaya, J.; Gonzalez, D.; Richard, R. L.; Berchem, J.; Hellinger, P.

2017-12-01

Spacecraft observations of Mercury's magnetosphere have shown that kinetic ion and electron particle effects play a major role in the transport, acceleration, and loss of plasma within the magnetospheric system. Kinetic processes include reconnection, the breakdown of particle adiabaticity and wave-particle interactions. Because of the vast range in spatial scales involved in magnetospheric dynamics, from local electron Debye length scales ( meters) to solar wind/planetary magnetic scale lengths (tens to hundreds of planetary radii), fully self-consistent kinetic simulations of a global planetary magnetosphere remain challenging. Most global simulations of Earth's and other planet's magnetosphere are carried out using MHD, enhanced MHD (e.g., Hall MHD), hybrid, or a combination of MHD and particle in cell (PIC) simulations. Here, 3D kinetic self-consistent hybrid (ion particle, electron fluid) and full PIC (ion and electron particle) simulations of the solar wind interaction with Mercury's magnetosphere are carried out. Using the implicit PIC and hybrid simulations, Mercury's relatively small, but highly kinetic magnetosphere will be examined to determine how the self-consistent inclusion of electrons affects magnetic reconnection, particle transport and acceleration of plasma at Mercury. Also the spatial and energy profiles of precipitating magnetospheric ions and electrons onto Mercury's surface, which can strongly affect the regolith in terms of space weathering and particle outflow, will be examined with the PIC and hybrid codes. MESSENGER spacecraft observations are used both to initiate and validate the global kinetic simulations to achieve a deeper understanding of the role kinetic physics play in magnetospheric dynamics.

Ion kinetic effects on the ignition and burn of inertial confinement fusion targets: A multi-scale approach

NASA Astrophysics Data System (ADS)

Peigney, B. E.; Larroche, O.; Tikhonchuk, V.

2014-12-01

In this article, we study the hydrodynamics and burn of the thermonuclear fuel in inertial confinement fusion pellets at the ion kinetic level. The analysis is based on a two-velocity-scale Vlasov-Fokker-Planck kinetic model that is specially tailored to treat fusion products (suprathermal α-particles) in a self-consistent manner with the thermal bulk. The model assumes spherical symmetry in configuration space and axial symmetry in velocity space around the mean flow velocity. A typical hot-spot ignition design is considered. Compared with fluid simulations where a multi-group diffusion scheme is applied to model α transport, the full ion-kinetic approach reveals significant non-local effects on the transport of energetic α-particles. This has a direct impact on hydrodynamic spatial profiles during combustion: the hot spot reactivity is reduced, while the inner dense fuel layers are pre-heated by the escaping α-suprathermal particles, which are transported farther out of the hot spot. We show how the kinetic transport enhancement of fusion products leads to a significant reduction of the fusion yield.

Ion kinetic effects on the ignition and burn of inertial confinement fusion targets: A multi-scale approach

DOE Office of Scientific and Technical Information (OSTI.GOV)

Peigney, B. E.; Larroche, O.; Tikhonchuk, V.

2014-12-15

In this article, we study the hydrodynamics and burn of the thermonuclear fuel in inertial confinement fusion pellets at the ion kinetic level. The analysis is based on a two-velocity-scale Vlasov-Fokker-Planck kinetic model that is specially tailored to treat fusion products (suprathermal α-particles) in a self-consistent manner with the thermal bulk. The model assumes spherical symmetry in configuration space and axial symmetry in velocity space around the mean flow velocity. A typical hot-spot ignition design is considered. Compared with fluid simulations where a multi-group diffusion scheme is applied to model α transport, the full ion-kinetic approach reveals significant non-local effectsmore » on the transport of energetic α-particles. This has a direct impact on hydrodynamic spatial profiles during combustion: the hot spot reactivity is reduced, while the inner dense fuel layers are pre-heated by the escaping α-suprathermal particles, which are transported farther out of the hot spot. We show how the kinetic transport enhancement of fusion products leads to a significant reduction of the fusion yield.« less

The Interaction of Coronal Mass Ejections with Alfvénic Turbulence

NASA Astrophysics Data System (ADS)

Manchester, Ward, IV; Van Der Holst, Bart

2017-09-01

We provide a first attempt to understand the interaction between Alfvén wave turbulence, kinetic instabilities and temperature anisotropies in the environment of a fast coronal mass ejection (CME) near the Sun. The impact of a fast CME on the solar corona causes turbulent energy, thermal energy and dissipative heating to increase by orders of magnitude, and produces conditions suitable for a host of kinetic instabilities. We study these CME-induced effects with the recently developed Alfvén Wave Solar Model, with which we are able to self-consistently simulate the turbulent energy transport and dissipation as well as isotropic electron heating and anisotropic proton heating. Furthermore, the model also offers the capability to address the effects of fire hose, mirror mode, and cyclotron kinetic instabilities on proton energy partitioning all in a global-scale numerical simulation. We find amplified turbulent energy in the CME sheath, along with strong wave reflection at the shock combine to cause wave dissipation rates to increase by more than a factor of 100. In contrast, wave energy is greatly diminished by adiabatic expansion in the flux rope. Finally, we find proton temperature anisotropies are limited by kinetic instabilities to a level consistent with solar wind observations.

The Interaction of Coronal Mass Ejections with Alfvenic Turbulence

NASA Astrophysics Data System (ADS)

Manchester, W.; van der Holst, B.

2017-12-01

We provide a first attempt to understand the interaction between Alfven wave turbulence, kinetic instabilities and temperature anisotropies in the environment of a fast coronal mass ejection (CME). The impact of a fast CME on the solar corona causes turbulent energy, thermal energy and dissipative heating to increase by orders of magnitude, and produces conditions suitable for a host of kinetic instabilities. We study these CME-induced effects with the recently developed Alfven Wave Solar Model, with which we are able to self-consistently simulate the turbulent energy transport and dissipation as well as isotropic electron heating and anisotropic proton heating. Furthermore, the model also offers the capability to address the effects of firehose, mirror mode, and cyclotron kinetic instabilities on proton energy partitioning, all in a global-scale numerical simulation. We find turbulent energy greatly enhanced in the CME sheath, strong wave reflection at the shock, which leads to wave dissipation rates increasing by more than a factor of 100. In contrast, wave energy is greatly diminished by adiabatic expansion in the flux rope. Finally, we find proton temperature anisotropies are limited by kinetic instabilities to a level consistent with solar wind observations.

A Theory for Self-consistent Acceleration of Energetic Charged Particles by Dynamic Small-scale Flux Ropes

NASA Astrophysics Data System (ADS)

le Roux, J. A.; Zank, G. P.; Khabarova, O.; Webb, G. M.

2016-12-01

Simulations of charged particle acceleration in turbulent plasma regions with numerous small-scale contracting and merging (reconnecting) magnetic islands/flux ropes emphasize the key role of temporary particle trapping in these structures for efficient acceleration that can result in power-law spectra. In response, a comprehensive kinetic transport theory framework was developed by Zank et al. and le Roux et al. to capture the essential physics of energetic particle acceleration in solar wind regions containing numerous dynamic small-scale flux ropes. Examples of test particle solutions exhibiting hard power-law spectra for energetic particles were presented in recent publications by both Zank et al. and le Roux et al.. However, the considerable pressure in the accelerated particles suggests the need for expanding the kinetic transport theory to enable a self-consistent description of energy exchange between energetic particles and small-scale flux ropes. We plan to present the equations of an expanded kinetic transport theory framework that will enable such a self-consistent description.

A Self-Consistent Model of the Interacting Ring Current Ions and Electromagnetic Ion Cyclotron Waves, Initial Results: Waves and Precipitating Fluxes

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Gamayunov, K. V.; Jordanova, V. K.; Krivorutsky, E. N.

2002-01-01

Initial results from a newly developed model of the interacting ring current ions and ion cyclotron waves are presented. The model is based on the system of two kinetic equations: one equation describes the ring current ion dynamics, and another equation describes wave evolution. The system gives a self-consistent description of the ring current ions and ion cyclotron waves in a quasilinear approach. These equations for the ion phase space distribution function and for the wave power spectral density were solved on aglobal magnetospheric scale undernonsteady state conditions during the 2-5 May 1998 storm. The structure and dynamics of the ring current proton precipitating flux regions and the ion cyclotron wave-active zones during extreme geomagnetic disturbances on 4 May 1998 are presented and discussed in detail.

A Self-Consistent Model of the Interacting Ring Current Ions and Electromagnetic ICWs. Initial Results: Waves and Precipitation Fluxes

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Gamayunov, K. V.; Jordanova, V. K.; Krivorutsky, E. N.; Whitaker, Ann F. (Technical Monitor)

2001-01-01

Initial results from the new developed model of the interacting ring current ions and ion cyclotron waves are presented. The model described by the system of two bound kinetic equations: one equation describes the ring current ion dynamics, and another one gives wave evolution. Such system gives a self-consistent description of the ring current ions and ion cyclotron waves in a quasilinear approach. Calculating ion-wave relationships, on a global scale under non steady-state conditions during May 2-5, 1998 storm, we presented the data at three time cuts around initial, main, and late recovery phases of May 4, 1998 storm phase. The structure and dynamics of the ring current proton precipitating flux regions and the wave active ones are discussed in detail.

A Self-Consistent Model of the Interacting Ring Current Ions with Electromagnetic ICWs

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Gamayunov, K. V.; Jordanova, V. K.; Krivorutsky, E. N.; Whitaker, Ann F. (Technical Monitor)

2001-01-01

Initial results from a newly developed model of the interacting ring current ions and ion cyclotron waves are presented. The model is based on the system of two bound kinetic equations: one equation describes the ring current ion dynamics, and another equation describes wave evolution. The system gives a self-consistent description of ring current ions and ion cyclotron waves in a quasilinear approach. These two equations were solved on a global scale under non steady-state conditions during the May 2-5, 1998 storm. The structure and dynamics of the ring current proton precipitating flux regions and the wave active zones at three time cuts around initial, main, and late recovery phases of the May 4, 1998 storm phase are presented and discussed in detail. Comparisons of the model wave-ion data with the Polar/HYDRA and Polar/MFE instruments results are presented..

Effect of EMIC Wave Normal Angle Distribution on Relativistic Electron Scattering Based on the Newly Developed Self-consistent RC/EMIC Waves Model by Khazanov et al. [2006

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Gallagher, D. L.; Gamayunov, K.

2007-01-01

It is well known that the effects of EMIC waves on RC ion and RB electron dynamics strongly depend on such particle/wave characteristics as the phase-space distribution function, frequency, wave-normal angle, wave energy, and the form of wave spectral energy density. Therefore, realistic characteristics of EMIC waves should be properly determined by modeling the RC-EMIC waves evolution self-consistently. Such a selfconsistent model progressively has been developing by Khaznnov et al. [2002-2006]. It solves a system of two coupled kinetic equations: one equation describes the RC ion dynamics and another equation describes the energy density evolution of EMIC waves. Using this model, we present the effectiveness of relativistic electron scattering and compare our results with previous work in this area of research.

Synthetic spectral analysis of a kinetic model for slow-magnetosonic waves in solar corona

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ruan, Wenzhi; He, Jiansen; Tu, Chuanyi

We propose a kinetic model of slow-magnetosonic waves to explain various observational features associated with the propagating intensity disturbances (PIDs) occurring in the solar corona. The characteristics of slow mode waves, e.g, inphase oscillations of density, velocity, and thermal speed, are reproduced in this kinetic model. Moreover, the red-blue (R-B) asymmetry of the velocity distribution as self-consistently generated in the model is found to be contributed from the beam component, as a result of the competition between Landau resonance and Coulomb collisions. Furthermore, we synthesize the spectral lines and make the spectral analysis, based on the kinetic simulation data ofmore » the flux tube plasmas and the hypothesis of the surrounding background plasmas. It is found that the fluctuations of parameters of the synthetic spectral lines are basically consistent with the observations: (1) the line intensity, Doppler shift, and line width are fluctuating in phase; (2) the R-B asymmetry usually oscillate out of phase with the former three parameters; (3) the blueward asymmetry is more evident than the redward asymmetry in the R-B fluctuations. The oscillations of line parameters become weakened for the case with denser surrounding background plasmas. Similar to the observations, there is no doubled-frequency oscillation of the line width for the case with flux-tube plasmas flowing bulkly upward among the static background plasmas. Therefore, we suggest that the “wave + beam flow” kinetic model may be a viable interpretation for the PIDs observed in the solar corona.« less

Transport in semiconductor nanowire superlattices described by coupled quantum mechanical and kinetic models.

PubMed

Alvaro, M; Bonilla, L L; Carretero, M; Melnik, R V N; Prabhakar, S

2013-08-21

In this paper we develop a kinetic model for the analysis of semiconductor superlattices, accounting for quantum effects. The model consists of a Boltzmann-Poisson type system of equations with simplified Bhatnagar-Gross-Krook collisions, obtained from the general time-dependent Schrödinger-Poisson model using Wigner functions. This system for superlattice transport is supplemented by the quantum mechanical part of the model based on the Ben-Daniel-Duke form of the Schrödinger equation for a cylindrical superlattice of finite radius. The resulting energy spectrum is used to characterize the Fermi-Dirac distribution that appears in the Bhatnagar-Gross-Krook collision, thereby coupling the quantum mechanical and kinetic parts of the model. The kinetic model uses the dispersion relation obtained by the generalized Kronig-Penney method, and allows us to estimate radii of quantum wire superlattices that have the same miniband widths as in experiments. It also allows us to determine more accurately the time-dependent characteristics of superlattices, in particular their current density. Results, for several experimentally grown superlattices, are discussed in the context of self-sustained coherent oscillations of the current density which are important in an increasing range of current and potential applications.

A quasilinear kinetic model for solar wind electrons and protons instabilities

NASA Astrophysics Data System (ADS)

Sarfraz, M.; Yoon, P. H.

2017-12-01

In situ measurements confirm the anisotropic behavior in temperatures of solar wind species. These anisotropies associated with charge particles are observed to be relaxed. In collionless limit, kinetic instabilities play a significant role to reshape particles distribution. The linear analysis results are encapsulated in inverse relationship between anisotropy and plasma beta based observations fittings techniques, simulations methods, or solution of linearized Vlasov equation. Here amacroscopic quasilinear technique is adopted to confirm inverse relationship through solutions of set of self-consistent kinetic equations. Firstly, for a homogeneous and non-collisional medium, quasilinear kinetic model is employed to display asymptotic variations of core and halo electrons temperatures and saturations of wave energy densities for electromagnetic electron cyclotron (EMEC) instability sourced by, T⊥}>T{∥ . It is shown that, in (β ∥ , T⊥}/T{∥ ) phase space, the saturations stages of anisotropies associated with core and halo electrons lined up on their respective marginal stability curves. Secondly, for case of electrons firehose instability ignited by excessive parallel temperature i.e T⊥}>T{∥ , both electrons and protons are allowed to dynamically evolve in time. It is also observed that, the trajectories of protons and electrons at saturation stages in phase space of anisotropy and plasma beta correspond to proton cyclotron and firehose marginal stability curves, respectively. Next, the outstanding issue that most of observed proton data resides in nearly isotropic state in phase space is interpreted. Here, in quasilinear frame-work of inhomogeneous solar wind system, a set of self-consistent quasilinear equations is formulated to show a dynamical variations of temperatures with spatial distributions. On choice of different initial parameters, it is shown that, interplay of electron and proton instabilities provides an counter-balancing force to slow down the protons away from marginal stability states. As we are dealing both, protons and electrons for radially expanding solar wind plasma, our present approach may eventually be incorporated in global-kinetic models of the solar wind species.

A model for the nonlocal transport and the associated distribution function deformation in magnetized laser-plasmas

NASA Astrophysics Data System (ADS)

Nicolaï, Ph.; Feugeas, J.-L.; Schurtz, G.

2006-06-01

We present a model of nonlocal transport for multidimensional radiation magneto hydrodynamic codes. In laser produced plasmas, it is now believed that the heat transfert can be strongly modified by the nonlocal nature of the electron conduction. Nevertheless other mechanisms as self generated magnetic fields may affect heat transport too. The model described in this work aims at extending the formula of G. Schurtz, Ph. Nicolaï and M. Busquet [1] to magnetized plasmas. A system of nonlocal equations is derived from kinetic equations with self-consistent electric and magnetic fields. These equations are analyzed and applied to a physical problem in order to demonstrate the main features of the model.

The role of chlorine chemistry in Antarctic ozone loss - Implications of new kinetic data

NASA Technical Reports Server (NTRS)

Rodriguez, Jose M.; Ko, Malcolm K. W.; Sze, Nien Dak

1990-01-01

New kinetic data yielding a slower formation rate and larger absorption cross sections of Cl2O2 are incorporated into a photochemical model to reassess the role of chlorine chemistry in accounting for the ozone reductions derived from TOMS observations in different regions of the Antarctic polar vortex during 1987. The model is further constrained by existing measurements from the Airborne Antarctic Ozone Experiment and the National Ozone Expedition II. Calculated concentrations of ClO based on the new kinetic data increase by almost a factor of two between the collar and core regions of the vortex during the second half of September. The calculated ozone reductions in the vortex core appear to be consistent with the TOMS observations in spite of the slower rate for the self-reaction of ClO.

Toroidal Ampere-Faraday Equations Solved Consistently with the CQL3D Fokker-Planck Time-Evolution

NASA Astrophysics Data System (ADS)

Harvey, R. W.; Petrov, Yu. V.

2013-10-01

A self-consistent, time-dependent toroidal electric field calculation is a key feature of a complete 3D Fokker-Planck kinetic distribution radial transport code for f(v,theta,rho,t). In the present CQL3D finite-difference model, the electric field E(rho,t) is either prescribed, or iteratively adjusted to obtain prescribed toroidal or parallel currents. We discuss first results of an implementation of the Ampere-Faraday equation for the self-consistent toroidal electric field, as applied to the runaway electron production in tokamaks due to rapid reduction of the plasma temperature as occurs in a plasma disruption. Our previous results assuming a constant current density (Lenz' Law) model showed that prompt ``hot-tail runaways'' dominated ``knock-on'' and Dreicer ``drizzle'' runaways; we will examine modifications due to the more complete Ampere-Faraday solution. Work supported by US DOE under DE-FG02-ER54744.

Torsion as a source of expansion in a Bianchi type-I universe in the self-consistent Einstein-Cartan theory of a perfect fluid with spin density

NASA Technical Reports Server (NTRS)

Bradas, James C.; Fennelly, Alphonsus J.; Smalley, Larry L.

1987-01-01

It is shown that a generalized (or 'power law') inflationary phase arises naturally and inevitably in a simple (Bianchi type-I) anisotropic cosmological model in the self-consistent Einstein-Cartan gravitation theory with the improved stress-energy-momentum tensor with the spin density of Ray and Smalley (1982, 1983). This is made explicit by an analytical solution of the field equations of motion of the fluid variables. The inflation is caused by the angular kinetic energy density due to spin. The model further elucidates the relationship between fluid vorticity, the angular velocity of the inertially dragged tetrads, and the precession of the principal axes of the shear ellipsoid. Shear is not effective in damping the inflation.

Kinetics of polyelectrolyte adsorption

NASA Astrophysics Data System (ADS)

Cohen Stuart, M. A.; Hoogendam, C. W.; de Keizer, A.

1997-09-01

The kinetics of polyelectrolyte adsorption has been investigated theoretically. In analogy with Kramers' rate theory for chemical reactions we present a model which is based on the assumption that a polyelectrolyte encounters a barrier in its motion towards an adsorbing surface. The height of the barrier, which is of electrostatic origin, is calculated with a self-consistent-field (SCF) model. The salt concentration strongly affects the height of the barrier. At moderate salt concentrations (0953-8984/9/37/009/img1) equilibrium in the adsorption is attained; at low salt concentration (0953-8984/9/37/009/img2) equilibrium is not reached on the time scale of experiments. The attachment process shows resemblances to the classical DLVO theory.

Ring Current Ion Coupling with Electromagnetic Ion Cyclotron Waves

NASA Technical Reports Server (NTRS)

Khazanov. G. V.; Gamayunov, K. V.; Jordanova, V. K.; Six, N. Frank (Technical Monitor)

2002-01-01

A new ring current global model has been developed that couples the system of two kinetic equations: one equation describes the ring current (RC) ion dynamic, and another equation describes wave evolution of electromagnetic ion cyclotron waves (EMIC). The coupled model is able to simulate, for the first time self-consistently calculated RC ion kinetic and evolution of EMIC waves that propagate along geomagnetic field lines and reflect from the ionosphere. Ionospheric properties affect the reflection index through the integral Pedersen and Hall conductivities. The structure and dynamics of the ring current proton precipitating flux regions, intensities of EMIC global RC energy balance, and some other parameters will be studied in detail for the selected geomagnetic storms.

Integrated fusion simulation with self-consistent core-pedestal coupling

DOE PAGES

Meneghini, O.; Snyder, P. B.; Smith, S. P.; ...

2016-04-20

In this study, accurate prediction of fusion performance in present and future tokamaks requires taking into account the strong interplay between core transport, pedestal structure, current profile and plasma equilibrium. An integrated modeling workflow capable of calculating the steady-state self- consistent solution to this strongly-coupled problem has been developed. The workflow leverages state-of-the-art components for collisional and turbulent core transport, equilibrium and pedestal stability. Validation against DIII-D discharges shows that the workflow is capable of robustly pre- dicting the kinetic profiles (electron and ion temperature and electron density) from the axis to the separatrix in good agreement with the experiments.more » An example application is presented, showing self-consistent optimization for the fusion performance of the 15 MA D-T ITER baseline scenario as functions of the pedestal density and ion effective charge Z eff.« less

Reply to "Comment on 'A Self-Consistent Model of the Interacting Ring Current Ions and Electromagnetic Ion Cyclotron Waves, Initial Results: Waves and Precipitation Fluxes' and 'Self-Consistent Model of the Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves in Multi-Ion Magnetosphere' by Khazanov et al. et al."

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Gamayunov, K. V.; Gallagher, D. L.; Kozyra, J. W.

2007-01-01

It is well-known that the effects of electromagnetic ion cyclotron (EMIC) waves on ring current (RC) ion and radiation belt (RB) electron dynamics strongly depend on such particle/wave characteristics as the phase-space distribution function, frequency, wavenormal angle, wave energy, and the form of wave spectral energy density. The consequence is that accurate modeling of EMIC waves and RC particles requires robust inclusion of the interdependent dynamics of wave growth/damping, wave propagation, and[ particles. Such a self-consistent model is being progressively developed by Khazanov et al. [2002, 2006, 2007]. This model is based on a system of coupled kinetic equations for the RC and EMIC wave power spectral density along with the ray tracing equations. Thome and Home [2007] (hereafter referred to as TH2007) call the Khazanov et al. [2002, 2006] results into question in their Comment. The points in contention can be summarized as follows. TH2007 claim that: (1) "the important damping of waves by thermal heavy ions is completely ignored", and Landau damping during resonant interaction with thermal electrons is not included in our model; (2) EMIC wave damping due to RC O + is not included in our simulation; (3) non-linear processes limiting EMIC wave amplitude are not included in our model; (4) growth of the background fluctuations to a physically significantamplitude"must occur during a single transit of the unstable region" with subsequent damping below bi-ion latitudes,and consequently"the bounce averaged wave kinetic equation employed in the code contains a physically erroneous 'assumption". Our reply will address each of these points as well as other criticisms mentioned in the Comment. TH2007 are focused on two of our papers that are separated by four years. Significant progress in the self-consistent treatment of the RC-EMIC wave system has been achieved during those years. The paper by Khazanov et al. [2006] presents the latest version of our model, and in this Reply we refer mostly to this paper.

Computational modelling of mesoscale dislocation patterning and plastic deformation of single crystals

NASA Astrophysics Data System (ADS)

Xia, Shengxu; El-Azab, Anter

2015-07-01

We present a continuum dislocation dynamics model that predicts the formation of dislocation cell structure in single crystals at low strains. The model features a set of kinetic equations of the curl type that govern the space and time evolution of the dislocation density in the crystal. These kinetic equations are coupled to stress equilibrium and deformation kinematics using the eigenstrain approach. A custom finite element method has been developed to solve the coupled system of equations of dislocation kinetics and crystal mechanics. The results show that, in general, dislocations self-organize in patterns under their mutual interactions. However, the famous dislocation cell structure has been found to form only when cross slip is implemented in the model. Cross slip is also found to lower the yield point, increase the hardening rate, and sustain an increase in the dislocation density over the hardening regime. Analysis of the cell structure evolution reveals that the average cell size decreases with the applied stress, which is consistent with the similitude principle.

Ion temperature effects on magnetotail Alfvén wave propagation and electron energization: ION TEMPERATURE EFFECTS ON ALFVÉN WAVES

DOE Office of Scientific and Technical Information (OSTI.GOV)

Damiano, P. A.; Johnson, J. R.; Chaston, C. C.

2015-07-01

A new 2-D self-consistent hybrid gyrofluid-kinetic electron model in dipolar coordinates is presented and used to simulate dispersive-scale Alfvén wave pulse propagation from the equator to the ionosphere along an L = 10 magnetic field line. The model is an extension of the hybrid MHD-kinetic electron model that incorporates ion Larmor radius corrections via the kinetic fluid model of Cheng and Johnson (1999). It is found that consideration of a realistic ion to electron temperature ratio decreases the propagation time of the wave from the plasma sheet to the ionosphere by several seconds relative to a ρi=0 case (which alsomore » implies shorter timing for a substorm onset signal) and leads to significant dispersion of wave energy perpendicular to the ambient magnetic field. Additionally, ion temperature effects reduce the parallel current and electron energization all along the field line for the same magnitude perpendicular electric field perturbation.« less

Simulations of anti-parallel reconnection using a nonlocal heat flux closure

DOE PAGES

Ng, Jonathan; Hakim, Ammar; Bhattacharjee, A.; ...

2017-08-08

The integration of kinetic effects in fluid models is important for global simulations of the Earth's magnetosphere. In particular, it has been shown that ion kinetics play a crucial role in the dynamics of large reconnecting systems, and that higher-order fluid moment models can account for some of these effects. Here, we use a ten-moment model for electrons and ions, which includes the off diagonal elements of the pressure tensor that are important for magnetic reconnection. Kinetic effects are recovered by using a nonlocal heat flux closure, which approximates linear Landau damping in the fluid framework. Moreover, the closure ismore » tested using the island coalescence problem, which is sensitive to ion dynamics. We also demonstrate that the nonlocal closure is able to self-consistently reproduce the structure of the ion diffusion region, pressure tensor, and ion velocity without the need for fine-tuning of relaxation coefficients present in earlier models.« less

SUPRATHERMAL SOLAR WIND ELECTRONS AND LANGMUIR TURBULENCE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kim, Sunjung; Yoon, Peter H.; Choe, G. S.

2016-09-01

The steady-state model recently put forth for the solar wind electron velocity distribution function during quiet time conditions, was originally composed of three population electrons (core, halo, and superhalo) with the core remaining nonresonant with any plasma waves while the halo and superhalo separately maintained steady-state resonance with whistler- and Langmuir-frequency range fluctuations, respectively. However, a recent paper demonstrates that whistler-range fluctuations in fact have no significant contribution. The present paper represents a consummation of the model in that a self-consistent model of the suprathermal electron population, which encompasses both the halo and the superhalo, is constructed solely on themore » basis of the Langmuir fluctuation spectrum. Numerical solutions to steady-state particle and wave kinetic equations are obtained on the basis of an initial trial electron distribution and Langmuir wave spectrum. Such a finding offers a self-consistent explanation for the observed steady-state electron distribution in the solar wind.« less

The nonlocal electron kinetics for a low-pressure glow discharge dusty plasma

NASA Astrophysics Data System (ADS)

Liang, Yonggan; Wang, Ying; Li, Hui; Tian, Ruihuan; Yuan, Chengxun; Kudryavtsev, A. A.; Rabadanov, K. M.; Wu, Jian; Zhou, Zhongxiang; Tian, Hao

2018-05-01

The nonlocal electron kinetic model based on the Boltzmann equation is developed in low-pressure argon glow discharge dusty plasmas. The additional electron-dust elastic and inelastic collision processes are considered when solving the kinetic equation numerically. The orbital motion limited theory and collision enhanced collection approximation are employed to calculate the dust surface potential. The electron energy distribution function (EEDF), effective electron temperature Teff, and dust surface potential are investigated under different plasma and dust conditions by solving the Boltzmann and the dust charging current balance equations self-consistently. A comparison of the calculation results obtained from nonlocal and local kinetic models is made. It is shown that the appearance of dust particles leads to a deviation of the EEDF from its original profile for both nonlocal and local kinetic models. With the increase in dust density and size, the effective electron temperature and dust surface potential decrease due to the high-energy electron loss on the dust surface. Meanwhile, the nonlocal and local results differ much from each other under the same calculation condition. It is concluded that, for low-pressure (PR ≤ 1 cm*Torr) glow discharge dusty plasmas, the existence of dust particles will amplify the difference of local and nonlocal EEDFs, which makes the local kinetic model more improper to determine the main parameters of the positive column. The nonlocal kinetic model should be used for the calculation of the EEDFs and dusty plasma parameters.

Self-Consistent Model of Magnetospheric Ring Current and Electromagnetic Ion Cyclotron Waves: The 2-7 May 1998 Storm

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Gamayunov, K. V.; Jordanova, V. K.

2003-01-01

A complete description of a self-consistent model of magnetospheric ring current interacting with electromagnetic ion cyclotron waves is presented. The model is based on the system of two kinetic equations; one equation describes the ring current ion dynamics, and another equation describes the wave evolution. The effects on ring current ions interacting with electromagnetic ion cyclotron waves and back on waves are considered self-consistently by solving both equations on a global magnetospheric scale under nonsteady state conditions. The developed model is employed to simulate the entire 2-7 May 1998 storm period. First, the trapped number fluxes of the ring current protons are calculated and presented along with comparison with the data measured by the three- dimensional hot plasma instrument Polar/HYDRA. Incorporating in the model the wave-particle interaction leads to much better agreement between the experimental data and the model results. Second, examining of the wave (MLT, L shell) distributions produced by the model during the storm progress reveals an essential intensification of the wave emission about 2 days after the main phase of the storm. This result is well consistent with the earlier ground-based observations. Finally, the theoretical shapes and the occurrence rates of the wave power spectral densities are studied. It is found that about 2 days after the storm s main phase on 4 May, mainly non-Gaussian shapes of power spectral densities are produced.

Controlled release of functional proteins through designer self-assembling peptide nanofiber hydrogel scaffold

PubMed Central

Koutsopoulos, Sotirios; Unsworth, Larry D.; Nagai, Yusuke; Zhang, Shuguang

2009-01-01

The release kinetics for a variety of proteins of a wide range of molecular mass, hydrodynamic radii, and isoelectric points through a nanofiber hydrogel scaffold consisting of designer self-assembling peptides were studied by using single-molecule fluorescence correlation spectroscopy (FCS). In contrast to classical diffusion experiments, the single-molecule approach allowed for the direct determination of diffusion coefficients for lysozyme, trypsin inhibitor, BSA, and IgG both inside the hydrogel and after being released into the solution. The results of the FCS analyses and the calculated pristine in-gel diffusion coefficients were compared with the values obtained from the Stokes–Einstein equation, Fickian diffusion models, and the literature. The release kinetics suggested that protein diffusion through nanofiber hydrogels depended primarily on the size of the protein. Protein diffusivities decreased, with increasing hydrogel nanofiber density providing a means of controlling the release kinetics. Secondary and tertiary structure analyses and biological assays of the released proteins showed that encapsulation and release did not affect the protein conformation and functionality. Our results show that this biocompatible and injectable designer self-assembling peptide hydrogel system may be useful as a carrier for therapeutic proteins for sustained release applications. PMID:19273853

Modeling Electric Field Influences on Plasmaspheric Refilling

NASA Technical Reports Server (NTRS)

Liemohn, M. W.; Kozyra, J. U.; Khazanov, G. V.; Craven, Paul D.

1998-01-01

We have a new model of ion transport that we have applied to the problem of plasmaspheric flux tube refilling after a geomagnetic disturbance. This model solves the Fokker-Planck kinetic equation by applying discrete difference numerical schemes to the various operators. Features of the model include a time-varying ionospheric source, self-consistent Coulomb collisions, field-aligned electric field, hot plasma interactions, and ion cyclotron wave heating. We see refilling rates similar to those of earlier observations and models, except when the electric field is included. In this case, the refilling rates can be quite different that previously predicted. Depending on the populations included and the values of relevant parameters, trap zone densities can increase or decrease. In particular, the inclusion of hot populations near the equatorial region (specifically warm pancake distributions and ring current ions) can dramatically alter the refilling rate. Results are compared with observations as well as previous hydrodynamic and kinetic particle model simulations.

Hydrodynamic Models for Multicomponent Plasmas with Collisional-Radiative Kinetics

DTIC Science & Technology

2014-12-01

16, 17]. The plasma, typically created by electric discharges , can deposit heat locally in the vicinity the flame, which quickly raises the gas...the corona layer of laser produced plasmas (LPP). Secondly, the self-consistent coupling of the plasma with the field gives rise to particle...excited species and reaction radicals; 7 n ncr solid transport layer (overdense) corona layer (underdense) temperature density shock wave ablation

Development of Self-Potential Tomography for Early Warning System of rainfall induced Landslides: Electro-kinetic Effects and Sandbox Experiments

NASA Astrophysics Data System (ADS)

Hattori, K.; Yamazaki, T.; Terajima, T.; Huang, Q.

2017-12-01

Electro-kinetic effects is one of the possible mechanism for ULF electromagnetic phenomena preceeding landlisdes and large earthquakes. To understand general tendencies of electromagnetic changes related to electro-kinetic effects, we struggle with the integrated research to clarify the coupling among hydrological, geotechnical, and electromagnetic changes. Our final goal is to understand the ULF elecromagnetic phenomena in order to develop a simple technology for earthquake monitoring/forecasting. So, in this paper, we first show the observed waveforms possibly related to the Boso slow slip events in 2003 and 2007 and indoor landslide experiments with artificial rainfall. Then, numerical computations on the self-potential variation by the simulated groundwater flow, and compare the results with those observed by laboratory experiments. In the result, the simulated self-potential variation is consistent with observed one. FInally, we developed self-potential tomography to estimate the ground water condition. And we also characterize the pressure from the self-potential data, and compare the result with observed pressure head that is measured by pore-pressure gauge and found that the inverted pressure head is consistent with observed one. In addition, we apply the self-potential data observed by the flume test. The estimated pressure head from observed self-potential data shows the consistency with observed pressure head. And estimated pressure head also show the characteristic distribution before the landslide occurred. These facts are highly suggestive in effectiveness of the self-potential tomography to monitor groundwater changes associated with landslide. The details will be given in our presentation.

3D Hall MHD-EPIC Simulations of Ganymede's Magnetosphere

NASA Astrophysics Data System (ADS)

Zhou, H.; Toth, G.; Jia, X.

2017-12-01

Fully kinetic modeling of a complete 3D magnetosphere is still computationally expensive and not feasible on current computers. While magnetohydrodynamic (MHD) models have been successfully applied to a wide range of plasma simulation, they cannot capture some important kinetic effects. We have recently developed a new modeling tool to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme (BATS-R-US) magnetohydrodynamic model. This results in a kinetic model of the regions where kinetic effects are important. In addition to the MHD-EPIC modeling of the magnetosphere, the improved model presented here is now able to represent the moon as a resistive body. We use a stretched spherical grid with adaptive mesh refinement (AMR) to capture the resistive body and its boundary. A semi-implicit scheme is employed for solving the magnetic induction equation to allow time steps that are not limited by the resistivity. We have applied the model to Ganymede, the only moon in the solar system known to possess a strong intrinsic magnetic field, and included finite resistivity beneath the moon`s surface to model the electrical properties of the interior in a self-consistent manner. The kinetic effects of electrons and ions on the dayside magnetopause and tail current sheet are captured with iPIC3D. Magnetic reconnections under different upstream background conditions of several Galileo flybys are simulated to study the global reconnection rate and the magnetospheric dynamics

The tight binding model study of the role of anisotropic AFM spin ordering in the charge ordered CMR manganites

NASA Astrophysics Data System (ADS)

Kar, J. K.; Panda, Saswati; Rout, G. C.

2017-05-01

We propose here a tight binding model study of the interplay between charge and spin orderings in the CMR manganites taking anisotropic effect due to electron hoppings and spin exchanges. The Hamiltonian consists of the kinetic energies of eg and t2g electrons of manganese ion. It further includes double exchange and Heisenberg interactions. The charge density wave interaction (CDW) describes an extra mechanism for the insulating character of the system. The CDW gap and spin parameters are calculated using Zubarev's Green's function technique and computed self-consistently. The results are reported in this communication.

Theoretical Technology Research for the International Solar Terrestrial Physics (ISTP) Program

NASA Technical Reports Server (NTRS)

Ashour-Abdalla, Maha; Curtis, Steve (Technical Monitor)

2002-01-01

During the last four years the UCLA (University of California, Los Angeles) IGPP (Institute of Geophysics and Planetary Physics) Space Plasma Simulation Group has continued its theoretical effort to develop a Mission Oriented Theory (MOT) for the International Solar Terrestrial Physics (ISTP) program. This effort has been based on a combination of approaches: analytical theory, large-scale kinetic (LSK) calculations, global magnetohydrodynamic (MHD) simulations and self-consistent plasma kinetic (SCK) simulations. These models have been used to formulate a global interpretation of local measurements made by the ISTP spacecraft. The regions of applications of the MOT cover most of the magnetosphere: solar wind, low- and high- latitude magnetospheric boundary, near-Earth and distant magnetotail, and auroral region. Most recent investigations include: plasma processes in the electron foreshock, response of the magnetospheric cusp, particle entry in the magnetosphere, sources of observed distribution functions in the magnetotail, transport of oxygen ions, self-consistent evolution of the magnetotail, substorm studies, effects of explosive reconnection, and auroral acceleration simulations. A complete list of the activities completed under the grant follow.

Role of stresses in annealing of ion-implantation damage in Si

DOE Office of Scientific and Technical Information (OSTI.GOV)

Seshan, K.; EerNisse, E.P.

Recent results showing a crystallographic orientation dependence of growth kinetics, secondary defects, and stress relief in annealing of ion-implanted Si are shown to be self-consistent if interpreted in terms of the influence of stresses upon annealing processes. The stress influence proposed is microplastic shear which is induced in (112) directions on (111) planes inclined to the implant surface by the biaxial stress created in the implant region by ion-implantation damage. The shear stresses are shown to be dependent on crystallographic orientation in a manner consistent with the model.

Numerical methods for solving moment equations in kinetic theory of neuronal network dynamics

NASA Astrophysics Data System (ADS)

Rangan, Aaditya V.; Cai, David; Tao, Louis

2007-02-01

Recently developed kinetic theory and related closures for neuronal network dynamics have been demonstrated to be a powerful theoretical framework for investigating coarse-grained dynamical properties of neuronal networks. The moment equations arising from the kinetic theory are a system of (1 + 1)-dimensional nonlinear partial differential equations (PDE) on a bounded domain with nonlinear boundary conditions. The PDEs themselves are self-consistently specified by parameters which are functions of the boundary values of the solution. The moment equations can be stiff in space and time. Numerical methods are presented here for efficiently and accurately solving these moment equations. The essential ingredients in our numerical methods include: (i) the system is discretized in time with an implicit Euler method within a spectral deferred correction framework, therefore, the PDEs of the kinetic theory are reduced to a sequence, in time, of boundary value problems (BVPs) with nonlinear boundary conditions; (ii) a set of auxiliary parameters is introduced to recast the original BVP with nonlinear boundary conditions as BVPs with linear boundary conditions - with additional algebraic constraints on the auxiliary parameters; (iii) a careful combination of two Newton's iterates for the nonlinear BVP with linear boundary condition, interlaced with a Newton's iterate for solving the associated algebraic constraints is constructed to achieve quadratic convergence for obtaining the solutions with self-consistent parameters. It is shown that a simple fixed-point iteration can only achieve a linear convergence for the self-consistent parameters. The practicability and efficiency of our numerical methods for solving the moment equations of the kinetic theory are illustrated with numerical examples. It is further demonstrated that the moment equations derived from the kinetic theory of neuronal network dynamics can very well capture the coarse-grained dynamical properties of integrate-and-fire neuronal networks.

Kinetic Simulations of Type II Radio Burst Emission Processes

NASA Astrophysics Data System (ADS)

Ganse, U.; Spanier, F. A.; Vainio, R. O.

2011-12-01

The fundamental emission process of Type II Radio Bursts has been under discussion for many decades. While analytic deliberations point to three wave interaction as the source for fundamental and harmonic radio emissions, sparse in-situ observational data and high computational demands for kinetic simulations have not allowed for a definite conclusion to be reached. A popular model puts the radio emission into the foreshock region of a coronal mass ejection's shock front, where shock drift acceleration can create eletrcon beam populations in the otherwise quiescent foreshock plasma. Beam-driven instabilities are then assumed to create waves, forming the starting point of three wave interaction processes. Using our kinetic particle-in-cell code, we have studied a number of emission scenarios based on electron beam populations in a CME foreshock, with focus on wave-interaction microphysics on kinetic scales. The self-consistent, fully kinetic simulations with completely physical mass-ratio show fundamental and harmonic emission of transverse electromagnetic waves and allow for detailled statistical analysis of all contributing wavemodes and their couplings.

Validation of a coupled core-transport, pedestal-structure, current-profile and equilibrium model

NASA Astrophysics Data System (ADS)

Meneghini, O.

2015-11-01

The first workflow capable of predicting the self-consistent solution to the coupled core-transport, pedestal structure, and equilibrium problems from first-principles and its experimental tests are presented. Validation with DIII-D discharges in high confinement regimes shows that the workflow is capable of robustly predicting the kinetic profiles from on axis to the separatrix and matching the experimental measurements to within their uncertainty, with no prior knowledge of the pedestal height nor of any measurement of the temperature or pressure. Self-consistent coupling has proven to be essential to match the experimental results, and capture the non-linear physics that governs the core and pedestal solutions. In particular, clear stabilization of the pedestal peeling ballooning instabilities by the global Shafranov shift and destabilization by additional edge bootstrap current, and subsequent effect on the core plasma profiles, have been clearly observed and documented. In our model, self-consistency is achieved by iterating between the TGYRO core transport solver (with NEO and TGLF for neoclassical and turbulent flux), and the pedestal structure predicted by the EPED model. A self-consistent equilibrium is calculated by EFIT, while the ONETWO transport package evolves the current profile and calculates the particle and energy sources. The capabilities of such workflow are shown to be critical for the design of future experiments such as ITER and FNSF, which operate in a regime where the equilibrium, the pedestal, and the core transport problems are strongly coupled, and for which none of these quantities can be assumed to be known. Self-consistent core-pedestal predictions for ITER, as well as initial optimizations, will be presented. Supported by the US Department of Energy under DE-FC02-04ER54698, DE-SC0012652.

Time-dependent Tonks-Langmuir model is unstable

NASA Astrophysics Data System (ADS)

Sheridan, T. E.; Baalrud, S. D.

2017-11-01

We investigate a time-dependent extension of the Tonks-Langmuir model for a one-dimensional plasma discharge with collisionless kinetic ions and Boltzmann electrons. Ions are created uniformly throughout the volume and flow from the center of the discharge to the boundary wall due to a self-consistent, zero-order electric field. Solving this model using a particle-in-cell simulation, we observe coherent low-frequency, long-wavelength unstable ion waves which move toward the boundary with a speed below both the ion acoustic speed and the average ion velocity. The maximum amplitude of the wave potential fluctuations peaks at ≈0.09 Te near the wall, where Te is the electron temperature in electron volts. Using linear kinetic theory, we identify this instability as slow ion-acoustic wave modes which are destabilized by the zero-order electric field.

Ring Current Ion Coupling with Electromagnetic Ion Cyclotron Waves

NASA Technical Reports Server (NTRS)

Khazanov, George V.

2002-01-01

A new ring current global model has been developed for the first time that couples the system of two kinetic equations: one equation describes the ring current (RC) ion dynamic, and another equation describes wave evolution of electromagnetic ion cyclotron waves (EMIC). The coupled model is able to simulate, for the first time self-consistently calculated RC ion kinetic and evolution of EMIC waves that propagate along geomagnetic field lines and reflect from the ionosphere. Ionospheric properties affect the reflection index through the integral Pedersen and Hall coductivities. The structure and dynamics of the ring current proton precipitating flux regions, intensities of EMIC, global RC energy balance, and some other parameters will be studied in detail for the selected geomagnetic storms. The space whether aspects of RC modelling and comparison with the data will also be discussed.

Mechanisms of kinetic trapping in self-assembly and phase transformation

PubMed Central

Hagan, Michael F.; Elrad, Oren M.; Jack, Robert L.

2011-01-01

In self-assembly processes, kinetic trapping effects often hinder the formation of thermodynamically stable ordered states. In a model of viral capsid assembly and in the phase transformation of a lattice gas, we show how simulations in a self-assembling steady state can be used to identify two distinct mechanisms of kinetic trapping. We argue that one of these mechanisms can be adequately captured by kinetic rate equations, while the other involves a breakdown of theories that rely on cluster size as a reaction coordinate. We discuss how these observations might be useful in designing and optimising self-assembly reactions. PMID:21932884

Self-consistent atmosphere modeling with cloud formation for low-mass stars and exoplanets

NASA Astrophysics Data System (ADS)

Juncher, Diana; Jørgensen, Uffe G.; Helling, Christiane

2017-12-01

Context. Low-mass stars and extrasolar planets have ultra-cool atmospheres where a rich chemistry occurs and clouds form. The increasing amount of spectroscopic observations for extrasolar planets requires self-consistent model atmosphere simulations to consistently include the formation processes that determine cloud formation and their feedback onto the atmosphere. Aims: Our aim is to complement the MARCS model atmosphere suit with simulations applicable to low-mass stars and exoplanets in preparation of E-ELT, JWST, PLATO and other upcoming facilities. Methods: The MARCS code calculates stellar atmosphere models, providing self-consistent solutions of the radiative transfer and the atmospheric structure and chemistry. We combine MARCS with a kinetic model that describes cloud formation in ultra-cool atmospheres (seed formation, growth/evaporation, gravitational settling, convective mixing, element depletion). Results: We present a small grid of self-consistently calculated atmosphere models for Teff = 2000-3000 K with solar initial abundances and log (g) = 4.5. Cloud formation in stellar and sub-stellar atmospheres appears for Teff < 2700 K and has a significant effect on the structure and the spectrum of the atmosphere for Teff < 2400 K. We have compared the synthetic spectra of our models with observed spectra and found that they fit the spectra of mid- to late-type M-dwarfs and early-type L-dwarfs well. The geometrical extension of the atmospheres (at τ = 1) changes with wavelength resulting in a flux variation of 10%. This translates into a change in geometrical extension of the atmosphere of about 50 km, which is the quantitative basis for exoplanetary transit spectroscopy. We also test DRIFT-MARCS for an example exoplanet and demonstrate that our simulations reproduce the Spitzer observations for WASP-19b rather well for Teff = 2600 K, log (g) = 3.2 and solar abundances. Our model points at an exoplanet with a deep cloud-free atmosphere with a substantial day-night energy transport and no temperature inversion.

A dynamical model of plasma turbulence in the solar wind

PubMed Central

Howes, G. G.

2015-01-01

A dynamical approach, rather than the usual statistical approach, is taken to explore the physical mechanisms underlying the nonlinear transfer of energy, the damping of the turbulent fluctuations, and the development of coherent structures in kinetic plasma turbulence. It is argued that the linear and nonlinear dynamics of Alfvén waves are responsible, at a very fundamental level, for some of the key qualitative features of plasma turbulence that distinguish it from hydrodynamic turbulence, including the anisotropic cascade of energy and the development of current sheets at small scales. The first dynamical model of kinetic turbulence in the weakly collisional solar wind plasma that combines self-consistently the physics of Alfvén waves with the development of small-scale current sheets is presented and its physical implications are discussed. This model leads to a simplified perspective on the nature of turbulence in a weakly collisional plasma: the nonlinear interactions responsible for the turbulent cascade of energy and the formation of current sheets are essentially fluid in nature, while the collisionless damping of the turbulent fluctuations and the energy injection by kinetic instabilities are essentially kinetic in nature. PMID:25848075

Conservative algorithms for non-Maxwellian plasma kinetics

DOE PAGES

Le, Hai P.; Cambier, Jean -Luc

2017-12-08

Here, we present a numerical model and a set of conservative algorithms for Non-Maxwellian plasma kinetics with inelastic collisions. These algorithms self-consistently solve for the time evolution of an isotropic electron energy distribution function interacting with an atomic state distribution function of an arbitrary number of levels through collisional excitation, deexcitation, as well as ionization and recombination. Electron-electron collisions, responsible for thermalization of the electron distribution, are also included in the model. The proposed algorithms guarantee mass/charge and energy conservation in a single step, and is applied to the case of non-uniform gridding of the energy axis in the phasemore » space of the electron distribution function. Numerical test cases are shown to demonstrate the accuracy of the method and its conservation properties.« less

Self-consistent perturbed equilibrium with neoclassical toroidal torque in tokamaks

DOE PAGES

Park, Jong-Kyu; Logan, Nikolas C.

2017-03-01

Toroidal torque is one of the most important consequences of non-axisymmetric fields in tokamaks. The well-known neoclassical toroidal viscosity (NTV) is due to the second-order toroidal force from anisotropic pressure tensor in the presence of these asymmetries. This work shows that the first-order toroidal force originating from the same anisotropic pressure tensor, despite having no flux surface average, can significantly modify the local perturbed force balance and thus must be included in perturbed equilibrium self-consistent with NTV. The force operator with an anisotropic pressure tensor is not self-adjoint when the NTV torque is finite and thus is solved directly formore » each component. This approach yields a modified, non-self-adjoint Euler-Lagrange equation that can be solved using a variety of common drift-kinetic models in generalized tokamak geometry. The resulting energy and torque integral provides a unique way to construct a torque response matrix, which contains all the information of self-consistent NTV torque profiles obtainable by applying non-axisymmetric fields to the plasma. This torque response matrix can then be used to systematically optimize non-axisymmetric field distributions for desired NTV profiles. Published by AIP Publishing.« less

Bayesian Calibration of Thermodynamic Databases and the Role of Kinetics

NASA Astrophysics Data System (ADS)

Wolf, A. S.; Ghiorso, M. S.

2017-12-01

Self-consistent thermodynamic databases of geologically relevant materials (like Berman, 1988; Holland and Powell, 1998, Stixrude & Lithgow-Bertelloni 2011) are crucial for simulating geological processes as well as interpreting rock samples from the field. These databases form the backbone of our understanding of how fluids and rocks interact at extreme planetary conditions. Considerable work is involved in their construction from experimental phase reaction data, as they must self-consistently describe the free energy surfaces (including relative offsets) of potentially hundreds of interacting phases. Standard database calibration methods typically utilize either linear programming or least squares regression. While both produce a viable model, they suffer from strong limitations on the training data (which must be filtered by hand), along with general ignorance of many of the sources of experimental uncertainty. We develop a new method for calibrating high P-T thermodynamic databases for use in geologic applications. The model is designed to handle pure solid endmember and free fluid phases and can be extended to include mixed solid solutions and melt phases. This new calibration effort utilizes Bayesian techniques to obtain optimal parameter values together with a full family of statistically acceptable models, summarized by the posterior. Unlike previous efforts, the Bayesian Logistic Uncertain Reaction (BLUR) model directly accounts for both measurement uncertainties and disequilibrium effects, by employing a kinetic reaction model whose parameters are empirically determined from the experiments themselves. Thus, along with the equilibrium free energy surfaces, we also provide rough estimates of the activation energies, entropies, and volumes for each reaction. As a first application, we demonstrate this new method on the three-phase aluminosilicate system, illustrating how it can produce superior estimates of the phase boundaries by incorporating constraints from all available data, while automatically handling variable data quality due to a combination of measurement errors and kinetic effects.

Nonequilibrium Saturation States and Fractional Kinetic Processes In The Turbulent Magnetotail

NASA Astrophysics Data System (ADS)

Milovanov, A. V.; Zelenyi, L. M.

Magnetotail regions with the considerably stretched and thinned magnetic lobe field offer a fertile playground for studying the fundamental properties of the self-organized turbulent systems. The focus of this report is on the turbulent nonequilibrium satu- ration states (NESS's) of the tail, where the plasma strongly couples with the self- organized magnetic and inductive electric fields. We advocate an unconventional de- scription of the NESS's, which brings together the ideas of fractal geometry, topology of manifolds, and fractional ("strange") kinetics. A self-consistent nonlinear fractional kinetic equation is proposed for the particle dynamics near the marginal NESS. We ar- gue that the inherent variability of the NESS is manifest in the low-frequency fluctu- ation spectrum f-1 often referred to as "flicker noise". The self-consistent plasma distribution function at the NESS is shown to reveal a power-law nonthermal tail (E) E-, where the slope 6 7 depends on the type of the spatiotem- poral correlations in the medium. Basic theoretical predictions are discussed against observations. This study was sponsored by the INTAS project 97-1612 and RFBR grants 00-02-17127 and 00-15-96631.

Three-dimensional drift kinetic response of high-β plasmas in the DIII-D tokamak.

PubMed

Wang, Z R; Lanctot, M J; Liu, Y Q; Park, J-K; Menard, J E

2015-04-10

A quantitative interpretation of the experimentally measured high-pressure plasma response to externally applied three-dimensional (3D) magnetic field perturbations, across the no-wall Troyon β limit, is achieved. The self-consistent inclusion of the drift kinetic effects in magnetohydrodynamic (MHD) modeling [Y. Q. Liu et al., Phys. Plasmas 15, 112503 (2008)] successfully resolves an outstanding issue of the ideal MHD model, which significantly overpredicts the plasma-induced field amplification near the no-wall limit, as compared to experiments. The model leads to quantitative agreement not only for the measured field amplitude and toroidal phase but also for the measured internal 3D displacement of the plasma. The results can be important to the prediction of the reliable plasma behavior in advanced fusion devices, such as ITER [K. Ikeda, Nucl. Fusion 47, S1 (2007)].

Kinetic studies of divertor heat fluxes in Alcator C-Mod

NASA Astrophysics Data System (ADS)

Pankin, A. Y.; Bateman, G.; Kritz, A. H.; Rafiq, T.; Park, G. Y.; Chang, C. S.; Brunner, D.; Hughes, J. W.; Labombard, B.; Terry, J.

2010-11-01

The kinetic XGC0 code [C.S. Chang et al, Phys. Plasmas 11 (2004) 2649] is used to model the H- mode pedestal and SOL regions in Alcator C-Mod discharges. The self-consistent simulations in this study include kinetic neoclassical physics and anomalous transport models along with the ExB flow shear effects. The heat fluxes on the divertor plates are computed and the fluxes to the outer plate are compared with experimental observations. The dynamics of the radial electric field near the separatrix and in the SOL region are computed with the XGC0 code, and the effect of the anomalous transport on the heat fluxes in the SOL region is investigated. In particular, the particle and thermal diffusivities obtained in the analysis mode are compared with predictions from the theory-based anomalous transport models such as MMM95 [G. Bateman et al, Phys. Plasmas 5 (1998) 1793] and DRIBM [T. Rafiq et al, to appear in Phys. Plasmas (2010)]. It is found that there is a notable pinch effect in the inner separatrix region. Possible physical mechanisms for the particle and thermal pinches are discussed.

Kinetic study on non-thermal volumetric plasma decay in the early afterglow of air discharge generated by a short pulse microwave or laser

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Wei, E-mail: yangwei861212@126.com; Zhou, Qianhong; Dong, Zhiwei

This paper reports a kinetic study on non-thermal plasma decay in the early afterglow of air discharge generated by short pulse microwave or laser. A global self-consistent model is based on the particle balance of complex plasma chemistry, electron energy equation, and gas thermal balance equation. Electron-ion Coulomb collision is included in the steady state Boltzmann equation solver to accurately describe the electron mobility and other transport coefficients. The model is used to simulate the afterglow of microsecond to nanosecond pulse microwave discharge in N{sub 2}, O{sub 2}, and air, as well as femtosecond laser filament discharge in dry andmore » humid air. The simulated results for electron density decay are in quantitative agreement with the available measured ones. The evolution of plasma decay under an external electric field is also investigated, and the effect of gas heating is considered. The underlying mechanism of plasma density decay is unveiled through the above kinetic modeling.« less

Self-Consistent Model of Magnetospheric Ring Current and Electromagnetic Ion Cyclotron Waves: The May 2-7, 1998, Storm

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Gamayunov, K. V.; Jordanova, V. K.

2003-01-01

Complete description of a self-consistent model for magnetospheric ring current interacting with electromagnetic ion cyclotron waves is presented. The model is based on the system of two kinetic equations; one equation describes the ring current ion dynamics, and another equation describes the wave evolution. The effects on ring current ions interacting with electromagnetic ion cyclotron waves, and back on waves, are considered self-consistently by solving both equations on a global magnetospheric scale under non steady-state conditions. In the paper by Khazanov et al. [2002] this self-consistent model has only been shortly outlined, and discussions of many the model related details have been omitted. For example, in present study for the first time a new algorithm for numerical finding of the resonant numbers for quasilinear wave-particle interaction is described, or it is demonstrated that in order to describe quasilinear interaction in a multi-ion thermal plasma correctly, both e and He(+) modes of electromagnetic ion cyclotron waves should be employed. The developed model is used to simulate the entire May 2-7, 1998 storm period. Trapped number fluxes of the ring current protons are calculated and presented along with their comparison with the data measured by the 3D hot plasma instrument Polar/HYDRA. Examining of the wave (MLT, L shell) distributions produced during the storm progress reveals an essential intensification of the wave emissions in about two days after main phase of storm. This result is well consistent with the earlier ground-based observations. Also the theoretical shapes and the occurrence rates for power spectral densities of electromagnetic ion cyclotron waves are studied. It is found that in about 2 days after the storm main phase on May 4, mainly non Gaussian shapes of power spectral densities are produced.

A practical nonlocal model for heat transport in magnetized laser plasmas

NASA Astrophysics Data System (ADS)

Nicolaï, Ph. D.; Feugeas, J.-L. A.; Schurtz, G. P.

2006-03-01

A model of nonlocal transport for multidimensional radiation magnetohydrodynamics codes is presented. In laser produced plasmas, it is now believed that the heat transport can be strongly modified by the nonlocal nature of the electron conduction. Other mechanisms, such as self-generated magnetic fields, may also affect the heat transport. The model described in this work, based on simplified Fokker-Planck equations aims at extending the model of G. Schurtz, Ph. Nicolaï, and M. Busquet [Phys. Plasmas 7, 4238 (2000)] to magnetized plasmas. A complete system of nonlocal equations is derived from kinetic equations with self-consistent electric and magnetic fields. These equations are analyzed and simplified in order to be implemented into large laser fusion codes and coupled to other relevant physics. The model is applied to two laser configurations that demonstrate the main features of the model and point out the nonlocal Righi-Leduc effect in a multidimensional case.

Nonequilibrium Self-Assembly of π-Conjugated Oligopeptides in Solution.

PubMed

Li, Bo; Li, Songsong; Zhou, Yuecheng; Ardoña, Herdeline Ann M; Valverde, Lawrence R; Wilson, William L; Tovar, John D; Schroeder, Charles M

2017-02-01

Supramolecular assembly is a powerful method that can be used to generate materials with well-defined structures across multiple length scales. Supramolecular assemblies consisting of biopolymer-synthetic polymer subunits are specifically known to exhibit exceptional structural and functional diversity as well as programmable control of noncovalent interactions through hydrogen bonding in biopolymer subunits. Despite recent progress, there is a need to control and quantitatively understand assembly under nonequilibrium conditions. In this work, we study the nonequilibrium self-assembly of π-conjugated synthetic oligopeptides using a combination of experiments and analytical modeling. By isolating an aqueous peptide solution droplet within an immiscible organic layer, the rate of peptide assembly in the aqueous solution can be controlled by tuning the transport rate of acid that is used to trigger assembly. Using this approach, peptides are guided to assemble under reaction-dominated and diffusion-dominated conditions, with results showing a transition from a diffusion-limited reaction front to spatially homogeneous assembly as the transport rate of acid decreases. Interestingly, our results show that the morphology of self-assembled peptide fibers is controlled by the assembly kinetics such that increasingly homogeneous structures of self-assembled synthetic oligopeptides were generally obtained using slower rates of assembly. We further developed an analytical reaction-diffusion model to describe oligopeptide assembly, and experimental results are compared to the reaction-diffusion model across a range of parameters. Overall, this work highlights the importance of molecular self-assembly under nonequilibrium conditions, specifically showing that oligopeptide assembly is governed by a delicate balance between reaction kinetics and transport processes.

Derivation of the density functional theory from the cluster expansion.

PubMed

Hsu, J Y

2003-09-26

The density functional theory is derived from a cluster expansion by truncating the higher-order correlations in one and only one term in the kinetic energy. The formulation allows self-consistent calculation of the exchange correlation effect without imposing additional assumptions to generalize the local density approximation. The pair correlation is described as a two-body collision of bound-state electrons, and modifies the electron- electron interaction energy as well as the kinetic energy. The theory admits excited states, and has no self-interaction energy.

Local rules simulation of the kinetics of virus capsid self-assembly.

PubMed

Schwartz, R; Shor, P W; Prevelige, P E; Berger, B

1998-12-01

A computer model is described for studying the kinetics of the self-assembly of icosahedral viral capsids. Solution of this problem is crucial to an understanding of the viral life cycle, which currently cannot be adequately addressed through laboratory techniques. The abstract simulation model employed to address this is based on the local rules theory of. Proc. Natl. Acad. Sci. USA. 91:7732-7736). It is shown that the principle of local rules, generalized with a model of kinetics and other extensions, can be used to simulate complicated problems in self-assembly. This approach allows for a computationally tractable molecular dynamics-like simulation of coat protein interactions while retaining many relevant features of capsid self-assembly. Three simple simulation experiments are presented to illustrate the use of this model. These show the dependence of growth and malformation rates on the energetics of binding interactions, the tolerance of errors in binding positions, and the concentration of subunits in the examples. These experiments demonstrate a tradeoff within the model between growth rate and fidelity of assembly for the three parameters. A detailed discussion of the computational model is also provided.

Interfacial mixing in high-energy-density matter with a multiphysics kinetic model

NASA Astrophysics Data System (ADS)

Haack, Jeffrey R.; Hauck, Cory D.; Murillo, Michael S.

2017-12-01

We have extended a recently developed multispecies, multitemperature Bhatnagar-Gross-Krook model [Haack et al., J. Stat. Phys. 168, 822 (2017), 10.1007/s10955-017-1824-9], to include multiphysics capabilities that enable modeling of a wider range of physical conditions. In terms of geometry, we have extended from the spatially homogeneous setting to one spatial dimension. In terms of the physics, we have included an atomic ionization model, accurate collision physics across coupling regimes, self-consistent electric fields, and degeneracy in the electronic screening. We apply the model to a warm dense matter scenario in which the ablator-fuel interface of an inertial confinement fusion target is heated, but for larger length and time scales and for much higher temperatures than can be simulated using molecular dynamics. Relative to molecular dynamics, the kinetic model greatly extends the temperature regime and the spatiotemporal scales over which we are able to model. In our numerical results we observe hydrogen from the ablator material jetting into the fuel during the early stages of the implosion and compare the relative size of various diffusion components (Fickean diffusion, electrodiffusion, and barodiffusion) that drive this process. We also examine kinetic effects, such as anisotropic distributions and velocity separation, in order to determine when this problem can be described with a hydrodynamic model.

Generalized Kinetic Description of Steady-State Collisionless Plasmas

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Liemohn, M. W.; Krivorutsky, E. N.

1997-01-01

We present a general solution to the collisionless Boltzmann (Vlasov) equation for a free-flowing plasma along a magnetic field line using Liouville's theorem, allowing for an arbitrary potential structure including non-monotonicities. The constraints of the existing collisionless kinetic transport models are explored, and the need for a more general approach to the problem of self- consistent potential energy calculations is described. Then a technique that handles an arbitrary potential energy distribution along the field line is presented and discussed. For precipitation of magnetospherically trapped hot plasma, this model yields moment calculations that vary by up to a factor of two for various potential energy structures with the same total potential drop. The differences are much greater for the high-latitude outflow scenario, giving order of magnitude variations depending on the shape of the potential energy distribution.

2-3D nonlocal transport model in magnetized laser plasmas.

NASA Astrophysics Data System (ADS)

Nicolaï, Philippe; Feugeas, Jean-Luc; Schurtz, Guy

2004-11-01

We present a model of nonlocal transport for multidimensional radiation magneto-hydrodynamics codes. This model, based on simplified Fokker-Planck equations, aims at extending the formulae of G Schurtz,Ph.Nicolaï and M. Busquet [Phys. Plasmas,7,4238 (2000)] to magnetized plasmas.The improvements concern various points as the electric field effects on nonlocal transport or conversely the kinetic effects on E field. However the main purpose of this work is to generalize the previous model by including magnetic field effects. A complete system of nonlocal equations is derived from kinetic equations with self-consistent E and B fields. These equations are analyzed and simplified in order to be implemented into large laser fusion codes and coupled to other relevent physics. Finally, our model allows to obtain the deformation of the electron distribution function due to nonlocal effects. This deformation leads to a non-maxwellian function which could be used to compute the influence on other physical processes.

Self-Assembly of Human Serum Albumin: A Simplex Phenomenon

PubMed Central

Thakur, Garima; Prashanthi, Kovur; Jiang, Keren; Thundat, Thomas

2017-01-01

Spontaneous self-assemblies of biomolecules can generate geometrical patterns. Our findings provide an insight into the mechanism of self-assembled ring pattern generation by human serum albumin (HSA). The self-assembly is a process guided by kinetic and thermodynamic parameters. The generated protein ring patterns display a behavior which is geometrically related to a n-simplex model and is explained through thermodynamics and chemical kinetics. PMID:28930179

A normalized plot as a novel and time-saving tool in complex enzyme kinetic analysis.

PubMed

Bravo, I G; Busto, F; De Arriaga, D; Ferrero, M A; Rodríguez-Aparicio, L B; Martínez-Blanco, H; Reglero, A

2001-09-15

A new data treatment is described for designing kinetic experiments and analysing kinetic results for multi-substrate enzymes. Normalized velocities are plotted against normalized substrate concentrations. Data are grouped into n + 1 families across the range of substrate or product tested, n being the number of substrates plus products assayed. It has the following advantages over traditional methods: (1) it reduces to less than a half the amount of data necessary for a proper description of the system; (2) it introduces a self-consistency checking parameter that ensures the 'scientific reliability' of the mathematical output; (3) it eliminates the need for a prior knowledge of Vmax; (4) the normalization of data allows the use of robust and fuzzy methods suitable for managing really 'noisy' data; (5) it is appropriate for analysing complex systems, as the complete general equation is used, and the actual influence of effectors can be typified; (6) it is amenable to being implemented as a software that incorporates testing and electing among rival kinetic models.

Non-equilibrium dynamics of artificial quantum matter

NASA Astrophysics Data System (ADS)

Babadi, Mehrtash

The rapid progress of the field of ultracold atoms during the past two decades has set new milestones in our control over matter. By cooling dilute atomic gases and molecules to nano-Kelvin temperatures, novel quantum mechanical states of matter can be realized and studied on a table-top experimental setup while bulk matter can be tailored to faithfully simulate abstract theoretical models. Two of such models which have witnessed significant experimental and theoretical attention are (1) the two-component Fermi gas with resonant s-wave interactions, and (2) the single-component Fermi gas with dipole-dipole interactions. This thesis is devoted to studying the non-equilibrium collective dynamics of these systems using the general framework of quantum kinetic theory. We present a concise review of the utilized mathematical methods in the first two chapters, including the Schwinger-Keldysh formalism of non-equilibrium quantum fields, two-particle irreducible (2PI) effective actions and the framework of quantum kinetic theory. We study the collective dynamics of the dipolar Fermi gas in a quasi-two-dimensional optical trap in chapter 3 and provide a detailed account of its dynamical crossover from the collisionless to the hydrodynamical regime. Chapter 4 is devoted to studying the dynamics of the attractive Fermi gas in the normal phase. Starting from the self-consistent T-matrix (pairing fluctuation) approximation, we systematically derive a set of quantum kinetic equations and show that they provide a globally valid description of the dynamics of the attractive Fermi gas, ranging from the weak-coupling Fermi liquid phase to the intermediate non-Fermi liquid pairing pseudogap regime and finally the strong-coupling Bose liquid phase. The shortcomings of the self-consistent T-matrix approximation in two spatial dimensions are discussed along with a proposal to overcome its unphysical behaviors. The developed kinetic formalism is finally utilized to reproduce and interpret the findings of a recent experiment done on the collective dynamics of trapped two-dimensional ultracold gases.

Active matter beyond mean-field: ring-kinetic theory for self-propelled particles.

PubMed

Chou, Yen-Liang; Ihle, Thomas

2015-02-01

Recently, Hanke et al. [Phys. Rev. E 88, 052309 (2013)] showed that mean-field kinetic theory fails to describe collective motion in soft active colloids and that correlations must not be neglected. Correlation effects are also expected to be essential in systems of biofilaments driven by molecular motors and in swarms of midges. To obtain correlations in an active matter system from first principles, we derive a ring-kinetic theory for Vicsek-style models of self-propelled agents from the exact N-particle evolution equation in phase space. The theory goes beyond mean-field and does not rely on Boltzmann's approximation of molecular chaos. It can handle precollisional correlations and cluster formation, which are both important to understand the phase transition to collective motion. We propose a diagrammatic technique to perform a small-density expansion of the collision operator and derive the first two equations of the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. An algorithm is presented that numerically solves the evolution equation for the two-particle correlations on a lattice. Agent-based simulations are performed and informative quantities such as orientational and density correlation functions are compared with those obtained by ring-kinetic theory. Excellent quantitative agreement between simulations and theory is found at not-too-small noises and mean free paths. This shows that there are parameter ranges in Vicsek-like models where the correlated closure of the BBGKY hierarchy gives correct and nontrivial results. We calculate the dependence of the orientational correlations on distance in the disordered phase and find that it seems to be consistent with a power law with an exponent around -1.8, followed by an exponential decay. General limitations of the kinetic theory and its numerical solution are discussed.

Approximate models for the ion-kinetic regime in inertial-confinement-fusion capsule implosions

DOE PAGES

Hoffman, Nelson M.; Zimmerman, George B.; Molvig, Kim; ...

2015-05-19

“Reduced” (i.e., simplified or approximate) ion-kinetic (RIK) models in radiation-hydrodynamic simulations permit a useful description of inertial-confinement-fusion (ICF) implosions where kinetic deviations from hydrodynamic behavior are important. For implosions in or near the kinetic regime (i.e., when ion mean free paths are comparable to the capsule size), simulations using a RIK model give a detailed picture of the time- and space-dependent structure of imploding capsules, allow an assessment of the relative importance of various kinetic processes during the implosion, enable explanations of past and current observations, and permit predictions of the results of future experiments. The RIK simulation method describedmore » here uses moment-based reduced kinetic models for transport of mass, momentum, and energy by long-mean-free-path ions, a model for the decrease of fusion reactivity owing to the associated modification of the ion distribution function, and a model of hydrodynamic turbulent mixing. Transport models are based on local gradient-diffusion approximations for the transport of moments of the ion distribution functions, with coefficients to impose flux limiting or account for transport modification. After calibration against a reference set of ICF implosions spanning the hydrodynamic-to-kinetic transition, the method has useful, quantifiable predictive ability over a broad range of capsule parameter space. Calibrated RIK simulations show that an important contributor to ion species separation in ICF capsule implosions is the preferential flux of longer-mean-free-path species out of the fuel and into the shell, leaving the fuel relatively enriched in species with shorter mean free paths. Also, the transport of ion thermal energy is enhanced in the kinetic regime, causing the fuel region to have a more uniform, lower ion temperature, extending over a larger volume, than implied by clean simulations. Furthermore, we expect that the success of our simple approach will motivate continued theoretical research into the development of first-principles-based, comprehensive, self-consistent, yet useable models of kinetic multispecies ion behavior in ICF plasmas.« less

Development of an Uncertainty Quantification Predictive Chemical Reaction Model for Syngas Combustion

DOE PAGES

Slavinskaya, N. A.; Abbasi, M.; Starcke, J. H.; ...

2017-01-24

An automated data-centric infrastructure, Process Informatics Model (PrIMe), was applied to validation and optimization of a syngas combustion model. The Bound-to-Bound Data Collaboration (B2BDC) module of PrIMe was employed to discover the limits of parameter modifications based on uncertainty quantification (UQ) and consistency analysis of the model–data system and experimental data, including shock-tube ignition delay times and laminar flame speeds. Existing syngas reaction models are reviewed, and the selected kinetic data are described in detail. Empirical rules were developed and applied to evaluate the uncertainty bounds of the literature experimental data. Here, the initial H 2/CO reaction model, assembled frommore » 73 reactions and 17 species, was subjected to a B2BDC analysis. For this purpose, a dataset was constructed that included a total of 167 experimental targets and 55 active model parameters. Consistency analysis of the composed dataset revealed disagreement between models and data. Further analysis suggested that removing 45 experimental targets, 8 of which were self-inconsistent, would lead to a consistent dataset. This dataset was subjected to a correlation analysis, which highlights possible directions for parameter modification and model improvement. Additionally, several methods of parameter optimization were applied, some of them unique to the B2BDC framework. The optimized models demonstrated improved agreement with experiments compared to the initially assembled model, and their predictions for experiments not included in the initial dataset (i.e., a blind prediction) were investigated. The results demonstrate benefits of applying the B2BDC methodology for developing predictive kinetic models.« less

Development of an Uncertainty Quantification Predictive Chemical Reaction Model for Syngas Combustion

DOE Office of Scientific and Technical Information (OSTI.GOV)

Slavinskaya, N. A.; Abbasi, M.; Starcke, J. H.

An automated data-centric infrastructure, Process Informatics Model (PrIMe), was applied to validation and optimization of a syngas combustion model. The Bound-to-Bound Data Collaboration (B2BDC) module of PrIMe was employed to discover the limits of parameter modifications based on uncertainty quantification (UQ) and consistency analysis of the model–data system and experimental data, including shock-tube ignition delay times and laminar flame speeds. Existing syngas reaction models are reviewed, and the selected kinetic data are described in detail. Empirical rules were developed and applied to evaluate the uncertainty bounds of the literature experimental data. Here, the initial H 2/CO reaction model, assembled frommore » 73 reactions and 17 species, was subjected to a B2BDC analysis. For this purpose, a dataset was constructed that included a total of 167 experimental targets and 55 active model parameters. Consistency analysis of the composed dataset revealed disagreement between models and data. Further analysis suggested that removing 45 experimental targets, 8 of which were self-inconsistent, would lead to a consistent dataset. This dataset was subjected to a correlation analysis, which highlights possible directions for parameter modification and model improvement. Additionally, several methods of parameter optimization were applied, some of them unique to the B2BDC framework. The optimized models demonstrated improved agreement with experiments compared to the initially assembled model, and their predictions for experiments not included in the initial dataset (i.e., a blind prediction) were investigated. The results demonstrate benefits of applying the B2BDC methodology for developing predictive kinetic models.« less

Kinetic and Structural Evolution of Self-gravitating, Magnetized Clouds: 2.5-dimensional Simulations of Decaying Turbulence

NASA Astrophysics Data System (ADS)

Ostriker, Eve C.; Gammie, Charles F.; Stone, James M.

1999-03-01

The molecular component of the Galaxy is comprised of turbulent, magnetized clouds, many of which are self-gravitating and form stars. To develop an understanding of how these clouds' kinetic and structural evolution may depend on their level of turbulence, mean magnetization, and degree of self-gravity, we perform a survey of direct numerical MHD simulations in which three parameters are independently varied. Our simulations consist of solutions to the time-dependent MHD equations on a two-dimensional grid with periodic boundary conditions; an additional ``half'' dimension is also incorporated as dependent variables in the third Cartesian direction. Two of our survey parameters, the mean magnetization parameter β≡c2sound/v2Alfven and the Jeans number nJ≡Lcloud/LJeans, allow us to model clouds that either meet or fail conditions for magneto-Jeans stability and magnetic criticality. Our third survey parameter, the sonic Mach number M≡σvelocity/csound, allows us to initiate turbulence of either sub- or super-Alfvénic amplitude; we employ an isothermal equation of state throughout. We evaluate the times for each cloud model to become gravitationally bound and measure each model's kinetic energy loss over the fluid-flow crossing time. We compare the evolution of density and magnetic field structural morphology and quantify the differences in the density contrast generated by internal stresses for models of differing mean magnetization. We find that the values of β and nJ, but not the initial Mach number M, determine the time for cloud gravitational binding and collapse: for mean cloud density nH2=100 cm-3, unmagnetized models collapse after ~5 Myr, and magnetically supercritical models generally collapse after 5-10 Myr (although the smallest magneto-Jeans stable clouds survive gravitational collapse until t~15 Myr), while magnetically subcritical clouds remain uncollapsed over the entire simulations; these cloud collapse times scale with the mean density as tg~n-1/2H2. We find, contrary to some previous expectations, less than a factor of 2 difference between turbulent decay times for models with varying magnetic field strength; the maximum decay time, for B~14 μG and nH2=100 cm-3, is 1.4 flow crossing times tcross=L/σvelocity (or 8 Myr for typical giant molecular cloud parameters). In all models we find turbulent amplification in the magnetic field strength up to at least the level βpert≡c2sound/δv2Alfven=0.1, with the turbulent magnetic energy between 25% and 60% of the turbulent kinetic energy after one flow crossing time. We find that for non-self-gravitating stages of evolution, when clouds have M=5-10, the mass-averaged density contrast magnitudes are in the range 0.2-0.5, with the contrast increasing both toward low and high β. Although our conclusions about density statistics may be affected by our isothermal assumption, we note that only the more strongly magnetized models appear to be consistent with estimates of clump/interclump density contrasts inferred in Galactic giant molecular clouds.

Experimental Study and Mathematical Modeling of Self-Sustained Kinetic Oscillations in Catalytic Oxidation of Methane over Nickel.

PubMed

Lashina, Elena A; Kaichev, Vasily V; Saraev, Andrey A; Vinokurov, Zakhar S; Chumakova, Nataliya A; Chumakov, Gennadii A; Bukhtiyarov, Valerii I

2017-09-21

The self-sustained kinetic oscillations in the oxidation of CH 4 over Ni foil have been studied at atmospheric pressure using an X-ray diffraction technique and mass spectrometry. It has been shown that the regular oscillations appear under oxygen-deficient conditions; CO, CO 2 , H 2 , and H 2 O are detected as the products. According to in situ X-ray diffraction measurements, nickel periodically oxidizes to NiO initiating the reaction-rate oscillations. To describe the oscillations, we have proposed a five-stage mechanism of the partial oxidation of methane over Ni and a corresponding three-variable kinetic model. The mechanism considers catalytic methane decomposition, dissociative adsorption of oxygen, transformation of chemisorbed oxygen to surface nickel oxide, and reaction of adsorbed carbon and oxygen species to form CO. Analysis of the kinetic model indicates that the competition of two processes, i.e., the oxidation and the carbonization of the catalyst surface, is the driving force of the self-sustained oscillations in the oxidation of methane. We have compared this mechanism with the detailed 18-stage mechanism described previously by Lashina et al. (Kinetics and Catalysis 2012, 53, 374-383). It has been shown that both kinetic mechanisms coupled with a continuous stirred-tank reactor model describe well the oscillatory behavior in the oxidation of methane under non-isothermal conditions.

Merging for Particle-Mesh Complex Particle Kinetic Modeling of the Multiple Plasma Beams

NASA Technical Reports Server (NTRS)

Lipatov, Alexander S.

2011-01-01

We suggest a merging procedure for the Particle-Mesh Complex Particle Kinetic (PMCPK) method in case of inter-penetrating flow (multiple plasma beams). We examine the standard particle-in-cell (PIC) and the PMCPK methods in the case of particle acceleration by shock surfing for a wide range of the control numerical parameters. The plasma dynamics is described by a hybrid (particle-ion-fluid-electron) model. Note that one may need a mesh if modeling with the computation of an electromagnetic field. Our calculations use specified, time-independent electromagnetic fields for the shock, rather than self-consistently generated fields. While a particle-mesh method is a well-verified approach, the CPK method seems to be a good approach for multiscale modeling that includes multiple regions with various particle/fluid plasma behavior. However, the CPK method is still in need of a verification for studying the basic plasma phenomena: particle heating and acceleration by collisionless shocks, magnetic field reconnection, beam dynamics, etc.

Self consistent hydrodynamic description of the plasma wake field excitation induced by a relativistic charged-particle beam in an unmagnetized plasma

NASA Astrophysics Data System (ADS)

Jovanović, Dušan; Fedele, Renato; De Nicola, Sergio; Akhter, Tamina; Belić, Milivoj

2017-12-01

A self-consistent nonlinear hydrodynamic theory is presented of the propagation of a long and thin relativistic electron beam, for a typical plasma wake field acceleration configuration in an unmagnetized and overdense plasma. The random component of the trajectories of the beam particles as well as of their velocity spread is modelled by an anisotropic temperature, allowing the beam dynamics to be approximated as a 3D adiabatic expansion/compression. It is shown that even in the absence of the nonlinear plasma wake force, the localisation of the beam in the transverse direction can be achieved owing to the nonlinearity associated with the adiabatic compression/rarefaction and a coherent stationary state is constructed. Numerical calculations reveal the possibility of the beam focussing and defocussing, but the lifetime of the beam can be significantly extended by the appropriate adjustments, so that transverse oscillations are observed, similar to those predicted within the thermal wave and Vlasov kinetic models.

Kinetics of motility-induced phase separation and swim pressure

NASA Astrophysics Data System (ADS)

Patch, Adam; Yllanes, David; Marchetti, M. Cristina

Active Brownian particles (ABPs) represent a minimal model of active matter consisting of self-propelled spheres with purely repulsive interactions and rotational noise. We correlate the time evolution of the mean pressure towards its steady state value with the kinetics of motility-induced phase separation. For parameter values corresponding to phase separated steady states, we identify two dynamical regimes. The pressure grows monotonically in time during the initial regime of rapid cluster formation, overshooting its steady state value and then quickly relaxing to it, and remains constant during the subsequent slower period of cluster coalescence and coarsening. The overshoot is a distinctive feature of active systems. NSF-DMR-1305184, NSF-DGE-1068780, ACI-1341006, FIS2015-65078-C02, BIFI-ZCAM.

Calculation of kinetic rate constants from thermodynamic data

NASA Technical Reports Server (NTRS)

Marek, C. John

1995-01-01

A new scheme for relating the absolute value for the kinetic rate constant k to the thermodynamic constant Kp is developed for gases. In this report the forward and reverse rate constants are individually related to the thermodynamic data. The kinetic rate constants computed from thermodynamics compare well with the current kinetic rate constants. This method is self consistent and does not have extensive rules. It is first demonstrated and calibrated by computing the HBr reaction from H2 and Br2. This method then is used on other reactions.

Self-Consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves in Multi-Ion Magnetosphere

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Gamayunov, K. V.; Gallagher, D. L.; Kozyra, J. U.

2006-01-01

The further development of a self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves (Khazanov et al., 2003) is presented In order to adequately take into account wave propagation and refraction in a multi-ion magnetosphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation. This is a major new feature of the present model and, to the best of our knowledge, the ray tracing equations for the first time are explicitly employed on a global magnetospheric scale in order to self-consistently simulate the spatial, temporal, and spectral evolution of the ring current and of electromagnetic ion cyclotron waves To demonstrate the effects of EMIC wave propagation and refraction on the wave energy distribution and evolution, we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. First, owing to the density gradient at the plasmapause, the net wave refraction is suppressed, and He+-mode grows preferably at the plasmapause. This result is in total agreement with previous ray tracing studies and is very clearly found in presented B field spectrograms. Second, comparison of global wave distributions with the results from another ring current model (Kozyra et al., 1997) reveals that this new model provides more intense and more highly plasmapause-organized wave distributions during the May 1998 storm period Finally, it is found that He(+)-mode energy distributions are not Gaussian distributions and most important that wave energy can occupy not only the region of generation, i.e., the region of small wave normal angles, but all wave normal angles, including those to near 90 . The latter is extremely crucial for energy transfer to thermal plasmaspheric electrons by resonant Landau damping and subsequent downward heat transport and excitation of stable auroral red arcs.

Self-Consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves. 1; Waves in Multi Ion Magnetosphere

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Gumayunov, K. V.; Gallagher, D. L.; Kozyra, J. U.

2006-01-01

The further development of a self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves [Khazanov et al., 2003] is presented. In order to adequately take into account the wave propagation and refraction in a multi-ion plasmasphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation. This is a major new feature of the present model and, to the best of our knowledge, the ray tracing equations for the first time are explicitly employed on a global magnetospheric scale in order to self-consistently simulate spatial, temporal, and spectral evolutions of the ring current and electromagnetic ion cyclotron waves. To demonstrate the effects of EMIC wave propagation and refraction on the EMIC wave energy distributions and evolution we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. First, due to the density gradient at the plasmapause, the net wave refraction is suppressed, and He(+)-mode grows preferably at plasmapause. This result is in a total agreement with the previous ray tracing studies, and very clear observed in presented B-field spectrograms. Second, comparison the global wave distributions with the results from other ring current model [Kozyra et al., 1997] reveals that our model provides more intense and higher plasmapause organized distributions during the May, 1998 storm period. Finally, the found He(+)-mode energy distributions are not Gaussian distributions, and most important that wave energy can occupy not only the region of generation, i. e. the region of small wave normal angles, but the entire wave normal angle region and even only the region near 90 degrees. The latter is extremely crucial for energy transfer to thermal plasmaspheric electrons by resonant Landau damping, and subsequent downward heat transport and excitation of stable auroral red arcs.

Kinetic Modeling of Slow Energy Release in Non-Ideal Carbon Rich Explosives

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vitello, P; Fried, L; Glaesemann, K

2006-06-20

We present here the first self-consistent kinetic based model for long time-scale energy release in detonation waves in the non-ideal explosive LX-17. Non-ideal, insensitive carbon rich explosives, such as those based on TATB, are believed to have significant late-time slow release in energy. One proposed source of this energy is diffusion-limited growth of carbon clusters. In this paper we consider the late-time energy release problem in detonation waves using the thermochemical code CHEETAH linked to a multidimensional ALE hydrodynamics model. The linked CHEETAH-ALE model dimensional treats slowly reacting chemical species using kinetic rate laws, with chemical equilibrium assumed for speciesmore » coupled via fast time-scale reactions. In the model presented here we include separate rate equations for the transformation of the un-reacted explosive to product gases and for the growth of a small particulate form of condensed graphite to a large particulate form. The small particulate graphite is assumed to be in chemical equilibrium with the gaseous species allowing for coupling between the instantaneous thermodynamic state and the production of graphite clusters. For the explosive burn rate a pressure dependent rate law was used. Low pressure freezing of the gas species mass fractions was also included to account for regions where the kinetic coupling rates become longer than the hydrodynamic time-scales. The model rate parameters were calibrated using cylinder and rate-stick experimental data. Excellent long time agreement and size effect results were achieved.« less

Constitutive relations of ferroelectric ceramics

NASA Astrophysics Data System (ADS)

Su, Yu

The objective of this thesis is to obtain a better understanding on the fundamental constitutive behavior of ferroelectric ceramics based on the physics of phase transition, micromechanics of heterogeneous materials, and principles of irreversible thermodynamics. Within this framework, a self-consistent model is developed to investigate the electromechanical responses of ferroelectric polycrystals under temperature change and electromechanical loading. Cooling of a paraelectric crystal below its curie temperature Tc would result in spontaneous polarization, whereas electromechanical loading on a poled crystal could lead to domain switch. Domain growth and reorientation inside ferroelectric crystals are studied in light of these phase transition and domain switch. In this process, the change of the effective elastic, dielectric and piezoelectric constants during the evolution of microstructures are examined. In addition, hysteresis loops for the electric displacement and other related phenomena are computed under cyclic electric load. On top of all methods implemented in this work, the kinetic equation derived from the irreversible thermodynamics is the key to study the domain evolution in ferroelectric crystals. The kinetic relation not only governs the growth of new domain in a ferroelectric crystal, but it also determines the onset of phase transition. This characteristic is used to study the effect of hydrostatic pressure on the shift of Curie temperature of a ferroelectric crystal. Based on the derived expressions, it is observed that the deriving force can increase or decrease upon applied hydrostatic mechanical loading, depending on the change of electromechanical moduli, eigenstrain and electro-polarization. Several typical cases are computed and it is found that the change of the electromechanical moduli during phase transformation plays the key role in the shift of Curie temperature. Since ferroelectric ceramics are in a polycrystal form, a self-consistent model is used to examine the issues involved. In this model, each grain is represented by a spherical inclusion embedded in an infinitely extended piezoelectric matrix, and the inclusion further possesses an eigenstrain and eigen polarization. Secant relations between the polycrystal-matrix and the embedded inclusion are established by extending Hill's [1] incremental relations. An iterative computational program is developed for this self-consistent model.

Clustering and heterogeneous dynamics in a kinetic Monte Carlo model of self-propelled hard disks

NASA Astrophysics Data System (ADS)

Levis, Demian; Berthier, Ludovic

2014-06-01

We introduce a kinetic Monte Carlo model for self-propelled hard disks to capture with minimal ingredients the interplay between thermal fluctuations, excluded volume, and self-propulsion in large assemblies of active particles. We analyze in detail the resulting (density, self-propulsion) nonequilibrium phase diagram over a broad range of parameters. We find that purely repulsive hard disks spontaneously aggregate into fractal clusters as self-propulsion is increased and rationalize the evolution of the average cluster size by developing a kinetic model of reversible aggregation. As density is increased, the nonequilibrium clusters percolate to form a ramified structure reminiscent of a physical gel. We show that the addition of a finite amount of noise is needed to trigger a nonequilibrium phase separation, showing that demixing in active Brownian particles results from a delicate balance between noise, interparticle interactions, and self-propulsion. We show that self-propulsion has a profound influence on the dynamics of the active fluid. We find that the diffusion constant has a nonmonotonic behavior as self-propulsion is increased at finite density and that activity produces strong deviations from Fickian diffusion that persist over large time scales and length scales, suggesting that systems of active particles generically behave as dynamically heterogeneous systems.

Orbital-free bond breaking via machine learning

NASA Astrophysics Data System (ADS)

Snyder, John C.; Rupp, Matthias; Hansen, Katja; Blooston, Leo; Müller, Klaus-Robert; Burke, Kieron

2013-12-01

Using a one-dimensional model, we explore the ability of machine learning to approximate the non-interacting kinetic energy density functional of diatomics. This nonlinear interpolation between Kohn-Sham reference calculations can (i) accurately dissociate a diatomic, (ii) be systematically improved with increased reference data and (iii) generate accurate self-consistent densities via a projection method that avoids directions with no data. With relatively few densities, the error due to the interpolation is smaller than typical errors in standard exchange-correlation functionals.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Park, Jong-Kyu; Logan, Nikolas C.

Toroidal torque is one of the most important consequences of non-axisymmetric fields in tokamaks. The well-known neoclassical toroidal viscosity (NTV) is due to the second-order toroidal force from anisotropic pressure tensor in the presence of these asymmetries. This work shows that the first-order toroidal force originating from the same anisotropic pressure tensor, despite having no flux surface average, can significantly modify the local perturbed force balance and thus must be included in perturbed equilibrium self-consistent with NTV. The force operator with an anisotropic pressure tensor is not self-adjoint when the NTV torque is finite and thus is solved directly formore » each component. This approach yields a modified, non-self-adjoint Euler-Lagrange equation that can be solved using a variety of common drift-kinetic models in generalized tokamak geometry. The resulting energy and torque integral provides a unique way to construct a torque response matrix, which contains all the information of self-consistent NTV torque profiles obtainable by applying non-axisymmetric fields to the plasma. This torque response matrix can then be used to systematically optimize non-axisymmetric field distributions for desired NTV profiles. Published by AIP Publishing.« less

Density functional tight-binding and infrequent metadynamics can capture entropic effects in intramolecular hydrogen transfer reactions

NASA Astrophysics Data System (ADS)

Oliveira, Luiz F. L.; Fu, Christopher D.; Pfaendtner, Jim

2018-04-01

Infrequent metadynamics uses biased simulations to estimate the unbiased kinetics of a system, facilitating the calculation of rates and barriers. Here the method is applied to study intramolecular hydrogen transfer reactions involving peroxy radicals, a class of reactions that is challenging to model due to the entropic contributions of the formation of ring structures in the transition state. Using the self-consistent charge density-functional based tight-binding (DFTB) method, we applied infrequent metadynamics to the study of four intramolecular H-transfer reactions, demonstrating that the method can qualitatively reproduce these high entropic contributions, as observed in experiments and those predicted by transition state theory modeled by higher levels of theory. We also show that infrequent metadynamics and DFTB are successful in describing the relationship between transition state ring size and kinetic coefficients (e.g., activation energies and the pre-exponential factors).

Peptide self-assembly: thermodynamics and kinetics.

PubMed

Wang, Juan; Liu, Kai; Xing, Ruirui; Yan, Xuehai

2016-10-21

Self-assembling systems play a significant role in physiological functions and have therefore attracted tremendous attention due to their great potential for applications in energy, biomedicine and nanotechnology. Peptides, consisting of amino acids, are among the most popular building blocks and programmable molecular motifs. Nanostructures and materials assembled using peptides exhibit important potential for green-life new technology and biomedical applications mostly because of their bio-friendliness and reversibility. The formation of these ordered nanostructures pertains to the synergistic effect of various intermolecular non-covalent interactions, including hydrogen-bonding, π-π stacking, electrostatic, hydrophobic, and van der Waals interactions. Therefore, the self-assembly process is mainly driven by thermodynamics; however, kinetics is also a critical factor in structural modulation and function integration. In this review, we focus on the influence of thermodynamic and kinetic factors on structural assembly and regulation based on different types of peptide building blocks, including aromatic dipeptides, amphiphilic peptides, polypeptides, and amyloid-relevant peptides.

Thermodynamic and kinetic modelling of fuel oxidation behaviour in operating defective fuel

NASA Astrophysics Data System (ADS)

Lewis, operating defective fuel B. J.; Thompson, W. T.; Akbari, F.; Thompson, D. M.; Thurgood, C.; Higgs, J.

2004-07-01

A theoretical treatment has been developed to predict the fuel oxidation behaviour in operating defective nuclear fuel elements. The equilibrium stoichiometry deviation in the hyper-stoichiometric fuel has been derived from thermodynamic considerations using a self-consistent set of thermodynamic properties for the U-O system, which emphasizes replication of solubilities and three-phase invariant conditions displayed in the U-O binary phase diagram. The kinetics model accounts for multi-phase transport including interstitial oxygen diffusion in the solid and gas-phase transport of hydrogen and steam in the fuel cracks. The fuel oxidation model is further coupled to a heat conduction model to account for the feedback effect of a reduced thermal conductivity in the hyper-stoichiometric fuel. A numerical solution has been developed using a finite-element technique with the FEMLAB software package. The model has been compared to available data from several in-reactor X-2 loop experiments with defective fuel conducted at the Chalk River Laboratories. The model has also been benchmarked against an O/U profile measurement for a spent defective fuel element discharged from a commercial reactor.

Continuum kinetic and multi-fluid simulations of classical sheaths

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cagas, P.; Hakim, A.; Juno, J.

The kinetic study of plasma sheaths is critical, among other things, to understand the deposition of heat on walls, the effect of sputtering, and contamination of the plasma with detrimental impurities. The plasma sheath also provides a boundary condition and can often have a significant global impact on the bulk plasma. In this paper, kinetic studies of classical sheaths are performed with the continuum kinetic code, Gkeyll, which directly solves the Vlasov-Maxwell equations. The code uses a novel version of the finite-element discontinuous Galerkin scheme that conserves energy in the continuous-time limit. The fields are computed using Maxwell equations. Ionizationmore » and scattering collisions are included; however, surface effects are neglected. The aim of this work is to introduce the continuum kinetic method and compare its results with those obtained from an already established finite-volume multi-fluid model also implemented in Gkeyll. Novel boundary conditions on the fluids allow the sheath to form without specifying wall fluxes, so the fluids and fields adjust self-consistently at the wall. Our work demonstrates that the kinetic and fluid results are in agreement for the momentum flux, showing that in certain regimes, a multifluid model can be a useful approximation for simulating the plasma boundary. There are differences in the electrostatic potential between the fluid and kinetic results. Further, the direct solutions of the distribution function presented here highlight the non-Maxwellian distribution of electrons in the sheath, emphasizing the need for a kinetic model. The densities, velocities, and the potential show a good agreement between the kinetic and fluid results. But, kinetic physics is highlighted through higher moments such as parallel and perpendicular temperatures which provide significant differences from the fluid results in which the temperature is assumed to be isotropic. Besides decompression cooling, the heat flux is shown to play a role in the temperature differences that are observed, especially inside the collisionless sheath. Published by AIP Publishing.« less

Continuum kinetic and multi-fluid simulations of classical sheaths

DOE PAGES

Cagas, P.; Hakim, A.; Juno, J.; ...

2017-02-21

The kinetic study of plasma sheaths is critical, among other things, to understand the deposition of heat on walls, the effect of sputtering, and contamination of the plasma with detrimental impurities. The plasma sheath also provides a boundary condition and can often have a significant global impact on the bulk plasma. In this paper, kinetic studies of classical sheaths are performed with the continuum kinetic code, Gkeyll, which directly solves the Vlasov-Maxwell equations. The code uses a novel version of the finite-element discontinuous Galerkin scheme that conserves energy in the continuous-time limit. The fields are computed using Maxwell equations. Ionizationmore » and scattering collisions are included; however, surface effects are neglected. The aim of this work is to introduce the continuum kinetic method and compare its results with those obtained from an already established finite-volume multi-fluid model also implemented in Gkeyll. Novel boundary conditions on the fluids allow the sheath to form without specifying wall fluxes, so the fluids and fields adjust self-consistently at the wall. Our work demonstrates that the kinetic and fluid results are in agreement for the momentum flux, showing that in certain regimes, a multifluid model can be a useful approximation for simulating the plasma boundary. There are differences in the electrostatic potential between the fluid and kinetic results. Further, the direct solutions of the distribution function presented here highlight the non-Maxwellian distribution of electrons in the sheath, emphasizing the need for a kinetic model. The densities, velocities, and the potential show a good agreement between the kinetic and fluid results. But, kinetic physics is highlighted through higher moments such as parallel and perpendicular temperatures which provide significant differences from the fluid results in which the temperature is assumed to be isotropic. Besides decompression cooling, the heat flux is shown to play a role in the temperature differences that are observed, especially inside the collisionless sheath. Published by AIP Publishing.« less

The mechanisms for nanoparticle surface diffusion and chain self-assembly determined from real-time nanoscale kinetics in liquid

DOE PAGES

Woehl, Taylor J.; Prozorov, Tanya

2015-08-20

The mechanisms for nanoparticle self-assembly are often inferred from the morphology of the final nanostructures in terms of attractive and repulsive interparticle interactions. Understanding how nanoparticle building blocks are pieced together during self-assembly is a key missing component needed to unlock new strategies and mechanistic understanding of this process. Here we use real-time nanoscale kinetics derived from liquid cell transmission electron microscopy investigation of nanoparticle self-assembly to show that nanoparticle mobility dictates the pathway for self-assembly and final nanostructure morphology. We describe a new method for modulating nanoparticle diffusion in a liquid cell, which we employ to systematically investigate themore » effect of mobility on self-assembly of nanoparticles. We interpret the observed diffusion in terms of electrostatically induced surface diffusion resulting from nanoparticle hopping on the liquid cell window surface. Slow-moving nanoparticles self-assemble predominantly into linear 1D chains by sequential attachment of nanoparticles to existing chains, while highly mobile nanoparticles self-assemble into chains and branched structures by chain–chain attachments. Self-assembly kinetics are consistent with a diffusion-driven mechanism; we attribute the change in self-assembly pathway to the increased self-assembly rate of highly mobile nanoparticles. Furthermore, these results indicate that nanoparticle mobility can dictate the self-assembly mechanism and final nanostructure morphology in a manner similar to interparticle interactions.« less

Self-consistent simulation of high-frequency driven plasma sheaths

NASA Astrophysics Data System (ADS)

Shihab, Mohammed; Eremin, Denis; Mussenbrock, Thomas; Brinkmann, Ralf

2011-10-01

Low pressure capacitively coupled plasmas are widely used in plasma processing and microelectronics industry. Understanding the dynamics of the boundary sheath is a fundamental problem. It controls the energy and angular distribution of ions bombarding the electrode, which in turn affects the surface reaction rate and the profile of microscopic features. In this contribution, we investigate the dynamics of plasma boundary sheaths by means of a kinetic self-consistent model, which is able to resolve the ion dynamics. Asymmetric sheath dynamics is observed for the intermediate RF regime, i.e., in the regime where the ion plasma frequency is equal to the driving frequency. The ion inertia causes an additional phase difference between the expansion and the contraction phase of the plasma sheath and an asymmetry for the ion energy distribution bimodal shape. A comparison with experimental results and particle in cell simulations is performed. Low pressure capacitively coupled plasmas are widely used in plasma processing and microelectronics industry. Understanding the dynamics of the boundary sheath is a fundamental problem. It controls the energy and angular distribution of ions bombarding the electrode, which in turn affects the surface reaction rate and the profile of microscopic features. In this contribution, we investigate the dynamics of plasma boundary sheaths by means of a kinetic self-consistent model, which is able to resolve the ion dynamics. Asymmetric sheath dynamics is observed for the intermediate RF regime, i.e., in the regime where the ion plasma frequency is equal to the driving frequency. The ion inertia causes an additional phase difference between the expansion and the contraction phase of the plasma sheath and an asymmetry for the ion energy distribution bimodal shape. A comparison with experimental results and particle in cell simulations is performed. The financial support from the Federal Ministry of Education and Research within the frame of the project ``Plasma-Technology-Grid'' and the support of the DFG via the collaborative research center SFB-TR87 is gratefully acknowledged.

One-dimensional model of interacting-step fluctuations on vicinal surfaces: Analytical formulas and kinetic Monte-Carlo simulations

NASA Astrophysics Data System (ADS)

Patrone, Paul; Einstein, T. L.; Margetis, Dionisios

2011-03-01

We study a 1+1D, stochastic, Burton-Cabrera-Frank (BCF) model of interacting steps fluctuating on a vicinal crystal. The step energy accounts for entropic and nearest-neighbor elastic-dipole interactions. Our goal is to formulate and validate a self-consistent mean-field (MF) formalism to approximately solve the system of coupled, nonlinear stochastic differential equations (SDEs) governing fluctuations in surface motion. We derive formulas for the time-dependent terrace width distribution (TWD) and its steady-state limit. By comparison with kinetic Monte-Carlo simulations, we show that our MF formalism improves upon models in which step interactions are linearized. We also indicate how fitting parameters of our steady state MF TWD may be used to determine the mass transport regime and step interaction energy of certain experimental systems. PP and TLE supported by NSF MRSEC under Grant DMR 05-20471 at U. of Maryland; DM supported by NSF under Grant DMS 08-47587.

Modeling viscosity and diffusion of plasma mixtures across coupling regimes

NASA Astrophysics Data System (ADS)

Arnault, Philippe

2014-10-01

Viscosity and diffusion of plasma for pure elements and multicomponent mixtures are modeled from the high-temperature low-density weakly coupled regime to the low-temperature high-density strongly coupled regime. Thanks to an atom in jellium modeling, the effect of electron screening on the ion-ion interaction is incorporated through a self-consistent definition of the ionization. This defines an effective One Component Plasma, or an effective Binary Ionic Mixture, that is representative of the strength of the interaction. For the viscosity and the interdiffusion of mixtures, approximate kinetic expressions are supplemented by mixing laws applied to the excess viscosity and self-diffusion of pure elements. The comparisons with classical and quantum molecular dynamics results reveal deviations in the range 20--40% on average with almost no predictions further than a factor of 2 over many decades of variation. Applications in the inertial confinement fusion context could help in predicting the growth of hydrodynamic instabilities.

Elucidating dominant pathways of the nano-particle self-assembly process.

PubMed

Zeng, Xiangze; Li, Bin; Qiao, Qin; Zhu, Lizhe; Lu, Zhong-Yuan; Huang, Xuhui

2016-09-14

Self-assembly processes play a key role in the fabrication of functional nano-structures with widespread application in drug delivery and micro-reactors. In addition to the thermodynamics, the kinetics of the self-assembled nano-structures also play an important role in determining the formed structures. However, as the self-assembly process is often highly heterogeneous, systematic elucidation of the dominant kinetic pathways of self-assembly is challenging. Here, based on mass flow, we developed a new method for the construction of kinetic network models and applied it to identify the dominant kinetic pathways for the self-assembly of star-like block copolymers. We found that the dominant pathways are controlled by two competing kinetic parameters: the encounter time Te, characterizing the frequency of collision and the transition time Tt for the aggregate morphology change from rod to sphere. Interestingly, two distinct self-assembly mechanisms, diffusion of an individual copolymer into the aggregate core and membrane closure, both appear at different stages (with different values of Tt) of a single self-assembly process. In particular, the diffusion mechanism dominates the middle-sized semi-vesicle formation stage (with large Tt), while the membrane closure mechanism dominates the large-sized vesicle formation stage (with small Tt). Through the rational design of the hydrophibicity of the copolymer, we successfully tuned the transition time Tt and altered the dominant self-assembly pathways.

Electron kinetics in atmospheric-pressure argon and nitrogen microwave microdischarges

DOE Office of Scientific and Technical Information (OSTI.GOV)

Levko, Dmitry; Raja, Laxminarayan L.

2016-04-28

Electron kinetics in atmospheric-pressure argon and nitrogen microwave (4 GHz) microdischarges is studied using a self-consistent one-dimensional Particle-in-Cell Monte Carlo Collisions model. The reversal of electric field (i.e., inverted sheath formation) is obtained in nitrogen and is not obtained in argon. This is explained by the different energy dependencies of electron-neutral collision cross sections in atomic and molecular gases and, as a consequence, different drag force acting on electrons. A non-local behavior of electron energy distribution function is obtained in both gases owing to electrons are generated in the plasma sheath. In both gases, electron energy relaxation length is comparable withmore » the interelectrode gap, and therefore, they penetrate the plasma bulk with large energies.« less

Ozone kinetics in low-pressure discharges: vibrationally excited ozone and molecule formation on surfaces

NASA Astrophysics Data System (ADS)

Marinov, Daniil; Guerra, Vasco; Guaitella, Olivier; Booth, Jean-Paul; Rousseau, Antoine

2013-10-01

A combined experimental and modeling investigation of the ozone kinetics in the afterglow of pulsed direct current discharges in oxygen is carried out. The discharge is generated in a cylindrical silica tube of radius 1 cm, with short pulse durations between 0.5 and 2 ms, pressures in the range 1-5 Torr and discharge currents ˜40-120 mA. Time-resolved absolute concentrations of ground-state atoms and ozone molecules were measured simultaneously in situ, by two-photon absorption laser-induced fluorescence and ultraviolet absorption, respectively. The experiments were complemented by a self-consistent model developed to interpret the results and, in particular, to evaluate the roles of vibrationally excited ozone and of ozone formation on surfaces. It is found that vibrationally excited ozone, O_3^{*} , plays an important role in the ozone kinetics, leading to a decrease in the ozone concentration and an increase in its formation time. In turn, the kinetics of O_3^{*} is strongly coupled with those of atomic oxygen and O2(a 1Δg) metastables. Ozone formation at the wall does not contribute significantly to the total ozone production under the present conditions. Upper limits for the effective heterogeneous recombination probability of O atoms into ozone are established.

Modeling plasma heating by ns laser pulse

NASA Astrophysics Data System (ADS)

Colonna, Gianpiero; Laricchiuta, Annarita; Pietanza, Lucia Daniela

2018-03-01

The transition to breakdown of a weakly ionized gas, considering inverse bremsstrahlung, has been investigated using a state-to-state self-consistent model for gas discharges, mimicking a ns laser pulse. The paper is focused on the role of the initial ionization on the plasma formation. The results give the hint that some anomalous behaviors, such as signal enhancement by metal nanoparticles, can be attributed to this feature. This approach has been applied to hydrogen gas regarded as a simplified model for LIBS plasmas, as a full kinetic scheme is available, including the collisional-radiative model for atoms and molecules. The model allows the influence of different parameters to be investigated, such as the initial electron molar fraction, on the ionization growth.

Self-Consistent Superthermal Electron Effects on Plasmaspheric Refilling

NASA Technical Reports Server (NTRS)

Liemohn, M. W.; Khazanov, G. V.; Moore, T. E.; Guiter, S. M.

1997-01-01

The effects of self-consistently including superthermal electrons in the definition of the ambipolar electric field are investigated for the case of plasmaspheric refilling after a geomagnetic storm. By using the total electron population in the hydrodynamic equations, a method for incorporating superthermal electron parameters in the electric field and electron temperature calculation is developed. Also, the ambipolar electric field is included in the kinetic equation for the superthermal electrons through a change of variables using the total energy and the first adiabatic invariant. Calculations based on these changes are performed by coupling time-dependent models of the thermal plasma and superthermal electrons. Results from this treatment of the electric field and the self-consistent development of the solution are discussed in detail. Specifically, there is a decreased thermal electron density in the plasmasphere during the first few minutes of refilling, a slightly accelerated proton shock front, and a decreased superthermal electron flux due to the deceleration by the electric field. The timescales of plasmaspheric refilling are discussed and determined to be somewhat shorter than previously calculated for the thermal plasma and superthermal electron population due to the effects of the field-aligned potential.

Self-Consistent Magnetosphere-Ionosphere Coupling and Associated Plasma Energization Processes

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Six, N. Frank (Technical Monitor)

2002-01-01

Magnetosphere-Ionosphere (MI) coupling and associated with this process electron and ion energization processes have interested scientists for decades and, in spite of experimental and theoretical research efforts, are still ones of the least well known dynamic processes in space plasma physics. The reason for this is that the numerous physical processes associated with MI coupling occur over multiple spatial lengths and temporal scales. One typical example of MI coupling is large scale ring current (RC) electrodynamic coupling that includes calculation of the magnetospheric electric field that is consistent with the ring current (RC) distribution. A general scheme for numerical simulation of such large-scale magnetosphere-ionosphere coupling processes has been presented earlier in many works. The mathematical formulation of these models are based on "modified frozen-in flux theorem" for an ensemble of adiabatically drifting particles in the magnetosphere. By tracking the flow of particles through the inner magnetosphere, the bounce-averaged phase space density of the hot ions and electrons can be reconstructed and the magnetospheric electric field can be calculated such that it is consistent with the particle distribution in the magnetosphere. The new a self-consistent ring current model has been developed that couples electron and ion magnetospheric dynamics with calculation of electric field. Two new features were taken into account in addition to the RC ions, we solve an electron kinetic equation in our model, self-consistently including these results in the solution. Second, using different analytical relationships, we calculate the height integrated ionospheric conductances as the function of precipitated high energy magnetospheric electrons and ions as produced by our model. This results in fundamental changes to the electric potential pattern in the inner magnetosphere, with a smaller Alfven boundary than previous potential formulations would predict but one consistent with recent satellite observations. This leads to deeper penetration of the plasma sheet ions and electrons into the inner magnetosphere and more effective ring current ions and electron energization.

Ion cyclotron instability at Io: Hybrid simulation results compared to in situ observations

NASA Astrophysics Data System (ADS)

Šebek, Ondřej; Trávníček, Pavel M.; Walker, Raymond J.; Hellinger, Petr

2016-08-01

We present analysis of global three-dimensional hybrid simulations of Io's interaction with Jovian magnetospheric plasma. We apply a single-species model with simplified neutral-plasma chemistry and downscale Io in order to resolve the ion kinetic scales. We consider charge exchange, electron impact ionization, and photoionization by using variable rates of these processes to investigate their impact. Our results are in a good qualitative agreement with the in situ magnetic field measurements for five Galileo flybys around Io. The hybrid model describes ion kinetics self-consistently. This allows us to assess the distribution of temperature anisotropies around Io and thereby determine the possible triggering mechanism for waves observed near Io. We compare simulated dynamic spectra of magnetic fluctuations with in situ observations made by Galileo. Our results are consistent with both the spatial distribution and local amplitude of magnetic fluctuations found in the observations. Cyclotron waves, triggered probably by the growth of ion cyclotron instability, are observed mainly downstream of Io and on the flanks in regions farther from Io where the ion pickup rate is relatively low. Growth of the ion cyclotron instability is governed mainly by the charge exchange rate.

Collisional radiative model of an argon atmospheric capillary surface-wave discharge

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yanguas-Gil, A.; Cotrino, J.; Gonzalez-Elipe, A.R.

2004-12-01

The characteristics of a microwave surface-wave sustained plasma operated at atmospheric pressure in an open-ended dielectric tube are investigated theoretically as a first step in the development of a self-consistent model for these discharges. The plasma column is sustained in flowing argon. A surface-wave discharge that fills the whole radial cross section of the discharge tube is considered. With experimental electron temperature profiles [Garcia et al., Spectrochim. Acta, Part B 55, 1733 (2000)] the numerical model is used to test the validity of the different approximations and to study the influence of the different kinetic processes and power loss mechanismsmore » on the discharge.« less

Understanding the Elementary Steps in DNA Tile-Based Self-Assembly.

PubMed

Jiang, Shuoxing; Hong, Fan; Hu, Huiyu; Yan, Hao; Liu, Yan

2017-09-26

Although many models have been developed to guide the design and implementation of DNA tile-based self-assembly systems with increasing complexity, the fundamental assumptions of the models have not been thoroughly tested. To expand the quantitative understanding of DNA tile-based self-assembly and to test the fundamental assumptions of self-assembly models, we investigated DNA tile attachment to preformed "multi-tile" arrays in real time and obtained the thermodynamic and kinetic parameters of single tile attachment in various sticky end association scenarios. With more sticky ends, tile attachment becomes more thermostable with an approximately linear decrease in the free energy change (more negative). The total binding free energy of sticky ends is partially compromised by a sequence-independent energy penalty when tile attachment forms a constrained configuration: "loop". The minimal loop is a 2 × 2 tetramer (Loop4). The energy penalty of loops of 4, 6, and 8 tiles was analyzed with the independent loop model assuming no interloop tension, which is generalizable to arbitrary tile configurations. More sticky ends also contribute to a faster on-rate under isothermal conditions when nucleation is the rate-limiting step. Incorrect sticky end contributes to neither the thermostability nor the kinetics. The thermodynamic and kinetic parameters of DNA tile attachment elucidated here will contribute to the future improvement and optimization of tile assembly modeling, precise control of experimental conditions, and structural design for error-free self-assembly.

Dechlorination kinetics of TCE at toxic TCE concentrations: Assessment of different models.

PubMed

Haest, P J; Springael, D; Smolders, E

2010-01-01

The reductive dechlorination of trichloroethene (TCE) in a TCE source zone can be self-inhibited by TCE toxicity. A study was set up to examine the toxicity of TCE in terms of species specific degradation kinetics and microbial growth and to evaluate models that describe this self-inhibition. A batch experiment was performed using the TCE dechlorinating KB-1 culture at initial TCE concentrations ranging from 0.04mM to saturation (8.4mM). Biodegradation activity was highest at 0.3mM TCE and no activity was found at concentrations from 4 to 8mM. Species specific TCE and cis-DCE (cis-dichloroethene) degradation rates and Dehalococcoides numbers were modeled with Monod kinetics combined with either Haldane inhibition or a log-logistic dose-response inhibition on these rates. The log-logistic toxicity model appeared the most appropriate model and predicts that the species specific degradation activities are reduced by a factor 2 at about 1mM TCE, respectively cis-DCE. However, the model showed that the inhibitive effects on the time for TCE to ethene degradation are a complex function of degradation kinetics and the initial cell densities of the dechlorinating species. Our analysis suggests that the self-inhibition on biodegradation cannot be predicted by a single concentration threshold without information on the cell densities.

Simulation of Mercury's magnetosheath with a combined hybrid-paraboloid model

NASA Astrophysics Data System (ADS)

Parunakian, David; Dyadechkin, Sergey; Alexeev, Igor; Belenkaya, Elena; Khodachenko, Maxim; Kallio, Esa; Alho, Markku

2017-08-01

In this paper we introduce a novel approach for modeling planetary magnetospheres that involves a combination of the hybrid model and the paraboloid magnetosphere model (PMM); we further refer to it as the combined hybrid model. While both of these individual models have been successfully applied in the past, their combination enables us both to overcome the traditional difficulties of hybrid models to develop a self-consistent magnetic field and to compensate the lack of plasma simulation in the PMM. We then use this combined model to simulate Mercury's magnetosphere and investigate the geometry and configuration of Mercury's magnetosheath controlled by various conditions in the interplanetary medium. The developed approach provides a unique comprehensive view of Mercury's magnetospheric environment for the first time. Using this setup, we compare the locations of the bow shock and the magnetopause as determined by simulations with the locations predicted by stand-alone PMM runs and also verify the magnetic and dynamic pressure balance at the magnetopause. We also compare the results produced by these simulations with observational data obtained by the magnetometer on board the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft along a dusk-dawn orbit and discuss the signatures of the magnetospheric features that appear in these simulations. Overall, our analysis suggests that combining the semiempirical PMM with a self-consistent global kinetic model creates new modeling possibilities which individual models cannot provide on their own.

Modeling of ion orbit loss and intrinsic toroidal rotation with the COGENT code

NASA Astrophysics Data System (ADS)

Dorf, M.; Dorr, M.; Cohen, R.; Rognlien, T.; Hittinger, J.

2014-10-01

We discuss recent advances in cross-separatrix neoclassical transport simulations with COGENT, a continuum gyro-kinetic code being developed by the Edge Simulation Laboratory (ESL) collaboration. The COGENT code models the axisymmetric transport properties of edge plasmas including the effects of nonlinear (Fokker-Planck) collisions and a self-consistent electrostatic potential. Our recent work has focused on studies of ion orbit loss and the associated toroidal rotation driven by this mechanism. The results of the COGENT simulations are discussed and analyzed for the parameters of the DIII-D experiment. Work performed for USDOE at LLNL under Contract DE-AC52-07NA27344.

Kinetics on cocondensation of phenol and urea

Treesearch

Bunichiro Tomita; Yasunori Yoshida; Chung-Yun. Hse

1993-01-01

The chemical kinetics on cocondensation between methylolphenols and urea under acidic condition were investigated using 2- and 4-hydroxybenzyl alcohols as well as 2,4,6-trimethylolphenol as model compounds. The reactivity of the cocondensation were compared between o- and p-methylol groups. Moreover, the kinetics on self-condensations of monomethylolphenols and...

Numerical modeling of the solar wind flow with observational boundary conditions

DOE PAGES

Pogorelov, N. V.; Borovikov, S. N.; Burlaga, L. F.; ...

2012-11-20

In this paper we describe our group efforts to develop a self-consistent, data-driven model of the solar wind (SW) interaction with the local interstellar medium. The motion of plasma in this model is described with the MHD approach, while the transport of neutral atoms is addressed by either kinetic or multi-fluid equations. The model and its implementation in the Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) are continuously tested and validated by comparing our results with other models and spacecraft measurements. In particular, it was successfully applied to explain an unusual SW behavior discovered by the Voyager 1 spacecraft, i.e., the developmentmore » of a substantial negative radial velocity component, flow turning in the transverse direction, while the latitudinal velocity component goes to very small values. We explain recent SW velocity measurements at Voyager 1 in the context of our 3-D, MHD modeling. We also present a comparison of different turbulence models in their ability to reproduce the SW temperature profile from Voyager 2 measurements. Lastly, the boundary conditions obtained at 50 solar radii from data-driven numerical simulations are used to model a CME event throughout the heliosphere.« less

Self-consistent Model of Magnetospheric Electric Field, RC and EMIC Waves

NASA Technical Reports Server (NTRS)

Gamayunov, K. V.; Khazanov, G. V.; Liemohn, M. W.; Fok, M.-C.

2007-01-01

Electromagnetic ion cyclotron (EMIC) waves are an important magnetospheric emission, which is excited near the magnetic equator with frequencies below the proton gyro-frequency. The source of bee energy for wave growth is provided by temperature anisotropy of ring current (RC) ions, which develops naturally during inward convection from the plasma sheet These waves strongly affect the dynamic s of resonant RC ions, thermal electrons and ions, and the outer radiation belt relativistic electrons, leading to non-adiabatic particle heating and/or pitch-angle scattering and loss to the atmosphere. The rate of ion and electron scattering/heating is strongly controlled by the Wave power spectral and spatial distributions, but unfortunately, the currently available observational information regarding EMIC wave power spectral density is poor. So combinations of reliable data and theoretical models should be utilized in order to obtain the power spectral density of EMIC waves over the entire magnetosphere throughout the different storm phases. In this study, we present the simulation results, which are based on two coupled RC models that our group has developed. The first model deals with the large-scale magnetosphere-ionosphere electrodynamic coupling, and provides a self-consistent description of RC ions/electrons and the magnetospheric electric field. The second model is based on a coupled system of two kinetic equations, one equation describes the RC ion dynamics and another equation describes the power spectral density evolution of EMIC waves, and self-consistently treats a micro-scale electrodynamic coupling of RC and EMIC waves. So far, these two models have been applied independently. However, the large-scale magnetosphere-ionosphere electrodynamics controls the convective patterns of both the RC ions and plasmasphere altering conditions for EMIC wave-particle interaction. In turn, the wave induced RC precipitation Changes the local field-aligned current distributions and the ionospheric conductances, which are crucial for a large-scale electrodynamics. The initial results from this new self-consistent model of the magnetospheric electric field, RC and EMIC waves will be shown in this presentation.

Self-Heating Dark Matter via Semiannihilation

NASA Astrophysics Data System (ADS)

Kamada, Ayuki; Kim, Hee Jung; Kim, Hyungjin; Sekiguchi, Toyokazu

2018-03-01

The freeze-out of dark matter (DM) depends on the evolution of the DM temperature. The DM temperature does not have to follow the standard model one, when the elastic scattering is not sufficient to maintain the kinetic equilibrium. We study the temperature evolution of the semiannihilating DM, where a pair of the DM particles annihilate into one DM particle and another particle coupled to the standard model sector. We find that the kinetic equilibrium is maintained solely via semiannihilation until the last stage of the freeze-out. After the freeze-out, semiannihilation converts the mass deficit to the kinetic energy of DM, which leads to nontrivial evolution of the DM temperature. We argue that the DM temperature redshifts like radiation as long as the DM self-interaction is efficient. We dub this novel temperature evolution as self-heating. Notably, the structure formation is suppressed at subgalactic scales like keV-scale warm DM but with GeV-scale self-heating DM if the self-heating lasts roughly until the matter-radiation equality. The long duration of the self-heating requires the large self-scattering cross section, which in turn flattens the DM density profile in inner halos. Consequently, self-heating DM can be a unified solution to apparent failures of cold DM to reproduce the observed subgalactic scale structure of the Universe.

Toroidal Ampere-Faraday Equations Solved Simultaneously with CQL3D Fokker-Planck Time-Evolution

NASA Astrophysics Data System (ADS)

Harvey, R. W. (Bob); Petrov, Yu. V. (Yuri); Forest, C. B.; La Haye, R. J.

2017-10-01

A self-consistent, time-dependent toroidal electric field calculation is a key feature of a complete 3D Fokker-Planck kinetic distribution radial transport code for f(v,theta,rho,t). We discuss benchmarking and first applications of an implementation of the Ampere-Faraday equation for the self-consistent toroidal electric field, as applied to (1) resistive turn on of applied electron cyclotron current in the DIII-D tokamak giving initial back current adjacent to the direct CD region and having possible NTM stabilization implications, and (2) runaway electron production in tokamaks due to rapid reduction of the plasma temperature as occurs in pellet injection, massive gas injection, or a plasma disruption. Our previous results assuming a constant current density (Lenz' Law) model showed that prompt ``hot-tail runaways'' dominated ``knock-on'' and Dreicer ``drizzle'' runaways; we perform full-radius modeling and examine modifications due to the more complete Ampere-Faraday solution. Presently, the implementation relies on a fixed shape eqdsk, and this limitation will be addressed in future work. Research supported by USDOE FES award ER54744.

Current-Voltage and Floating-Potential characteristics of cylindrical emissive probes from a full-kinetic model based on the orbital motion theory

NASA Astrophysics Data System (ADS)

Chen, Xin; Sánchez-Arriaga, Gonzalo

2018-02-01

To model the sheath structure around an emissive probe with cylindrical geometry, the Orbital-Motion theory takes advantage of three conserved quantities (distribution function, transverse energy, and angular momentum) to transform the stationary Vlasov-Poisson system into a single integro-differential equation. For a stationary collisionless unmagnetized plasma, this equation describes self-consistently the probe characteristics. By solving such an equation numerically, parametric analyses for the current-voltage (IV) and floating-potential (FP) characteristics can be performed, which show that: (a) for strong emission, the space-charge effects increase with probe radius; (b) the probe can float at a positive potential relative to the plasma; (c) a smaller probe radius is preferred for the FP method to determine the plasma potential; (d) the work function of the emitting material and the plasma-ion properties do not influence the reliability of the floating-potential method. Analytical analysis demonstrates that the inflection point of an IV curve for non-emitting probes occurs at the plasma potential. The flat potential is not a self-consistent solution for emissive probes.

Improving the Ionospheric Auroral Conductance in a Global Ring Current Model and the Effects on the Ionospheric Electrodynamics

NASA Astrophysics Data System (ADS)

Yu, Y.; Jordanova, V. K.; McGranaghan, R. M.; Solomon, S. C.

2017-12-01

The ionospheric conductance, height-integrated electric conductivity, can regulate both the ionospheric electrodynamics and the magnetospheric dynamics because of its key role in determining the electric field within the coupled magnetosphere-ionosphere system. State-of-the-art global magnetosphere models commonly adopt empirical conductance calculators to obtain the auroral conductance. Such specification can bypass the complexity of the ionosphere-thermosphere chemistry but on the other hand breaks the self-consistent link within the coupled system. In this study, we couple a kinetic ring current model RAM-SCB-E that solves for anisotropic particle distributions with a two-stream electron transport code (GLOW) to more self-consistently compute the height-dependent electric conductivity, provided the auroral electron precipitation from the ring current model. Comparisons with the traditional empirical formula are carried out. It is found that the newly coupled modeling framework reveals smaller Hall and Pedersen conductance, resulting in a larger electric field. As a consequence, the subauroral polarization streams demonstrate a better agreement with observations from DMSP satellites. It is further found that the commonly assumed Maxwellian spectrum of the particle precipitation is not globally appropriate. Instead, a full precipitation spectrum resulted from wave particle interactions in the ring current accounts for a more comprehensive precipitation spectrum.

Theory verification and numerical benchmarking on neoclassical toroidal viscosity

NASA Astrophysics Data System (ADS)

Wang, Z. R.; Park, J.-K.; Liu, Y. Q.; Logan, N. C.; Menard, J. E.

2013-10-01

Systematic verification and numerical benchmarking has been successfully carried out among three different approaches of neoclassical toroidal viscosity (NTV) theory and the corresponding codes: IPEC-PENT is developed based on the combined NTV theory but without geometric simplifications; MARS-K originally calculating the kinetic energy is upgraded to calculate the NTV torque based on the equivalence between kinetic energy and NTV torque; MARS-Q includes smoothly connected NTV formula. The derivation and numerical results both indicate that the imaginary part of kinetic energy calculated by MARS-K is equivalent to the NTV torque in IPEC-PENT. In the benchmark of precession resonance between MARS-Q and MARS-K/IPEC-PENT, it is first time to show the agreement and the correlation between the connected NTV formula and the combined NTV theory in different collisional region. Additionally, both IPEC-PENT and MARS-K indicates the importance of the bounce harmonic resonance which could greatly enhance the NTV torque when E cross B drift frequency reaches the bounce resonance condition. Since MARS-K also has the capability to calculate the plasma response including the kinetic effect self-consistently, the self-consistent NTV torque calculations have also been tested. This work is supported by DOE Contract No. DE-AC02-09CH11466.

Demonstration of Anisotropic Fluid Closure Capturing the Kinetic Structure of Magnetic Reconnection

NASA Astrophysics Data System (ADS)

Ohia, Obioma

2012-10-01

Magnetic reconnection in collisionless plasmas plays an important role in space and laboratory plasmas. Allowing magnetic stress to be reduced by a rearrangement of magnetic line topology, this process is often accompanied by a large release of magnetic field energy, which can heat the plasma, drive large scale flows, or accelerate particles. Reconnection has been widely studied through fluid models and kinetic simulations. While two-fluid models often reproduce the fast reconnection that is observed in nature and seen in kinetic simulations, it is found that the structure surrounding the electron diffusion region and the electron current layer differ vastly between fluid models and kinetic simulations [1]. Recently, using an adiabatic solution of the Vlasov equation, a new fluid closure has been obtained for electrons that relate parallel and perpendicular pressures to the density and magnetic field [2]. Here we present the results of fluid simulation, developed using the HiFi framework [3], that implements new equations of state for guide-field reconnection. The new fluid closure accurately accounts for the anisotropic electron pressure that builds in the reconnection region due to electric and magnetic trapping of electrons. In contrast to previous fluid models, our fluid simulation reproduces the detailed reconnection region as observed in fully kinetic simulations [4]. We hereby demonstrate that the new fluid closure self-consistently captures all the physics relevant to the structure of the reconnection region, providing a gateway to a renewed and deeper theoretical understanding for reconnection in weakly collisional regimes.[4pt] [1] Daughton W et al., Phys. Plasmas 13, 072101 (2006).[0pt] [2] Le A et al., Phys. Rev. Lett. 102, 085001 (2009). [0pt] [3] Lukin VS, Linton MG, Nonlinear Proc. Geoph. 18, 871 (2011). [0pt] [4] Ohia O, et al., Phys. Rev. Lett. In Press (2012).

Ozone kinetics in low-pressure discharges

NASA Astrophysics Data System (ADS)

Guerra, Vasco; Marinov, Daniil; Guaitella, Olivier; Rousseau, Antoine

2012-10-01

Ozone kinetics is quite well established at atmospheric pressure, due to the importance of ozone in atmospheric chemistry and to the development of industrial ozone reactors. However, as the pressure is decreased and the dominant three-body reactions lose importance, the main mechanisms involved in the creation and destruction of ozone are still surrounded by important uncertainties. In this work we develop a self-consistent model for a pulsed discharge and its afterglow operating in a Pyrex reactor with inner radius 1 cm, at pressures in the range 1-5 Torr and discharge currents of 40-120 mA. The model couples the electron Boltzmann equation with a system of equations for the time evolution of the heavy particles. The calculations are compared with time-dependent measurements of ozone and atomic oxygen. Parametric studies are performed in order to clarify the role of vibrationally excited ozone in the overall kinetics and to establish the conditions where ozone production on the surface may become important. It is shown that vibrationally excited ozone does play a significant role, by increasing the time constants of ozone formation. Moreover, an upper limit for the ozone formation at the wall in these conditions is set at 10(-4).

Parameter estimation in large-scale systems biology models: a parallel and self-adaptive cooperative strategy.

PubMed

Penas, David R; González, Patricia; Egea, Jose A; Doallo, Ramón; Banga, Julio R

2017-01-21

The development of large-scale kinetic models is one of the current key issues in computational systems biology and bioinformatics. Here we consider the problem of parameter estimation in nonlinear dynamic models. Global optimization methods can be used to solve this type of problems but the associated computational cost is very large. Moreover, many of these methods need the tuning of a number of adjustable search parameters, requiring a number of initial exploratory runs and therefore further increasing the computation times. Here we present a novel parallel method, self-adaptive cooperative enhanced scatter search (saCeSS), to accelerate the solution of this class of problems. The method is based on the scatter search optimization metaheuristic and incorporates several key new mechanisms: (i) asynchronous cooperation between parallel processes, (ii) coarse and fine-grained parallelism, and (iii) self-tuning strategies. The performance and robustness of saCeSS is illustrated by solving a set of challenging parameter estimation problems, including medium and large-scale kinetic models of the bacterium E. coli, bakerés yeast S. cerevisiae, the vinegar fly D. melanogaster, Chinese Hamster Ovary cells, and a generic signal transduction network. The results consistently show that saCeSS is a robust and efficient method, allowing very significant reduction of computation times with respect to several previous state of the art methods (from days to minutes, in several cases) even when only a small number of processors is used. The new parallel cooperative method presented here allows the solution of medium and large scale parameter estimation problems in reasonable computation times and with small hardware requirements. Further, the method includes self-tuning mechanisms which facilitate its use by non-experts. We believe that this new method can play a key role in the development of large-scale and even whole-cell dynamic models.

Development and application of kinetic model on biological anoxic/aerobic filter.

PubMed

Kim, Youngnoh; Tanaka, Kazuhiro; Lee, Yong-Woo; Chung, Jinwook

2008-01-01

An up-flow biological anoxic filter (BANF) has been developed to achieve high removal performance of suspended solids and BOD removal as well as nitrogen. With a view to understand treatment mechanisms, we developed a filtration model that incorporates filtration, deposit scoring and biological reactions simultaneously. The biological reactions consist of four types of reaction; dissolution of organic particles; utilization of dissolved organic matter; denitrification; and self-degradation of bacteria. Whereas the reactor is generally assumed to be a plug flow reactor in the filtration model, it is assumed a continuous-flow stirred tank reactor (CSTR) in the model of biological reactions. The hydrodynamics is supposed that the filter bottom (the portion sludge settled) is a CSTR and the filter bed (the portion filled with filter media) consists of number of CSTR of equal size arranged in series. The model obtained in this study was verified and simulated using experimental results taken from a pilot-scale plant and predicted the experimental data well, applying to design and operate BANF.

Coupling Fluid and Kineitc Effects in Space Weather: an interdisciplinary task

NASA Astrophysics Data System (ADS)

Lapenta, Giovanni; González-Herrero, Diego; Boella, Elisabetta; Siddi, Lorenzo; Cazzola, Emanuele

2017-04-01

Two agents are key to space weather: electromagentic fields and energetic particles. Magnetic fields carried by plasmas in the solar wind interact with the Earth magnetosphere and solar energetic particles produced by solar events or in cosmic rays affect the space environment. Describing both is challenging. Magnetized plasmas are most effectively described by magneto-hydrodynamics, MHD, a fluid theory based on describing some fields defined in space: electromagnetic fields, density, velocity and temperature of the plasma. High energy particles, instead need a more detailed approach , kinetic theory, where statistical distributions of particles are governed by the Boltzmann equation. While fluid models are based on the ordinary space and time, kinetic models require a six dimensional space, called phase space, besides time. The two methods are not separated, the processes leading to the production of energetic particles are the same that involve space plasamas and fields. Arriving at a single self-consistent model has been the goal of the Swiff project funded by the EC in FP7 and it is now a key goal of the ongoing DEEP-ER project. We present a new approach developed with the goal of extending the reach of kinetic models to the fluid scales. Kinetic models are a higher order description and all fluid effects are included in them. However, the cost in terms of computing power is much higher and it has been so far prohibitively expensive to treat space weather events fully kinetically. We have now designed a new method capable of reducing that cost by several orders of magnitude making it possible for kinetic models to study space weather events [1,2]. We will report the new methodology and show its application to space weather mdeling. [1] Giovanni Lapenta,Exactly Energy Conserving Semi-Implicit Particle in Cell Formulation, to appear, JCP, arXiv:1602.06326 [2] Giovanni Lapenta, Diego Gonzalez-Herrero, Elisabetta Boella, Multiple scale kinetic simulations with the energy conserving semi implicit particle in cell (PIC) method, submitted JPP, arXiv:1612.08289

Coarse-grained molecular dynamics modeling of the kinetics of lamellar block copolymer defect annealing

NASA Astrophysics Data System (ADS)

Peters, Andrew J.; Lawson, Richard A.; Nation, Benjamin D.; Ludovice, Peter J.; Henderson, Clifford L.

2016-01-01

State-of-the-art block copolymer (BCP)-directed self-assembly (DSA) methods still yield defect densities orders of magnitude higher than is necessary in semiconductor fabrication despite free-energy calculations that suggest equilibrium defect densities are much lower than is necessary for economic fabrication. This disparity suggests that the main problem may lie in the kinetics of defect removal. This work uses a coarse-grained model to study the rates, pathways, and dependencies of healing a common defect to give insight into the fundamental processes that control defect healing and give guidance on optimal process conditions for BCP-DSA. It is found that bulk simulations yield an exponential drop in defect heal rate above χN˜30. Thin films show no change in rate associated with the energy barrier below χN˜50, significantly higher than the χN values found previously for self-consistent field theory studies that neglect fluctuations. Above χN˜50, the simulations show an increase in energy barrier scaling with 1/2 to 1/3 of the bulk systems. This is because thin films always begin healing at the free interface or the BCP-underlayer interface, where the increased A-B contact area associated with the transition state is minimized, while the infinitely thick films cannot begin healing at an interface.

A degradation function consistent with Cocks–Ashby porosity kinetics

DOE PAGES

Moore, John A.

2017-10-14

Here, the load carrying capacity of ductile materials degrades as a function of porosity, stress state and strain-rate. The effect of these variables on porosity kinetics is captured by the Cocks–Ashby model; however, the Cocks–Ashby model does not account for material degradation directly. This work uses a yield criteria to form a degradation function that is consistent with Cocks–Ashby porosity kinetics and is a function of porosity, stress state and strain-rate dependence. Approximations of this degradation function for pure hydrostatic stress states are also explored.

A degradation function consistent with Cocks–Ashby porosity kinetics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Moore, John A.

Here, the load carrying capacity of ductile materials degrades as a function of porosity, stress state and strain-rate. The effect of these variables on porosity kinetics is captured by the Cocks–Ashby model; however, the Cocks–Ashby model does not account for material degradation directly. This work uses a yield criteria to form a degradation function that is consistent with Cocks–Ashby porosity kinetics and is a function of porosity, stress state and strain-rate dependence. Approximations of this degradation function for pure hydrostatic stress states are also explored.

Kinetic theory approach to modeling of cellular repair mechanisms under genome stress.

PubMed

Qi, Jinpeng; Ding, Yongsheng; Zhu, Ying; Wu, Yizhi

2011-01-01

Under acute perturbations from outer environment, a normal cell can trigger cellular self-defense mechanism in response to genome stress. To investigate the kinetics of cellular self-repair process at single cell level further, a model of DNA damage generating and repair is proposed under acute Ion Radiation (IR) by using mathematical framework of kinetic theory of active particles (KTAP). Firstly, we focus on illustrating the profile of Cellular Repair System (CRS) instituted by two sub-populations, each of which is made up of the active particles with different discrete states. Then, we implement the mathematical framework of cellular self-repair mechanism, and illustrate the dynamic processes of Double Strand Breaks (DSBs) and Repair Protein (RP) generating, DSB-protein complexes (DSBCs) synthesizing, and toxins accumulating. Finally, we roughly analyze the capability of cellular self-repair mechanism, cellular activity of transferring DNA damage, and genome stability, especially the different fates of a certain cell before and after the time thresholds of IR perturbations that a cell can tolerate maximally under different IR perturbation circumstances.

Integrated modeling applications for tokamak experiments with OMFIT

NASA Astrophysics Data System (ADS)

Meneghini, O.; Smith, S. P.; Lao, L. L.; Izacard, O.; Ren, Q.; Park, J. M.; Candy, J.; Wang, Z.; Luna, C. J.; Izzo, V. A.; Grierson, B. A.; Snyder, P. B.; Holland, C.; Penna, J.; Lu, G.; Raum, P.; McCubbin, A.; Orlov, D. M.; Belli, E. A.; Ferraro, N. M.; Prater, R.; Osborne, T. H.; Turnbull, A. D.; Staebler, G. M.

2015-08-01

One modeling framework for integrated tasks (OMFIT) is a comprehensive integrated modeling framework which has been developed to enable physics codes to interact in complicated workflows, and support scientists at all stages of the modeling cycle. The OMFIT development follows a unique bottom-up approach, where the framework design and capabilities organically evolve to support progressive integration of the components that are required to accomplish physics goals of increasing complexity. OMFIT provides a workflow for easily generating full kinetic equilibrium reconstructions that are constrained by magnetic and motional Stark effect measurements, and kinetic profile information that includes fast-ion pressure modeled by a transport code. It was found that magnetic measurements can be used to quantify the amount of anomalous fast-ion diffusion that is present in DIII-D discharges, and provide an estimate that is consistent with what would be needed for transport simulations to match the measured neutron rates. OMFIT was used to streamline edge-stability analyses, and evaluate the effect of resonant magnetic perturbation (RMP) on the pedestal stability, which have been found to be consistent with the experimental observations. The development of a five-dimensional numerical fluid model for estimating the effects of the interaction between magnetohydrodynamic (MHD) and microturbulence, and its systematic verification against analytic models was also supported by the framework. OMFIT was used for optimizing an innovative high-harmonic fast wave system proposed for DIII-D. For a parallel refractive index {{n}\\parallel}>3 , the conditions for strong electron-Landau damping were found to be independent of launched {{n}\\parallel} and poloidal angle. OMFIT has been the platform of choice for developing a neural-network based approach to efficiently perform a non-linear multivariate regression of local transport fluxes as a function of local dimensionless parameters. Transport predictions for thousands of DIII-D discharges showed excellent agreement with the power balance calculations across the whole plasma radius and over a broad range of operating regimes. Concerning predictive transport simulations, the framework made possible the design and automation of a workflow that enables self-consistent predictions of kinetic profiles and the plasma equilibrium. It is found that the feedback between the transport fluxes and plasma equilibrium can significantly affect the kinetic profiles predictions. Such a rich set of results provide tangible evidence of how bottom-up approaches can potentially provide a fast track to integrated modeling solutions that are functional, cost-effective, and in sync with the research effort of the community.

LIGKA: A linear gyrokinetic code for the description of background kinetic and fast particle effects on the MHD stability in tokamaks

NASA Astrophysics Data System (ADS)

Lauber, Ph.; Günter, S.; Könies, A.; Pinches, S. D.

2007-09-01

In a plasma with a population of super-thermal particles generated by heating or fusion processes, kinetic effects can lead to the additional destabilisation of MHD modes or even to additional energetic particle modes. In order to describe these modes, a new linear gyrokinetic MHD code has been developed and tested, LIGKA (linear gyrokinetic shear Alfvén physics) [Ph. Lauber, Linear gyrokinetic description of fast particle effects on the MHD stability in tokamaks, Ph.D. Thesis, TU München, 2003; Ph. Lauber, S. Günter, S.D. Pinches, Phys. Plasmas 12 (2005) 122501], based on a gyrokinetic model [H. Qin, Gyrokinetic theory and computational methods for electromagnetic perturbations in tokamaks, Ph.D. Thesis, Princeton University, 1998]. A finite Larmor radius expansion together with the construction of some fluid moments and specification to the shear Alfvén regime results in a self-consistent, electromagnetic, non-perturbative model, that allows not only for growing or damped eigenvalues but also for a change in mode-structure of the magnetic perturbation due to the energetic particles and background kinetic effects. Compared to previous implementations [H. Qin, mentioned above], this model is coded in a more general and comprehensive way. LIGKA uses a Fourier decomposition in the poloidal coordinate and a finite element discretisation in the radial direction. Both analytical and numerical equilibria can be treated. Integration over the unperturbed particle orbits is performed with the drift-kinetic HAGIS code [S.D. Pinches, Ph.D. Thesis, The University of Nottingham, 1996; S.D. Pinches et al., CPC 111 (1998) 131] which accurately describes the particles' trajectories. This allows finite-banana-width effects to be implemented in a rigorous way since the linear formulation of the model allows the exchange of the unperturbed orbit integration and the discretisation of the perturbed potentials in the radial direction. Successful benchmarks for toroidal Alfvén eigenmodes (TAEs) and kinetic Alfvén waves (KAWs) with analytical results, ideal MHD codes, drift-kinetic codes and other codes based on kinetic models are reported.

Modeling the heating and atomic kinetics of a photoionized neon plasma experiment

NASA Astrophysics Data System (ADS)

Lockard, Tom E.

Motivated by gas cell photoionized plasma experiments performed by our group at the Z facility of Sandia National Laboratories, we discuss in this dissertation a modeling study of the heating and ionization of the plasma for conditions characteristic of these experiments. Photoionized plasmas are non-equilibrium systems driven by a broadband x-ray radiation flux. They are commonly found in astrophysics but rarely seen in the laboratory. Several modeling tools have been employed: (1) a view-factor computer code constrained with side x-ray power and gated monochromatic image measurements of the z-pinch radiation, to model the time-history of the photon-energy resolved x-ray flux driving the photoionized plasma, (2) a Boltzmann self-consistent electron and atomic kinetics model to simulate the electron distribution function and configuration-averaged atomic kinetics, (3) a radiation-hydrodynamics code with inline non-equilibrium atomic kinetics to perform a comprehensive numerical simulation of the experiment and plasma heating, and (4) steady-state and time-dependent collisional-radiative atomic kinetics calculations with fine-structure energy level description to assess transient effects in the ionization and charge state distribution of the plasma. The results indicate that the photon-energy resolved x-ray flux impinging on the front window of the gas cell is very well approximated by a linear combination of three geometrically-diluted Planckian distributions. Knowledge of the spectral details of the x-ray drive turned out to be important for the heating and ionization of the plasma. The free electrons in the plasma thermalize quickly relative to the timescales associated with the time-history of the x-ray drive and the plasma atomic kinetics. Hence, electrons are well described by a Maxwellian energy distribution of a single temperature. This finding is important to support the application of a radiation-hydrodynamic model to simulate the experiment. It is found that the computed plasma heating compares well with experimental observation when the effects of the windows, hydrodynamics, and non-equilbirium neon emissivity and opacity are employed. The atomic kinetics shows significant time-dependent effects because the timescale of the x-ray drive is too short compared to that of the photoionization process. These modeling and simulation results are important to test theory and modeling assumptions and approximations, and also to provide guidance on data interpretation and analysis.

Analysis and Modeling of DIII-D Experiments With OMFIT and Neural Networks

NASA Astrophysics Data System (ADS)

Meneghini, O.; Luna, C.; Smith, S. P.; Lao, L. L.; GA Theory Team

2013-10-01

The OMFIT integrated modeling framework is designed to facilitate experimental data analysis and enable integrated simulations. This talk introduces this framework and presents a selection of its applications to the DIII-D experiment. Examples include kinetic equilibrium reconstruction analysis; evaluation of MHD stability in the core and in the edge; and self-consistent predictive steady-state transport modeling. The OMFIT framework also provides the platform for an innovative approach based on neural networks to predict electron and ion energy fluxes. In our study a multi-layer feed-forward back-propagation neural network is built and trained over a database of DIII-D data. It is found that given the same parameters that the highest fidelity models use, the neural network model is able to predict to a large degree the heat transport profiles observed in the DIII-D experiments. Once the network is built, the numerical cost of evaluating the transport coefficients is virtually nonexistent, thus making the neural network model particularly well suited for plasma control and quick exploration of operational scenarios. The implementation of the neural network model and benchmark with experimental results and gyro-kinetic models will be discussed. Work supported in part by the US DOE under DE-FG02-95ER54309.

Kinetics of motility-induced phase separation and swim pressure

NASA Astrophysics Data System (ADS)

Patch, Adam; Yllanes, David; Marchetti, M. Cristina

2017-01-01

Active Brownian particles (ABPs) represent a minimal model of active matter consisting of self-propelled spheres with purely repulsive interactions and rotational noise. Here we examine the pressure of ABPs in two dimensions in both closed boxes and systems with periodic boundary conditions and show that its nonmonotonic behavior with density is a general property of ABPs and is not the result of finite-size effects. We correlate the time evolution of the mean pressure towards its steady-state value with the kinetics of motility-induced phase separation. For parameter values corresponding to phase-separated steady states, we identify two dynamical regimes. The pressure grows monotonically in time during the initial regime of rapid cluster formation, overshooting its steady-state value and then quickly relaxing to it, and remains constant during the subsequent slower period of cluster coalescence and coarsening. The overshoot is a distinctive feature of active systems.

From Excessive Journal Self-Cites to Citation Stacking: Analysis of Journal Self-Citation Kinetics in Search for Journals, Which Boost Their Scientometric Indicators.

PubMed

Heneberg, Petr

2016-01-01

Bibliometric indicators increasingly affect careers, funding, and reputation of individuals, their institutions and journals themselves. In contrast to author self-citations, little is known about kinetics of journal self-citations. Here we hypothesized that they may show a generalizable pattern within particular research fields or across multiple fields. We thus analyzed self-cites to 60 journals from three research fields (multidisciplinary sciences, parasitology, and information science). We also hypothesized that the kinetics of journal self-citations and citations received from other journals of the same publisher may differ from foreign citations. We analyzed the journals published the American Association for the Advancement of Science, Nature Publishing Group, and Editura Academiei Române. We found that although the kinetics of journal self-cites is generally faster compared to foreign cites, it shows some field-specific characteristics. Particularly in information science journals, the initial increase in a share of journal self-citations during post-publication year 0 was completely absent. Self-promoting journal self-citations of top-tier journals have rather indirect but negligible direct effects on bibliometric indicators, affecting just the immediacy index and marginally increasing the impact factor itself as long as the affected journals are well established in their fields. In contrast, other forms of journal self-citations and citation stacking may severely affect the impact factor, or other citation-based indices. We identified here a network consisting of three Romanian physics journals Proceedings of the Romanian Academy, Series A, Romanian Journal of Physics, and Romanian Reports in Physics, which displayed low to moderate ratio of journal self-citations, but which multiplied recently their impact factors, and were mutually responsible for 55.9%, 64.7% and 63.3% of citations within the impact factor calculation window to the three journals, respectively. They did not receive nearly any network self-cites prior impact factor calculation window, and their network self-cites decreased sharply after the impact factor calculation window. Journal self-citations and citation stacking requires increased attention and elimination from citation indices.

From Excessive Journal Self-Cites to Citation Stacking: Analysis of Journal Self-Citation Kinetics in Search for Journals, Which Boost Their Scientometric Indicators

PubMed Central

2016-01-01

Bibliometric indicators increasingly affect careers, funding, and reputation of individuals, their institutions and journals themselves. In contrast to author self-citations, little is known about kinetics of journal self-citations. Here we hypothesized that they may show a generalizable pattern within particular research fields or across multiple fields. We thus analyzed self-cites to 60 journals from three research fields (multidisciplinary sciences, parasitology, and information science). We also hypothesized that the kinetics of journal self-citations and citations received from other journals of the same publisher may differ from foreign citations. We analyzed the journals published the American Association for the Advancement of Science, Nature Publishing Group, and Editura Academiei Române. We found that although the kinetics of journal self-cites is generally faster compared to foreign cites, it shows some field-specific characteristics. Particularly in information science journals, the initial increase in a share of journal self-citations during post-publication year 0 was completely absent. Self-promoting journal self-citations of top-tier journals have rather indirect but negligible direct effects on bibliometric indicators, affecting just the immediacy index and marginally increasing the impact factor itself as long as the affected journals are well established in their fields. In contrast, other forms of journal self-citations and citation stacking may severely affect the impact factor, or other citation-based indices. We identified here a network consisting of three Romanian physics journals Proceedings of the Romanian Academy, Series A, Romanian Journal of Physics, and Romanian Reports in Physics, which displayed low to moderate ratio of journal self-citations, but which multiplied recently their impact factors, and were mutually responsible for 55.9%, 64.7% and 63.3% of citations within the impact factor calculation window to the three journals, respectively. They did not receive nearly any network self-cites prior impact factor calculation window, and their network self-cites decreased sharply after the impact factor calculation window. Journal self-citations and citation stacking requires increased attention and elimination from citation indices. PMID:27088862

Signatures of a macroscopic switching transition for a dynamic microtubule

NASA Astrophysics Data System (ADS)

Aparna, J. S.; Padinhateeri, Ranjith; Das, Dibyendu

2017-04-01

Characterising complex kinetics of non-equilibrium self-assembly of bio-filaments is of general interest. Dynamic instability in microtubules, consisting of successive catastrophes and rescues, is observed to occur as a result of the non-equilibrium conversion of GTP-tubulin to GDP-tubulin. We study this phenomenon using a model for microtubule kinetics with GTP/GDP state-dependent polymerisation, depolymerisation and hydrolysis of subunits. Our results reveal a sharp switch-like transition in the mean velocity of the filaments, from a growth phase to a shrinkage phase, with an associated co-existence of the two phases. This transition is reminiscent of the discontinuous phase transition across the liquid-gas boundary. We probe the extent of discontinuity in the transition quantitatively using characteristic signatures such as bimodality in velocity distribution, variance and Binder cumulant, and also hysteresis behaviour of the system. We further investigate ageing behaviour in catastrophes of the filament, and find that the multi-step nature of catastrophes is intensified in the vicinity of the switching transition. This assumes importance in the context of Microtubule Associated Proteins which have the potential of altering kinetic parameter values.

Effects of Different PER Translational Kinetics on the Dynamics of a Core Circadian Clock Model

PubMed Central

Nieto, Paula S.; Revelli, Jorge A.; Garbarino-Pico, Eduardo; Condat, Carlos A.; Guido, Mario E.; Tamarit, Francisco A.

2015-01-01

Living beings display self-sustained daily rhythms in multiple biological processes, which persist in the absence of external cues since they are generated by endogenous circadian clocks. The period (per) gene is a central player within the core molecular mechanism for keeping circadian time in most animals. Recently, the modulation PER translation has been reported, both in mammals and flies, suggesting that translational regulation of clock components is important for the proper clock gene expression and molecular clock performance. Because translational regulation ultimately implies changes in the kinetics of translation and, therefore, in the circadian clock dynamics, we sought to study how and to what extent the molecular clock dynamics is affected by the kinetics of PER translation. With this objective, we used a minimal mathematical model of the molecular circadian clock to qualitatively characterize the dynamical changes derived from kinetically different PER translational mechanisms. We found that the emergence of self-sustained oscillations with characteristic period, amplitude, and phase lag (time delays) between per mRNA and protein expression depends on the kinetic parameters related to PER translation. Interestingly, under certain conditions, a PER translation mechanism with saturable kinetics introduces longer time delays than a mechanism ruled by a first-order kinetics. In addition, the kinetic laws of PER translation significantly changed the sensitivity of our model to parameters related to the synthesis and degradation of per mRNA and PER degradation. Lastly, we found a set of parameters, with realistic values, for which our model reproduces some experimental results reported recently for Drosophila melanogaster and we present some predictions derived from our analysis. PMID:25607544

Effects of different per translational kinetics on the dynamics of a core circadian clock model.

PubMed

Nieto, Paula S; Revelli, Jorge A; Garbarino-Pico, Eduardo; Condat, Carlos A; Guido, Mario E; Tamarit, Francisco A

2015-01-01

Living beings display self-sustained daily rhythms in multiple biological processes, which persist in the absence of external cues since they are generated by endogenous circadian clocks. The period (per) gene is a central player within the core molecular mechanism for keeping circadian time in most animals. Recently, the modulation PER translation has been reported, both in mammals and flies, suggesting that translational regulation of clock components is important for the proper clock gene expression and molecular clock performance. Because translational regulation ultimately implies changes in the kinetics of translation and, therefore, in the circadian clock dynamics, we sought to study how and to what extent the molecular clock dynamics is affected by the kinetics of PER translation. With this objective, we used a minimal mathematical model of the molecular circadian clock to qualitatively characterize the dynamical changes derived from kinetically different PER translational mechanisms. We found that the emergence of self-sustained oscillations with characteristic period, amplitude, and phase lag (time delays) between per mRNA and protein expression depends on the kinetic parameters related to PER translation. Interestingly, under certain conditions, a PER translation mechanism with saturable kinetics introduces longer time delays than a mechanism ruled by a first-order kinetics. In addition, the kinetic laws of PER translation significantly changed the sensitivity of our model to parameters related to the synthesis and degradation of per mRNA and PER degradation. Lastly, we found a set of parameters, with realistic values, for which our model reproduces some experimental results reported recently for Drosophila melanogaster and we present some predictions derived from our analysis.

Growth of vertically aligned nanowires in metal-oxide nanocomposites: kinetic Monte-Carlo modeling versus experiments.

PubMed

Hennes, M; Schuler, V; Weng, X; Buchwald, J; Demaille, D; Zheng, Y; Vidal, F

2018-04-26

We employ kinetic Monte-Carlo simulations to study the growth process of metal-oxide nanocomposites obtained via sequential pulsed laser deposition. Using Ni-SrTiO3 (Ni-STO) as a model system, we reduce the complexity of the computational problem by choosing a coarse-grained approach mapping Sr, Ti and O atoms onto a single effective STO pseudo-atom species. With this ansatz, we scrutinize the kinetics of the sequential synthesis process, governed by alternating deposition and relaxation steps, and analyze the self-organization propensity of Ni atoms into straight vertically aligned nanowires embedded in the surrounding STO matrix. We finally compare the predictions of our binary toy model with experiments and demonstrate that our computational approach captures fundamental aspects of self-assembled nanowire synthesis. Despite its simplicity, our modeling strategy successfully describes the impact of relevant parameters like the concentration or laser frequency on the final nanoarchitecture of metal-oxide thin films grown via pulsed laser deposition.

On the Voyager 1 Zero Radial Velocity Measurements in the Inner Heliosheath

NASA Astrophysics Data System (ADS)

Pogorelov, N. V.; Borovikov, S. N.; Kryukov, I.; Zank, G. P.

2011-12-01

Theoretical analysis of the Voyager 1 data revealed a very small, or even negative, value of the solar wind (SW) radial velocity component. This should not be surprising if we take into account time-dependent processes that take place in the inner heliosheath (IHS). We analyze solar cycle modeling of the SW interaction with the local interstellar medium (LISM) and demonstrate the existence of small and negative values of the SW radial velocity. It shown that, in reality, a similar picture can be observed in the outer heliosheath (OHS), where on the contrary, extended regions of the positive radial velocity are observed. Another scenario discussed in this talk is related to effects of transients, such as global merged interaction regions and corotating interacting regions. Numerical results are obtained with the SW-LISM interaction model developed in the UAHuntsville and implemented in the Multi-Scale Fluid-Kinetic Simulation Suite. This model treats ions magnetohydrodynamically while the transport of neutral atoms is performed kinetically by solving the Boltzmann equation with a Monte Carlo approach or using a multi-fluid approach developed in mid-90's by G. P. Zank. Pickup ions (PUIs) can be treated as a separate fluid. The evolution of the PUI-generated turbulence is addressed on the differential level by adding three additional equations, which are solved self-consistently with the MHD-kinetic system of equations.

Three-dimensional Kinetic Pulsar Magnetosphere Models: Connecting to Gamma-Ray Observations

NASA Astrophysics Data System (ADS)

Kalapotharakos, Constantinos; Brambilla, Gabriele; Timokhin, Andrey; Harding, Alice K.; Kazanas, Demosthenes

2018-04-01

We present three-dimensional (3D) global kinetic pulsar magnetosphere models, where the charged particle trajectories and the corresponding electromagnetic fields are treated self-consistently. For our study, we have developed a Cartesian 3D relativistic particle-in-cell code that incorporates radiation reaction forces. We describe our code and discuss the related technical issues, treatments, and assumptions. Injecting particles up to large distances in the magnetosphere, we apply arbitrarily low to high particle injection rates, and obtain an entire spectrum of solutions from close to the vacuum-retarded dipole to close to the force-free (FF) solution, respectively. For high particle injection rates (close to FF solutions), significant accelerating electric field components are confined only near the equatorial current sheet outside the light cylinder. A judicious interpretation of our models allows the particle emission to be calculated, and consequently, the corresponding realistic high-energy sky maps and spectra to be derived. Using model parameters that cover the entire range of spin-down powers of Fermi young and millisecond pulsars, we compare the corresponding model γ-ray light curves, cutoff energies, and total γ-ray luminosities with those observed by Fermi to discover a dependence of the particle injection rate, { \\mathcal F }, on the spin-down power, \\dot{{ \\mathcal E }}, indicating an increase of { \\mathcal F } with \\dot{{ \\mathcal E }}. Our models, guided by Fermi observations, provide field structures and particle distributions that are not only consistent with each other but also able to reproduce a broad range of the observed γ-ray phenomenologies of both young and millisecond pulsars.

ATOMIC-SCALE DESIGN OF IRON FISCHER-TROPSCH CATALYSTS: A COMBINED COMPUTATIONAL CHEMISTRY, EXPERIMENTAL, AND MICROKINETIC MODELING APPROACH

DOE Office of Scientific and Technical Information (OSTI.GOV)

Manos Mavrikakis; James A. Dumesic; Amit A. Gokhale

2005-03-22

Efforts during this first year focused on four areas: (1) searching/summarizing published FTS mechanistic and kinetic studies of FTS reactions on iron catalysts; (2) construction of mass spectrometer-TPD and Berty CSTR reactor systems; (3) preparation and characterization of unsupported iron and alumina-supported iron catalysts at various iron loadings (4) Determination of thermochemical parameters such as binding energies of reactive intermediates, heat of FTS elementary reaction steps, and kinetic parameters such as activation energies, and frequency factors of FTS elementary reaction steps on a number of model surfaces. Literature describing mechanistic and kinetic studies of Fischer-Tropsch synthesis on iron catalysts wasmore » compiled in a draft review. Construction of the mass spectrometer-TPD system is 90% complete and of a Berty CSTR reactor system 98% complete. Three unsupported iron catalysts and three alumina-supported iron catalysts were prepared by nonaqueous-evaporative deposition (NED) or aqueous impregnation (AI) and characterized by chemisorption, BET, extent-of-reduction, XRD, and TEM methods. These catalysts, covering a wide range of dispersions and metal loadings, are well-reduced and relatively thermally stable up to 500-600 C in H{sub 2}, thus ideal for kinetic and mechanistic studies. The alumina-supported iron catalysts will be used for kinetic and mechanistic studies. In the coming year, adsorption/desorption properties, rates of elementary steps, and global reaction rates will be measured for these catalysts, with and without promoters, providing a database for understanding effects of dispersion, metal loading, and support on elementary kinetic parameters and for validation of computational models that incorporate effects of surface structure and promoters. Furthermore, using state-of-the-art self-consistent Density Functional Theory (DFT) methods, we have extensively studied the thermochemistry and kinetics of various elementary steps on three different model surfaces: (1) Fe(110), (2) Fe(110) modified by subsurface C, and (3) Fe surface modified with Pt adatoms. These studies have yielded valuable insights into the reactivity of Fe surfaces for FTS, and provided accurate estimates for the effect of Fe modifiers such as subsurface C and surface Pt.« less

Thermodynamic and dynamic behaviors of self-organizing polymeric systems

NASA Astrophysics Data System (ADS)

Zhao, Yiqiang

Two topics of self-organizing polymeric systems are explored in this work: thermodynamic and dynamic properties of liquid crystal polymers in solutions and rheological behaviors of self-organizing gels. For dilute nematic solutions of end-on side-chain liquid crystal polysiloxanes (SCLCP) dissolved in 5CB, the chain anisotropies R∥/R ⊥, obtained from electrorheological(ER) analysis based on the Brochard model, are consistent with independent measurements of Rg∥/R g⊥ via small-angle neutron scattering (SANS), which unambiguously demonstrating a slightly prolate SCLCP chain conformation. Dissolution of this prolate SCLCP in flow-aligning 5CB produces a tumbling flow, clearly indicating a discrepancy with the Brochard hydrodynamic theory which predicts such a transition only for oblate conformation. A numerical comparison using a modified version of the Brochard model leads to improved self-consistent agreement between SANS, ER and shear transient experiments. The molecular weight dependence of the chain conformational relaxation time it indicates an extended SCLCP chain conformation in 5CB. SANS analysis suggests that the SCLCP conformation is sensitive to the solvent interaction, i.e. a more extended conformation is observed in isotropic acetone-d6 than in nematic 5CB. A SANS conformational study of SCLCCs with methoxyphenylbenzoate mesogenic side group in CDC13 demonstrates that the form factor of a single comb-like SCLCP chain is well described by a wormlike chain model with finite cross-sectional thickness over the entire q range, taking into account the molecular weight polydispersity. Consistent with measurement of a large R g from low q analysis, the resulting persistence length lp is in the range 28˜32 A, substantially larger than that of unsubstituted polydimethylsiloxane (PDMS) chain (l p =5.8 A), which suggests a relatively rigid SCLCP chain due to the influence of densely attached mesogenic groups. For nematic mixtures of copolysiloxane SCLCP in 5OCB, a metastably extended miscible nematic range is observed at low SCLCP concentration upon cooling. Onset of an induced smectic phase occurs upon cooling at 60%wt SCLCP concentration which corresponds to 48:52 molar ratio of mesogens. Dielectric spectra of these mixtures over a wide concentration range exhibit two distinct regimes of relaxation behavior reflecting the crossover from dilute and semidilute to concentrated regime. Rheological behavior of a metallo-supramolecular gel with thixotropic feature is explored to understand the viscoelastic behaviors of self-assembling networks consisting of "living polymers". A well-defined yield point and non-linear viscoelasticity at small strain are probed via the controlled-stress and controlled-strain measurements, respectively. The self-assembled network is readily presheared into a Newtonian sol and displays a three-stage kinetic recovery process, closely associated with the metal ion-ligand binding kinetics and related phase behaviors. Finally, we investigate the viscoelastic properties of a novel colloidal gel in which macrocycles self-assemble into interconnected self-organized clusters. A series of rheological experiments are combined to reveal the shear responsive nature as well as linear and nonlinear viscoelasticity of this gel. Certain features of observed viscoelastic properties demonstrate the characteristics of the behavior of colloidal gels which show slow glassy dynamics. The negative temperature dependence of the storage modulus at low frequency suggests that enthalpic contributions to elasticity need to be considered, presumably due to internal energy changes upon deformation.

Thermodynamics versus Kinetics Dichotomy in the Linear Self-Assembly of Mixed Nanoblocks.

PubMed

Ruiz, L; Keten, S

2014-06-05

We report classical and replica exchange molecular dynamics simulations that establish the mechanisms underpinning the growth kinetics of a binary mix of nanorings that form striped nanotubes via self-assembly. A step-growth coalescence model captures the growth process of the nanotubes, which suggests that high aspect ratio nanostructures can grow by obeying the universal laws of self-similar coarsening, contrary to systems that grow through nucleation and elongation. Notably, striped patterns do not depend on specific growth mechanisms, but are governed by tempering conditions that control the likelihood of depropagation and fragmentation.

Kinetic Simulations of Dense Plasma Focus Breakdown

NASA Astrophysics Data System (ADS)

Schmidt, A.; Higginson, D. P.; Jiang, S.; Link, A.; Povilus, A.; Sears, J.; Bennett, N.; Rose, D. V.; Welch, D. R.

2015-11-01

A dense plasma focus (DPF) device is a type of plasma gun that drives current through a set of coaxial electrodes to assemble gas inside the device and then implode that gas on axis to form a Z-pinch. This implosion drives hydrodynamic and kinetic instabilities that generate strong electric fields, which produces a short intense pulse of x-rays, high-energy (>100 keV) electrons and ions, and (in deuterium gas) neutrons. A strong factor in pinch performance is the initial breakdown and ionization of the gas along the insulator surface separating the two electrodes. The smoothness and isotropy of this ionized sheath are imprinted on the current sheath that travels along the electrodes, thus making it an important portion of the DPF to both understand and optimize. Here we use kinetic simulations in the Particle-in-cell code LSP to model the breakdown. Simulations are initiated with neutral gas and the breakdown modeled self-consistently as driven by a charged capacitor system. We also investigate novel geometries for the insulator and electrodes to attempt to control the electric field profile. The initial ionization fraction of gas is explored computationally to gauge possible advantages of pre-ionization which could be created experimentally via lasers or a glow-discharge. Prepared by LLNL under Contract DE-AC52-07NA27344.

1D kinetic simulations of a short glow discharge in helium

NASA Astrophysics Data System (ADS)

Yuan, Chengxun; Bogdanov, E. A.; Eliseev, S. I.; Kudryavtsev, A. A.

2017-07-01

This paper presents a 1D model of a direct current glow discharge based on the solution of the kinetic Boltzmann equation in the two-term approximation. The model takes into account electron-electron coulomb collisions, the corresponding collision integral is written in both detailed and simplified forms. The Boltzmann equation for electrons is coupled with continuity equations for ions and metastable atoms and the Poisson equation for electric potential. Simulations are carried out self-consistently for the whole length of discharge in helium (from cathode to anode) for cases p = 1 Torr, L = 3.6 cm and p = 20 Torr, L = 1.8 mm, so that pL = 3.6 cm.Torr in both cases. It is shown that simulations based on the kinetic approach give lower values of electron temperature in plasma than fluid simulations. Peaks in spatial differential flux corresponding to the electrons originating from superelastic collisions and Penning ionization were observed in simulations. Different approaches of taking coulomb collisions into account give significantly different values of electron density and electron temperature in plasma. Analysis showed that using a simplified approach gives a non-zero contribution to the electron energy balance, which is comparable to energy losses on elastic and inelastic collisions and leads to significant errors and thus is not recommended.

Self-consistent implementation of meta-GGA functionals for the ONETEP linear-scaling electronic structure package.

PubMed

Womack, James C; Mardirossian, Narbe; Head-Gordon, Martin; Skylaris, Chris-Kriton

2016-11-28

Accurate and computationally efficient exchange-correlation functionals are critical to the successful application of linear-scaling density functional theory (DFT). Local and semi-local functionals of the density are naturally compatible with linear-scaling approaches, having a general form which assumes the locality of electronic interactions and which can be efficiently evaluated by numerical quadrature. Presently, the most sophisticated and flexible semi-local functionals are members of the meta-generalized-gradient approximation (meta-GGA) family, and depend upon the kinetic energy density, τ, in addition to the charge density and its gradient. In order to extend the theoretical and computational advantages of τ-dependent meta-GGA functionals to large-scale DFT calculations on thousands of atoms, we have implemented support for τ-dependent meta-GGA functionals in the ONETEP program. In this paper we lay out the theoretical innovations necessary to implement τ-dependent meta-GGA functionals within ONETEP's linear-scaling formalism. We present expressions for the gradient of the τ-dependent exchange-correlation energy, necessary for direct energy minimization. We also derive the forms of the τ-dependent exchange-correlation potential and kinetic energy density in terms of the strictly localized, self-consistently optimized orbitals used by ONETEP. To validate the numerical accuracy of our self-consistent meta-GGA implementation, we performed calculations using the B97M-V and PKZB meta-GGAs on a variety of small molecules. Using only a minimal basis set of self-consistently optimized local orbitals, we obtain energies in excellent agreement with large basis set calculations performed using other codes. Finally, to establish the linear-scaling computational cost and applicability of our approach to large-scale calculations, we present the outcome of self-consistent meta-GGA calculations on amyloid fibrils of increasing size, up to tens of thousands of atoms.

Self-consistent implementation of meta-GGA functionals for the ONETEP linear-scaling electronic structure package

NASA Astrophysics Data System (ADS)

Womack, James C.; Mardirossian, Narbe; Head-Gordon, Martin; Skylaris, Chris-Kriton

2016-11-01

Accurate and computationally efficient exchange-correlation functionals are critical to the successful application of linear-scaling density functional theory (DFT). Local and semi-local functionals of the density are naturally compatible with linear-scaling approaches, having a general form which assumes the locality of electronic interactions and which can be efficiently evaluated by numerical quadrature. Presently, the most sophisticated and flexible semi-local functionals are members of the meta-generalized-gradient approximation (meta-GGA) family, and depend upon the kinetic energy density, τ, in addition to the charge density and its gradient. In order to extend the theoretical and computational advantages of τ-dependent meta-GGA functionals to large-scale DFT calculations on thousands of atoms, we have implemented support for τ-dependent meta-GGA functionals in the ONETEP program. In this paper we lay out the theoretical innovations necessary to implement τ-dependent meta-GGA functionals within ONETEP's linear-scaling formalism. We present expressions for the gradient of the τ-dependent exchange-correlation energy, necessary for direct energy minimization. We also derive the forms of the τ-dependent exchange-correlation potential and kinetic energy density in terms of the strictly localized, self-consistently optimized orbitals used by ONETEP. To validate the numerical accuracy of our self-consistent meta-GGA implementation, we performed calculations using the B97M-V and PKZB meta-GGAs on a variety of small molecules. Using only a minimal basis set of self-consistently optimized local orbitals, we obtain energies in excellent agreement with large basis set calculations performed using other codes. Finally, to establish the linear-scaling computational cost and applicability of our approach to large-scale calculations, we present the outcome of self-consistent meta-GGA calculations on amyloid fibrils of increasing size, up to tens of thousands of atoms.

Cometabolic degradation kinetics of TCE and phenol by Pseudomonas putida.

PubMed

Chen, Yan-Min; Lin, Tsair-Fuh; Huang, Chih; Lin, Jui-Che

2008-08-01

Modeling of cometabolic kinetics is important for better understanding of degradation reaction and in situ application of bio-remediation. In this study, a model incorporated cell growth and decay, loss of transformation activity, competitive inhibition between growth substrate and non-growth substrate and self-inhibition of non-growth substrate was proposed to simulate the degradation kinetics of phenol and trichloroethylene (TCE) by Pseudomonas putida. All the intrinsic parameters employed in this study were measured independently, and were then used for predicting the batch experimental data. The model predictions conformed well to the observed data at different phenol and TCE concentrations. At low TCE concentrations (<2 mg l(-1)), the models with or without self-inhibition of non-growth substrate both simulated the experimental data well. However, at higher TCE concentrations (>6 mg l(-1)), only the model considering self-inhibition can describe the experimental data, suggesting that a self-inhibition of TCE was present in the system. The proposed model was also employed in predicting the experimental data conducted in a repeated batch reactor, and good agreements were observed between model predictions and experimental data. The results also indicated that the biomass loss in the degradation of TCE below 2 mg l(-1) can be totally recovered in the absence of TCE for the next cycle, and it could be used for the next batch experiment for the degradation of phenol and TCE. However, for higher concentration of TCE (>6 mg l(-1)), the recovery of biomass may not be as good as that at lower TCE concentrations.

Active Galactic Nucleus Feedback in an Isolated Elliptical Galaxy: The Effect of Strong Radiative Feedback in the Kinetic Mode

NASA Astrophysics Data System (ADS)

Gan, Zhaoming; Yuan, Feng; Ostriker, Jeremiah P.; Ciotti, Luca; Novak, Gregory S.

2014-07-01

Based on two-dimensional high-resolution hydrodynamic numerical simulation, we study the mechanical and radiative feedback effects from the central active galactic nucleus (AGN) on the cosmological evolution of an isolated elliptical galaxy. The inner boundary of the simulation domain is carefully chosen so that the fiducial Bondi radius is resolved and the accretion rate of the black hole is determined self-consistently. It is well known that when the accretion rates are high and low, the central AGNs will be in cold and hot accretion modes, which correspond to the radiative and kinetic feedback modes, respectively. The emitted spectrum from the hot accretion flows is harder than that from the cold accretion flows, which could result in a higher Compton temperature accompanied by a more efficient radiative heating, according to previous theoretical works. Such a difference of the Compton temperature between the two feedback modes, the focus of this study, has been neglected in previous works. Significant differences in the kinetic feedback mode are found as a result of the stronger Compton heating. More importantly, if we constrain models to correctly predict black hole growth and AGN duty cycle after cosmological evolution, we find that the favored model parameters are constrained: mechanical feedback efficiency diminishes with decreasing luminosity (the maximum efficiency being ~= 10-3.5), and X-ray Compton temperature increases with decreasing luminosity, although models with fixed mechanical efficiency and Compton temperature can be found that are satisfactory as well. We conclude that radiative feedback in the kinetic mode is much more important than previously thought.

Growing Actin Networks Form Lamellipodium and Lamellum by Self-Assembly

PubMed Central

Huber, Florian; Käs, Josef; Stuhrmann, Björn

2008-01-01

Many different cell types are able to migrate by formation of a thin actin-based cytoskeletal extension. Recently, it became evident that this extension consists of two distinct substructures, designated lamellipodium and lamellum, which differ significantly in their kinetic and kinematic properties as well as their biochemical composition. We developed a stochastic two-dimensional computer simulation that includes chemical reaction kinetics, G-actin diffusion, and filament transport to investigate the formation of growing actin networks in migrating cells. Model parameters were chosen based on experimental data or theoretical considerations. In this work, we demonstrate the system's ability to form two distinct networks by self-organization. We found a characteristic transition in mean filament length as well as a distinct maximum in depolymerization flux, both within the first 1–2 μm. The separation into two distinct substructures was found to be extremely robust with respect to initial conditions and variation of model parameters. We quantitatively investigated the complex interplay between ADF/cofilin and tropomyosin and propose a plausible mechanism that leads to spatial separation of, respectively, ADF/cofilin- or tropomyosin-dominated compartments. Tropomyosin was found to play an important role in stabilizing the lamellar actin network. Furthermore, the influence of filament severing and annealing on the network properties is explored, and simulation data are compared to existing experimental data. PMID:18708450

Power spectra in warm G inflation and its consistency: Stochastic approach

NASA Astrophysics Data System (ADS)

Motaharfar, Meysam; Massaeli, Erfan; Sepangi, Hamid Reza

2017-11-01

Recently, it has been realized that the so-called G -inflation model inspired by supplementing a generalized covariant Galileon-like nonlinear derivative self-interacting term to the standard kinetic term should be ruled out from inflationary models. This is due to the fact that it suffers from lack of an oscillatory phase at the end of the inflationary regime which is typically accompanied by the appearance of a negative squared propagation speed of the scalar mode leading to instabilities of small-scale perturbations. In this regard, the warm G -inflation scenario is proposed where for G inflation to survive, the Galileon scalar field is coupled to the radiation field through a dissipation term which results in removing the reheating period due to the characteristics of warm inflationary scenario. In so doing, a linear stability analysis is first performed to obtain the appropriate slow-roll conditions in such a proposal. Cosmological perturbations of the model are then investigated by utilizing a fluctuation-dissipation theorem and analytical expressions are derived for observable quantities; the power spectrum, tilt spectral index, and tensor-to-scalar ratio in terms of P S R parameters and Galileon flow functions. Finally, the model is solved for chaotic self-interacting potentials, particularly the renormalizable Higgs potential λ/4 ϕ4, and shown to be consistent with observations in the weak dissipation Q ≪1 +3 δ/GXδX and G -dominant 3 δ/GXδX≫1 regime despite its large self-coupling, since the energy scale at the horizon crossing is depressed by the synergy of Galileon and thermal effects.

KINETIC EVOLUTION OF CORONAL HOLE PROTONS BY IMBALANCED ION-CYCLOTRON WAVES: IMPLICATIONS FOR MEASUREMENTS BY SOLAR PROBE PLUS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Isenberg, Philip A.; Vasquez, Bernard J.

We extend the kinetic guiding-center model of collisionless coronal hole protons presented in Isenberg and Vasquez to consider driving by imbalanced spectra of obliquely propagating ion-cyclotron waves. These waves are assumed to be a small by-product of the imbalanced turbulent cascade to high perpendicular wavenumber, and their total intensity is taken to be 1% of the total fluctuation energy. We also extend the kinetic solutions for the proton distribution function in the resulting fast solar wind to heliocentric distances of 20 solar radii, which will be attainable by the Solar Probe Plus spacecraft. We consider three ratios of outward-propagating tomore » inward-propagating resonant intensities: 1, 4, and 9. The self-consistent bulk flow speed reaches fast solar wind values in all cases, and these speeds are basically independent of the intensity ratio. The steady-state proton distribution is highly organized into nested constant-density shells by the resonant wave-particle interaction. The radial evolution of this kinetic distribution as the coronal hole plasma flows outward is understood as a competition between the inward- and outward-directed large-scale forces, causing an effective circulation of particles through the (v{sub ∥}, v{sub ⊥}) phase space and a characteristic asymmetric shape to the distribution. These asymmetries are substantial and persist to the outer limit of the model computation, where they should be observable by the Solar Probe Plus instruments.« less

DOE Office of Scientific and Technical Information (OSTI.GOV)

Le, Hai P.; Cambier, Jean -Luc

Here, we present a numerical model and a set of conservative algorithms for Non-Maxwellian plasma kinetics with inelastic collisions. These algorithms self-consistently solve for the time evolution of an isotropic electron energy distribution function interacting with an atomic state distribution function of an arbitrary number of levels through collisional excitation, deexcitation, as well as ionization and recombination. Electron-electron collisions, responsible for thermalization of the electron distribution, are also included in the model. The proposed algorithms guarantee mass/charge and energy conservation in a single step, and is applied to the case of non-uniform gridding of the energy axis in the phasemore » space of the electron distribution function. Numerical test cases are shown to demonstrate the accuracy of the method and its conservation properties.« less

A three-dimensional phase field model for nanowire growth by the vapor-liquid-solid mechanism

NASA Astrophysics Data System (ADS)

Wang, Yanming; Ryu, Seunghwa; McIntyre, Paul C.; Cai, Wei

2014-07-01

We present a three-dimensional multi-phase field model for catalyzed nanowire (NW) growth by the vapor-liquid-solid (VLS) mechanism. The equation of motion contains both a Ginzburg-Landau term for deposition and a diffusion (Cahn-Hilliard) term for interface relaxation without deposition. Direct deposition from vapor to solid, which competes with NW crystal growth through the molten catalyst droplet, is suppressed by assigning a very small kinetic coefficient at the solid-vapor interface. The thermodynamic self-consistency of the model is demonstrated by its ability to reproduce the equilibrium contact angles at the VLS junction. The incorporation of orientation dependent gradient energy leads to faceting of the solid-liquid and solid-vapor interfaces. The model successfully captures the curved shape of the NW base and the Gibbs-Thomson effect on growth velocity.

A physics-based crystallographic modeling framework for describing the thermal creep behavior of Fe-Cr alloys

DOE PAGES

Wen, Wei; Capolungo, Laurent; Patra, Anirban; ...

2017-02-23

In this work, a physics-based thermal creep model is developed based on the understanding of the microstructure in Fe-Cr alloys. This model is associated with a transition state theory based framework that considers the distribution of internal stresses at sub-material point level. The thermally activated dislocation glide and climb mechanisms are coupled in the obstacle-bypass processes for both dislocation and precipitate-type barriers. A kinetic law is proposed to track the dislocation densities evolution in the subgrain interior and in the cell wall. The predicted results show that this model, embedded in the visco-plastic self-consistent (VPSC) framework, captures well the creepmore » behaviors for primary and steady-state stages under various loading conditions. We also discuss the roles of the mechanisms involved.« less

Distribution of nitrogen and defects in SiO(x)N(y)/Si structures formed by the thermal nitridation of SiO2/Si

NASA Technical Reports Server (NTRS)

Vasquez, R. P.; Madhukar, A.; Grunthaner, F. J.; Naiman, M. L.

1986-01-01

Previously reported nitrogen distributions in SiO2 films on Si which have been thermally nitrided at 1000 C have been explained by a kinetic model of the nitridation process which rests upon the effects of interfacial strain. A critical test of this kinetic model is the validity of the predictions regarding nitrogen distributions obtained at other nitridation temperatures. In this work, nitrogen distributions determined via X-ray photoelectron spectroscopy are reported for samples nitrided at 800 and 1150 C, and are shown to be consistent with the kinetic model. In addition, the intensity of a fluorine marker is found to correlate with the nitrogen distribution, and is postulated to be related to kinetically generated defects in the dielectric film, consistent with the strain-dependent energy of formation of defects proposed recently to explain electrical data.

Understanding the Pulsar High Energy Emission: Macroscopic and Kinetic Models

NASA Astrophysics Data System (ADS)

Kalapotharakos, Constantinos; Brambilla, Gabriele; Timokhin, Andrey; Kust Harding, Alice; Kazanas, Demos

2017-08-01

Pulsars are extraordinary objects powered by the rotation of magnetic fields of order 10^8, 10^12G anchored onto neutron stars and rotating with periods 10^(-3)-10s. These fields mediate the conversion of their rotational energy into MHD winds and at the same time accelerate particles to energies sufficiently high to produce GeV photons. Fermi, since its launch in 2008, has established several trends among the observed gamma-ray pulsar properties playing a catalytic role in the current modeling of the high energy emission in pulsar magnetospheres. We judiciously use the guidance provided by the Fermi data to yield meaningful constraints on the macroscopic parameters of our global dissipative pulsar magnetosphere models. Our FIDO (Force-Free Inside, Dissipative Outside) models indicate that the dissipative regions lie outside the light cylinder near the equatorial current sheet. Our models reproduce the light-curve phenomenology while a detailed comparison of the model spectral properties with those observed by Fermi reveals the dependence of the macroscopic conductivity parameter on the spin-down rate providing a unique insight into the understanding of the physical mechanisms behind the high-energy emission in pulsar magnetospheres. Finally, we further exploit these important results by building self-consistent 3D global kinetic particle-in-cell (PIC) models which, eventually, provide the dependence of the macroscopic parameter behavior (e.g. conductivity) on the microphysical properties (e.g. particle multiplicities, particle injection rates). Our PIC models provide field structures and particle distributions that are not only consistent with each other but also able to reproduce a broad range of the observed gamma-ray phenomenology (light curves and spectral properties) of both young and millisecond pulsars.

Nucleation and growth in cluster dynamics: A quantitative test of the classical kinetic approach

NASA Astrophysics Data System (ADS)

Gránásy, László; James, Peter F.

2000-12-01

Nucleation and size dependent growth of nanometer sized crystalline particles in glassy media have been studied by numerically solving the Turnbull-Fisher master equations that describe the time evolution of cluster population. Time dependencies of the formation rate and number density are determined for large clusters (built of up to 2×105 formula units, containing 1.8×106 atoms). We demonstrate that the formation rate and number density of such clusters are well approximated by Shneidman's asymptotically exact analytical solution. A quantitative test of the kinetic Turnbull-Fisher model has been performed: Evaluating the kinetic coefficients and interfacial parameters from the transient time and steady-state nucleation rates measured on six stoichiometric oxide glass compositions (lithium-disilicate, barium-disilicate, lithium-diborate, wollastonite, 1:2:3 and 2:1:3 soda-lime-silica glass compositions), we calculated the macroscopic growth rates and compared with experiments. For wollastonite, lithium-diborate and the 1:2:3 soda-lime-silica glass, differences of 2 to 4 orders of magnitude have been observed between theory and experiment. This inadequacy of the microscopic kinetic parameters in describing macroscopic growth cannot be explained by either the curvature effect on the interfacial free energy or the self-consistency correction for the cluster free energy. The origin of the discrepancy is discussed.

Operational advances in ring current modeling using RAM-SCB

DOE Office of Scientific and Technical Information (OSTI.GOV)

Welling, Daniel T; Jordanova, Vania K; Zaharia, Sorin G

The Ring current Atmosphere interaction Model with Self-Consistently calculated 3D Magnetic field (RAM-SCB) combines a kinetic model of the ring current with a force-balanced model of the magnetospheric magnetic field to create an inner magnetospheric model that is magnetically self consistent. RAM-SCB produces a wealth of outputs that are valuable to space weather applications. For example, the anisotropic particle distribution of the KeV-energy population calculated by the code is key for predicting surface charging on spacecraft. Furthermore, radiation belt codes stand to benefit substantially from RAM-SCB calculated magnetic field values and plasma wave growth rates - both important for determiningmore » the evolution of relativistic electron populations. RAM-SCB is undergoing development to bring these benefits to the space weather community. Data-model validation efforts are underway to assess the performance of the system. 'Virtual Satellite' capability has been added to yield satellite-specific particle distribution and magnetic field output. The code's outer boundary is being expanded to 10 Earth Radii to encompass previously neglected geosynchronous orbits and allow the code to be driven completely by either empirical or first-principles based inputs. These advances are culminating towards a new, real-time version of the code, rtRAM-SCB, that can monitor the inner magnetosphere conditions on both a global and spacecraft-specific level. This paper summarizes these new features as well as the benefits they provide the space weather community.« less

Operational Advances in Ring Current Modeling Using RAM-SCB

NASA Astrophysics Data System (ADS)

Morley, S.; Welling, D. T.; Zaharia, S. G.; Jordanova, V. K.

2010-12-01

The Ring current Atmosphere interaction Model with Self-Consistently calculated 3D Magnetic field (RAM-SCB) combines a kinetic model of the ring current with a force-balanced model of the magnetospheric magnetic field to create an inner magnetospheric model that is magnetically self consistent. RAM-SCB produces a wealth of outputs that are valuable to space weather applications. For example, the anisotropic particle distribution of the KeV-energy population calculated by the code is key for predicting surface charging on spacecraft. Furthermore, radiation belt codes stand to benefit substantially from RAM-SCB calculated magnetic field values and plasma wave growth rates - both important for determining the evolution of relativistic electron populations. RAM-SCB is undergoing development to bring these benefits to the space weather community. Data-model validation efforts are underway to assess the performance of the system. “Virtual Satellite” capability has been added to yield satellite-specific particle distribution and magnetic field output. The code’s outer boundary is being expanded to 10 Earth Radii to encompass previously neglected geosynchronous orbits and allow the code to be driven completely by either empirical or first-principles based inputs. These advances are culminating towards a new, real-time version of the code, rtRAM-SCB, that can monitor the inner magnetosphere conditions on both a global and spacecraft-specific level. This paper summarizes these new features as well as the benefits they provide the space weather community.

Cluster dynamics and cluster size distributions in systems of self-propelled particles

NASA Astrophysics Data System (ADS)

Peruani, F.; Schimansky-Geier, L.; Bär, M.

2010-12-01

Systems of self-propelled particles (SPP) interacting by a velocity alignment mechanism in the presence of noise exhibit rich clustering dynamics. Often, clusters are responsible for the distribution of (local) information in these systems. Here, we investigate the properties of individual clusters in SPP systems, in particular the asymmetric spreading behavior of clusters with respect to their direction of motion. In addition, we formulate a Smoluchowski-type kinetic model to describe the evolution of the cluster size distribution (CSD). This model predicts the emergence of steady-state CSDs in SPP systems. We test our theoretical predictions in simulations of SPP with nematic interactions and find that our simple kinetic model reproduces qualitatively the transition to aggregation observed in simulations.

Transport properties in dilute UN (X ) solid solutions (X =Xe ,Kr )

NASA Astrophysics Data System (ADS)

Claisse, Antoine; Schuler, Thomas; Lopes, Denise Adorno; Olsson, Pär

2016-11-01

Uranium nitride (UN) is a candidate fuel for current GEN III fission reactors, for which it is investigated as an accident-tolerant fuel, as well as for future GEN IV reactors. In this study, we investigate the kinetic properties of gas fission products (Xe and Kr) in UN. Binding and migration energies are obtained using density functional theory, with an added Hubbard correlation to model f electrons, and the occupation matrix control scheme to avoid metastable states. These energies are then used as input for the self-consistent mean field method which enables to determine transport coefficients for vacancy-mediated diffusion of Xe and Kr on the U sublattice. The magnetic ordering of the UN structure is explicitly taken into account, for both energetic and transport properties. Solute diffusivities are compared with experimental measurements and the effect of various parameters on the theoretical model is carefully investigated. We find that kinetic correlations are very strong in this system, and that despite atomic migration anisotropy, macroscopic solute diffusivities show limited anisotropy. Our model indicates that the discrepancy between experimental measurements probably results from different irradiation conditions, and hence different defect concentrations.

Simulating complex atomistic processes: On-the-fly kinetic Monte Carlo scheme with selective active volumes

NASA Astrophysics Data System (ADS)

Xu, Haixuan; Osetsky, Yury N.; Stoller, Roger E.

2011-10-01

An accelerated atomistic kinetic Monte Carlo (KMC) approach for evolving complex atomistic structures has been developed. The method incorporates on-the-fly calculations of transition states (TSs) with a scheme for defining active volumes (AVs) in an off-lattice (relaxed) system. In contrast to conventional KMC models that require all reactions to be predetermined, this approach is self-evolving and any physically relevant motion or reaction may occur. Application of this self-evolving atomistic kinetic Monte Carlo (SEAK-MC) approach is illustrated by predicting the evolution of a complex defect configuration obtained in a molecular dynamics (MD) simulation of a displacement cascade in Fe. Over much longer times, it was shown that interstitial clusters interacting with other defects may change their structure, e.g., from glissile to sessile configuration. The direct comparison with MD modeling confirms the atomistic fidelity of the approach, while the longer time simulation demonstrates the unique capability of the model.

Kinetic neoclassical transport in the H-mode pedestal

DOE PAGES

Battaglia, D. J.; Burrell, K. H.; Chang, C. S.; ...

2014-07-16

Multi-species kinetic neoclassical transport through the QH-mode pedestal and scrapeoff layer on DIII-D is calculated using XGC0, a 5D full-f particle-in-cell drift-kinetic solver with self-consistent neutral recycling and sheath potentials. We achieved quantitative agreement between the fluxdriven simulation and the experimental electron density, impurity density and orthogonal measurements of impurity temperature and flow profiles by adding random-walk particle diffusion to the guiding-center drift motion. Furthermore, we computed the radial electric field (Er) that maintains ambipolar transport across flux surfaces and to the wall self-consistently on closed and open magnetic field lines, and is in excellent agreement with experiment. The Ermore » inside the separatrix is the unique solution that balances the outward flux of thermal tail deuterium ions against the outward neoclassical electron flux and inward pinch of impurity and colder deuterium ions. Particle transport in the pedestal is primarily due to anomalous transport, while the ion heat and momentum transport is primarily due to the neoclassical transport. The full-f treatment quantifies the non-Maxwellian energy distributions that describe a number of experimental observations in low-collisionallity pedestals on DIII-D, including intrinsic co-Ip parallel flows in the pedestal, ion temperature anisotropy and large impurity temperatures in the scrape-off layer.« less

Monodisperse self-assembly in a model with protein-like interactions

NASA Astrophysics Data System (ADS)

Wilber, Alex W.; Doye, Jonathan P. K.; Louis, Ard A.; Lewis, Anna C. F.

2009-11-01

We study the self-assembly behavior of patchy particles with "proteinlike" interactions that can be considered as a minimal model for the assembly of viral capsids and other shell-like protein complexes. We thoroughly explore the thermodynamics and dynamics of self-assembly as a function of the parameters of the model and find robust assembly of all target structures considered. Optimal assembly occurs in the region of parameter space where a free energy barrier regulates the rate of nucleation, thus preventing the premature exhaustion of the supply of monomers that can lead to the formation of incomplete shells. The interactions also need to be specific enough to prevent the assembly of malformed shells, but while maintaining kinetic accessibility. Free energy landscapes computed for our model have a funnel-like topography guiding the system to form the target structure and show that the torsional component of the interparticle interactions prevents the formation of disordered aggregates that would otherwise act as kinetic traps.

Ergodic model for the expansion of spherical nanoplasmas.

PubMed

Peano, F; Coppa, G; Peinetti, F; Mulas, R; Silva, L O

2007-06-01

Recently, the collisionless expansion of spherical nanoplasmas has been analyzed with a new ergodic model, clarifying the transition from hydrodynamiclike to Coulomb-explosion regimes, and providing accurate laws for the relevant features of the phenomenon. A complete derivation of the model is presented here. The important issue of the self-consistent initial conditions is addressed by analyzing the initial charging transient due to the electron expansion, in the approximation of immobile ions. A comparison among different kinetic models for the expansion is presented, showing that the ergodic model provides a simplified description, which retains the essential information on the electron distribution, in particular, the energy spectrum. Results are presented for a wide range of initial conditions (determined from a single dimensionless parameter), in excellent agreement with calculations from the exact Vlasov-Poisson theory, thus providing a complete and detailed characterization of all the stages of the expansion.

Study on the kinetic self-assembly of type I collagen from tilapia (Oreochromis niloticus) skin using the fluorescence probe thioflavin T.

PubMed

Yan, Mingyan; Wang, Xinping

2018-05-27

The kinetic self-assembly of type I collagen from tilapia (Oreochromis niloticus) skin was characterized by the fluorescence method based on thioflavin T (ThT). The fluorescence probe could bind to the active monomeric collagen with a higher ordered degree of molecule, which displayed the pH and ionic strength dependence, the binding constant higher at neutral pH and proportional to the NaCl concentration. Compared to the turbidity method, ThT was more suitable to characterize the nucleation phase of collagen self-assembly. The nucleus size was determined through the ThT fluorescence and linear-polymerization model. At various pH and ionic strength, the nucleus size was nearly identical, either one or two monomers, demonstrating that one or two active monomeric collagen formed into the nucleus and different pH and ionic strength didn't alter the self-assembly mechanism of collagen. This approach was beneficial to advance the understanding of the kinetic self-assembly of the fish-sourced collagen in vitro. Copyright © 2018 Elsevier B.V. All rights reserved.

Extracting surface diffusion coefficients from batch adsorption measurement data: application of the classic Langmuir kinetics model.

PubMed

Chu, Khim Hoong

2017-11-09

Surface diffusion coefficients may be estimated by fitting solutions of a diffusion model to batch kinetic data. For non-linear systems, a numerical solution of the diffusion model's governing equations is generally required. We report here the application of the classic Langmuir kinetics model to extract surface diffusion coefficients from batch kinetic data. The use of the Langmuir kinetics model in lieu of the conventional surface diffusion model allows derivation of an analytical expression. The parameter estimation procedure requires determining the Langmuir rate coefficient from which the pertinent surface diffusion coefficient is calculated. Surface diffusion coefficients within the 10 -9 to 10 -6  cm 2 /s range obtained by fitting the Langmuir kinetics model to experimental kinetic data taken from the literature are found to be consistent with the corresponding values obtained from the traditional surface diffusion model. The virtue of this simplified parameter estimation method is that it reduces the computational complexity as the analytical expression involves only an algebraic equation in closed form which is easily evaluated by spreadsheet computation.

Particle-In-Cell Simulations of the Solar Wind Interaction with Lunar Crustal Magnetic Anomalies: Magnetic Cusp Regions

NASA Technical Reports Server (NTRS)

Poppe, A. R.; Halekas, J. S.; Delory, G. T.; Farrell, W. M.

2012-01-01

As the solar wind is incident upon the lunar surface, it will occasionally encounter lunar crustal remanent magnetic fields. These magnetic fields are small-scale, highly non-dipolar, have strengths up to hundreds of nanotesla, and typically interact with the solar wind in a kinetic fashion. Simulations, theoretical analyses, and spacecraft observations have shown that crustal fields can reflect solar wind protons via a combination of magnetic and electrostatic reflection; however, analyses of surface properties have suggested that protons may still access the lunar surface in the cusp regions of crustal magnetic fields. In this first report from a planned series of studies, we use a 1 1/2-dimensional, electrostatic particle-in-cell code to model the self-consistent interaction between the solar wind, the cusp regions of lunar crustal remanent magnetic fields, and the lunar surface. We describe the self-consistent electrostatic environment within crustal cusp regions and discuss the implications of this work for the role that crustal fields may play regulating space weathering of the lunar surface via proton bombardment.

Multi-Fluid Moment Simulations of Ganymede using the Next-Generation OpenGGCM

NASA Astrophysics Data System (ADS)

Wang, L.; Germaschewski, K.; Hakim, A.; Bhattacharjee, A.; Raeder, J.

2015-12-01

We coupled the multi-fluid moment code Gkeyll[1,2] to the next-generation OpenGGCM[3], and studied the reconnection dynamics at the Ganymede. This work is part of our effort to tackle the grand challenge of integrating kinetic effects into global fluid models. The multi-fluid moment model integrates kinetic effects in that it can capture crucial kinetic physics like pressure tensor effects by evolving moments of the Vlasov equations for each species. This approach has advantages over previous models: desired kinetic effects, together with other important effects like the Hall effect, are self-consistently embedded in the moment equations, and can be efficiently implemented, while not suffering from severe time-step restriction due to plasma oscillation nor artificial whistler modes. This model also handles multiple ion species naturally, which opens up opportunties in investigating the role of oxygen in magnetospheric reconnection and improved coupling to ionosphere models. In this work, the multi-fluid moment solver in Gkeyll was wrapped as a time-stepping module for the high performance, highly flexible next-generation OpenGGCM. Gkeyll is only used to provide the local plasma solver, while computational aspects like parallelization and boundary conditions are handled entirely by OpenGGCM, including interfacing to other models like ionospheric boundary conditions provided by coupling with CTIM [3]. The coupled code is used to study the dynamics near Ganymede, and the results are compared with MHD and Hall MHD results by Dorelli et al. [4]. Hakim, A. (2008). Journal of Fusion Energy, 27, 36-43. Hakim, A., Loverich, J., & Shumlak, U. (2006). Journal of Computational Physics, 219, 418-442. Raeder, J., Larson, D., Li, W., Kepko, E. L., & Fuller-Rowell, T. (2008). Space Science Reviews, 141(1-4), 535-555. Dorelli, J. C., Glocer, A., Collinson, G., & Tóth, G. (2015). Journal of Geophysical Research: Space Physics, 120.

The basis for cosmic ray feedback: Written on the wind

PubMed Central

Zweibel, Ellen G.

2017-01-01

Star formation and supermassive black hole growth in galaxies appear to be self-limiting. The mechanisms for self-regulation are known as feedback. Cosmic rays, the relativistic particle component of interstellar and intergalactic plasma, are among the agents of feedback. Because cosmic rays are virtually collisionless in the plasma environments of interest, their interaction with the ambient medium is primarily mediated by large scale magnetic fields and kinetic scale plasma waves. Because kinetic scales are much smaller than global scales, this interaction is most conveniently described by fluid models. In this paper, I discuss the kinetic theory and the classical theory of cosmic ray hydrodynamics (CCRH) which follows from assuming cosmic rays interact only with self-excited waves. I generalize CCRH to generalized cosmic ray hydrodynamics, which accommodates interactions with extrinsic turbulence, present examples of cosmic ray feedback, and assess where progress is needed. PMID:28579734

The basis for cosmic ray feedback: Written on the wind

NASA Astrophysics Data System (ADS)

Zweibel, Ellen G.

2017-05-01

Star formation and supermassive black hole growth in galaxies appear to be self-limiting. The mechanisms for self-regulation are known as feedback. Cosmic rays, the relativistic particle component of interstellar and intergalactic plasma, are among the agents of feedback. Because cosmic rays are virtually collisionless in the plasma environments of interest, their interaction with the ambient medium is primarily mediated by large scale magnetic fields and kinetic scale plasma waves. Because kinetic scales are much smaller than global scales, this interaction is most conveniently described by fluid models. In this paper, I discuss the kinetic theory and the classical theory of cosmic ray hydrodynamics (CCRH) which follows from assuming cosmic rays interact only with self-excited waves. I generalize CCRH to generalized cosmic ray hydrodynamics, which accommodates interactions with extrinsic turbulence, present examples of cosmic ray feedback, and assess where progress is needed.

The basis for cosmic ray feedback: Written on the wind.

PubMed

Zweibel, Ellen G

2017-05-01

Star formation and supermassive black hole growth in galaxies appear to be self-limiting. The mechanisms for self-regulation are known as feedback . Cosmic rays, the relativistic particle component of interstellar and intergalactic plasma, are among the agents of feedback. Because cosmic rays are virtually collisionless in the plasma environments of interest, their interaction with the ambient medium is primarily mediated by large scale magnetic fields and kinetic scale plasma waves. Because kinetic scales are much smaller than global scales, this interaction is most conveniently described by fluid models. In this paper, I discuss the kinetic theory and the classical theory of cosmic ray hydrodynamics (CCRH) which follows from assuming cosmic rays interact only with self-excited waves. I generalize CCRH to generalized cosmic ray hydrodynamics, which accommodates interactions with extrinsic turbulence, present examples of cosmic ray feedback, and assess where progress is needed.

Sensitivity of equilibrium profile reconstruction to motional Stark effect measurements

NASA Astrophysics Data System (ADS)

Batha, S. H.; Levinton, F. M.; Hirshman, S. P.; Bell, M. G.; Wieland, R. M.

1996-09-01

The magnetic-field pitch-angle profile, gamma p(R) identical to tan-1(Bpol/Btor), is measured on TFTR using a motional Stark effect (MSE) polarimeter. Measured pitch angle profiles, along with kinetic profiles and external magnetic measurements, are used to compute a self-consistent equilibrium using the free-boundary variational moments equilibrium code VMEC. Uncertainties in the q profile due to uncertainties in gamma P(R), magnetic measurements and kinetic measurements are found to be small. Subsequent uncertainties in the VMEC-calculated current density and shear profiles are also small

The Harrison Diffusion Kinetics Regimes in Solute Grain Boundary Diffusion

DOE Office of Scientific and Technical Information (OSTI.GOV)

Belova, Irina; Fiedler, T; Kulkarni, Nagraj S

2012-01-01

Knowledge of the limits of the principal Harrison kinetics regimes (Type-A, B and C) for grain boundary diffusion is very important for the correct analysis of the depth profiles in a tracer diffusion experiment. These regimes for self-diffusion have been extensively studied in the past by making use of the phenomenological Lattice Monte Carlo (LMC) method with the result that the limits are now well established. The relationship of those self-diffusion limits to the corresponding ones for solute diffusion in the presence of solute segregation to the grain boundaries remains unclear. In the present study, the influence of solute segregationmore » on the limits is investigated with the LMC method for the well-known parallel grain boundary slab model by showing the equivalence of two diffusion models. It is shown which diffusion parameters are useful for identifying the limits of the Harrison kinetics regimes for solute grain boundary diffusion. It is also shown how the measured segregation factor from the diffusion experiment in the Harrison Type-B kinetics regime may differ from the global segregation factor.« less

Kinetic Approaches to Shear-Driven Magnetic Reconnection for Multi-Scale Modeling of CME Initiation

NASA Astrophysics Data System (ADS)

Black, C.; Antiochos, S. K.; DeVore, C.; Germaschewski, K.; Karpen, J. T.

2013-12-01

In the standard model for coronal mass ejections (CME) and/or solar flares, the free energy for the event resides in the strongly sheared magnetic field of a filament channel. The pre-eruption force balance, consisting of an upward force due to the magnetic pressure of the sheared field balanced by a downward tension due to overlying un-sheared field, is widely believed to be disrupted by magnetic reconnection. Therefore, understanding initiation of solar explosive phenomena requires a true multi-scale model of reconnection onset driven by the buildup of magnetic shear. While the application of magnetic-field shear is a trivial matter in MHD simulations, it is a significant challenge in a PIC code. The driver must be implemented in a self-consistent manner and with boundary conditions that avoid the generation of waves that destroy the applied shear. In this work, we describe drivers for 2.5D, aperiodic, PIC systems and discuss the implementation of driver-consistent boundary conditions that allow a net electric current to flow through the walls. Preliminary tests of these boundaries with a MHD equilibrium are shown. This work was supported, in part, by the NASA Living With a Star TR&T Program.

The Einstein-Vlasov System/Kinetic Theory.

PubMed

Andréasson, Håkan

2005-01-01

The main purpose of this article is to provide a guide to theorems on global properties of solutions to the Einstein-Vlasov system. This system couples Einstein's equations to a kinetic matter model. Kinetic theory has been an important field of research during several decades in which the main focus has been on nonrelativistic and special relativistic physics, i.e. to model the dynamics of neutral gases, plasmas, and Newtonian self-gravitating systems. In 1990, Rendall and Rein initiated a mathematical study of the Einstein-Vlasov system. Since then many theorems on global properties of solutions to this system have been established. The Vlasov equation describes matter phenomenologically, and it should be stressed that most of the theorems presented in this article are not presently known for other such matter models (i.e. fluid models). This paper gives introductions to kinetic theory in non-curved spacetimes and then the Einstein-Vlasov system is introduced. We believe that a good understanding of kinetic theory in non-curved spacetimes is fundamental to good comprehension of kinetic theory in general relativity.

The Einstein-Vlasov System/Kinetic Theory.

PubMed

Andréasson, Håkan

2002-01-01

The main purpose of this article is to provide a guide to theorems on global properties of solutions to the Einstein-Vlasov system. This system couples Einstein's equations to a kinetic matter model. Kinetic theory has been an important field of research during several decades in which the main focus has been on nonrelativistic and special relativistic physics, i.e. to model the dynamics of neutral gases, plasmas, and Newtonian self-gravitating systems. In 1990, Rendall and Rein initiated a mathematical study of the Einstein-Vlasov system. Since then many theorems on global properties of solutions to this system have been established. The Vlasov equation describes matter phenomenologically, and it should be stressed that most of the theorems presented in this article are not presently known for other such matter models (i.e. fluid models). This paper gives introductions to kinetic theory in non-curved spacetimes and then the Einstein-Vlasov system is introduced. We believe that a good understanding of kinetic theory in non-curved spacetimes is fundamental to good comprehension of kinetic theory in general relativity.

Active galactic nucleus feedback in an isolated elliptical galaxy: The effect of strong radiative feedback in the kinetic mode

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gan, Zhaoming; Yuan, Feng; Ostriker, Jeremiah P.

2014-07-10

Based on two-dimensional high-resolution hydrodynamic numerical simulation, we study the mechanical and radiative feedback effects from the central active galactic nucleus (AGN) on the cosmological evolution of an isolated elliptical galaxy. The inner boundary of the simulation domain is carefully chosen so that the fiducial Bondi radius is resolved and the accretion rate of the black hole is determined self-consistently. It is well known that when the accretion rates are high and low, the central AGNs will be in cold and hot accretion modes, which correspond to the radiative and kinetic feedback modes, respectively. The emitted spectrum from the hotmore » accretion flows is harder than that from the cold accretion flows, which could result in a higher Compton temperature accompanied by a more efficient radiative heating, according to previous theoretical works. Such a difference of the Compton temperature between the two feedback modes, the focus of this study, has been neglected in previous works. Significant differences in the kinetic feedback mode are found as a result of the stronger Compton heating. More importantly, if we constrain models to correctly predict black hole growth and AGN duty cycle after cosmological evolution, we find that the favored model parameters are constrained: mechanical feedback efficiency diminishes with decreasing luminosity (the maximum efficiency being ≅ 10{sup –3.5}), and X-ray Compton temperature increases with decreasing luminosity, although models with fixed mechanical efficiency and Compton temperature can be found that are satisfactory as well. We conclude that radiative feedback in the kinetic mode is much more important than previously thought.« less

Self-consistent particle-in-cell simulations of fundamental and harmonic radio plasma emission mechanisms

NASA Astrophysics Data System (ADS)

Tsiklauri, D.; Thurgood, J. O.

2015-12-01

first co-author Jonathan O. Thurgood (QMUL) The simulation of three-wave interaction based plasma emission, an underlying mechanism for type III solar radio bursts, is a challenging task requiring fully-kinetic, multi-dimensional models. This paper aims to resolve a contradiction in past attempts, whereby some authors report that no such processes occur and others draw conflicting conclusions, by using 2D, fully kinetic, particle-in-cell simulations of relaxing electron beams. Here we present the results of particle-in-cell simulations which for different physical parameters permit or prohibit the plasma emission. We show that the possibility of plasma emission is contingent upon the frequency of the initial electrostatic waves generated by the bump-in-tail instability, and that these waves may be prohibited from participating in the necessary three-wave interactions due to the frequency beat requirements. We caution against simulating astrophysical radio bursts using unrealistically dense beams (a common approach which reduces run time), as the resulting non-Langmuir characteristics of the initial wave modes significantly suppresses the emission. Comparison of our results indicates that, contrary to the suggestions of previous authors, a plasma emission mechanism based on two counter-propagating beams is unnecessary in astrophysical context. Finally, we also consider the action of the Weibel instability, which generates an electromagnetic beam mode. As this provides a stronger contribution to electromagnetic energy than the emission, we stress that evidence of plasma emission in simulations must disentangle the two contributions and not simply interpret changes in total electromagnetic energy as the evidence of plasma emission. In summary, we present the first self-consistent demonstration of fundamental and harmonic plasma emission from a single-beam system via fully kinetic numerical simulation. Pre-print can be found at http://astro.qmul.ac.uk/~tsiklauri/jtdt1.pdf

The Einstein-Vlasov System/Kinetic Theory.

PubMed

Andréasson, Håkan

2011-01-01

The main purpose of this article is to provide a guide to theorems on global properties of solutions to the Einstein-Vlasov system. This system couples Einstein's equations to a kinetic matter model. Kinetic theory has been an important field of research during several decades in which the main focus has been on non-relativistic and special relativistic physics, i.e., to model the dynamics of neutral gases, plasmas, and Newtonian self-gravitating systems. In 1990, Rendall and Rein initiated a mathematical study of the Einstein-Vlasov system. Since then many theorems on global properties of solutions to this system have been established. This paper gives introductions to kinetic theory in non-curved spacetimes and then the Einstein-Vlasov system is introduced. We believe that a good understanding of kinetic theory in non-curved spacetimes is fundamental to a good comprehension of kinetic theory in general relativity.

The alternative strategy for designing covalent drugs through kinetic effects of pi-stacking on the self-assembled nanoparticles: a model study with antibiotics

NASA Astrophysics Data System (ADS)

Du, Libo; Suo, Siqingaowa; Zhang, Han; Jia, Hongying; Liu, Ke Jian; Zhang, Xue Ji; Liu, Yang

2016-11-01

It is still a huge challenge to find a new strategy for rationally designing covalent drugs because most of them are discovered by serendipity. Considering that the effect of covalent drugs is closely associated with the kinetics of the reaction between drug molecule and its target protein, here we first demonstrate an example of the kinetic effect of pi-stacking of drug molecules on covalent antimicrobial drug design. When PEGylated 7-aminocephalosporanic acid (PEG-ACA) is used as a substrate drug, pi-stacking of the ACA group via the self-assembly of PEG-ACA on the surface of gold nanoparticles (i.e. Au@ACA) exhibits antibacterial activity against E. coli fourfold higher than a PEG-ACA monomer does. The reason can be reasonably attributed to the kinetic rate enhancement for the covalent reaction between Au@ACA and penicillin binding proteins. We believe that the self-assembly of functional groups onto the surface of gold nanoparticles represents a new strategy for covalent drug design.

A new energy transfer model for turbulent free shear flow

NASA Technical Reports Server (NTRS)

Liou, William W.-W.

1992-01-01

A new model for the energy transfer mechanism in the large-scale turbulent kinetic energy equation is proposed. An estimate of the characteristic length scale of the energy containing large structures is obtained from the wavelength associated with the structures predicted by a weakly nonlinear analysis for turbulent free shear flows. With the inclusion of the proposed energy transfer model, the weakly nonlinear wave models for the turbulent large-scale structures are self-contained and are likely to be independent flow geometries. The model is tested against a plane mixing layer. Reasonably good agreement is achieved. Finally, it is shown by using the Liapunov function method, the balance between the production and the drainage of the kinetic energy of the turbulent large-scale structures is asymptotically stable as their amplitude saturates. The saturation of the wave amplitude provides an alternative indicator for flow self-similarity.

Class of self-limiting growth models in the presence of nonlinear diffusion

NASA Astrophysics Data System (ADS)

Kar, Sandip; Banik, Suman Kumar; Ray, Deb Shankar

2002-06-01

The source term in a reaction-diffusion system, in general, does not involve explicit time dependence. A class of self-limiting growth models dealing with animal and tumor growth and bacterial population in a culture, on the other hand, are described by kinetics with explicit functions of time. We analyze a reaction-diffusion system to study the propagation of spatial front for these models.

Breakdown parameter for kinetic modeling of multiscale gas flows.

PubMed

Meng, Jianping; Dongari, Nishanth; Reese, Jason M; Zhang, Yonghao

2014-06-01

Multiscale methods built purely on the kinetic theory of gases provide information about the molecular velocity distribution function. It is therefore both important and feasible to establish new breakdown parameters for assessing the appropriateness of a fluid description at the continuum level by utilizing kinetic information rather than macroscopic flow quantities alone. We propose a new kinetic criterion to indirectly assess the errors introduced by a continuum-level description of the gas flow. The analysis, which includes numerical demonstrations, focuses on the validity of the Navier-Stokes-Fourier equations and corresponding kinetic models and reveals that the new criterion can consistently indicate the validity of continuum-level modeling in both low-speed and high-speed flows at different Knudsen numbers.

A kinetic approach to magnetospheric modeling

NASA Technical Reports Server (NTRS)

Whipple, E. C., Jr.

1979-01-01

The earth's magnetosphere is caused by the interaction between the flowing solar wind and the earth's magnetic dipole, with the distorted magnetic field in the outer parts of the magnetosphere due to the current systems resulting from this interaction. It is surprising that even the conceptually simple problem of the collisionless interaction of a flowing plasma with a dipole magnetic field has not been solved. A kinetic approach is essential if one is to take into account the dispersion of particles with different energies and pitch angles and the fact that particles on different trajectories have different histories and may come from different sources. Solving the interaction problem involves finding the various types of possible trajectories, populating them with particles appropriately, and then treating the electric and magnetic fields self-consistently with the resulting particle densities and currents. This approach is illustrated by formulating a procedure for solving the collisionless interaction problem on open field lines in the case of a slowly flowing magnetized plasma interacting with a magnetic dipole.

Electrical conductivity of a methane-air burning plasma under the action of weak electric fields

NASA Astrophysics Data System (ADS)

Colonna, G.; Pietanza, L. D.; D'Angola, A.; Laricchiuta, A.; Di Vita, A.

2017-02-01

This paper focuses on the calculation of the electrical conductivity of a methane-air flame in the presence of weak electric fields, solving the Boltzmann equation for free electrons self-consistently coupled with chemical kinetics. The chemical model GRI-Mech 3.0 has been completed with chemi-ionization reactions to model ionization in the absence of fields, and a database of cross sections for electron-impact-induced processes to account for reactions and transitions activated in the flame during discharge. The dependence of plasma properties on the frequency of an oscillating field has been studied under different pressure and gas temperature conditions. Fitting expressions of the electrical conductivity as a function of gas temperature and methane consumption are provided for different operational conditions in the Ansaldo Energia burner.

Kinetic modeling of x-ray laser-driven solid Al plasmas via particle-in-cell simulation

NASA Astrophysics Data System (ADS)

Royle, R.; Sentoku, Y.; Mancini, R. C.; Paraschiv, I.; Johzaki, T.

2017-06-01

Solid-density plasmas driven by intense x-ray free-electron laser (XFEL) radiation are seeded by sources of nonthermal photoelectrons and Auger electrons that ionize and heat the target via collisions. Simulation codes that are commonly used to model such plasmas, such as collisional-radiative (CR) codes, typically assume a Maxwellian distribution and thus instantaneous thermalization of the source electrons. In this study, we present a detailed description and initial applications of a collisional particle-in-cell code, picls, that has been extended with a self-consistent radiation transport model and Monte Carlo models for photoionization and K L L Auger ionization, enabling the fully kinetic simulation of XFEL-driven plasmas. The code is used to simulate two experiments previously performed at the Linac Coherent Light Source investigating XFEL-driven solid-density Al plasmas. It is shown that picls-simulated pulse transmissions using the Ecker-Kröll continuum-lowering model agree much better with measurements than do simulations using the Stewart-Pyatt model. Good quantitative agreement is also found between the time-dependent picls results and those of analogous simulations by the CR code scfly, which was used in the analysis of the experiments to accurately reproduce the observed K α emissions and pulse transmissions. Finally, it is shown that the effects of the nonthermal electrons are negligible for the conditions of the particular experiments under investigation.

Complete protein-protein association kinetics in atomic detail revealed by molecular dynamics simulations and Markov modelling

NASA Astrophysics Data System (ADS)

Plattner, Nuria; Doerr, Stefan; de Fabritiis, Gianni; Noé, Frank

2017-10-01

Protein-protein association is fundamental to many life processes. However, a microscopic model describing the structures and kinetics during association and dissociation is lacking on account of the long lifetimes of associated states, which have prevented efficient sampling by direct molecular dynamics (MD) simulations. Here we demonstrate protein-protein association and dissociation in atomistic resolution for the ribonuclease barnase and its inhibitor barstar by combining adaptive high-throughput MD simulations and hidden Markov modelling. The model reveals experimentally consistent intermediate structures, energetics and kinetics on timescales from microseconds to hours. A variety of flexibly attached intermediates and misbound states funnel down to a transition state and a native basin consisting of the loosely bound near-native state and the tightly bound crystallographic state. These results offer a deeper level of insight into macromolecular recognition and our approach opens the door for understanding and manipulating a wide range of macromolecular association processes.

Nonequilibrium fixed points in longitudinally expanding scalar theories: Infrared cascade, Bose condensation and a challenge for kinetic theory

DOE PAGES

Berges, J.; Schlichting, S.; Boguslavski, K.; ...

2015-11-05

In [Phys. Rev. Lett. 114, 061601 (2015)], we reported on a new universality class for longitudinally expanding systems, encompassing strongly correlated non-Abelian plasmas and N-component self-interacting scalar field theories. Using classical-statistical methods, we showed that these systems share the same self-similar scaling properties for a wide range of momenta in a limit where particles are weakly coupled but their occupancy is high. Here we significantly expand on our previous work and delineate two further self-similar regimes. One of these occurs in the deep infrared (IR) regime of very high occupancies, where the nonequilibrium dynamics leads to the formation of amore » Bose-Einstein condensate. The universal IR scaling exponents and the spectral index characterizing the isotropic IR distributions are described by an effective theory derived from a systematic large-N expansion at next-to-leading order. Remarkably, this effective theory can be cast as a vertex-resummed kinetic theory. The other novel self-similar regime occurs close to the hard physical scale of the theory, and sets in only at later times. In this study, we argue that the important role of the infrared dynamics ensures that key features of our results for scalar and gauge theories cannot be reproduced consistently in conventional kinetic theory frameworks.« less

Nonequilibrium fixed points in longitudinally expanding scalar theories: Infrared cascade, Bose condensation and a challenge for kinetic theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Berges, J.; Schlichting, S.; Boguslavski, K.

In [Phys. Rev. Lett. 114, 061601 (2015)], we reported on a new universality class for longitudinally expanding systems, encompassing strongly correlated non-Abelian plasmas and N-component self-interacting scalar field theories. Using classical-statistical methods, we showed that these systems share the same self-similar scaling properties for a wide range of momenta in a limit where particles are weakly coupled but their occupancy is high. Here we significantly expand on our previous work and delineate two further self-similar regimes. One of these occurs in the deep infrared (IR) regime of very high occupancies, where the nonequilibrium dynamics leads to the formation of amore » Bose-Einstein condensate. The universal IR scaling exponents and the spectral index characterizing the isotropic IR distributions are described by an effective theory derived from a systematic large-N expansion at next-to-leading order. Remarkably, this effective theory can be cast as a vertex-resummed kinetic theory. The other novel self-similar regime occurs close to the hard physical scale of the theory, and sets in only at later times. In this study, we argue that the important role of the infrared dynamics ensures that key features of our results for scalar and gauge theories cannot be reproduced consistently in conventional kinetic theory frameworks.« less

Results from the OH-PT model: a Kinetic-MHD Model of the Outer Heliosphere within SWMF

NASA Astrophysics Data System (ADS)

Michael, A.; Opher, M.; Tenishev, V.; Borovikov, D.; Toth, G.

2017-12-01

We present an update of the OH-PT model, a kinetic-MHD model of the outer heliosphere. The OH-PT model couples the Outer Heliosphere (OH) and Particle Tracker (PT) components within the Space Weather Modeling Framework (SWMF). The OH component utilizes the Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) MHD code, a highly parallel, 3D, and block-adaptive solver. As a stand-alone model, the OH component solves the ideal MHD equations for the plasma and a separate set of Euler's equations for the different populations of neutral atoms. The neutrals and plasma in the outer heliosphere are coupled through charge-exchange. While this provides an accurate solution for the plasma, it is an inaccurate description of the neutrals. The charge-exchange mean free path is on the order of the size of the heliosphere; therefore the neutrals cannot be described as a fluid. The PT component is based on the Adaptive Mesh Particle Simulator (AMPS) model, a 3D, direct simulation Monte Carlo model that solves the Boltzmann equation for the motion and interaction of multi-species plasma and is used to model the neutral distribution functions throughout the domain. The charge-exchange process occurs within AMPS, which handles each event on a particle-by-particle basis and calculates the resulting source terms to the MHD equations. The OH-PT model combines the MHD solution for the plasma with the kinetic solution for the neutrals to form a self-consistent model of the heliosphere. In this work, we present verification and validation of the model as well as demonstrate the codes capabilities. Furthermore we provide a comparison of the OH-PT model to our multi-fluid approximation and detail the differences between the models in both the plasma solution and neutral distribution functions.

Self-consistent nonlocal feedback theory for electrocatalytic swimmers with heterogeneous surface chemical kinetics

NASA Astrophysics Data System (ADS)

Nourhani, Amir; Crespi, Vincent H.; Lammert, Paul E.

2015-06-01

We present a self-consistent nonlocal feedback theory for the phoretic propulsion mechanisms of electrocatalytic micromotors or nanomotors. These swimmers, such as bimetallic platinum and gold rods catalyzing decomposition of hydrogen peroxide in aqueous solution, have received considerable theoretical attention. In contrast, the heterogeneous electrochemical processes with nonlocal feedback that are the actual "engines" of such motors are relatively neglected. We present a flexible approach to these processes using bias potential as a control parameter field and a locally-open-circuit reference state, carried through in detail for a spherical motor. While the phenomenological flavor makes meaningful contact with experiment easier, required inputs can also conceivably come from, e.g., Frumkin-Butler-Volmer kinetics. Previously obtained results are recovered in the weak-heterogeneity limit and improved small-basis approximations tailored to structural heterogeneity are presented. Under the assumption of weak inhomogeneity, a scaling form is deduced for motor speed as a function of fuel concentration and swimmer size. We argue that this form should be robust and demonstrate a good fit to experimental data.

A Rate-Theory-Phase-Field Model of Irradiation-Induced Recrystallization in UMo Nuclear Fuels

NASA Astrophysics Data System (ADS)

Hu, Shenyang; Joshi, Vineet; Lavender, Curt A.

2017-12-01

In this work, we developed a recrystallization model to study the effect of microstructures and radiation conditions on recrystallization kinetics in UMo fuels. The model integrates the rate theory of intragranular gas bubble and interstitial loop evolutions and a phase-field model of recrystallization zone evolution. A first passage method is employed to describe one-dimensional diffusion of interstitials with a diffusivity value several orders of magnitude larger than that of fission gas xenons. With the model, the effect of grain sizes on recrystallization kinetics is simulated. The results show that (1) recrystallization in large grains starts earlier than that in small grains, (2) the recrystallization kinetics (recrystallization volume fraction) decrease as the grain size increases, (3) the predicted recrystallization kinetics are consistent with the experimental results, and (4) the recrystallization kinetics can be described by the modified Avrami equation, but the parameters of the Avrami equation strongly depend on the grain size.

Modelling DNA origami self-assembly at the domain level.

PubMed

Dannenberg, Frits; Dunn, Katherine E; Bath, Jonathan; Kwiatkowska, Marta; Turberfield, Andrew J; Ouldridge, Thomas E

2015-10-28

We present a modelling framework, and basic model parameterization, for the study of DNA origami folding at the level of DNA domains. Our approach is explicitly kinetic and does not assume a specific folding pathway. The binding of each staple is associated with a free-energy change that depends on staple sequence, the possibility of coaxial stacking with neighbouring domains, and the entropic cost of constraining the scaffold by inserting staple crossovers. A rigorous thermodynamic model is difficult to implement as a result of the complex, multiply connected geometry of the scaffold: we present a solution to this problem for planar origami. Coaxial stacking of helices and entropic terms, particularly when loop closure exponents are taken to be larger than those for ideal chains, introduce interactions between staples. These cooperative interactions lead to the prediction of sharp assembly transitions with notable hysteresis that are consistent with experimental observations. We show that the model reproduces the experimentally observed consequences of reducing staple concentration, accelerated cooling, and absent staples. We also present a simpler methodology that gives consistent results and can be used to study a wider range of systems including non-planar origami.

Modelling DNA origami self-assembly at the domain level

NASA Astrophysics Data System (ADS)

Dannenberg, Frits; Dunn, Katherine E.; Bath, Jonathan; Kwiatkowska, Marta; Turberfield, Andrew J.; Ouldridge, Thomas E.

2015-10-01

We present a modelling framework, and basic model parameterization, for the study of DNA origami folding at the level of DNA domains. Our approach is explicitly kinetic and does not assume a specific folding pathway. The binding of each staple is associated with a free-energy change that depends on staple sequence, the possibility of coaxial stacking with neighbouring domains, and the entropic cost of constraining the scaffold by inserting staple crossovers. A rigorous thermodynamic model is difficult to implement as a result of the complex, multiply connected geometry of the scaffold: we present a solution to this problem for planar origami. Coaxial stacking of helices and entropic terms, particularly when loop closure exponents are taken to be larger than those for ideal chains, introduce interactions between staples. These cooperative interactions lead to the prediction of sharp assembly transitions with notable hysteresis that are consistent with experimental observations. We show that the model reproduces the experimentally observed consequences of reducing staple concentration, accelerated cooling, and absent staples. We also present a simpler methodology that gives consistent results and can be used to study a wider range of systems including non-planar origami.

Consistent Chemical Mechanism from Collaborative Data Processing

DOE PAGES

Slavinskaya, Nadezda; Starcke, Jan-Hendrik; Abbasi, Mehdi; ...

2016-04-01

Numerical tool of Process Informatics Model (PrIMe) is mathematically rigorous and numerically efficient approach for analysis and optimization of chemical systems. It handles heterogeneous data and is scalable to a large number of parameters. The Boundto-Bound Data Collaboration module of the automated data-centric infrastructure of PrIMe was used for the systematic uncertainty and data consistency analyses of the H 2/CO reaction model (73/17) and 94 experimental targets (ignition delay times). The empirical rule for evaluation of the shock tube experimental data is proposed. The initial results demonstrate clear benefits of the PrIMe methods for an evaluation of the kinetic datamore » quality and data consistency and for developing predictive kinetic models.« less

Kinetic and radiation-hydrodynamic modeling of x-ray heating in laboratory photoionized plasmas

NASA Astrophysics Data System (ADS)

Mancini, Roberto

2017-06-01

In experiments performed at the Z facility of Sandia National Laboratories a cm-scale cell filled with neon gas was driven by the burst of broadband x-rays emitted at the collapse of a wire-array z-pinch turning the gas into a photoionized plasma. Transmission spectroscopy of a narrowband portion of the x-ray flux was used to diagnose the plasma. The data show a highly-ionized neon plasma with a rich line absorption spectrum that permits the extraction of the ionization distribution among Be-, Li-, He- and H-like ions. Analysis of the spectra produced atomic ground and low excited state areal densities in these ions, and from the ratio of first-excited to ground state populations in Li-like neon a temperature of 19±4eV was extracted to characterize the x-ray heating of the plasma. To interpret this observation, we have performed data-constrained view-factor calculations of the spectral distribution of the x-ray drive, self-consistent modeling of electron and atomic kinetics, and radiation-hydrodynamic simulations. For the conditions of the experiment, the electron distribution thermalizes quickly, has a negligible high-energy tail, and is very well approximated by a single Maxwellian distribution. Radiation-hydrodynamic simulations with either LTE or NLTE (i.e. non-equilibrium) atomic physics provide a more complete modeling of the experiment. We found that in order to compute electron temperatures consistent with observation inline non-equilibrium collisional-radiative neon atomic kinetics needs to be taken into account. We discuss the details of LTE and NLTE simulations, and the impact of atomic physics on the radiation heating and cooling rates that determine the plasma temperature. This work was sponsored in part by DOE Office of Science Grant DE-SC0014451, and the Z Facility Fundamental Science Program of SNL.

Theory of low transitions in CO discharge lasers

NASA Technical Reports Server (NTRS)

Sidney, B. D.; Mcinuille, R. M.; Smith, N. S.; Hassan, H. A.

1976-01-01

A self consistent theoretical model which couples the electron and heavy particle finite rate kinetics with the optical and fluid dynamic processes has been employed to identify the various parameters and explain the mechanism responsible for producing low lying transitions in slow flowing CO lasers. It is found that lasing on low lying transitions can be achieved at low temperatures for low pressures (or low flow rates) together with high partial pressures of the He and N2. The role of N2 has been identified as an additive responsible for reducing the electron temperature to a range where the transfer of electrical power to the lower vibrational modes of CO is optimum.

Kinetics of Substrate Biodegradation under the Cumulative Effects of Bioavailability and Self-Inhibition.

PubMed

Gharasoo, Mehdi; Centler, Florian; Van Cappellen, Philippe; Wick, Lukas Y; Thullner, Martin

2015-05-05

Microbial degradation is an important process in many environments controlling for instance the cycling of nutrients or the biodegradation of contaminants. At high substrate concentrations toxic effects may inhibit the degradation process. Bioavailability limitations of a degradable substrate can therefore either improve the overall dynamics of degradation by softening the contaminant toxicity effects to microorganisms, or slow down the biodegradation by reducing the microbial access to the substrate. Many studies on biodegradation kinetics of a self-inhibitive substrate have mainly focused on physiological responses of the bacteria to substrate concentration levels without considering the substrate bioavailability limitations rising from different geophysical and geochemical dynamics at pore-scale. In this regard, the role of bioavailability effects on the kinetics of self-inhibiting substrates is poorly understood. In this study, we theoretically analyze this role and assess the interactions between self-inhibition and mass transfer-limitations using analytical/numerical solutions, and show the findings practical relevance for a simple model scenario. Although individually self-inhibition and mass-transfer limitations negatively impact biodegradation, their combined effect may enhance biodegradation rates above a concentration threshold. To our knowledge, this is the first theoretical study describing the cumulative effects of the two mechanisms together.

Quasilinear diffusion coefficients in a finite Larmor radius expansion for ion cyclotron heated plasmas

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lee, Jungpyo; Wright, John; Bertelli, Nicola

In this study, a reduced model of quasilinear velocity diffusion by a small Larmor radius approximation is derived to couple the Maxwell’s equations and the Fokker Planck equation self-consistently for the ion cyclotron range of frequency waves in a tokamak. The reduced model ensures the important properties of the full model by Kennel-Engelmann diffusion, such as diffusion directions, wave polarizations, and H-theorem. The kinetic energy change (Wdot ) is used to derive the reduced model diffusion coefficients for the fundamental damping (n = 1) and the second harmonic damping (n = 2) to the lowest order of the finite Larmormore » radius expansion. The quasilinear diffusion coefficients are implemented in a coupled code (TORIC-CQL3D) with the equivalent reduced model of the dielectric tensor. We also present the simulations of the ITER minority heating scenario, in which the reduced model is verified within the allowable errors from the full model results.« less

Quasilinear diffusion coefficients in a finite Larmor radius expansion for ion cyclotron heated plasmas

DOE PAGES

Lee, Jungpyo; Wright, John; Bertelli, Nicola; ...

2017-04-24

In this study, a reduced model of quasilinear velocity diffusion by a small Larmor radius approximation is derived to couple the Maxwell’s equations and the Fokker Planck equation self-consistently for the ion cyclotron range of frequency waves in a tokamak. The reduced model ensures the important properties of the full model by Kennel-Engelmann diffusion, such as diffusion directions, wave polarizations, and H-theorem. The kinetic energy change (Wdot ) is used to derive the reduced model diffusion coefficients for the fundamental damping (n = 1) and the second harmonic damping (n = 2) to the lowest order of the finite Larmormore » radius expansion. The quasilinear diffusion coefficients are implemented in a coupled code (TORIC-CQL3D) with the equivalent reduced model of the dielectric tensor. We also present the simulations of the ITER minority heating scenario, in which the reduced model is verified within the allowable errors from the full model results.« less

Self-consistent current sheet structures in the quiet-time magnetotail

NASA Technical Reports Server (NTRS)

Holland, Daniel L.; Chen, James

1993-01-01

The structure of the quiet-time magnetotail is studied using a test particle simulation. Vlasov equilibria are obtained in the regime where v(D) = E(y) c/B(z) is much less than the ion thermal velocity and are self-consistent in that the current and magnetic field satisfy Ampere's law. Force balance between the plasma and magnetic field is satisfied everywhere. The global structure of the current sheet is found to be critically dependent on the source distribution function. The pressure tensor is nondiagonal in the current sheet with anisotropic temperature. A kinetic mechanism is proposed whereby changes in the source distribution results in a thinning of the current sheet.

Kinetic Modeling of Ultraintense X-Ray Laser-Matter Interactions

NASA Astrophysics Data System (ADS)

Royle, Ryan; Sentoku, Yasuhiko; Mancini, Roberto; Johzaki, Tomoyuki

2015-11-01

High-intensity XFELs have become a novel way of creating and studying hot dense plasmas. The LCLS at Stanford can deliver a millijoule of energy with more than 1012 photons in a ~ 100 femtosecond pulse. By tightly focusing the beam to a micron-scale spot size, the XFEL can be intensified to more than 1018 W/cm2, making it possible to heat solid matter isochorically beyond a million degrees (>100 eV). Such extreme states of matter are of considerable interest due to their relevance to astrophysical plasmas. Additionally, they will allow novel ways of studying equation-of-state and opacity physics under Gbar pressure and strong fields. Photoionization is the dominant x-ray absorption mechanism and triggers the heating processes. A photoionization model that takes into account the subshell cross-sections has been developed in a kinetic plasma simulation code, PICLS, that solves the x-ray transport self-consistently. The XFEL-matter interaction with several elements, including solid carbon, aluminum, and iron, is studied with the code, and the results are compared with recent LCLS experiments. This work was supported by the DOE/OFES under Contract No. DE-SC0008827.

Kinetic and morphology study of alginate-vineyard pruning waste biocomposite vs. non modified vineyard pruning waste for dye removal.

PubMed

Vecino, Xanel; Devesa-Rey, Rosa; Villagrasa, Salvador; Cruz, Jose M; Moldes, Ana B

2015-12-01

In this work a comparative bioadsorption study between a biocomposite consisting of hydrolysed vineyard pruning waste entrapped in calcium alginate spheres and non entrapped vineyard residue was carried out. Results have demonstrated that the biocomposite based on lignocellulose-calcium alginate spheres removed 77.3% of dyes, while non entrapped lignocellulose eliminated only removed 27.8% of colour compounds. The experimental data were fitted to several kinetic models (pseudo-first order, pseudo-second order, Chien-Clayton model, intraparticle diffusion model and Bangham model); being pseudo-second order the kinetic model that better described the adsorption of dyes onto both bioadsorbents. In addition, a morphological study (roughness and shape) of alginate-vineyard biocomposite was established under extreme conditions, observing significant differences between hydrated and dehydrated alginate-vineyard biocomposite. The techniques used to carry out this morphological study consisted of scanning electron microscopy (SEM), perfilometry and 3D surface analysis. Copyright © 2015. Published by Elsevier B.V.

A double-layer based model of ion confinement in electron cyclotron resonance ion source.

PubMed

Mascali, D; Neri, L; Celona, L; Castro, G; Torrisi, G; Gammino, S; Sorbello, G; Ciavola, G

2014-02-01

The paper proposes a new model of ion confinement in ECRIS, which can be easily generalized to any magnetic configuration characterized by closed magnetic surfaces. Traditionally, ion confinement in B-min configurations is ascribed to a negative potential dip due to superhot electrons, adiabatically confined by the magneto-static field. However, kinetic simulations including RF heating affected by cavity modes structures indicate that high energy electrons populate just a thin slab overlapping the ECR layer, while their density drops down of more than one order of magnitude outside. Ions, instead, diffuse across the electron layer due to their high collisionality. This is the proper physical condition to establish a double-layer (DL) configuration which self-consistently originates a potential barrier; this "barrier" confines the ions inside the plasma core surrounded by the ECR surface. The paper will describe a simplified ion confinement model based on plasma density non-homogeneity and DL formation.

One-dimensional model of interacting-step fluctuations on vicinal surfaces: Analytical formulas and kinetic Monte Carlo simulations

NASA Astrophysics Data System (ADS)

Patrone, Paul N.; Einstein, T. L.; Margetis, Dionisios

2010-12-01

We study analytically and numerically a one-dimensional model of interacting line defects (steps) fluctuating on a vicinal crystal. Our goal is to formulate and validate analytical techniques for approximately solving systems of coupled nonlinear stochastic differential equations (SDEs) governing fluctuations in surface motion. In our analytical approach, the starting point is the Burton-Cabrera-Frank (BCF) model by which step motion is driven by diffusion of adsorbed atoms on terraces and atom attachment-detachment at steps. The step energy accounts for entropic and nearest-neighbor elastic-dipole interactions. By including Gaussian white noise to the equations of motion for terrace widths, we formulate large systems of SDEs under different choices of diffusion coefficients for the noise. We simplify this description via (i) perturbation theory and linearization of the step interactions and, alternatively, (ii) a mean-field (MF) approximation whereby widths of adjacent terraces are replaced by a self-consistent field but nonlinearities in step interactions are retained. We derive simplified formulas for the time-dependent terrace-width distribution (TWD) and its steady-state limit. Our MF analytical predictions for the TWD compare favorably with kinetic Monte Carlo simulations under the addition of a suitably conservative white noise in the BCF equations.

Forced unfolding of single-chain polymeric nanoparticles.

PubMed

Hosono, Nobuhiko; Kushner, Aaron M; Chung, Jaeyoon; Palmans, Anja R A; Guan, Zhibin; Meijer, E W

2015-06-03

Atomic force microscopy (AFM)-based single-molecule force spectroscopy (SMFS) is applied to single-chain polymeric nanoparticles (SCPNs) to acquire information about the internal folding structure of SCPNs and inherent kinetic parameters of supramolecular self-assembling motifs embedded into the SCPNs. The SCPNs used here are polyacrylate-based polymers carrying 2-ureido-4-[1H]-pyrimidinone (UPy) or benzene-1,3,5-tricarboxamide (BTA) pendants that induce an intramolecular chain collapse into nanoparticles consisting of one polymer chain only via internal supramolecular cross-linking. The SCPN is stretched by an AFM cantilever to unfold mechanically, which allows measuring of force-extension profiles of the SCPNs. Consecutive peaks observed in the force profiles are attributed to rupture events of self-assembled UPy/BTA units in the SCPNs. The force profiles have been analyzed statistically for a series of polymers with different UPy/BTA incorporation densities. The results provide insights into the internal conformation of SCPNs, where the folding structure can be changed with the incorporation density of UPy/BTA. In addition, dynamic loading rate analysis allows the determination of kinetic parameters of BTA self-assembly, which has not been accessible by any other method. This study offers a rational tool for understanding the folding structure, kinetics, and pathway of two series of SCPNs.

Analysis of a kinetic multi-segment foot model part II: kinetics and clinical implications.

PubMed

Bruening, Dustin A; Cooney, Kevin M; Buczek, Frank L

2012-04-01

Kinematic multi-segment foot models have seen increased use in clinical and research settings, but the addition of kinetics has been limited and hampered by measurement limitations and modeling assumptions. In this second of two companion papers, we complete the presentation and analysis of a three segment kinetic foot model by incorporating kinetic parameters and calculating joint moments and powers. The model was tested on 17 pediatric subjects (ages 7-18 years) during normal gait. Ground reaction forces were measured using two adjacent force platforms, requiring targeted walking and the creation of two sub-models to analyze ankle, midtarsal, and 1st metatarsophalangeal joints. Targeted walking resulted in only minimal kinematic and kinetic differences compared with walking at self selected speeds. Joint moments and powers were calculated and ensemble averages are presented as a normative database for comparison purposes. Ankle joint powers are shown to be overestimated when using a traditional single-segment foot model, as substantial angular velocities are attributed to the mid-tarsal joint. Power transfer is apparent between the 1st metatarsophalangeal and mid-tarsal joints in terminal stance/pre-swing. While the measurement approach presented here is limited to clinical populations with only minimal impairments, some elements of the model can also be incorporated into routine clinical gait analysis. Copyright © 2011 Elsevier B.V. All rights reserved.

Towards an evolutionary theory of the origin of life based on kinetics and thermodynamics.

PubMed

Pascal, Robert; Pross, Addy; Sutherland, John D

2013-11-06

A sudden transition in a system from an inanimate state to the living state-defined on the basis of present day living organisms-would constitute a highly unlikely event hardly predictable from physical laws. From this uncontroversial idea, a self-consistent representation of the origin of life process is built up, which is based on the possibility of a series of intermediate stages. This approach requires a particular kind of stability for these stages-dynamic kinetic stability (DKS)-which is not usually observed in regular chemistry, and which is reflected in the persistence of entities capable of self-reproduction. The necessary connection of this kinetic behaviour with far-from-equilibrium thermodynamic conditions is emphasized and this leads to an evolutionary view for the origin of life in which multiplying entities must be associated with the dissipation of free energy. Any kind of entity involved in this process has to pay the energetic cost of irreversibility, but, by doing so, the contingent emergence of new functions is made feasible. The consequences of these views on the studies of processes by which life can emerge are inferred.

A General Framework for Thermodynamically Consistent Parameterization and Efficient Sampling of Enzymatic Reactions

PubMed Central

Saa, Pedro; Nielsen, Lars K.

2015-01-01

Kinetic models provide the means to understand and predict the dynamic behaviour of enzymes upon different perturbations. Despite their obvious advantages, classical parameterizations require large amounts of data to fit their parameters. Particularly, enzymes displaying complex reaction and regulatory (allosteric) mechanisms require a great number of parameters and are therefore often represented by approximate formulae, thereby facilitating the fitting but ignoring many real kinetic behaviours. Here, we show that full exploration of the plausible kinetic space for any enzyme can be achieved using sampling strategies provided a thermodynamically feasible parameterization is used. To this end, we developed a General Reaction Assembly and Sampling Platform (GRASP) capable of consistently parameterizing and sampling accurate kinetic models using minimal reference data. The former integrates the generalized MWC model and the elementary reaction formalism. By formulating the appropriate thermodynamic constraints, our framework enables parameterization of any oligomeric enzyme kinetics without sacrificing complexity or using simplifying assumptions. This thermodynamically safe parameterization relies on the definition of a reference state upon which feasible parameter sets can be efficiently sampled. Uniform sampling of the kinetics space enabled dissecting enzyme catalysis and revealing the impact of thermodynamics on reaction kinetics. Our analysis distinguished three reaction elasticity regions for common biochemical reactions: a steep linear region (0> ΔGr >-2 kJ/mol), a transition region (-2> ΔGr >-20 kJ/mol) and a constant elasticity region (ΔGr <-20 kJ/mol). We also applied this framework to model more complex kinetic behaviours such as the monomeric cooperativity of the mammalian glucokinase and the ultrasensitive response of the phosphoenolpyruvate carboxylase of Escherichia coli. In both cases, our approach described appropriately not only the kinetic behaviour of these enzymes, but it also provided insights about the particular features underpinning the observed kinetics. Overall, this framework will enable systematic parameterization and sampling of enzymatic reactions. PMID:25874556

The kinetics and dynamics of the coma of Halley's comet

NASA Technical Reports Server (NTRS)

Combi, Michael R.

1994-01-01

This grant to the University of Michigan supported the efforts of Michael R. Combi to serve as a co-investigator in collaboration with a larger effort by the principal investigator, William Smyth of Atmospheric and Environmental Research, Inc. The overall objective of this project was to analyze in a self-consistent manner unique optical O((sup 1)D) and NH2 ultra-high resolution line profile data of excellent quality and other supporting lower-resolution spectral data for the coma of comet P/Halley by using highly developed and physically-based cometary coma models in order to determine and explain in terms of physical processes the actual dynamics and photochemical kinetics that occur in the coma. The justification for this work is that it provides a valuable and underlying physical base from which to interpret significantly different types of coma observations in a self-consistent manner and hence bring into agreement (or avoid) apparent inconsistencies that arise from non-physically based interpretations. The level of effort for the Michigan component amounted to less than three person-months over a planned period of three years. The period had been extended at no extra cost to four years because the Michigan grant and the AER contract did not have coincident time periods. An effort of somewhat larger scope was undertaken by the PI. The importance of the O((sup 1)D) profiles is that they provide a direct trace of the water distribution in comets. The line profile shape is produced by the convolution of the outflow velocity and thermal dispersion of the parent water molecules with the photokinetic ejection of the oxygen atoms upon photodissociation of the parent water molecules. Our understanding of the NH2 and its precursor ammonia are important for comet-to-comet composition variations as they relate to the cosmo-chemistry of the early solar nebula. Modeling of the distribution of NH2 is necessary in order to infer the ammonia production rates from NH2 observations.

A Reduced-order NLTE Kinetic Model for Radiating Plasmas of Outer Envelopes of Stellar Atmospheres

DOE Office of Scientific and Technical Information (OSTI.GOV)

Munafò, Alessandro; Mansour, Nagi N.; Panesi, Marco, E-mail: munafo@illinois.edu, E-mail: nagi.n.mansour@nasa.gov, E-mail: m.panesi@illinois.edu

The present work proposes a self-consistent reduced-order NLTE kinetic model for radiating plasmas found in the outer layers of stellar atmospheres. A detailed collisional-radiative kinetic mechanism is constructed by leveraging the most up-to-date set of ab initio and experimental data available in the literature. This constitutes the starting point for the derivation of a reduced-order model, obtained by lumping the bound energy states into groups. In order to determine the needed thermo-physical group properties, uniform and Maxwell–Boltzmann energy distributions are used to reconstruct the energy population of each group. Finally, the reduced set of governing equations for the material gasmore » and the radiation field is obtained based on the moment method. Applications consider the steady flow across a shock wave in partially ionized hydrogen. The results clearly demonstrate that adopting a Maxwell–Boltzmann grouping allows, on the one hand, for a substantial reduction of the number of unknowns and, on the other, to maintain accuracy for both gas and radiation quantities. Also, it is observed that, when neglecting line radiation, the use of two groups already leads to a very accurate resolution of the photo-ionization precursor, internal relaxation, and radiative cooling regions. The inclusion of line radiation requires adopting just one additional group to account for optically thin losses in the α , β , and γ lines of the Balmer and Paschen series. This trend has been observed for a wide range of shock wave velocities.« less

Polymer damage mitigation---predictive lifetime models of polymer insulation degradation and biorenewable thermosets through cationic polymerization for self-healing applications

NASA Astrophysics Data System (ADS)

Hondred, Peter Raymond

Over the past 50 years, the industrial development and applications for polymers and polymer composites has become expansive. However, as with any young technology, the techniques for predicting material damage and resolving material failure are in need of continued development and refinement. This thesis work takes two approaches to polymer damage mitigation---material lifetime prediction and spontaneous damage repair through self-healing while incorporating bio-renewable feedstock. First, material lifetime prediction offers the benefit of identifying and isolating material failures before the effects of damage results in catastrophic failure. Second, self-healing provides a systematic approach to repairing damaged polymer composites, specifically in applications where a hands-on approach or removing the part from service are not feasible. With regard to lifetime prediction, we investigated three specific polymeric materials---polytetrafluoroethylene (PTFE), poly(ethylene-alt-tetrafluoroethylene) (ETFE), and Kapton. All three have been utilized extensively in the aerospace field as a wire insulation coating. Because of the vast amount of electrical wiring used in aerospace constructions and the potential for electrical and thermal failure, this work develops mathematical models for both the thermal degradation kinetics as well as a lifetime prediction model for electrothermal breakdown. Isoconversional kinetic methods, which plot activation energy as a function of the extent of degradation, present insight into the development each kinetic model. The models for PTFE, ETFE, and Kapton are one step, consecutive three-step, and competitive and consecutive five-step respectively. Statistical analysis shows that an nth order autocatalytic reaction best defined the reaction kinetics for each polymer's degradation. Self-healing polymers arrest crack propagation through the use of an imbedded adhesive that reacts when cracks form. This form of damage mitigation focuses on repairing damage before the damage causes a failure in the polymer's function. In this work, the healing agent (adhesive) is developed using bio-renewable oils instead of solely relying on petroleum based feedstocks. Several bio-renewable thermosetting polymers were successfully prepared from tung oil through cationic polymerization for the use as the healing agent in self-healing microencapsulated applications. Modifications to both the monomers in the resin and the catalyst for polymerization were made and the subsequent changes to mechanical, thermal, and structural properties were identified. Furthermore, compressive lap shear testing was used to confirm that the adhesive properties would be beneficial for self-healing applications. Finally, scanning electron microscopy of the crack plane was used to study the fracture mechanism of the crack.

Dynamic reassembly of peptide RADA16 nanofiber scaffold

NASA Astrophysics Data System (ADS)

Yokoi, Hidenori; Kinoshita, Takatoshi; Zhang, Shuguang

2005-06-01

Nanofiber structures of some peptides and proteins as biological materials have been studied extensively, but their molecular mechanism of self-assembly and reassembly still remains unclear. We report here the reassembly of an ionic self-complementary peptide RADARADARADARADA (RADA16-I) that forms a well defined nanofiber scaffold. The 16-residue peptide forms stable -sheet structure and undergoes molecular self-assembly into nanofibers and eventually a scaffold hydrogel consisting of >99.5% water. In this study, the nanofiber scaffold was sonicated into smaller fragments. Circular dichroism, atomic force microscopy, and rheology were used to follow the kinetics of the reassembly. These sonicated fragments not only quickly reassemble into nanofibers that were indistinguishable from the original material, but their reassembly also correlated with the rheological analyses showing an increase of scaffold rigidity as a function of nanofiber length. The disassembly and reassembly processes were repeated four times and, each time, the reassembly reached the original length. We proposed a plausible sliding diffusion model to interpret the reassembly involving complementary nanofiber cohesive ends. This reassembly process is important for fabrication of new scaffolds for 3D cell culture, tissue repair, and regenerative medicine. atomic force microscopy | circular dichroism | dynamic behaviors | ionic self-complementary peptides | nanofiber hydrogels

Hydrodynamics beyond Navier-Stokes: the slip flow model.

PubMed

Yudistiawan, Wahyu P; Ansumali, Santosh; Karlin, Iliya V

2008-07-01

Recently, analytical solutions for the nonlinear Couette flow demonstrated the relevance of the lattice Boltzmann (LB) models to hydrodynamics beyond the continuum limit [S. Ansumali, Phys. Rev. Lett. 98, 124502 (2007)]. In this paper, we present a systematic study of the simplest LB kinetic equation-the nine-bit model in two dimensions--in order to quantify it as a slip flow approximation. Details of the aforementioned analytical solution are presented, and results are extended to include a general shear- and force-driven unidirectional flow in confined geometry. Exact solutions for the velocity, as well as for pertinent higher-order moments of the distribution functions, are obtained in both Couette and Poiseuille steady-state flows for all values of rarefaction parameter (Knudsen number). Results are compared with the slip flow solution by Cercignani, and a good quantitative agreement is found for both flow situations. Thus, the standard nine-bit LB model is characterized as a valid and self-consistent slip flow model for simulations beyond the Navier-Stokes approximation.

Flow-Induced New Channels of Energy Exchange in Multi-Scale Plasma Dynamics - Revisiting Perturbative Hybrid Kinetic-MHD Theory.

PubMed

Shiraishi, Junya; Miyato, Naoaki; Matsunaga, Go

2016-05-10

It is found that new channels of energy exchange between macro- and microscopic dynamics exist in plasmas. They are induced by macroscopic plasma flow. This finding is based on the kinetic-magnetohydrodynamic (MHD) theory, which analyses interaction between macroscopic (MHD-scale) motion and microscopic (particle-scale) dynamics. The kinetic-MHD theory is extended to include effects of macroscopic plasma flow self-consistently. The extension is realised by generalising an energy exchange term due to wave-particle resonance, denoted by δ WK. The first extension is generalisation of the particle's Lagrangian, and the second one stems from modification to the particle distribution function due to flow. These extensions lead to a generalised expression of δ WK, which affects the MHD stability of plasmas.

Supermassive black holes and their feedback effects in the IllustrisTNG simulation

NASA Astrophysics Data System (ADS)

Weinberger, Rainer; Springel, Volker; Pakmor, Rüdiger; Nelson, Dylan; Genel, Shy; Pillepich, Annalisa; Vogelsberger, Mark; Marinacci, Federico; Naiman, Jill; Torrey, Paul; Hernquist, Lars

2018-06-01

We study the population of supermassive black holes (SMBHs) and their effects on massive central galaxies in the IllustrisTNG cosmological hydrodynamical simulations of galaxy formation. The employed model for SMBH growth and feedback assumes a two-mode scenario in which the feedback from active galactic nuclei occurs through a kinetic, comparatively efficient mode at low accretion rates relative to the Eddington limit, and in the form of a thermal, less efficient mode at high accretion rates. We show that the quenching of massive central galaxies happens coincidently with kinetic-mode feedback, consistent with the notion that active supermassive black cause the low specific star formation rates observed in massive galaxies. However, major galaxy mergers are not responsible for initiating most of the quenching events in our model. Up to black hole masses of about 108.5 M⊙, the dominant growth channel for SMBHs is in the thermal mode. Higher mass black holes stay mainly in the kinetic mode and gas accretion is self-regulated via their feedback, which causes their Eddington ratios to drop, with SMBH mergers becoming the main channel for residual mass growth. As a consequence, the quasar luminosity function is dominated by rapidly accreting, moderately massive black holes in the thermal mode. We show that the associated growth history of SMBHs produces a low-redshift quasar luminosity function and a redshift zero black hole mass - stellar bulge mass relation in good agreement with observations, whereas the simulation tends to over-predict the high-redshift quasar luminosity function.

On the interpretation of energy and energy fluxes of nonlinear internal waves: An example from Massachusetts Bay

USGS Publications Warehouse

Scotti, A.; Beardsley, R.; Butman, B.

2006-01-01

A self-consistent formalism to estimate baroclinic energy densities and fluxes resulting from the propagation of internal waves of arbitrary amplitude is derived using the concept of available potential energy. The method can be applied to numerical, laboratory or field data. The total energy flux is shown to be the sum of the linear energy flux ??? u??? p??? dz (primes denote baroclinic quantities), plus contributions from the non-hydrostatic pressure anomaly and the self-advection of kinetic and available potential energy. Using highly resolved observations in Massachusetts Bay, it is shown that due to the presence of nonlinear internal waves periodically propagating in the area, ??? u??? p??? dz accounts for only half of the total flux. The same data show that equipartition of available potential and kinetic energy can be violated, especially when the nonlinear waves begin to interact with the bottom. ?? 2006 Cambridge University Press.

A comparison of observed and numerically predicted eddy kinetic energy budgets for a developing extratropical cyclone

NASA Technical Reports Server (NTRS)

Dare, P. M.; Smith, P. J.

1983-01-01

The eddy kinetic energy budget is calculated for a 48-hour forecast of an intense occluding winter cyclone associated with a strong well-developed jet stream. The model output consists of the initialized (1200 GMT January 9, 1975) and the 12, 24, 36, and 48 hour forecast fields from the Drexel/NCAR Limited Area Mesoscale Prediction System (LAMPS) model. The LAMPS forecast compares well with observations for the first 24 hours, but then overdevelops the low-level cyclone while inadequately developing the upper-air wave and jet. Eddy kinetic energy was found to be concentrated in the upper-troposphere with maxima flanking the primary trough. The increases in kinetic energy were found to be due to an excess of the primary source term of kinetic energy content, which is the horizontal flux of eddy kinetic energy over the primary sinks, and the generation and dissipation of eddy kinetic energy.

Thermodynamic Analysis of Chemically Reacting Mixtures-Comparison of First and Second Order Models.

PubMed

Pekař, Miloslav

2018-01-01

Recently, a method based on non-equilibrium continuum thermodynamics which derives thermodynamically consistent reaction rate models together with thermodynamic constraints on their parameters was analyzed using a triangular reaction scheme. The scheme was kinetically of the first order. Here, the analysis is further developed for several first and second order schemes to gain a deeper insight into the thermodynamic consistency of rate equations and relationships between chemical thermodynamic and kinetics. It is shown that the thermodynamic constraints on the so-called proper rate coefficient are usually simple sign restrictions consistent with the supposed reaction directions. Constraints on the so-called coupling rate coefficients are more complex and weaker. This means more freedom in kinetic coupling between reaction steps in a scheme, i.e., in the kinetic effects of other reactions on the rate of some reaction in a reacting system. When compared with traditional mass-action rate equations, the method allows a reduction in the number of traditional rate constants to be evaluated from data, i.e., a reduction in the dimensionality of the parameter estimation problem. This is due to identifying relationships between mass-action rate constants (relationships which also include thermodynamic equilibrium constants) which have so far been unknown.

The interaction between fishbone modes and shear Alfvén waves in tokamak plasmas

NASA Astrophysics Data System (ADS)

He, Hongda; Liu, Yueqiang; Dong, J. Q.; Hao, G. Z.; Wu, Tingting; He, Zhixiong; Zhao, K.

2016-05-01

The resonant interaction between the energetic particle triggered fishbone mode and the shear Alfvén waves is computationally investigated and firmly demonstrated based on a tokamak plasma equilibrium, using the self-consistent MHD-kinetic hybrid code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503). This type of continuum resonance, occurring critically due to the mode’s toroidal rotation in the plasma frame, significantly modifies the eigenmode structure of the fishbone instability, by introducing two large peaks of the perturbed parallel current density near but offside the q  =  1 rational surface (q is the safety factor). The self-consistently computed radial plasma displacement substantially differs from that being assumed in the conventional fishbone theory.

Synoptic, Global Mhd Model For The Solar Corona

NASA Astrophysics Data System (ADS)

Cohen, Ofer; Sokolov, I. V.; Roussev, I. I.; Gombosi, T. I.

2007-05-01

The common techniques for mimic the solar corona heating and the solar wind acceleration in global MHD models are as follow. 1) Additional terms in the momentum and energy equations derived from the WKB approximation for the Alfv’en wave turbulence; 2) some empirical heat source in the energy equation; 3) a non-uniform distribution of the polytropic index, γ, used in the energy equation. In our model, we choose the latter approach. However, in order to get a more realistic distribution of γ, we use the empirical Wang-Sheeley-Arge (WSA) model to constrain the MHD solution. The WSA model provides the distribution of the asymptotic solar wind speed from the potential field approximation; therefore it also provides the distribution of the kinetic energy. Assuming that far from the Sun the total energy is dominated by the energy of the bulk motion and assuming the conservation of the Bernoulli integral, we can trace the total energy along a magnetic field line to the solar surface. On the surface the gravity is known and the kinetic energy is negligible. Therefore, we can get the surface distribution of γ as a function of the final speed originating from this point. By interpolation γ to spherically uniform value on the source surface, we use this spatial distribution of γ in the energy equation to obtain a self-consistent, steady state MHD solution for the solar corona. We present the model result for different Carrington Rotations.

The value of swarm data for practical modeling of plasma devices

NASA Astrophysics Data System (ADS)

Napartovich, A. P.; Kochetov, I. V.

2011-04-01

The non-thermal plasma is a key component in gas lasers, waste gas cleaners, ozone generators, plasma igniters, flame holders, flow control in high-speed aerodynamics and other applications. The specific feature of the non-thermal plasma is its high sensitivity to variations in governing parameters (gas composition, pressure, pulse duration, E/N parameter). The reactivity of the plasma is due to the appearance of atoms and chemical radicals. For the efficient production of chemically active species high average electron energy is required, which is controlled by the balance of gain from the electric field and loss in inelastic collisions. In low-ionized plasma the electron energy distribution function is far from Maxwellian and must be found numerically for specified conditions. Numerical modeling of processes in plasma technologies requires vast databases on electron scattering cross sections to be available. The only reliable criterion for evaluations of validity of a set of cross sections for a particular species is a correct prediction of electron transport and kinetic coefficients measured in swarm experiments. This criterion is used traditionally to improve experimentally measured cross sections, as was suggested earlier by Phelps. The set of cross sections subjected to this procedure is called a self-consistent set. Nowadays, such reliable self-consistent sets are known for many species. Problems encountered in implementation of the fitting procedure and examples of its successful applications are described in the paper. .

Active colloids in the context of chemical kinetics

NASA Astrophysics Data System (ADS)

Oshanin, G.; Popescu, M. N.; Dietrich, S.

2017-03-01

We study a mesoscopic model of a chemically active colloidal particle which on certain parts of its surface promotes chemical reactions in the surrounding solution. For reasons of simplicity and conceptual clarity, we focus on the case in which only electrically neutral species are present in the solution and on chemical reactions which are described by first order kinetics. Within a self-consistent approach we explicitly determine the steady state product and reactant number density fields around the colloid as functionals of the interaction potentials of the various molecular species in solution with the colloid. By using a reciprocal theorem, this allows us to compute and to interpret—in a transparent way in terms of the classical Smoluchowski theory of chemical kinetics—the external force needed to keep such a catalytically active colloid at rest (stall force) or, equivalently, the corresponding velocity of the colloid if it is free to move. We use the particular case of triangular-well interaction potentials as a benchmark example for applying the general theoretical framework developed here. For this latter case, we derive explicit expressions for the dependences of the quantities of interest on the diffusion coefficients of the chemical species, the reaction rate constant, the coverage by catalyst, the size of the colloid, as well as on the parameters of the interaction potentials. These expressions provide a detailed picture of the phenomenology associated with catalytically-active colloids and self-diffusiophoresis.

Peierls-Nabarro barrier and protein loop propagation

NASA Astrophysics Data System (ADS)

Sieradzan, Adam K.; Niemi, Antti; Peng, Xubiao

2014-12-01

When a self-localized quasiparticle excitation propagates along a discrete one-dimensional lattice, it becomes subject to a dissipation that converts the kinetic energy into lattice vibrations. Eventually the kinetic energy no longer enables the excitation to cross over the minimum energy barrier between neighboring sites, and the excitation becomes localized within a lattice cell. In the case of a protein, the lattice structure consists of the Cα backbone. The self-localized quasiparticle excitation is the elemental building block of loops. It can be modeled by a kink that solves a variant of the discrete nonlinear Schrödinger equation. We study the propagation of such a kink in the case of the protein G related albumin-binding domain, using the united residue coarse-grained molecular-dynamics force field. We estimate the height of the energy barriers that the kink needs to cross over in order to propagate along the backbone lattice. We analyze how these barriers give rise to both stresses and reliefs, which control the kink movement. For this, we deform a natively folded protein structure by parallel translating the kink along the backbone away from its native position. We release the transposed kink, and we follow how it propagates along the backbone toward the native location. We observe that the dissipative forces that are exerted on the kink by the various energy barriers have a pivotal role in determining how a protein folds toward its native state.

Swarming in viscous fluids: three-dimensional patterns in swimmer- and force-induced flows

NASA Astrophysics Data System (ADS)

Chuang, Yao-Li; D'Orsogna, Maria R.; Chou, Tom

Mathematical models of self-propelled interacting particles have reproduced various fascinating ``swarming'' patterns observed in natural and artificial systems. The formulation of such models usually ignores the influence of the surrounding medium in which the particles swarm. Here we develop from first principles a three-dimensional theory of swarming particles in a viscous fluid environment and investigate how the hydrodynamic coupling among the particles may affect their collective behavior. Specifically, we examine the hydrodynamic coupling among self-propelled particles interacting through ``social'' or ``mechanical'' forces. We discover that new patterns arise as a consequence of different interactions and self-propulsion mechanisms. Examples include flocks with prolate or oblate shapes, intermittent mills, recirculating peloton-like structures, and jet-like fluid flows that kinetically destabilize mill-like structures. Our results reveal possible mechanisms for three-dimensional swarms to kinetically control their collective behaviors in fluids. Supported by NSF DMS 1021818 & 1021850, ARO W1911NF-14-1-0472, ARO MURI W1911NF-11-10332.

Microscopic Modeling of Intersubband Optical Processes in Type II Semiconductor Quantum Wells: Linear Absorption

NASA Technical Reports Server (NTRS)

Li, Jian-Zhong; Kolokolov, Kanstantin I.; Ning, Cun-Zheng

2003-01-01

Linear absorption spectra arising from intersubband transitions in semiconductor quantum well heterostructures are analyzed using quantum kinetic theory by treating correlations to the first order within Hartree-Fock approximation. The resulting intersubband semiconductor Bloch equations take into account extrinsic dephasing contributions, carrier-longitudinal optical phonon interaction and carrier-interface roughness interaction which is considered with Ando s theory. As input for resonance lineshape calculation, a spurious-states-free 8-band kp Hamiltonian is used, in conjunction with the envelop function approximation, to compute self-consistently the energy subband structure of electrons in type II InAs/AlSb single quantum well structures. We demonstrate the interplay of nonparabolicity and many-body effects in the mid-infrared frequency range for such heterostructures.

Numerical Solution of the Gyrokinetic Poisson Equation in TEMPEST

NASA Astrophysics Data System (ADS)

Dorr, Milo; Cohen, Bruce; Cohen, Ronald; Dimits, Andris; Hittinger, Jeffrey; Kerbel, Gary; Nevins, William; Rognlien, Thomas; Umansky, Maxim; Xiong, Andrew; Xu, Xueqiao

2006-10-01

The gyrokinetic Poisson (GKP) model in the TEMPEST continuum gyrokinetic edge plasma code yields the electrostatic potential due to the charge density of electrons and an arbitrary number of ion species including the effects of gyroaveraging in the limit kρ1. The TEMPEST equations are integrated as a differential algebraic system involving a nonlinear system solve via Newton-Krylov iteration. The GKP preconditioner block is inverted using a multigrid preconditioned conjugate gradient (CG) algorithm. Electrons are treated as kinetic or adiabatic. The Boltzmann relation in the adiabatic option employs flux surface averaging to maintain neutrality within field lines and is solved self-consistently with the GKP equation. A decomposition procedure circumvents the near singularity of the GKP Jacobian block that otherwise degrades CG convergence.

Impact of velocity space distribution on hybrid kinetic-magnetohydrodynamic simulation of the (1,1) mode

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kim, Charlson C.

2008-07-15

Numeric studies of the impact of the velocity space distribution on the stabilization of (1,1) internal kink mode and excitation of the fishbone mode are performed with a hybrid kinetic-magnetohydrodynamic model. These simulations demonstrate an extension of the physics capabilities of NIMROD[C. R. Sovinec et al., J. Comput. Phys. 195, 355 (2004)], a three-dimensional extended magnetohydrodynamic (MHD) code, to include the kinetic effects of an energetic minority ion species. Kinetic effects are captured by a modification of the usual MHD momentum equation to include a pressure tensor calculated from the {delta}f particle-in-cell method [S. E. Parker and W. W. Lee,more » Phys. Fluids B 5, 77 (1993)]. The particles are advanced in the self-consistent NIMROD fields. We outline the implementation and present simulation results of energetic minority ion stabilization of the (1,1) internal kink mode and excitation of the fishbone mode. A benchmark of the linear growth rate and real frequency is shown to agree well with another code. The impact of the details of the velocity space distribution is examined; particularly extending the velocity space cutoff of the simulation particles. Modestly increasing the cutoff strongly impacts the (1,1) mode. Numeric experiments are performed to study the impact of passing versus trapped particles. Observations of these numeric experiments suggest that assumptions of energetic particle effects should be re-examined.« less

Smouldering (thermal) remediation of soil contaminated with industrial organic liquids: novel insights into heat transfer and kinetics uncovered by integrating experiments and modelling

NASA Astrophysics Data System (ADS)

Gerhard, J.; Zanoni, M. A. B.; Torero, J. L.

2017-12-01

Smouldering (i.e., flameless combustion) underpins the technology Self-sustaining Treatment for Active Remediation (STAR). STAR achieves the in situ destruction of nonaqueous phase liquids (NAPLs) by generating a self-sustained smouldering reaction that propagates through the source zone. This research explores the nature of the travelling reaction and the influence of key in situ and engineered characteristics. A novel one-dimensional numerical model was developed (in COMSOL) to simulate the smouldering remediation of bitumen-contaminated sand. This model was validated against laboratory column experiments. Achieving model validation depended on correctly simulating the energy balance at the reaction front, including properly accounting for heat transfer, smouldering kinetics, and heat losses. Heat transfer between soil and air was demonstrated to be generally not at equilibrium. Moreover, existing heat transfer correlations were found to be inappropriate for the low air flow Reynold's numbers (Re < 30) relevant in this and similar thermal remediation systems. Therefore, a suite of experiments were conducted to generate a new heat transfer correlation, which generated correct simulations of convective heat flow through soil. Moreover, it was found that, for most cases of interest, a simple two-step pyrolysis/oxidation set of kinetic reactions was sufficient. Arrhenius parameters, calculated independently from thermogravimetric experiments, allowed the reaction kinetics to be validated in the smouldering model. Furthermore, a simple heat loss term sufficiently accounted for radial heat losses from the column. Altogether, these advances allow this simple model to reasonably predict the self-sustaining process including the peak reaction temperature, the reaction velocity, and the complete destruction of bitumen behind the front. Simulations with the validated model revealed numerous unique insights, including how the system inherently recycles energy, how air flow rate and NAPL saturation dictate contaminant destruction rates, and the extremes that lead to extinction. Overall, this research provides unique insights into the complex interplay of thermochemical processes that govern the success of smouldering as well as other thermal remediation approaches.

A Monster CME Obscuring a Demon Star Flare

NASA Astrophysics Data System (ADS)

Moschou, Sofia-Paraskevi; Drake, Jeremy J.; Cohen, Ofer; Alvarado-Gomez, Julian D.; Garraffo, Cecilia

2017-12-01

We explore the scenario of a coronal mass ejection (CME) being the cause of the observed continuous X-ray absorption of the 1997 August 30 superflare on the eclipsing binary Algol (the Demon Star). The temporal decay of the absorption is consistent with absorption by a CME undergoing self-similar evolution with uniform expansion velocity. We investigate the kinematic and energetic properties of the CME using the ice cream cone model for its three-dimensional structure in combination with the observed profile of the hydrogen column density decline with time. Different physically justified length scales were used that allowed us to estimate lower and upper limits of the possible CME characteristics. Further consideration of the maximum available magnetic energy in starspots leads us to quantify its mass as likely lying in the range 2× {10}21 {--} 2× {10}22 g and kinetic energy in the range 7× {10}35 {--} 3× {10}38 erg. The results are in reasonable agreement with extrapolated relations between flare X-ray fluence and CME mass and kinetic energy derived for solar CMEs.

Variational formulation of macroparticle models for electromagnetic plasma simulations

DOE PAGES

Stamm, Alexander B.; Shadwick, Bradley A.; Evstatiev, Evstati G.

2014-06-01

A variational method is used to derive a self-consistent macroparticle model for relativistic electromagnetic kinetic plasma simulations. Extending earlier work, discretization of the electromagnetic Low Lagrangian is performed via a reduction of the phase-space distribution function onto a collection of finite-sized macroparticles of arbitrary shape and discretization of field quantities onto a spatial grid. This approach may be used with lab frame coordinates or moving window coordinates; the latter can greatly improve computational efficiency for studying some types of laser-plasma interactions. The primary advantage of the variational approach is the preservation of Lagrangian symmetries, which in our case leads tomore » energy conservation and thus avoids difficulties with grid heating. In addition, this approach decouples particle size from grid spacing and relaxes restrictions on particle shape, leading to low numerical noise. The variational approach also guarantees consistent approximations in the equations of motion and is amenable to higher order methods in both space and time. We restrict our attention to the 1.5-D case (one coordinate and two momenta). Lastly, simulations are performed with the new models and demonstrate energy conservation and low noise.« less

Mechanisms of kinetic stabilization by the drugs paclitaxel and vinblastine

PubMed Central

Castle, Brian T.; McCubbin, Seth; Prahl, Louis S.; Bernens, Jordan N.; Sept, David; Odde, David J.

2017-01-01

Microtubule-targeting agents (MTAs), widely used as biological probes and chemotherapeutic drugs, bind directly to tubulin subunits and “kinetically stabilize” microtubules, suppressing the characteristic self-assembly process of dynamic instability. However, the molecular-level mechanisms of kinetic stabilization are unclear, and the fundamental thermodynamic and kinetic requirements for dynamic instability and its elimination by MTAs have yet to be defined. Here we integrate a computational model for microtubule assembly with nanometer-scale fluorescence microscopy measurements to identify the kinetic and thermodynamic basis of kinetic stabilization by the MTAs paclitaxel, an assembly promoter, and vinblastine, a disassembly promoter. We identify two distinct modes of kinetic stabilization in live cells, one that truly suppresses on-off kinetics, characteristic of vinblastine, and the other a “pseudo” kinetic stabilization, characteristic of paclitaxel, that nearly eliminates the energy difference between the GTP- and GDP-tubulin thermodynamic states. By either mechanism, the main effect of both MTAs is to effectively stabilize the microtubule against disassembly in the absence of a robust GTP cap. PMID:28298489

Numerical simulation of flood inundation using a well-balanced kinetic scheme for the shallow water equations with bulk recharge and discharge

NASA Astrophysics Data System (ADS)

Ersoy, Mehmet; Lakkis, Omar; Townsend, Philip

2016-04-01

The flow of water in rivers and oceans can, under general assumptions, be efficiently modelled using Saint-Venant's shallow water system of equations (SWE). SWE is a hyperbolic system of conservation laws (HSCL) which can be derived from a starting point of incompressible Navier-Stokes. A common difficulty in the numerical simulation of HSCLs is the conservation of physical entropy. Work by Audusse, Bristeau, Perthame (2000) and Perthame, Simeoni (2001), proposed numerical SWE solvers known as kinetic schemes (KSs), which can be shown to have desirable entropy-consistent properties, and are thus called well-balanced schemes. A KS is derived from kinetic equations that can be integrated into the SWE. In flood risk assessment models the SWE must be coupled with other equations describing interacting meteorological and hydrogeological phenomena such as rain and groundwater flows. The SWE must therefore be appropriately modified to accommodate source and sink terms, so kinetic schemes are no longer valid. While modifications of SWE in this direction have been recently proposed, e.g., Delestre (2010), we depart from the extant literature by proposing a novel model that is "entropy-consistent" and naturally extends the SWE by respecting its kinetic formulation connections. This allows us to derive a system of partial differential equations modelling flow of a one-dimensional river with both a precipitation term and a groundwater flow model to account for potential infiltration and recharge. We exhibit numerical simulations of the corresponding kinetic schemes. These simulations can be applied to both real world flood prediction and the tackling of wider issues on how climate and societal change are affecting flood risk.

Critical evaluation of Jet-A spray combustion using propane chemical kinetics in gas turbine combustion simulated by KIVA-2

NASA Technical Reports Server (NTRS)

Nguyen, H. L.; Ying, S.-J.

1990-01-01

Jet-A spray combustion has been evaluated in gas turbine combustion with the use of propane chemical kinetics as the first approximation for the chemical reactions. Here, the numerical solutions are obtained by using the KIVA-2 computer code. The KIVA-2 code is the most developed of the available multidimensional combustion computer programs for application of the in-cylinder combustion dynamics of internal combustion engines. The released version of KIVA-2 assumes that 12 chemical species are present; the code uses an Arrhenius kinetic-controlled combustion model governed by a four-step global chemical reaction and six equilibrium reactions. Researchers efforts involve the addition of Jet-A thermophysical properties and the implementation of detailed reaction mechanisms for propane oxidation. Three different detailed reaction mechanism models are considered. The first model consists of 131 reactions and 45 species. This is considered as the full mechanism which is developed through the study of chemical kinetics of propane combustion in an enclosed chamber. The full mechanism is evaluated by comparing calculated ignition delay times with available shock tube data. However, these detailed reactions occupy too much computer memory and CPU time for the computation. Therefore, it only serves as a benchmark case by which to evaluate other simplified models. Two possible simplified models were tested in the existing computer code KIVA-2 for the same conditions as used with the full mechanism. One model is obtained through a sensitivity analysis using LSENS, the general kinetics and sensitivity analysis program code of D. A. Bittker and K. Radhakrishnan. This model consists of 45 chemical reactions and 27 species. The other model is based on the work published by C. K. Westbrook and F. L. Dryer.

Development of an Austenitization Kinetics Model for 22MnB5 Steel

NASA Astrophysics Data System (ADS)

Di Ciano, M.; Field, N.; Wells, M. A.; Daun, K. J.

2018-03-01

This paper presents a first-order austenitization kinetics model for 22MnB5 steel, commonly used in hot forming die quenching. Model parameters are derived from constant heating rate dilatometry measurements. Vickers hardness measurements made on coupons that were quenched at intermediate stages of the process were used to verify the model, and the Ac 1 and Ac 3 temperatures inferred from dilatometry are consistent with correlations found in the literature. The austenitization model was extended to consider non-constant heating rates typical of industrial furnaces and again showed reasonable agreement between predictions and measurements. Finally, the model is used to predict latent heat evolution during industrial heating and is shown to be consistent with values inferred from thermocouple measurements of furnace-heated 22MnB5 coupons reported in the literature.

Thermodynamic criteria for estimating the kinetic parameters of catalytic reactions

NASA Astrophysics Data System (ADS)

Mitrichev, I. I.; Zhensa, A. V.; Kol'tsova, E. M.

2017-01-01

Kinetic parameters are estimated using two criteria in addition to the traditional criterion that considers the consistency between experimental and modeled conversion data: thermodynamic consistency and the consistency with entropy production (i.e., the absolute rate of the change in entropy due to exchange with the environment is consistent with the rate of entropy production in the steady state). A special procedure is developed and executed on a computer to achieve the thermodynamic consistency of a set of kinetic parameters with respect to both the standard entropy of a reaction and the standard enthalpy of a reaction. A problem of multi-criterion optimization, reduced to a single-criterion problem by summing weighted values of the three criteria listed above, is solved. Using the reaction of NO reduction with CO on a platinum catalyst as an example, it is shown that the set of parameters proposed by D.B. Mantri and P. Aghalayam gives much worse agreement with experimental values than the set obtained on the basis of three criteria: the sum of the squares of deviations for conversion, the thermodynamic consistency, and the consistency with entropy production.

The Roles of Biofilm Conductivity and Donor Substrate Kinetics in a Mixed-Culture Biofilm Anod

EPA Science Inventory

We experimentally assessed kinetics and thermodynamics of electron transfer (ET) from the donor substrate (acetate) to the anode for a mixed-culture biofilm anode. We interpreted the results with a modified biofilm-conduction model consisting of three ET steps: (1) intracellular...

Influence of thermodynamically unfavorable secondary structures on DNA hybridization kinetics

PubMed Central

Hata, Hiroaki; Kitajima, Tetsuro

2018-01-01

Abstract Nucleic acid secondary structure plays an important role in nucleic acid–nucleic acid recognition/hybridization processes, and is also a vital consideration in DNA nanotechnology. Although the influence of stable secondary structures on hybridization kinetics has been characterized, unstable secondary structures, which show positive ΔG° with self-folding, can also form, and their effects have not been systematically investigated. Such thermodynamically unfavorable secondary structures should not be ignored in DNA hybridization kinetics, especially under isothermal conditions. Here, we report that positive ΔG° secondary structures can change the hybridization rate by two-orders of magnitude, despite the fact that their hybridization obeyed second-order reaction kinetics. The temperature dependence of hybridization rates showed non-Arrhenius behavior; thus, their hybridization is considered to be nucleation limited. We derived a model describing how ΔG° positive secondary structures affect hybridization kinetics in stopped-flow experiments with 47 pairs of oligonucleotides. The calculated hybridization rates, which were based on the model, quantitatively agreed with the experimental rate constant. PMID:29220504

Monte carlo simulations of enzyme reactions in two dimensions: fractal kinetics and spatial segregation.

PubMed

Berry, Hugues

2002-10-01

Conventional equations for enzyme kinetics are based on mass-action laws, that may fail in low-dimensional and disordered media such as biological membranes. We present Monte Carlo simulations of an isolated Michaelis-Menten enzyme reaction on two-dimensional lattices with varying obstacle densities, as models of biological membranes. The model predicts that, as a result of anomalous diffusion on these low-dimensional media, the kinetics are of the fractal type. Consequently, the conventional equations for enzyme kinetics fail to describe the reaction. In particular, we show that the quasi-stationary-state assumption can hardly be retained in these conditions. Moreover, the fractal characteristics of the kinetics are increasingly pronounced as obstacle density and initial substrate concentration increase. The simulations indicate that these two influences are mainly additive. Finally, the simulations show pronounced S-P segregation over the lattice at obstacle densities compatible with in vivo conditions. This phenomenon could be a source of spatial self organization in biological membranes.

Monte carlo simulations of enzyme reactions in two dimensions: fractal kinetics and spatial segregation.

PubMed Central

Berry, Hugues

2002-01-01

Conventional equations for enzyme kinetics are based on mass-action laws, that may fail in low-dimensional and disordered media such as biological membranes. We present Monte Carlo simulations of an isolated Michaelis-Menten enzyme reaction on two-dimensional lattices with varying obstacle densities, as models of biological membranes. The model predicts that, as a result of anomalous diffusion on these low-dimensional media, the kinetics are of the fractal type. Consequently, the conventional equations for enzyme kinetics fail to describe the reaction. In particular, we show that the quasi-stationary-state assumption can hardly be retained in these conditions. Moreover, the fractal characteristics of the kinetics are increasingly pronounced as obstacle density and initial substrate concentration increase. The simulations indicate that these two influences are mainly additive. Finally, the simulations show pronounced S-P segregation over the lattice at obstacle densities compatible with in vivo conditions. This phenomenon could be a source of spatial self organization in biological membranes. PMID:12324410

Flow-Induced New Channels of Energy Exchange in Multi-Scale Plasma Dynamics – Revisiting Perturbative Hybrid Kinetic-MHD Theory

PubMed Central

Shiraishi, Junya; Miyato, Naoaki; Matsunaga, Go

2016-01-01

It is found that new channels of energy exchange between macro- and microscopic dynamics exist in plasmas. They are induced by macroscopic plasma flow. This finding is based on the kinetic-magnetohydrodynamic (MHD) theory, which analyses interaction between macroscopic (MHD-scale) motion and microscopic (particle-scale) dynamics. The kinetic-MHD theory is extended to include effects of macroscopic plasma flow self-consistently. The extension is realised by generalising an energy exchange term due to wave-particle resonance, denoted by δ WK. The first extension is generalisation of the particle’s Lagrangian, and the second one stems from modification to the particle distribution function due to flow. These extensions lead to a generalised expression of δ WK, which affects the MHD stability of plasmas. PMID:27160346

Branch-Based Model for the Diameters of the Pulmonary Airways: Accounting for Departures From Self-Consistency and Registration Errors

PubMed Central

Neradilek, Moni B.; Polissar, Nayak L.; Einstein, Daniel R.; Glenny, Robb W.; Minard, Kevin R.; Carson, James P.; Jiao, Xiangmin; Jacob, Richard E.; Cox, Timothy C.; Postlethwait, Edward M.; Corley, Richard A.

2017-01-01

We examine a previously published branch-based approach for modeling airway diameters that is predicated on the assumption of self-consistency across all levels of the tree. We mathematically formulate this assumption, propose a method to test it and develop a more general model to be used when the assumption is violated. We discuss the effect of measurement error on the estimated models and propose methods that take account of error. The methods are illustrated on data from MRI and CT images of silicone casts of two rats, two normal monkeys, and one ozone-exposed monkey. Our results showed substantial departures from self-consistency in all five subjects. When departures from self-consistency exist, we do not recommend using the self-consistency model, even as an approximation, as we have shown that it may likely lead to an incorrect representation of the diameter geometry. The new variance model can be used instead. Measurement error has an important impact on the estimated morphometry models and needs to be addressed in the analysis. PMID:22528468

Controlling Self-Assembly in Al(110) Homoepitaxy

NASA Astrophysics Data System (ADS)

Tiwary, Yogesh; Fichthorn, Kristen

2010-03-01

Homoepitaxial growth on Al(110) exhibits nanoscale self-assembly into huts with well-defined (100) and (111) facets [1]. Although some of the diffusion mechanisms underlying this kinetic self-assembly were identified and incorporated into a two-dimensional model [2], we used density-functional theory (DFT) to identify many other mechanisms that are needed to describe the three-dimensional assembly seen experimentally [3]. We developed a three-dimensional kinetic Monte Carlo (KMC) model of Al(110) homoepitaxy. The inputs to the model were obtained from DFT [3,4]. Our model is in agreement with experimentally observed trends for this system. We used KMC to predict self-assembly under various growth conditions. To achieve precise placement of Al nanohuts, we simulated thermal-field-directed assembly [5]. Our results indicate that this technique can be used to create uniform arrays of nanostructures. [1] F. Buatier de Mongeot, W. Zhu, A. Molle, R. Buzio, C. Boragno, U. Valbusa, E. Wang, and Z. Zhang, Phys. Rev. Lett. 91, 016102 (2003). [2] W. Zhu, F. Buatier de Mongeot, U. Valbusa, E. G. Wang, and Z. Y. Zhang, Phys. Rev. Lett. 92, 106102 (2004). [3] Y. Tiwary and K. A. Fichthorn, submitted to Phys. Rev. B. [4] Y. Tiwary and K. A. Fichthorn, Phys. Rev. B 78, 205418 (2008). [5] C. Zhang and R. Kalyanaraman, Appl. Phys. Lett. 83, 4827 (2003).

Growth Kinetics of Intracellular RNA/Protein Droplets: Signature of a Liquid-Liquid Phase Transition?

NASA Astrophysics Data System (ADS)

Berry, Joel; Weber, Stephanie C.; Vaidya, Nilesh; Zhu, Lian; Haataja, Mikko; Brangwynne, Clifford P.

2015-03-01

Nonmembrane-bound organelles are functional, dynamic assemblies of RNA and/or protein that can self-assemble and disassemble within the cytoplasm or nucleoplasm. The possibility that underlying intracellular phase transitions may drive and mediate the morphological evolution of some membrane-less organelles has been supported by several recent studies. In this talk, results from a collaborative experimental-theoretical study of the growth and dissolution kinetics of nucleoli and extranucleolar droplets (ENDs) in C. elegans embryos will be presented. We have employed Flory-Huggins solution theory, reaction-diffusion kinetics, and quantitative statistical dynamic scaling analysis to characterize the specific growth mechanisms at work. Our findings indicate that both in vivo and in vitro droplet scaling and growth kinetics are consistent with those resulting from an equilibrium liquid-liquid phase transition mediated by passive nonequilibrium growth mechanisms - simultaneous Brownian coalescence and Ostwald ripening. This supports a view in which cells can employ phase transitions to drive structural organization, while utilizing active processes, such as local transcriptional activity, to fine tune the kinetics of these phase transitions in response to given conditions.

Global Energetics of Solar Flares. VI. Refined Energetics of Coronal Mass Ejections

NASA Astrophysics Data System (ADS)

Aschwanden, Markus J.

2017-09-01

In this study, we refine the coronal mass ejection (CME) model that was presented in an earlier study of the global energetics of solar flares and associated CMEs and apply it to all (860) GOES M- and X-class flare events observed during the first seven years (2010-2016) of the Solar Dynamics Observatory (SDO) mission. The model refinements include (1) the CME geometry in terms of a 3D volume undergoing self-similar adiabatic expansion, (2) the solar gravitational deceleration during the propagation of the CME, which discriminates between eruptive and confined CMEs, (3) a self-consistent relationship between the CME center-of-mass motion detected during EUV dimming and the leading-edge motion observed in white-light coronagraphs, (4) the equipartition of the CME’s kinetic and thermal energies, and (5) the Rosner-Tucker-Vaiana scaling law. The refined CME model is entirely based on EUV-dimming observations (using Atmospheric Imager Assembly (AIA)/SDO data) and complements the traditional white-light scattering model (using Large-Angle and Spectrometric Coronagraph Experiment (LASCO)/Solar and Heliospheric Observatory data), and both models are independently capable of determining fundamental CME parameters. Comparing the two methods, we find that (1) LASCO is less sensitive than AIA in detecting CMEs (in 24% of the cases), (2) CME masses below {m}{cme}≲ {10}14 g are underestimated by LASCO, (3) AIA and LASCO masses, speeds, and energies agree closely in the statistical mean after the elimination of outliers, and (4) the CME parameters speed v, emission measure-weighted flare peak temperature T e , and length scale L are consistent with the following scaling laws: v\\propto {T}e1/2, v\\propto {({m}{cme})}1/4, and {m}{cme}\\propto {L}2.

Detailed kinetics and regulation of mammalian 2-oxoglutarate dehydrogenase

PubMed Central

2011-01-01

Background Mitochondrial 2-oxoglutarate (α-ketoglutarate) dehydrogenase complex (OGDHC), a key regulatory point of tricarboxylic acid (TCA) cycle, plays vital roles in multiple pathways of energy metabolism and biosynthesis. The catalytic mechanism and allosteric regulation of this large enzyme complex are not fully understood. Here computer simulation is used to test possible catalytic mechanisms and mechanisms of allosteric regulation of the enzyme by nucleotides (ATP, ADP), pH, and metal ion cofactors (Ca2+ and Mg2+). Results A model was developed based on an ordered ter-ter enzyme kinetic mechanism combined with con-formational changes that involve rotation of one lipoic acid between three catalytic sites inside the enzyme complex. The model was parameterized using a large number of kinetic data sets on the activity of OGDHC, and validated by comparison of model predictions to independent data. Conclusions The developed model suggests a hybrid rapid-equilibrium ping-pong random mechanism for the kinetics of OGDHC, consistent with previously reported mechanisms, and accurately describes the experimentally observed regulatory effects of cofactors on the OGDHC activity. This analysis provides a single consistent theoretical explanation for a number of apparently contradictory results on the roles of phosphorylation potential, NAD (H) oxidation-reduction state ratio, as well as the regulatory effects of metal ions on ODGHC function. PMID:21943256

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hochberg, Yonit; Pyle, Matt; Zhao, Yue

We examine in greater detail the recent proposal of using superconductors for detecting dark matter as light as the warm dark matter limit of O(keV). Detection of suc light dark matter is possible if the entire kinetic energy of the dark matter is extracted in the scattering, and if the experiment is sensitive to O(meV) energy depositions. This is the case for Fermi-degenerate materials in which the Fermi velocity exceeds the dark matter velocity dispersion in the Milky Way of ~10 –3. We focus on a concrete experimental proposal using a superconducting target with a transition edge sensor in ordermore » to detect the small energy deposits from the dark matter scatterings. Considering a wide variety of constraints, from dark matter self-interactions to the cosmic microwave background, we show that models consistent with cosmological/astrophysical and terrestrial constraints are observable with such detectors. A wider range of viable models with dark matter mass below an MeV is available if dark matter or mediator properties (such as couplings or masses) differ at BBN epoch or in stellar interiors from those in superconductors. We also show that metal targets pay a strong in-medium suppression for kinetically mixed mediators; this suppression is alleviated with insulating targets.« less

Detecting superlight dark matter with Fermi-degenerate materials

DOE PAGES

Hochberg, Yonit; Pyle, Matt; Zhao, Yue; ...

2016-08-08

We examine in greater detail the recent proposal of using superconductors for detecting dark matter as light as the warm dark matter limit of O(keV). Detection of suc light dark matter is possible if the entire kinetic energy of the dark matter is extracted in the scattering, and if the experiment is sensitive to O(meV) energy depositions. This is the case for Fermi-degenerate materials in which the Fermi velocity exceeds the dark matter velocity dispersion in the Milky Way of ~10 –3. We focus on a concrete experimental proposal using a superconducting target with a transition edge sensor in ordermore » to detect the small energy deposits from the dark matter scatterings. Considering a wide variety of constraints, from dark matter self-interactions to the cosmic microwave background, we show that models consistent with cosmological/astrophysical and terrestrial constraints are observable with such detectors. A wider range of viable models with dark matter mass below an MeV is available if dark matter or mediator properties (such as couplings or masses) differ at BBN epoch or in stellar interiors from those in superconductors. We also show that metal targets pay a strong in-medium suppression for kinetically mixed mediators; this suppression is alleviated with insulating targets.« less

Transition in Electron Physics of Magnetic Reconnection in Weakly Collisional Plasma

NASA Astrophysics Data System (ADS)

Le, A.; Roytershteyn, V.; Karimabadi, H.; Daughton, W. S.; Egedal, J.; Forest, C.

2013-12-01

Using self-consistent fully kinetic simulations with a Monte-Carlo treatment of the Coulomb collision operator, we explore the transition between collisional and kinetic regimes of magnetic reconnection in high-Lundquist-number current sheets. Recent research in collisionless reconnection has shown that electron kinetic physics plays a key role in the evolution. Large-scale electron current sheets may form, leading to secondary island formation and turbulent flux rope interactions in 3D. The new collisional simulations demonstrate how increasing collisionality modifies or eliminates these electron structures in the kinetic regimes. Additional basic questions that are addressed include how the reconnection rate and the release of magnetic energy into electrons and ions vary with collisionality. The numerical study provides insight into reconnection in dense regions of the solar corona, the solar wind, and upcoming laboratory experiments at MRX (Princeton) and MPDX (UW-Madison). The implications of these results for studies of turbulence dissipation in weakly collisional plasmas are discussed.

Molecular dynamics of reversible self-healing materials

NASA Astrophysics Data System (ADS)

Madden, Ian; Luijten, Erik

Hydrolyzable polymers have numerous industrial applications as degradable materials. Recent experimental work by Cheng and co-workers has introduced the concept of hindered urea bond (HUB) chemistry to design self-healing systems. Important control parameters are the steric hindrance of the HUB structures, which is used to tune the hydrolytic degradation kinetics, and their density. We employ molecular dynamics simulations of polymeric interfaces to systematically explore the role of these properties in a coarse-grained model, and make direct comparison to experimental data. Our model provides direct insight into the self-healing process, permitting optimization of the control parameters.

Multicomponent lattice Boltzmann model from continuum kinetic theory.

PubMed

Shan, Xiaowen

2010-04-01

We derive from the continuum kinetic theory a multicomponent lattice Boltzmann model with intermolecular interaction. The resulting model is found to be consistent with the model previously derived from a lattice-gas cellular automaton [X. Shan and H. Chen, Phys. Rev. E 47, 1815 (1993)] but applies in a much broader domain. A number of important insights are gained from the kinetic theory perspective. First, it is shown that even in the isothermal case, the energy equipartition principle dictates the form of the equilibrium distribution function. Second, thermal diffusion is shown to exist and the corresponding diffusivities are given in terms of macroscopic parameters. Third, the ordinary diffusion is shown to satisfy the Maxwell-Stefan equation at the ideal-gas limit.

Kinetics of hairpin ribozyme cleavage in yeast.

PubMed Central

Donahue, C P; Fedor, M J

1997-01-01

Hairpin ribozymes catalyze a self-cleavage reaction that provides a simple model for quantitative analyses of intracellular mechanisms of RNA catalysis. Decay rates of chimeric mRNAs containing self-cleaving ribozymes give a direct measure of intracellular cleavage kinetics in yeast. Intracellular ribozyme-mediated cleavage occurs at similar rates and shows similar inhibition by ribozyme mutations as ribozyme-mediated reactions in vitro, but only when ribozymes are located in a favorable mRNA sequence context. The impact of cleavage on mRNA abundance is shown to depend directly on intrinsic mRNA stability. Surprisingly, cleavage products are no more labile than uncleaved mRNAs despite the loss of terminal cap structures or poly (A). PMID:9292496

3-Dimensional Modeling of Capacitively and Inductively Coupled Plasma Etching Systems

NASA Astrophysics Data System (ADS)

Rauf, Shahid

2008-10-01

Low temperature plasmas are widely used for thin film etching during micro and nano-electronic device fabrication. Fluid and hybrid plasma models were developed 15-20 years ago to understand the fundamentals of these plasmas and plasma etching. These models have significantly evolved since then, and are now a major tool used for new plasma hardware design and problem resolution. Plasma etching is a complex physical phenomenon, where inter-coupled plasma, electromagnetic, fluid dynamics, and thermal effects all have a major influence. The next frontier in the evolution of fluid-based plasma models is where these models are able to self-consistently treat the inter-coupling of plasma physics with fluid dynamics, electromagnetics, heat transfer and magnetostatics. We describe one such model in this paper and illustrate its use in solving engineering problems of interest for next generation plasma etcher design. Our 3-dimensional plasma model includes the full set of Maxwell equations, transport equations for all charged and neutral species in the plasma, the Navier-Stokes equation for fluid flow, and Kirchhoff's equations for the lumped external circuit. This model also includes Monte Carlo based kinetic models for secondary electrons and stochastic heating, and can take account of plasma chemistry. This modeling formalism allows us to self-consistently treat the dynamics in commercial inductively and capacitively coupled plasma etching reactors with realistic plasma chemistries, magnetic fields, and reactor geometries. We are also able to investigate the influence of the distributed electromagnetic circuit at very high frequencies (VHF) on the plasma dynamics. The model is used to assess the impact of azimuthal asymmetries in plasma reactor design (e.g., off-center pump, 3D magnetic field, slit valve, flow restrictor) on plasma characteristics at frequencies from 2 -- 180 MHz. With Jason Kenney, Ankur Agarwal, Ajit Balakrishna, Kallol Bera, and Ken Collins.

INTERDISCIPLINARY PHYSICS AND RELATED AREAS OF SCIENCE AND TECHNOLOGY: Competition Between Self-birth and Catalyzed Death in Aggregation Growth with Catalysis Injection

NASA Astrophysics Data System (ADS)

Chen, Dan; Lin, Zhen-Quan; Sun, Yun-Fei; Ke, Jian-Hong

2009-12-01

We propose two irreversible aggregation growth models of aggregates of two distinct species (A and B) to study the interactions between virus aggregates and medicine efficacy aggregates in the virus-medicine cooperative evolution system. The A-species aggregates evolve driven by self monomer birth and B-species aggregate-catalyzed monomer death in model I and by self birth, catalyzed death, and self monomer exchange reactions in model II, while the catalyst B-species aggregates are assumed to be injected into the system sustainedly or at a periodic time-dependent rate. The kinetic behaviors of the A-species aggregates are investigated by the rate equation approach based on the mean-field theory with the self birth rate kernel IA(k) = Ik, catalyzed death rate kernel JAB(k) = Jk and self exchange rate kernel KA (k, l) = Kkl. The kinetic behaviors of the A-species aggregates are mainly dominated by the competition between the two effects of the self birth (with the effective rate I) and the catalyzed death (with the effective rate JB0), while the effects of the self exchanges of the A-species aggregates which appear in an effective rate KA0 play important roles in the cases of I > JB0 and I = JB0. The evolution behaviors of the total mass MA1(t) and the total aggregate number MA0 (t) are obtained, and the aggregate size distribution αk (t) of species A is found to approach a generalized scaling form in the case of I >= JB0 and a special modified scaling form in the case of I < JB0. The periodical evolution of the B-monomers concentration plays an exponential form of the periodic modulation.

Testing the nonlocal kinetic energy functional of an inhomogeneous, two-dimensional degenerate Fermi gas within the average density approximation

NASA Astrophysics Data System (ADS)

Towers, J.; van Zyl, B. P.; Kirkby, W.

2015-08-01

In a recent paper [B. P. van Zyl et al., Phys. Rev. A 89, 022503 (2014), 10.1103/PhysRevA.89.022503], the average density approximation (ADA) was implemented to develop a parameter-free, nonlocal kinetic energy functional to be used in the orbital-free density functional theory of an inhomogeneous, two-dimensional (2D) Fermi gas. In this work, we provide a detailed comparison of self-consistent calculations within the ADA with the exact results of the Kohn-Sham density functional theory and the elementary Thomas-Fermi (TF) approximation. We demonstrate that the ADA for the 2D kinetic energy functional works very well under a wide variety of confinement potentials, even for relatively small particle numbers. Remarkably, the TF approximation for the kinetic energy functional, without any gradient corrections, also yields good agreement with the exact kinetic energy for all confining potentials considered, although at the expense of the spatial and kinetic energy densities exhibiting poor pointwise agreement, particularly near the TF radius. Our findings illustrate that the ADA kinetic energy functional yields accurate results for both the local and global equilibrium properties of an inhomogeneous 2D Fermi gas, without the need for any fitting parameters.

Atomic-Scale Design of Iron Fischer-Tropsch Catalysts; A Combined Computational Chemistry, Experimental, and Microkinetic Modeling Approach

DOE Office of Scientific and Technical Information (OSTI.GOV)

Manos Mavrikakis; James Dumesic; Rahul Nabar

2008-09-29

This work focuses on (1) searching/summarizing published Fischer-Tropsch synthesis (FTS) mechanistic and kinetic studies of FTS reactions on iron catalysts; (2) preparation and characterization of unsupported iron catalysts with/without potassium/platinum promoters; (3) measurement of H{sub 2} and CO adsorption/dissociation kinetics on iron catalysts using transient methods; (3) analysis of the transient rate data to calculate kinetic parameters of early elementary steps in FTS; (4) construction of a microkinetic model of FTS on iron, and (5) validation of the model from collection of steady-state rate data for FTS on iron catalysts. Three unsupported iron catalysts and three alumina-supported iron catalysts weremore » prepared by non-aqueous-evaporative deposition (NED) or aqueous impregnation (AI) and characterized by chemisorption, BET, temperature-programmed reduction (TPR), extent-of-reduction, XRD, and TEM methods. These catalysts, covering a wide range of dispersions and metal loadings, are well-reduced and relatively thermally stable up to 500-600 C in H{sub 2} and thus ideal for kinetic and mechanistic studies. Kinetic parameters for CO adsorption, CO dissociation, and surface carbon hydrogenation on these catalysts were determined from temperature-programmed desorption (TPD) of CO and temperature programmed surface hydrogenation (TPSR), temperature-programmed hydrogenation (TPH), and isothermal, transient hydrogenation (ITH). A microkinetic model was constructed for the early steps in FTS on polycrystalline iron from the kinetic parameters of elementary steps determined experimentally in this work and from literature values. Steady-state rate data were collected in a Berty reactor and used for validation of the microkinetic model. These rate data were fitted to 'smart' Langmuir-Hinshelwood rate expressions derived from a sequence of elementary steps and using a combination of fitted steady-state parameters and parameters specified from the transient measurements. The results provide a platform for further development of microkinetic models of FTS on Fe and a basis for more precise modeling of FTS activity of Fe catalysts. Calculations using periodic, self-consistent Density Functional Theory (DFT) methods were performed on various realistic models of industrial, Fe-based FTS catalysts. Close-packed, most stable Fe(110) facet was analyzed and subsequently carbide formation was found to be facile leading to the choice of the FeC(110) model representing a Fe facet with a sub-surface C atom. The Pt adatom (Fe{sup Pt}(110)) was found to be the most stable model for our studies into Pt promotion and finally the role of steps was elucidated by recourse to the defected Fe(211) facet. Binding Energies(BEs), preferred adsorption sites and geometries for all FTS relevant stable species and intermediates were evaluated on each model catalyst facet. A mechanistic model (comprising of 32 elementary steps involving 19 species) was constructed and each elementary step therein was fully characterized with respect to its thermochemistry and kinetics. Kinetic calculations involved evaluation of the Minimum Energy Pathways (MEPs) and activation energies (barriers) for each step. Vibrational frequencies were evaluated for the preferred adsorption configuration of each species with the aim of evaluating entropy-changes, pre exponential factors and serving as a useful connection with experimental surface science techniques. Comparative analysis among these four facets revealed important trends in their relative behavior and roles in FTS catalysis. Overall the First Principles Calculations afforded us a new insight into FTS catalysis on Fe and modified-Fe catalysts.« less

Discovery of a Significant Acetone•Hydroperoxy Adduct Chaperone Effect and Its Impact on the Determination of Room Temperature Rate Constants for Acetonylperoxy/Hydroperoxy Self-Reactions and Cross Reaction Via Infrared Kinetic Spectroscopy.

NASA Astrophysics Data System (ADS)

Grieman, F. J.; Hui, A. O.; Okumura, M.; Sander, S. P.

2017-12-01

In order to model the upper troposphere/lower stratosphere in regions containing acetone properly, the kinetics of the acetonylperoxy/hydroperoxy self-reactions and cross reaction have been studied over a wide temperature range using Infrared Kinetic Spectroscopy. We report here the determination of different rate constants for the acetonylperoxy chemistry that we obtained at 298 K compared to currently accepted values. A considerable increase in the observed HO2 self-reaction rate constant due to rate enhancement via the chaperone effect from the reaction between HO2 and the (CH3)2CO•HO2 hydrogen-bonded adduct, even at room temperature, was discovered that was previously ignored. Correct determination of the acetonylperoxy and hydroperoxy kinetics must include this dependence of the HO2 self-reaction rate on acetone concentration. Via excimer laser flash photolysis to create the radical reactants, HO2 absorption was monitored in the infrared by diode laser wavelength modulation detection simultaneously with CH3C(O)CH2O2absorption monitored in the ultraviolet at 300 nm as a function of time. Resulting decay curves were fit concurrently first over a short time scale to obtain the rate constants minimizing subsequent product reactions. Modeling/fitting with a complete reaction scheme was then performed to refine the rate constants and test their veracity. Experiments were carried out over a variety of concentrations of acetone and methanol. Although no effect due to methanol concentration was found at room temperature, the rate constant for the hydroperoxy self-reaction was found to increase linearly with acetone concentration which is interpreted as the adduct being formed and resulting in a chaperone mechanism that enhances the self-reaction rate: (CH3)2CO·HO2 + HO2 → H2O2 + O2 + (CH3)2CO Including this effect, the resulting room temperature rate constants for the cross reaction and the acetonylperoxy self-reaction were found to be 2-3 times smaller than previously reported. This complex formation/chaperone mechanism is similar to that found for methanol, but different in that it occurs at room temperature. No precursor concentration dependence was found for the acetonylperoxy radical reactions. The equilibrium constant for the complex formation will also be presented.

On the biophysics and kinetics of toehold-mediated DNA strand displacement

PubMed Central

Srinivas, Niranjan; Ouldridge, Thomas E.; Šulc, Petr; Schaeffer, Joseph M.; Yurke, Bernard; Louis, Ard A.; Doye, Jonathan P. K.; Winfree, Erik

2013-01-01

Dynamic DNA nanotechnology often uses toehold-mediated strand displacement for controlling reaction kinetics. Although the dependence of strand displacement kinetics on toehold length has been experimentally characterized and phenomenologically modeled, detailed biophysical understanding has remained elusive. Here, we study strand displacement at multiple levels of detail, using an intuitive model of a random walk on a 1D energy landscape, a secondary structure kinetics model with single base-pair steps and a coarse-grained molecular model that incorporates 3D geometric and steric effects. Further, we experimentally investigate the thermodynamics of three-way branch migration. Two factors explain the dependence of strand displacement kinetics on toehold length: (i) the physical process by which a single step of branch migration occurs is significantly slower than the fraying of a single base pair and (ii) initiating branch migration incurs a thermodynamic penalty, not captured by state-of-the-art nearest neighbor models of DNA, due to the additional overhang it engenders at the junction. Our findings are consistent with previously measured or inferred rates for hybridization, fraying and branch migration, and they provide a biophysical explanation of strand displacement kinetics. Our work paves the way for accurate modeling of strand displacement cascades, which would facilitate the simulation and construction of more complex molecular systems. PMID:24019238

On the biophysics and kinetics of toehold-mediated DNA strand displacement.

PubMed

Srinivas, Niranjan; Ouldridge, Thomas E; Sulc, Petr; Schaeffer, Joseph M; Yurke, Bernard; Louis, Ard A; Doye, Jonathan P K; Winfree, Erik

2013-12-01

Dynamic DNA nanotechnology often uses toehold-mediated strand displacement for controlling reaction kinetics. Although the dependence of strand displacement kinetics on toehold length has been experimentally characterized and phenomenologically modeled, detailed biophysical understanding has remained elusive. Here, we study strand displacement at multiple levels of detail, using an intuitive model of a random walk on a 1D energy landscape, a secondary structure kinetics model with single base-pair steps and a coarse-grained molecular model that incorporates 3D geometric and steric effects. Further, we experimentally investigate the thermodynamics of three-way branch migration. Two factors explain the dependence of strand displacement kinetics on toehold length: (i) the physical process by which a single step of branch migration occurs is significantly slower than the fraying of a single base pair and (ii) initiating branch migration incurs a thermodynamic penalty, not captured by state-of-the-art nearest neighbor models of DNA, due to the additional overhang it engenders at the junction. Our findings are consistent with previously measured or inferred rates for hybridization, fraying and branch migration, and they provide a biophysical explanation of strand displacement kinetics. Our work paves the way for accurate modeling of strand displacement cascades, which would facilitate the simulation and construction of more complex molecular systems.

Phase-field simulations of velocity selection in rapidly solidified binary alloys

NASA Astrophysics Data System (ADS)

Fan, Jun; Greenwood, Michael; Haataja, Mikko; Provatas, Nikolas

2006-09-01

Time-dependent simulations of two-dimensional isothermal Ni-Cu dendrites are simulated using a phase-field model solved with a finite-difference adaptive mesh refinement technique. Dendrite tip velocity selection is examined and found to exhibit a transition between two markedly different regimes as undercooling is increased. At low undercooling, the dendrite tip growth rate is consistent with the kinetics of the classical Stefan problem, where the interface is assume to be in local equilibrium. At high undercooling, the growth velocity selected approaches a linear dependence on melt undercooling, consistent with the continuous growth kinetics of Aziz and with a one-dimensional steady-state phase-field asymptotic analysis of Ahmad [Phys. Rev. E 58, 3436 (1998)]. Our simulations are also consistent with other previously observed behaviors of dendritic growth as undercooling is increased. These include the transition of dendritic morphology to absolute stability and nonequilibrium solute partitioning. Our results show that phase-field models of solidification, which inherently contain a nonzero interface width, can be used to study the dynamics of complex solidification phenomena involving both equilibrium and nonequilibrium interface growth kinetics.

Self-discharge of electrochemical capacitors based on soluble or grafted quinone.

PubMed

Shul, Galyna; Bélanger, Daniel

2016-07-28

The self-discharge of hybrid electrochemical capacitors based on the redox activity of electrolyte additives or grafted species to the electrode material is investigated simultaneously for the cell and each individual electrode. Electrochemical capacitors using a redox-active electrolyte consisting in hydroquinone added to the electrolyte solution and a redox-active electrode based on anthraquinone-grafted carbon as a negative electrode are investigated. The results are analyzed by using Conway kinetic models and compared to those of a common electrochemical double layer capacitor. The self-discharge investigation is complemented by charge/discharge cycling and it is shown that processes affecting galvanostatic charge/discharge cycling and the self-discharge rate occurring at each electrode of an electrochemical capacitor are different but related to each other. The electrochemical capacitor containing hydroquinone in the electrolyte exhibits a much quicker self-discharge rate than that using a negative electrode based on grafted anthraquinone with a 50% decay of the cell voltage of the fully charged device in 0.6 and 6 h, respectively. The fast self-discharge of the former is due to the diffusion of benzoquinone molecules (formed at the positive electrode during charging) to the negative electrode, where they are reduced, causing a quick depolarization. The grafting of anthraquinone molecules on the carbon material of the negative electrode led to a much slower self-discharge, which nonetheless occurred, by the reaction of the reduced form of the grafted species with electrolyte species.

Toroidal modelling of resistive internal kink and fishbone instabilities

NASA Astrophysics Data System (ADS)

Wu, Tingting; He, Hongda; Liu, Yueqiang; Liu, Yue; Hao, G. Z.; Zhu, Jinxia

2018-05-01

The influence of energetic particles and plasma resistivity on the n =1 ( n is the toroidal mode number) internal kink and fishbone modes in tokamak plasmas is numerically investigated, using the full toroidal, resistive magnetohydrodynamic-kinetic hybrid stability code MARS-K [Liu et al., Phys. Plasmas 15 112503 (2008)]. The results show that energetic particles can either stabilize or destabilize the ideal internal kink mode, depending on the radial profiles of the particles' density and pressure. Resistive fishbones with and without an ideal wall are investigated. It is found that, in the presence of energetic particles as well as plasma resistivity, two branches of unstable roots exist, for a plasma which is ideally stable to the internal kink instability. One is the resistive internal kink mode. The other is the resistive fishbone mode. These two-branch solutions show similar behaviors, independent of whether the initial ideal kink stability is due to an ideal wall stabilization for high-beta plasmas, or due to a stable equilibrium below the Bussac pressure limit. For a realistic toroidal plasma, the resistive internal kink is the dominant instability, which grows much faster than the resistive fishbone. The plasma resistivity destabilizes the resistive internal kink while stabilizes the resistive fishbone. Systematic comparison with an analytic model qualitatively confirms the MARS-K results. Compared to analytic models based on the perturbative approach, MARS-K offers an improved physics model via self-consistent treatment of coupling between the fluid and kinetic effects due to energetic particles.

Application of kinetic flux vector splitting scheme for solving multi-dimensional hydrodynamical models of semiconductor devices

NASA Astrophysics Data System (ADS)

Nisar, Ubaid Ahmed; Ashraf, Waqas; Qamar, Shamsul

In this article, one and two-dimensional hydrodynamical models of semiconductor devices are numerically investigated. The models treat the propagation of electrons in a semiconductor device as the flow of a charged compressible fluid. It plays an important role in predicting the behavior of electron flow in semiconductor devices. Mathematically, the governing equations form a convection-diffusion type system with a right hand side describing the relaxation effects and interaction with a self consistent electric field. The proposed numerical scheme is a splitting scheme based on the kinetic flux-vector splitting (KFVS) method for the hyperbolic step, and a semi-implicit Runge-Kutta method for the relaxation step. The KFVS method is based on the direct splitting of macroscopic flux functions of the system on the cell interfaces. The second order accuracy of the scheme is achieved by using MUSCL-type initial reconstruction and Runge-Kutta time stepping method. Several case studies are considered. For validation, the results of current scheme are compared with those obtained from the splitting scheme based on the NT central scheme. The effects of various parameters such as low field mobility, device length, lattice temperature and voltage are analyzed. The accuracy, efficiency and simplicity of the proposed KFVS scheme validates its generic applicability to the given model equations. A two dimensional simulation is also performed by KFVS method for a MESFET device, producing results in good agreement with those obtained by NT-central scheme.

SUPECA kinetics for scaling redox reactions in networks of mixed substrates and consumers and an example application to aerobic soil respiration

DOE PAGES

Tang, Jin-Yun; Riley, William J.

2017-09-05

Several land biogeochemical models used for studying carbon–climate feedbacks have begun explicitly representing microbial dynamics. However, to our knowledge, there has been no theoretical work on how to achieve a consistent scaling of the complex biogeochemical reactions from microbial individuals to populations, communities, and interactions with plants and mineral soils. We focus here on developing a mathematical formulation of the substrate–consumer relationships for consumer-mediated redox reactions of the form A +  B E→  products, where products could be, e.g., microbial biomass or bioproducts. Under the quasi-steady-state approximation, these substrate–consumer relationships can be formulated as the computationally difficult full equilibrium chemistrymore » problem or approximated analytically with the dual Monod (DM) or synthesizing unit (SU) kinetics. We find that DM kinetics is scaling inconsistently for reaction networks because (1) substrate limitations are not considered, (2) contradictory assumptions are made regarding the substrate processing rate when transitioning from single- to multi-substrate redox reactions, and (3) the product generation rate cannot be scaled from one to multiple substrates. In contrast, SU kinetics consistently scales the product generation rate from one to multiple substrates but predicts unrealistic results as consumer abundances reach large values with respect to their substrates. We attribute this deficit to SU's failure to incorporate substrate limitation in its derivation. To address these issues, we propose SUPECA (SU plus the equilibrium chemistry approximation – ECA) kinetics, which consistently imposes substrate and consumer mass balance constraints. We show that SUPECA kinetics satisfies the partition principle, i.e., scaling invariance across a network of an arbitrary number of reactions (e.g., as in Newton's law of motion and Dalton's law of partial pressures). We tested SUPECA kinetics with the equilibrium chemistry solution for some simple problems and found SUPECA outperformed SU kinetics. As an example application, we show that a steady-state SUPECA-based approach predicted an aerobic soil respiration moisture response function that agreed well with laboratory observations. We conclude that, as an extension to SU and ECA kinetics, SUPECA provides a robust mathematical representation of complex soil substrate–consumer interactions and can be applied to improve Earth system model (ESM) land models.« less

SUPECA kinetics for scaling redox reactions in networks of mixed substrates and consumers and an example application to aerobic soil respiration

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tang, Jin-Yun; Riley, William J.

Several land biogeochemical models used for studying carbon–climate feedbacks have begun explicitly representing microbial dynamics. However, to our knowledge, there has been no theoretical work on how to achieve a consistent scaling of the complex biogeochemical reactions from microbial individuals to populations, communities, and interactions with plants and mineral soils. We focus here on developing a mathematical formulation of the substrate–consumer relationships for consumer-mediated redox reactions of the form A +  B E→  products, where products could be, e.g., microbial biomass or bioproducts. Under the quasi-steady-state approximation, these substrate–consumer relationships can be formulated as the computationally difficult full equilibrium chemistrymore » problem or approximated analytically with the dual Monod (DM) or synthesizing unit (SU) kinetics. We find that DM kinetics is scaling inconsistently for reaction networks because (1) substrate limitations are not considered, (2) contradictory assumptions are made regarding the substrate processing rate when transitioning from single- to multi-substrate redox reactions, and (3) the product generation rate cannot be scaled from one to multiple substrates. In contrast, SU kinetics consistently scales the product generation rate from one to multiple substrates but predicts unrealistic results as consumer abundances reach large values with respect to their substrates. We attribute this deficit to SU's failure to incorporate substrate limitation in its derivation. To address these issues, we propose SUPECA (SU plus the equilibrium chemistry approximation – ECA) kinetics, which consistently imposes substrate and consumer mass balance constraints. We show that SUPECA kinetics satisfies the partition principle, i.e., scaling invariance across a network of an arbitrary number of reactions (e.g., as in Newton's law of motion and Dalton's law of partial pressures). We tested SUPECA kinetics with the equilibrium chemistry solution for some simple problems and found SUPECA outperformed SU kinetics. As an example application, we show that a steady-state SUPECA-based approach predicted an aerobic soil respiration moisture response function that agreed well with laboratory observations. We conclude that, as an extension to SU and ECA kinetics, SUPECA provides a robust mathematical representation of complex soil substrate–consumer interactions and can be applied to improve Earth system model (ESM) land models.« less

SUPECA kinetics for scaling redox reactions in networks of mixed substrates and consumers and an example application to aerobic soil respiration

NASA Astrophysics Data System (ADS)

Tang, Jin-Yun; Riley, William J.

2017-09-01

Several land biogeochemical models used for studying carbon-climate feedbacks have begun explicitly representing microbial dynamics. However, to our knowledge, there has been no theoretical work on how to achieve a consistent scaling of the complex biogeochemical reactions from microbial individuals to populations, communities, and interactions with plants and mineral soils. We focus here on developing a mathematical formulation of the substrate-consumer relationships for consumer-mediated redox reactions of the form A + BE→ products, where products could be, e.g., microbial biomass or bioproducts. Under the quasi-steady-state approximation, these substrate-consumer relationships can be formulated as the computationally difficult full equilibrium chemistry problem or approximated analytically with the dual Monod (DM) or synthesizing unit (SU) kinetics. We find that DM kinetics is scaling inconsistently for reaction networks because (1) substrate limitations are not considered, (2) contradictory assumptions are made regarding the substrate processing rate when transitioning from single- to multi-substrate redox reactions, and (3) the product generation rate cannot be scaled from one to multiple substrates. In contrast, SU kinetics consistently scales the product generation rate from one to multiple substrates but predicts unrealistic results as consumer abundances reach large values with respect to their substrates. We attribute this deficit to SU's failure to incorporate substrate limitation in its derivation. To address these issues, we propose SUPECA (SU plus the equilibrium chemistry approximation - ECA) kinetics, which consistently imposes substrate and consumer mass balance constraints. We show that SUPECA kinetics satisfies the partition principle, i.e., scaling invariance across a network of an arbitrary number of reactions (e.g., as in Newton's law of motion and Dalton's law of partial pressures). We tested SUPECA kinetics with the equilibrium chemistry solution for some simple problems and found SUPECA outperformed SU kinetics. As an example application, we show that a steady-state SUPECA-based approach predicted an aerobic soil respiration moisture response function that agreed well with laboratory observations. We conclude that, as an extension to SU and ECA kinetics, SUPECA provides a robust mathematical representation of complex soil substrate-consumer interactions and can be applied to improve Earth system model (ESM) land models.