US20060158183A1 - System and method for characterizing a sample by low-frequency spectra - Google Patents
System and method for characterizing a sample by low-frequency spectra Download PDFInfo
- Publication number
- US20060158183A1 US20060158183A1 US11/349,857 US34985706A US2006158183A1 US 20060158183 A1 US20060158183 A1 US 20060158183A1 US 34985706 A US34985706 A US 34985706A US 2006158183 A1 US2006158183 A1 US 2006158183A1
- Authority
- US
- United States
- Prior art keywords
- sample
- noise
- signal
- frequency
- spectral
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N37/00—Details not covered by any other group of this subclass
- G01N37/005—Measurement methods not based on established scientific theories
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/032—Measuring direction or magnitude of magnetic fields or magnetic flux using magneto-optic devices, e.g. Faraday or Cotton-Mouton effect
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
- Measuring Magnetic Variables (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
A method and apparatus for interrogating a sample that exhibits molecular rotation are disclosed. In practicing the method, the sample is placed in a container having both magnetic and electromagnetic shielding, and Gaussian noise is injected into the sample. An electromagnetic time-domain signal composed of sample source radiation superimposed on the injected Guassian noise is detected, and this signal is used to generate a spectral plot that displays, at a selected power setting of the Gaussian noise source, low-frequency spectral components characteristic of the sample in a selected frequency range between DC and 50 kHz. In one embodiment, the spectral plot that is generated is a histogram of stochastic resonance events over the selected frequency range. From this spectrum, one or more low-frequency signal components that are characteristic of the sample being interrogated are identified.
Description
- This application is a continuation-in-part, and claims the benefit of, U.S. patent application Ser. No. 10/112,927, filed Mar. 29, 2002, entitled APPARATUS AND METHOD FOR MEASURING MOLECULAR ELECTROMAGNETIC SIGNALS WITH A SQUID DEVICE AND STOCHASTIC RESONANCE TO MEASURE LOW-THRESHOLD SIGNALS (Attorney Docket No. 38547.8007.US00), International Patent Application No. PCT/US03/09544, filed Mar. 28, 2003, entitled SYSTEM AND METHOD FOR CHARACTERIZING A SAMPLE BY LOW-FREQUENCY SPECTRA (Attorney Docket No. 38547.8005.WO00), and International Patent Application No. PCT/US03/11834, filed Apr. 18, 2003, entitled SYSTEM AND METHOD FOR SAMPLE DETECTION BASED ON LOW-FREQUENCY SPECTRAL COMPONENTS (Attorney Docket No. 38547.8006.WO00), all incorporated herein by reference.
- There are a variety of spectroscopic tools for characterizing atomic or molecular compound. These include, but are not limited to, x-ray, UV, visible-light, infrared and microwave spectroscopy, and nuclear and electron spin resonance (NMR and ESR) spectroscopy. In general, spectroscopic tools are useful for at least four different type of chemical-analytical problems: first, to characterize an atomic and molecular compound according to its spectrographic features, e.g., spectral components; second, to determine the atomic composition of a compound, according to the spectral characteristics of atoms making up the compound; third, to determine 2-D or 3-D conformation of a molecular compound according to the spectral characteristic of atom-atom interactions in the compound; and fourth, to detect and identify components, e.g., contaminants, in a sample according to the distinguishing spectral characteristics of the compound being detected.
- Most existing spectroscopic tools provide some unique advantage(s) in terms of sensitivity, the information gained, ease of measurement and cost. Because each tool provides information not otherwise available, it is generally advantageous to be able to bring to bear on any chemical-analytical, as many pertinent spectroscopic tools as possible.
- The invention includes, in one aspect, a method of characterizing spectral emission features of a sample material, e.g., low-frequency emissions related to molecular motion within the sample. The method uses a time-domain signal of the sample over a sample-duration time T, and a sampling rate F for sampling the time domain signal, where F*T is the total sample count S, F is approximately twice the frequency domain resolution f of a Real Fast Fourier Transform of the time-domain signal sampled at sampling rate F, and S>f*n, where n is at least 10. The program selects S/n samples from the stored time domain signal and performs a Real Fast Fourier Transform (RFFT) on the samples. The RFFT is then normalized (e.g., by setting the highest value to 1), and an average power for the signal is calculated from the normalized signal. Next, the program places an event count in each of f selected-frequency event bins, where the measured power at the corresponding selected frequency≧average power*ε obtains, where 0<ε<1, and ε is chosen such that the total number of counts placed in an event bin is between about 20-50% of the maximum possible bin counts in that bin. These steps are repeated n times, generating a histogram that shows, for each event bin f over a selected frequency range, the number of event counts in each bin.
- The method may further include the step of placing the normalized power value from the RFFT in f corresponding-frequency power bins, and, after n cycles of program operation, (a) dividing the accumulated values placed in each of the f power bins by n, to yield an average power in each bin, and (b) displaying on the histogram, the average power in each bin. The method may further include identifying those bins in the histogram that have an event count above a given threshold and an average power.
- Also disclosed is a low-frequency spectral signature associated with a material of interest comprising a list of frequency components in the DC-50 kHz frequency range that are generated by the above method. The frequencies in the list may be identified from a histogram of the number of sample-dependent stochastic events occurring at each of a plurality of spectral increments within a selected frequency range between DC and 50 kHz.
- In another aspect, the invention includes an apparatus for interrogating a sample that exhibits low-frequency molecular motion. The apparatus includes a magnetically and electromagnetically shielded container adapted for receiving the sample, an adjustable-power source of Gaussian noise for injection into the sample, with the sample in said container, and a detector for detecting an electromagnetic time-domain signal composed of sample source radiation superimposed on the injected Gaussian noise. An electronic computer in the apparatus receives the time-domain signal from the detector, and processes the signal to generate a spectral plot that displays, at a selected power setting of the Gaussian noise source, low-frequency spectral components characteristic of the sample in a selected frequency range between DC and 50 kHz.
- In one general embodiment, the electronic computer includes a signal analyzer that functions to (i) calculate a series of Fourier spectra of the time-domain signal over each of a plurality of defined time periods, in a selected frequency range between DC and 50 kHz, and (ii) average the Fourier spectra.
- Typically, at least five Fourier spectra are calculated, each taken over a 1-5 second time-domain interval.
- In another general embodiment, the electronic computer includes machine-readable code operable to carry out the method described above for generating a histogram of spectral events.
- The source of Gaussian noise in the apparatus may be an adjustable-power Gaussian noise generator and a Helmholz coil which is contained within the magnetic cage and the Faraday cage, and which receives a selected noise output signal from the noise generator in the
range 100 mV to 1 V. The injector is designed to inject Gaussian noise into the sample at a frequency, for example, between DC and 8 kHz. - The detector in the apparatus may be a first-derivative superconducting gradiometer which outputs a current signal, and a SQUID operatively connected to the gradiometer to convert the current signal to an amplified voltage signal.
- The container in the apparatus may include an attenuation tube having a sample-holding region, a magnetic shielding cage surrounding the region, and a Faraday cage contained within the magnetic shielding cage and also surrounding the region. In this embodiment, the source of Gaussian noise may include a Gaussian noise generator and a Helmholz coil which is contained within the magnetic cage and the Faraday cage, and which receives a noise output signal from the noise generator, and which further includes, for use in removing stationary noise components in the time-dependent signal, a signal inverter operatively connected to the said noise source and to said SQUID, for receiving Gaussian noise from the noise source and outputting into said SQUID, Gaussian noise in inverted form with respect to the Gaussian noise injected into the sample.
- In still another aspect, the invention includes a method for interrogating a sample that exhibits low-frequency molecular motion. In practicing the method, the sample is placed in a container having both magnetic and electromagnetic shielding, and Gaussian noise is injected into the sample at a selected noise amplitude. An electromagnetic time-domain signal composed of sample source radiation superimposed on the injected Gaussian noise, is recorded, and from this, a spectral plot that contains, at a selected power setting of the Gaussian noise source, low-frequency, sample-dependent spectral components characteristic of the sample in a selected frequency range between DC and 50 kHz is generated. The steps are repeated at different selected noise amplitudes until a plot showing a maximum or near maximum number of spectral components characteristic of the sample are generated.
- In one embodiment, the spectral plot is generated by (i) calculating a series of Fourier spectra of the time-domain signal over each of a plurality of defined time periods, in a selected frequency range between DC and 50 kHz, and (ii) averaging the Fourier spectra.
- In another general embodiment, the spectral plot is generated by the histogram method above.
- These and other objects and features of the invention will become more fully apparent when the following detailed description of the invention is read in conjunction with the accompanying drawings.
-
FIG. 1 is an isometric view of one embodiment of a molecular electromagnetic signaling detection apparatus formed in accordance with one embodiment of the present invention; -
FIG. 2 is an enlarged, detail view of the faraday cage and its contents shown inFIG. 1 ; and -
FIG. 3 is an enlarged, cross sectional view of one of the attenuation tubes shown inFIGS. 1 and 2 . -
FIG. 4 is a cross-section view of the faraday cage and its contents shown inFIG. 2 . -
FIG. 5 is a cross-section view of an alternative embodiment of the invention shown inFIGS. 1 through 4 . -
FIG. 6 is an enlarged, detail view of the frames supporting the coils of the Helmholtz transformer described herein. -
FIG. 7 is a diagram of an alternative electromagnetic emission detection system. -
FIG. 8 diagram of the processing unit included in the detection system of the above Figures. -
FIG. 9 is a diagram of an alternative processing unit to that ofFIG. 8 . -
FIG. 10 is a flow diagram of the signal detection and processing performed by the present system. -
FIG. 11A is a spectral plot of the emissions of a first sample. -
FIG. 11B is spectral plot of the emissions of a second sample. -
FIGS. 12A and 12B are spectral plots, in the spectral region between 500-530 Hz, for a sample of saturated NaCl, generated by Fourier transforming a non-correlated time-domain sample signal (12A), and Fourier transforming a cross-correlated sample spectrum (12B). -
FIGS. 13A and 13B are spectral plots, in the spectral region between 500-530 Hz, for a sample of alkyl ether sulfate, generated by Fourier transforming a non-correlated time-domain sample signal (13A), and Fourier transforming a cross-correlated sample spectrum (13B). -
FIGS. 14A-14F are spectral plots, in the spectral region between 500-530 Hz, for samples of deionized water (14A), a saturated NaCl solution (14B), a solution of 1% NaCl in deionized water (14C); a saturated NaBr sample (14D), alkyl ether sulfate in deionized water (14E), and no sample (14F). -
FIGS. 15A-15F are spectral plots, in the spectral region between 500 and 535 Hz, of a sample of an amino acid at a 1:100 wt/volume solution (15A) and at increasing w/v dilutions of 1:10,000 (15B), 1:1 million (15C), 1:100 million (15D), 1:10 billion (15E and 15F), where the spectra inFIGS. 15A-15E were generated with 50 second recordings and 40 minute correlations, and the spectrum ofFIG. 15F was generated with a 4:25 minute recording with a 12 hour correlation. -
FIG. 16 is a schematic diagram illustrating an alternative embodiment of a molecular electromagnetic signaling detection apparatus. -
FIG. 17A is a cross-sectional view of the alternative embodiment ofFIG. 16 . -
FIG. 17B is an enlargement of a portion ofFIG. 17A . -
FIG. 18 is a cross-sectional isometric view ofFIG. 17B . -
FIG. 19 is a diagram of an alternative processing unit to that ofFIG. 9 . -
FIG. 20 shows a high-level flow diagram of data flow for the histogram spectral plot method of the invention; -
FIG. 21 is a flow diagram of the algorithm for generating a spectral plot histogram, in accordance with the invention, and -
FIGS. 22A-22D are histogram spectra of a sample taken at four different noise power levels. -
FIGS. 23A-23C are computer screen shots displaying a user interface for generating and displaying a spectral plot histogram, withFIG. 23C being in color. - In the drawings, identical reference numbers identify identical or substantially similar elements or acts. To easily identify the disc
- Discussion of any particular element or art, the most significant digit or digits in a reference number refer to the Figure number in which that element is first introduced.
- I. Definitions
- The terms below have the following definitions unless indicated otherwise.
- “Sample that exhibits molecular rotation” refers to a sample material, which may be in gaseous, liquid or solid form (other than a solid metal) in which one or more of the molecular compounds or atomic ions making up or present in the sample exhibit rotation.
- “Magnetic shielding” refers to shielding that inhibits or prevents passage of magnetic flux as a result of the magnetic permeability of the shielding material.
- “Electromagnetic shielding” refers to, e.g., standard Faraday electromagnetic shielding.
- “Time-domain signal” or “time-series signal” refers to a signal with transient signal properties that change over time.
- “Sample-source radiation” refers to magnetic flux emissions resulting from molecular motion of a sample, such as the rotation of a molecular dipole in a magnetic field.
- “Gaussian noise” means random noise having a Gaussian power distribution.
- “Stationary white Gaussian noise” means random Gaussian noise that has no predictable future components.
- “Frequency-domain spectrum” refers to a Fourier frequency plot of a time-domain signal.
- “Spectral components” refer to singular or repeating qualities within a time-domain signal that can be measured in the frequency, amplitude, and/or phase domains. Spectral components will typically refer to signals present in the frequency domain.
- “Similar sample,” with reference to a first sample, refers to the same sample or a sample having substantially the same sample components as the first sample.
- “Faraday cage” refers to an electromagnetic shielding configuration that provides an electrical path to ground for unwanted electromagnetic radiation, thereby quieting an electromagnetic environment.
- II. Apparatus
- Described in detail below is a system and method for detecting, processing, and presenting low frequency electromagnetic emissions or signals of a sample of interest. In one embodiment, a known white or Gaussian noise signal is introduced to the sample. The Gaussian noise is configured to permit the electromagnetic emissions from the sample to be sufficiently detected by a signal detection system. Sets of detected signals are processed together to ensure repeatability and statistical relevance. The resulting emission pattern or spectrum can be displayed, stored, and/or identified as a particular substance.
- The following description provides specific details for a thorough understanding of, and enabling description for, embodiments of the invention. However, one skilled in the art will understand that the invention may be practiced without these details. In other instances, well-known structures and functions have not been shown or described in detail to avoid unnecessarily obscuring the description of embodiments of the invention.
- As explained in detail below, embodiments of the present invention are directed to providing an apparatus and method for the repeatable detection and recording of low-threshold molecular electromagnetic signals. A magnetically shielded faraday cage shields the sample material and detection apparatus from extraneous electromagnetic signals. Within the magnetically shielded faraday cage, a coil injects white or Gaussian noise, a nonferrous tray holds the sample, and a gradiometer detects low-threshold molecular electromagnetic signals. The apparatus further includes a superconducting quantum interference device (“SQUID”) and a preamplifier.
- The apparatus is used by placing a sample within the magnetically shielded faraday cage in close proximity to the noise coil and gradiometer. White noise is injected through the noise coil and modulated until the molecular electromagnetic signal is enhanced through stochastic resonance. The enhanced molecular electromagnetic signal, shielded from external interference by the faraday cage and the field generated by the noise coil, is then detected and measured by the gradiometer and SQUID. The signal is then amplified and transmitted to any appropriate recording or measuring equipment.
- Referring to
FIG. 1 , there is shown a shieldingstructure 10 which includes, in an outer to inner direction, aconductive wire cage 16 which is a magnetic shield and innerconductive wire cages - Referring to
FIG. 2 , thefaraday cage 10 is open at the top, and includesside openings faraday cage 10 is further comprised of threecopper mesh cages copper mesh cages -
Side openings attenuation tubes faraday cage 10 while isolating the interior of the cage from external sources of interference. Referring toFIG. 3 ,attenuation tube 24 is comprised of threecopper mesh tubes copper mesh cages copper mesh tubes Attenuation tube 24 is further capped withcap 32, with thecap having hole 34.Attenuation tube 22 is similarly comprised ofcopper mesh tubes cap 32. - Referring again to
FIG. 2 , a low-densitynonferrous sample tray 50 is mounted in the interior of thefaraday cage 10. Thesample tray 50 is mounted so that it may be removed from thefaraday cage 10 through theattenuation tube 22 andside opening 12. Threerods 52, each of which is greater in length than the distance from the center vertical axis of thefaraday cage 10 to the outermost edge of theattenuation tube 22, are attached to thesample tray 50. The threerods 52 are adapted to conform to the interior curve of theattenuation tube 22, so that thesample tray 50 may be positioned in the center of thefaraday cage 10 by resting the rods in the attenuation tube. In the illustrated embodiment, thesample tray 50 androds 52 are made of glass fiber epoxy. It will be readily apparent to those skilled in the art that thesample tray 50 androds 52 may be made of other nonferrous materials, and the tray may be mounted in thefaraday cage 10 by other means, such as by a single rod. - Referring again to
FIG. 2 , mounted within thefaraday cage 10 and above thesample tray 50 is acryogenic dewar 100. In the disclosed embodiment, thedewar 100 is adapted to fit within the opening at the top offaraday cage 10 and is a Model BMD-6 Liquid Helium Dewar manufactured by Tristan Technologies, Inc. Thedewar 100 is constructed of a glass-fiber epoxy composite. Agradiometer 110 with a very narrow field of view is mounted within thedewar 100 in position so that its field of view encompasses thesample tray 50. In the illustrated embodiment, thegradiometer 110 is a first order axial detection coil, nominally 1 centimeter in diameter, with a 2% balance, and is formed from a superconductor. The gradiometer can be any form of gradiometer excluding a planar gradiometer. Thegradiometer 110 is connected to the input coil of one low temperature direct current superconducting quantum interference device (“SQUID”) 120. In the disclosed embodiment, the SQUID is a Model LSQ/20 LTS dc SQUID manufactured by Tristan Technologies, Inc. It will be recognized by those skilled in the art that high temperature or alternating current SQUIDs can be used without departing from the spirit and scope of the invention. In an alternative embodiment, theSQUID 120 includes anoise suppression coil 124. - The disclosed combination of
gradiometer 110 andSQUID 120 have a sensitivity of 5 microTesla/√Hz when measuring magnetic fields. - The output of
SQUID 120 is connected to a ModelSP Cryogenic Cable 130 manufactured by Tristan Technologies, Inc. TheCryogenic Cable 130 is capable of withstanding the temperatures within and without thedewar 100 and transfers the signal from theSQUID 120 to Flux-LockedLoop 140, which is mounted externally to thefaraday cage 10 anddewar 100. The Flux-LockedLoop 140 in the disclosed embodiment is an iFL-301-L Flux Locked Loop manufactured by Tristan Technologies, Inc. - Referring to
FIG. 1 , the Flux LockedLoop 140 further amplifies and outputs the signal received from theSQUID 120 via high-level output circuit 142 to an iMC-303 iMAG® SQUID controller 150. The Flux-LockedLoop 140 is also connected via a model CC-60 six-meter fiber-opticcomposite connecting cable 144 to theSQUID controller 150. The fiber-optic connecting cable 144 andSQUID controller 150 are manufactured by Tristan Technologies, Inc. Thecontroller 150 is mounted externally to themagnetic shielding cage 40. The fiber-optic connecting cable 144 carriers control signals from theSQUID controller 150 to the Flux LockedLoop 140, further reducing the possibility of electromagnetic interference with the signal to be measured. It will be apparent to those skilled in the art that other Flux-Locked Loops, connecting cables, and Squid controllers can be used without departing from the spirit and scope of the invention. - The
SQUID controller 150 further comprises high resolution analog todigital converters 152, a standard GP-IB bus 154 to output digitalized signals, andBNC connectors 156 to output analog signals. In the illustrated embodiment, the BNC connectors are connected to adual trace oscilloscope 160 throughpatch cord 162. - Referring to
FIG. 2 , a two-element Helmholtz transformer 60 is installed to either side of thesample tray 50 when the sample tray is fully inserted within thefaraday cage 10. In the illustrated embodiment, thecoil windings Helmholtz transformer 60 are designed to operate in the direct current to 50 kilohertz range, with a center frequency of 25 kilohertz and self-resonant frequency of 8.8 megahertz. In the illustrated embodiment, thecoil windings gradiometer 110 andsample tray 50 are positioned within the field produced by the Helmholtz coil. Each ofcoil windings nonferrous frames frames legs 70.Frames legs 70 to permit vertical movement of the frames in relation to the lower portion ofdewar 100. Movement of the frames permits adjustment of thecoil windings Helmholtz transformer 60 to vary the amplitude of white noise received atgradiometer 110. Thelegs 70 rest on or are epoxied onto the bottom of thefaraday cage 10. In the illustrated embodiment, theframes legs 70 are made of glass fiber epoxy. Other arrangements of transformers or coils may be used around thesample tray 50 without departing from the spirit and scope of the invention. - Referring to
FIG. 4 , there is shown a cross-sectional view of the faraday cage and its contents, showingwindings 62 ofHelmholtz transformer 60 in relation todewar 100 andfaraday cage 10. Note also inFIG. 4 the positioning ofsample tray 50 andsample 200. - Referring to
FIG. 5 , there is shown an alternative embodiment in which theHelmholtz coil windings additional noise coil 300 is positioned belowsample tray 50. The windings of theadditional noise coil 300 are substantially perpendicular to thevertical windings Helmholtz transformer 60, and the windings of theadditional noise coil 300 are thus substantially in parallel orientation to the bottom offaraday cage 10. - In this alternative embodiment, noise would be fed to
noise coil 300 from an identical twisted pair wire (not shown) as that supplying the Helmholtz coil. The noise source would originate with the same noise generator used to supply noise to the Helmholtz coil. Noise would be sampled either at the noise generator via an additional noise output connection, or via a balanced splitter from an output connection to the noise generator. Attenuation of the noise signal atadditional noise coil 300 would be through an adjustable RF signal attenuation circuit, of which many are available commercially, or via a suitable series of fixed value RF attenuation filters. - Referring to
FIG. 6 , a detail of the frames supporting the coils ofHelmholtz transformer 60 may be seen; the reference point ofFIG. 6 is 90 degrees from the view ofFIG. 4 , and omits thefaraday cage 10.Frames Frames 66′ and 68′ illustrate the rotation of said frames about the axis of the hinged connection joining said frames, so as to dispose the coil windings of the Helmholtz transformer in an non-parallel relationship with one another. - Referring again to
FIG. 1 , an amplitude adjustablewhite noise generator 80 is external tomagnetic shielding cage 40, and is electrically connected to theHelmholtz transformer 60 throughfilter 90 byelectrical cable 82. Referring toFIG. 3 ,cable 82 is run throughside opening 12,attenuation tube 24, and throughcap 32 viahole 34.Cable 82 is a co-axial cable further comprising a twisted pair ofcopper conductors 84 surrounded by interior and exterior magnetic shielding 86 and 88, respectively. In other embodiments, the conductors can be any nonmagnetic electrically conductive material, such as silver or gold. The interior and exterior magnetic shielding 86 and 88 terminates atcap 32, leaving thetwisted pair 84 to span the remaining distance from the end cap to theHelmholtz transformer 60 shown inFIG. 1 . The interior magnetic shielding 86 is electrically connected toFaraday cage 16 throughcap 32, while the exterior magnetic shielding is electrically connected to the magnetically shieldedcage 40 shown inFIG. 1 . - Referring to
FIG. 1 , thewhite noise generator 80 can generate nearly uniform noise across a frequency spectrum from zero to 100 kilohertz. In the illustrated embodiment, thefilter 90 filters out noise above 50 kilohertz, but other frequency ranges may be used without departing from the spirit and scope of the invention. -
White noise generator 80 is also electrically connected to the other input ofdual trace oscilloscope 160 throughpatch cord 164. - Referring to
FIGS. 1, 2 and 3, a sample of thesubstance 200 to be measured is placed on thesample tray 50 and the sample tray is placed within thefaraday cage 10. In the first embodiment, thewhite noise generator 80 is used to inject white noise through theHelmholtz transformer 60. The noise signal creates an induced voltage in thegradiometer 110. The induced voltage in thegradiometer 110 is then detected and amplified by theSQUID 120, the output from the SQUID is further amplified by the flux lockedloop 140 and sent to theSQUID controller 150, and then sent to thedual trace oscilloscope 160. Thedual trace oscilloscope 160 is also used to display the signal generated bywhite noise generator 80. - The white noise signal is adjusted by altering the output of the
white noise generator 80 and by rotating theHelmholtz transformer 60 around thesample 200, shown inFIG. 2 . Rotation of theHelmholtz transformer 60 about the axis of the hinged connection offrames gradiometer 110. Depending upon the desired phase alteration, the hinged connection offrames sample tray 50. The hinged connection also permitswindings Helmholtz transformer 60 with respect togradiometer 110. The typical adjustment of phase will include this out-of-parallel orientation, although the other orientation may be preferred in certain circumstances, to accommodate an irregularly shapedsample 200, for example. Noise is applied and adjusted until the noise is 30 to 35 decibels above the molecular electromagnetic emissions sought to be detected. At this noise level, the noise takes on the characteristics of the molecular electromagnetic signal through the well-known phenomenon of stochastic resonance. The stochastic product sought is observed when the oscilloscope trace reflecting the signal detected bygradiometer 110 varies from the trace reflecting the signal directly fromwhite noise generator 80. In alternative embodiments, the signal can be recorded and or processed by any commercially available equipment. - In an alternative embodiment, the method of detecting the molecular electromagnetic signals further comprises injecting noise 180° out of phase with the original noise signal applied at the
Helmholtz transformer 60 through thenoise suppression coil 124 of theSQUID 120. The stochastic product sought can then be observed when the oscilloscope trace reflecting the signal detected bygradiometer 110 becomes non-random. - Regardless of how the noise is injected and adjusted, the stochastic product can also be determined by observing when an increase in spectral peaks occurs. The spectral peaks can be observed as either a line plot on
oscilloscope 160 or as numerical values, or by other well known measuring devices. - Embodiments of the present invention provide a method and apparatus for detecting extremely low-threshold molecular electromagnetic signals without external interference. They further provide for the output of those signals in a format readily usable by a wide variety of signal recording and processing equipment.
- Referring now to
FIG. 7 , an alternative embodiment to the molecular electromagnetic emission detection and processing system of the above Figures is shown. Asystem 700 includes adetection unit 702 coupled to aprocessing unit 704. Although theprocessing unit 704 is shown external to thedetection unit 702, at least a part of the processing unit can be located within the detection unit. - The
detection unit 702, which is shown in a cross-sectional view inFIG. 7 , includes a plurality of components nested or concentric with each other. A sample chamber orfaraday cage 706 is nested within ametal cage 708. Each of thesample chamber 706 and themetal cage 708 can be comprised of aluminum material. Thesample chamber 706 can be maintained in a vacuum and may be temperature controlled to a preset temperature. Themetal cage 708 is configured to function as a low pass filter. - Between the
sample chamber 706 and themetal cage 708 and encircling thesample chamber 706 are a set of parallel heating coils orelements 710. One ormore temperature sensor 711 is also located proximate to theheating elements 710 and thesample chamber 706. For example, four temperature sensors may be positioned at different locations around the exterior of thesample chamber 706. Theheating elements 710 and the temperature sensor(s) 711 are configured to maintain a certain temperature inside thesample chamber 706. - A
shield 712 encircles themetal cage 708. Theshield 712 is configured to provide additional magnetic field shielding or isolation for thesample chamber 706. Theshield 712 can be comprised of lead or other magnetic shielding materials. Theshield 712 is optional when sufficient shielding is provided by thesample chamber 706 and/or themetal cage 708. - Surrounding the
shield 712 is acryogen layer 716 with G10 insulation. The cryogen may be liquid helium. The cryogen layer 716 (also referred to as a cryogenic Dewar) is at an operating temperature of 4 degrees Kelvin. Surrounding thecryogen layer 716 is anouter shield 718. Theouter shield 718 is comprised of nickel alloy and is configured to be a magnetic shield. The total amount of magnetic shielding provided by thedetection unit 702 is approximately −100 dB, −100 dB, and −120 dB along the three orthogonal planes of a Cartesian coordinate system. - The various elements described above are electrically isolated from each other by air gaps or dielectric barriers (not shown). It should also be understood that the elements are not shown to scale relative to each other for ease of description.
- A
sample holder 720 can be manually or mechanically positioned within thesample chamber 706. Thesample holder 720 may be lowered, raised, or removed from the top of thesample chamber 706. Thesample holder 720 is comprised of a material that will not introduce Eddy currents and exhibits little or no inherent molecular rotation. As an example, thesample holder 720 can be comprised of high quality glass or Pyrex. - The
detection unit 702 is configured to handle solid, liquid, or gas samples. Various sample holders may be utilized in thedetection unit 702. For example, depending on the size of the sample, a larger sample holder may be utilized. As another example, when the sample is reactive to air, the sample holder can be configured to encapsulate or form an airtight seal around the sample. In still another example, when the sample is in a gaseous state, the sample can be introduced inside thesample chamber 706 without thesample holder 720. For such samples, thesample chamber 706 is held at a vacuum. Avacuum seal 721 at the top of thesample chamber 706 aids in maintaining a vacuum and/or accommodating thesample holder 720. - A
sense coil 722 and asense coil 724, also referred to as detection coils, are provided above and below thesample holder 720, respectively. The coil windings of the sense coils 722, 724 are configured to operate in the direct current (DC) to approximately 50 kilohertz (kHz) range, with a center frequency of 25 kHz and a self-resonant frequency of 8.8 MHz. The sense coils 722, 724 are in the second derivative form and are configured to achieve approximately 100% coupling. In one embodiment, thecoils coils -
Helmholtz coils shield 712 and themetal cage 708, as explained herein. Each of thecoils coils coils - The degree of coupling between the emissions from the sample and the
coils sample holder 720 relative to thecoils coils sample holder 720. - The
processing unit 704 is electrically coupled to thecoils processing unit 704 specifies the white or Gaussian noise to be injected by thecoils coils - Referring to
FIG. 8 , a processing unit employing aspects of the invention includes asample tray 840 that permits asample 842 to be inserted into, and removed from, aFaraday cage 844 and Helmholtz coil 746. A SQUID/gradiometer detector assembly 848 is positioned within acryogenic dewar 850. A flux-lockedloop 852 is coupled between the SQUID/gradiometer detector assembly 848 and aSQUID controller 854. TheSQUID controller 854 may be a model iMC-303 iMAG multichannel controller provided by Tristan. - An
analog noise generator 856 provides a noise signal (as noted above) to aphase lock loop 858. The x-axis output of the phase lock loop is provided to theHelmholtz coil 846, and may be attenuated, such as by 20 dB. The y-axis output of the phase lock loop is split by asignal splitter 860. One portion of the y-axis output is input the noise cancellation coil at the SQUID, which has a separate input for the gradiometer. The other portion of the y-axis signal isinput oscilloscope 862, such as an analog/digital oscilloscope having Fourier functions like the Tektronix TDS 3000b (e.g., model 3032b). That is, the x-axis output of the phase lock loop drives the Helmholz coil, and the y-axis output, which is in inverted form, is split to input the SQUID and the oscilloscope. Thus, the phase lock loop functions as a signal inverter. The oscilloscope trace is used to monitor the analog noise signal, for example, for determining when a sufficient level of noise for producing non-stationary spectral components is achieved. An analog tape recorder orrecording device 864, coupled to thecontroller 854, records signals output from the device, and is preferably a wideband (e.g. 50 kHz) recorder. APC controller 866 may be an MS Windows based PC interfacing with thecontroller 854 via, for example, anRS 232 port. - In
FIG. 9 , a block diagram of another embodiment of the processing unit is shown. A dual phase lock-inamplifier 202 is configured to provide a first signal (e.g., “x” or noise signal) to thecoils amplifier 202 is configured to lock without an external reference and may be a Perkins Elmer model 7265 DSP lock-in amplifier. This amplifier works in a “virtual mode,” where it locks to an initial reference frequency, and then removes the reference frequency to allow it to run freely and lock to “noise.” - An
analog noise generator 200 is electrically coupled to theamplifier 202. Thegenerator 200 is configured to generate or induce an analog white Gaussian noise at thecoils amplifier 202. As an example, thegenerator 200 may be a model 1380 manufactured by General Radio. - An
impedance transformer 204 is electrically coupled between theSQUID 206 and theamplifier 202. Theimpedance transformer 204 is configured to provide impedance matching between theSQUID 206 andamplifier 202. - The noise cancellation feature of the
SQUID 206 can be turned on or off. When the noise cancellation feature is turned on, theSQUID 206 is capable of canceling or nullifying the injected noise component from the detected emissions. To provide the noise cancellation, the first signal to thecoils SQUID 206 is a noise cancellation signal and is inverted from the first signal at an amplitude sufficient to null the noise at the SQUID output (e.g., 180 degrees out of phase with respect to the first signal). - The
SQUID 206 is a low temperature direct element SQUID. As an example, theSQUID 206 may be a model LSQ/20 LTS dC SQUID manufactured by Tristan Technologies, Inc. Alternatively, a high temperature or alternating current SQUID can be used. Thecoils 722, 724 (e.g., gradiometer) and the SQUID 206 (collectively referred to as the SQUID/gradiometer detector assembly) combined has a magnetic field measuring sensitivity of approximately 5 microTesla/√Hz. The induced voltage in thecoils SQUID 206. The output of theSQUID 206 is a voltage approximately in the range of 0.2-0.8 microVolts. - The output of the
SQUID 206 is the input to aSQUID controller 208. TheSQUID controller 208 is configured to control the operational state of theSQUID 206 and further condition the detected signal. As an example, theSQUID controller 208 may be an iMC-303 iMAG multi-channel SQUID controller manufactured by Tristan Technologies, Inc. - The output of the
SQUID controller 208 is inputted to anamplifier 210. Theamplifier 210 is configured to provide a gain in the range of 0-100 dB. A gain of approximately 20 dB is provided when noise cancellation node is turned on at theSQUID 206. A gain of approximately 50 dB is provided when theSQUID 206 is providing no noise cancellation. - The amplified signal is inputted to a recorder or
storage device 212. Therecorder 212 is configured to convert the analog amplified signal to a digital signal and store the digital signal. In one embodiment, therecorder 212 stores 8600 data points per Hz and can handle 2.46 Mbits/sec. As an example, therecorder 212 may be a Sony digital audiotape (DAT) recorder. Using a DAT recorder, the raw signals or data sets can be sent to a third party for display or specific processing as desired. - A
lowpass filter 214 filters the digitized data set from therecorder 212. Thelowpass filter 214 is an analog filter and may be a Butterworth filter. The cutoff frequency is at approximately 50 kHz. - A
bandpass filter 216 next filters the filtered data sets. Thebandpass filter 216 is configured to be a digital filter with a bandwidth between DC to 50 kHz. Thebandpass filter 216 can be adjusted for different bandwidths. - The output of the
bandpass filter 216 is the input to aFourier transformer processor 218. TheFourier transform processor 218 is configured to convert the data set, which is in the time domain, to a data set in the frequency domain. TheFourier transform processor 218 performs a Fast Fourier Transform (FFT) type of transform. - The Fourier transformed data sets are the input to a correlation and
comparison processor 220. The output of therecorder 212 is also an input to theprocessor 220. Theprocessor 220 is configured to correlate the data set with previously recorded data sets, determine thresholds, and perform noise cancellation (when no noise cancellation is provided by the SQUID 206). The output of theprocessor 220 is a final data set representative of the spectrum of the sample's molecular low frequency electromagnetic emissions. - A user interface (UI) 222, such as a graphical user interface (GUI), may also be connected to at least the
filter 216 and theprocessor 220 to specify signal processing parameters. Thefilter 216,processor 218, and theprocessor 220 can be implemented as hardware, software, or firmware. For example, thefilter 216 and theprocessor 218 may be implemented in one or more semiconductor chips. Theprocessor 220 may be software implemented in a computing device. - This amplifier works in a “virtual mode,” where it locks to an initial reference frequency, and then removes the reference frequency to allow it to run freely and lock to “noise.” The analog noise generator (which is produced by General Radio, a truly analog noise generator) requires 20 dB and 45-dB attenuation for the Helmholz and noise cancellation coil, respectively.
- The Helmholz coil may have a sweet spot of about one cubic inch with a balance of 1/100th of a percent. In an alternative embodiments, the Helmholtz coil may move both vertically, rotationally (about the vertical access), and from a parallel to spread apart in a pie shape. In one embodiment, the SQUID, gradiometer, and driving transformer (controller) have values of 1.8, 1.5 and 0.3 micro-Henrys, respectively. The Helmholtz coil may have a sensitivity of 0.5 Gauss per amp at the sweet spot.
- Approximately 10 to 15 microvolts may be needed for a stochastic response. By injecting noise, the system has raised the sensitivity of the SQUID device. The SQUID device had a sensitivity of about 5 femtotesla without the noise. This system has been able to improve the sensitivity by 25 to 35 dB by injecting noise and using this stochastic resonance response, which amounts to nearly a 1,500% increase.
- After receiving and recording signals from the system, a computer, such as a mainframe computer, supercomputer or high-performance computer does both pre and post processing, such by employing the Autosignal software product by Systat Software of Richmond Calif., for the pre-processing, while Flexpro software product does the post-processing. Flexpro is a data (statistical) analysis software supplied by Dewetron, Inc. The following equations or options may be used in the Autosignal and Flexpro products.
FFT Algorithm:
Best Exact N using Temperton's Prime Factor FFT (C. Temperton, “Implementation of a Self-Sorting In-Place Prime Factor FFT Algorithm, Journal of Computation Physics, v. 58, p. 283, 1985).
Data Tapering Windows: - [cs4 BHarris min] 0.35875−0.48829*cos(2*Pi*i/(n−1))+0.14128*cos(4*Pi*i/(n−1))−0.01168*(6*Pi*i/(n−1)), i=0.n−1
- [Rectangular] No fixed shape tapering available (Oscilloscope)
- Magnitude: sqrt(Re*Re+Im*Im) [Re=real component, Im=imaginary component]
- Amplitude: 2.0*sqrt(Re*Re+Im*Im)/n
- db, decibels: 10.0*log10(Re*Re+Im*Im)
Averaging Replicates: - Replicates are based on the X-values coinciding to within 1e-8 fractional precision.
- Reference Subtraction:
- Reference Signal Subtraction (baseline noise) is performed on Y axis (amplitude) at each point (channel) along the X (time) axis. Negative Y values are then zeroed.
- Cross-Correlation:
- The function calculates the cross correlation function using summation and integration. Since the signal is transient, the correlation function is calculated using direct multiplication and integration. All of the values required for the calculation which lie outside the source channels (data series) are taken to be 0. The points for which t<0 are also calculated.
- Fourier Significance Levels:
- Monte Carlo data is fitted to parametric models. Where data size N is the only factor, univariate TableCurve 2D parametric models are used. For a segmented FFT where segment size and overlap are additional influences, trivariate Chebyshev polynominals are implemented. These are options selected under Autosignal. One could have data sets that analyze individually, or could be analyzed in an overlapping fashion where data set one would be analyzed, then the second half of data set one and the first half of data set two, then data set two, then the second half.
- A flow diagram of the signal detection and processing performed by the
system 100 is shown inFIG. 10 . When a sample is of interest, at least four signal detections or data runs are performed: a first data run at a time t1 without the sample, a second data run at a time t2 with the sample, a third data run at a time t3 with the sample, and a fourth data run at a time t4 without the sample. Performing and collecting data sets from more than one data run increases accuracy of the final (e.g., correlated) data set. In the four data runs, the parameters and conditions of thesystem 100 are held constant (e.g., temperature, amount of amplification, position of the coils, the noise signal, etc.). - At a
block 300, the appropriate sample (or if it's a first or fourth data run, no sample), is placed in thesystem 100. A given sample, without injected noise, emits electromagnetic emissions in the DC-50 kHz range at an amplitude equal to or less than approximately 0.001 microTesla. To capture such low emissions, a white Gaussian noise is injected at ablock 301. - At a
block 302, thecoils - At a
block 304, the injected noise is cancelled as the induced voltage is being detected. This block is omitted when the noise cancellation feature of theSQUID 206 is turned off. - At a
block 306, the voltage values of the data set are amplified by 20-50 dB, depending on whether noise cancellation occurred at theblock 304. And at ablock 308, the amplified data set undergoes analog to digital (A/D) conversion and is stored in therecorder 212. A digitized data set can comprise millions of rows of data. - After the acquired data set is stored, at a block 310 a check is performed to see whether at least four data runs for the sample have occurred (e.g., have acquired at least four data sets). If four data sets for a given sample have been obtained, then lowpass filtering occurs at a
block 312. Otherwise, the next data run is initiated (return to the block 300). - After lowpass filtering (block 312) and bandpass filtering (at a block 314) the digitized data sets, the data sets are converted to the frequency domain at a
Fourier transform block 316. - Next, at a
block 318, like data sets are correlated with each other at each data point. For example, the first data set corresponding to the first data run (e.g., a baseline or ambient noise data run) and the fourth data set corresponding to the fourth data run (e.g., another noise data run) are correlated to each other. IF the amplitude value of the first data set at a given frequency is the same as the amplitude value of the fourth data set at that given frequency, then the correlation value or number for that given frequency would be 1.0. Alternatively, the range of correlation values may be set at between 0-100. Such correlation or comparison also occurs for the second and third data runs (e.g., the sample data runs). Because the acquired data sets are stored, they can be accessed at a later time as the remaining data runs are completed. - When the
SQUID 206 provides no noise cancellation, then predetermined threshold levels are applied to each correlated data set to eliminate statistically irrelevant correlation values. A variety of threshold values may be used, depending on the length of the data runs (the longer the data runs, greater the accuracy of the acquired data) and the likely similarity of the sample's actual emission spectrum to other types of samples. In addition to the threshold levels, the correlations are averaged. Use of thresholds and averaging correlation results in the injected noise component becoming very small in the resulting correlated data set. - If noise cancellation is provided at the
SQUID 206, then the use of thresholds and averaging correlations are not necessary. - Once the two sample data sets have been refined to a correlated sample data set and the two noise data sets have been refined to a correlated noise data set, the correlated noise data set is subtracted from the correlated sample data set. The resulting data set is the final data set (e.g., a data set representative of the emission spectrum of the sample) (block 320).
- Since there can be 8600 data points per Hz and the final data set can have data points for a frequency range of DC-50 kHz, the final data set can comprise several hundred million rows of data. Each row of data can include the frequency, amplitude, phase, and a correlation value.
- In
FIGS. 11A and 11B , there are shown examples of sample emission spectrums. AFourier plot 400 shown inFIG. 11A corresponds to a spectrum of a sample of saturated sodium chloride solution. AFourier plot 500 shown inFIG. 11B corresponds to a spectrum of a sample of an enzyme. - Referring to
FIG. 16 , another alternative embodiment to the systems described above will now be described as asystem 1600. In general, alternatives alternative embodiments described herein are substantially similar to previously described embodiments, and the same reference numbers often identify common elements and functions. Only significant differences in construction or operation are described in detail. - A second derivative gradiometer is shown as 1602, where the target sample is positioned between upper and lower pairs of coils. Two inner coils on opposite sides of the sample complement each other, while two outer coils (top and bottom coils) each complement each other, and oppose the two inner coils. Such an arrangement allows for greater signal extraction from the sample and improved noise rejection.
- While shown in the Figures and described in greater detail below, the
system 1600 employs a concentric series of elements and an arrangement along a central axis extending into the dewar. Astepper motor 1604 allows the sample to be positioned axially within this arrangement of concentric elements. In particular, the sample may be positioned at a desired location within a middle of thegradiometer 1602. - Likewise, a
micrometer adjustment mechanism 1606, such as a mechanical micrometer or stepper motor, allows the Helmholtz coils to be aligned with respect to elements in the system (such as the sample and gradiometer). Such an adjustment of the Helmholtz coil aids in manufacture and calibration of thesystem 1600, as well as allowing precise alignment of fields within the system, such as providing a uniform field with respect to thegradiometer 1602. It may be useful to also provide a field off set or change in field gradient to produce a better stochastic result, to offset noise in the system, or to provide other benefits. -
FIGS. 17A, 17B , and 18 show more clearly the concentric arrangement of elements within thesystem 1600, wherein the sample tube extends axially through a center of a low pass filtering metal shield 1802 (such as a stainless steel alloy) to pass signals below 2 kHz. An outer magnetic (MU) shield surrounds the gradiometer, Helmholtz coils and sample. The arrangement ofsystem 1600 is generally self-explanatory with respect to the Figures. - The random white noise generator, model 1381, manufactured by General Radio and described above, may be replaced by a programmable Gaussian white noise generator manufactured by Noise/Com. Such a generator employs two outputs, one inverted from the other. One output may be connected to the Helmholtz coil, with the other (inverted) output connected to the SQUID noise cancellation coil noted above.
- Likewise, as shown in
FIG. 19 , the Tektronix digital oscilliscope noted above may be replaced by a two-channeleddynamic signal analyzer 1902, model SR 785, manufactured by Stanford Research Systems. Such a signal analyzer may process incoming signals by sampling multiple time domain signals and averaging them across multiple frequency domain FFT's. This may result in a full spectrum frequency domain record of all non-random signal components. Other changes that may be made include replacing the digital audio tape storage system with a digital versatile disk (DVD) recorder 1904. Further, adata acquisition board 1906 manufactured by Keithley,model 3801, may be used, which works with software for generating histograms, as described below. - In the alternative embodiment shown in
FIG. 19 , anoise cancellation coil 1908 is connected between the gradiometer and SQUID. (While a first derivative gradiometer is shown, a second derivative gradiometer, such as that shown inFIG. 16 , may be used.) While not shown inFIG. 19 , an inverted noise channel (inverted with respect to noise applied to the Helmholtz coils) may be applied to the noise cancellation coil 1908 (and may first pass through an impedance transformer that attenuates the noise signal by, for example, 45 dB). In an alternative embodiment, not shown, the noise cancellation coil may be positioned within theSQUID 120, between the SQUID input and output coils. - III. Histogram Method of Generating Spectral Information
-
FIG. 20 is a high level data flow diagram in the histogram method for generating spectral information. Data acquired from the SQUID (box 2002) or stored data (box 2004) is saved as 16 bit WAV data (box 2006), and converted into double-precision floating point data (box 2008). The converted data may be saved (box 2010) or displayed as a raw waveform (box 2012). The converted data is then passed to the algorithm described below with respect toFIG. 21 , and indicated by thebox 2014 labeled Fourier Analysis. The histogram can be displayed at 2016. - With reference to
FIG. 21 , the general flow of the histogram algorithm is to take a discrete sampled time-domain signal and use Fourier analysis to convert it to a frequency domain spectrum for further analysis. The time-domain signals are acquired from an ADC (analog/digital converter) and stored in the buffer indicated at 2102. This sample is SampleDuration seconds long, and is sampled at SampleRate samples per second, thus providing SampleCount (SampleDuration*SampleRate) samples. The FrequencyRange that can be recovered from the signal is defined as half the SampleRate, as defined by Nyquist. Thus, if a time-series signal is sampled at 10,000 samples per second, the FrequencyRange will be 0 Hz to 5 kHz. One Fourier algorithm that may be used is aRadix 2 Real Fast Fourier Transform (RFFT), which has a selectable frequency domain resolution (FFTSize) of powers of two up to 216. An FFTSize of 8192 is selected, to provide provides enough resolution to have at least one spectrum bin per Hertz as long as the FrequencyRange stays at or below 8 kHz. The SampleDuration should be long enough such that SampleCount>(2*)FFTSize* 10 to ensure reliable results. - Since this FFT can only act on FFTSize samples at a time, the program must perform the FFT on the samples sequentially and average the results together to get the final spectrum. If one chooses to skip FFTSize samples for each FFT, a statistical error of 1/FFTSizeˆ0.5 is introduced. If, however, one chooses to overlap the FFT input by half the FFTSize, this error is reduced to 1/(0.81*2*FFTSize)ˆ0.5. This reduces the error from 0.0110485435 to 0.0086805556. Additional information about errors and correlation analyses in general, consult Bendat & Piersol, “Engineering Applications of Correlation and Spectral Analysis”, 1993.
- Prior to performing the FFT on a given window, a data tapering filter may be applied to avoid spectral leakage due to sampling aliasing. This filter can be chosen from among Rectangular (no filter), Hamming, Hanning, Bartlett, Blackman and Blackman/Harris, as examples.
- In an exemplary method, and as shown in
box 2104, we have chosen 8192 for the variable FFTSize, which will be the number of time-domain samples we operate on at a time, as well as the number of discrete frequencies output by the FFT. Note that FFTSize=8192 is the resolution, or number of bins in the range which is dictated by the sampling rate. The variable n, which dictates how many discrete RFFT's (Real FFT's) performed, is set by dividing the SampleCount by FFTSize*2, the number of FFT bins. In order for the algorithm to generate sensible results, this number n should be at least 10 to 20 (although other valves are possible), where more may be preferred to pick up weaker signals. This implies that for a given SampleRate and FFTSize, the SampleDuration must be long enough. A counter m, which counts from 0 to n, is initialized to zero, also as shown inbox 2104. - The program first establishes three buffers:
buffer 2108 for FFTSize histogram bins, that will accumulate counts at each bin frequency;buffer 2110 for average power at each bin frequency, and abuffer 2112 containing the FFTSize copied samples for each m. - The program initializes the histograms and arrays (box 2113) and copies FFTSize samples of the wave data into
buffer 2112, at 2114, and performs an RFFT on the wave data (box 2115). The FFT is normalized so that the highest amplitude is 1 (box 2116) and the average power for all FFTSize bins is determined from the normalized signal (box 2117). For each bin frequency, the normalized value from the FFT at that frequency is added to each bin in buffer 2108 (box 2118). - In
box 2119 the program then looks at the power at each bin frequency, relative to the average power calculated from above. If the power is within a certain factor epsilon (between 0 and 1) of the average power, then it is counted and the corresponding bin is incremented in the histogram buffer at 16. Otherwise it is discarded. - Note that the average power it is comparing to is for this FFT instance only. An enhanced, albeit slower algorithm might take two passes through the data and compute the average over all time before setting histogram levels. The comparison to epsilon helps to represent a power value that is significant enough for a frequency bin. Or in broader terms, the equation employing epsilon helps answer the question, “is there a signal at this frequency at this time?” If the answer is yes, it could due be one of two things: (1) stationary noise which is landing in this bin just this one time, or (2) a real low level periodic signal which will occur nearly every time. Thus, the histogram counts will weed out the noise hits, and enhance the low level signal hits. So, the averaging and epsilon factor allow one to select the smallest power level considered significant.
- Counter m is incremented at
box 2120, and the above process is repeated for each n set of WAV data until m is equal to n (box 2121). At each cycle, the average power for each bin is added to the associated bin at 2118, and each histogram bin is incremented by one when the power amplitude condition at 2114 is met. - When all n cycles of data have been considered, the average power in each bin is determined by dividing the total accumulated average power in each bin by n, the total number of cycles (box 2122) and the results displayed (box 2123). Except where structured noise exists, e.g., DC=0 or at multiples of 60 Hz, the average power in each bin will be some relatively low number. This is indicated in the plots shown at FIGS. 22A-D (the histograms generated at 400, 600, 700, and 900 mV). The plots of
FIGS. 22A-22D show only a portion of the histogram bins, namely a spectrum from 7953 Hz through 8533 Hz. As shown inFIGS. 22A and 22B , no stochastic event is visible at 400 mV or 600 mV of injected noise, respectively. However, as shown inFIG. 22C , at 700 mV, a visible stochastic event is evident. Thereafter, as shown inFIG. 22D , at 900 mV, the stochastic event is lost. - The histogram produced by the above steps contains, in each bin, a count between 0 and n of the number of times that the power at that frequency was above (epsilon*the average power for that whole FFT output). If a bin count is incremented due to unstructured noise, that noise will be distributed across all the frequency bins over time, thus not adding up to much in a given bin. If there is consistent signal at a given frequency, it will be present at each of the n time slices and thus have a bin count approaching n. Large amplitude noise, such as sixty hertz and its harmonics have both high bin counts as well as high average power. We can differentiate between these frequencies, and the ones we are interested in that have low average power, but high bin counts.
-
FIGS. 22A-22D show histograms generated by the method at four different noise power inputs. As shown, the program may display average power at each frequency as a vertical bar. The histogram bin counts may be represented as a connected upper line. If the power is considered “low” (e.g. less than average/3), and the histogram has a certain count, then a connecting line may become observable between the peak of a power bar and a peak of a histogram bar. Bins highlighted by the connecting lines are likely candidates for low energy molecular spectra. - It can be appreciated from
FIGS. 22A-22D and from the above, that there are two settings of note used in generating a meaningful histogram, that is, a histogram that shows stochastic resonance effects related to a sample being interrogated. The first is the power level of Gaussian white noise supplied to the sample. If this level is too low, the noise level is not sufficient to create stochastic resonance and the bin histogram reflects noise only. If the power input is too high, the average power level calculated for each bin is high and stochastic events cannot be distinguished. - The second setting is the value of epsilon. This value determines a power value that will be used to distinguish an event over average value. At a value of 1, no events will be detected, since power will never be greater than average power. As epsilon approaches zero, virtually every value will be placed in a bin. Between 0 and 1, and typically at a value that gives a number of bin counts between about 20-50% of total bin counts for structured noise, epsilon will have a maximum “spectral character,” meaning the stochastic resonance events will be most highly favored over pure noise.
- Therefore, in practicing the invention, one can systematically increase the power gain on the noise input, e.g., in 100 mV increments between 0 and 1 V, and at each power setting, adjust epsilon until a histogram having well defined peaks is observed. Where, for example, the sample being processed represents a 20 second time interval, total processing time for each different power and epsilon will be about 25 seconds. When a well-defined signal is observed, either the power setting or epsilon or both can be refined until an optimal histogram, meaning one with the largest number of identifiable peaks, is produced.
- Under this algorithm, numerous bins may be filled and associated histogram rendered for low frequencies due to the general occurrence of noise (such as environmental noise) at the low frequencies. Thus, the system may simply ignore bins below a given frequency (e.g., below 1 kHz), but still render sufficient bin values at higher frequencies to determine unique signal signatures between samples.
- Alternatively, since a purpose of the epsilon variable is to accommodate different average power levels determined in each cycle, the program could itself automatically adjust epsilon using a predefined function relating average power level to an optimal value of epsilon.
- Similarly, the program could compare peak heights at each power setting, and automatically adjust the noise power setting until optimal peak heights or character is observed in the histograms.
- Although the value of epsilon may be a fixed value for all frequencies, it is also contemplated to employ a frequency-dependent value for epsilon, to adjust for the higher value average energies that may be observed at low frequencies, e.g., DC to 1,000. A frequency-dependent epsilon factor could be determined, for example, by averaging a large number of low-frequency FFT regions, and determining a value of epsilon that “adjusts” average values to values comparable to those observed at higher frequencies.
- Referring to
FIGS. 23A-23C , an example of a user interface for generating histograms is shown. Aslider bar 2302 determines the length of a sample waveform segment, such as up to 300-600 seconds, and allows a user to effectively scroll within a waveform. Abox 2304 allows the user to set a Nyquist frequency, such as 5, 10 or 20 kHz, and also provided is an adjacent reset button. Aslider bar 2306 allows the user to move the baseline for histograms, while a 60Hz checkbox 2308 allows the user to identify the 60 Hz bin and all related 60 Hz harmonics with vertical lines (as shown inFIG. 23C ). When anacquire button 2312 is selected, the software generates or acquires a waveform from a sample, such as that shown inFIG. 23B . When anfft button 2310 is selected, the software generates a histogram plot, such as that shown inFIG. 23C . - IV. Methods and Applications
- This section describes the use of the apparatus described above for interrogating a sample, and a variety of applications of the apparatus in characterizing a sample and in detecting sample components. Also disclosed, in accordance with the invention, is a low-frequency spectroscopic signature or data set by which a sample can be characterized, and a time-domain signal of a sample, used, for example, in generating the sample spectroscopic signature.
- A. Method of Interrogating a Sample
- An objective of the method of the invention is to generate spectroscopic information relating a sample being interrogated. As will be seen, the information may be in the form of a spectral plot, in a selected low-frequency spectral range, or a data set which identifies low-frequency spectral components characterizing the sample, or actual identification of a sample or sample components, based on the characteristic frequencies identified for the sample.
- The sample may be any material having atomic or molecular components, e.g., ionic salt components or molecular compound in ionized or nonionized form, or protonated or non-protonated form, that has molecular rotation, and preferably a dipole moment such that molecular rotation in a magnetic field, e.g., the earth's magnetic field, is effective to produce a low-frequency electromagnetic emission. The sample is typically a liquid sample, but may be gaseous or solid or semi-solid as well, as long as at least one component of the sample has one or more rotational degrees of freedom. Typical samples are aqueous or organic solutions having one or more solute components, which may be the sample material of interest, dissolved in the solvent.
- The sample is placed in a suitable vessel, preferably one such as Pyrex glass that has little observable low-frequency spectral components, and the vessel is then positioned in the apparatus container as described in Section II. With the sample positioned in the apparatus container, the Gaussian noise generator is activated to inject Gaussian noise into the sample. The amplitude (mean amplitude) of Gaussian noise injected is preferably sufficient to produce non-stationary composite time-domain signal components. This may be done, for example, using an oscilloscope with a Fourier transform capability, and observing the frequency-domain signal in a suitable range, e.g., 200-800 Hz window. A suitable noise level is selected when detectable frequency components are first observed.
- During noise injection, the recording device records a time-domain electromagnetic signal from the detector over a preset time interval. The recording interval may be relatively short, e.g., 30-60 seconds, or may be several minutes or more, depending on the final spectral resolution required. The signals recorded are stored in a suitable signal storage device, e.g., a tape or hard disc, for use in later signal processing operations now to be described.
- In general, it is desirable to enhance sample signal components by cross-correlating the sample time-domain signal recorded with a second time-domain signal of the same sample or, less preferably, an identical sample or a sample having the same sample components of interest. The recording time for the second signal is preferably the same as for the first signal. The two signals are cross-correlated using a standard cross-correlation algorithm in the time domain. This results in a spreadsheet or spectrum identifying the signal spectral components that are common in both signals that hold up over time, and a correlation value for each component which measures the relationship between spectral components common to both signals.
- The improvement in spectral resolution obtained by the signal cross-correlation is seen in the
FIGS. 12A and 12B , andFIGS. 13A and 13B . The Figures are Fast Fourier transforms of a first time-dependent signal in the frequency domain (FIGS. 12A and 13A ) or a fast Fourier transform of the first and second cross-correlated frequency-domain spectra (the spreadsheet referred to above) (FIGS. 13B and 13B ) to plot the spectral components in the frequency domain, and in the spectral range of 500-530 Hz. - Comparing
FIGS. 12A and 12B for a NaCl sample, it is seen that the cross-correlation signal processing significantly enhances signal-to-noise ratio, brings out much more detail in the sample-specific spectral component whose peak is at 522.5 Hz, and also produces a significantly refined peak position. Similar (and exemplary) results were observed for the alkyl ether sulfate sample, whose spectral features in the 500-530 Hz range are seen inFIGS. 13A and 13B for uncorrelated and correlated signals, respectively. As with the NaCl sample, the spectrum derived from the correlated signals gave much lower signal-to-noise ratio, much more detail and information as to sample-specific spectral components. The signal correlation can also be applied, conventionally, to produce a spreadsheet relating frequency and phase (rather than frequency and amplitude). - The correlated time-domain spectrum from above is plotted in the frequency domain by applying a fast Fourier transform to the spectrum, where the spectral correlations values are represented in the y axis as amplitudes. The plot is within the frequency range DC to 50 kHz, preferably in the region DC to 6.5 kHz. As will be seen below, the dominant spectral features of many samples are found in the 100-1,500 Hz range, particularly the 500-550 Hz range; and thus the spectrum generated may be confined accordingly, e.g., in the 500-530 Hz range. The FTT is carried out by a well known FTT algorithm. The correlated time-domain signal may also, or alternatively, be transformed to phase-domain or amplitude or magnitude domain signals, to extract signal information related to phase or amplitude components of the sample spectrum.
- Once the frequency-domain spectrum is generated, either by the cross-correlating or FFT step, the spectrum is used to identify one or more low-frequency signal components that are characteristic of the sample being interrogated. This step may be performed by the user from direct viewing or by computer analysis of the spectra.
-
FIGS. 14A-14F show spectral features for the samples deionized water (FIG. 14A ), saturated NaCl (FIG. 14B ), 1% NaCl in deionized water (FIG. 14C ), saturated NaBr (FIG. 14D ), alkyl ether sulfate (FIG. 14E ), and empty sample vessel (FIG. 14F ), all in the spectral range of about 500-530. As seen, each sample has distinctive spectral components characterized by one or more peaks at well defined frequencies. - B. Characterizing a Sample
- Accordance with another aspect of the invention, the method above is used to generate a data set of low-frequency spectral components of a given sample, also referred to herein as a low-frequency signature signal of the material.
- The 500-530 spectral range shown for the several samples above was selected to illustrate having prominent spectral features in the various samples. In order to obtain a more complete data set of spectral frequency components of a sample, spectral components over a wider frequency range, e.g., 100-1,500 Hz, should be determined. In one aspect, the invention includes a data set of spectral components associated with a given sample material, e.g., a solvent, gas, or solute component of a solution. The data set includes a list of the low-frequency spectral components of the sample, e.g., in the 100-1,500 Hz range, whose cross-spectral correlations have a selected statistical measure above background spectral noise, or selected ones of these components that are unique to the sample.
- A variety of signal-analysis methods may be employed in generating the low-frequency data set for a given sample. In one exemplary method, a cross-correlated sample signal spectrum is compared with a cross-correlated noise (no sample) signal. The algorithm next advances incremental, e.g., in 0.1 Hz intervals across the cross-correlated sample spectrum and the cross-correlated noise spectrum, looking at the correlation value at each frequency point, and subtracts the noise correlation from the sample correlation at that point, to yield a frequency plot of corrected correlation values. These values will be relative to a particular sample, and depend, for example, on the relative amplitude of any noise component.
- In general, frequency components having a higher correlation value (relative to the other values in that sample) will tend to hold up (be observed) over many interrogations of the same sample. To identify those that do hold up, the frequency components observed for the sample over two or more sample sets, each obtained as above, are compared, and only those that are seen in two (or more, if available) sets are taken as valid components of the data set for that sample. In the tables below, data sets for several samples (as identified in the tables) are given along with the correlations determined from a single sample interrogation. Those values indicated in italics (typically having the smaller correlation values) were found not to hold up in multiple data sets from the same sample material.
- Thus, for example, for the saturated NaCl sample in Table 1, spectral components at 522.58, 523.12, 523.47, and 523.85 Hz correlate from sample to sample, and would form a data set for the sample in the frequency range 500-530 Hz. Additional members of the data set may be included in an expanded frequency range.
- Similarly, for the amino acid sample of Table 3, the data set would include components at 262.93, 257.81, 257.23, 536.68, 448.05, 531.37, 528.80, 593.44, 588.68, 583.74, 578.61, 769.59, and 744.14 in the frequency range of between about 250 and 1,400 Hz. The greater spectral composition of the amino acid sample, relative to NaCl, presumably reflects in part, the greater complexity of the sample molecule.
Reversing Noise NaCl (Sat) NaCl (1%) NaBr (Sat) Frequency / Frequency / Frequency / Correlation Correlation Correlation √ A √ A √ A 522.58 .3762 521.12 1.4161 520.57 2.0847 523.12 .1732 521.48 1.4100 519.84 2.0704 523.47 .1235 515.99 1.3865 509.37 2.0304 523.85 .1021 520.75 1.3641 513.45 2.0155 507.38 .0832 514.34 1.3735 516.35 1.9950 524.43 .0768 525.86 1.3440 519.46 1.9950 512.71 .0753 523.70 1.3400 518.33 1.9929 — — 526.61 1.3364 522.78 1.9635 -
Reversing Noise DI Water Alkyl Ether Sulfate Noise Frequency / Frequency / Frequency / Correlation Correlation Correlation √ A √ A √ A 521.12 1.5324 517.81 .3376 514.34 .0734 521.67 1.0818 516.50 .3375 513.79 .0432 520.20 1.0630 517.08 .2776 506.28 .0326 511.23 1.0502 515.46 .2749 512.70 .0277 515.44 1.0457 518.37 .2508 522.58 .0220 513.06 1.0451 519.47 .2425 525.15 .0177 525.51 1.0371 515.44 .2400 516.36 .0149 520.75 1.0301 519.84 .2383 523.13 .0140 -
Spectra Frequency Correlation Number (Hz) Factor 1 262.93 .139 2 340.39 .134 3 257.81 .126 4 357.23 .114 5 417.48 .110 6 536.68 .101 7 448.05 .096 8 531.37 .096 9 528.80 .077 10 593.44 .071 11 588.68 .065 12 583.74 .058 13 1408.99 .052 14 840.08 .050 15 1393.99 .048 16 578.61 .045 17 1348.99 .044 18 769.59 .042 19 1042.96 .042 20 1238.52 .042 21 1472.16 .042 22 1062.92 .041 23 1281.73 .041 24 744.14 .039 - The data above demonstrates that both simple and more complex molecular samples can be characterized in terms of unique low-frequency spectral compounds. The data set associated with a given sample material may also include (as shown in the tables) the associated correlations values of the spectral components. The data set may be used for example, in identifying components in an unknown sample and/or for estimating the relative concentrations of a material in a sample. The use of the method for identifying low-concentration components in a sample is discussed in the next section.
- C. Identifying Components in a Sample
- It is often desirable to detect sample components, e.g., trace contaminants, present in a multi-component sample material, such as a liquid sample with unknown contaminants, or other samples capable of holding or supporting a contaminant that it is desired to detect.
- An analytical method for detecting a component of a sample, in accordance with another aspect of the invention, includes first identifying the low-frequency sample spectral components of a sample (i) in a selected frequency range between DC and 50 kHz (ii) whose cross-spectral correlations have a selected statistical measure above background spectral noise, as described above.
- The sample spectral components are then compared with characteristic low-frequency spectral components of known compounds suspected of being present in the sample. In a typical example, the sample components are compared against the data set of each of the components suspected of being in the sample and which one desires to detect. A components, e.g., compound is identified as being present in the sample if one or more of its characteristic low-frequency spectral components correspond to one or more low-frequency spectral components of a known sample.
- As shown in the set of
FIGS. 11A-11F , detection of a compound (an amino acid) can occur at very low levels, e.g., in the parts per billion range or lower. In particular, even at a dilution of 1:10 billion w/v, a characteristic spectral component at about 531 Hz is observed. The Figures demonstrate that signal amplitude, corresponding spectral component correlation, does decline with increasing compound dilution. However, the loss in signal amplitude at low concentration can be compensated for by extending the recording time, in this example, from 50 seconds for the first group of Figures to 4.25 minutes for the most dilute sample (FIG. 11F ). - Where, as in the above example, the spectral component amplitude declines with decreasing concentration, the amount of compound can be estimated on the basis of signal amplitude, assuming that the data set for the compound also includes concentration dependent amplitude information.
- It has also been observed in some cases that the frequency of the characteristic spectral components may shift by as much as 3 Hz in a systematic fashion with changes in concentration. For such compounds, the amount of material present in a sample can be estimated by changes in amplitude and/or frequency shift in one or more of the spectral components. It will be appreciated for materials showing a concentration dependent frequency shift that a data set for that compound could include concentration-dependent frequencies as well as concentration-dependent amplitudes for particular components.
- D. Time-Domain Signals
- In still another aspect, the invention includes a time-domain signal associated with a material of interest. The time-domain signal, and its method of production, have been discussed above. Briefly, the signal is produced by placing the sample of interest in a container having both magnetic and electromagnetic shielding, injecting Gaussian noise into the sample; and recording an electromagnetic time-domain signal composed of sample source radiation superimposed on the injected Gaussian noise.
- The signal may be used to characterize a sample, much as a spectral component data set is used to characterize a material. Alternatively, the signal may be used for generating a low frequency signal signature of spectral components associated with a material of interest. The signal signature can be generated, also as described above, by (i) cross correlating the time-domain signal recorded with a second time domain signal separately recorded from the same or similar sample, to produce a frequency domain spectrum in a frequency range within DC to 50 kHz.
- Unless the context clearly requires otherwise, throughout the description and the claims, the words “comprise,” “comprising,” and the like are to be construed in an inclusive sense as opposed to an exclusive or exhaustive sense; that is to say, in the sense of “including, but not limited to.” Words in the above detailed description using the singular or plural number may also include the plural or singular number respectively. Additionally, the words “herein,” “above,” “below” words of similar import, when used in this application, shall refer to this application as a whole and not to any particular portions of this application. When the claims use the word “or” in reference to a list of two or more items, that word covers all of the following interpretations of the word: any of the items in the list, all of the items in the list and any combination of the items in the list.
- The above detailed descriptions of embodiments of the invention are not intended to be exhaustive or to limit the invention to the precise form disclosed above. While specific embodiments of, and examples for, the invention are described above for illustrative purposes, various equivalent modifications are possible within the scope of the invention, as those skilled in the relevant art will recognize. For example, while processes or steps are presented in a given order, alternative embodiments may perform routines having steps in a different order, and some steps may be deleted, moved, added, subdivided, combined, and/or modified. Each of these steps may be implemented in a variety of different ways. Also, while these steps are shown as being performed in series, these steps may instead be performed in parallel, or may be performed at different times.
- The teachings of the invention provided herein can be applied to other systems, not necessarily the system described herein. These and other changes can be made to the invention in light of the detailed description. The elements and acts of the various embodiments described above can be combined to provide further embodiments.
- All of the above patents and applications and other references, including any that may be listed in accompanying filing papers, are incorporated herein by reference. Aspects of the invention can be modified, if necessary, to employ the systems, functions and concepts of the various references described above to provide yet further embodiments of the invention.
- These and other changes can be made to the invention in light of the above detailed description. While the above description details certain embodiments of the invention and describes the best mode contemplated, no matter how detailed the above appears in text, the invention can be practiced in many ways. Details of the signal acquiring and analysis system may vary considerably in its implementation details, while still be encompassed by the invention disclosed herein. As noted above, particular terminology used when describing certain features or aspects of the invention should not be taken to imply that the terminology is being re-defined herein to be restricted to any specific characteristics, features or aspects of the invention with which that terminology is associated. In general, the terms used in the following claims should not be construed to limit the invention to the specific embodiments disclosed in the specification, unless the above Detailed Description section explicitly defines such terms. Accordingly, the actual scope of the invention encompasses not only the disclosed embodiments, but also all equivalent ways of practicing or implementing the invention under the claims.
- While certain aspects of the invention are presented below in certain claim forms, the inventors contemplate the various aspects of the invention in any number of claim forms. For example, while only one aspect of the invention is recited as embodied in a method claim format, it may likewise be embodied in a computer-readable medium claim format. Accordingly, the inventors reserve the right to add additional claims after filing the application to pursue such additional claim forms for other aspects of the invention.
Claims (17)
1-22. (canceled)
23. Apparatus for interrogating a sample that exhibits low-frequency molecular motion, comprising:
a container adapted for receiving the sample, the container having both magnetic and electromagnetic shielding;
an adjustable-power source of Gaussian noise for directing Gaussian noise to the sample, with the sample in the container;
a detector for detecting an electromagnetic time-domain signal composed of sample source radiation superimposed with the directed Gaussian noise; and
an electronic computer adapted to receive the time-domain signal from the detector, and to process the signal to generate a spectral plot that displays, at a selected power setting of the Gaussian noise source, low-frequency spectral components characteristic of the sample in a selected frequency range between DC and 50 kHz.
24. The apparatus of claim 23 , wherein the electronic computer includes a signal analyzer that functions to (i) calculate a series of Fourier spectra of the time-domain signal over each of a plurality of defined time periods, in a selected frequency range between 100 Hz and 50 kHz, and (ii) average the Fourier spectra.
25. The apparatus of claim 24 , wherein the calculating includes calculating at least five Fourier spectra, each taken over a 1-5 second time-domain interval.
26. The apparatus of claim 23 , wherein the source of Gaussian noise includes an adjustable-power Gaussian noise generator and a Helmholz coil which is contained within the magnetic electromagnetic shielding, and which receives a selected noise output signal from the noise generator in a range 100 mV to 1 V.
27. The apparatus of claim 26 , wherein the generator is designed to inject Gaussian noise into the sample at a frequency between DC and 2 kHz.
28. The apparatus of claim 23 , wherein the detector is a second-derivative gradiometer which outputs a current signal, and a SQUID operatively connected to the gradiometer to convert the current signal to an amplified voltage signal.
29. A method for interrogating a sample that exhibits low-frequency molecular motion, comprising:
placing the sample in a container having both magnetic and electromagnetic shielding,
injecting noise into the sample at a selected noise amplitude;
recording an electromagnetic time-domain signal composed of sample source radiation superimposed on the injected noise,
generating a spectral plot that contains, at a selected power setting of the noise source, low-frequency, sample-dependent spectral components characteristic of the sample in a selected frequency range between 100 and 50 kHz, and
repeating the injecting, recording and generating at different selected noise amplitudes until a plot showing a maximum or near maximum number of spectral components characteristic of the sample is generated.
30. The method of claim 29 , wherein the generating includes (i) calculating a series of Fourier spectra of the time-domain signal over each of a plurality of defined time periods, in a selected frequency range between 100 Hz and 50 kHz, and (ii) averaging the Fourier spectra.
31. A method for interrogating a sample that exhibits low-frequency molecular motion, comprising:
receiving the sample;
magnetically and electromagnetically shielding the sample;
directing adjustable-power Gaussian noise to the sample;
detecting an electromagnetic time-domain signal composed of sample source radiation superimposed with the directed Gaussian noise; and
receiving the time-domain signal, and processing the signal to generate a spectral plot that displays, at a selected power setting of the Gaussian noise source, low-frequency spectral components characteristic of the sample in a selected frequency range between DC and 50 kHz.
32. The method of claim 31 , wherein the processing includes (i) calculating a series of Fourier spectra of the time-domain signal over each of a plurality of defined time periods, in a selected frequency range between 100 Hz and 50 kHz, and (ii) averaging the Fourier spectra.
33. The method of claim 32 , wherein the calculating includes calculating at least five Fourier spectra, each taken over a 1-5 second time-domain interval.
34. The method of claim 31 , wherein the directing Gaussian noise includes providing a Helmholz coil contained within magnetic electromagnetic shielding, which receives a selected noise output signal in a range 100 mV to 1 V.
35. The method of claim 34 , wherein the directing Gaussian noise includes injecting Gaussian noise into the sample at a frequency between DC and 2 kHz.
36. The method of claim 31 , wherein the detecting is performed using a second-derivative gradiometer which outputs a current signal, and a SQUID operatively connected to the gradiometer to convert the current signal to an amplified voltage signal.
37. An apparatus for interrogating a sample that exhibits low-frequency molecular motion, comprising:
means for placing the sample in a container having both magnetic and electromagnetic shielding;
means for injecting noise into the sample at a selected noise amplitude;
means for recording an electromagnetic time-domain signal composed of sample source radiation superimposed on the injected noise;
means for generating a spectral plot that contains, at a selected power setting of the noise source, low-frequency, sample-dependent spectral components characteristic of the sample in a selected frequency range between 100 and 50 kHz; and,
means for repeating the injecting, recording and generating at different selected noise amplitudes until a plot showing a maximum or near maximum number of spectral components characteristic of the sample is generated.
38. The apparatus of claim 37 , wherein the means for generating includes means for calculating a series of Fourier spectra of the time-domain signal over each of a plurality of defined time periods, in a selected frequency range between 100 Hz and 50 kHz, and means for averaging the Fourier spectra.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/349,857 US20060158183A1 (en) | 2002-03-29 | 2006-02-07 | System and method for characterizing a sample by low-frequency spectra |
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/112,927 US6724188B2 (en) | 2002-03-29 | 2002-03-29 | Apparatus and method for measuring molecular electromagnetic signals with a squid device and stochastic resonance to measure low-threshold signals |
PCT/US2003/009544 WO2003083439A2 (en) | 2002-03-29 | 2003-03-28 | System and method for characterizing a sample by low-frequency spectra |
PCT/US2003/011834 WO2003102566A2 (en) | 2002-04-19 | 2003-04-18 | System and method for sample detection based on low-frequency spectral components |
US10/683,875 US6995558B2 (en) | 2002-03-29 | 2003-10-09 | System and method for characterizing a sample by low-frequency spectra |
US11/349,857 US20060158183A1 (en) | 2002-03-29 | 2006-02-07 | System and method for characterizing a sample by low-frequency spectra |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/683,875 Continuation US6995558B2 (en) | 2002-03-29 | 2003-10-09 | System and method for characterizing a sample by low-frequency spectra |
Publications (1)
Publication Number | Publication Date |
---|---|
US20060158183A1 true US20060158183A1 (en) | 2006-07-20 |
Family
ID=34435399
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/683,875 Expired - Lifetime US6995558B2 (en) | 2002-03-29 | 2003-10-09 | System and method for characterizing a sample by low-frequency spectra |
US10/575,182 Abandoned US20070210790A1 (en) | 2003-10-09 | 2004-10-08 | System And Method For Characterizing A Sample By Low-Frequency Spectra |
US11/349,857 Abandoned US20060158183A1 (en) | 2002-03-29 | 2006-02-07 | System and method for characterizing a sample by low-frequency spectra |
Family Applications Before (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/683,875 Expired - Lifetime US6995558B2 (en) | 2002-03-29 | 2003-10-09 | System and method for characterizing a sample by low-frequency spectra |
US10/575,182 Abandoned US20070210790A1 (en) | 2003-10-09 | 2004-10-08 | System And Method For Characterizing A Sample By Low-Frequency Spectra |
Country Status (8)
Country | Link |
---|---|
US (3) | US6995558B2 (en) |
EP (1) | EP1671144A4 (en) |
JP (2) | JP4425922B2 (en) |
CN (1) | CN1864073B (en) |
AU (1) | AU2004280998B2 (en) |
BR (1) | BRPI0415235B1 (en) |
CA (1) | CA2538988C (en) |
WO (1) | WO2005036131A2 (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100277164A1 (en) * | 2006-09-01 | 2010-11-04 | Commonwealth Scientific And Industrial Research Organisation | Method and apparatus for signal recovery |
WO2014145114A1 (en) * | 2013-03-15 | 2014-09-18 | Aliphcom | Dynamic control of sampling rate of motion to modify power consumption |
US9417257B2 (en) | 2004-07-27 | 2016-08-16 | Nativis, Inc. | System and method for collecting, storing, processing, transmitting and presenting very low amplitude signals |
US10046172B2 (en) | 2013-03-15 | 2018-08-14 | Nativis, Inc. | Controller and flexible coils for administering therapy, such as for cancer therapy |
CN109686528A (en) * | 2018-12-18 | 2019-04-26 | 武汉船用电力推进装置研究所(中国船舶重工集团公司第七一二研究所) | A kind of HTS-SMES magnet device |
US11097623B2 (en) | 2014-09-24 | 2021-08-24 | Rohm Co., Ltd. | Current mode control type switching power supply device |
Families Citing this family (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6995558B2 (en) * | 2002-03-29 | 2006-02-07 | Wavbank, Inc. | System and method for characterizing a sample by low-frequency spectra |
US6724188B2 (en) * | 2002-03-29 | 2004-04-20 | Wavbank, Inc. | Apparatus and method for measuring molecular electromagnetic signals with a squid device and stochastic resonance to measure low-threshold signals |
CA2460794C (en) | 2002-04-19 | 2005-02-08 | Bennett M. Butters | System and method for sample detection based on low-frequency spectral components |
GB0424717D0 (en) * | 2004-11-09 | 2004-12-08 | Metis Instr & Equipment N V | Sensor for measuring magnetic flux |
US7620509B2 (en) * | 2005-11-03 | 2009-11-17 | Tektronix, Inc. | Detection of time-frequency codes using a spectrogram |
FR2894673B1 (en) * | 2005-12-14 | 2014-10-31 | Luc Montagnier | METHOD FOR CHARACTERIZING A BIOCHEMICAL ELEMENT HAVING BIOLOGICAL ACTIVITY BY ANALYZING LOW FREQUENCY ELECTROMAGNETIC SIGNALS |
JP4823726B2 (en) * | 2006-03-17 | 2011-11-24 | 国立大学法人豊橋技術科学大学 | Biopolymer detection method and detection apparatus thereof |
EP2030209B1 (en) * | 2006-05-19 | 2013-06-26 | ABB Technology Ltd. | Reactor shield |
JP4399610B2 (en) * | 2006-09-26 | 2010-01-20 | 国立大学法人九州工業大学 | Magnetic property measuring method and system |
US8150037B2 (en) * | 2007-02-20 | 2012-04-03 | Carnegie Mellon University | Apparatus and method for secure, user-friendly deployment of information |
JP2010525892A (en) * | 2007-05-04 | 2010-07-29 | カリフォルニア インスティテュート オブ テクノロジー | Low magnetic field SQUID-MRI apparatus, component and method |
US8647273B2 (en) * | 2007-06-21 | 2014-02-11 | RF Science & Technology, Inc. | Non-invasive weight and performance management |
US10264993B2 (en) * | 2007-06-21 | 2019-04-23 | Rf Science & Technology Inc. | Sample scanning and analysis system and methods for using the same |
US8259299B2 (en) * | 2007-06-21 | 2012-09-04 | Rf Science & Technology Inc. | Gas scanning and analysis |
US8382668B2 (en) * | 2007-06-21 | 2013-02-26 | Rf Science & Technology Inc. | Non-invasive determination of characteristics of a sample |
US8647272B2 (en) * | 2007-06-21 | 2014-02-11 | Rf Science & Technology Inc | Non-invasive scanning apparatuses |
US8022864B2 (en) * | 2007-11-08 | 2011-09-20 | The United States Of America As Represented By The Secretary Of Commerce | Detection of transient signals in doppler spectra |
US8405379B1 (en) | 2008-09-18 | 2013-03-26 | Luc Montagnier | System and method for the analysis of DNA sequences in biological fluids |
KR101012107B1 (en) | 2009-04-22 | 2011-02-07 | 한국표준과학연구원 | Data acquisition system of multi-channel SQUID signal |
TWI407121B (en) * | 2010-06-15 | 2013-09-01 | Jiann Fuh Chen | Acoustic emission detection device for discharge in high-pressure cable |
US8942783B2 (en) * | 2011-06-09 | 2015-01-27 | Element 1 Systems Llc | Stochastic resonance and brownian motion for the reduction of sudden infant death syndrome (SIDS) |
US9965250B2 (en) * | 2011-12-07 | 2018-05-08 | Quintessencelabs Pty Ltd. | Integrated quantum-random noise generator using quantum vacuum states of light |
BR112015013573B1 (en) * | 2012-12-10 | 2020-11-03 | Arcelormittal Investigacion Y Desarrollo, S.L. | test method of a austenitic steel reformer tube. |
US10094897B2 (en) * | 2014-07-25 | 2018-10-09 | Gersh Z. Taicher | Nuclear magnetic resonance apparatus and methods |
JP6357391B2 (en) * | 2014-09-22 | 2018-07-11 | 日本電子株式会社 | Information processing apparatus and information processing method |
CN104300226B (en) * | 2014-09-29 | 2017-05-03 | 北京智谷睿拓技术服务有限公司 | Method and device for realizing directional antenna and directional antenna |
CN105911487B (en) * | 2016-04-14 | 2018-12-14 | 中国科学院上海微系统与信息技术研究所 | A kind of superconducting magnetic sensor search coil and detector |
CN106199680A (en) * | 2016-09-13 | 2016-12-07 | 清华大学 | Slow neutron detection device |
JP7188907B2 (en) * | 2018-04-27 | 2022-12-13 | ラピスセミコンダクタ株式会社 | OBJECT DETERMINATION DEVICE, PROGRAM, OBJECT DETERMINATION METHOD AND SEMICONDUCTOR DEVICE |
US10088535B1 (en) | 2018-06-06 | 2018-10-02 | QuSpin, Inc. | System and method for measuring a magnetic gradient field |
EP4153296A1 (en) | 2020-05-22 | 2023-03-29 | Nearfield Atomics Inc. | Systems for therapy delivery using near field magnetic induction devices |
RU2747100C1 (en) * | 2020-07-07 | 2021-04-26 | Федеральное государственное бюджетное научное учреждение "Федеральный исследовательский центр "Красноярский научный центр Сибирского отделения Российской академии наук" (ФИЦ КНЦ СО РАН, КНЦ СО РАН) | Ferromagnetic resonance scanning spectrometer |
US20220274604A1 (en) | 2021-02-26 | 2022-09-01 | NearField Atomics Inc. | Systems and methods to keep drivers, pilots, or others alert or focused by electromagnetic delivery of signals |
WO2022203904A1 (en) | 2021-03-22 | 2022-09-29 | NearField Atomics Inc. | Systems and methods for measuring magnetic fields from solvated target molecules using a magnetoresistive sensor |
CN117272022A (en) * | 2023-09-19 | 2023-12-22 | 小谷粒(广州)母婴用品有限公司 | Detection method of MEMS oscillator |
Citations (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4031462A (en) * | 1975-07-07 | 1977-06-21 | Motorola, Inc. | Frequency spectrum analyzer |
US4095168A (en) * | 1977-02-22 | 1978-06-13 | Varian Associates, Inc. | Rf pick-up coil circuit for a wide tuning range nuclear magnetic resonance probe |
US4365303A (en) * | 1980-02-07 | 1982-12-21 | The Perkin-Elmer Corporation | Method and apparatus for determining the nature of an unknown chemical substance |
US4682027A (en) * | 1986-04-25 | 1987-07-21 | Varian Associates, Inc. | Method and apparatus for sample confirmation in gas chromatography |
US4692685A (en) * | 1984-03-14 | 1987-09-08 | Blaze Kevin L | Electrical measuring apparatus, and methods for determining the condition or identity of biological material |
US4751515A (en) * | 1980-07-09 | 1988-06-14 | Corum James F | Electromagnetic structure and method |
US4822169A (en) * | 1984-09-25 | 1989-04-18 | Richard Distl | Measuring assembly for analyzing electromagnetic radiation |
US5254950A (en) * | 1991-09-13 | 1993-10-19 | The Regents, University Of California | DC superconducting quantum interference device usable in nuclear quadrupole resonance and zero field nuclear magnetic spectrometers |
US5339811A (en) * | 1991-11-06 | 1994-08-23 | Mitsui Mining & Smelting Co., Ltd. | Magnetoencephalograph |
US5343147A (en) * | 1992-09-08 | 1994-08-30 | Quantum Magnetics, Inc. | Method and apparatus for using stochastic excitation and a superconducting quantum interference device (SAUID) to perform wideband frequency response measurements |
US5446681A (en) * | 1990-10-12 | 1995-08-29 | Exxon Research And Engineering Company | Method of estimating property and/or composition data of a test sample |
US5458142A (en) * | 1993-03-19 | 1995-10-17 | Farmer; Edward J. | Device for monitoring a magnetic field emanating from an organism |
US5465049A (en) * | 1992-04-20 | 1995-11-07 | Sumitomo Electric Industries, Ltd. | Integrated type planar magnetic sensor having SQUID and flux transformer formed of oxide superconductor |
US5508203A (en) * | 1993-08-06 | 1996-04-16 | Fuller; Milton E. | Apparatus and method for radio frequency spectroscopy using spectral analysis |
US5541413A (en) * | 1992-04-24 | 1996-07-30 | Thiokol Corporation | Acousto-optic tunable filter-based surface scanning system and process |
US5574369A (en) * | 1994-08-19 | 1996-11-12 | Hibbs; Andrew D. | Detection and communications device employing stochastic resonance |
US5583432A (en) * | 1994-04-11 | 1996-12-10 | Sci-Nostics Limited | Electrical method and apparatus for non-contact determination of physical and/or chemical properties of a sample, particularly of blood |
US5656937A (en) * | 1995-06-07 | 1997-08-12 | Conductus, Inc. | Low-noise symmetric dc SQUID system having two pairs of washer coils and a pair of Josephson junctions connected in series |
US5696691A (en) * | 1995-06-29 | 1997-12-09 | The United States Of America As Represented By The Secretary Of The Navy | Self-adjusting statistical noise analyzer with interference suppression |
US5734353A (en) * | 1995-08-14 | 1998-03-31 | Vortekx P.C. | Contrawound toroidal helical antenna |
US5752514A (en) * | 1995-08-31 | 1998-05-19 | Shimadzu Corporation | Biomagnetism measuring method and apparatus |
US5789961A (en) * | 1996-06-28 | 1998-08-04 | The United States Of America As Represented By The Secretary Of The Navy | Noise- and coupling-tuned signal processor with arrays of nonlinear dynamic elements |
US5944782A (en) * | 1996-10-16 | 1999-08-31 | Veritas Software Corporation | Event management system for distributed computing environment |
US5955400A (en) * | 1996-08-02 | 1999-09-21 | Hitachi, Ltd. | SQUID integrated with pickup coils |
US5959548A (en) * | 1997-10-31 | 1999-09-28 | Halliburton Energy Services, Inc. | Electromagnetic signal pickup device |
US6020782A (en) * | 1994-05-25 | 2000-02-01 | The United States Of America As Represented By The Secretary Of The Navy | Noise assisted signal processor with nonlinearly coupled arrays of nonlinear dynamic elements |
US6028558A (en) * | 1992-12-15 | 2000-02-22 | Van Voorhies; Kurt L. | Toroidal antenna |
US6084399A (en) * | 1997-03-26 | 2000-07-04 | Sumitomo Electric Industries, Ltd. | Method and apparatus for determining concentration of magnetic substances in a non-magnetic substance using a SQUID |
US6136541A (en) * | 1999-02-22 | 2000-10-24 | Vialogy Corporation | Method and apparatus for analyzing hybridized biochip patterns using resonance interactions employing quantum expressor functions |
US6142681A (en) * | 1999-02-22 | 2000-11-07 | Vialogy Corporation | Method and apparatus for interpreting hybridized bioelectronic DNA microarray patterns using self-scaling convergent reverberant dynamics |
US6150812A (en) * | 1994-02-21 | 2000-11-21 | Fermion Inc. | Detection of electromagnetic fields as a determinant of an event |
US6159444A (en) * | 1998-09-11 | 2000-12-12 | The Regents Of The University Of California | NMR/MRI with hyperpolarized gas and high Tc SQUID |
US6196057B1 (en) * | 1998-04-02 | 2001-03-06 | Reliance Electric Technologies, Llc | Integrated multi-element lubrication sensor and lubricant health assessment |
US6285249B1 (en) * | 2000-01-21 | 2001-09-04 | The United States Of America As Represented By The Secretary Of The Navy | Controlled stochastic resonance circuit |
US6294911B1 (en) * | 1998-03-12 | 2001-09-25 | Tdk Corporation | Measurement method of magnetization direction of magnetoresistive effect devices and measurement apparatus based on the method of TDK corporation |
US6320369B1 (en) * | 1998-10-07 | 2001-11-20 | Nec Corporation | Superconducting current measuring circuit having detection loop |
US6323632B1 (en) * | 1999-08-13 | 2001-11-27 | Coulter International Corp. | Solid state RF oscillator-detector for flow cytometer |
US20020158631A1 (en) * | 2001-04-27 | 2002-10-31 | Akihiko Kandori | Biomagnetic field measuring apparatus |
US6541978B1 (en) * | 1998-09-23 | 2003-04-01 | Digibio | Method, system and device for producing signals from a substance biological and/or chemical activity |
US6724188B2 (en) * | 2002-03-29 | 2004-04-20 | Wavbank, Inc. | Apparatus and method for measuring molecular electromagnetic signals with a squid device and stochastic resonance to measure low-threshold signals |
US6760674B2 (en) * | 2001-10-08 | 2004-07-06 | Microchip Technology Incorporated | Audio spectrum analyzer implemented with a minimum number of multiply operations |
US6815949B2 (en) * | 2001-07-19 | 2004-11-09 | Hitachi, Ltd. | Apparatus for measuring a magnetic field |
US6885192B2 (en) * | 2002-02-06 | 2005-04-26 | The Regents Of The University Of California | SQUID detected NMR and MRI at ultralow fields |
US6952652B2 (en) * | 2002-04-19 | 2005-10-04 | Wavbank, Inc. | System and method for sample detection based on low-frequency spectral components |
US6995558B2 (en) * | 2002-03-29 | 2006-02-07 | Wavbank, Inc. | System and method for characterizing a sample by low-frequency spectra |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2598414B1 (en) | 1985-11-12 | 1988-10-14 | Centre Nat Rech Scient | PHARMACEUTICAL COMPOSITION CONTAINING HISPIDULIN OR A DERIVATIVE AND USE OF SUCH COMPOUNDS IN THE PREPARATION OF ANTI-ASTHMATIC COMPOSITIONS |
JP2624696B2 (en) * | 1987-07-31 | 1997-06-25 | 日本分光株式会社 | Spectrum estimator |
CA2008009C (en) * | 1989-01-20 | 1994-05-03 | Hajime Hayashi | Apparatus for measuring magnetic field |
FR2659740B1 (en) | 1990-03-13 | 1992-07-10 | Inst Nat Sante Rech Med | METHOD FOR CONTROLLING THE DILUTION AND CONTAMINATION OF HIGH DILUTION SOLUTIONS. |
FR2659553B1 (en) | 1990-03-13 | 1992-07-10 | Inst Nat Sante Rech Med | PROCESS FOR THE PREPARATION OF HOMEOPATHIC COMPOSITIONS FROM AN INITIAL SOLUTION CONTAINING AN ACTIVE SUBSTANCE. |
JPH0464059A (en) * | 1990-07-03 | 1992-02-28 | Kuwabara Yasunaga | Processing device of analysis data |
US5311125A (en) * | 1992-03-18 | 1994-05-10 | Lake Shore Cryotronics, Inc. | Magnetic property characterization system employing a single sensing coil arrangement to measure AC susceptibility and DC moment of a sample |
US6084309A (en) * | 1992-10-20 | 2000-07-04 | Fujitsu Limited | Semiconductor device and semiconductor device mounting structure |
FR2700628B1 (en) | 1993-01-21 | 1995-03-24 | Benvenistre Jacques | Method and device for transmitting in the form of a signal the biological activity of a carrier material to another carrier material, and for processing such a signal, and product obtained with such a method. |
US5606691A (en) * | 1994-11-10 | 1997-02-25 | Harman Interactive Inc. | Method of searching a database using selected criterion having implicit logical operation |
US5737342A (en) * | 1996-05-31 | 1998-04-07 | Quantum Corporation | Method for in-chip testing of digital circuits of a synchronously sampled data detection channel |
US6133734A (en) * | 1997-12-30 | 2000-10-17 | Schlumberger Technology Corporation | Method and apparatus for evaluating an earth formation using nuclear magnetic resonance techiques |
FR2777656B1 (en) | 1998-04-20 | 2000-06-16 | Digibio | METHOD FOR AMPLIFYING THE FORMATION OF LIGAND-RECEPTOR COMPLEXES AND USES THEREOF |
FR2780651B1 (en) | 1998-07-01 | 2001-07-20 | Digibio | METHOD FOR ACTIVATING AN INACTIVE AND VERY LOW CONCENTRATION SOLUTION OF A BIOLOGICAL AND / OR CHEMICAL SUBSDETERMINE IN A SOLVENT |
FR2783606B1 (en) | 1998-09-23 | 2000-11-24 | Digibio | METHOD AND SYSTEM FOR PRODUCING A SUBSTANCE OR SIGNAL HAVING A COAGULATING OR ANTICOAGULATING EFFECT. THERAPEUTIC APPLICATIONS OF SAID SUBSTANCE OR SIGNAL |
EP1137951A1 (en) * | 1998-12-10 | 2001-10-04 | University Of Maryland | Scanning single electron transistor microscope for imaging ambient temperature objects |
US6411108B1 (en) * | 1999-11-05 | 2002-06-25 | Sensor Technologies, Inc. | Noncontact signal analyzer |
JP2001237901A (en) * | 2000-02-23 | 2001-08-31 | Ando Electric Co Ltd | Waveform analysis method and device therefor |
US6201821B1 (en) | 2000-06-05 | 2001-03-13 | Agilent Technologies, Inc. | Coherent population trapping-based frequency standard having a reduced magnitude of total a.c. stark shift |
JP3877131B2 (en) * | 2000-11-10 | 2007-02-07 | ホーチキ株式会社 | Fire detector and fire detection method |
US6586931B2 (en) * | 2001-04-20 | 2003-07-01 | Baker Hughes Incorporated | NMR logging in the earth's magnetic field |
US6825655B2 (en) * | 2001-08-31 | 2004-11-30 | Imego Ab | Method and arrangement for detecting changes of a magnetic response in magnetic particles |
-
2003
- 2003-10-09 US US10/683,875 patent/US6995558B2/en not_active Expired - Lifetime
-
2004
- 2004-10-08 JP JP2006534425A patent/JP4425922B2/en active Active
- 2004-10-08 AU AU2004280998A patent/AU2004280998B2/en active Active
- 2004-10-08 CA CA2538988A patent/CA2538988C/en active Active
- 2004-10-08 EP EP04785407A patent/EP1671144A4/en not_active Withdrawn
- 2004-10-08 CN CN2004800294902A patent/CN1864073B/en active Active
- 2004-10-08 BR BRPI0415235-2A patent/BRPI0415235B1/en not_active IP Right Cessation
- 2004-10-08 US US10/575,182 patent/US20070210790A1/en not_active Abandoned
- 2004-10-08 WO PCT/US2004/033383 patent/WO2005036131A2/en active Application Filing
-
2006
- 2006-02-07 US US11/349,857 patent/US20060158183A1/en not_active Abandoned
-
2009
- 2009-11-12 JP JP2009258652A patent/JP2010032551A/en active Pending
Patent Citations (47)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4031462A (en) * | 1975-07-07 | 1977-06-21 | Motorola, Inc. | Frequency spectrum analyzer |
US4095168A (en) * | 1977-02-22 | 1978-06-13 | Varian Associates, Inc. | Rf pick-up coil circuit for a wide tuning range nuclear magnetic resonance probe |
US4365303A (en) * | 1980-02-07 | 1982-12-21 | The Perkin-Elmer Corporation | Method and apparatus for determining the nature of an unknown chemical substance |
US4751515A (en) * | 1980-07-09 | 1988-06-14 | Corum James F | Electromagnetic structure and method |
US4692685A (en) * | 1984-03-14 | 1987-09-08 | Blaze Kevin L | Electrical measuring apparatus, and methods for determining the condition or identity of biological material |
US4822169A (en) * | 1984-09-25 | 1989-04-18 | Richard Distl | Measuring assembly for analyzing electromagnetic radiation |
US4682027A (en) * | 1986-04-25 | 1987-07-21 | Varian Associates, Inc. | Method and apparatus for sample confirmation in gas chromatography |
US5446681A (en) * | 1990-10-12 | 1995-08-29 | Exxon Research And Engineering Company | Method of estimating property and/or composition data of a test sample |
US5254950A (en) * | 1991-09-13 | 1993-10-19 | The Regents, University Of California | DC superconducting quantum interference device usable in nuclear quadrupole resonance and zero field nuclear magnetic spectrometers |
US5339811A (en) * | 1991-11-06 | 1994-08-23 | Mitsui Mining & Smelting Co., Ltd. | Magnetoencephalograph |
US5465049A (en) * | 1992-04-20 | 1995-11-07 | Sumitomo Electric Industries, Ltd. | Integrated type planar magnetic sensor having SQUID and flux transformer formed of oxide superconductor |
US5541413A (en) * | 1992-04-24 | 1996-07-30 | Thiokol Corporation | Acousto-optic tunable filter-based surface scanning system and process |
US5343147A (en) * | 1992-09-08 | 1994-08-30 | Quantum Magnetics, Inc. | Method and apparatus for using stochastic excitation and a superconducting quantum interference device (SAUID) to perform wideband frequency response measurements |
US6204821B1 (en) * | 1992-12-15 | 2001-03-20 | West Virginia University | Toroidal antenna |
US6028558A (en) * | 1992-12-15 | 2000-02-22 | Van Voorhies; Kurt L. | Toroidal antenna |
US5458142A (en) * | 1993-03-19 | 1995-10-17 | Farmer; Edward J. | Device for monitoring a magnetic field emanating from an organism |
US5508203A (en) * | 1993-08-06 | 1996-04-16 | Fuller; Milton E. | Apparatus and method for radio frequency spectroscopy using spectral analysis |
US6150812A (en) * | 1994-02-21 | 2000-11-21 | Fermion Inc. | Detection of electromagnetic fields as a determinant of an event |
US5583432A (en) * | 1994-04-11 | 1996-12-10 | Sci-Nostics Limited | Electrical method and apparatus for non-contact determination of physical and/or chemical properties of a sample, particularly of blood |
US6020782A (en) * | 1994-05-25 | 2000-02-01 | The United States Of America As Represented By The Secretary Of The Navy | Noise assisted signal processor with nonlinearly coupled arrays of nonlinear dynamic elements |
US5574369A (en) * | 1994-08-19 | 1996-11-12 | Hibbs; Andrew D. | Detection and communications device employing stochastic resonance |
US5656937A (en) * | 1995-06-07 | 1997-08-12 | Conductus, Inc. | Low-noise symmetric dc SQUID system having two pairs of washer coils and a pair of Josephson junctions connected in series |
US5696691A (en) * | 1995-06-29 | 1997-12-09 | The United States Of America As Represented By The Secretary Of The Navy | Self-adjusting statistical noise analyzer with interference suppression |
US5734353A (en) * | 1995-08-14 | 1998-03-31 | Vortekx P.C. | Contrawound toroidal helical antenna |
US5952978A (en) * | 1995-08-14 | 1999-09-14 | Vortekx, Inc. | Contrawound toroidal antenna |
US5752514A (en) * | 1995-08-31 | 1998-05-19 | Shimadzu Corporation | Biomagnetism measuring method and apparatus |
US5789961A (en) * | 1996-06-28 | 1998-08-04 | The United States Of America As Represented By The Secretary Of The Navy | Noise- and coupling-tuned signal processor with arrays of nonlinear dynamic elements |
US5955400A (en) * | 1996-08-02 | 1999-09-21 | Hitachi, Ltd. | SQUID integrated with pickup coils |
US5944782A (en) * | 1996-10-16 | 1999-08-31 | Veritas Software Corporation | Event management system for distributed computing environment |
US6084399A (en) * | 1997-03-26 | 2000-07-04 | Sumitomo Electric Industries, Ltd. | Method and apparatus for determining concentration of magnetic substances in a non-magnetic substance using a SQUID |
US5959548A (en) * | 1997-10-31 | 1999-09-28 | Halliburton Energy Services, Inc. | Electromagnetic signal pickup device |
US6294911B1 (en) * | 1998-03-12 | 2001-09-25 | Tdk Corporation | Measurement method of magnetization direction of magnetoresistive effect devices and measurement apparatus based on the method of TDK corporation |
US6196057B1 (en) * | 1998-04-02 | 2001-03-06 | Reliance Electric Technologies, Llc | Integrated multi-element lubrication sensor and lubricant health assessment |
US6159444A (en) * | 1998-09-11 | 2000-12-12 | The Regents Of The University Of California | NMR/MRI with hyperpolarized gas and high Tc SQUID |
US6541978B1 (en) * | 1998-09-23 | 2003-04-01 | Digibio | Method, system and device for producing signals from a substance biological and/or chemical activity |
US6320369B1 (en) * | 1998-10-07 | 2001-11-20 | Nec Corporation | Superconducting current measuring circuit having detection loop |
US6136541A (en) * | 1999-02-22 | 2000-10-24 | Vialogy Corporation | Method and apparatus for analyzing hybridized biochip patterns using resonance interactions employing quantum expressor functions |
US6142681A (en) * | 1999-02-22 | 2000-11-07 | Vialogy Corporation | Method and apparatus for interpreting hybridized bioelectronic DNA microarray patterns using self-scaling convergent reverberant dynamics |
US6323632B1 (en) * | 1999-08-13 | 2001-11-27 | Coulter International Corp. | Solid state RF oscillator-detector for flow cytometer |
US6285249B1 (en) * | 2000-01-21 | 2001-09-04 | The United States Of America As Represented By The Secretary Of The Navy | Controlled stochastic resonance circuit |
US20020158631A1 (en) * | 2001-04-27 | 2002-10-31 | Akihiko Kandori | Biomagnetic field measuring apparatus |
US6815949B2 (en) * | 2001-07-19 | 2004-11-09 | Hitachi, Ltd. | Apparatus for measuring a magnetic field |
US6760674B2 (en) * | 2001-10-08 | 2004-07-06 | Microchip Technology Incorporated | Audio spectrum analyzer implemented with a minimum number of multiply operations |
US6885192B2 (en) * | 2002-02-06 | 2005-04-26 | The Regents Of The University Of California | SQUID detected NMR and MRI at ultralow fields |
US6724188B2 (en) * | 2002-03-29 | 2004-04-20 | Wavbank, Inc. | Apparatus and method for measuring molecular electromagnetic signals with a squid device and stochastic resonance to measure low-threshold signals |
US6995558B2 (en) * | 2002-03-29 | 2006-02-07 | Wavbank, Inc. | System and method for characterizing a sample by low-frequency spectra |
US6952652B2 (en) * | 2002-04-19 | 2005-10-04 | Wavbank, Inc. | System and method for sample detection based on low-frequency spectral components |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9417257B2 (en) | 2004-07-27 | 2016-08-16 | Nativis, Inc. | System and method for collecting, storing, processing, transmitting and presenting very low amplitude signals |
US20100277164A1 (en) * | 2006-09-01 | 2010-11-04 | Commonwealth Scientific And Industrial Research Organisation | Method and apparatus for signal recovery |
WO2014145114A1 (en) * | 2013-03-15 | 2014-09-18 | Aliphcom | Dynamic control of sampling rate of motion to modify power consumption |
US10046172B2 (en) | 2013-03-15 | 2018-08-14 | Nativis, Inc. | Controller and flexible coils for administering therapy, such as for cancer therapy |
US11103721B2 (en) | 2013-03-15 | 2021-08-31 | Natives, Inc. | Controller and flexible coils for administering therapy, such as for cancer therapy |
US11097623B2 (en) | 2014-09-24 | 2021-08-24 | Rohm Co., Ltd. | Current mode control type switching power supply device |
CN109686528A (en) * | 2018-12-18 | 2019-04-26 | 武汉船用电力推进装置研究所(中国船舶重工集团公司第七一二研究所) | A kind of HTS-SMES magnet device |
Also Published As
Publication number | Publication date |
---|---|
AU2004280998B2 (en) | 2008-04-10 |
CA2538988A1 (en) | 2005-04-21 |
WO2005036131A2 (en) | 2005-04-21 |
US20070210790A1 (en) | 2007-09-13 |
AU2004280998A1 (en) | 2005-04-21 |
BRPI0415235B1 (en) | 2017-07-04 |
CN1864073A (en) | 2006-11-15 |
JP4425922B2 (en) | 2010-03-03 |
CA2538988C (en) | 2011-02-15 |
JP2010032551A (en) | 2010-02-12 |
EP1671144A4 (en) | 2007-04-18 |
WO2005036131A3 (en) | 2005-07-07 |
CN1864073B (en) | 2010-05-05 |
US20040183530A1 (en) | 2004-09-23 |
US6995558B2 (en) | 2006-02-07 |
EP1671144A2 (en) | 2006-06-21 |
JP2007533972A (en) | 2007-11-22 |
BRPI0415235A (en) | 2006-12-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6995558B2 (en) | System and method for characterizing a sample by low-frequency spectra | |
CA2473142C (en) | System and method for characterizing a sample by low-frequency spectra | |
AU2013290020B2 (en) | Miniaturized molecular interrogation and data system | |
US6952652B2 (en) | System and method for sample detection based on low-frequency spectral components | |
CA2573350C (en) | System and method for collecting, storing, processing, transmitting and presenting very low amplitude signals | |
AU2007321907A1 (en) | Apparatus and method for transducing an in vitro or mammalian system with a low-frequency signal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: NATIVIS, INC., CALIFORNIA Free format text: CHANGE OF NAME;ASSIGNOR:WAVBANK, INC.;REEL/FRAME:017766/0105 Effective date: 20060213 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |