Laser and environmental monitoring system

Description of equivalent WO 2007001308 A2

LASER AND ENVIRONMENTAL MONITORING SYSTEM

CROSS REFERENCE TO RELATED APPLICATIONS [0001] This is a continuation of 10/884,798, filed July 2, 2004, which is a continuation-in-part of 10/791 ,377, filed March 2, 2004, which is a continuation- in-part of 10/265,211 , filed October 4, 2002, which is a continuation-in-part of PCT/US02/02548, filed January 28, 2002, which claims priority to provisional application 60/265,133, filed January 30, 2001. Furthermore, this is continuation- in-part of 10/265,211 , filed October 4, 2002, which is a continuation-in-part of PCT/US02/02548, filed January 28, 2002, which claims priority to provisional application 60/265,133, filed January 30, 2001. This is also a continuation-in- part of 10/628,874, filed July 28, 2003, which is a continuation of PCT/US02/02548, filed January 28, 2002, which claims priority to provisional application 60/265,133, filed January 30, 2001. All of these applications are incorporated by reference herein.

BACKGROUND [0002] The present invention generally relates to a monitoring system and more particularly to a laser and environmental monitoring system.

[0003] Conventionally, laser desorption mass spectrometry has been used with a fixed laser beam pulse shape and computers for simple chemical analysis processes on purified molecules. The laser beam pulse shape was not considered an important parameter and was not modified; whatever fixed shape was set by the manufacturer for the ultraviolet laser was used in the tests. The general concept of typically laser selective ion formation from molecules in a molecular beam is disclosed in the following publication: Assion et al., "Control of Chemical Reactions by Feedback-Optimized Phase-Shaped Femtosecond Laser Pulses," Science, Vol. 282, page 919 (October 30, 1998). The pulse shaping process with a learning algorithm is disclosed in Judson et al., "Teaching Lasers to Control Molecules," Physical Review Letters, Vol. 68, No. 10, page 1500 (March 9, 1992). It is noteworthy, however, that the Assion article discloses use of an 80 femtosecond laser pulse and requires molecules to be isolated in a molecular beam, while the Judson article discloses use of a one nanosecond laser pulse and is purely conceptual as it does not include experimental results.

[0004] Commercially practical femtosecond lasers have been unavailable until recently. For example, lasers which can generate 10 femtosecond or less laser pulse durations have traditionally been extremely expensive, required unrealistically high electrical energy consumption (for extensive cooling, by way of example) and depended on laser dyes that had to be replenished every month thereby leading to commercial impracticality.

[0005] Ultrashort pulses are prone to suffer phase distortions as they propagate through or reflect from optics because of their broad bandwidth. There has been significant progress in correcting these unwanted phase distortions. There have been recent experimental attempts to purposely shape the phase of ultrashort pulses since shaped pulses have been shown to increase the yield of certain chemical reactions and multiphoton excitation, although the mechanism for the observed changes remains unknown in most cases. As usually practiced, the output waveform is determined by the Fourier transform (FT) of a spatial pattern transferred by a mask or a modulator array onto the dispersed optical spectrum. The introduction of liquid crystal modulator arrays and acousto-optic (A/O) modulators into FT pulse shapers led to computer programmable pulse shaping, with millisecond and microsecond reprogramming times, respectively, and widespread adoption of this technique. These shaped pulses require a very large data set and in many cases, complex learning calculations for determining the pulse shaping characteristics for a particular application. The optimal pulse for the particular application is not known in advance. Since the variation shape of the possible pulse shapes is huge, scanning the entire parameter space is impossible and as such the optimized pulse shape could not have been predicted by theory. For a pulse shaper with Λ/ pixels, one can generate (P *A)^N shaped pulses, where P and A are the number of different phases and amplitudes a pixel can take. If it is assumed 100 pixels, each taking 10 different amplitude values and 100 different phase values, the number of different pulses is of order of magnitude 10³⁰⁰. This dataset is extremely large, therefore, while in principle, the field exists to achieve the desired photonic transformation or excitation, finding it is a great challenge. It would be desirable for a system to control ultrashort pulses with a smaller dataset, operable to generate very complex pulse shapes that are optimal for the particular application and are highly reproducible.

[0006] Additionally, monitoring the environment for chemical and biological agents, including explosives, from terrorist threats or from industrial contamination has become a necessity for reasons of national security and the well being of humans. Conventional devices are only designed for use to detect a single known agent or are inaccurate. Accordingly, to avoid a costly false positive or false negative identification, it would be desirable to employ femtosecond laser and control technology to environmental monitoring.

SUMMARY OF THE INVENTION [0007] In accordance with the present invention, a laser and monitoring system is provided. In another aspect of the present invention, the system includes a laser, pulse shaper and detection device. A further aspect of the present invention employs a femtosecond laser and binary pulse shaping (BPS). Still another aspect of the present invention uses a laser beam pulse, a pulse shaper and a SHG crystal. In yet another aspect of the present invention, the multiphoton intrapulse interference phase scan (hereinafter "MIIPS") method is used to characterize the spectral phase of femtosecond laser pulses and to correct them. A further aspect of the system of the present invention is employed to monitor environmental chemical and biological agents, including toxins, explosives, and diseases, among others.

[0008] The laser and environmental monitoring system of the present invention is advantageous over conventional constructions since the present invention allows for automated analysis and identification of molecules in complex mixtures. It is envisioned that the present invention is capable of fast (for example, one second), accurate (for example, even in a chemically complex environment), robust (for example, stand alone, closed-loop and portable) and reproducible sensing. Operationally, the system employs a computer controlled pulse shaping module which interfaces with a commercially available femtosecond pulsed laser and mass spectrometry module thereby using cost effective instruments. The system operates on the principle of molecular control based on shaped laser fields and using genetic and evolutionary, learning search methods, whereby a series of laser fields are chosen to unequivocally identify each chemical or biological agent of interest based on their electronic and nuclear structure. In case of a positive identification, the unit automatically contacts a command center and/or changes the state of a building heating, ventilating and air conditioning system in order to advantageously minimize spreading of the agent or to neutralize the agent, in a very fast and real-time analysis and feedback manner. Additional advantages and features of the present invention will become apparent from the following description and appended claims, taken in conjunction with the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

[0009] Figure 1 is a diagrammatic perspective view showing a first preferred embodiment building employing the laser and environmental monitoring system of the present invention;

[0010] Figure 2 is a diagrammatic side view showing a second preferred embodiment building employing the system;

[0011] Figures 3 and 4 are diagrammatic views showing the preferred embodiment system;

[0012] Figure 5 is a diagrammatic view showing the preferred embodiment system employing MIIPS;

[0013] Figures 6A-C are flow charts showing the operation of the preferred embodiment system in the evolutionary learning, library mode;

[0014] Figures 7, 8A and 8B are flow charts for the method and computer software operation employed with the preferred embodiment system;

[0015] Figure 9 is a schematic representation of an SLM-256 (CRI, Inc.) amplitude/phase mask employed with one variation of the preferred embodiment system; [0016] Figure 10 shows a diagrammatic representation of the direction of electric field rotation which is different for the two liquid crystal masks employed with the preferred embodiment system of Figure 9;

[0017] Figure 11 shows an experimentally measured transmission as a function of voltage used for calibration of the SLM mask of Figure 9 proposed to be employed with the present invention;

[0018] Figure 12 shows experimental demonstration of phase characterization followed by compensation proposed to be employed with the present invention system (the transform limited pulse width retrieved from FROG was 28-fs, the MIIPS traces are shown on the left column and the SHG-FROG traces are shown on the right column);

[0019] Figure 13 shows schematic representations of three binary phase masks that have been optimized by a genetic learning program to maximize the indicated intensity ratio of the fragment ions labeled a, b and c, as proposed to be employed with the present invention system;

[0020] Figure 14 is a flow chart showing the operation of the preferred embodiment system used in the field to monitor the environment;

[0021] Figure 15 is a perspective view showing a first alternate embodiment of the present invention system; and [0022] Figure 16 is a perspective view showing a second alternate embodiment of the present invention system;

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS [0023] A laser and environmental monitoring system 1 is provided to monitor airborne chemical and biological agents in the desired environmental specimen or target area. This unit, because of its sensitivity, can be used to sense smoke from a fire. Referring to Figure 1 , a first preferred embodiment system 1 of the present invention includes a multi-story building 2, having a plurality of wings with interconnecting hallways, with each wing having separate ventilating systems, such as heating/ventilating/air conditioning climate control units (hereinafter "HVAC" units) 3. A laser and detection apparatus 21 is housed in a self-contained unit 4 with inlet valve controls 5. Furthermore, an HVAC control unit 6, a fire door/alarm/sprinkler emergency control unit 7 and a remote command center 8, such as a public safety or fire department 8, are also provided. Laser and detection apparatus 21 includes a femtosecond laser, pulse shaping optics, a mass spectrometer and a computer controller 37 as will be discussed in greater detail hereinafter. Computer controller 37 is electrically connected to HVAC controller 6, emergency controller 7, and command center 8 to allow for automatic two way communications. Computer controller 37 of unit 4 further automatically controls valve settings of a manifold or the like coupled to ducts 9, flexible hoses or other conduits which are coupled to an exhaust plenum in each discrete room, floor and/or wing subject to environmental monitoring.

[0024] Figure 2 discloses a second preferred embodiment building employing the system 1 of the present invention. In this embodiment, air analysis and monitoring unit 4 is directly coupled within the exhaust conduits for each of the three floors of the building shown. This embodiment is employed separately with each building wing or separate building structure, and multiple communication lines can optionally couple the computer controllers of each unit 4 with a main command station computer or the like.

[0025] Referring to Figures 1 , 2 and 14, for both preferred embodiments, monitoring and air analysis unit 4 is initially set up at its manufacturing facility with a library of pre-stored data of acceptable background data, unacceptable chemical and biological agent data, and the corresponding pulse shaping control data. The self-contained unit 4 is subsequently installed at the facility or building 2. Unit 4 is acting as an "artificial nose" in that it receives exhausted air from inside the target environmental area (depending upon the HVAC ducting and zone configuration), as selectively regulated by the valve controls 5 in coordination with the laser and detection computer controller 37. Upon installation and set up, and optionally at regular intervals of known-safe conditions, the laser and detection computer controller 37 senses or monitors acceptable background chemicals in the air, such as oxygen, carbon dioxide, chlorine (if adjacent a swimming pool), benzene (if adjacent a bus station or high volume street), or hexachlorabenzene (if the building is a hospital), or the like. These acceptable background chemicals are automatically stored in the memory of computer controller 37 and matched with pulse shapes and mass spectrums previously stored in the library of data. Active baseline subtraction is thereafter employed during regular environmental monitoring so these benign chemicals and their fragments are not mistakenly considered part of a dangerous compound, thereby reducing the risk of creating a false detection and warning.

[0026] Subsequently, the monitoring and analysis unit 4 automatically samples air from each environmental target area at a repetition rate of approximately one per minute, a rate which is then increased to approximately one thousand times per second if an unusual background chemical or undesirable agent is detected. The background and regularly monitored chemical analyses are stored in the computer memory for later automatic and/or manual comparisons in the event of agent detection. The varying sample rate minimizes the memory and historical trend information from becoming overwhelmed with unnecessary results, and the variable repetition rate saves laser power costs and prolongs the life of laser components and detection equipment.

[0027] If an undesirably harmful chemical or biological agent is sensed, then the laser and detection computer controller will automatically override and dictate action of the HVAC controller 6 and/or emergency controller 7, while transmitting warning and other calculated information to the command center 8. Examples of such undesirable agents include and are not limited to: methane, propane, mustard gas, nerve agents, phosgene, vinyl chloride, cyanide, anthrax and other pathogens and diseases. Depending upon the agent sensed and calculated, the computer controller will initiate low impact counter-measures such as causing closure of the HVAC dampers or exhaust doors in order to isolate the infected room and environment while further analysis is ongoing at a very high rate. Monitoring unit 4 will also recalibrate itself and use multiple shaped pulses to reduce the likelihood of error, and to insure that there is no mass spec overlap of results between undesired agents and acceptable chemicals. The unit will further calculate a confidence percentage of the test results compared to an unacceptable agent believed to have been identified based on previously stored library data. If a low confidence level is determined then the unit will rapidly conduct one or more further tests on the environment within seconds of the initial noteworthy reading.

[0028] If the confidence level calculated of an undesirable agent is statistically significant then laser and detection computer controller 37 will enact serious counter-measures such as turning off and isolating all HVAC units, causing closure of fire doors, and/or sending alarms to evacuate humans from the building. The laser and detection computer will also send the appropriate warnings and information to the command center and department 8 including currently calculated values and historical trends. Furthermore, monitoring unit 4 will continue sensing all environmental areas within the building so as to track movement and cross-contamination between environmental areas for use by emergency personnel. As an additional option, the computer controller of monitoring unit 4 will enact neutralizing actions depending upon the agent determination, such as opening external building exhaust vents and/or windows in order to flush out the contaminant agent to the external atmosphere, turn on the water sprinklers within the building, and/or inject direct countermeasures such as the introduction of an antidote chemical into the inlet side of the HVAC system so as to neutralize the agent.

[0029] On a more detailed level within monitoring unit 4, operational parameters such as ultimate sensitivity are initially obtained: once six predetermined shaped laser pulses (as more fully described hereinafter for the library) are found for a given molecule, together with the TL pulse (also as more fully described hereinafter for the library), the system is ready to make a positive identification of that particular chemical. In fact, every time the measurement is performed with a different pulse shape, the accuracy is improved exponentially with the number of independent measurements. This multidimensional analysis of the present invention is believed to make it a million times more accurate than presently used mass spectrometry units. The goal is to build a library of multidimensional spectra that result from shaped laser pulses that have been optimized for maximum contrast for each given chemical or biological threat agent (and their analogs). It is important to optimize the pulses for each molecule, because each molecule has a distinct electronic and nuclear structure that determines its susceptibility for ionization and fragmentation. It is noteworthy again that the evolutionary learning program searches will only be used in a laboratory setting and not in the field with the preferred embodiment system.

[0030] The present invention pulse shaper monitoring of the environment will be carried out as follows. Every second, the unit will obtain a mass spectrum using TL pulses. Under these conditions, the system is at its highest sensitivity. If a mass spectrum indicates that a possible chemical agent is present, the system immediately narrows the list of possibilities based on the available TL mass spectrum and performs tests with the first three fields defined in the library for the suspected chemical agents. These subsequent tests are run at the maximum repetition rate. Within seconds the unit should have confirmation of the suspected chemical agent(s) present. At this point the unit runs a calibration test to make sure that the laser, shaper and mass spectrometer are operating within specifications. The unit performs three additional measurements with the final three shaped laser pulses for a final and definitive determination. If the identification is positive for a known threat, the unit contacts a command center and uplinks the experimental data used to make the determination. The command center has the required information to make the best strategic decision based on solid information within a minute of detection. Additional rounds of tests can be requested from the command center, with each full cycle taking less than one minute to complete. The speed and efficiency of the method is achieved by having previously determined the six best shaped laser pulses for each molecule. The pulse shapes in the library are determined using the automated evolutionary learning program in a laboratory setting. Updated pulse shape libraries are periodically uploaded to the computer controller of the monitoring device.

[0031 ] One of the advantages of the system and method of the present invention is that it is robust enough to operate in the presence of a complex chemical environment. For example, if the building environment being monitored is next to many automotive vehicles or a bus stop, the unit will detect diesel, gasoline and exhaust fumes. Under these conditions, it will be imperative to take advantage of the multidimensional properties of the present invention. TL pulses will yield mass spectra with a number of suspect ion masses. However, upon further analysis, the unit will determine that no unacceptable, dangerous chemicals are present. If these conditions persist, the computer will automatically define a background level of contaminants and actively subtract it. The power of multidimensional analysis is extremely valuable if a mixture of chemical agents is present. The present invention system would have no problems determining which chemicals are present even under these adverse conditions. Finally, if a new chemical threat is developed by a terrorist group, the command center can update all the units in the field with a new library that contains the spectrometric signature that is obtained by multidimensional laser interrogation of that compound. The addition to the library would take less than a day to generate in the laboratory, and would be available for immediate use by the monitoring unit upon upload.

[0032] One laser and detector unit creates the library. That unit is in a laboratory and is capable of handling nasty chemicals. The other unit, is compact and automated and uses the libraries prepared by the first unit. The second unit can be compact and is field transportable. Both units have fs laser, shaper, miips, binary shaping, computer controller. The only difference is that one is used to generate the library in a laboratory setting. The other one is fully automated and field deployable. It is made to be a closed box that requires minimum maintenance.

Pulse Shaping System

[0033] The preferred embodiment of a control system and apparatus 21 of the present invention for use with the laser and environmental monitoring system is generally shown in Figures 3 and 4. Apparatus 21 includes a femtosecond laser 23, an upstream grating 25, an upstream convex mirror 27, a spatial light modulator 29, a downstream concave mirror 31 , a downstream grating 33, a detection device 35, and a personal computer 37. The pulse shaper is the entire optical setup from grating to grating. As will be discussed later, the gratings can be replaced by a prism, and the spatial light modulator can be a deformable mirror, a liquid crystal spatial light modulator or a microelectronic micromechanical system MEMS. The pulse shaper can also be reflective instead of transmissive. Finally, the pulse shaper can also be an accousto optic shaper.

[0034] Personal computer 37 has a microprocessor based electrical control system, memory, an output screen, a data storage device, an input keyboard, and a storage disk. More specifically, the detection device is a mass spectrometer 39. Bursts or pulses of a laser beam 43 are emitted from laser 23, through the optics 25, 27, 31 and 33, as well as through the spatial light modulator 29 for detection and sensing by spectrometer 39 for further evaluation, analysis, comparison and subsequent control by personal computer 37.

[0035] The laser is preferably an ultra-fast femtosecond laser capable of high peak intensity (with a typical peak greater than 10¹⁴ watts/cm²) which preferably emits laser beam pulses of less than 100 femtosecond duration, and more preferably at or less than 50 femtoseconds, and for the environmental monitoring applications even more preferably as short as 10 femtosecond duration, for each pulse burst or shot. The intense optical pulses are formed in a Kerr-Lens mode locked titanium sapphire oscillator. Such lasers are capable of producing hundreds of nanometers of coherent bandwidth, although only about 50 nm are typically used. The output is amplified in a 1 kHz regenerative chirped pulsed amplifier. The output pulse is typically 100 fs long with a central wavelength of 800 nm and total pulse energy of 0.1 to 1 mJ. Preferred lasers include: the Kapteyn and Murnane femtosecond laser oscillator, which can produce less than 15 fs pulses at 100 MHz; the Hurricane model from Spectra Physics Inc., which is diode pumped and gives 0.8 mJ per pulse with sub-50 fs pulses at 1 kHz; and the CPA-2001+ model from Clark-MXR Inc., which gives 1.3 mJ per pulse with sub-150 fs pulses at 1 kHz, pumping a Clark-MXR Inc. non-collinear parametric amplifier (hereinafter "NOPA") which produces 0.2 mJ per pulse, and is capable of generating sub-20 fs pulses. A Ytterbium Tungstenate laser may also be desirable. A time of flight mass spectrometer (TOF-MS), and a compact Ocean Optics spectrometer, are preferred for detection purposes. The entire system is preferably mounted on a vibration- damping optical table. [0036] A Fourier plane pulse shaper is preferably used with the present invention for the transmissive construction illustrated with this embodiment. Ultra-short laser pulses contain from one to fifty optical cycles, and last only a few femtoseconds. This is much faster than most current electronics and therefore shaping with fast time gates is very difficult. On the other hand, because of the uncertainty principle, the optical spectrum spans tens to hundreds of nanometers. Such a large bandwidth is relatively easy to measure and to filter, and there are several techniques to shape the spectrum in the frequency domain, and thereby shape the temporal pulse upon recompression. [0037] In order to access the frequency domain and the individual frequency components that comprise the pulse, a geometric arrangement is employed, using two back-to-back spectrometers. The spectrometers are especially designed to introduce no net temporal dispersion: that is, all colors pass through the spectrometers within the same amount of time. The first spectrometer (including grating 25 and mirror 27) spreads the unshaped pulse spectrum along a line according to its dispersion function y(α). The light intercepts spatial amplitude and phase mask spatial light modulator 29 at this point. The mask output then forms the entrance to a second spectrometer (including grating 33 and mirror 31) which recombines the colors into a single shaped pulse.

[0038] The heart of the pulse shaper is the programmable 256 pixel liquid-crystal mask (consisting of two overlapping 128 pixel liquid crystal arrays) that is placed at the Fourier plane 29. For the applications envisioned herein, the mask must be capable of shifting the phase of individual frequencies. For alternate embodiment pulse shapers, a different electronically programmable mask that is capable of controlling phase has been demonstrated: a liquid crystal display (hereinafter "LCD"), an acousto-optic modulator (hereinafter "AOM"), a deformable mirror, and a permanently deformed mirror. A LCD pulse shaper can be obtained from CRI Co. and has a modulator electronic driver. [0039] The AOM consists of an anti-reflection coated Tellurium Dioxide

(TeO₂) crystal with a piezo electric transducer glued onto one end. The central frequency of the acoustic wave is αc/2π = 200 MHz. The acoustic velocity vs in the crystal is 4.2 km/s and the light pulse spends less than 10 ps in the crystal, so the acoustic wave moves less than 0.002 λ acoustic during the transit of the light field through the crystal. Since the acoustic wave is essentially frozen as the optical pulse travels through the crystal, the complex amplitude of the acoustic wave traveling through the crystal in the y direction, A(t)cosαct =

A(y/vs)cosαct, is mapped onto the optical field E(α) as it passes through the

AOM. If some of the dispersed optical field encounters a weak acoustic wave, that frequency is attenuated; if the acoustic wave carrier is shifted by phase angle ø, that phase shift is imposed on the optical field. This pulse shaper has a total efficiency of about 20% including the diffraction efficiency of the AOM and the diffraction efficiency of the gratings. The diffracted light is used and the undiffracted "zero order" beam is blocked, to allow full modulation of both amplitude and phase in the shaped beam. The shaped beam then has the form

E_shaped {ω) ^ E_input {ω)xa{ω)xeⁱ^^> where a(ω)e^m<0) - A[y(ω)/v_s]; a is the frequency, and e is a constant.

Fixed pulse shaping optics, such as chirped mirrors or permanently etched reflective masks, can also be employed. The laser pulses are fed into the pulse shaper where the grating 25 causes dispersion. Curved mirror 27 focuses the spectrum onto Fourier plane 29. Changes in the phase ø of the spectral components indicated by the computer are used to tailor the laser pulse before reconstruction with second curved mirror 31 and grating 33. Once shaped, a small portion of the pulse is directed to spectrometer 39 for evaluation and calibration using MIIPS. The greater percentage of the laser, for example 95%, is sent towards the mass spectrometer for monitoring purposes.

[0040] In this embodiment, the phase and amplitude masks of the pulse shaper are controlled by the computer wherein the laser pulse shape takes a dynamic role. The microprocessor within personal computer 37 will then control laser 23, receive an essentially real time feedback input signal from spectrometer 39, and then perform calculations, comparisons and evaluations, and possibly automatic variation of subsequent pulse shapes. These automated steps can be substituted with manual user calculations and decisions if desired based on personal computer outputs.

[0041] As applied to the environmental monitoring applications herein, selective control of multiphoton processes in large molecules, including proteins, is possible using a simple pulse shaping method that is based on taking maximum advantage of the multiphoton intrapulse interference caused in short pulses with large bandwidths, and the resulting fluctuations in the electromagnetic pulse. The results show an extraordinary level of control that is robust and sample independent, with contrast ratios near two orders of magnitude. Such large contrast ratios allow for more precise cancellation control of undesired photons and other laser beam characteristics, such that nonlinear transitions induced by each pulse are controlled. Moreover, a fs-pulse shaper arrangement can be used to achieve selective fragmentation of the sample molecules, and this property can be exploited for identification. The details of the underlying technology will be described in greater detail hereinafter. Alternately, instead of the LCD-SLM used in Figure 5 to introduce phase functions prestored in the memory unit of the controller, phase functions can be incorporated into a passive optical component such as a turret with different static phase masks manufactured on a transparent substrate. The turret can toggle the static phase masks in front of the mirror (see Figure 5). A system with such a turret would reduce the complexity and expense of computer controlled LCD-SLM. Nevertheless, the ability to run a MIIPS test is valuable to make sure that the laser is operating properly, thus, a computer controller pulse shaper will be preferred for this application.

Multiphoton lntrapulse Interference

[0042] A multiphoton intrapulse interference phase scan (hereinafter "MIIPS") system and method are preferably employed with the environmental monitoring system of the present invention to characterize the spectral phase of the femtosecond laser pulses. The phase across the spectrum of an ultrafast pulse can affect the multiphoton process in a number of ways. Phase can increase the pulse length and hence reduce the peak intensity of the pulse, thereby preventing saturation, a common result under high intensity excitation. Phase can also be used to synchronize changes in the electric field with intramolecular wave packet dynamics. Finally, phase can be used to cause interference in the way multiple frequencies combine to achieve multiphoton excitation. This process is known as multiphoton intrapulse interference (hereinafter "Mil").

[0043] The technique of Mil and its application to control multiphoton processes is based on rationally designing an electric field required to achieve a particular target with a minimum number of parameters. The method is based on calculating the amplitude of the nth-order electric field and comparing it to the absorption spectrum of the molecules being controlled. This provides a strong physical understanding of the control process, which can be very useful in the interpretation of experiments where the field is optimized by computer programs based on evolutionary learning or similar methods. For strong field photofragmentation and ionization, required for environmental monitoring, no such calculations are possible. Therefore, determination of the optimum pulse shapes to identify each molecule will be obtained experimentally using a learning method.

[0044] Multiphoton intrapulse interference phase scan is capable of both pulse characterization and compensation of subsequent pulses. Within minutes, the pulses are characterized and compensated to yield transform- limited (hereinafter "TL") or user-specified shaped pulses at the sample. This capability is important for the present invention to make sure the pulse shaper and laser are operating within specifications. [0045] MIIPS is a single-beam method that does not require an interferometer. To make a precise and accurate measurement of the spectral phase using MIIPS, a known phase delay is imposed on the frequencies that make up the pulse using a calibrated pulse shaper. The pulse shaper 121 essentially behaves as two back-to-back spectrometers. In one embodiment, shown in Figure 5, the pulse emitted by laser 123 is dispersed with a grating and collimated with a 200-mm cylindrical lens 133. At the Fourier plane, where all the frequencies are isolated, their phases are manipulated by a computer- controlled LCD^* spatial light modulator (hereinafter "SLM"). The SLM applies the reference phase function to the input pulse. The resulting pulse is then retro reflected and reconstituted to the time domain and the output is picked by mirror 131. The SLM can be updated every pulse (presently limited to 1 kHz). The LCD has a 250-ms response time, so in principle it can be updated at 4 kHz. A small percent of the output beam is reflected by a partially reflective mirror 513 and is analyzed by focusing (with lens 515) onto a 0.10-mm-thick beta barium borate crystal 507 for second-harmonic generation (hereinafter "SHG") in its path, usually at the place where optimum pulses are required. The use of the second harmonic is important to the method. The average SHG output is collimated by second lens 515 and directed to a dispersive spectrometer 503 and detected by a detector array 511. For each reference phase function that is introduced by the computer-controlled SLM, a different spectrum is recorded and stored in computer controller 531. In a sense, the pulse autocorrelates itself at the SHG crystal.

[0046] Pulse characterization involves the introduction of a reference phase-modulation function of the form Φ= ∞ cos {γ Ω - δ), where ∞. is the magnitude of the phase delay, γ is the periodicity Ω is the frequency detuning from the carrier frequency of the pulse, and δ is the position in the spectrum at which the cosine function is equal to one. The reference phase function, with typical values ∞ = 2π, and γ- pulse duration, is programmed into the SLM and scanned for different values of δ ranging from 0 to 2 π. For each value of δ, the spectrum of the frequency- doubled pulse changes, achieving a maximum in the spectral region over which the SLM compensates for the phase distortions. [0047] Qualitatively, the distance between the diagonal features determines linear chirp while the angle between the features determines the quadratic chirp. The full quantitative determination of the spectral phase by integration can be obtained. Once the MIIPS system has characterized the pulse and retrieved the phase distortions inherent to the pulses, it can use that information to drive the SLM such that it compensates for the distortions. The first step in compensation is to take the phase determined from the first scan and program it into the SLM with a negative sign so that it subtracts the distortions. The system carries out a new phase scan to determine the remaining spectral phase modulation (usually about 10% of the original). Typically, three such iterations will yield transform-limited pulses. Because the laser is not focused in the pulse shaper, the method can be used with pulses that are relatively high in energy. Pulses ranging from about 10 pJ to about 30 mJ and pulse durations from less than 5 to about 500 fs can be used. Once the pulses are compensated (transform-limited), the laser can be focused to produce peak intensities from about 10¹² to about 10¹⁸ W/cm², depending on the input energy.

[0048] This single beam method is capable of retrieving the magnitude and sign of second and third order phase modulation (in other words, linear and quadratic chirp) directly, without iteration or inversion procedures. MIIPS achieves accurate phase retrieval from chirped ultrashort pulses. For MIIPS, no synchronous autocorrelation, beam splitting, or time delays are required because the second harmonic spectrum depends on the relative phases of all frequencies within the pulse. The amplitude of the pulse is obtained directly from a spectrometer in a communications receiver. Referring to Figure 5, in order to precisely determine of the phase of all frequency components in a pulse from a femtosecond laser 123, a pulse shaper, such as the one described in A. M. Weiner, "Femtosecond pulse shaping using spatial light modulators," Rev. Sci. Instrum. 71 , pp. 1929-1960 (2000), is employed to introduce a reference phase function designed to yield this information directly, as further described hereinafter. The shaped pulses are frequency doubled by a thin SHG crystal 507 and the output is directed to spectrometer 503. [0049] The MIIPS method is based on the principle that second harmonic generation, as well as other nonlinear optical processes, depend on the phase function φ (ω) across the spectrum of the laser pulse. Higher order phase distortions, such as self-phase modulation and quadratic phase components can be obtained from the curvature of the line defined by the maximum SHG response. The MIIPS can be programmed to find the phase distortions on the laser pulses directly by integration and to introduce a compensation phase function that eliminates the distortions. This mode of operation can be used to find arbitrary phase deformations and yield transform limited pulses, which in a MIIPS scan, look like straight parallel lines separated by π.

[0050] The present invention provides a system and method to characterize the spectral phase of femtosecond pulses. This single beam method is capable of retrieving the magnitude and sign of linear and quadratic chirp with high resolution. Pulse retrieval is based on analytical expressions that yield the phase distortion, without iteration or inversion procedures. Linear and quadratic chirp values, and to some extent cubic chirp values, are important because there are knobs on the laser that can be used to correct for this distortion by mechanically adjusting the grating spacing in the laser beam amplifier compressor. The method can be used with very short pulses. This adjustment can be automatically controlled with the computer controlled software. The method is very versatile, and can be used with high or very low intensity pulses for any wavelength for which low cost, off-the-shelf SHG crystals exist. MIIPS can also be used by obtaining third or higher order harmonics in gases. The maximum signal makes the method useful for the characterization of pulses in wavelength regions for which SHG crystals are not available.

Second - Harmonic Generation with, Powders

[0051] For some variations of the present invention, chemical powders, adhered onto transparent quartz carriers, are employed in place of thin SHG crystals. The powder embodiment is presently preferred to significantly reduce cost in high energy (for example, one nanojoule or greater) applications such as for MIIPS, nonlinear optical characterizations and FROG analysis. The chemical powder is preferably Potassium Dihydrogen Phosphate (KDP or KD^*P) or Beta Barium Borate, which is glued onto a glass, microscope slide using a silicone- rubber or cyanoacrylate (C₅H₅NO2) adhesive. Type I phase matching is used. The power particle size is preferably between about 0.5 to 20 microns, depending on the focusing length, intensity of the laser and the sensitivity of the detector. Furthermore, the powder approach is advantageous by having large variety of random crystal orientations therein which creates improved average results for laser pulses shorter than 20 fs.

Binary Phase Shaping [0052] Laser control is dominated by interference between different nonlinear optical pathways connecting the initial and final states. The challenge is finding the proper phase for each frequency within the pulse to achieve constructive interference at the desired pathway and destructive interference elsewhere. To a very good approximation, it is sufficient to set the phase of different frequency components of femtosecond pulse to two values separated by π. To obtain a very good approximation, it is sufficient to set the phase of different frequency components of a femtosecond pulse to two values separated by π. The phase between photons of different frequencies takes only two values, preferably 0 or π, to maximize or minimize a given pathway. Any two values whose difference is π work equivalently well. This method is defined as binary phase shaping (hereinafter "BPS"). BPS is preferably used to solve the problem of selective multiphoton excitation with ultrashort laser pulses. The use of a small number of phase values between 3 and 10 is regarded as a trivial extension of our binary approach. [0053] In order to control the behavior of molecules under the intense radiation of a femtosecond laser, it important to modulate the phase of the frequencies within the bandwidth of the pulse. For the invention being described, this implies, finding the best binary value for each of the spectral components that traverses a discrete region in the shaper corresponding to a pixel. In some shapers, the pixel may be discrete as in a liquid crystal, or a moving MEMS component. In other cases the pixel may be continuous as a region in an AO crystal or a deformable mirror. The appropriate phases that generate a specific desired outcome that is useful in the identification of a chemical, is found using the learning algorithm illustrated by the flowcharts in Figure 6.

[0054] In another variation of the present invention, specific phase functions designed to produce a specific time-domain fluctuation of the electric field in the time domain will be programmed in the analyzer, and the behavior of all chemical agents of interest evaluated for their behavior when interrogated by these specific phase functions. These functions could be introduced as fixed optics produced by etching or other micromachining methods, or programmed on 128-pixel SLM or a SLM with greater than 256 pixels. [0055] A titanium-sapphire regeneratively amplified laser is preferably employed, which can be obtained from Coherent Lasers, and it is capable of generating pulses as short as 30 fs. The spectral phase of the pulse is tailored using a computer-controlled pulse shaper. Preferably, the pulses are centered near 800 nm. A Ytterbium laser, centered at 1040 nm, may alternately be used. The spectral phase of each pulse is corrected using the Mil phase-scan (MIIPS) method, which compensates phase distortions to obtain transform-limited (TL) pulses. The binary phase is introduced as an addition to the compensation phase. About 1% of the shaped laser pulse intensity, with energy ~ 0.003 mJ per pulse and 1 kHz repetition rate, is focused, to a spot size of ~ 20 microns in diameter, on a 100 micron thin beta barium borate (βBBO) type I SHG crystal. The frequency-doubled light is collected with an optical fiber and dispersed on a compact spectrometer, preferably obtainable from Ocean Optics. This setup is used at regular intervals to ascertain that the laser system and pulse shaper are operating properly. Most of the intensity of the laser beam with energy -0.3 mJ is directed towards the TOF-MS where it interacts with the sampled air, causing photofragmentation and ionization.

[0056] The advantage of BPS is that computational redundancies are greatly reduced. For BPS and 128 active pixels, the search space is reduced by hundreds of orders of magnitude compared to arbitrary phase and amplitude pulse shaping as discussed above. The resulting space is small enough that a greater percentage of the search space can be evaluated experimentally. A learning feedback method or simple evolutionary learning computer program can quickly converge towards significantly improved solutions. BPS is expected to have significant technological advantages. A retardation equivalent to π is easy and fast to obtain and calibrate. Permanently etched masks can be made in advance and used for specific applications. Library Evolutionary Learning Computer Program

[0057] Initially, the pulse shaper is programmed for monitoring the environment based on requirements of a library of data stored in memory of the computer. Fast, accurate and reproducible monitoring of the environment for chemical and biological agents, will not usually permit the pulse shaper to run in an evolutionary learning mode in the field, although such field use still falls within the scope of the present invention. Instead, a much more efficient mode of operation is used, whereby the unit is continuously monitoring the environment with transform-limited (TL) pulses. Under these circumstances, the unit takes a fraction of a second to make preliminary, but highly sensitive measurements. In the event that a suspected molecular ion fragment is detected, then the unit activates a library search for the pre-stored suspected agents and employs a series of pre-programmed shaped pulses to make an absolute identification. While in the search mode, the unit can speed by three to six orders of magnitude in repetition rate to provide the fastest analysis possible. Because it is using a pre-programmed series of shaped pulses, the whole identification process can be completed in under a minute.

[0058] The evolutionary learning program, shown in flow chart 6B and used in a laboratory setting to define a library of pulses, begins with an initial random population of pulse shapes, each characterized by binary code that specifies the phase for each pixel in the pulse shaper device. Each pulse shape is tested for its fitness (ability to generate the result that most resembles a target selected in advance). New sets of pulse parameters are created by modifying (mutating) and recombining (via crossover) elements of the best of the previous pulse shapes to create new ones (survival of the fittest). This basic series of processes is iterated and the fitness converges toward a "best" value.

[0059] In simple terms, the role of the pulse shaper is to advance or retard individual frequencies within a laser pulse by a factor of π. For example, in a TL pulse, all the frequencies are locked, and have net zero retardation. In this case, the spectral phase is flat. The pulse shaper can be used on a TL pulse to make some frequencies arrive before others. In these general terms, a pulse shaper can be defined by a number of parameters: input bandwidth, frequency resolution, and maximum retardance. The spectral resolution of a pulse-shaper setup is determined at the Fourier plane, where the pulse is resolved in the frequency domain. The phase retardation must be calibrated and checked for accuracy and reproducibility. This step will be done by MIIPS as described earlier. [0060] For example, the pulse shaper shown in Figure 5 is in the 2f reflective configuration is provided in combination with a dispersion grating, a cylindrical mirror, the phase-retardation unit, which, in this case, is a liquid crystal spatial light modulator (LC-SLM), and a retro reflector mirror. The SLM unit is located at the Fourier plane of the shaper. As mentioned above, at the Fourier plane, the pulses are spread into their frequency components where each can be delayed separately. Experimentally, spatial light modulators have a finite resolution that is determined by the number of individual liquid crystal elements (pixels) and the optics that disperse the spectrum over the SLM. The liquid crystal is birefringent, therefore, depending on the polarization of the incoming light, a voltage can introduce pure phase retardation or a combination of phase retardation and polarization rotation. The key parameter for determining the optical resolution is the focal spot size that is projected by the lens across the SLM. The optical setup must take into account the effective confocal parameter. The SLM unit is about 1 cm thick such that a setup with a sufficiently long confocal parameter is needed to preserve the spectral resolution.

[0061] The frequency resolution of the pulse shaper is determined by the number of pixels in the SLM, as described above. Frequency resolution plays two different roles. First, one can think of frequency resolution in the time domain. Because of the Fourier transform relation between spectral and time resolution, the higher the spectral resolution, the longer the pulses that can be produced by the pulse shaper and the higher the finesse with which one can control narrow frequency resonances in the sample. This is of particular importance for samples that have narrow resonances, like gas-phase systems. The number of active pixels can be related to the overall complexity of the resulting shaped pulse. 128 active pixels is adequate for effecting molecular control and influencing ionization and fragmentation. [0062] When the optical axis of the liquid crystal mask is oriented at a

45° angle with respect to the polarization of the incident electric field, polarization rotation is introduced in addition to retardance. When two such SLM units are lapped back-to-back and in opposite angles of rotation, and they are flanked by input and output polarizers, one can control phase and amplitude of the transmitted light. The double-mask SLM unit operation is illustrated in Figures 9- 11. Referring to Figure 10, a light incident from the left goes through a polarizer and then transmits through the first liquid crystal plate, where it undergoes retardation and polarization rotation. When the light reaches the second liquid crystal plate, it undergoes retardation, but this time polarization rotation takes place in the opposite direction. Finally, a second polarizer is used for amplitude control. When the polarization is rotated away from the incident, horizontal direction, it is attenuated by the second polarizer. A rotation of 90 degrees results in zero transmission. This polarization dependence is very useful for calibration of the pulse shaper. Ramping the voltage of one of the liquid crystal plates while maintaining the other at a constant voltage results in the transmission function, shown in Figure 11 , which can be used to accurately calibrate the dependence of retardance on voltage. The total retardation φ is determined experimentally taking advantage of changes in the transmission given by: T= cos²[π(f?i ( V_Λ )-R₂( V₂))! λ],

where R(V) is the retardance as a function of voltage introduced by each SLM unit. By fixing V₂ and scanning V_\, T(V^) is measured and Ri(Vi) is calculated by the computer (except for a constant). R₂(V₂) is obtained by measuring T(V₂) while keeping Vi constant. Knowing /?i(Vi) and R₂(V₂), the phase delay φ is calculated according to: φ = π(Ri(Vi)+R₂(V₂))/λ). It is possible to program the pulse shaper through the computer to provide amplitude-only or phase-only pulse shaping, using the formulas given above. The pulse shaper can be programmed to automatically perform a full calibration, a process that usually takes a matter of minutes. When the LCD in the elements is not oriented at 45 degree angle, a half-wavelength plate can be used to rotate the polarization of the beam and perform the calibration.

[0063] Once the shaping unit is calibrated, it is imperative that the pulse shaper compensates the phase deformations in the femtosecond laser system. Even under ideal conditions, a commercial femtosecond laser produces pulses with phase deformations that are primarily quadratic or cubic in the frequency domain. Multiphoton lntrapulse interference phase scan (MIIPS) automatically determines the spectral phase deformations and compensates for them. Within a minute, the phase distortions are eliminated and TL pulses are obtained (see Figure 12) where distorted pulses are shown (before) and compensated pulses are obtained after MIIPS. TL pulses are identified by MIIPS by producing straight parallel lines separated by π. The MIIPS method is believed to be at least one order of magnitude more accurate than alternative methods commercially available. The shaper calibration described above, together with accurate phase characterization and compensation provided by MIIPS, are important steps that are believed to ensure the robustness and reproducibility sought by the system of the present invention.

[0064] The monitoring system depends on the identification of certain shaped laser fields that will produce different uniquely identifiable fingerprints from each chemical or biological compound. The search for these shaped laser fields requires the search of a large parameter space using the evolutionary learning program. BPS reduces the search space and drastically increases the reproducibility. Each chemical will be entrained in a molecular beam, which will be interrogated at right angles by the focused shaped laser field. The mass spectrum resulting from ionization and fragmentation will be recorded and compared to others. The evolutionary learning program will have a search target and will determine the shape that best approaches the target. [0065] The ability of strong shaped laser fields to influence the ionization and fragmentation of molecules is utilized in order to find uniquely identifying pulses shapes for each molecule. Because laser-molecule control depends on the electronic and nuclear structure of the molecule, no two compounds would yield the same results. Reference should be made to Figure

13. Specifically, the peak with the largest mass (a), and the subsequent two highest peaks (b) and (c), respectively, are labeled. Six targets are optimized.

The first three are those that increase the following peak to valley ratios (a/b),

(a/c), and (b/c). The final three laser shapes are those that generate the lowest values for the same ratios. In order to maximize the identification capabilities of the method, the six pulses are required to produce identifiably different results with the six different pulse shapes and different regarding the ionization fragmentation pattern obtained by excitation with TL pulses. This scheme is outlined as an illustration for the strategy to be used for identifying the pulses to be included in a library, but other strategies could be conceived.

[0066] An evolutionary learning computer program (ELC), also known as a learning feedback method, is initially employed in the environmental monitoring system of the present invention to determine and pre-store a library of the undesired agents and the acceptable background chemicals. The phase masks of the pulse shaper are controlled by computer 37 (see Figure 3) and 531 (see Figure 5). The adaptive laser source is part of a learning feedback method or ELC that modifies the laser pulse shape based on its success at optimizing the yield of charged agents, which may include chemicals or proteins. In the present application, the laser pulse shape takes a dynamic role. Pulse shapes are envisioned which include sequences of pulses where each pulse in the sequence plays a different role, for example, melting, excitation, selective fragmentation, proton transfer and evaporation.

[0067] The physical process runs itself by an intelligent "feedback" method by means of an intelligent loop. The learning method tries various pulse shapes, assesses their success in achieving the desired target excitation, and uses the knowledge gained in this way to improve the pulse shapes on subsequent laser shots, all with only minimal intervention of the researcher or system user. Changing conditions are automatically corrected within the learning method or feedback loop.

[0068] Reference should now be made to Figures 6-8B. The feedback software is implemented with an initial random population of pulse shapes. Each pulse shape is characterized by a series of numbers that specify the spectral phase and amplitude in each of the wavelengths or frequencies within the pulse. This particular parameterization for the pulse shape is itself subject to optimization; this makes the algorithm adaptive. The most important part of the method is the test for fitness of a given pulse shape. Each pulse shape is tested for its ability to generate the result that most resembles a target selected in advance. For example, the program will calculate the ratio between the amplitude at the desired protein signal and the background. Once relative success is quantified, a new generation of pulse shapes is produced by mating different parts of the amplitude and phase information from pairs of the most fit pulse shapes from the current generation. In addition, the method prescribes a small probability (15%) of random changes or mutations in the successful pulse shapes. Furthermore, a new set (10%) of pulses is generated by setting random portions of the phase mask to zero, thus generating new phase masks. Once the new set of individual phases is generated they are evaluated experimentally. This basic series of processes is iterated until the fitness converges on a "best" value.

[0069] The convergence and robustness of the feedback method solutions can be measured in two different ways. First, the variance in the amplitude and phase information itself can be monitored. As the feedback method converges on a solution, the values fall into a narrow range that produces the best result. Over the course of many generations, some of the parameters become very stable, which is an indication that those spectral phases and amplitudes are particularly important for driving the process that determines fitness. Secondly, the information for different initial conditions is monitored. If the feedback method is working it should converge to a similar result. [0070] New sets of parameters are created by modifying and combining elements of previous pulse shapes to create new ones. This is carried out by statistical operators that act on the phases and amplitudes of the pulse shapes. Operators that can be used include multi-point crossover, mutation, as discussed above.

[0071] A well-chosen set of operators can greatly enhance the performance of the feedback method and lend additional physical insight. However, the proper choice is usually far from obvious, so the method is allowed to adapt itself by letting it choose how often to use a given operator to produce new pulse shapes. The use of adaptive operators helps speed up convergence, and, perhaps more importantly, it helps shed light on the control mechanism at work. For example, crossover is more effective in the beginning of the algorithm when there is maximal uncertainty, since it does a good job of mixing up the information of the initial pulse shapes. It becomes less effective as the feedback method converges to the best solutions, since at this point there is much less change in the parameters, so there is no longer a need to drastically change the information. Ideally, the learning program learns from its past mistakes and does not test possible pulse shapes which it now knows will fail, which saves a considerable amount of computing time.

[0072] For each pulse shape, a number of spectra will be obtained. The number of laser shots that are averaged per pulse shape will depend on achievement of a statistically significant spectrum. At first, when the pulse shapes are the result of random phases and amplitudes we imagine that up to 1000 repetitions may be needed to distinguish the more efficient pulse shapes. This will allow one pulse shape per second. As the selection process proceeds large gains in efficiency can be expected. The final stages of the optimization may be carried out at a rate of 100 different pulse shapes per second. The goal is to reach single pulse, femto-mol sensitivity. The second goal is to define a series of such pulses that generates uniquely different spectra from the same compound. It is the predictable different behavior that a compound exhibits when interrogated by each of the optimized shaped pulses, the response patter, that will be used for positive identification. A given response that will stand out even in the presence of a chemically complex environment.

[0073] The learning feedback software employed in the present invention control system and apparatus is more fully described as follows. The preliminary investigation method and computer software steps for analyzing a pre-test unknown sample or specimen can be observed in Figure 7. For any new system, the test should start with pre-defined pulse shapes in order to obtain a basic understanding of the system. Among the pre-defined pulses, the shortest pulse is expected to ionize molecules on the surface of the sample with minimum decomposition, the longest pulse is expected to mimic the nanosecond experiments where the chemical undergoes substantial fragmentation. The microprocessor within personal computer 37 (see Figures 3 and 4) will then control laser 23 in accordance with step D, receive an essentially real time feedback input signal from mass spectrometer 39 in accordance with step F and then perform calculations, comparisons and evaluations in accordance with steps G, H and I. These automated steps can be substituted with manual user calculations and decisions if desired based on personal computer outputs.

[0074] The objective of the software routine of Figures 6B, 8A and 8B is to aid in the selection of unique identification shaped pulses that will become integral part of the library for field monitoring. An optional alternate embodiment subroutine includes shooting a long nanosecond laser beam pulse on a solid target to cause a layer of the solid to evaporate. The evaporated compounds are then analyzed by the shaped laser pulses. The short pulse of approximately 50 femtoseconds is shaped and focused on the plume of gases in order to look for fragmentation and sample identification. The long and short pulse combinations can be used in addition to or without the benefit of pulse shaping. Otherwise, the control system and apparatus are the same as discussed herein. This combination is specially useful for the study of non-volatile compounds such as spores and other biological agents. Another alternate variation is where particulate matter in the air is deposited on a sticky matrix and then a nanosecond laser is used to desorb it and then the fs laser is used to analyze the particles for biological agents such as spores. This would pertain to the case where two pulses are used, the first one is a nanosecond or femtosecond pulse and is only used for desorption. The second pulse after an adjustable delay time, is used for fragmentation and further ionization.

[0075] The real time learning feedback method and computer software will now be described in greater detail. This method and software are employed to statistically optimize the repetitive identification of molecularly complex and unknown samples or specimens in a highly automated and relatively quick manner. The data that is obtained in the mass spectrometer for each laser pulse consists of a two-dimensional array of numbers which are stored in random access memory of the personal computer. The first column of the array contains arrival times for the ions, when the data is obtained from the time-of-flight mass spectrometer. Equivalent numbers can be obtained from different mass spectrometry units such as quadrupole or ion-cyclotron spectrometers. This numbers can be converted to a mass assuming a charge for the species. The second number identifies the shaped pulse that caused the specific mass spectrum. The data acquisition involves collecting a data set that unequivocally identifies a chemical compound even when it is a small percent of a complex mixture.

[0076] Each pulse shape is defined by a one-dimensional array of numbers which are stored in random access memory of the personal computer or otherwise accessible on a read only basis by the personal computer from a disk or the like. The length of the array determines the resolution with which the spectrum of the laser pulse is being sculpted. For a liquid spatial-light modulator, the number of pixels typically determines this number. The array of numbers determines the phase delay for every pixel. The entire array of phases determines the final shape of the output pulse. The shortest pulse is pre-defined as the shortest duration possible for the laser system will the longest pulse is pre-defined as the longest pulse that can be made with the pulse shaper. A two- pulse combination is pre-defined as the combination of an unshaped pulse with a shaped pulse. Pre-defined ultraviolet or infrared pulses at 400 or 800 nm, for example, can be used. [0077] Step R allows the microprocessor to determine the best pulse shape with the highest success value and store it as the optimum value in the random access memory of the central processing unit. The computer will then pick approximately the ten percent best pulse shapes based on the highest success values and then reset the discarded values and automatically generate new laser pulse shapes in step V for subsequent testing iterations on the same specimen. The generation of new seed pulses is important for the success of the feedback method. The goal is to arrive at the optimum pulse in the shortest number of iterations making sure that one has searched the entire range of parameters, the global maximum. The "cost functional" refers to the statistical pressure that is placed on the optimum pulse shape in order to simplify it. For example, once an optical pulse shape or other characteristic is found, it may be important to determine how sensitive the outcome is to each of its amplitude and phase components. Perhaps a much lower resolution can produce the results. The simpler the pulse the easier it is to reproduce and interpret the results in terms of physical concepts. For certain cases, the shape can be simple enough that it can be prepared without a shaper which would allow for a less expensive alternative to the preferred pulse shaping; for example, a combination of two or three different 800 nm pulses, or a combination of infra-red and ultra-violet pulses could be employed as a modified or optimized pulse. Once statistical convergence has been determined by the personal computer, then the test is complete by determining the optimum pulse characteristics (whether they be pulse shape, pulse duration or any other such variable laser beam characteristic) for the corresponding and now post-test identified specimen.

[0078] Subsequently, Figure 6C illustrates the computer program and method for inclusion of a desired phase shape for acceptable background items and unacceptable agents. The ten statistically best pulse/phase shaping solutions are used which identify the chemical compound. These are then automatically ranked by the computer controller based on recognizable differences with previously stored phase shapes. These are then, in turn, ranked based on differences to other previously stored mass spectrum data. Next, the controller automatically selects for storage in memory, the new phase and pulse shaping characteristics for a specific chemical compound or agent that do not overlap with previously stored data. This method significantly reduces the chance of erroneous identifications and determinations.

Alternate Applications

[0079] Figure 15 illustrates a first alternate embodiment use of a laser and environmental monitoring system 701 of the present invention. In this embodiment, a monitoring unit 703, having a femtosecond laser, pulse shaper, mass spectrometer and computer, are fluidically coupled to various inlet ports 705 within walls of a metal detector, security checkpoint unit 707. Such a device is installed within an airport, building entry way, train station or other location of public access. The monitoring unit is set up to sense and determine if explosives or other undesired chemicals are present on the person, luggage or other item being moved adjacent inlet ports 705. A read out of the sensed undesirable chemical or biological (including a disease) agent is projected upon a CRT output screen 709 and warning lights or alarms can be activated. Active baseline subtraction and library comparisons are employed as was previously disclosed.

[0080] A second alternate embodiment is illustrated in Figure 16. In this exemplary embodiment, an environmental monitoring unit 731 is employed to sense an undesired disease or illness agent such as diabetes, hepatitis, SARS, or for drug and alcohol testing. The person first breathes into a balloon 733 or directly into an inlet tube 735. If the balloon variation is employed, the balloon is attached to inlet tube 735 which then flows the breathed out air containing various chemicals into the mass spectrometer contained within monitoring unit 731. Through use of a femtosecond laser, pulse shaper, computer and previously stored library of data, in a manner previously disclosed herein, an output display 737 will automatically disclose if any of the undesired agents are present.

[0081] While various embodiments have been disclosed herein, it should be appreciated that other modifications may be made that are covered by the system and methods of the present invention. For example, alternate lasers, chemicals, optics, computer controllers and HVAC devices can be employed as long as they function as described. The description of the invention is merely exemplary in nature and, thus, variations that do not depart from the gist of the invention are intended to be within the scope of the invention. Such variations are not to be regarded as a departure from the spirit and scope of the invention.