Dipolar and J encoded DQ MAS spectra under rotational resonance

A two-dimensional (2D) double-quantum (DQ) experiment under rotational resonance (R(2)) conditions is introduced for evaluating dipolar couplings in rotating solids. The contributions from the R(2)-recoupled dipolar interaction and the J coupling can be conveniently separated in the resulting 2D R(2)-DQ spectrum, so that the unknown dipolar coupling can readily be extracted, provided that the values of the involved J coupling constants are known. Since the measured parameters are integral intensity ratios between suitably chosen absorption peaks in the 2D spectrum, the proposed method is characterized by a reduced sensitivity to relaxation parameters. The effect of rotor-modulated terms, including chemical shift anisotropy, is efficiently averaged out by synchronizing the excitation/reconversion time with the rotor period. All of these features are demonstrated theoretically by the example of two model systems, namely, isolated spin-pairs and a three-spin system. The results of the theoretical models are applied to both (13)C and (1)H nuclei to extract dipolar couplings in uniformly (13)C labeled L-alanine and a crosslinked natural rubber.     

Wave Simulation in Frozen Porous Media

We propose a numerical algorithm for simulation of wave propagation in frozen porous media, where the pore space is filled with ice and water. The model, based on a Biot-type three-phase theory, predicts three compressional waves and two shear waves and models the attenuation level observed in rocks. Attenuation is modeled with exponential relaxation functions which allow a differential formulation based on memory variables. The wavefield is obtained using a grid method based on the Fourier differential operator and a Runge–Kutta time-integration algorithm. Since the presence of slow quasistatic modes makes the differential equations stiff, a time-splitting integration algorithm is used to solve the stiff part analytically. The modeling is second-order accurate in the time discretization and has spectral accuracy in the calculation of the spatial derivatives.     

Linewidth-Resolved N HSQC, a Simple 3D Method to Measure N Relaxation Times from T1 and T2 Linewidths

A three-dimensional approach for measuring ¹⁵N relaxation times is described. Instead of selecting particular values for the relaxation period, in the proposed method the relaxation period is incremented periodically in order to create a 3D spectrum. This additional frequency domain of the transformed spectrum contains the relaxation time information in the T₁ and T₂ linewidths, and thus the longitudinal and transverse ¹⁵N relaxation times can be measured without determination of 2D cross peak volumes/intensities and subsequent curve fitting procedures.     

Changes in Porcine Muscle Water Characteristics during Growth—An in Vitro Low-Field NMR Relaxation Study

This study investigates the effects of developmental stage and muscle type on the mobility and distribution of water within skeletal muscles, using low-field ¹H-NMR transverse relaxation measurements in vitro on four different porcine muscles (M. longissimus dorsi, M. semitendinosus, M. biceps femoris, M. vastus intermedius) from a total of 48 pigs slaughtered at various weight classes between 25 kg and 150 kg. Principal component analysis (PCA) revealed effects of both slaughter weight and muscle type on the transverse relaxation decay. Independent of developmental stage and muscle type, distributed exponential analysis of the NMR T₂ relaxation data imparted the existence of three distinct water populations, T_2b, T₂₁, and T₂₂, with relaxation times of approximately 1–10, 45–120, and 200–500 ms, respectively. The most profound change during muscle growth was a shift toward faster relaxation in the intermediate time constant, T₂₁. It decreased by approx. 24% in all four muscle types during the period from 25 to 150 kg live weight. Determination of dry matter, fat, and protein content in the muscles showed that the changes in relaxation time of the intermediate time constant, T₂₁, during growth should be ascribed mainly to a change in protein content, as the protein content explained 77% of the variation in the T₂₁ time constant. Partial least squares (PLS) regression revealed validated correlations in the region of 0.58 to 0.77 between NMR transverse relaxation data and muscle development for all the four muscle types, which indicates that NMR relaxation measurements may be used in the prediction of muscle developmental stage.     

Relaxation Behavior of Laser-Polarized Xe Gas: Size Dependency and Wall Effect of the T1 Relaxation Time in Glass and Gelatin Bulbs

Size dependency of the relaxation time T₁ was measured for laser-polarized ¹²⁹Xe gas encapsulated in different sized cavities made by glass bulbs or gelatin capsules. The use of laser-polarized gas enhances the sensitivity a great deal, making it possible to measure the longer ¹²⁹Xe relaxation time in quite a short time. The size dependency is analyzed on the basis of the kinetic theory of gases and a relationship is derived in which the relaxation rate is connected with the square inverse of the diameter of the cavity. Such an analysis provides a novel parameter which denotes the wall effect on the relaxation rate when a gas molecule collides with the surface once in a second. The relaxation time of ¹²⁹Xe gas is also dependent on the material which forms the cavity. This dependency is large and the relaxation study using polarized ¹²⁹Xe gas is expected to offer important information about the state of the matter of the cavity wall.     

Measuring the longitudinal NMR relaxation rates of fast relaxing nuclei using a signal eliminating relaxation filter

A new experiment for selective determination of the relaxation rates of fast relaxing NMR signals is presented. The experiment is derived from the conventional inversion recovery experiment by substituting the 180 degrees inversion pulse of this experiment with a signal eliminating relaxation filter (SERF) consisting of three 180 degrees pulses separated by two variable delays, Delta1 and Delta2. The SERF experiment allows a selective suppression of signals with relaxation rates below a given limit while monitoring the relaxation of faster relaxing signals. The experiment was tested on a sample of 20% oxidized plastocyanin from Anabaena variabilis, where the fast exchange of an electron between the reduced (diamagnetic) and the oxidized (paramagnetic) form results in a series of average signals with widely different relaxation rates. To ensure an optimum extraction of information from the experimental data, the relaxation rates were obtained from the SERF experiment by a simultaneous analysis of all the FIDs of the experiment using a fast linear prediction model method developed previously. The reliability of the relaxation rates obtained from the SERF experiment was confirmed by a comparison of the rates with the corresponding rates obtained from a conventional inversion recovery experiment.     

The effects of dissolved oxygen upon amide proton relaxation and chemical shift in a perdeuterated protein

The effects of dissolved molecular oxygen upon amide proton ((1)H(N)) longitudinal and transverse relaxation rates and chemical shifts were studied for a small protein domain, the second type 2 module of fibronectin ((2)F2)-isotopically enriched to 99% (2)H, 98% (15)N. Longitudinal relaxation rate enhancements, R(O(2))((1)H(N)), of individual backbone (1)H(N) nuclei varied up to 14 fold between a degassed and oxygenated (1 bar) solution, indicating that the oxygen distribution within the protein is inhomogeneous. On average, smaller relaxation rate enhancements were observed for (1)H(N) nuclei associated with the core of the protein compared to (1)H(N) nuclei closer to the surface, suggesting restricted oxygen accessibility to some regions. In agreement with an O(2)-(1)H(N) hyperfine interaction in the extreme narrowing limit, the (1)H(N) transverse relaxation rates showed no significant change, up to an oxygen pressure of 9.5 bar (the maximum pressure used in this study). For most (1)H(N) resonances, small deltadelta(O(2))((1)H(N)) hyperfine chemical shifts could be detected between oxygen pressures of 1 bar and 9.5 bar.     

NMR Relaxation and Diffusion Measurements on Iron(III)-Doped Kaolin Clay

Kaolin clay samples were mixed with various amounts of Fe₂O₃ powder. The influence of this magnetic impurity on NMR relaxation and diffusion measurements on the water in this porous material was investigated. The NMR relaxation measurements showed a nearly mono-exponential decay, leading to the conclusion that the pore size distribution of the clay samples is either narrow and/or that the pores are interconnected very well. Both the longitudinal and the transverse relaxation rate depend linearly on the concentration of the Fe₂O₃ impurity. The NMR diffusion measurements revealed that the Fe₂O₃ causes internal magnetic field gradients that largely exceed the maximum external gradient that could be applied by our NMR apparatus (0.3 T/m). Additional SQUID measurements yielded the magnetization and magnetic susceptibility of the samples at the magnetic field strength used in the NMR measurements (0.8 T). A theoretical estimate of the internal magnetic field gradients leads to the conclusion that the water in the porous clay samples cannot be described by the commonly observed motional averaging regime. Probably an intermediate or a localization regime is induced by the large internal gradients, which are estimated to be on the order of 1 to 10 T/m in the pore volume and may exceed 1000 T/m at the pore surface.     

Time-optimal coherence-order-selective transfer of in-phase coherence in heteronuclear IS spin systems

Based on principles of geometric optimal control theory, coherence transfer building blocks can be derived which achieve optimal sensitivity. Here, experimental pulse sequences are presented that achieve the best possible coherence-order-selective in-phase transfer (S(-)-->I(-)) for a heteronuclear 2-spin system for any given mixing time in the absence of relaxation. For short mixing times, the optimal experiment improves the sensitivity of isotropic mixing by up to 12.5%.     

Implementation Aspects of 3D Lattice-BGK: Boundaries, Accuracy, and a New Fast Relaxation Method

In many realistic fluid-dynamical simulations the specification of the boundary conditions, the error sources, and the number of time steps to reach a steady state are important practical considerations. In this paper we study these issues in the case of the lattice-BGK model. The objective is to present a comprehensive overview of some pitfalls and shortcomings of the lattice-BGK method and to introduce some new ideas useful in practical simulations. We begin with an evaluation of the widely used bounce-back boundary condition in staircase geometries by simulating flow in an inclined tube. It is shown that the bounce-back scheme is first-order accurate in space when the location of the non-slip wall is assumed to be at the boundary nodes. Moreover, for a specific inclination angle of 45 degrees, the scheme is found to be second-order accurate when the location of the non-slip velocity is fitted halfway between the last fluid nodes and the first solid nodes. The error as a function of the relaxation parameter is in that case qualitatively similar to that of flat walls. Next, a comparison of simulations of fluid flow by means of pressure boundaries and by means of body force is presented. A good agreement between these two boundary conditions has been found in the creeping-flow regime. For higher Reynolds numbers differences have been found that are probably caused by problems associated with the pressure boundaries. Furthermore, two widely used 3D models, namelyD₃Q₁₅andD₃Q₁₉, are analysed. It is shown that theD₃Q₁₅model may induce artificial checkerboard invariants due to the connectivity of the lattice. Finally, a new iterative method, which significantly reduces the saturation time, is presented and validated on different benchmark problems.     

A New Two-Dimensional Pulse Sequence for T2* Measurements of Protons in C Isotopomers

A new two-dimensional pulse sequence for T₂* measurement of protons directly coupled to ¹³C spins is proposed. The sequence measures the tranverse relaxation time of heteronuclear proton single-quantum coherence under conditions of free precession and is therefore well suited to evaluate relaxation losses of proton magnetization during preparation delays of heteronuclear pulse experiments in analytical NMR. The relevant part of the pulse sequence can be inserted as a “building block” into any direct or inverse detecting H,C correlation pulse sequence if proton spin–spin relaxation is to be investigated. In this contribution, the building block is inserted into a HETCOR as well as into a HMQC pulse sequence. Experimental results for the HETCOR-based sequence are given.     

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.     

Particle Transport through Scattering Regions with Clear Layers and Inclusions

This paper introduces generalized diffusion models for the transport of particles in scattering media with nonscattering inclusions. Classical diffusion is known as a good approximation of transport only in scattering media. Based on asymptotic expansions and the coupling of transport and diffusion models, generalized diffusion equations with nonlocal interface conditions are proposed which offer a computationally cheap, yet accurate, alternative to solving the full phase-space transport equations. The paper shows which computational model should be used depending on the size and shape of the nonscattering inclusions in the simplified setting of two space dimensions. An important application is the treatment of clear layers in near-infrared (NIR) spectroscopy, an imaging technique based on the propagation of NIR photons in human tissues.     

Improved narrowband dipolar recoupling for homonuclear distance measurements in rotating solids

Recovery of the magnetic dipolar interaction between nuclei bearing the same gyromagnetic ratio in rotating solids can be promoted by synchronous rf irradiation. Determination of the dipolar interaction strength can serve as a tool for structural elucidation in polycrystalline powders. Spinning frequency dependent narrow-band (nb) RFDR and SEDRA experiments are utilized as simple techniques for the determination of dipolar interactions between the nuclei in coupled homonuclear spin pairs. The magnetization exchange and coherence dephasing due to a fixed number of rotor-synchronously applied pi-pulses is monitored at spinning frequencies in the vicinity of the rotational resonance (R(2)) conditions. The powder nbRFDR and nbSEDRA decay curves of spin magnetizations and coherences, respectively, as a function of the spinning frequency can be measured and analyzed using simple rate equations providing a quantitative measure of the dipolar coupling. The effects of the phenomenological relaxation parameters in these rate equations are discussed and an improved methodology is suggested for analyzing nbRFDR data for small dipolar couplings. The distance between the labeled nuclei in the 1,3-(13)C(2)-hydroxybutyric acid molecule is rederived using existing nbRFDR results and the new simulation procedure. A nbSEDRA experiment has been performed successfully on a powder sample of singly labeled 1-(13)C-L-leucine measuring the dipolar interaction between the labeled carboxyl carbon and the natural abundant beta-carbon. Both narrowband techniques are employed for the determination of the nuclear distances between the side-chain carbons of leucine and its carbonyl carbon in a tripeptide Leu-Gly-Phe that is singly (13)C-labeled at the leucine carbonyl carbon position.     

Multigrid Solution of the Incompressible Navier–Stokes Equations for Density-Stratified Flow past Three-Dimensional Obstacles

The steady incompressible Navier–Stokes equations in three dimensions are solved for neutral and stably stratified flow past three-dimensional obstacles of increasing spanwise width. The continuous equations are approximated using a finite volume discretisation on staggered grids with a flux-limited monotonic scheme for the advective terms. The discrete equations which arise are solved using a nonlinear multigrid algorithm with up to four grid levels using the SIMPLE pressure correction method as smoother. When at its most effective the multigrid algorithm is demonstrated to yield convergence rates which are independent of the grid density. However, it is found that the asymptotic convergence rate depends on the choice of the limiter used for the advective terms of the density equation, and some commonly used schemes are investigated. The variation with obstacle width of the influence of the stratification on the flow field is described and the results of the three-dimensional computations are compared with those of the corresponding computation of flow over a two-dimensional obstacle (of effectively infinite width). Also given are the results of time-dependent computations for three-dimensional flows under conditions of strong static stability when lee-wave propagation is present and the multigrid algorithm is used to compute the flow at each time step.     

Dynamics of Spin I=3/2 under Spin-Locking Conditions in an Ordered Environment

We have derived approximate analytic solutions to the master equation describing the evolution of the spin I=3/2 density operator in the presence of a radio-frequency (RF) field and both static and fluctuating quadrupolar interactions. Spectra resulting from Fourier transformation of the evolutions of the on-resonance spin-locked magnetization into the various coherences display two satellite pairs and, in some cases, a central line. The central line is generally trimodal, consisting of a narrow component related to a slowly relaxing mode and two broad components pertaining to two faster relaxing modes. The rates of the fast modes are sensitive to slow molecular motion. Neither the amplitude nor the width of the narrow component is affected by the magnitude of the static coupling, whereas the corresponding features of the broad components depend in a rather complicated manner on the spin-lock field strength and static quadrupolar interaction. Under certain experimental conditions, the dependencies of the amplitudes on the dynamics are seen to vanish and the relaxation rates reduce to relatively simple expressions. One of the promising emerging features is the fact that the evolutions into the selectively detected quadrupolar spin polarization order and the rank-two double-quantum coherence do not exhibit a slowly relaxing mode and are particularly sensitive to slow molecular motion. Furthermore, these coherences can only be excited in the presence of a static coupling and this makes it possible to discern nuclei in anisotropic from those in isotropic environment. The feasibility of the spin-lock pulse sequences with limited RF power and a nonvanishing average electric field gradient has been demonstrated through experiments on sodium in a dense lyotropic DNA liquid crystal.     

Convection-Diffusion Lattice Boltzmann Scheme for Irregular Lattices

In this paper, a lattice Boltzmann (LB) scheme for convection diffusion on irregular lattices is presented, which is free of any interpolation or coarse graining step. The scheme is derived using the axioma that the velocity moments of the equilibrium distribution equal those of the Maxwell–Boltzmann distribution. The axioma holds for both Bravais and irregular lattices, implying a single framework for LB schemes for all lattice types. By solving benchmark problems we have shown that the scheme is indeed consistent with convection diffusion. Furthermore, we have compared the performance of the LB schemes with that of finite difference and finite element schemes. The comparison shows that the LB scheme has a similar performance as the one-step second-order Lax–Wendroff scheme: it has little numerical diffusion, but has a slight dispersion error. By changing the relaxation parameter ω the dispersion error can be balanced by a small increase of the numerical diffusion.     

Comparison of Several Spatial Discretizations for the Navier–Stokes Equations

Grid convergence studies for subsonic and transonic flows over airfoils are presented in order to compare the accuracy of several spatial discretizations for the compressible Navier–Stokes equations. The discretizations include the following schemes for the inviscid fluxes: (1) second-order-accurate centered differences with third-order matrix numerical dissipation, (2) the second-order convective upstream split pressure scheme (CUSP), (3) third-order upwind-biased differencing with Roe's flux-difference splitting, and (4) fourth-order centered differences with third-order matrix numerical dissipation. The first three are combined with second-order differencing for the grid metrics and viscous terms. The fourth discretization uses fourth-order differencing for the grid metrics and viscous terms, as well as higher-order approximations near boundaries and for the numerical integration used to calculate forces and moments. The results indicate that the discretization using higher-order approximations for all terms is substantially more accurate than the others, producing less than two percent numerical error in lift and drag components on grids with less than 13,000 nodes for subsonic cases and less than 18,000 nodes for transonic cases. Since the cost per grid node of all of the discretizations studied is comparable, the higher-order discretization produces solutions of a given accuracy much more efficiently than the others.