QENS DATA ANALYSIS Workbook, 1st Edition

 Hilary James
 8 months ago
 Views:
Transcription
1 Workbook, 1st Edition This workbook is designed to instruct those neutron users new to the QuasiElastic Neutron Scattering technique about basic QENS data reduction. The learning approach adopted is guided analysis of real data to illustrate some of the challenges posed by actual problems and datasets. The worked examples chosen also illustrate the kind of information one can extract from a QENS experiment once data reduction is complete and analysis begins. The reduction and analysis process is specifically underpinned by the Mantid package. Mark Telling NBIAW, November 2014
2 Introduction This workbook is designed to instruct those neutron users new to the QuasiElastic Neutron Scattering technique about basic QENS data reduction. The learning approach adopted is guided analysis of real data to illustrate the challenges posed by actual problems and datasets. The worked examples chosen (translational diffusion of water molecules and side group rotation in peptides) also illustrate the kind of information one can extract from a QENS experiment once data reduction is complete and analysis begins. The reduction and analysis process is specifically underpinned by the Mantid package. Comments and suggestions about this workbook are most welcome and should be sent to: All the data files referred to in this workbook can be downloaded as a.zip file from: under Related Links The Mantid project provides a framework that supports highperformance computing and visualization of scientific data. Mantid has been created to manipulate and analyse neutron and muon data, but could be applied to many other techniques. The framework is open source and supported on multiple target platforms (Windows, Linux, Mac). Tailor made data reduction and analysis procedures can be formulated using Python scripting. However, this tutorial focuses on the reduction and analysis of QENS data collected from any neutron instrument via a Graphical User Interface (GUI). This interface has been specifically designed to guide the user through the data reduction and analysis procedure. NB: All procedures have been tested using the Windows 7 installation and Mantid 3.2 Mark Telling NBIAW, th November
3 Mantid can be downloaded from, download.mantidproject.org For reliability choose the Latest Stable Build. However, as the interface is developing rapidly, the Nightly Development Build may contain additional features and bug fixes that are useful for data analysis. The page, contains a wealth of information and supplementary information regarding the Mantid project and its functionality. Mark Telling NBIAW, th November
4 Getting Started To launch the main Mantid interface type: MantidPlot (Unix) or click on the Windows desktop icon The worked examples presented in this workbook will use two dedicated GUI s. These are launched by following: Interfaces > Indirect > Convert To Energy or Interfaces > Indirect > Indirect Data Analysis from the tool bar at the top of the main MantidPlot window. Before starting you should ensure that Mantid is configured for your needs. In particular you should add paths to directories which contain the worked example data. To do this, click on File > Manage User Directories. The window below will appear. Data can be added as a path by browsing or typing the directory name in the field below Data Search Directories and then clicking Add Directory. Once a path has been added, data stored in that location can be accessed simply by entering the run number. The Manage User Directories window Mark Telling NBIAW, th November
5 Basic Mantid Data Manipulation This section explains: the concept of workspaces how to export data how to use QtiPlot fitting how to overlay data I. Workspaces Mantid utilises workspaces to hold data. These come in several forms, but the most common is the Matrix Workspace which represents XYE data for one, or several, spectra. The data contained in a workspace can viewed as a matrix or as a table, and graphed in many ways including line graphs, contour plots and 3D graphs. Interaction with workspaces is typically through a Graphical User interface. Matrix Workspaces are typically created by executing one of Mantid's 'Load' algorithms or are the output of algorithms which take a Matrix Workspace as input. In addition to data, workspaces hold a workspace history which lists the algorithms that have been applied to the data. Two types of workspace are listed in the Workspace List pane in the figure below NB: In the following discussion workspace IRIS has been created by loading raw data file via the LOAD button at the top of the Workspaces window and IRS57252_graphite002_red has been created using the Interfaces > Indirect > Convert To Energy GUI. These two methods of loading data will be discussed later Mark Telling NBIAW, th November
6 The workspace list window a) Workspace: IRIS Workspace, IRIS , holds within it raw neutron counts / s for each individual IRIS detector. To examine the data collected in a single detector, right click on workspace IRIS and select Plot Spectrum. The following input fields appear. Raw data plotting options Mark Telling NBIAW, th November
7 On IRIS, as way of example, there are 114 detectors hence ID numbers: (1,2 = neutron monitors, 353 = graphite analyser detectors, = mica analyser detectors, = diffraction detectors). Enter a detector (ID) number and click OK to plot the associated raw data. The raw counts / s for detector 5 are shown below. Alternatively, entering '5,10,20' plots spectra 5,10 and 20 on one graph Raw data (counts / s) collected in IRIS detector no. 5 (left) and in detectors 5,10,20 (right) NB: Click on the small arrows beside the workspace name to expand the entry and see specific workspace information. b) Workspace: IRS57252_graphite002_red Contains corrected or _REDuced data converted from counts/s to counts/energy transfer. Workspace IRS57252_graphite002_red contains only 10 spectra/histograms since the raw data in workspace IRIS has been reduced and corrected such that a) only detectors viewing the graphite analyser are included (i.e. 3 to 53) in the analysis and b) the raw counts in every 5 detectors have been added together to give 10 output spectra. Mark Telling NBIAW, th November
8 II. Exporting Data To export the data contained within any listed workspace click on the Algorithms tab at the bottom of the workspace list pane. The options shown below should appear. Select SaveAscii from the dropdown menu (beside Execute) and click Execute. The SaveAscii Input Dialog box should appear. Select a directory (for the written data file) and specify a file name. The workspace to be exported can be selected from the dropdown list next to InputWorkspace. The user can also, for example, define which workspace spectra i.e. indices (Min and Max) to export, choose the type of data separator used in the file (CSV is usually a reliable option), add comments and decide whether to WriteXError. The Algorithms options Mark Telling NBIAW, th November
9 The SaveAscii dialog box III. QtiPlot Enabling QtiPlot allows data sets and workspaces to be analysed using a suite of inbuilt Origin fitting algorithms. To enable QtiPlot follow View > Preferences > Fitting. The following window should appear. Select Enable QtiPlot Fitting The Choose Default Settings window Mark Telling NBIAW, th November
10 When QtiPlot fitting has been enabled, and a data file has been plotted, Analysis appears in the list of options above the main MantidPlot window toolbar. Tool bar options at the top of MantidPlot. Analysis appears when QtiPlot is enabled The fitting options found under Analysis are shown below The QtiPlot analysis options Mark Telling NBIAW, th November
11 IV Overlaying Data To overlay data, select the plot window in which you wish to compare data. Right click on the plot window and select Add/Remove Curve. A list of available data sets will appear. As an example, in the figure below we have the option of adding the data contained in the workspaces IRIS_57252_H20_290K_002sp10 and IRIS_57252_H20_290K_002sp20 to the plot window of IRIS_57252_H20_290K_ To do so, first select the two spectra in the left hand pane and then click the right pointing arrow. Reversing this process removes data from a plot window. Add / Remove Curves options window The plot style (line, scatter, line + scatter etc) of the added data set can be selected using the New Curves Style drop down menu at the top of the Add/Remove Curve window. Mark Telling NBIAW, th November
12 Example One: translational diffusion Aims To reduce and correct (using Mantid s inbuilt analytical and numerical integration algorithms) QENS data collected from water (H K) To perform line width analysis and determine the diffusion coefficient of H 2 0 at 290K To compare the result with literature To consider the limitations imposed by the experimental data on the analysis Data files File Name Sample Proton Current IRIS_57242_H20_280K_002 Water 280 K uamps IRIS_57246_H20_290K_002 Water 290 K uamps IRIS_57249_H20_300K_002 Water 300 K uamps IRIS_57252_H20_310K_002 Water 310 K uamps IRIS_57081_Van_RT_002 Van Cylinder uamps IRIS_56068_Empty_Annular_RT_002 Empty Annular Can uamps Mark Telling NBIAW, th November
13 Task One: optimal sample thickness To minimise multiple scattering effects (and hence complex correction procedures which will not be discussed here) a sample s thickness should be limited such that it scatters only % of the incident beam. For more info about multiple scattering see: Bée M., Quasi Elastic Neutron Scattering, 1988, Adam Hilger, Chapter 4, p 107) Using the BeerLambert relation below, calculate the optimal thickness of a water sample I trans = I incident exp(nt) n = number density = total scattering cross section / formula unit t = sample thickness 1. Calculate the thickness of the water sample Task Two: Detector efficiency correction file Different detectors will have slightly different efficiencies. The result is an observable fluctuation in measured neutron intensity across a detector bank. Ideally, a wholly elastic scattering sample should have the same measured intensity in all detectors. To correct for such discrepancy, data is collected from a sample that only scatters neutrons elastically; usually a thin walled vanadium cylinder ( ~ 10% scatterer). This measurement should be performed before the actual experiment begins and then used to create a detector calibration file that is applied to all subsequent measurements during analysis. Mark Telling NBIAW, th November
14 To create a detector calibration file using Mantid: 1. Launch the GUI: Interfaces > Indirect > Convert To Energy 2. Select the correct spectrometer 3. Click on the tab Calibration 4. Browse for the Vanadium Cylinder data file and Open 5. Select Create RES File (default limits are loaded but ensure Start/End < 0.5) 6. The calibration algorithm works by subtracting a background level and integrating over the elastic peak for each of detector. The resulting integrals, or areas, are then compared and a correction factor generated. 7. To see the resulting detector calibration and RES files instrument Run_No _AnalyserReflection_calib instrument Run_No _AnalyserReflection_res check Plot result and Save Result and then Run Calibration The IRIS57081_graphite002_calib file Mark Telling NBIAW, th November
15 Q. How might the calibration plot differ if you created a detector calibration file using data collected from a 2mm thick flat vanadium plate oriented at 90 o to the incident beam? A note on Create Res File : Choosing this option creates a workspace labelled Instrument Run_No _AnalyserReflection_res which will be used later for QENS line width analysis. The algorithm takes the experimental vanadium calibration data in each detector, efficiency corrects the raw counts and sums all the spectra. The result / output is a single high statistic instrument resolution file. Q. What are the possible assumptions / limitations imposed by creating a RES file in this manner? Task Three: Applying the detector efficiency correction file to the data 1. Select tab Energy Transfer 2. Load data file IRIS_57246_H20_290K_002 in Run Files input field 3. Load the detector calibration file and check Use Calib File 4. In Mapping select individual the mapping option allows you to define how you group the detectors. The option individual treats all 51 IRIS spectra individually (surprise surprise!) 5. Select contour from the Plot Output options 6. Click Run Energy Transfer and the contour plot below should appear (a new workspace with detector efficiency corrected H20 290K data in it should appear with the name ****_red ) 7. Reduced the empty sample can data (56068) in the same way Mark Telling NBIAW, th November
16 Efficiency corrected data : individual spectra : IRIS57246_graphite002_red Tip: to visualise a workspace in 3D double click on the workspace name and the option 3D plot will appear in the main taskbar Q. How could you test that the detector efficiency corrections were being applied correctly? Try it! Task Four: Absorption corrections and empty can subtraction The main correction to be applied to neutron scattering data is that for neutron absorption both in the sample and its container, when present. For flat plate geometry, theses corrections can be analytical and have been discussed, for example, by Carlile [C J Carlile, Rutherford Laboratory report, RL (1974)]. The situation for cylindrical geometry is more complex and requires numerical integration. These techniques are well known, used in liquid and amorphous diffraction and are described in the ATLAS manual [A K Soper, W S Howells & A C Hannon, RAL Report RAL (1989) H H Paalman & C J Pings, J Appl Phys Mark Telling NBIAW, th November
17 (1962)]. The routines used here have been developed from the corrections programs in the ATLAS suite and take into account the wavelength variation of both the absorption and the scattering crosssections for the inelastic flight paths.the application of theoretical absorption corrections to experimental data is a 2stage process. Stage 1: Generate correction files based on user input sample and samplecan geometry and composition information: 1. Launch the Indirect Data Analysis GUI 2. Click on Calculate Corrections 3. Input type File and select detector efficiency corrected H 2 0 file 4. Select Use Can this option will incorporate the effect of scattering from the aluminium sample can. Deselect if you just want to treat the sample 5. The water data was collected using an aluminium annular can. Enter the 3 radii and step size (see For Reference below). NB: The calculation cuts the cylinder into several annuli where stepsize is the radial increment. There are a minimum number of annuli, n, required for reliable results (20) where n is determined from n = (r2r1)/step. If n < 20 then an error will occur. The smaller the step the better the result but computing time increases. 6. Calculate number density and absorption / scattering cross sections for sample and sample can. Enter calculated parameters. Number density in atoms/å 3. NB In is possible to generate cross section information by entering the samples chemical formula. For example, H2O is entered: H2O but check the values generated in the Result Log 7. Select Plot Output wavelength, tick Save Result and click Run 8. Workspace ****_cyl_abs should be created. This work space contains neutron scattering and absorption data for both sample and samplecan as a function of wavelength. Mark Telling NBIAW, th November
18 Stage 2: Apply correction file to experimental data and subtract scattering intensity from empty samplecell 1. Select tab Apply Corrections 2. Tick all relevant check boxes: a. Tick Use Corrections to apply the ****_cyl_abs correction file to the data b. Tick Use Can to subtract empty samplecan scattering intensity c. Tick Scale Can By if exact number densities / neutron scattering properties are unknown. The Scale value will adjust the intensity of the empty cell data to avoid over correction (which may result in an inverted peak at E = 0) 3. Enter empty samplecan file or select workspace 4. Select Plot Contributions to see relative intensities and Plot Output > contour 5. Select Save Result and click Run 6. Two new workspaces are created. Depending on the level of correction performed the workspaces will have different labels. For fully corrected data (i.e absorption and empty cell subtraction) the extensions will be: ****_correct_runno_rqw, ****_correct_runno_red ****_correct_runno_result NB: For just empty can subtracted data correct will be replaced by subtract For Reference: Aluminium has a Face Centre Cubic crystal structure with a lattice constant of Å and a density of 2.7 gcm 3 Mark Telling NBIAW, th November
19 Total scattering and absorption cross sections can be found at: For annular geometry enter values of r1, r2 and r3 where: r1 = inner sample radius r2 = outer sample radius r3 = outer sample radius + (2 x can wall thickness) For the sample can use: Inner sample radius = 1cm Sample thickness = use calculated value from Task One Sample can wall thickness = 0.06cm Beam width = 2 cm Task Five: Linewidth analysis in energy There are several approaches to QENS line width analysis; the method chosen is really a matter of personal preference and can include basic least squares fitting or Bayesian probability (D.Sivia, Physica B, 202, 332, 1994). Here we will use the former, and the ConvFit option, to convolve N lorentzian functions with the instrument resolution file (and a background if necessary) to get the best fit to the experimental data. To perform line width analysis: Mark Telling NBIAW, th November
20 1. In the Indirect Data Analysis GUI click on ConvFit 2. Select your corrected water data file 3. Select the _res file create in Task Two 4. Select Fit Type One Lorentzian Q. Do you need a background term? Include one for the purpose of this demonstration 5. Adjust your starting parameters  fitting range (blue lines), background level etc. If you want to Fix a parameter for all fits, right click on the Property name and select Fix 6. Manually adjusting the red lines allows you to guess a starting line width and amplitude for your fit. To activate them, and visualise the effect of broadening/narrowing the lorentzian, click Plot Guess you can also adjust the parameters of the lorentzian by using the Lorentzian 1 or 2 Property box 7. Select Plot Output (Single Run) and click Run this will plot the fit for just one spectrum in a new plot window with the residual; the spectrum shown being determined by the value in the first Spectra Range input field. Suggestion: choose spectrum If the resulting fit is poor adjust your fit parameters accordingly 9. When satisfied with the fit of one spectrum select the complete spectral range for fitting, Plot Output : FWHM, check Save Result and click Run Sequential Fit 10. A plot of line width vs. Q will appear. The plot data (data, fits, result, residual) for a sequential 1L fit will be written to workspace: *****_graphite002_conv_1lfitl_ start _to_ end _Workspace 11. All fit parameters will be written to table *****_graphite002_conv_1lfitl_0_to_50_parameters Mark Telling NBIAW, th November
21 12. Graphs of line width and amplitude will be written to *****_graphite002_conv_1lfitl_0_to_50_results Q. From a data analysis point of view, what are the limitations of the data that has been collected? Q. From an experimental point of view, what could you do to improve the results of your analysis? Task Six: Fitting By modelling the Q dependence of the QENS line width, diffusion coefficients, residence times and jump lengths associated with diffusive motions of molecules can be extracted. Several models are proposed in the literature and each considers deviations from the continuous diffusion model, or Fick s Law, for specific environments when the mechanism of the diffusion needs to be considered at smaller and smaller length scales i.e. larger and larger Q. For a comprehensive overview see Bée M., QuasiElastic Neutron Scattering, 1988, Adam Hilger, Chapter 5, p 148. As way of examples, Model Function Fick s law ( hwhm)( Q) DQ l Q 6 The ChudleyElliot Jump diffusion model (1961) 1 sin( Ql) ( hwhm)( Q) 1 Ql Mark Telling NBIAW, th November
22 Diffusive modes in water: i) Singwi and Sjölander (1960) exp( Q D o R l ( hwhm)( Q) 1 2 o 1 Q D o 2 1 ) ii) Teixeira (1985) 2 2 DQ l ( hwhm)( Q) ; D 2 1 DQ 6 o o av Use the Teixeira model, fit your line width data and extract root mean jump length, L = <l 2 > 0.5, the residence time, o, and ultimately the translational diffusion constant, D. To code and fit a function in Mantid: 1. Click on the plot of line width vs. Q 2. Click on Analysis in the main Mantid window tool bar 3. Select Fit Wizard 4. Select Category > User defined 5. Enter a name for your new function 6. Enter the equation into the blank input field. For example, y=mx+c might be coded: Name: Linear Fit Parameters: m, c Entry on equation field: (m*x)+c 7. Save the file and click on Fit to launch the Fit Wizard dialogue box 8. Click Fit Mark Telling NBIAW, th November
23 Q. How do your values compare with values reported in the literature? Q. A single lorentzian (and background) fit to S(Q,) is a greatly oversimplified description of the water diffusion problem. Why? What does such a simple model neglect? Suggestions: Further Analysis Process the 280, 300 and 310K data sets and extract mean squared jump lengths and residence times. Do the results follow the trends presented in the Teixeira paper? Create a plot to compare of line width vs. Q vs. temperature and associated fits. Mark Telling NBIAW, th November
24 Example Two: localised diffusion While QENS line width analysis can tell us a lot about the motion, or diffusive properties, of molecules the Q dependence of the intensity of only those neutrons scattered elastically (i.e. the intensity at E = 0) also contains a wealth of information; from the geometry of a localised motion to transition temperatures. For more detailed information about the science driving the following worked analysis please see: pubs.rsc.org/en/content/articlelanding/2011/sm/c1sm05603d Aims To reduce and correct elastic incoherent scattering data collected as a function of temperature from the lyophilised protein, apoferritin. To identify transition temperatures via analysis of the elastic scattering intensity To investigate the effect of hydration upon the protein via analysis of the mean squared displacement To investigate the Qdependence of the line width associated with a localised mode To determine the elastic incoherent structure factor (EISF) associated with the dynamic process Mark Telling NBIAW, th November
25 Data files File Name Temp Sample Start Temp / Increment Proton Current IRS K Dry Apo IRS38114 to IRS38156 Dry Apo Start Temp 15K T=5K 30.1 IRS K Dry Apo IRS K Empty Flat Plate Cell OSI K Hydrated Apo h= OSI65562 to OSI65608 Hydrated Apo h=0.14 Start Temp 10 K T=10K 65 OSI65610 to OSI65636 Hydrated Apo h=0.14 Start Temp 245 K T=5K 65 OSI K Empty Flat Plate Cell 65 OSI K Hydrated Apo h= OSI65439 to OSI65495 Hydrated Apo h=0.31 Start Temp 20 K T=10K 65.1 OSI K Empty Flat Plate Cell Task One: Detector efficiency correction file 1. Taking the DRY (aka LYOPHILISED) PROTEIN data only! Follow the procedure set out in Example 1: Task 2. However, this time create a detector efficiency file from the sample itself measured at low temperature. NB: This approach should ONLY be used if the temperature of the sample is low enough that you are sure that ALL scattering is elastic! The advantage of this approach, however, is that Mark Telling NBIAW, th November
26 you have a direct measure of absorption effects from both sample and sample can as well as correction for the orientation of the sample can should a flat plate cell be used. Q. What might be the disadvantages of this method when creating a RES file from experimental flat samplecan data? Task Two: Applying the detector efficiency correction file to the data 1. Create a suitable detector mapping file (****.map) to process individual detectors BUT such that the dead detectors seen in the efficiency plot are removed from all subsequent analysis. A mapping file should have the following structure: 3 ; number of spectra to be created 1 ; spectrum 1 13 ; number of detectors to be added ; detector nos 2 ; spectrum The mapping file above will create 3 spectra. These spectra will be the sum of counts in detectors 315, and Mark Telling NBIAW, th November
27 3. Check that the mapping file is applied correctly 4. Apply the detector efficiency correction and mapping file to all the lyophilised data in the Data Files list above so you have a _red file for > and also for the empty cell Task Three: Empty sample can subtraction 1. Follow the procedure set out in Stage 2, Example 1: Task 4 to subtract the empty cell from all detector efficiency corrected data files Task Four: Analysis of the elastic scattering intensity Since energy must be conserved during the scattering process, any onset of QENS broadening must result in an associated decrease in elastic scattering intensity. Monitoring the elastic scattering intensity, I(Q,T, 0), as a function of temperature and momentum transfer therefore allows transition temperatures to be identified. To analyse the intensity of the elastic peak from our lyophilised protein: 1. Select the Indirect Data Analysis GUI 2. Click on the Elwin tab 3. Select the _subtract_ data files which correspond to temperatures 15K to 290K 4. The blue lines show the range over which integration of the elastic line will be performed Q. Does this range seem suitable? Remember ONLY neutrons scattered elastically should be included in the analysis Mark Telling NBIAW, th November
28 7. Check Save Result and click Run 8. Three workspaces are written: Inst First Run No _AnalyserReflection_to_ Last Run No _eq1 = intensity vs. Q vs. temp Inst First Run No _ AnalyserReflection _to_ Last Run No _eq2 = ln(intensity) vs. Q 2 vs. temp (for mean squared displacement analysis) Inst First Run No _ AnalyserReflection _to_ Last Run No _elf = intensity vs. Q vs. temp NB: if only one data file is processed using ELWIN the output workspaces will be named inst Run No _ AnalyserReflection _eq'n' and inst Run No _ AnalyserReflection _elf NB: depending on when a data file was collected, the xaxis in the.elf workspace will be either the last 3 digits of the run number (i.e > 113) OR the temperature associated with the measurement. Sample temps will be used if the raw data file has a temperature log associated with it. For example, when loading 38113, Elwin looks for SE Log Name: sample. However, this log name does not exist in the raw data log information (right click on the irs38113_graphite002_red and open Sample Logs) In contrast, irs27252_graphite002_red does contain a temperature log Name called sample 9. Select the _elf file > right click > Plot Spectrum With Errors > enter 3 spectrum values of your choice (suggestion: 10,30,40) and click OK. Alternatively, plot the _elf file in 3D Q. What do you notice? Is the elastic intensity constant as a function of run number? Q. Is the data too noisy? If so, how could you reduce the noise? Mark Telling NBIAW, th November
29 Task Five: Mean square displacement, <r(t) 2 >, analysis The mean square displacement, <r(t) 2 >, of atoms in a material can be extracted from elastic neutron scattering intensities, I(Q,T, 0), by fitting, S inc ( Q, T, 0) I( Q, T, 0) exp( Q r( T) ) (1) 3 Q Fit exp 2 r 3 2 Strictly speaking, this form is valid for harmonic oscillations or equivalent atoms. Any deviation from harmonic behaviour is noted as the activation of anharmonic behaviour. Experimentally, harmonic to anharmonic behaviour is usually identified by an inflexion in the temperature dependence of <r(t) 2 >. Equation (1) is only valid in the lowq regime. However, <r(t) 2 > still proves itself to be a useful tool for relative changes in <r(t) 2 > as a function of, for example, hydration level, temperature and / or sample. Mark Telling NBIAW, th November
30 To extract <r(t) 2 > from the elastic scattering created in Task Four: 1. Select tab MSD Fit 2. Browse for and load the _eq2 file 3. Select one of the higher temperature spectra to display using Plot Spectrum (suggestion : 30) 4. Eq (1) above will be fitted to data points delineated by the blue lines. Q. Does the fit range in Q look suitable? If not, adjust the range by modifying StartX and EndX in the MSD tab Property field or drag the blue lines on the plot 5. Check Plot Result and Save Result 6. Click on Run to fit the data on screen OR Run Sequential to fit all temperatures 7. The fit of Eqn (1) to the data on screen is shown 8. Three workspaces are created: Inst First Run No _graphite002_to_ Last Run No _msd_workspaces contains individual workspaces with raw data, associated fits and residuals Inst First Run No _graphite002_to_ Last Run No _msd contains the initial amplitude (A0) and gradient (A1) vs run no/temp Inst First Run No _graphite002_to_ Last Run No _msd_parameters table of A0 and A1 vs run number / temperature 16. To see <r 2 > vs. temperature  right click on Inst First Run No _graphite002_to_ Last Run No _msd and select Plot Spectrum With Errors. Right click on the plot window that appear, select Add/Remove Curve and remove workspace _a0 Mark Telling NBIAW, th November
31 Q. Is there a change in the slope of <r(t) 2 >? If so, does it broadly correspond with the temperature at which the elastic intensity starts to decrease? NB: if no temperature information is included in the output files then the temps will have to be entered manually by viewing the table in workspace Inst First Run No _graphite002_to_ Last Run No _msd_parametersand changing the xaxis values Task Six: <r(t) 2 > and the effect of hydration 1. Now repeat the data reduction, ELWIN and MSD Fit analysis detailed above using the hydrated protein data sets. Plot the resulting <r(t) 2 > from the hydrated materials on top of the lyophilised protein response. Q. What do you notice? Q. Should you actually compare the three experimental data sets? Q. If so, what is this result telling you about the protein? Task Seven: Line width analysis in time Experimentally, the measured scattering function, S meas inc(q,), is a convolution of S inc (Q,) and the resolution function of the neutron instrument, R(Q,). For spectrometers operating in Q space, S meas inc(q,w) = S inc (Q,) x R(Q,). In its simplest form the instrument resolution approximates to a Gaussian or Lorentzian function of finite width, res (usually quoted as full width at half maximum). As we saw in Example 1, using either a measured or theoretical R(Q,), least squares fitting or Bayesian analysis routines can be used to isolate the intensities and widths of the spectral contributions to S inc (Q,). Here, however, we will Mark Telling NBIAW, th November
32 adopt an analysis method which relies on the Fourier transform of the measured scattering function. The merits of fitting in the time regime are discussed in V. Arrighi, J. S. Higgins, A. N. Burgess and W. S. Howells, Macromolecules, 1995, 28, Using Fast Fourier Transform (FFT) methods, the measured QENS and resolution spectra are converted to the time domain. Deconvolution of R(Q,) and S meas inc(q,), is achieved by simply dividing the Fourier response of the sample by that of the resolution. The result is the timedependent intermediate scattering function, I(Q,t). In the simplest case, a single relaxation process will manifest itself in the time domain as a simple exponential, I(Q,t) = A o (Q) + [1  A o (Q)]exp((t/). Here, is the relaxation time and A o (Q) is the Elastic Incoherent Structure Factor (EISF). A system that exhibits a distribution of relaxation rates, however, may be better described using the Kohlrausch Williams Watt (KWW), or stretched exponential, form i.e. I(Q,t) = A o (Q) + [1  A o (Q)]exp((t/ KWW ) ) It should be noted that here KWW is an effective relaxation time which is dependent upon both T and, or more correctly the temperature dependence of the spectral shape of the distribution. Nonexponential behaviour is immediately apparent should the stretching parameter fall below unity. Another advantage of converting S(Q,) to I(Q,t) is for direct comparison with data collected using NSE. Such comparison can allow relaxation phenomena to be followed out to much longer time scales and thus modelled with much greater certainty; at least at certain momentum transfer vectors. In this task we will look to analyse data taken at 300K in the time domain. Analysis is a 2 stage process i.e. first we FFT the spectra, changing energy to time, and then we fit. Mark Telling NBIAW, th November
33 Stage 1: 1. Using the Convert To Energy GUI, reduce data file (lyophilised protein, 300K) such that a) run (lyophilised protein, 10K) is used to create the detector efficiency file and b) a suitable mapping file is produced that collates every 5 detectors. 2. Reduce the empty samplecan data and the base temp sample file as above 3. Using the Indirect Data Analysis GUI, subtract the empty samplecell data from both reduced sample files 4. Select Indirect Data Analysis > Fury 5. Browsed for the 300K '_subtract_' data file and Plot Input 7. For Resolution, browse for the 10K '_subtract_' data file 8. Adjust E low, E high and E width to be 0.5, 0.01 and 0.5 respectively 9. Check Save Result and click Run 10. The resulting workspace from Fury analysis is named 'inst'run No'_AnalyserReflection_iqt' I(Q,t) data for spectra 1,4,7, 10 from IRS38157_graphite002_iqt Mark Telling NBIAW, th November
34 Stage 2: There are several options to fit I(Q,t) curves a) the FuryFit GUI, b) creating a user defined fit function in Fit Wizard or c) using an external program such as DAVE FuryFit: The following details how to perform I(Q,t) analysis using the FuryFit GUI a GUI very similar in operation to the ConvFit analysis tool described earlier. 1. Select File or workspace IRS38157_graphite002_iqt. A plot of the data in spectrum N will appear (default for N is spectrum 0). Select spectrum 5 2. Select Fit Type: 1 Stretched Exponential 3. Click Constrain Intensities to insure that the sum of the initial amplitudes of all individual fit components is 1 4. Define the range of the fit by either entering values in the StartX / End X fields or dragging the blue lines on the plot 5. Define the background level by either entering values in the LinearBackground field or by dragging the green base line up to the required level 6. Selecting Plot Guess automatically adds a fit to the data based only the initial values associated with each fit function Property. 7. Once initial values are set, click Run to fit the spectrum displayed in the GUI only NB Selecting Plot Output launches a new window (when Run is pressed) showing the single fit, data and residuals in greater detail and only over the x and y fit range entered 8. Function parameters can be fixed by right clicking on the parameter name and selecting Fix. Fixed parameters are unfixed by following the same procedure Mark Telling NBIAW, th November
35 9. Once happy with the fit of a single spectrum, select which fit parameter to plot (e.g. tau, intensity or beta), select Save Result and click Run Sequential Fit to analyse all the spectra 10. Three workspaces are created: Inst Run No _graphite002_fury_1s_s0_to_9_result contains fit parameter values for all spectra in WORKSPACE form which can be plotted individually Inst Run No _graphite002_fury_1s_s0_to_9_parameters contains fit parameter values for all spectra in TABLE form Inst Run No _graphite002_fury_1s_s0_to_9_workspace contains data (ID 0), fits (ID 1) and residuals (ID 2) for each spectrum. Expand the workspace (click on the side arrow) to see individual spectrum workspaces. To create the multiplot below, select four workspaces > right click > Plot Spectrum > enter, 01 (for data and fit only) I(Q,t) data and sequential fits for spectra 1,4,7, 10 from IRS38157_graphite002_iqt Mark Telling NBIAW, th November
36 11. Examine the fits and fit parameters. Are they reasonable? Are any fit parametrs constant (within error) as a function of Q? Could any be fixed at a single value to reduce uncertainty in the other parameters? NB If beta is constant, within error, then it can be constrained at a mean value by either a) FIXING the parameter OR b) clicking Constrain Beta Over All Q. The latter creates the output files detailed above BUT with the additional description _1Smult 12. Analyse the I(Q,t) data using FuryFit and ensure you can broadly replicate the plots in the reference SoftMatter paper. NB: for consistency with the data in the SoftMatter paper fix beta at 0.64 Fit Wizard: Analysis can also be performed using the Fit Wizard GUI described earlier. As such only additional information required to analyse the data will be given below. 1. Using irs38157_graphite002_iqt, create a plot showing all I(Q,t) curves > right click, Plot Spectrum, enter spectrum ID numbers : Adjust axis values to: x = 0 to 0.1. y = 0.82 to Click on the graph that appears and select Main Menu > Analysis > Fit Wizard 3. Create and save a new User function with the form : (1A)*exp((x/t)^B))+A 4. Select Fit With Selected User Function and click Fit >> 5. In Custom Output >> click One Table For All Fits, click Parameter Table and then << Fit. A new table named Parameters should appear 6. Select Curve ***sp1, change fit range to x = 0 to 0.1 and click Fit. The fit parameters should be appended to the Parameter table 7. Repeat for Curves ****sp2 to ****sp Manipulate the fit parameters and plot the Parameter table values to replicate the plots in the reference SoftMatter paper Mark Telling NBIAW, th November
37 DAVE: For completeness, the following details how to perform I(Q,t) analysis using the very versatile IDL based package, DAVE, which can be downloaded here:www.ncnr.nist.gov/dave/ Taking a monkeyseemonkeydo approach to speed up analysis 1 In Mantid, click on the Algorithms tab 2 In the Execute input field, type SaveDaveGrp and click Execute 3 A new dialogue box will appear which will allow us to convert the Mantid 'inst'run No'_AnalyserReflection_iqt' file into a format DAVE can read. Select the relevant '_iqt' file, check ToMicroEv and enter a filename for the written dataset 4 Click Run 5 Launch DAVE by typing Dave 6 In the DAVE Gui > File > Preferences : Set Data And Working Directory 7 Load your I(Q,t) file by following: Data Input/Output > Read Dataset From > ASCII 8 Click on OK when the Specify Axes Labels For Data window appears 9 Data Input/Output > Write Dataset As > DAVE  a file with the extension.dave is created in your working directory 10 Analysis > 1D Peak Fitting (PAN). The PAN GUI will launch and the data in the.dave file should be visible in the plot window 11 In the PAN GUI, define the fit range by following File > Preferences and check / enter the following: Mark Telling NBIAW, th November
38 9. Click Accept and the Dismiss 10. To write a function: Select Function > Userfunction. Create a simple exponential decay function. Save your expression Q. Do you need to add a background to your equation? What does any background (or asymptotic I(Q,t) level) intensity tell you? 11. Click Accept > enter suitable start values > click Apply And Dismiss 12. Click Fit Current Group Q. Does a simple exponential provide a suitable description of the data? Note: linear axes can be deceptive! Examine the data with logx and logy axes by following Plot_Options > XLog and YLog. Focus the plot range by following File > Preferences. Suggestion: x = 0.01 to 0.1 and y = 0.87 to Click Fit All Groups and examine the goodness of fit at different Q values. 14. Click Clear All Curves and create a new stretched exponential function. Refit the data Mark Telling NBIAW, th November
39 Q. Is the stretched exponential description more appropriate? What does the stretched exponential form tell you about the internal mobility of the protein? 15. To examine the resulting fit parameters click > Plot Fit Parameter and select the parameter to view. For example, for the simple exponential model Parameter #0 will probably be the background level and Parameter #1 will probably be the decay time, Q. Compared to the water example would you say that the relaxation rate,, is Q dependent? Q. What does a Q independent relaxation rate suggest about the motion being probed? Q. To improve the fit / reduce errors could you fix one of the parameters in your model? Refit the data with or fixed. To perform a new fit, first click on Clear All Curves and then follow the procedure above Task Eight: Fitting the Elastic Incoherent Structure Factor (EISF) The EISF, A o (Q), provides information about the geometry of a particular localised motion and can be extracted from I(Q,t) if the decay curve reaches a plateau at 'long' decay times. From the I(Q,t) analysis above, the Q dependence of the background level gives reasonable access to the EISF; or, in the case of a partially mobile sample as is the case here, an effective structure factor. The theoretical Q dependence of elastic scattering intensity has been computed for many different types of localised motions. For a comprehensive list see Bée M., QuasiElastic Neutron Scattering, 1988, Adam Hilger, Chapter 6, p 200). As examples: Mark Telling NBIAW, th November
40 A o,ch3 (Q) = 1/3[1 + 2j o (Qr 3)] the theoretical EISF expected for a proton undergoing a localised 3fold jump rotation. r = Angstrom and j o is the zeroorder Bessel function. A o,jump (Q) = 1/2[1+j o (Qd)] the theoretical EISF response predicted for protons undergoing 2site jumps of distance, d A o,diff (Q) = (3j 1 (Qr))/Qr 2 the theoretical EISF response predicted for protons diffusing on a sphere of radius, r, Use the equations above to model the EISF response generated from your I(Q,t) analysis. NB: when fitting the model EISF functions remember that your experimental EISF is in fact an effective EISF since only a certain percentage of the sample is actually mobile the theoretical EISF functions assume the entire sample is mobile! The theoretical expressions therefore have to be modified to reflect this : 1. From your Mantid I(Q,t) FuryFit analysis, plot the Qdependence of your background parameter 2. Fit the data using FitWIzard and the three EISF expressions above; suitably modified to describe an effective EISF Q. Which of the three models provide the best description of the experimental data? Q. From your fits, what percentage of the protons in the material is mobile? Q. What would you need to do to truly ascertain the model which best describes the experimental data? Mark Telling NBIAW, th November
41 Suggestions: Further Analysis By modelling the temperature dependence of for the lyophilised material, determine a mean activation energy, E a, ave, for the protein dynamics observed. How does this compare to E a, ave determined from other dry proteins? Mark Telling NBIAW, th November
42 Example Three: Dental Cements Right...now we have all the tools to explore QENS data we will focus on data collected from dental cements; data which is so far is unanalysed and the 'result' unknown. The rational behind this work is given in the proposal (raw data directory) submitted to gain access to beam time on the IRIS instrument. Use the various data analysis methods explored in Examples One and Two to gain insight into the role of water in the solidification of dental cements. Don t forget to make use of advice from all the experimental team! Data Files File Name Temp Description K 28 day old HEMA cement K K Elastic Window Scan Data 54435, 54498, K Liquid POLY 54436, K 54437, K K Elastic Window Scan Data K Liquid HEMA K K K Vanadium Mark Telling NBIAW, th November
Experiment 5: Exploring Resolution, Signal, and Noise using an FTIR CH3400: Instrumental Analysis, Plymouth State University, Fall 2013
Experiment 5: Exploring Resolution, Signal, and Noise using an FTIR CH3400: Instrumental Analysis, Plymouth State University, Fall 2013 Adapted from JP Blitz and DG Klarup, "SignaltoNoise Ratio, Signal
More informationTutorial: Crystal structure refinement of oxalic acid dihydrate using GSAS
Tutorial: Crystal structure refinement of oxalic acid dihydrate using GSAS The aim of this tutorial is to use GSAS to locate hydrogen in oxalic acid dihydrate and refine the crystal structure. By no means
More informationUsing Excel for Graphical Analysis of Data
Using Excel for Graphical Analysis of Data Introduction In several upcoming labs, a primary goal will be to determine the mathematical relationship between two variable physical parameters. Graphs are
More informationQuick Instruction. For EXAFS Data Analysis
Quick Instruction For EXAFS Data Analysis Using WinXAS, FEFF 7 and ATOMS Programs Farideh Jalilehvand Department of Chemistry University of Calgary email: faridehj@ucalgary.ca Last updated: December 2002
More informationFundamentals of Rietveld Refinement III. Additional Examples
Fundamentals of Rietveld Refinement III. Additional Examples An Introduction to Rietveld Refinement using PANalytical X Pert HighScore Plus v3.0d Scott A Speakman, Ph.D. MIT Center for Materials Science
More informationModeling Evaporating Liquid Spray
Tutorial 17. Modeling Evaporating Liquid Spray Introduction In this tutorial, the airblast atomizer model in ANSYS FLUENT is used to predict the behavior of an evaporating methanol spray. Initially, the
More informationFundamentals of Rietveld Refinement III. Refinement of a Mixture
Fundamentals of Rietveld Refinement III. Refinement of a Mixture An Introduction to Rietveld Refinement using PANalytical X Pert HighScore Plus v3.0e Scott A Speakman, Ph.D. MIT Center for Materials Science
More informationFundamentals of Rietveld Refinement II. Refinement of a Single Phase
Fundamentals of Rietveld Refinement II. Refinement of a Single Phase An Introduction to Rietveld Refinement using PANalytical X Pert HighScore Plus v3.0a Scott A Speakman, Ph.D. MIT Center for Materials
More informationksa 400 Growth Rate Analysis Routines
kspace Associates, Inc., 2182 Bishop Circle East, Dexter, MI 48130 USA ksa 400 Growth Rate Analysis Routines Table of Contents ksa 400 Growth Rate Analysis Routines... 2 1. Introduction... 2 1.1. Scan
More informationChromatography Software Training Materials. Contents
Chromatography Software Training Materials This document contains information on how to build a method, start the instrument to acquire data, and then process the data using the Galaxie Program. You will
More informationTutorial 7: Automated Peak Picking in Skyline
Tutorial 7: Automated Peak Picking in Skyline Skyline now supports the ability to create custom advanced peak picking and scoring models for both selected reaction monitoring (SRM) and dataindependent
More informationDynamics and Vibration. Tutorial
Dynamics and Vibration Tutorial Startup To use Dynamics and Vibration Analysis (DVA), you must first start TK Solver. Once in TK, select Dynamics & Vibration from the Applications Menu. The DVA Menu will
More informationUsing a Single Rotating Reference Frame
Tutorial 9. Using a Single Rotating Reference Frame Introduction This tutorial considers the flow within a 2D, axisymmetric, corotating disk cavity system. Understanding the behavior of such flows is
More informationssrna Example SASSIE Workflow Data Interpolation
NOTE: This PDF file is for reference purposes only. This lab should be accessed directly from the web at https://sassieweb.chem.utk.edu/sassie2/docs/sample_work_flows/ssrna_example/ssrna_example.html.
More informationLecture 7 Notes: 07 / 11. Reflection and refraction
Lecture 7 Notes: 07 / 11 Reflection and refraction When an electromagnetic wave, such as light, encounters the surface of a medium, some of it is reflected off the surface, while some crosses the boundary
More informationGraphical Analysis of Data using Microsoft Excel [2016 Version]
Graphical Analysis of Data using Microsoft Excel [2016 Version] Introduction In several upcoming labs, a primary goal will be to determine the mathematical relationship between two variable physical parameters.
More informationANSYS Workbench Guide
ANSYS Workbench Guide Introduction This document serves as a stepbystep guide for conducting a Finite Element Analysis (FEA) using ANSYS Workbench. It will cover the use of the simulation package through
More informationHeat Exchanger Efficiency
6 Heat Exchanger Efficiency Flow Simulation can be used to study the fluid flow and heat transfer for a wide variety of engineering equipment. In this example we use Flow Simulation to determine the efficiency
More informationSkyline Targeted Method Editing
Skyline Targeted Method Editing This tutorial will cover many of the features available in the Skyline Targeted Proteomics Environment for creating new instrument methods for Selected Reaction Monitoring
More informationSteady Flow: LidDriven Cavity Flow
STARCCM+ User Guide Steady Flow: LidDriven Cavity Flow 2 Steady Flow: LidDriven Cavity Flow This tutorial demonstrates the performance of STARCCM+ in solving a traditional square liddriven cavity
More informationAppendix C. Vernier Tutorial
C1. Vernier Tutorial Introduction: In this lab course, you will collect, analyze and interpret data. The purpose of this tutorial is to teach you how to use the Vernier System to collect and transfer
More informationLogger Pro 3. Quick Reference
Logger Pro 3 Quick Reference Getting Started Logger Pro Requirements To use Logger Pro, you must have the following equipment: Windows 98, 2000, ME, NT, or XP on a Pentium processor or equivalent, 133
More informationTraceFinder Analysis Quick Reference Guide
TraceFinder Analysis Quick Reference Guide This quick reference guide describes the Analysis mode tasks assigned to the Technician role in the Thermo TraceFinder 3.0 analytical software. For detailed descriptions
More informationCurve Correction in Atomic Absorption
Curve Correction in Atomic Absorption Application Note Atomic Absorption Authors B. E. Limbek C. J. Rowe Introduction The Atomic Absorption technique ultimately produces an output measured in optical units
More informationChapter 3: Rate Laws Excel Tutorial on Fitting logarithmic data
Chapter 3: Rate Laws Excel Tutorial on Fitting logarithmic data The following table shows the raw data which you need to fit to an appropriate equation k (s 1 ) T (K) 0.00043 312.5 0.00103 318.47 0.0018
More informationLesson 1: Routine ITC Data Analysis and Fitting. Routine ITC Data Analysis. Lesson 1:Routine ITC Data Analysis and Fitting
Lesson 1:Routine ITC Data Analysis and Fitting Lesson 1: Routine ITC Data Analysis and Fitting In this lesson you will learn to perform routine analysis of ITC data. For reasonably good data, Origin makes
More informationData integration and scaling
Data integration and scaling Harry Powell MRC Laboratory of Molecular Biology 3rd February 2009 Abstract Processing diffraction images involves three basic steps, which are indexing the images, refinement
More informationHarris: Quantitative Chemical Analysis, Eight Edition CHAPTER 05: QUALITY ASSURANCE AND CALIBRATION METHODS
Harris: Quantitative Chemical Analysis, Eight Edition CHAPTER 05: QUALITY ASSURANCE AND CALIBRATION METHODS 50. International Measurement Evaluation Program Sample: Pb in river water (blind sample) :
More informationCapstone Appendix. A guide to your lab computer software
Capstone Appendix A guide to your lab computer software Important Notes Many of the Images will look slightly different from what you will see in lab. This is because each lab setup is different and so
More informationPCI1C. User Guide. Port City Instruments, LLC Market Street, Suite A271 Wilmington, NC
User Guide Port City Instruments, LLC 8209 Market Street, Suite A271 Wilmington, NC 28411 8664562488 www.portcityinstruments.com Copyright 2014 Port City Instruments, LLC. All Rights Reserved. This document
More informationGeneral Program Description
General Program Description This program is designed to interpret the results of a sampling inspection, for the purpose of judging compliance with chosen limits. It may also be used to identify outlying
More information2IN35 VLSI Programming Lab Work Assignment 1: Hardware design using Verilog
2IN35 VLSI Programming Lab Work Assignment 1: Hardware design using Verilog Hrishikesh Salunkhe, h.l.salunkhe@tue.nl, Alok Lele, a.lele@tue.nl April 28, 2015 1 Contents 1 Introduction 3 2 Hardware design
More informationLab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.
More informationData and Function Plotting with MATLAB (Linux10)
Data and Function Plotting with MATLAB (Linux10) This tutorial describes the use of MATLAB for general plotting of experimental data and equations and for special plots like histograms. (Astronomers 
More information7 Fractions. Number Sense and Numeration Measurement Geometry and Spatial Sense Patterning and Algebra Data Management and Probability
7 Fractions GRADE 7 FRACTIONS continue to develop proficiency by using fractions in mental strategies and in selecting and justifying use; develop proficiency in adding and subtracting simple fractions;
More informationError Analysis, Statistics and Graphing
Error Analysis, Statistics and Graphing This semester, most of labs we require us to calculate a numerical answer based on the data we obtain. A hard question to answer in most cases is how good is your
More informationModel FP6500 Spectrofluorometer Instruction Manual. FP6500 for Windows
Model FP6500 Spectrofluorometer Instruction Manual FP6500 for Windows P/N: 03029999 April 2000 Contents Safety Considerations...i Regulatory Statements... iii Preface... iv Installation Conditions...v
More informationcief Data Analysis Chapter Overview Chapter 12:
page 285 Chapter 12: cief Data Analysis Chapter Overview Analysis Screen Overview Opening Run Files How Run Data is Displayed Viewing Run Data Data Notifications and Warnings Checking Your Results Group
More informationIntroduction to Hermes
Introduction to Hermes Version 2.0 November 2017 Hermes v1.9 Table of Contents Introduction... 2 Visualising and Editing the MLL1 fusion protein... 2 Opening Files in Hermes... 3 Setting Style Preferences...
More informationFollowing a tour is the easiest way to learn Prism.
Page 1 of 25 A tour of Prism Following a tour is the easiest way to learn Prism. View a movie Watch and listen to a ten minute introductory movie from Prism's Welcome dialog. Or view it on the web. Read
More informationEffective Medium Theory, Rough Surfaces, and Moth s Eyes
Effective Medium Theory, Rough Surfaces, and Moth s Eyes R. Steven Turley, David Allred, Anthony Willey, Joseph Muhlestein, and Zephne Larsen Brigham Young University, Provo, Utah Abstract Optics in the
More informationT 1 Relaxation Measurement: The InversionRecovery Experiment (Using IconNMR) Revised
T 1 Relaxation Measurement: The InversionRecovery Experiment (Using IconNMR) Revised 4212015 Relaxation times are divided into two types: longitudinal, which concerns change in magnetization along the
More informationLaboratory 6: Light and the Laser
Laboratory 6: Light and the Laser WARNING NEVER LOOK DIRECTLY AT LASER LIGHT Index of Refraction: Snell's Law 1. Read the section on physical optics in some introductory physics text. 2. Set the semicircular
More informationImporting and processing a DGGE gel image
BioNumerics Tutorial: Importing and processing a DGGE gel image 1 Aim Comprehensive tools for the processing of electrophoresis fingerprints, both from slab gels and capillary sequencers are incorporated
More informationENVI Classic Tutorial: Multispectral Analysis of MASTER HDF Data 2
ENVI Classic Tutorial: Multispectral Analysis of MASTER HDF Data Multispectral Analysis of MASTER HDF Data 2 Files Used in This Tutorial 2 Background 2 Shortwave Infrared (SWIR) Analysis 3 Opening the
More informationspecular diffuse reflection.
Lesson 8 Light and Optics The Nature of Light Properties of Light: Reflection Refraction Interference Diffraction Polarization Dispersion and Prisms Total Internal Reflection Huygens s Principle The Nature
More informationScaling and Power Spectra of Natural Images
Scaling and Power Spectra of Natural Images R. P. Millane, S. Alzaidi and W. H. Hsiao Department of Electrical and Computer Engineering University of Canterbury Private Bag 4800, Christchurch, New Zealand
More informationXray Powder Diffraction
Xray Powder Diffraction Chemistry 754 Solid State Chemistry Lecture #8 April 15, 2004 Single Crystal Diffraction Diffracted Beam Incident Beam Powder Diffraction Diffracted Beam Incident Beam In powder
More informationUsing LoggerPro. Nothing is more terrible than to see ignorance in action. J. W. Goethe ( )
Using LoggerPro Nothing is more terrible than to see ignorance in action. J. W. Goethe (17491832) LoggerPro is a generalpurpose program for acquiring, graphing and analyzing data. It can accept input
More information= 3 + (5*4) + (1/2)*(4/2)^2.
Physics 100 Lab 1: Use of a Spreadsheet to Analyze Data by Kenneth Hahn and Michael Goggin In this lab you will learn how to enter data into a spreadsheet and to manipulate the data in meaningful ways.
More informationExperiment 6. Snell s Law. Use Snell s Law to determine the index of refraction of Lucite.
Experiment 6 Snell s Law 6.1 Objectives Use Snell s Law to determine the index of refraction of Lucite. Observe total internal reflection and calculate the critical angle. Explain the basis of how optical
More informationDiffraction. Singleslit diffraction. Diffraction by a circular aperture. Chapter 38. In the forward direction, the intensity is maximal.
Diffraction Chapter 38 Huygens construction may be used to find the wave observed on the downstream side of an aperture of any shape. Diffraction The interference pattern encodes the shape as a Fourier
More informationImport and preprocessing of raw spectrum data
BioNumerics Tutorial: Import and preprocessing of raw spectrum data 1 Aim Comprehensive tools for the import of spectrum data, both raw spectrum data as processed spectrum data are incorporated into BioNumerics.
More informationThe Mathcad Workspace 7
For information on system requirements and how to install Mathcad on your computer, refer to Chapter 1, Welcome to Mathcad. When you start Mathcad, you ll see a window like that shown in Figure 21. By
More informationModeling Evaporating Liquid Spray
Tutorial 16. Modeling Evaporating Liquid Spray Introduction In this tutorial, FLUENT s airblast atomizer model is used to predict the behavior of an evaporating methanol spray. Initially, the air flow
More informationMAN0314/2.00. Software update summary
Filename: Software Update Notification for DTS v4.20.doc Doct. Version No: 1 Date: Feb 2006 Author: M.T.Connah Product Software suite Product Zetasizer Nano series Firmware Manuals Zetasizer Nano series
More informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More informationLab Activity #2 Statistics and Graphing
Lab Activity #2 Statistics and Graphing Graphical Representation of Data and the Use of Google Sheets : Scientists answer posed questions by performing experiments which provide information about a given
More informationIntroduction to Raman spectroscopy measurement data processing using Igor Pro
Introduction to Raman spectroscopy measurement data processing using Igor Pro This introduction is intended to minimally guide beginners to processing Raman spectroscopy measurement data, which includes
More informationIntroduction to C omputational F luid Dynamics. D. Murrin
Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena
More informationTo Measure a Constant Velocity. Enter.
To Measure a Constant Velocity Apparatus calculator, black lead, calculator based ranger (cbr, shown), Physics application this text, the use of the program becomes second nature. At the Vernier Software
More informationUsing HyperChem for CH 463. Chun Park
Using HyperChem for CH 463 Chun Park Using HyperChem for CH 463 1 st Edition Written by Chun Park Oregon State University Table of Contents Introduction 1 Chapter 1 Getting Started 2 Starting HyperChem...
More informationEquipment Support Structures
Page 1 Equipment Support Structures Preface Using This Guide Where to Find More Information Conventions What's New? Getting Started Setting Up Your Session Creating a Simple Structural Frame Creating Nonuniform
More informationCHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force
CHAPTER 4 Numerical Models This chapter presents the development of numerical models for sandwich beams/plates subjected to fourpoint bending and the hydromat test system. Detailed descriptions of the
More informationISO INTERNATIONAL STANDARD. Particle size analysis Laser diffraction methods. Analyse granulométrique Méthodes par diffraction laser
INTERNATIONAL STANDARD ISO 13320 First edition 20091001 Corrected version 20091201 Particle size analysis Laser diffraction methods Analyse granulométrique Méthodes par diffraction laser Reference
More informationChromQuest 5.0 Quick Reference Guide
ChromQuest 5.0 Quick Reference Guide This guide contains an overview of the ChromQuest chromatography data system, with topics organized by workflow. For more information, refer to the ChromQuest User
More informationScientific Graphing in Excel 2007
Scientific Graphing in Excel 2007 When you start Excel, you will see the screen below. Various parts of the display are labelled in red, with arrows, to define the terms used in the remainder of this overview.
More informationMonte Carlo simulation of photon and electron transport
First Barcelona Techno Week Course on semiconductor detectors ICCUB, 1115th July 2016 Monte Carlo simulation of photon and electron transport Francesc Salvat Monte Carlo 1 Simulations performed with the
More informationAdvanced ANSYS FLUENT Acoustics
Workshop Modeling FlowInduced (Aeroacoustic) Noise 14.5 Release Advanced ANSYS FLUENT Acoustics 2011 ANSYS, Inc. November 7, 2012 1 Introduction This tutorial demonstrates how to model 2D turbulent flow
More informationChapter 4 Determining Cell Size
Chapter 4 Determining Cell Size Chapter 4 Determining Cell Size The third tutorial is designed to give you a demonstration in using the Cell Size Calculator to obtain the optimal cell size for your circuit
More informationSolution Recording and Playback: Vortex Shedding
STARCCM+ User Guide 6663 Solution Recording and Playback: Vortex Shedding This tutorial demonstrates how to use the solution recording and playback module for capturing the results of transient phenomena.
More informationLead Discovery 5.2. User Guide. Powered by TIBCO Spotfire
User Guide Powered by TIBCO Spotfire Last Modified: July 26, 2013 Table of Contents 1. Introduction... 5 2. Loading Data... 6 2.1. Opening an SDFile... 6 2.2. Importing a ChemDraw for Excel File... 6 2.3.
More informationDetermination of the aperture of the LHCb VELO RF foil
LHCbPUB21412 April 1, 214 Determination of the aperture of the LHCb VELO RF foil M. FerroLuzzi 1, T. Latham 2, C. Wallace 2. 1 CERN, Geneva, Switzerland 2 University of Warwick, United Kingdom LHCbPUB21412
More informationOperating instructions  Cary 100 Bio UVvisible Spectrophotometer
Operating instructions  Cary 100 Bio UVvisible Spectrophotometer This document describes how to power up the spectrophotometer, set its measurement parameters, insert one or more of samples for analysis
More informationIntroduction to Solid Modeling Using SolidWorks 2008 COSMOSMotion Tutorial Page 1
Introduction to Solid Modeling Using SolidWorks 2008 COSMOSMotion Tutorial Page 1 In this tutorial, we will learn the basics of performing motion analysis using COSMOSMotion. Although the tutorial can
More informationContents Contents Creating a Simulation Example: A Dipole Antenna AMDS User s Guide
Contents Contents 1 Creating a Simulation 7 Introduction 8 Data Files for Examples 8 Software Organization 9 Constructing the Geometry 10 Creating the Mesh 11 Defining Run Parameters 13 Requesting Results
More informationFigure 1. Figure 2. The BOOTSTRAP
The BOOTSTRAP Normal Errors The definition of error of a fitted variable from the variancecovariance method relies on one assumption that the source of the error is such that the noise measured has a
More informationISOCS Characterization of Sodium Iodide Detectors for GammaRay Spectrometry
ISOCS Characterization of Sodium Iodide Detectors for GammaRay Spectrometry Sasha A. Philips, Frazier Bronson, Ram Venkataraman, Brian M. Young AbstractActivity measurements require knowledge of the
More informationPhysics 101, Lab 1: LINEAR KINEMATICS PREDICTION SHEET
Physics 101, Lab 1: LINEAR KINEMATICS PREDICTION SHEET After reading through the Introduction, Purpose and Principles sections of the lab manual (and skimming through the procedures), answer the following
More informationMuon Data Analysis Workshop: An Introduction to the Practicalities
Muon Data Analysis Workshop: An Introduction to the Practicalities Francis Pratt (ISIS) ISIS Muon Training Course March 2012 Outline A. Acquiring the data B. Location and format of the data C. Data analysis
More informationc Fluent Inc. May 16,
Tutorial 1. Office Ventilation Introduction: This tutorial demonstrates how to model an office shared by two people working at computers, using Airpak. In this tutorial, you will learn how to: Open a new
More informationReflection and Refraction of Light
PC1222 Fundamentals of Physics II Reflection and Refraction of Light 1 Objectives Investigate for reflection of rays from a plane surface, the dependence of the angle of reflection on the angle of incidence.
More informationLesson 1: Creating T Spline Forms. In Samples section of your Data Panel, browse to: Fusion 101 Training > 03 Sculpt > 03_Sculpting_Introduction.
3.1: Sculpting Sculpting in Fusion 360 allows for the intuitive freeform creation of organic solid bodies and surfaces by leveraging the T Splines technology. In the Sculpt Workspace, you can rapidly
More informationExcel 2013 Intermediate
Instructor s Excel 2013 Tutorial 2  Charts Excel 2013 Intermediate 103124 Unit 2  Charts Quick Links Chart Concepts Page EX197 EX199 EX200 Selecting Source Data Pages EX198 EX234 EX237 Creating a Chart
More informationTransactions on Modelling and Simulation vol 10, 1995 WIT Press, ISSN X
Hydrodynamic coefficients and motions due to a floating cylinder in waves D.D. Bhatta, M. Rahman Department of Applied Mathematics, Technical University of Nova Scotia, Halifax, Nova Scotia, Canada B3J
More informationUMASIS, AN ANALYSIS AND VISUALIZATION TOOL FOR DEVELOPING AND OPTIMIZING ULTRASONIC INSPECTION TECHNIQUES
UMASIS, AN ANALYSIS AND VISUALIZATION TOOL FOR DEVELOPING AND OPTIMIZING ULTRASONIC INSPECTION TECHNIQUES A.W.F. Volker, J. G.P. Bloom TNO Science & Industry, Stieltjesweg 1, 2628CK Delft, The Netherlands
More informationVersion 1.0 November2016 Hermes V1.8.2
Hermes in a Nutshell Version 1.0 November2016 Hermes V1.8.2 Table of Contents Hermes in a Nutshell... 1 Introduction... 2 Example 1. Visualizing and Editing the MLL1 fusion protein... 3 Setting Your Display...
More informationSpreadsheet View and Basic Statistics Concepts
Spreadsheet View and Basic Statistics Concepts GeoGebra 3.2 Workshop Handout 9 Judith and Markus Hohenwarter www.geogebra.org Table of Contents 1. Introduction to GeoGebra s Spreadsheet View 2 2. Record
More informationU90 Ladder Software Manual. Version 3.50, 6/03
U90 Ladder Software Manual Version 3.50, 6/03 Table Of Contents Welcome to U90 Ladder... 1 Program Editors... 1 Project Navigation Tree...1 Browse Sequences...1 Printing Documentation...2 Interface Language...
More informationControl of Light. Emmett Ientilucci Digital Imaging and Remote Sensing Laboratory Chester F. Carlson Center for Imaging Science 8 May 2007
Control of Light Emmett Ientilucci Digital Imaging and Remote Sensing Laboratory Chester F. Carlson Center for Imaging Science 8 May 007 Spectroradiometry Spectral Considerations Chromatic dispersion
More informationTorsionallateral buckling large displacement analysis with a simple beam using Abaqus 6.10
Torsionallateral buckling large displacement analysis with a simple beam using Abaqus 6.10 This document contains an Abaqus tutorial for performing a buckling analysis using the finite element program
More informationANSYS AIM Tutorial Structural Analysis of a Plate with Hole
ANSYS AIM Tutorial Structural Analysis of a Plate with Hole Author(s): Sebastian Vecchi, ANSYS Created using ANSYS AIM 18.1 Problem Specification PreAnalysis & Start Up Analytical vs. Numerical Approaches
More informationINVESTIGATIONS OF CROSSCORRELATION AND EUCLIDEAN DISTANCE TARGET MATCHING TECHNIQUES IN THE MPEF ENVIRONMENT. Greg And Ken Holmlund # ABSTRACT
INVESTIGATIONS OF CROSSCORRELATION AND EUCLIDEAN DISTANCE TARGET MATCHING TECHNIQUES IN THE MPEF ENVIRONME Greg Dew @ And Ken Holmlund # @ Logica # EUMETSAT ABSTRACT CrossCorrelation and Euclidean Distance
More informationGeostatistics 2D GMS 7.0 TUTORIALS. 1 Introduction. 1.1 Contents
GMS 7.0 TUTORIALS 1 Introduction Twodimensional geostatistics (interpolation) can be performed in GMS using the 2D Scatter Point module. The module is used to interpolate from sets of 2D scatter points
More information[1] CURVE FITTING WITH EXCEL
1 Lecture 04 February 9, 2010 Tuesday Today is our third Excel lecture. Our two central themes are: (1) curvefitting, and (2) linear algebra (matrices). We will have a 4 th lecture on Excel to further
More informationModeField Diameter and Spot Size Measurements of Lensed and Tapered Specialty Fibers
ModeField Diameter and Spot Size Measurements of Lensed and Tapered Specialty Fibers By Jeffrey L. Guttman, Ph.D., Director of Engineering, OphirSpiricon Abstract: The ModeField Diameter (MFD) and spot
More informationInteractive comment on Quantification and mitigation of the impact of scene inhomogeneity on Sentinel4 UVN UVVIS retrievals by S. Noël et al.
Atmos. Meas. Tech. Discuss., www.atmosmeastechdiscuss.net/5/c741/2012/ Author(s) 2012. This work is distributed under the Creative Commons Attribute 3.0 License. Atmospheric Measurement Techniques Discussions
More informationInteractive comment on Quantification and mitigation of the impact of scene inhomogeneity on Sentinel4 UVN UVVIS retrievals by S. Noël et al.
Atmos. Meas. Tech. Discuss., 5, C741 C750, 2012 www.atmosmeastechdiscuss.net/5/c741/2012/ Author(s) 2012. This work is distributed under the Creative Commons Attribute 3.0 License. Atmospheric Measurement
More informationProbabilistic Analysis Tutorial
Probabilistic Analysis Tutorial 21 Probabilistic Analysis Tutorial This tutorial will familiarize the user with the Probabilistic Analysis features of Swedge. In a Probabilistic Analysis, you can define
More information