Speed of reaction-transport processes

We present an approach to determining the speed of wave-front solutions to reaction-transport processes. This method is more accurate than previous ones. This is explicitly shown for several cases of practical interest: (i) the anomalous diffusion reaction, (ii) reaction diffusion in an advective field, and (iii) time-delayed reaction diffusion. There is good agreement with the results of numerical simulations

Femtosecond dynamics of molecules and clusters

The real-time dynamics of multiphoton ionization and fragmentation of molecules - Na_2 , Na_3 - and clusters - Na_n, Hg_n - has been studied in molecular beam experiments employing ion and electron spectroscopy together with femtosecond pump-probe techniques. Experiments with Na_2 and Na_3 reveal unexpected features of the dynamics of the absorption of several photons as seen in the one- and three dimensional vibrational wave packet motion in different potential surfaces and in high laser fields. Cluster size dependent studies of physical properties such as absorption resonances, lifetimes and decay channels have been performed using tunable femtosecond light pulses in resonance enhanced multiphoton ionization (REMPI) of the cluster size under investigation. This method failed in ns-laser experiments due to the ultrafast decay of the studied cluster. For Na_n, cluster we find that for cluster sizes n \le 21 molecular excitations and properties prevail over collective excitations of plasmon-like resonances. In the case of Hg_n cluster prompt formation of singly and doubly charged cluster are observed up to n \approx 60. The transient multiphoton ionization spectra show a 'short' time wave packet dynamics, which is identical for singly and doubly charged mercury clusters while the 'long' time fragmentation dynamics is different.

Structure, Dynamics and Reactivity in the Organic Solid State: Anthracene Derivatives and Charge Transfer Crystals

The work presented in this thesis tackles some important points concerning the collective properties of two typical categories of molecular crystals, i.e., anthracene derivatives and charge transfer crystals. Anthracene derivatives have constituted the class of materials from which systematical investigations of crystal-to-crystal photodimerization reactions started, developed and have been the subject of a new awakening in the recent years. In this work some of these compounds, namely, 9-cyanoanthacene, 9-anthacenecarboxylic acid and 9-methylanthracene, have been selected as model systems for a phenomenological approach to some key properties of the solid state, investigated by spectroscopic methods. The present results show that, on the basis of the solid state organization and the chemical nature of each compound, photo-reaction dynamics and kinetics display distinctive behaviors, which allows for a classification of the various processes in topochemical, non topochemical, reversible or topophysical. The second part of the thesis was focused on charge transfer crystals, binary systems formed by stoichiometric combinations of the charge donating perylene (D) and the charge accepting tetracyano-quinodimethane (A), this latter also in its fluorinated derivatives. The work was focused on the growth of single crystals, some of which not yet reported in the literature, by PVT technique. Structural and spectroscopic characterizations have been performed, with the aim of determining the degree of charge transfer between donor and acceptor in the co-crystals. An interesting outcome of the systematic search performed in this work is the definition of the experimental conditions which drive the crystal growth of the binary systems either towards the low (1:1) or the high ratio (3:1 or 3:2) stoichiometries.

Maximizing Product Formation and Minimizing Degradation: A Look at Buffering Conditions and Post-Reaction Degradation for the In-Line Jaffe Reaction with Capillary Electrophoresis

Creatinine levels in blood serum are typically used to assess renal function. Clinical determination of creatinine is often based on the Jaffe reaction, in which creatinine in the serum reacts with sodium picrate, resulting in a spectrophotometrically quantifiable product. Previous work from our lab has introduced an electrophoretically mediated initiation of this reaction, in which nanoliter plugs of individual reagent solutions can be added to the capillary and then mixed and reacted. Following electrophoretic separation of the product from excess reactant(s), the product can be directly determined on column. This work aims to gain a detailed understanding of the in-capillary reagent mixing dynamics, in-line reaction yield, and product degradation during electrophoresis, with an overall goal of improving assay sensitivity. One set of experiments focuses on maximizing product formation through manipulation of various conditions such as pH, voltage applied, and timing of the applied voltage, in addition to manipulations in the identity, concentration, and pH of the background electrolyte. Through this work, it was determined that dramatic changes in local voltage fields within the various reagent zones lead to ineffective reagent overlapping. Use of the software simulation program Simul 5 enabled visualization of the reaction dynamics within the capillary, specifically the wide variance between the electric field intensities within the creatinine and picrate zones. Because of this simulation work, the experimental method was modified to increase the ionic strength of the creatinine reagent zone to lower the local voltage field, thus producing more predictable and effective overlap conditions for the reagents and allowing the formation of more Jaffe product. As second set of experiments focuses on controlling the post-reaction product degradation. In that vein, we have systematically explored the importance of the identity, concentration, and pH of the background electrolyte on the post-reaction degradation rate of the product. Although prior work with borate background electrolytes indicated that product degradation was probably a function of the ionic strength of the background electrolyte, this work with a glycine background electrolyte demonstrates that degradation is in fact not a function of ionic strength of the background electrolyte. As the concentration and pH of the glycine background increased, the rate of degradation of product did not change dramatically, whereas in borate-buffered systems, the rate of Jaffe product degradation increased linearly with background electrolyte concentration above 100.0 mM borate. Similarly, increasing pH of the glycine background electrolyte did not result in a corresponding increase in product degradation, as it had with the borate background electrolyte. Other general trends that were observed include: increasing background electrolyte concentration increases peak efficiency and higher pH favors product formation; thus, it appears that use of a background electrolyte other than borate, such as glycine, the rate of degradation of the Jaffe product can be slowed, increasing the sensitivity of this in-line assay.

A subgrid viscosity Lagrance-Galerkin method for convection-diffusion problems

We present and analyze a subgrid viscosity Lagrange-Galerk in method that combines the subgrid eddy viscosity method proposed in W. Layton, A connection between subgrid scale eddy viscosity and mixed methods. Appl. Math. Comp., 133: 14 7-157, 2002, and a conventional Lagrange-Galerkin method in the framework of P1⊕ cubic bubble finite elements. This results in an efficient and easy to implement stabilized method for convection dominated convection diffusion reaction problems. Numerical experiments support the numerical analysis results and show that the new method is more accurate than the conventional Lagrange-Galerkin one.

Diffusion of linear paraffins in nanoporous silica

The diffusion of hexane, heptane, octane, and decane in nanoporous MCM-41 silica at various temperatures is investigated by the zero-length-column method. The diffusion coefficients are derived by a complete-time-range analysis of desorption curves at different purge flow rates and temperatures. The results show that the calculated low-coverage diffusivity values decrease monotonically, and the derived Henry's law constants increase, as the carbon number of paraffins increases. The study reveals that transport is strongly influenced by intracrystalline diffusion and dominated by the sorbate-sorbent interaction. The diffusion activation energy and adsorption isosteric heat at zero loading increase monotonically with the carbon number of linear paraffins, but their ratio is essentially constant for each adsorbate compound.

Diffusion of n-decane in mesoporous MCM-41 silicas

We investigate here the diffusion of n-decane in nanoporous MCM-41 silicas with pore diameters between 3.0 and 4.3 nm, and at various temperatures and purge flow rates, by the Zero Length Column method. A complete-time-range analysis of desorption curves is used to derive the diffusion coefficient, and the effect of pore size, purge flow rate and temperature on the diffusion character is systematically studied. The results show that the calculated low-coverage diffusivity values are strongly dependent on temperature but only weakly dependent on pore size. The study reveals that transport is controlled by intracrystalline diffusion and dominated by sorbate-sorbent interaction, with the experimental isosteric heat matching the potential energy of flat-lying n-decane molecules on the surface, determined using a united atom model. The diffusion activation energy and adsorption isosteric heat at zero loading for the different pore size MCM-41 samples vary in a narrow range respectively, and their ratio is essentially constant over the pore size range studied. The study shows that the ZLC method is an effective tool to investigate the diffusion kinetics of hydrocarbons in mesoporous MCM-41 materials. (c) 2005 Elsevier Inc. All rights reserved.

Reactive clusters on a membrane

We investigate the reaction dynamics of diffusive molecules with immobile binding partners. The fixed reactants build clusters that comprise just a few tens of molecules, which leads to small cluster sizes. These molecules participate in the reaction only if they are activated. The dynamics of activation is mapped to a time-dependent size of an active region within the cluster. We focus on the deterministic description of the dynamics of a single cluster. The spatial setup accounts for one of the most important determinants of the dynamics of a cluster, i.e. diffusional transport of reaction partners toward or away from the active region of the cluster. We provide numerical and analytical evidence that diffusion influences decisively the dynamic regimes of the reactions. The application of our methods to intracellular Ca²⁺ dynamics shows that large local concentrations saturate the Ca²⁺ feedback to the channel state control. That eliminates oscillations depending on this feedback.

Optimization of Heat Treatment Profiles Applied to Nanometric-Scale Nb(3)Sn Wires With Cu-Sn Artificial Pinning Centers

Nb(3)Sn is one of the most used superconducting materials for applications in high magnetic fields. The improvement of the critical current densities (J(c)) is important, and must be analyzed together with the optimization of the flux pinning acting in the material. For Nb(3)Sn, it is known that the grain boundaries are the most effective pinning centers. However, the introduction of artificial pinning centers (APCs) with different superconducting properties has been proved to be beneficial for J(c). As these APCs are normally in the nanometric-scale, the conventional heat treatment profiles used for Nb(3)Sn wires cannot be directly applied, leading to excessive grain growth and/or increase of the APCs cross sections. In this work, the heat treatment profiles for Nb(3)Sn superconductor wires with Cu(Sn) artificial pinning centers in nanometric-scale were analyzed in an attempt to improve J(c) . It is described a methodology to optimize the heat treatment profiles in respect to diffusion, reaction and formation of the superconducting phases. Microstructural, transport and magnetic characterization were performed in an attempt to find the pinning mechanisms acting in the samples. It was concluded that the maximum current densities were found when normal phases (due to the introduction of the APCs) are acting as main pinning centers in the global behavior of the Nb(3)Sn superconducting wire.

Heat treatment influence on the superconducting properties of nanometric-scale Nb(3)Sn wires with Cu-Sn artificial pinning centers

Since the discovery of Nb(3)Sn superconductors many efforts have been expended to improve the transport properties in these materials. In this work, the heat treatment profiles for Nb(3)Sn superconductor wires with Cu(Sn) artificial pinning centers (APCs) with nanometric-scale sizes were analyzed in an attempt to improve the critical current densities and upper critical magnetic field. The methodology to optimize the heat treatment profiles in respect to the diffusion, reaction and formation of the superconducting phases is described. Microstructural characterization, transport and magnetic measurements were performed in an attempt to relate the microstructure to the pinning mechanisms acting in the samples. It was concluded that the maximum current densities occur due to normal phases (APCs) that act as the main pinning centers in the global behavior of the Nb(3)Sn superconducting wire. The APC technique was shown to be very powerful because it permitted mixing of the pinning mechanism. This achievement was not possible in other studies in Nb(3)Sn wires reported up to now.

Solution of a heterogeneous modeling of a horizontal-flow anaerobic immobilized biomass (HAIB) reactor by the sequencing method

A modeling study was completed to develop a methodology that combines the sequencing and finite difference methods for the simulation of a heterogeneous model of a tubular reactor applied in the treatment of wastewater. The system included a liquid phase (convection diffusion transport) and a solid phase (diffusion reaction) that was obtained by completing a mass balance in the reactor and in the particle, respectively. The model was solved using a pilot-scale horizontal-flow anaerobic immobilized biomass (HAIB) reactor to treat domestic sewage, with the concentration results compared with the experimental data. A comparison of the behavior of the liquid phase concentration profile and the experimental results indicated that both the numerical methods offer a good description of the behavior of the concentration along the reactor. The advantage of the sequencing method over the finite difference method is that it is easier to apply and requires less computational time to model the dynamic simulation of outlet response of HAIB.

Useful representations of series solutions for transport problems

Simple techniques are presented for rearrangement of an infinite series in a systematic way such that the convergence of the resulting expression is accelerated. These procedures also allow calculation of required boundary derivatives. Several examples of conduction and diffusion-reaction problems illustrate the methods.

Lanczos subspace time-independent wave packet calculations of S(D-1)+H-2 reactive scattering

In this paper. we present the results of quantum dynamical simulations of the S (D-1) + H-2 insertion reaction on a newly developed potential energy surface (J. Chem. Phys. 2001, 114, 320). State-to-state reaction probabilities. product state distributions, and initial-state resolved cumulative reaction probabilities from a given incoming reactant channel are obtained from a time-independent wave packet analysis, performed within a single Lanczos subspace. Integral reaction cross sections are then estimated by J-shifting method and compared with the results from molecular beam experiment and QCT calculations.

Calculation of bound and resonance states of HO2 for nonzero total angular momentum

Bound and resonance states of HO2 have been calculated quantum mechanically by the Lanczos homogeneous filter diagonalization method [Zhang and Smith, Phys. Chem. Chem. Phys. 3, 2282 (2001); J. Chem. Phys. 115, 5751 (2001)] for nonzero total angular momentum J = 1,2,3. For lower bound states, agreement between the results in this paper and previous work is quite satisfactory; while for high lying bound states and resonances these are the first reported results. A helicity quantum number V assignment (within the helicity conserving approximation) is performed and the results indicate that for lower bound states it is possible to assign the V quantum numbers unambiguously, but for resonances it is impossible to assign the V helicity quantum numbers due to strong mixing. In fact, for the high-lying bound states, the mixing has already appeared. These results indicate that the helicity conserving approximation is not good for the resonance state calculations and exact quantum calculations are needed to accurately describe the reaction dynamics for HO2 system. Analysis of the resonance widths shows that most of the resonances are overlapping and the interferences between them lead to large fluctuations from one resonance to another. In accord with the conclusions from earlier J = 0 calculations, this indicates that the dissociation of HO2 is essentially irregular. (C) 2003 American Institute of Physics.

Structure-function relationships in a glycosyltransferase, a phosphatase and an oxidoreductase

Dissertation presented to obtain the Ph.D degree in Biochemistry