An energy-transport model for semiconductors derived from the Boltzmann equation

An energy-transport model is rigorously derived from the Boltzmann transport equation of semiconductors under the hypothesis that the energy gain or loss of the electrons by the phonon collisions is weak. Retaining at leading order electron-electron collisions and elastic collisions (i.e., impurity scattering and the elastic part of phonon collisions), a rigorous diffusion limit of the Boltzmann equation can be carried over, which leads to a set of diffusion equations for the electron density and temperature. The derivation is given in both the degenerate and nondegenerate cases.     

Diffusion in Lorentz lattice gas automata with backscattering

The probability of first return to the initial intervalx and the diffusion tensorD _x are calculated exactly for a ballistic Lorentz gas on a Bethe lattice or Cayley tree. It consists of a moving particle and a fixed array of scatterers, located at the nodes, and the lengths of the intervals between scatterers are determined by a geometric distribution. The same values forx andD _x apply also to a regular space lattice with a fraction of sites occupied by a scatterer in the limit of a small concentration of scatterers. If backscattering occurs, the results are very different from the Boltzmann approximation. The theory is applied to different types of lattices and different types of scatterers having rotational or mirror symmetries.     

Multi-element Saha–Boltzmann and Boltzmann plots in laser-induced plasmas

A multi-element Saha–Boltzmann plot method is proposed for the determination of the temperature and the relative number density in laser-induced plasmas, assuming local thermodynamic equilibrium and stoichiometry conservation. The method has been applied to the characterization of a plasma generated with a Cu–Fe–Ni–Mn alloy, using a Nd:YAG laser in air at atmospheric pressure. Spectra of the local emissivity have been obtained by spatial deconvolution of the intensity spectra, obtained with spatial resolution. Saha–Boltzmann plots obtained from the emissivities of 58 spectral lines of Fe I, Fe II, Ni I, Ni II, Mn I and Mn II have been fitted to linear behavior with high correlation, which shows the validity of the equation proposed. Radial distributions of the temperature and number densities of neutral atoms and ions have been determined. The results obtained reinforce the initial considerations of local thermodynamic equilibrium and conservation of stoichiometry. The proposed equation can also be applied to only one ionization species (multi-element Boltzmann plot). Spatially-integrated measurements of the plasma emission have also been performed to show that, in this case, the application of the method to the line intensities provides the two different apparent temperatures for neutral atoms and ions.     

Molecular mechanics and dynamics of biomolecules using a solvent continuum model

An easy implementation of molecular mechanics and molecular dynamics simulation using a continuum solvent model is presented that is particularly suitable for biomolecular simulations. The computation of solvation forces is made using the linear Poisson-Boltzmann equation (polar contribution) and the solvent-accessible surface area approach (nonpolar contribution). The feasibility of the methodology is demonstrated on a small protein and a small DNA hairpin. Although the parameters employed in this model must be refined to gain reliability, the performance of the method, with a standard choice of parameters, is comparable with results obtained by explicit water simulations. Copyright 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1830-1842, 2001     

Origin of the Landau-Lifshitz hydrodynamic fluctuations in nonequilibrium systems and a new method for reducing the Boltzmann equation

The Landau-Lifshitz fluctuating fluxes in fluctuating hydrodynamics are derived from the deterministic Boltzmann equation with the aid of a reduction method developed by Fujisaka and Mori. Thus it is shown that the hydrodynamic fluctuations innonequilibrium systems are generated by the reduction of variables from the-space distribution function to its five momentum moments, i.e., the hydrodynamic variables. This differs from the Bixon-Zwanzig and Fox-Uhlenbeck theories, in which the Landau-Lifshitz fluctuating fluxes are derived from the molecular fluctuating force in thestochastic Boltzmann-Langevin equation, which is, however, negligible in nonequilibrium systems. Thus the present method improves the Chapman-Enskog reduction method so as to include the hydrodynamic fluctuations generated by the reduction of variables.Supported in part by the Scientific Research Fund of the Ministry of Education.     

The stochastic Boltzmann equation and hydrodynamic fluctuations

Based on the assumption of a kinetic equation in space, a stochastic differential equation of the one-particle distribution is derived without the use of the linear approximation. It is just the Boltzmann equation with a Langevin-fluctuating force term. The result is the general form of the linearized Boltzmann equation with fluctuations found by Bixon and Zwanzig and by Fox and Uhlenbeck. It reduces to the general Landau-Lifshitz equations of fluid dynamics in the presence of fluctuations in a similar hydrodynamic approximation to that used by Chapman and Enskog with respect to the Boltzmann equation.This work received financial support from the Alexander von Humboldt Foundation.     

Action of HCl on 3-hydroxypyrazolo(iso)quinolines to give 1-chloropyrazoles: evidence for an addition–elimination mechanism by ab initio calculations in gas phase and water

 Addition–elimination reactions involving a nucleophile and a remote leaving group [S^H _N(AE)^tele] are well-known under basic conditions, especially amongst electron-poor six-membered heterocycles, but are less commonly encountered for five-membered heterocycles and are rare under acidic conditions. Concentrated HCl converts 1-hydroxy-1H-pyrazolo[3,4-c] isoquinoline and 1-hydroxy-1H-pyrazolo[3,4-c]quinoline into 3-chloro-1H-pyrazolo[3,4-c]isoquinoline and 3-chloro-1H-pyrazolo[3,4-c]quinoline, respectively. However, apparently neither the isomeric 1-hydroxy-1H-pyrazolo[4,3-c](iso)-quinolines nor the parent 1-hydroxypyrazole undergo this reaction. Additionally, all these systems are refractory under basic conditions. We present a plausible mechanism for the reaction, involving the 3-addition of Cl^- to the diprotonated heterocycle, followed by the elimination of water. Calculations of the initial transition states and intermediates, using optimisation at B3LYP/6-311+G(d,p), including thermochemistry [HF/6-31+G(d)], and single-point Poisson–Boltzmann self-consistent reaction field determination of the free energy of solvation (Jaguar Poisson–Boltzmann self-consistent reaction field), support this mechanism and reproduce the observed order of reactivity, the addition step being 2–4 kcal less favourable for the isomeric 1-hydroxy-1H-pyrazolo[4,3-c](iso)quinolines and provide a rationalisation for the role of strong acid. Received: 27 June 2002 / Accepted: 6 September 2002 / Published online: 14 February 2003     

Spatial relaxation of electrons in nonisothermal plasmas

The spatial relaxation of electrons to homogeneous states under the action of space-independent electric fields is investigated in helium, krypton, and N₂ plasmas for various electric field strengths. These investigations are based on a new method recently developed for solving the one-dimensional inhomogeneous electron Boltzmann equation in weakly ionized, collision-dominated plasmas. Elastic as well as conservative inelastic collisions of electrons with gas atoms have been included in the kinetic treatment. The spatial relaxation is caused by an imposed direct disturbance in the velocity distribution of the electrons on a spatial boundary. A pronounced dependence of the relaxation structure and the resultant relaxation length on the atomic data of the electron collision processes in different gases has been found. Furthermore the relaxation process sensitively depends on the electric field strength in the region of medium field values.     

Recent development of non-faradaic potentiometry

Non-faradaic potentiometry has been plagued by a great many fundamental errors and a lack of conceptualization. Of greatest concern is the second Nernst equation hiatus. Potentiometry may be generally classified as faradaic and non-faradaic. The former deals with the redox reactions using the Nernst equation to explain the potential origin. The latter deals with the non-redox reactions using the Boltzmann and modified Boltzmann equations to explain the origin of electrode potential. Redox faradaic potentiometry has been well described in the textbooks. However, non-faradaic potentiometry has been almost completely neglected in the literature. Many well-known electrodes, such as the pH glass electrode, common reference electrodes, and ion selective electrodes (ISE) have been mistakenly interpreted as redox reactions or ion exchange reactions. New theories and experimental results show their mechanisms to be non-faradaic in nature. Furthermore, the reaction mechanisms for ISE have been confused in textbooks with redox reactions and the Nernst equation. The ISE potentials originating from adsorption of ions or charged particles based on surface charge density will be explained using the double and counterion triple layers concept. The new counterion triple layer concept may be applied to the potential development of sensors. The reason for a new concept, theory, or mechanism is to better explain the phenomena. Examples will be given of how our new concept explains the capacitor, counterion triple layer, surface adsorbed layers interactions, and the interface structure. We will also discuss the new sensor development based on the new adsorption concept. For the first time a new type of Ag/AgCl reference electrode for non-faradaic potentiometry will be presented, one without a liquid junction and with a Pt wire instead of a salt bridge. They will help open up a new horizon for electrochemical sensor research and may be used under unusual conditions, such as high temperature and high pressure, stirring, etc.     

Semiglobal simplex optimization and its application to determining the preferred solvation sites of proteins

The classical simplex method is extended into the Semiglobal Simplex (SGS) algorithm. Although SGS does not guarantee finding the global minimum, it affords a much more thorough exploration of the local minima than any traditional minimization method. The basic idea of SGS is to perform a local minimization in each step of the simplex algorithm, and thus, similarly to the Convex Global Underestimator (CGU) method, the search is carried out on a surface spanned by local minima. The SGS and CGU methods are compared by minimizing a set of test functions of increasing complexity, each with a known global minimum and many local minima. Although CGU delivers substantially better success rates in simple problems, the two methods become comparable as the complexity of the problems increases. Because SGS is generally faster than CGU, it is the method of choice for solving optimization problems in which function evaluation is computationally inexpensive and the search region is large. The extreme simplicity of the method is also a factor. The SGS method is applied here to the problem of finding the most preferred (i.e., minimum free energy) solvation sites on a streptavidin monomer. It is shown that the SGS method locates the same lowest free energy positions as an exhaustive multistart Simplex search of the protein surface, with less than one-tenth the number of minizations. The combination of the two methods, i.e.. multistart simplex and SGS, provides a reliable procedure for predicting all potential solvation sites of a protein.