Sparse space-time boundary element methods for the heat equation

The goal of this work is the efficient solution of the heat equation with Dirichlet or Neumann boundary conditions using the Boundary Elements Method (BEM). Efficiently solving the heat equation is useful, as it is a simple model problem for other types of parabolic problems. In complicated spatial domains as often found in engineering, BEM can be beneficial since only the boundary of the domain has to be discretised. This makes BEM easier than domain methods such as finite elements and finite differences, conventionally combined with time-stepping schemes to solve this problem. The contribution of this work is to further decrease the complexity of solving the heat equation, leading both to speed gains (in CPU time) as well as requiring smaller amounts of memory to solve the same problem. To do this we will combine the complexity gains of boundary reduction by integral equation formulations with a discretisation using wavelet bases. This reduces the total work to O(h

Heat Conduction Equation Solution in the Presence of a Change of State in a Bounded Axisymmetric Cylindrical Domain

The heat conduction problem, in the presence of a change of state, was solved for the case of an indefinitely long cylindrical layer cavity. As boundary conditions, it is imposed that the internal surface of the cavity is maintained below the fusion temperature of the infilling substance and the external surface is kept above it. The solution, obtained in nondimensional variables, consists in two closed form heat conduction equation solutions for the solidified and liquid regions, which formally depend of the, at first, unknown position of the phase change front. The energy balance through the phase change front furnishes the equation for time dependence of the front position, which is numerically solved. Substitution of the front position for a particular instant in the heat conduction equation solutions gives the temperature distribution inside the cavity at that moment. The solution is illustrated with numerical examples. [DOI: 10.1115/1.4003542]

Determination of Very Slow Diffusion of Paramagnetic Ions by NMR Spectroscopy

In this paper, we propose a new method of measuring the very slow paramagnetic ion diffusion coefficient using a commercial high-resolution spectrometer. If there are distinct paramagnetic ions influencing the hydrogen nuclear magnetic relaxation time differently, their diffusion coefficients can be measured separately. A cylindrical phantom filled with Fricke xylenol gel solution and irradiated with gamma rays was used to validate the method. The Fricke xylenol gel solution was prepared with 270 Bloom porcine gelatin, the phantom was irradiated with gamma rays originated from a (60)Co source and a high-resolution 200 MHz nuclear magnetic resonance (NMR) spectrometer was used to obtain the phantom (1)H profile in the presence of a linear magnetic field gradient. By observing the temporal evolution of the phantom NMR profile, an apparent ferric ion diffusion coefficient of 0.50 mu m(2)/ms due to ferric ions diffusion was obtained. In any medical process where the ionizing radiation is used, the dose planning and the dose delivery are the key elements for the patient safety and success of treatment. These points become even more important in modern conformal radio therapy techniques, such as stereotactic radiosurgery, where the delivered dose in a single session of treatment can be an order of magnitude higher than the regular doses of radiotherapy. Several methods have been proposed to obtain the three-dimensional (3-D) dose distribution. Recently, we proposed an alternative method for the 3-D radiation dose mapping, where the ionizing radiation modifies the local relative concentration of Fe(2+)/Fe(3+) in a phantom containing Fricke gel and this variation is associated to the MR image intensity. The smearing of the intensity gradient is proportional to the diffusion coefficient of the Fe(3+) and Fe(2+) in the phantom. There are several methods for measurement of the ionic diffusion using NMR, however, they are applicable when the diffusion is not very slow.

Structural Characterization of Rare-Earth Doped Yttrium Aluminoborate Laser Glasses Using Solid State NMR

The structure of laser glasses in the system (B(2)O(3))(0.6){(Al(2)O(3))(0.4-x)(Y(2)O(3))(x)} (0.1 <= x <= 0.25) has been investigated by means of (11)B, (27)Al, and (89)Y solid state NMR as well as Y-3d core-level X-ray photoelectron spectroscopy, (11)B magic-angle spinning (MAS) NMR spectra reveal that the majority of the boron atoms are three-coordinated, and a slight increase of four-coordinated boron content with increasing x can be noticed. (27)Al MAS NMR spectra show that the alumina species are present in the coordination states four, five and six. All of them are in intimate contact with both the three- and the four-coordinate boron species and vice versa, as indicated by (11)B/(27)Al rotational echo double resonance (REDOR) data. These results are consistent with the formation of a homogeneous, nonsegregated glass structure. For the first time, (89)Y solid state NMR has been used to probe the local environment of Y(3+) ions in a glass-forming system. The intrinsic sensitivity problem associated with (89)Y NMR has been overcome by combining the benefits of paramagnetic doping with those of signal accumulation via Carr-Purcell spin echo trains. Both the (89)Y chemical shifts and the Y-3d core level binding energies are found to be rather sensitive to the yttrium bonding state and reveal that the bonding properties of the yttrium atoms in these glasses are similar to those found in the model compounds YBO(3) and YAl(3)(BO(3))(4), Based on charge balance considerations as well as (11)B NMR line shape analyses, the dominant borate species are concluded to be meta- and pyroborate anions.

A study of the relaxation dynamics in a quadrupolar NMR system using Quantum State Tomography

This article reports a relaxation study in an oriented system containing spin 3/2 nuclei using quantum state tomography (QST). The use of QST allowed evaluating the time evolution of all density matrix elements starting from several initial states. Using an appropriated treatment based on the Redfield theory, the relaxation rate of each density matrix element was measured and the reduced spectral densities that describe the system relaxation were determined. All the experimental data could be well described assuming pure quadrupolar relaxation and reduced spectral densities corresponding to a superposition of slow and fast motions. The data were also analyzed in the context of Quantum Information Processing, where the coherence loss of each qubit of the system was determined using the partial trace operation. (C) 2008 Elsevier Inc. All rights reserved.

Ruthenium(II) complexes containing N(4)-tolyl-2-acetylpyridine thiosemicarbazones and phosphine ligands: NMR and electrochemical studies of cis-trans isomerization

[Ru(HL)(PPh3)(2)Cl]Cl complexes have been obtained in which HL = N(4)-ortho (complex 1), N(4)-meta (complex 2) and N(4) pctratolyl 2-acetylpyridine thiosemicarbazone (complex 3). NMR and electrochemical studies indicate that both cis and trans isomers exist in solution, and that the cis isomers are converted into the trans isomers with time. Crystal structure determination of (1) reveals that the traps isomer is formed in the solid state. (c) 2007 Elsevier B.V. All rights reserved.

CONTRA: Improving the performance of dynamic investigations in natural abundance organic solids by mirror-symmetric constant-time CODEX

We present a minor but essential modification to the CODEX 1D-MAS exchange experiment. The new CONTRA method, which requires minor changes of the original sequence only, has advantages over the previously introduced S-CODEX, since it is less sensitive to artefacts caused by finite pulse lengths. The performance of this variant, including the finite pulse effect, was confirmed by SIMPSON calculations and demonstrated on a number of dynamic systems. (C) 2007 Elsevier Inc. All rights reserved.

NMR study of slow motions of HDA hydrocarbons chains inside lamellar structures

Measurements of H-1 and C-13 Nuclear Magnetic Resonance (NMR) for the nano-composite materials formed by the intercalation of hexadecylamine (HDA) in metal oxides (TiO2, V2O5 and MoO3), are reported. The H-1 NMR spin-lattice relaxation in the rotating frame was described by using the spectral density due to Davidson and Cole, which incorporates a distribution of correlation times characterized by a width parameter epsilon. The fitting of the data was obtained for epsilon = 0.74, indicating that the correlation times are distributed over a narrow range in this system. High-resolution C-13 NMR techniques were used to resolve the NMR lines of middle-chain methylene groups in the spectra and variable contact time cross-polarization {H-1-}C-13 experiments were employed to analyze the reorientation dynamics of the CH3 and CH2 groups in the HDA chains.

Normalization procedure for relaxation studies in NMR quantum information processing

NMR quantum information processing studies rely on the reconstruction of the density matrix representing the so-called pseudo-pure states (PPS). An initially pure part of a PPS state undergoes unitary and non-unitary (relaxation) transformations during a computation process, causing a ""loss of purity"" until the equilibrium is reached. Besides, upon relaxation, the nuclear polarization varies in time, a fact which must be taken into account when comparing density matrices at different instants. Attempting to use time-fixed normalization procedures when relaxation is present, leads to various anomalies on matrices populations. On this paper we propose a method which takes into account the time-dependence of the normalization factor. From a generic form for the deviation density matrix an expression for the relaxing initial pure state is deduced. The method is exemplified with an experiment of relaxation of the concurrence of a pseudo-entangled state, which exhibits the phenomenon of sudden death, and the relaxation of the Wigner function of a pseudo-cat state.

Billiards in a General Domain with Random Reflections

We study stochastic billiards on general tables: a particle moves according to its constant velocity inside some domain D R(d) until it hits the boundary and bounces randomly inside, according to some reflection law. We assume that the boundary of the domain is locally Lipschitz and almost everywhere continuously differentiable. The angle of the outgoing velocity with the inner normal vector has a specified, absolutely continuous density. We construct the discrete time and the continuous time processes recording the sequence of hitting points on the boundary and the pair location/velocity. We mainly focus on the case of bounded domains. Then, we prove exponential ergodicity of these two Markov processes, we study their invariant distribution and their normal (Gaussian) fluctuations. Of particular interest is the case of the cosine reflection law: the stationary distributions for the two processes are uniform in this case, the discrete time chain is reversible though the continuous time process is quasi-reversible. Also in this case, we give a natural construction of a chord ""picked at random"" in D, and we study the angle of intersection of the process with a (d - 1) -dimensional manifold contained in D.

The Aedes aegypti larval transcriptome: a comparative perspective with emphasis on trypsins and the domain structure of peritrophins

The genome sequence of Aedes aegypti was recently reported. A significant amount of Expressed Sequence Tags (ESTs) were sequenced to aid in the gene prediction process. In the present work we describe an integrated analysis of the genomic and EST data, focusing on genes with preferential expression in larvae (LG), adults (AG) and in both stages (SG). A total of 913 genes (5.4% of the transcript complement) are LG, including ion transporters and cuticle proteins that are important for ion homeostasis and defense. From a starting set of 245 genes encoding the trypsin domain, we identified 66 putative LG, AG, and SG trypsins by manual curation. Phylogenetic analyses showed that AG trypsins are divergent from their larval counterparts (LG), grouping with blood-induced trypsins from Anopheles gambiae and Simulium vittatum. These results support the hypothesis that blood-feeding arose only once, in the ancestral Culicomorpha. Peritrophins are proteins that interlock chitin fibrils to form the peritrophic membrane (PM) that compartmentalizes the food in the midgut. These proteins are recognized by having chitin-binding domains with 6 conserved Cys and may also present mucin-like domains (regions expected to be highly O-glycosylated). PM may be formed by a ring of cells (type 2, seen in Ae. aegypti larvae and Drosophila melanogaster) or by most midgut cells (type 1, found in Ae. aegypti adult and Tribolium castaneum). LG and D. melanogaster peritrophins have more complex domain structures than AG and T. castaneum peritrophins. Furthermore, mucin-like domains of peritrophins from T. castaneum (feeding on rough food) are lengthier than those of adult Ae. aegypti (blood-feeding). This suggests, for the first time, that type 1 and type 2 PM may have variable molecular architectures determined by different peritrophins and/or ancillary proteins, which may be partly modulated by diet.

PILZ Protein Structure and Interactions with PILB and the FIMX EAL Domain: Implications for Control of Type IV Pilus Biogenesis

The PilZ protein was originally identified as necessary for type IV pilus (T4P) biogenesis. Since then, a large and diverse family of bacterial PilZ homology domains have been identified, some of which have been implicated in signaling pathways that control important processes, including motility, virulence and biofilm formation. Furthermore, many PilZ homology domains, though not PilZ itself, have been shown to bind the important bacterial second messenger bis(3`-> 5`)cyclic diGMP (c-diGMP). The crystal structures of the PilZ orthologs from Xanthomonas axonopodis pv Citri (PilZ(XAC1133), this work) and from Xanthomonas campestris pv campestris (XC1028) present significant structural differences to other PilZ homologs that explain its failure to bind c-diGMP. NMR analysis of PilZ(XAC1133) shows that these structural differences are maintained in solution. In spite of their emerging importance in bacterial signaling, the means by which NZ proteins regulate specific processes is not clear. In this study, we show that PilZ(XAC1133) binds to PilB, an ATPase required for TV polymerization, and to the EAL domain of FiMX(XAC2398), which regulates TV biogenesis and localization in other bacterial species. These interactions were confirmed in NMR, two-hybrid and far-Western blot assays and are the first interactions observed between any PilZ domain and a target protein. While we were unable to detect phosphodiesterase activity for FimXX(AC2398) in vitro, we show that it binds c-diGMP both in the presence and in the absence of PilZ(XAC1133). Site-directed mutagenesis studies for conserved and exposed residues suggest that PilZ(XAC1133) interactions with FimX(XAC2398) and PilB(XAC3239) are mediated through a hydrophobic surface and an unstructured C-terminal extension conserved only in PilZ orthologs. The FimX-PilZ-PilB interactions involve a full set of ""degenerate"" GGDEF, EAL and PilZ domains and provide the first evidence of the means by which PilZ orthologs and FimX interact directly with the TP4 machinery. (C) 2009 Elsevier Ltd. All rights reserved.

(1)H NMR spectroscopy as a tool to determine accurate photoisomerization quantum yields of stilbene-like ligands coordinated to rhenium(I) polypyridyl complexes

In this work, the use of proton nuclear magnetic resonance, (1)H NMR, was fully described as a powerful tool to follow a photoreaction and to determine accurate quantum yields, so called true quantum yields (Phi(true)), when a reactant and photoproduct absorption overlap. For this, Phi(true) for the trans-cis photoisomerization process were determined for rhenium(I) polypyridyl complexes, fac-[Re(CO)(3)(NN)(trans-L)](+) (NN = 1,10-phenanthroline, phen, or 4,7-diphenyl-1,10-phenanthroline, ph(2)phen, and L = 1,2-bis(4-pyridyl) ethylene, bpe, or 4-styrylpyridine, stpy). The true values determined at 365 nm irradiation (e. g. Phi(NMR) = 0.80 for fac-[Re(CO)(3)(phen)(trans-bpe)](+)) were much higher than those determined by absorption spectral changes (Phi(UV-Vis) = 0.39 for fac-[Re(CO)(3)(phen)(trans-bpe)](+)). Phi(NMR) are more accurate in these cases due to the distinct proton signals of trans and cis-isomers, which allow the actual determination of each component concentration under given irradiation time. Nevertheless when the photoproduct or reactant contribution at the probe wavelength is negligible, one can determine Phi(true) by regular absorption spectral changes. For instance, Phi(313) nm for free ligand photoisomerization determined both by absorption and (1)H NMR variation are equal within the experimental error (bpe: Phi(UV-Vis) = 0.27, Phi(NMR) = 0.26; stpy: Phi(UV-Vis) = 0.49, Phi(NMR) = 0.49). Moreover, (1)H NMR data combined with electronic spectra allowed molar absorptivity determination of difficult to isolate cis-complexes. (C) 2009 Elsevier B. V. All rights reserved.

Catalytic ethanolysis of soybean oil with immobilized lipase from Candida antarctica and (1)H NMR and GC quantification of the ethyl esters (biodiesel) produced

The catalytic ethanolysis of soybean oil with commercial immobilized lipase type B from Candida antarctica to yield ethyl esters (biodiesel) has been investigated. Transesterification was monitored with respect to the following parameters: quantity of biocatalyst, reaction time, amount of water added and turnover of lipase. The highest yields of biodiesel (87% by (1)H NMR; 82.9% by GC) were obtained after a reaction time of 24 h at 32 degrees C in the presence of lipase equivalent to 5.0% (w/w) of the amount of soybean oil present. The production of ethyl esters by enzymatic ethanolysis was not influenced by the addition of water up to 4.0% (v/v) of the alcohol indicating that it is possible to use hydrated ethanol in the production of biodiesel catalyzed by lipase. The immobilized enzyme showed high stability under moderate reaction conditions and retained its activity after five production cycles. The (1)H NMR methodology elaborated for the quantification of biodiesel in unpurified reaction mixtures showed good correlations between the signal areas of peaks associated with the alpha-methylene groups of the ethyl esters and those of the triacyl-glycerides in residual soybean oil. Monoacylglycerides, diacylglycerides and triglycerides could also be detected and quantified in the crude biodiesel using (1)H NMR spectroscopic and GC-FID chromatographic methods. The biodiesel production by enzymatic catalysis was promising. In this case, was produced a low concentration of glycerol (0.74%) and easily removed by water extraction. (C) 2010 Elsevier B.V. All rights reserved.

Nuclear magnetic resonance water relaxation time changes in bananas during ripening: a new mechanism

BACKGROUND: Nuclear magnetic resonance studies of banana fragments during ripening show an increase on the water transverse relaxation time (T(2)) and a decrease in water self-diffusion coefficient (D). As T(2) and D are normally directly correlated, we studied these two properties in intact bananas during ripening, in an attempt to rule out the effect of injury on the apparent discrepancies in the behavior of T(2) and D. RESULTS: The results show that injury in bananas causes a decrease in T(2) of the water in vacuoles (T(2vac)). They also show that T(2vac) increased and D decreased during ripening, ruling out the injury effect. To explain the apparent discrepancies, we propose a new hypothesis for the increase in T(2) values, based on the reduction of Fe(3+) ions to Fe(2+) by galacturonic acid, produced by the hydrolysis of pectin and a decrease in internal oxygen concentration during ripening. CONCLUSION: As injury alters T(2) values it is necessary to use intact bananas to study relaxation times during ripening. The novel interpretation for the increase in T(2vac) based on reduction of Fe(+3) and O(2) concentration is an alternative mechanism to that based on the hydrolysis of starch in amyloplasts. (C) 2010 Society of Chemical Industry