Homonuclear J coupling effects in volume localized NMR spectroscopy: pitfalls and solutions

Design of phase-modulated broadband refocusing pulses

Broadband linear-phase refocusing pulses were designed with the Shinnar-Le Roux (SLR) transformation and verified experimentally. The design works in several steps: initially, a linear-phase B polynomial is created with the Parks-McClellan/Remez exchange algorithm. The complementary A polynomial required for the SLR transformation is generated with the Hilbert transformation, yielding the minimum-phase response. The phase response of the A polynomial is altered by zero-flipping, which changes the overall pulse shape while retaining its refocusing profile. Optimal pulses in terms of minimal B(1max) and hence broadest bandwidth were found with non-linear optimisation of the zero-flipping pattern. These pulses are generally phase modulated with a time-symmetric amplitude and anti-symmetric phase modulation. In this work, a whole range of pulses were designed to demonstrate the underlying relationships. Five exemplary pulses were implemented into a PRESS sequence and validated by acquiring images of a water-oil phantom and lactate spectra at TE = 144 ms.

A detailed analysis of localized J-difference GABA editing: theoretical and experimental study at 4 T

The problem of low signal-to-noise ratio for gamma-aminobutyric acid (GABA) in vivo is exacerbated by inefficient detection schemes and non-optimal experimental parameters. To analyze the mechanisms for GABA signal loss of a MEGA-PRESS J-difference sequence at 4 T, numerical simulations were performed ranging from ideal to realistic experimental implementation, including volume selection and experimental radio frequency (RF) pulse shapes with a macromolecular minimization scheme. The simulations were found to be in good agreement with phantom and in vivo data from human brain. The overall GABA signal intensity for the simulations with realistic conditions for the MEGA-PRESS difference spectrum was calculated to be almost half of the signal simulated under ideal conditions (~43% signal loss). In contrast, creatine was reduced significantly less then GABA (~19% signal loss). The 'four-compartment' distribution due to J-coupling in the PRESS-based localization was one of the most significant sources of GABA signal loss, in addition to imperfect RF profiles for volume selection and editing. An alternative strategy that reduces signal loss due to the four-compartment distribution is suggested. In summary, a detailed analysis of J-difference editing is provided with estimates of the relative amounts of GABA signal losses due to various mechanisms. The numerical simulations presented in this study should facilitate both implementation of the more efficient acquisition and quantification process of J-coupled systems.

Incorporating homonuclear polarization transfer into PRESS for proton spectral editing: illustration with lactate and glutathione

A proton spectral editing pulse sequence for the detection of metabolites with spin systems that involve weak coupling is presented. The sequence is based on homonuclear polarization transfer incorporated into the standard PRESS (Point RESolved Spectroscopy) sequence, which is a volume-selective double spin echo method, to enable spatial localization. All peaks in the region of interest are initially suppressed whether they are peaks from the target metabolite or from contaminating background. The target signal is then restored by polarization transfer from a proton that has a resonance outside the suppressed region and to which the target spins are weakly coupled. This is achieved by the application of a 90 degrees hard pulse with phase orthogonal to that of the PRESS excitation pulse at the location of the first echo in PRESS and by optimizing the two PRESS timings, TE(1) and TE(2), for most efficient yield. Background signal not coupled to any protons outside the initially saturated region remains suppressed. The advantage of this sequence compared to multiple quantum filters is that signal from singlet peaks outside the suppressed area are preserved and can thus be used as a reference. The efficacy of the sequence was verified experimentally on phantom solutions of lactate and glutathione at 3.0 T. For the AX(3) spin system of lactate, the sequence timings were optimized by product operator calculations whereas for the ABX spin system of the cysteinyl group of glutathione numerical calculations were performed for sequence timing optimization.

Lactate quantification by proton magnetic resonance spectroscopy using a clinical MRI machine: A basic study

The purpose of this study was to establish quantification method of lactate concentration by proton magnetic resonance spectroscopy (MRS) carried out using a conventional 1.5-T MRI machine. We used a lactate phantom with known concentrations (1, 1.5, 3, 6, 12 and 14 mmol/L). As a clinical example, a patient with mitochondrial myopathy, encephalopathy, lactic acidosis and stroke-like episodes (MELAS) was evaluated. Proton MRS was carried out using a clinical 1.5-T super-conducting magnetic resonance whole-body system. Data were acquired by point resolved spectroscopy. A coupling constant of J = 7.35 Hz (2/J = 272 ms) and two long in-phase echo time of 272 ms and 544 ms were used to calculate the T2 relaxation time. The tissue water signal was used as an internal standard to quantify lactate. The correlation coefficient R between the calculated lactate concentrations and the known concentration of lactate was 0.99 with a constant factor of 0.32 (1/3.14). In patients with MELAS, the lactate concentration measured by MRS was 6.2 mmol/kg wet weight, which is similar to the value obtained in previous studies. In the present study, we have established a reliable method for lactate quantification in a phantom study and have shown a sample of clinical case of MELAS.

Direct noninvasive quantification of lactate and high energy phosphates simultaneously in exercising human skeletal muscle by localized magnetic resonance spectroscopy

A novel method based on interleaved localized 31P- and 1H MRS is presented, by which lactate accumulation and the accompanying changes in high energy phosphates in human skeletal muscle can be monitored simultaneously during exercise and recovery. Lactate is quantified using a localized double quantum filter suppressing the abundant lipid signals while taking into account orientation dependent signal modulations. Lactate concentration after ischemic exercise directly quantified by DQF 1H spectroscopy was 24 +/- 3 mmol/L cell water, while 22 +/- 3 mmol/L was expected on the basis of 31P MRS acquired simultaneously. Lactate concentration in a sample of porcine meat was estimated to be 40 +/- 7 mmol/L by means of DQF quantitation, versus 39 +/- 5 mmol/L by biochemical methods. Excellent agreement is shown between lactate concentrations measured noninvasively by 1H MRS, measured biochemically ex vivo, and inferred indirectly in vivo from changes in pH, P(i), and PCr as obtained from 31P MRS data.

Elimination of spatial interference in PRESS-localized editing spectroscopy

Unambiguous detection of gamma-amino butyric acid (GABA) in the human brain is hindered by low concentration and spectral overlap with other metabolites. The popular MEGA-PRESS (PRESS: point-resolved spectroscopic sequence) method allows spectral separation of GABA from other metabolites, but suffers from a significant signal-to-noise ratio (SNR) reduction due to the 4-compartment artifact. An alternative PRESS localization technique (PRESS+4) was investigated and compared to MEGA-PRESS using numerical simulations, phantom, and in vivo experiments. It was shown that while the MEGA-PRESS method suffers significant signal loss ( approximately equal 20% for the difference spectrum), GABA signal intensity in PRESS+4 is reduced by only 2% compared to the nonlocalized condition at 4T. The improved method retains important features of the popular MEGA-PRESS such as additional water suppression and macromolecular elimination as demonstrated in human brain experiments. This method is not limited to GABA J-difference editing, but can be applied in any PRESS-based experiments. It should prove particularly useful at higher field, where the 4-compartment artifact is especially detrimental.

MR spectroscopy and spectroscopic imaging: comparing 3.0 T versus 1.5 T

In vivo magnetic resonance spectroscopy (MR spectroscopy) offers the unique possibility to monitor human brain metabolism in a noninvasive way. At 3.0 T, MR spectroscopy not only profits from higher available signal compared with 1.5 T, but from increased chemical shift dispersion as well. These gains may be exchanged into increased spatial resolution or speed in MR spectroscopic imaging. However, some adverse effects related to the higher field strength, such as increased field inhomogeneities and sequence restrictions caused by safety limitations need to be considered. These require protocol adaptations and technical advances that have not yet fully found their way onto the clinical platform. If neglected, effects such as chemical shift misregistration at higher field strength can lead to wrong localizations or loss of signals of certain metabolites, which can intervene with the diagnostic value of a spectrum. This article tries to give an understanding of the potentials and challenges of MR spectroscopy at the higher field strength of 3.0 T, and to give insight into new techniques that hopefully soon will become available in daily clinical routine to fully exploit all benefits of the higher field strength.

Spatial localization in nuclear magnetic resonance spectroscopy

The ability to select a discrete region within the body for signal acquisition is a fundamental requirement of in vivo NMR spectroscopy. Ideally, it should be possible to tailor the selected volume to coincide exactly with the lesion or tissue of interest, without loss of signal from within this volume or contamination with extraneous signals. Many techniques have been developed over the past 25 years employing a combination of RF coil properties, static magnetic field gradients and pulse sequence design in an attempt to meet these goals. This review presents a comprehensive survey of these techniques, their various advantages and disadvantages, and implications for clinical applications. Particular emphasis is placed on the reliability of the techniques in terms of signal loss, contamination and the effect of nuclear relaxation and J-coupling. The survey includes techniques based on RF coil and pulse design alone, those using static magnetic field gradients, and magnetic resonance spectroscopic imaging. Although there is an emphasis on techniques currently in widespread use (PRESS, STEAM, ISIS and MRSI), the review also includes earlier techniques, in order to provide historical context, and techniques that are promising for future use in clinical and biomedical applications.

Optimized detection of lactate at high fields using inner volume saturation

In localized proton MR spectroscopy ((1)H-MRS) in vivo, the detection of lactate (Lac) is affected by modulation of its resonances due to homonuclear scalar couplings (J). A simple and convenient way to distinguish Lac from lipids is to set the TE to 1/J so that the Lac signal is inverted while other resonances (such as lipid) remain in-phase. However, at high field strengths, such as 3 Tesla or above, the modulation of the Lac signal is complicated by chemical shift effects that cause modulation patterns to vary within different subregions of the localized volume. Under some conditions the Lac signal may even disappear completely. In this note we introduce the concept of inner volume saturation (IVS), which makes use of high bandwidth spatial pulses to remove the signal corresponding to the regions of the localized volume that contribute unwanted modulation patterns. The method is described theoretically and demonstrated experimentally at 3 Tesla in a phantom and a patient with acute stroke. The phantom measurements indicate that virtually 100% of the Lac signal can be recovered using this method. The method should be feasible at magnetic fields above 3 Tesla, and may also be applied to other coupled spin systems in which modulation effects are important.

Non-invasive quantification of lactate by proton MR spectroscopy and its clinical applications

Lactate is an important metabolite in clinical cases indicating the status of metabolic impairment. We applied a clinically relevant simple method for lactate quantification using magnetic resonance spectroscopy (MRS). We used two long in-phase echo time (TE=272,544 ms) to calculate T2 relaxation time and the absolute concentration of lactate. This method was optimized using phantom study and applied to clinical cases. This technique does not require complicated processing and could be applied in daily clinical practice. Moreover, this technique enables lactate quantification in cases (e.g. tumor) where lipid peak is overlapped with the lactate peak at short echo times.

Spin-echo MRS in humans at high field: LASER localisation using FOCI pulses

Significant improvements in spin-echo MRS are possible when voxel localisation is performed using high bandwidth frequency offset corrected inversion (FOCI) pulses as opposed to more conventional lower bandwidth pulses. The reduced chemical shift displacement errors result in a spectrum that more accurately reflects the actual metabolite distribution within any region of interest that is selected graphically on a series of scout images, and can lead to improved metabolite detection in the case of homonuclear J-coupled spins. At 4.7T, FOCI pulses with a 20 kHz bandwidth result in extremely sharp and uniform selection profiles, and negligible contamination from outside of the voxel of interest, for all signals in the 1H spectral range that is normally studied. A 'FOCI' adiabatic half-passage is observed to provide good excitation over the 1H spectral range. Single shot performance with echo-time (TE)48 ms is reported using a four-port drive birdcage head coil. GAMMA simulations show that, for many detectable metabolites at 4.7 T, LASER localisation using FOCI pulses with TE=48 ms results in 1H anti-phase spectral components that are the same order as would be obtained from a symmetric PRESS sequence with TE=32 ms. Timing schemes are proposed to enable good measurement of lactate with very little signal loss arising from chemical shift displacement errors at TE=144 and 288 ms.

Numerical simulation of PRESS localized MR spectroscopy

Numerical simulations of NMR spectra can provide a rapid and convenient method for optimizing acquisition sequence parameters and generating prior spectral information required for parametric spectral analysis. For spatially resolved spectroscopy, spatially dependent variables affect the resultant spectral amplitudes and phases, which must therefore be taken into account in any spectral simulation model. In this study, methods for numerical simulation of spectra obtained using the PRESS localization pulse sequence are examined. A comparison is made between three different simulation models that include different levels of detail regarding the spatial distributions of the excitation functions, and spin evolution during application of the pulses. These methods were evaluated for measurement of spectra from J-coupled spin systems that are of interest for in vivo proton spectroscopy and results compared with experimental data. It is demonstrated that for optimized refocusing pulses it is sufficient to account for chemical shift effects only, although there is some advantage to implementing a more general numerical simulation approach that includes information on RF pulse excitation profiles, which provides sufficient accuracy while maintaining moderate computational requirements and flexibility to handle different spin systems.

Brain γ-aminobutyric acid measurement by proton double-quantum filtering with selectiveJ rewinding

An optimized single-shot proton double-quantum (DQ) filter for the quantification of gamma-aminobutyric acid (GABA) levels in human brain is reported. It is demonstrated that creation of DQ coherences following dual-resonance-selective refocusing gives a theoretical editing efficiency of 50% for the detection of the GABA resonance at 3.01 ppm. The sequence times are optimized with both numerical and experimental analyses of the editing performance, giving an experimental editing efficiency of 42%. It is acknowledged that homocarnosine is partially coedited, leading to a 20% contribution to the edited signal; however, macromolecule contamination is negligible in vivo under these experimental conditions. The GABA concentration in human prefrontal cortex is estimated to be 0.8 +/- 0.1 micromol/g (mean +/- SD, n = 6), with reference to the internal standard creatine at 9 micromol/g.

Multiple spin-echo spectroscopic imaging for rapid quantitative assessment of N-acetylaspartate and lactate in acute stroke

Monitoring the signal levels of lactate (Lac) and N-acetylaspartate (NAA) by chemical shift imaging can provide additional knowledge about tissue damage in acute stroke. Despite the need for this metabolic information, spectroscopic imaging (SI) has not been used routinely for acute stroke patients, mainly due to the long acquisition time required. The presented data demonstrate that the application of a fast multiple spin-echo (MSE) SI sequence can reduce the measurement time to 6 min (four spin echoes per echo train, 32 x 32 matrix). Quantification of Lac and NAA in terms of absolute concentrations (i.e., mmol/l) can be achieved by means of the phantom replacement approach, with correction terms for the longitudinal and transversal relaxation adapted to the multiple spin-echo sequence. In this pilot study of 10 stroke patients (symptom onset < 24 hr), metabolite concentrations obtained from MSE-SI add important information regarding tissue viability that is not provided by other sequences (e.g., diffusion-weighted imaging (DWI) and perfusion-weighted imaging (PWI)). Metabolic changes extended beyond the borders of the apparent diffusion coefficient (ADC) lesion in nine of the 10 patients, showing a rise in Lac concentrations up to 18 mmol/l, while NAA levels sometimes dropped below the detection level. Considerable differences among the patients in terms of the Lac concentrations and the size of the SI-ADC mismatch were observed.

Slice-selective J-coupled coherence transfer using symmetric linear phase pulses: applications to localized GABA spectroscopy

Symmetric, linear phase, slice-selective RF pulses were analyzed theoretically for performing slice-selective coherence transfer. It was shown using numerical simulations of product operators that, when a prefocusing gradient of the same area as that of the refocusing gradient is added, these pulses become slice-selective universal rotator pulses, therefore, capable of performing slice-selective coherence transfer. As an example, a slice-selective universal rotator pulse based on a seven-lobe hamming-filtered sinc pulse was applied to in vivo single-shot simultaneous spectral editing and spatial localization of neurotransmitter GABA in the human brain.

1H magnetic resonance spectroscopy in human hydrocephalus

PURPOSE

To evaluate cerebral metabolism in clinical hydrocephalus with (1)H magnetic resonance spectroscopy (MRS).

MATERIALS AND METHODS

In 24 children and adults with progressive, arrested, or normal pressure hydrocephalus, long-echo time (1)H MR spectra were acquired from periventricular white matter and intraventricular cerebrospinal fluid (CSF). Metabolite ratios, and the presence of lactate, were compared with 38 age-matched controls.

RESULTS

Metabolite ratios of patients were within the 95% confidence interval (CI) of controls. A small lactate resonance was detected in 20% of control and hydrocephalic subjects. Lactate was consistently visible in CSF spectra, though lactate concentrations were normal. The CSF lactate T(2) was long in comparison with the known intracellular metabolite T(2) relaxation times. In three neonates with hydrocephalus and spina bifida, 3-hydroxybutyrate was detected in CSF in vivo.

CONCLUSION

Within the limits of the present methods, (1)H MRS could not detect cerebral metabolic abnormalities in human hydrocephalus and provided no additional diagnostic information. The long T(2) of lactate in CSF explains its high visibility. Hence, the detection of lactate in spectra acquired from voxels that contain CSF does not necessarily imply cerebral ischemia.

Proton magnetic resonance spectroscopy in the brain: report of AAPM MR Task Group #9

AAPM Magnetic Resonance Task Group #9 on proton magnetic resonance spectroscopy (MRS) in the brain was formed to provide a reference document for acquiring and processing proton (1H) MRS acquired from brain tissue. MRS is becoming a common adjunct to magnetic resonance imaging (MRI), especially for the differential diagnosis of tumors in the brain. Even though MR imaging is an offshoot of MR spectroscopy, clinical medical physicists familiar with MRI may not be familiar with many of the common practical issues regarding MRS. Numerous research laboratories perform in vivo MRS on other magnetic nuclei, such as 31P, 13C, and 19F. However, most commercial MR scanners are generally only capable of spectroscopy using the signals from protons. Therefore this paper is of limited scope, giving an overview of technical issues that are important to clinical proton MRS, discussing some common clinical MRS problems, and suggesting how they might be resolved. Some fundamental issues covered in this paper are common to many forms of magnetic resonance spectroscopy and are written as an introduction for the reader to these methods. These topics include shimming, eddy currents, spatial localization, solvent saturation, and post-processing methods. The document also provides an extensive review of the literature to guide the practicing medical physicist to resources that may be useful for dealing with issues not covered in the current article.

Quantitative evaluation of the lactate signal loss and its spatial dependence in press localized (1)H NMR spectroscopy

Localized (1)H NMR spectroscopy using the 90 degrees -t(1)-180 degrees -t(1)+t(2)-180 degrees -t(2)-Acq. PRESS sequence can lead to a signal loss for the lactate doublet compared with signals from uncoupled nuclei which is dependent on the choice of t(1) and t(2). The most striking signal loss of up to 78% of the total signal occurs with the symmetrical PRESS sequence (t(1)=t(2)) at an echo time of 2/J (approximately 290 ms). Calculations have shown that this signal loss is related to the pulse angle distributions produced by the two refocusing pulses which leads to the creation of single quantum polarization transfer (PT) as well as to not directly observable states (NDOS) of the lactate AX(3) spin system: zero- and multiple-quantum coherences, and longitudinal spin orders. In addition, the chemical shift dependent voxel displacement (VOD) leads to further signal loss. By calculating the density operator for various of the echo times TE=n/J, n=1, 2, 3,..., we calculated quantitatively the contributions of these effects to the signal loss as well as their spatial distribution. A maximum signal loss of 75% can be expected from theory for the symmetrical PRESS sequence and TE=2/J for Hamming filtered sinc pulses, whereby 47% are due to the creation of NDOS and up to 28% arise from PT. Taking also the VOD effect into account (2 mT/m slice selection gradients, 20-mm slices) leads to 54% signal loss from NDOS and up to 24% from PT, leading to a maximum signal loss of 78%. Using RE-BURP pulses with their more rectangular pulse angle distributions reduces the maximum signal loss to 44%. Experiments at 1.5 T using a lactate solution demonstrated a maximum lactate signal loss for sinc pulses of 82% (52% NDOS, 30% PT) at TE=290 ms using the symmetrical PRESS sequence. The great signal loss and its spatial distribution is of importance for investigations using a symmetrical PRESS sequence at TE=2/J.

The effects of slice-selective excitation/refocusing in localized spectral editing with gradient-selected double-quantum coherence transfer

Spectral editing using gradient-selected double-quantum filtering (DQF) with PRESS localization has been used for selective observation of metabolites in vivo. In previous studies using localized DQF sequences, it is generally assumed that the slice-selective pulses used in the sequence have no roles in coherence transfer, and do not interfere with DQF. To validate this assumption, the effects of slice-selective excitation/refocusing on DQF were investigated in DQF lactate editing sequences combined with PRESS localization. Contrary to the previous assumption, the results show that, due to chemical shift displacement artifact and J coupling, slice selection in DQF does interfere with coherence transfer, affecting both the accuracy of spatial localization and the detection sensitivity adversely. In the case of lactate editing, the effects of this interference can be accounted for simply by adjusting the strength of the slice-selection gradients and by using narrowband slice-selective refocusing pulses.

Two-voxel localization sequence for in vivo two-dimensional homonuclear correlation spectroscopy

The combination of localized 2D 1H MR correlation spectroscopy and Hadamard encoding allows the simultaneous acquisition of multiple volumes of interest without an increase in the experimental duration, compared to single-voxel acquisition. In the present study, 2D correlation spectra were acquired simultaneously within 20 to 40 min in two voxels located in each hemisphere of the rat brain. An intervoxel distance of 20% of the voxel size was sufficient to limit spatial contamination. The following cerebral metabolites gave detectable crosspeaks: N-acetylaspartate, the glutamate/glutamine pool, aspartate, phosphoethanolamine, glucose, glutathione, taurine, myo-inositols, lactate, threonine, gamma-aminobutyric acid, and alanine. Most of the metabolites were measured without contamination of other resonances.

Proton NMR chemical shifts and coupling constants for brain metabolites

Proton NMR chemical shift and J-coupling values are presented for 35 metabolites that can be detected by in vivo or in vitro NMR studies of mammalian brain. Measurements were obtained using high-field NMR spectra of metabolites in solution, under conditions typical for normal physiological temperature and pH. This information is presented with an accuracy that is suitable for computer simulation of metabolite spectra to be used as basis functions of a parametric spectral analysis procedure. This procedure is verified by the analysis of a rat brain extract spectrum, using the measured spectral parameters. In addition, the metabolite structures and example spectra are presented, and clinical applications and MR spectroscopic measurements of these metabolites are reviewed.

Sources of variability in the response of coupled spins to the PRESS sequence and their potential impact on metabolite quantification

Using a numerical method of solving the equation of motion of the density matrix, an evaluation is presented of the sources of the marked variability in the response to the point resolved spectroscopy (PRESS) pulse sequence of coupled proton spin systems. The consequences of an inappropriate 180 degrees pulse design and of the limitations on radiofrequency power are demonstrated for a weakly coupled example, lactate. The dominating role of strong coupling, which is present in most brain metabolites, is demonstrated for glutamate, in which 160 terms in the density operator were tracked to monitor the gross changes in lineshape and signal intensity as a function of the two echo times. The predictions of the numerical solutions were confirmed by experiments on phantoms of aqueous metabolite solutions.

Lactate detection at 3T: compensating J coupling effects with BASING

Detection of lactate by in vivo 1H magnetic resonance spectroscopy may provide a means of identifying regions of metabolic stress in brain and other human tissue, potentially identifying regional ischemia in stroke or necrosis in tumors. At higher field strengths (3 and 4 T), which have recently become available for whole-body human studies, the chemical shift difference between the doublet from the methyl protons and the quartet from the methine proton becomes comparable to the available radiofrequency (RF) pulse bandwidth. In this case "anomalous" J modulation occurs in PRESS and STEAM because the coupling partner of the observed resonance may or may not be refocused by the RF pulses depending on the position of the molecule within the voxel and the size of the chemical shift misregistration artifact. These anomalies lead to signal cancellation for echo times near odd multiples of 1/J (often used to highlight the inverted lactate doublet against nearby lipid peaks) in single voxel studies, and spatial variation of the doublet lineshape in chemical shift imaging studies, producing erroneous determination of relative lactate concentrations. While increasing the band-width of the RF pulses can reduce this effect by reducing the signal cancellation, some cancellation will always remain. A means of eliminating this effect using BASING/ MEGA (Mescher M et al. Solvent suppression using selective echo dephasing J Magn Reson A 1996;123:226-229; Star-Lack J et al. Improved water and lipid suppression for 3D PRESS CSI using RF band selective inversion with gradient dephasing (BASING). Magn Reson Med 1997;38: 311-321) water suppression pulses will be described, along with some of its limitations.

Automated spectral analysis I: formation of a priori information by spectral simulation

A spectral simulation method is described for generating a priori information for use in parametric spectral analysis. The method makes use of GAMMA (S. A. Smith, T. O. Levante, B. H. Meier, R. R. Ernst, J. Magn. Reson., 106A, 75-105, 1994), a programming environment that facilitates simulation of magnetic resonance phenomena. The input parameters consist of the chemical shifts and scalar spin-coupling constants for the compounds to be analyzed, the acquisition pulse sequence, and the field strength used. The resultant spectral information consists of the relative amplitude, frequency, and phase of all resonances, which are stored in a spectral database. This procedure can be rapidly and conveniently modified to reflect different acquisition parameters and data analysis requirements.