A Kinetic Vlasov Model for Plasma Simulation Using Discontinuous Galerkin Method on Many-Core Architectures

Thesis (Ph.D.)--University of Washington, 2016-03

Interpretation of the temperature dependence of rate constants in biosensor studies

A comparison is made between Arrhenius and transition-state analyses of the temperature dependence of rate constants reported in four published biosensor studies. Although the Eyring transition-state theory seemingly affords a more definitive solution to the problem of characterizing the activation energetics, the analysis is equivocal because of inherent assumptions about reaction mechanism and the magnitude of the transmission coefficient. In view of those uncertainties it is suggested that a preferable course of action entails reversion to the empirical Arrhenius analysis with regard to the energy of activation and a preexponential factor. The former is essentially equivalent to the enthalpy of activation, whereas the magnitude of the latter indicates directly the extent of disparity between the frequency of product formation and the universal frequency factor (temperature multiplied by the ratio of the Boltzmann and Planck constants) and hence the likelihood of a more complicated kinetic mechanism than that encompassed by the Eyring transition-state theory. (C) 2004 Elsevier Inc. All rights reserved.

Influence of lattice interactions on the jahn-teller distortion of the [Cu(H2O)(6)](2+) ion: Dependence of the crystal structure of K-2[Cu(H2O)(6)](SO4)(2x)(SeO4)(2-2x) upon the sulfate/selenate ratio

The temperature dependence of the structure of the mixed-anion Tutton salt K-2[Cu(H2O)(6)](SO4)(2x)(SeO4)(2-2x) has been determined for crystals with 0, 17, 25, 68, 78, and 100% sulfate over the temperature range of 85-320 K. In every case, the [Cu(H2O)(6)](2+) ion adopts a tetragonally elongated coordination geometry with an orthorhombic distortion. However, for the compounds with 0, 17, and 25% sulfate, the long and intermediate bonds occur on a different pair of water molecules from those with 68, 78, and 100% sulfate. A thermal equilibrium between the two forms is observed for each crystal, with this developing more readily as the proportions of the two counterions become more similar. Attempts to prepare a crystal with approximately equal amounts of sulfate and selenate were unsuccessful. The temperature dependence of the bond lengths has been analyzed using a model in which the Jahn-Teller potential surface of the [Cu(H2O)(6)](2+) ion is perturbed by a lattice-strain interaction. The magnitude and sign of the orthorhombic component of this strain interaction depends on the proportion of sulfate to selenate. Significant deviations from Boltzmann statistics are observed for those crystals exhibiting a large temperature dependence of the average bond lengths, and this may be explained by cooperative interactions between neighboring complexes.

Noisy fitness evaluation in genetic algorithms and the dynamics of learning

A theoretical model is presented which describes selection in a genetic algorithm (GA) under a stochastic fitness measure and correctly accounts for finite population effects. Although this model describes a number of selection schemes, we only consider Boltzmann selection in detail here as results for this form of selection are particularly transparent when fitness is corrupted by additive Gaussian noise. Finite population effects are shown to be of fundamental importance in this case, as the noise has no effect in the infinite population limit. In the limit of weak selection we show how the effects of any Gaussian noise can be removed by increasing the population size appropriately. The theory is tested on two closely related problems: the one-max problem corrupted by Gaussian noise and generalization in a perceptron with binary weights. The averaged dynamics can be accurately modelled for both problems using a formalism which describes the dynamics of the GA using methods from statistical mechanics. The second problem is a simple example of a learning problem and by considering this problem we show how the accurate characterization of noise in the fitness evaluation may be relevant in machine learning. The training error (negative fitness) is the number of misclassified training examples in a batch and can be considered as a noisy version of the generalization error if an independent batch is used for each evaluation. The noise is due to the finite batch size and in the limit of large problem size and weak selection we show how the effect of this noise can be removed by increasing the population size. This allows the optimal batch size to be determined, which minimizes computation time as well as the total number of training examples required.

The charge effect on the hindrance factors for diffusion and convection of a solute in pores II

The diffusion and convection of a solute suspended in a fluid across porous membranes are known to be reduced compared to those in a bulk solution, owing to the fluid mechanical interaction between the solute and the pore wall as well as steric restriction. If the solute and the pore wall are electrically charged, the electrostatic interaction between them could affect the hindrance to diffusion and convection. In this study, the transport of charged spherical solutes through charged circular cylindrical pores filled with an electrolyte solution containing small ions was studied numerically by using a fluid mechanical and electrostatic model. Based on a mean field theory, the electrostatic interaction energy between the solute and the pore wall was estimated from the Poisson-Boltzmann equation, and the charge effect on the solute transport was examined for the solute and pore wall of like charge. The results were compared with those obtained from the linearized form of the Poisson-Boltzmann equation, i.e.the Debye-Hückel equation. © 2012 The Japan Society of Fluid Mechanics and IOP Publishing Ltd.

Fluid mechanical and electrostatic study on the osmotic flow through cylindrical pores

When two solutions differing in solute concentration are separated by a porous membrane, the osmotic pressure will generate a net volume flux of the suspending fluid across the membrane; this is termed osmotic flow. We consider the osmotic flow across a membrane with circular cylindrical pores when the solute and the pore walls are electrically charged, and the suspending fluid is an electrolytic solution containing small cations and anions. Under the condition in which the radius of the pores and that of the solute molecules greatly exceed those of the solvent as well as the ions, a fluid mechanical and electrostatic theory is introduced to describe the osmotic flow in the presence of electric charge. The interaction energy, including the electrostatic interaction between the solute and the pore wall, plays a key role in determining the osmotic flow. We examine the electrostatic effect on the osmotic flow and discuss the difference in the interaction energy determined from the nonlinear Poisson-Boltzmann equation and from its linearized equation (the Debye-Hückel equation).

Electrical charge effect on osmotic flow through pores

An electrostatic model for osmotic flow through circular cylindrical pores is developed to describe the reflection coefficient for the membrane transport in the presence of surface charges on the pore wall and the solute. For a spherical solute placed at an arbitrary radial position in the pore, the electrical potential was computed by a spectral element method applied to the Poisson-Boltzmann equation together with the condition of electrical neutrality. The interaction energy between the surface charges was used to estimate the osmotic reflection coefficient. The proposed model predicts that even for a small Debye length compared to the pore radius, the repulsive electrostatic interaction between the surface charges could significantly increase the osmotic flow through the pore.

Assessment of physical and hydrological properties of biological soil crusts using x-ray microtomography and modelling

Biological soil crusts (BSCs) are formed by aggregates of soil particles and communities of microbial organisms and are common in all drylands. The role of BSCs on infiltration remains uncertain due to the lack of data on their role in affecting soil physical properties such as porosity and structure. Quantitative assessment of these properties is primarily hindered by the fragile nature of the crusts. Here we show how the use of a combination of non-destructive imaging X-ray microtomography (XMT) and Lattice Boltzmann method (LBM) enables quantification of key soil physical parameters and the modeling of water flow through BSCs samples from Kalahari Sands, Botswana. We quantify porosity and flow changes as a result of mechanical disturbance of such a fragile cyanobacteria-dominated crust. Results show significant variations in porosity between different types of crusts and how they affect the flow and that disturbance of a cyanobacteria-dominated crust results in the breakdown of larger pore spaces and reduces flow rates through the surface layer. We conclude that the XMT–LBM approach is well suited for study of fragile surface crust samples where physical and hydraulic properties cannot be easily quantified using conventional methods.

High-frequency conductivity of optically excited charge carriers in hydrogenated nanocrystalline silicon investigated by spectroscopic femtosecond pump-probe reflectivity measurements

We report an investigation into the high-frequency conductivity of optically excited charge carriers far from equilibrium with the lattice. The investigated samples consist of hydrogenated nanocrystalline silicon films grown on a thin film of silicon oxide on top of a silicon substrate. For the investigation, we used an optical femtosecond pump-probe setup to measure the reflectance change of a probe beam. The pump beam ranged between 580 and 820nm, whereas the probe wavelength spanned 770 to 810nm. The pump fluence was fixed at 0.6mJ/cm2. We show that at a fixed delay time of 300fs, the conductivity of the excited electron-hole plasma is described well by a classical conductivity model of a hot charge carrier gas found at Maxwell-Boltzmann distribution, while Fermi-Dirac statics is not suitable. This is corroborated by values retrieved from pump-probe reflectance measurements of the conductivity and its dependence on the excitation wavelength and carrier temperature. The conductivity decreases monotonically as a function of the excitation wavelength, as expected for a nondegenerate charge carrier gas.

Fractional Calculus of P-transforms

MSC 2010: 44A20, 33C60, 44A10, 26A33, 33C20, 85A99

Numerical investigation on the heat transfer enhancement using micro/nano phase-change particulate flow

The introduction of phase change material fluid and nanofluid in micro-channel heat sink design can significantly increase the cooling capacity of the heat sink because of the unique features of these two kinds of fluids. To better assist the design of a high performance micro-channel heat sink using phase change fluid and nanofluid, the heat transfer enhancement mechanism behind the flow with such fluids must be completely understood. ^ A detailed parametric study is conducted to further investigate the heat transfer enhancement of the phase change material particle suspension flow, by using the two-phase non-thermal-equilibrium model developed by Hao and Tao (2004). The parametric study is conducted under normal conditions with Reynolds numbers of Re = 90–600 and phase change material particle concentrations of ϵp ≤ 0.25, as well as extreme conditions of very low Reynolds numbers (Re < 50) and high phase change material particle concentration (ϵp = 50%–70%) slurry flow. By using the two newly-defined parameters, named effectiveness factor ϵeff and performance index PI, respectively, it is found that there exists an optimal relation between the channel design parameters L and D, particle volume fraction ϵp, Reynolds number Re, and the wall heat flux qw. The influence of the particle volume fraction ϵp, particle size dp, and the particle viscosity μ p, to the phase change material suspension flow, are investigated and discussed. The model was validated by available experimental data. The conclusions will assist designers in making their decisions that relate to the design or selection of a micro-pump suitable for micro or mini scale heat transfer devices. ^ To understand the heat transfer enhancement mechanism of the nanofluid flow from the particle level, the lattice Boltzmann method is used because of its mesoscopic feature and its many numerical advantages. By using a two-component lattice Boltzmann model, the heat transfer enhancement of the nanofluid is analyzed, through incorporating the different forces acting on the nanoparticles to the two-component lattice Boltzmann model. It is found that the nanofluid has better heat transfer enhancement at low Reynolds numbers, and the Brownian motion effect of the nanoparticles will be weakened by the increase of flow speed. ^

Quantifying protein folding pathways

A deep understanding of the proteins folding dynamics can be get quantifying folding landscape by calculating how the number of microscopic configurations (entropy) varies with the energy of the chain, Ω=Ω(E). Because of the incredibly large number of microstates available to a protein, direct enumeration of Ω(E) is not possible on realistic computer simulations. An estimate of Ω(E) can be obtained by use of a combination of statistical mechanics and thermodynamics. By combining different definitions of entropy that are valid for a system whose probability for occupying a state is given by the canonical Boltzmann probability, computers allow the determination of Ω(E). ^ The energy landscapes of two similar, but not identical model proteins were studied. One protein contains no kinetic tracks. Results show a smooth funnel for the folding landscape. That allows the contour determination of the folding funnel. Also it was presented results for the folding landscape for a modified protein with kinetic traps. Final results show that the computational approach is able to distinguish and explore regions of the folding landscape that are due to kinetic traps from the native state folding funnel.^

Nonlinear quantum transport in low-dimensional electronic devices

The study of transport processes in low-dimensional semiconductors requires a rigorous quantum mechanical treatment. However, a full-fledged quantum transport theory of electrons (or holes) in semiconductors of small scale, applicable in the presence of external fields of arbitrary strength, is still not available. In the literature, different approaches have been proposed, including: (a) the semiclassical Boltzmann equation, (b) perturbation theory based on Keldysh's Green functions, and (c) the Quantum Boltzmann Equation (QBE), previously derived by Van Vliet and coworkers, applicable in the realm of Kubo's Linear Response Theory (LRT). ^ In the present work, we follow the method originally proposed by Van Wet in LRT. The Hamiltonian in this approach is of the form: H = H 0(E, B) + λV, where H0 contains the externally applied fields, and λV includes many-body interactions. This Hamiltonian differs from the LRT Hamiltonian, H = H0 - AF(t) + λV, which contains the external field in the field-response part, -AF(t). For the nonlinear problem, the eigenfunctions of the system Hamiltonian, H0(E, B), include the external fields without any limitation on strength. ^ In Part A of this dissertation, both the diagonal and nondiagonal Master equations are obtained after applying projection operators to the von Neumann equation for the density operator in the interaction picture, and taking the Van Hove limit, (λ → 0, t → ∞, so that (λ2 t)n remains finite). Similarly, the many-body current operator J is obtained from the Heisenberg equation of motion. ^ In Part B, the Quantum Boltzmann Equation is obtained in the occupation-number representation for an electron gas, interacting with phonons or impurities. On the one-body level, the current operator obtained in Part A leads to the Generalized Calecki current for electric and magnetic fields of arbitrary strength. Furthermore, in this part, the LRT results for the current and conductance are recovered in the limit of small electric fields. ^ In Part C, we apply the above results to the study of both linear and nonlinear longitudinal magneto-conductance in quasi one-dimensional quantum wires (1D QW). We have thus been able to quantitatively explain the experimental results, recently published by C. Brick, et al., on these novel frontier-type devices. ^

Nonlinear quantum transport in low-dimensional electronic devices

The study of transport processes in low-dimensional semiconductors requires a rigorous quantum mechanical treatment. However, a full-fledged quantum transport theory of electrons (or holes) in semiconductors of small scale, applicable in the presence of external fields of arbitrary strength, is still not available. In the literature, different approaches have been proposed, including: (a) the semiclassical Boltzmann equation, (b) perturbation theory based on Keldysh's Green functions, and (c) the Quantum Boltzmann Equation (QBE), previously derived by Van Vliet and coworkers, applicable in the realm of Kubo's Linear Response Theory (LRT). In the present work, we follow the method originally proposed by Van Vliet in LRT. The Hamiltonian in this approach is of the form: H = H°(E, B) + λV, where H0 contains the externally applied fields, and λV includes many-body interactions. This Hamiltonian differs from the LRT Hamiltonian, H = H° - AF(t) + λV, which contains the external field in the field-response part, -AF(t). For the nonlinear problem, the eigenfunctions of the system Hamiltonian, H°(E, B) , include the external fields without any limitation on strength. In Part A of this dissertation, both the diagonal and nondiagonal Master equations are obtained after applying projection operators to the von Neumann equation for the density operator in the interaction picture, and taking the Van Hove limit, (λ → 0 , t → ∞ , so that (λ2 t)n remains finite). Similarly, the many-body current operator J is obtained from the Heisenberg equation of motion. In Part B, the Quantum Boltzmann Equation is obtained in the occupation-number representation for an electron gas, interacting with phonons or impurities. On the one-body level, the current operator obtained in Part A leads to the Generalized Calecki current for electric and magnetic fields of arbitrary strength. Furthermore, in this part, the LRT results for the current and conductance are recovered in the limit of small electric fields. In Part C, we apply the above results to the study of both linear and nonlinear longitudinal magneto-conductance in quasi one-dimensional quantum wires (1D QW). We have thus been able to quantitatively explain the experimental results, recently published by C. Brick, et al., on these novel frontier-type devices.

Física do meio ambiente como disciplina nas licenciaturas em física: proposta de uma unidade didática para ensinar efeito estufa e aquecimento global

In this dissertation we propose a Teaching Unit of Physics to teach content through environmental discussions of the greenhouse effect and global warming. This teaching unit is based on a problem-methodological intervention from the application of the method of the Arch of Charles Maguerez. The methodological foundations of the thesis are embedded in action research and this is structured in five chapters: the first chapter deals with the Physical Environment (FMA) as a subject in Degree Courses in Physics in Brazil, bringing the concern of how this discipline has been taught. We started the first chapter explaining the reasons behind the inclusion of the discipline of Physical Environment in a Physics Degree Courses. Then we did a search on the websites of Institutions of Higher Education, to know of the existence or not of this discipline on curricular. We then analyzed the menus to see what bibliographies are being adopted and what content of Physics are being worked, and how it has been done. The courses surveyed were those of Federal and Federal Institutes Universities. Thus ended the first chapter. Given the inseparability between studies in Physics Teaching and studies on competencies, skills and significant learning, wrote the second chapter. In this chapter we discuss the challenge of converting information into knowledge. Initially on initial teacher training, because even if this is not our focus, the study is a discipline on the upper reaches, therefore, offered to future teachers. Then we talked about the culture of knowledge, where we emphasize the use of a teaching approach that promotes meanings taught by content and make sense to the student. We finished the third chapter, making some considerations on skills and abilities, in order to identify what skills and competencies were developed and worked during and after the implementation of Curriculum Unit. The third chapter is the result of a literature review and study of the radioactive EarthSun interaction. The subjects researched approach from the generation of energy in the sun to topics stain solar coronal mass ejections, solar wind, black body radiation, Wien displacement law, Stefan-Boltzmann Law, greenhouse effect and global warming. This chapter deals with material support for the teacher of the aforementioned discipline. The fourth chapter talks about the arc method of Charles Maguerez; Here we explain the structure of each of the five steps of the Arc and how to use them in teaching. We also show another version of this method adapted by Bordenave. In the fifth and final chapter brought a description of how the method of Arc was used in physics classes of Environment, with students majoring in Physics IFRN Campus Santa Cruz. Here, in this chapter, a transcript of classes to show how was the application of a problem-based methodology in the teaching of content proposed for Physics Teaching Unit from the environmental discussion about the greenhouse effect and global warming phenomena