Signal enhancement in 5QMAS spectra of spin-5/2 quadrupolar nuclei

5QMAS experiments on spin-5/2 systems display a low sensitivity compared with their 3QMAS counterparts. Nevertheless, the superior resolution of 5QMAS over 3QMAS makes these experiments a favorable choice for many materials. We report an enhancement scheme for the 5QMAS experiment, using an improved five-quantum excitation pulse scheme combined with a FAM-II conversion pulse. The results are verified experimentally on a polycrystalline sample of gamma-(27)Al(2)O(3), showing an enhancement factor of 2.4 over the simple two-pulse (CW) 5QMAS scheme. Numerical computations of the efficiency parameter epsilon support these results.     

Phosporus-31 Solid-State NMR in High-Field Gradients: Prospects for Imaging Bone Using the Long Echo-Train Summation Technique (LETS)

Stray-field techniques are reported for ³¹P studies of solids for a variety of compounds including bone, bone meal and calcium hydroxyapatite. Long Hahn echo trains produced by the application of many pulses were used as in the long echo-train summation technique. Double-resonance enhancements of ³¹P by use of both direct and indirect experiments were attempted on a sample of NH₄PF₆:³¹P{¹⁹F} double resonance produced, at most, a 26% increase in the initial level of the ³¹P echo signal.     

An open magnet utilizing ferro-refraction current magnification

Ferro-refraction is the field magnification that is obtained when a current segment is near a high magnetic permeable boundary. It is shown that ferro-refraction may be used in the design of magnets for NMR or MRI to increase the efficiency of these magnets. The field may be modeled analytically with the Biot--Savart law and the inclusion of mirror image currents. Ferro-refraction is particularly useful in the design of monohedral magnets, magnets producing a remote homogeneous region which have the magnetic sources arranged to one side. These magnets have also been called planar magnets. Two designs for a monohedral magnet which produce good agreement between experimental and analytic results are presented.     

Use of the Carbonyl Chemical Shift to Relieve Degeneracies in Triple-Resonance Assignment Experiments

We illustrate an approach that uses the backbone carbonyl chemical shift to relieve resonance overlaps in triple-resonance assignment experiments conducted on protein samples. We apply this approach to two cases of simultaneous overlaps: those of ((1)H(N), (15)N) spin pairs and those of ((1)H(alpha), (13)C(alpha)) spin pairs in residues preceding prolines. For these cases we employed respectively CBCACO(N)H and H(CA)CON experiments, simple variants of the commonly used CBCA(CO)NH and HCA(CO)N experiments obtained by replacing one of the indirect dimensions with a carbonyl dimension. We present data collected on ribosomal protein S4 using these experiments, along with overlap statistics for four other polypeptides ranging in size from 76 to 263 residues. These data indicate that the CBCACO(N)H, in combination with the CBCA(CO)NH, can relieve >83% of the ((1)H(N), (15)N) and ((1)H(N), (13)C') overlaps for these proteins. The data also reveal how the H(CA)CON experiment successfully completed the assignment of triply and quadruply degenerate X-Pro spin systems in a mobile, proline-rich region of S4, even when X was a glycine. Finally, we discuss the relative sensitivities of these experiments compared to those of existing sequences, an analysis that reinforces the usefulness of these experiments in assigning extensively overlapped and/or proline-rich sequences in proteins.     

Numerical simulation of one- and two-dimensional ESEEM experiments

Numerical simulation has become an indispensable tool for the interpretation of pulse EPR experiments. In this work it is shown how automatic orientation selection, grouping of operator factors, and direct selection and elimination of coherences can be used to improve the efficiency of time-domain simulations of one- and two-dimensional electron spin echo envelope modulation (ESEEM) spectra. The program allows for the computation of magnetic interactions of any symmetry and can be used to simulate spin systems with an arbitrary number of nuclei with any spin quantum number. Experimental restrictions due to finite microwave pulse lengths are addressed and the enhancement of forbidden coherences by microwave pulse matching is illustrated. A comparison of simulated and experimental HYSCORE (hyperfine sublevel correlation) spectra of ordered and disordered systems with varying complexity shows good qualitative agreement.     

Mass Flux Schemes and Connection to Shock Instability

We analyze numerical mass fluxes with an emphasis on their capability for accurately capturing shock and contact discontinuities. The study of mass flux is useful because it is the term common to all conservation equations and the numerical diffusivity introduced in it bears a direct consequence to the prediction of contact (stationary and moving) discontinuities, which are considered to be the limiting case of the boundary layer. We examine several prominent numerical flux schemes and analyze the structure of numerical diffusivity. This leads to a detailed investigation into the cause of certain catastrophic breakdowns by some numerical flux schemes. In particular, we identify the dissipative terms that are responsible for shock instabilities, such as the odd–even decoupling and the so-called “carbuncle phenomenon”. As a result, we propose a conjecture stating the connection of the pressure difference term to these multidimensional shock instabilities and hence a cure to those difficulties. The validity of this conjecture has been confirmed by examining a wide class of upwind schemes. The conjecture is useful to the flux function development, for it indicates whether the flux scheme under consideration will be afflicted with these kinds of failings. Thus, a class of shock-stable schemes can be identified. Interestingly, a shock-stable scheme's self-correcting capability is demonstrated with respect to carbuncle-contaminated profiles for flows at both low supersonic and high Mach numbers.     

Design and implementation of a protection system for NMR spectrometers

We have implemented a scheme, SPECMON, for monitoring various parameters of a spectrometer, such as nitrogen pressure and sample temperature, and taking corrective action. The scheme is based on considerations of protection management which are of general application. Evaluation of the spectrometer state is incorporated in macros of the application software (VNMR) and is therefore very flexible. In contrast, corrective action is limited to the single one which is deemed fully safe: complete shutdown of the spectrometer and logging. Shutdown is implemented by a minor hardware modification of the spectrometer: the introduction of a second input to a relay already present for protection of the spectrometer power supply. Monitoring is handled by the host computer, and the shutdown command is transmitted via control lines of its series port, independent of the standard connection between the host computer and the NMR system console. The monitoring system (software and hardware) is unobtrusive in normal conditions, and it can be tested without affecting the operation of the spectrometer.     

H NMR Characterization of the Product from Single Solid-Phase Resin Beads Using Capillary NMR Flow Probes

A capillary NMR flow probe was designed to generate high-resolution ¹H NMR spectra at 600 MHz from the cleaved product of individual 160-μm Tentagel combinatorial chemistry beads. By injecting a dissolved sample sandwiched between an immiscible, perfluorinated organic liquid directly into the probe, NMR spectra of the product cleaved from single beads were acquired in just 1 h of spectrometer time without diffusional dilution. Sample handling efficiency on the single bead scale was comparable to that obtained with a bulk sample. Using the relative intensity of the DMSO-d₅H versus the analyte signals in a fully relaxed CPMG spectrum, the amount of product cleaved from a single bead was determined to be 540±170 pmol in one of the samples. Following the NMR data collection, the samples were examined with electrospray ionization mass spectrometry to provide additional structural information. By coupling with microliter-volume fluidic capabilities, the capillary flow probe described here will enable multidimensional characterization of single solid-phase resin products in an online manner.     

A rotational approach to localized SPAMM 1-1 tagging

Magnetic resonance tagging usually relies on controlling the phase dispersion of the transverse magnetization component. Phase dispersion is, however, affected by the inherent phase of selective excitation pulses, thus limiting their combination with tagging sequences to the application of refocusable pulses, as in the localized spatial modulation of magnetization (L-SPAMM) technique. In this study, we examine the effect of selective excitation pulses on a L-SPAMM 1-1 sequence, showing that in the case of two identical pulses the phase component is canceled out, and thus preemphasis and refocus gradients are not needed, allowing us to take advantage of a constant gradient throughout the tagging sequence, and also that one might choose nonrefocusable maximum and minimum phase pulses.     

Determination of membrane Peptide orientation by 1H-detected 2H NMR spectroscopy

We demonstrate the application of the proton inverse detected deuteron (PRIDE) NMR technique to the measurement of the orientation of membrane-bound peptides with enhanced sensitivity. Gramicidin D, a transmembrane peptide, and ovispirin, a surface-bound peptide, were used as model systems. The peptides were 2H-labeled by 1H/2H exchange and oriented uniaxially on glass plates. The directly detected 2H spectra of both peptides showed only a strong D(2)O signal and no large quadrupolar splittings. In contrast, the PRIDE spectrum of gramicidin exhibited quadrupolar splittings as large as 281 kHz, consistent with its transmembrane orientation. Moreover, the large D(2)O signal in the directly detected 2H spectra was cleanly suppressed in the PRIDE spectrum. For ovispirin, the 1H indirectly detected 2H spectrum revealed a 104 kHz splitting and a zero-frequency peak. The former reflects the in-plane orientation of most of the helix axis, while the latter results from residues with a magic-angle orientation of the N-D bonds. These are consistent with previous 15N NMR results on ovispirin. The combination of PRIDE and exchange labeling provides an economical and sensitive method of studying membrane peptide orientations in lipid bilayers without the influence of D(2)O and with the ability to detect N-D bonds at the magic angle from the bilayer normal.     

An Extension of Alias Sampling Method for Parametrized Probability Distributions

An extension of the alias sampling technique for distribution functions depending on a number of parameters was developed. It takes advantage of modern computer architectures with large amounts of cheap memory, by using discrete representations of probability distribution functions. The sampling is done by fast interpolation techniques involving only elementary logical and arithmetical operations, allowing one to keep a higher degree of accuracy as the grids spacing is controlled by the user. By this method it is possible to obtain the value of interest by direct interpolation between the sampled values obtained with the same set of random numbers for the grid values of the parameters adjacent to the values of interest. Sampling tests carried for the case of Molière electron multi-scatter angle distribution show that this method can be successfully used in Monte Carlo codes for sampling complex probability distributions.     

Application of high-resolution magic-angle spinning NMR spectroscopy to define the cell uptake of MRI contrast agents

A new method, based on proton high-resolution magic-angle spinning ((1)H HR-MAS) NMR spectroscopy, has been employed to study the cell uptake of magnetic resonance imaging contrast agents (MRI-CAs). The method was tested on human red blood cells (HRBC) and white blood cells (HWBC) by using three gadolinium complexes, widely used in diagnostics, Gd-BOPTA, Gd-DTPA, and Gd-DOTA, and the analogous complexes obtained by replacing Gd(III) with Dy(III), Nd(III), and Tb(III) (i.e., complexes isostructural to the ones of gadolinium but acting as shift agents). The method is based on the evaluation of the magnetic effects, line broadening, or induced lanthanide shift (LIS) caused by these complexes on NMR signals of intra- and extracellular water. Since magnetic effects are directly linked to permeability, this method is direct. In all the tests, these magnetic effects were detected for the extracellular water signal only, providing a direct proof that these complexes are not able to cross the cell membrane. Line broadening effects (i.e., the use of gadolinium complexes) only allow qualitative evaluations. On the contrary, LIS effects can be measured with high precision and they can be related to the concentration of the paramagnetic species in the cellular compartments. This is possible because the HR-MAS technique provides the complete elimination of bulk magnetic susceptibility (BMS) shift and the differentiation of extra- and intracellular water signals. Thus with this method, the rapid quantification of the MRI-CA amount inside and outside the cells is actually feasible.     

IR Spectrum of C8H2: Integrated Band Intensities and Some Observational Implications

Polyynes are of astrophysical interest since they appear to be involved in organic chemistry in very different mediums. In Titan's atmosphere, the lightest polyyne, C₄H₂, was detected by Voyager. Recently C₄H₂ and C₆H₂ have been discovered in a protoplanetary nebula, suggesting polyynes as a possible chemical pathway to PAH (polycyclic aromatic hydrocarbons). Moreover, several experimental simulations and modeling imply their production from the photochemistry of methane and their involvement in the formation of organic aerosols. After the study of C₄H₂ and C₆H₂ spectra in the UV and IR wavelength range, we report here the first spectrum of gaseous C₈H₂ in the range 400–4000 cm⁻¹ at room temperature and low resolution. The task was hardly achieved because of the high instability of this molecule with temperature and pressure. To avoid exothermic polymerization, the compound as mixed with a solvent. We have performed a separate spectroscopic study of the solvent to determine C₈H₂ partial pressure within the mixture. This allowed us to calculate C₈H₂ integrated band intensities. In the studied wavelength range, C₈H₂ presents three main bands similar to those of C₆H₂ in terms of vibrational type, position, and relative intensity. To study the possible identification of these polyynes by spatial observatories (Cassini–Huygens, ISO), we have also measured the C₆H₂ and C₈H₂ infrared spectra in the range 400–1500 cm⁻¹ at 0.35 cm⁻¹ resolution.     

H Spectroscopy without Solvent Suppression: Characterization of Signal Modulations at Short Echo Times

While most proton (¹H) spectra acquired in vivo utilize selective suppression of the solvent signal for more sensitive detection of signals from the dilute metabolites, recent reports have demonstrated the feasibility and advantages of collecting in vivo data without solvent attenuation. When these acquisitions are performed at short echo times, the presence of frequency modulations of the water resonance may become an obstacle to the identification and quantitation of metabolite resonances. The present report addresses the characteristics, origin, and elimination of these sidebands. Sideband amplitudes were measured as a function of delay time between gradient pulse and data collection, as a function of gradient pulse amplitude, and as a function of spatial location of the sample for each of the three orthogonal gradient sets. Acoustic acquisitions were performed to demonstrate the correlation between mechanical vibration resonances and the frequencies of MR sidebands. A mathematical framework is developed and compared with the experimental results. This derivation is based on the theory that these frequency modulations are induced by magnetic field fluctuations generated by the transient oscillations of gradient coils.     

The Lattice Solid Model to Simulate the Physics of Rocks and Earthquakes: Incorporation of Friction

The particle-based lattice solid model developed to study the physics of rocks and the nonlinear dynamics of earthquakes is refined by incorporating intrinsic friction between particles. The model provides a means for studying the causes of seismic wave attenuation, as well as frictional heat generation, fault zone evolution, and localisation phenomena. A modified velocity–Verlat scheme that allows friction to be precisely modelled is developed. This is a difficult computational problem given that a discontinuity must be accurately simulated by the numerical approach (i.e., the transition from static to dynamical frictional behaviour). This is achieved using a half time step integration scheme. At each half time step, a nonlinear system is solved to compute the static frictional forces and states of touching particle-pairs. Improved efficiency is achieved by adaptively adjusting the time step increment, depending on the particle velocities in the system. The total energy is calculated and verified to remain constant to a high precision during simulations. Numerical experiments show that the model can be applied to the study of earthquake dynamics, the stick–slip instability, heat generation, and fault zone evolution. Such experiments may lead to a conclusive resolution of the heat flow paradox and improved understanding of earthquake precursory phenomena and dynamics.     

A subspace time-domain algorithm for automated NMR spectral normalization

Recently, two methods have been proposed for quantitatively comparing NMR spectra of control and treated samples, in order to examine the possible occurring variations in cell metabolism and/or structure in response to numerous physical, chemical, and biological agents. These methods are the maximum superposition normalization algorithm (MaSNAl) and the minimum rank normalization algorithm (MiRaNAl). In this paper a new subspace-based time-domain normalization algorithm, denoted by SuTdNAl (subspace time-domain normalization algorithm), is presented. By the determination of the intersection of the column spaces of two Hankel matrices, the common signal poles and further on the components having proportionally varying amplitudes are detected. The method has the advantage that it is computationally less intensive than the MaSNAl and the MiRaNAl. Furthermore, no approximate estimate of the normalization factor is required. The algorithm was tested by Monte Carlo simulations on a set of simulation signals. It was shown that the SuTdNAl has a statistical performance similar to that of the MiRaNAl, which itself is an improvement over the MaSNAl. Furthermore, two samples of known contents are compared with the MiRaNAl, the SuTdNAl, and an older method using a standard. Finally, the SuTdNAl is tested on a realistic simulation example derived from an in vitro measurement on cells.     

Probe--sample coupling in the magnetic resonance force microscope

The magnetic resonance force microscope (MRFM) provides a route to achieving scanned probe magnetic resonance imaging with extremely high spatial resolution. Achieving this capability will require understanding the force exerted on a microscopic magnetic probe by a spatially extended sample over which the probe is scanned. Here we present a detailed analysis of this interaction between probe and sample. We focus on understanding the situation where the micromagnet mounted on the mechanical resonator generates a very inhomogeneous magnetic field and is scanned over a sample with at least one spatial dimension much larger than that of the micromagnet. This situation differs quite significantly from the conditions under which most MRFM experiments have been carried out where the sample is mounted on the mechanical resonator and placed in a rather weak magnetic field gradient. In addition to the concept of a sensitive slice (the spatial region where the magnetic resonance condition is met) it is valuable to map the forces exerted on the probe by spins at various locations; this leads to the concept of the force slice (the region in which spins exert force on the resonator). Results of this analysis, obtained both analytically and numerically, will be qualitatively compared with an initial experimental finding from an EPR-MRFM experiment carried out on DPPH at 4 K.     

A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows

We propose a new model and a solution method for two-phase compressible flows. The model involves six equations obtained from conservation principles applied to each phase, completed by a seventh equation for the evolution of the volume fraction. This equation is necessary to close the overall system. The model is valid for fluid mixtures, as well as for pure fluids. The system of partial differential equations is hyperbolic. Hyperbolicity is obtained because each phase is considered to be compressible. Two difficulties arise for the solution: one of the equations is written in non-conservative form; non-conservative terms exist in the momentum and energy equations. We propose robust and accurate discretisation of these terms. The method solves the same system at each mesh point with the same algorithm. It allows the simulation of interface problems between pure fluids as well as multiphase mixtures. Several test cases where fluids have compressible behavior are shown as well as some other test problems where one of the phases is incompressible. The method provides reliable results, is able to compute strong shock waves, and deals with complex equations of state.     

Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy

In this paper we design a class of numerical schemes that are higher-order extensions of the weighted essentially non-oscillatory (WENO) schemes of G.-S. Jiang and C.-W. Shu (1996) and X.-D. Liu, S. Osher, and T. Chan (1994). Used by themselves, the schemes may not always be monotonicity preserving but coupled with the monotonicity preserving bounds of A. Suresh and H. T. Huynh (1997) they perform very well. The resulting monotonicity preserving weighted essentially non-oscillatory (MPWENO) schemes have high phase accuracy and high order of accuracy. The higher-order members of this family are almost spectrally accurate for smooth problems. Nevertheless, they, have robust shock capturing ability. The schemes are stable under normal CFL numbers. They are also efficient and do not have a computational complexity that is substantially greater than that of the lower-order members of this same family of schemes. The higher accuracy that these schemes offer coupled with their relatively low computational complexity makes them viable competitors to lower-order schemes, such as the older total variation diminishing schemes, for problems containing both discontinuities and rich smooth region structure. We describe the MPWENO schemes here as well as show their ability to reach their designed accuracies for smooth flow. We also examine the role of steepening algorithms such as the artificial compression method in the design of very high order schemes. Several test problems in one and two dimensions are presented. For multidimensional problems where the flow is not aligned with any of the grid directions it is shown that the present schemes have a substantial advantage over lower-order schemes. It is argued that the methods designed here have great utility for direct numerical simulations and large eddy simulations of compressible turbulence. The methodology developed here is applicable to other hyperbolic systems, which is demonstrated by showing that the MPWENO schemes also work very well on magnetohydrodynamical test problems.     

Solution of Time-Dependent Diffusion Equations with Variable Coefficients Using Multiwavelets

A new numerical algorithm is developed for the solution of time-dependent differential equations of diffusion type. It allows for an accurate and efficient treatment of multidimensional problems with variable coefficients, nonlinearities, and general boundary conditions. For space discretization we use the multiwavelet bases introduced by Alpert (1993,SIAM J. Math. Anal.24, 246–262), and then applied to the representation of differential operators and functions of operators presented by Alpert, Beylkin, and Vozovoi (Representation of operators in the multiwavelet basis, in preparation). An important advantage of multiwavelet basis functions is the fact that they are supported only on non-overlapping subdomains. Thus multiwavelet bases are attractive for solving problems in finite (non periodic) domains. Boundary conditions are imposed with a penalty technique of Hesthaven and Gottlieb (1996,SIAM J. Sci. Comput., 579–612) which can be used to impose rather general boundary conditions. The penalty approach was extended to a procedure for ensuring the continuity of the solution and its first derivative across interior boundaries between neighboring subdomains while time stepping the solution of a time dependent problem. This penalty procedure on the interfaces allows for a simplification and sparsification of the representation of differential operators by discarding the elements responsible for interactions between neighboring subdomains. Consequently the matrices representing the differential operators (on the finest scale) have block-diagonal structure. For a fixed order of multiwavelets (i.e., a fixed number of vanishing moments) the computational complexity of the present algorithm is proportional to the number of subdomains. The time discretization method of Beylkin, Keiser, and Vozovoi (1998, PAM Report 347) is used in view of its favorable stability properties. Numerical results are presented for evolution equations with variable coefficients in one and two dimensions.