Ionization and dissociation of H2 in intense laser fields at 1.064 microm, 532 nm, and 355 nm

Efficient Control of Electron Localization and Probability Modulation with Synthesized Two-Color Intense Laser Pulses

A coupled electron-nuclear dynamical study at attosecond time scale is performed on the HD+ and H2+ molecular ions under the influence of synthesized intense two-color electric fields. We have employed ω - 2ω and also, ω - 3ω two-color fields in the infrared/mid-infrared regime to study the different fragmentation processes originating from the interference of n - (n + i) (i = 1, 2) photon absorption pathways. The branching ratios corresponding to different photofragments are controlled by tuning the relative phase as well as intensity of the two-color pulses, while the effect of the initial nuclear wave function is also studied by taking an individual vibrational eigenstate or a coherent superposition of several eigenstates of HD+ and H2+. By comprehensive analysis, the efficacy of the two different types of synthesized two-color pulses (ω - 2ω and ω - 3ω) are analyzed with respect to one-color intense pulses in terms of controlling the probability modulation and electron localization asymmetry and compared with previous theoretical calculations and experimental findings. Through the detailed investigation, we have addressed which one is the major controlling knob to have better electron localization as well as probability modulation.

Electro-Nuclear Dynamics of Single and Double Ionization of H2 in Ultrafast Intense Laser Pulses

We present an efficient method for modeling the single and double ionization dynamics of the H2 molecule in ultrashort, intense laser fields. This method is based on a semianalytical approach to calculate the time-dependent single and double molecular ionization rates and on a numerical approach to describe the vibrational motion that takes place in the intermediate molecular ion H2+. This model allows for the prediction of the single and double ionization probabilities of the H2 molecule to be made over a wide range of frequencies and laser intensities with limited computational time while providing a realistic estimate of the energy of the products of the dissociative ionization and of the Coulomb explosion of the H2 molecule. The effect of vibrational dynamics on ionization yields and proton kinetic energy release spectra is demonstrated and, in the case of the latter, is discussed in terms of basic strong-field molecular fragmentation mechanisms.

Rabi oscillations in a stretching molecule

Rabi oscillation is an elementary laser-driven physical process in atoms and artificial atoms from solid-state systems, while it is rarely demonstrated in molecules. Here, we investigate the bond-length-dependent Rabi oscillations with varying Rabi frequencies in strong-laser-field dissociation of H2+. The coupling of the bond stretching and Rabi oscillations makes the nuclei gain different kinetic energies while the electron is alternatively absorbing and emitting photons. The resulting proton kinetic energy spectra show rich structures beyond the prediction of the Floquet theorem and the well-accepted resonant one-photon dissociation pathway. Our study shows that the laser-driven Rabi oscillations accompanied by nuclear motions are essential to understanding the bond-breaking mechanism and provide a time-resolved perspective to manipulate rich dynamics of the strong-laser-field dissociation of molecules.

Isolating Attosecond Electron Dynamics in Molecules where Nuclei Move Fast

Capturing electronic dynamics in real time has been the ultimate goal of attosecond science since its beginning. While for atomic targets the existing measurement techniques have been thoroughly validated, in molecules there are open questions due to the inevitable copresence of moving nuclei, which are not always mere spectators of the phototriggered electron dynamics. Previous work has shown that not only can nuclear motion affect the way electrons move in a molecule, but it can also lead to contradictory interpretations depending on the chosen experimental approach. In this Letter we investigate how nuclear motion affects and eventually distorts the electronic dynamics measured by using two of the most popular attosecond techniques, reconstruction of attosecond beating by interference of two-photon transitions and attosecond streaking. Both methods are employed, in combination with ab initio theoretical calculations, to retrieve photoionization delays in the dissociative ionization of H_{2}, H_{2}→H^{+}+H+e^{-}, in the region of the Q_{1} series of autoionizing states, where nuclear motion plays a prominent role. We find that the experimental reconstruction of attosecond beating by interference of two-photon transitions results are very sensitive to bond softening around the Q_{1} threshold (27.8 eV), even at relatively low infrared (IR) intensity (I_{0}∼1.4×10^{11}  W/cm^{2}), due to the long duration of the probe pulse that is inherent to this technique. Streaking, on the other hand, seems to be a better choice to isolate attosecond electron dynamics, since shorter pulses can be used, thus reducing the role of bond softening. This conclusion is supported by very good agreement between our streaking measurements and the results of accurate theoretical calculations. Additionally, the streaking technique offers the necessary energy resolution to accurately retrieve the fast-oscillating phase of the photoionization matrix elements, an essential requirement for extending this technique to even more complicated molecular targets.

Circularly polarized light-induced potentials and the demise of excited states

In the presence of strong electric fields, the excited states of single-electron molecules and molecules with large transient dipoles become unstable because of anti-alignment, the rotation of the molecular axis perpendicular to the field vector, where bond hardening is not possible. We show how to overcome this problem by using circularly polarized electromagnetic fields. Using a full quantum description of the electronic, vibrational, and rotational degrees of freedom, we characterize the excited electronic state dressed by the field and analyze its dependence on the bond length and angle and the stability of its vibro-rotational eigenstates. Although the dynamics is metastable, most of the population remains trapped in this excited state for hundreds of femtoseconds, allowing quantum control. Contrary to what happens with linearly polarized fields, the photodissociation occurs along the initial molecular axis, not perpendicular to it.

Anti-alignment driven dynamics in the excited states of molecules under strong fields

We develop two novel models of the H2+ molecule and its isotopes from which we assess quantum-mechanically and semi-classically whether the molecule anti-aligns with the field in the first excited electronic state. The results from both models allow us to predict anti-alignment dynamics even for the HD+ isotope, which possesses a permanent dipole moment. The molecule dissociates at angles perpendicular to the field polarization in both the excited and the ground electronic state, as the population is exchanged through a conical intersection. The quantum mechanical dispersion of the initial state is sufficient to cause full dissociation. We conclude that the stabilization of these molecules in the excited state through bond-hardening under a strong field is highly unlikely.

Clocking Enhanced Ionization of Hydrogen Molecules with Rotational Wave Packets

Laser-induced rotational wave packets of H_{2} and D_{2} molecules were experimentally measured in real time by using two sequential 25-fs laser pulses and a reaction microscope. By measuring the time-dependent yields of the above-threshold dissociation and the enhanced ionization of the molecule, we observed a few-femtosecond time delay between the two dissociation channels for both H_{2} and D_{2}. The delay was interpreted and reproduced by a classical model that considers enhanced ionization and thus additional interaction within the laser pulse. We demonstrate that by accurately measuring the phase of the rotational wave packet in hydrogen molecules we can resolve dissociation dynamics which is occurring within a fraction of a molecular rotation. Such a rotational clock is a general concept applicable to sequential fragmentation processes in other molecules.

Generalized Phase Sensitivity of Directional Bond Breaking in the Laser-Molecule Interaction

We establish a generalized picture of the phase sensitivity of laser-induced directional bond breaking using the H_{2} molecule as the example. We show that the well-known proton ejection anisotropy measured with few-cycle pulses as a function of their carrier-envelope phases arises as an amplitude modulation of an intrinsic anisotropy that is sensitive to the laser phase at the ionization time and determined by the molecule's electronic structure. Our work furthermore reveals a strong electron-proton correlation that may open up a new approach to experimentally accessing the laser-sub-cycle intramolecular electron dynamics also in larger molecules.

Probing multiphoton light-induced molecular potentials

The strong coupling between intense laser fields and valence electrons in molecules causes distortions of the potential energy hypersurfaces which determine the motion of the nuclei and influence possible reaction pathways. The coupling strength varies with the angle between the light electric field and valence orbital, and thereby adds another dimension to the effective molecular potential energy surface, leading to the emergence of light-induced conical intersections. Here, we demonstrate that multiphoton couplings can give rise to complex light-induced potential energy surfaces that govern molecular behavior. In the laser-induced dissociation of H2+, the simplest of molecules, we measure a strongly modulated angular distribution of protons which has escaped prior observation. Using two-color Floquet theory, we show that the modulations result from ultrafast dynamics on light-induced molecular potentials. These potentials are shaped by the amplitude, duration and phase of the dressing fields, allowing for manipulating the dissociation dynamics of small molecules.

Imaging an isolated water molecule using a single electron wave packet

Observing changes in molecular structure requires atomic-scale Ångstrom and femtosecond spatio-temporal resolution. We use the Fourier transform (FT) variant of laser-induced electron diffraction (LIED), FT-LIED, to directly retrieve the molecular structure of H₂O⁺ with picometer and femtosecond resolution without a priori knowledge of the molecular structure nor the use of retrieval algorithms or ab initio calculations. We identify a symmetrically stretched H₂O⁺ field-dressed structure that is most likely in the ground electronic state. We subsequently study the nuclear response of an isolated water molecule to an external laser field at four different field strengths. We show that upon increasing the laser field strength from 2.5 to 3.8 V/Å, the O-H bond is further stretched and the molecule slightly bends. The observed ultrafast structural changes lead to an increase in the dipole moment of water and, in turn, a stronger dipole interaction between the nuclear framework of the molecule and the intense laser field. Our results provide important insights into the coupling of the nuclear framework to a laser field as the molecular geometry of H₂O⁺ is altered in the presence of an external field.

Conical Intersections Induced by Quantum Light: Field-Dressed Spectra from the Weak to the Ultrastrong Coupling Regimes

In classical laser fields with frequencies resonant with the electronic excitation in molecules, it is by now known that conical intersections are induced by the field and are called light-induced conical intersections (LICIs). As optical cavities have become accessible, the question arises whether their quantized modes could also lead to the appearance of LICIs. A theoretical framework is formulated for the investigation of LICIs of diatomics in such quantum light. The eigenvalue spectrum of the dressed states in the cavity is studied, putting particular emphasis on the investigation of absorption spectra of the Na₂ molecule, that is, on the transitions between dressed states, measured by employing a weak probe pulse. The dependence of the spectra on the light-matter coupling strength in the cavity and on the frequency of the cavity mode is studied in detail. The computations demonstrate strong nonadiabatic effects caused by the appearing LICI.

Strong-Field Fourier Transform Vibrational Spectroscopy of D_{2}^{+} Using Few-Cycle Near-Infrared Laser Pulses

The photoionization of D_{2} and dissociation of the resultant D_{2}^{+} are monitored by pump-probe measurements using intense near-infrared few-cycle laser pulses. The yields of D_{2}^{+} and D^{+} recorded up to the pump-probe delay time of 527 ps exhibit oscillatory structures reflecting the motion of the created vibrational wave packet of D_{2}^{+}, and the Fourier transform of the data in time domain reveals the vibrational level separations with uncertainties of 0.0002-0.0097  cm^{-1}, showing a potential application of the strong-field pump-probe measurements to high-resolution spectroscopy of molecular ions.

Direct Signatures of Light-Induced Conical Intersections on the Field-Dressed Spectrum of Na2

Rovibronic spectra of the field-dressed homonuclear diatomic Na₂ molecule are investigated to identify direct signatures of the light-induced conical intersection (LICI) on the spectrum. The theoretical framework formulated allows the computation of the (1) field-dressed rovibronic states induced by a medium-intensity continuous-wave laser light and the (2) transition amplitudes between these field-dressed states with respect to an additional weak probe pulse. The field-dressed spectrum features absorption peaks resembling the field-free spectrum as well as stimulated emission peaks corresponding to transitions not visible in the field-free case. By investigating the dependence of the field-dressed spectra on the dressing-field wavelength, in both full- and reduced-dimensional simulations, direct signatures of the LICI can be identified. These signatures include (1) the appearance of new peaks and the splitting of peaks for both absorption and stimulated emission and (2) the manifestation of an intensity-borrowing effect in the field-dressed spectrum.

Insight into the role of UV-irradiation in photothermal catalytic Fischer–Tropsch synthesis over TiO2 nanotube-supported cobalt nanoparticles

To explore an efficient catalytic system with high activity and selectivity is the key to improve Fischer–Tropsch synthesis (FTS) technology and the main focus in the academic field.

The importance of Rydberg orbitals in dissociative ionization of small hydrocarbon molecules in intense laser fields

Much of our intuition about strong-field processes is built upon studies of diatomic molecules, which typically have electronic states that are relatively well separated in energy. In polyatomic molecules, however, the electronic states are closer together, leading to more complex interactions. A combined experimental and theoretical investigation of strong-field ionization followed by hydrogen elimination in the hydrocarbon series C₂D₂, C₂D₄ and C₂D₆ reveals that the photofragment angular distributions can only be understood when the field-dressed orbitals rather than the field-free orbitals are considered. Our measured angular distributions and intensity dependence show that these field-dressed orbitals can have strong Rydberg character for certain orientations of the molecule relative to the laser polarization and that they may contribute significantly to the hydrogen elimination dissociative ionization yield. These findings suggest that Rydberg contributions to field-dressed orbitals should be routinely considered when studying polyatomic molecules in intense laser fields.

Toward the Full Quantum Dynamical Description of Photon Induced Processes in D2. J Phys Chem A 2016;120:9411-9421. [PMID: 27934332 DOI: 10.1021/acs.jpca.6b09623]

The dissociative ionization (multiphoton regime) of the D₂⁺ ion by ultrashort laser pulses has been studied theoretically using ab initio calculations. The combined ionization and dissociation spectrum was explored for fixed molecular axis orientations. In accordance with previous investigations, the dominant features in the obtained joint energy spectrum were multiphoton peaks. In addition to this, in the present work, photoelectron angular distributions were analyzed as well. By performing a partial wave analysis for each multiphoton peak, we have identified the number of absorbed photons. Moreover, we also found that the angular distribution can significantly change inside a multiphoton peak as a function of electron and nuclear kinetic energy.

Molecular vibrational trapping revisited: a case study with D2

The present theoretical study is concerned with the vibrational trapping or bond hardening, which is a well-known phenomenon predicted by a dressed state representation of small molecules like and in an intense laser field. This phenomenon is associated with a condition where the energy of the light induced, vibrational level coincides with one of the vibrational levels on the field-free potential curve, which at the same time maximizes the wave function overlap between these two levels. One-dimensional numerical simulations were performed to investigate this phenomenon in a more quantitative way than has been done previously by calculating the photodissociation probability of for a wide range of photon energy. The obtained results undoubtedly show that the nodal structure of the field-free vibrational wave functions plays a decisive role in the vibrational trapping, in addition to the current understanding of this phenomenon.

Isotope Effect in Tunneling Ionization of Neutral Hydrogen Molecules

It has been recently predicted theoretically that due to nuclear motion light and heavy hydrogen molecules exposed to strong electric field should exhibit substantially different tunneling ionization rates [O. I. Tolstikhin, H. J. Worner, and T. Morishita, Phys. Rev. A 87, 041401(R) (2013)]. We studied that isotope effect experimentally by measuring relative ionization yields for each species in a mixed H_{2}/D_{2} gas jet interacting with intense femtosecond laser pulses. In a reaction microscope apparatus, we detected ionic fragments from all contributing channels (single ionization, dissociation, and sequential double ionization) and determined the ratio of total single ionization yields for H_{2} and D_{2}. The measured ratio agrees quantitatively with the prediction of the generalized weak-field asymptotic theory in an apparent failure of the frozen-nuclei approximation.

Observation of Quantum Interferences via Light-Induced Conical Intersections in Diatomic Molecules

We observe energy-dependent angle-resolved diffraction patterns in protons from strong-field dissociation of the molecular hydrogen ion H_{2}^{+}. The interference is a characteristic of dissociation around a laser-induced conical intersection (LICI), which is a point of contact between two surfaces in the dressed two-dimensional Born-Oppenheimer potential energy landscape of a diatomic molecule in a strong laser field. The interference magnitude and angular period depend strongly on the energy difference between the initial state and the LICI, consistent with coherent diffraction around a cone-shaped potential barrier whose width and thickness depend on the relative energy of the initial state and the cone apex. These findings are supported by numerical solutions of the time-dependent Schrödinger equation for similar experimental conditions.

Direct Signature of Light-Induced Conical Intersections in Diatomics

Nonadiabatic effects are ubiquitous in physics, chemistry, and biology. They are strongly amplified by conical intersections (CIs), which are degeneracies between electronic states of triatomic or larger molecules. A few years ago it was revealed that CIs in molecular systems can be formed by laser light, even in diatomics. Because of the prevailing strong nonadiabatic couplings, the existence of such laser-induced conical intersections (LICIs) may considerably change the dynamical behavior of molecular systems. By analyzing the photodissociation process of the D2+ molecule carefully, we found a robust effect in the angular distribution of the photofragments that serves as a direct signature of the LICI, providing undoubted evidence of its existence.

Electron-nuclear energy sharing in above-threshold multiphoton dissociative ionization of H2

We report experimental observation of the energy sharing between electron and nuclei in above-threshold multiphoton dissociative ionization of H2 by strong laser fields. The absorbed photon energy is shared between the ejected electron and nuclei in a correlated fashion, resulting in multiple diagonal lines in their joint energy spectrum governed by the energy conservation of all fragment particles.

Spectral encoding of x-ray/optical relative delay

We present a new technique for measuring the relative delay between a soft x-ray FEL pulse and an optical laser that indicates a sub 25 fs RMS measurement error. An ultra-short x-ray pulse photo-ionizes a semiconductor (Si(3)N(4)) membrane and changes the optical transmission. An optical continuum pulse with a temporally chirped bandwidth spanning 630 nm-710 nm interacts with the membrane such that the timing of the x-ray pulse can be determined from the onset of the spectral modulation of the transmitted optical pulse. This experiment demonstrates a nearly in situ single-shot measurement of the x-ray pulse arrival time relative to the ultra-short optical pulse.

A time-independent Hermitian Floquet approach for high harmonic generation in H(2)+ and HD+: effect of nonadiabatic interaction in HD+

We have theoretically investigated the high harmonic generation (HHG) in H(2)(+) and HD(+) using a time-independent Hermitian nonperturbative three-dimensional Floquet approach for continuous wave monochromatic lasers of intensities of 2.59x10(13), 4.0x10(13), and 5.6x10(13) W/cm(2), and wavelengths of 1064, 532, and 355 nm. For the moderate intensities and the wavelengths used, the Keldysh parameter gamma > 1 and no tunnel ionization occurs. We have endeavored to explain the dynamics of HHG in H(2)(+) and HD(+) within the framework of transitions due to electronic as well as intrinsic (for HD(+)) dipole moments and the nuclear motion on the field coupled ground and the first excited electronic states of these single-electron molecular ions, without considering any ionization. To evaluate the HHG spectra, the resonance Floquet quasienergy and the Fourier components of the Floquet state corresponding to the initial vibrational-rotational level v = 0, J = 0 have been calculated by solving the time-independent close-coupled Schrödinger equation following the Floquet ansatz. For HD(+), we have taken into account the symmetry breaking effect of the nonadiabatic interaction or coupling [beyond the Born-Oppenheimer (BO) dynamics] in our calculations of the HHG spectra and have compared the results with the HHG spectra of HD(+) obtained with the BO approximation.

Ion-energy dependence of asymmetric dissociation of D2 by a two-color laser field

Two-color (800 and 400 nm) short (45 fs) linearly polarized pulses are used to ionize and dissociate D2 into a neutral deuterium atom and a deuteron. The yields and energies of the ions are measured left and right along the polarization vector. As the relative phase of the two colors is varied, strong yield asymmetries are found in the ion-energy regions traditionally identified as bond softening, above-threshold dissociation and rescattering. The asymmetries in these regions are quite different. A model based on the dynamic coupling by the laser field of the gerade and ungerade states in the molecular ion accounts for many of the observed features.

Large-area imager of hydrogen leaks in fuel cells using laser-induced breakdown spectroscopy

We constructed a simple device, which utilized laser-induced breakdown spectroscopy to image H2 gas leaking from the surfaces of hydrogen fuel cells to ambient air. Nanosecond laser pulses of wavelength lambda=532 nm emitted from a neodymium-doped yttrium aluminum garnet laser were first compressed to a pulse length Deltat<1 ns using a stimulated Brillouin backscattering cell. Relay-imaging optics then focused this beam onto the H(2) leak and initiated the breakdown plasma. The Balmer-alpha (H-alpha) emission that emerged from this was collected with a 2-m-long macrolens assembly with a 90-mm-diameter image area, which covered a solid angle of approximately 1 x 10(-3)pi steradians seen from the plasma. The H-alpha light was isolated by two 100-mm-diameter interference filters with a 2 nm bandpass, and imaged by a thermoelectrically cooled charge-coupled device camera. By scanning the position of the laser focus, the spatial distribution of H2 gas over a 90-mm-diameter area was photographed with a spatial resolution of < or = 5 mm. Photoionization of the water vapor in the air caused a strong H-alpha background. By using pure N2 as a buffer gas, H2 leaks with rates of <1 cc/min were imaged. We also studied the possibilities of detecting He, Ne, or Xe gas leaks.

Master equation and control of an open quantum system with leakage

Given a multilevel system coupled to a bath, we use a Feshbach P, Q partitioning technique to derive an exact trace-nonpreserving master equation for a subspace S_{i} of the system. The resultant equation properly treats the leakage effect from S_{i} into the remainder of the system space. Focusing on a second-order approximation, we show that a one-dimensional master equation is sufficient to study problems of quantum state storage and is a good approximation, or exact, for several analytical models. It allows a natural definition of a leakage function and its control and provides a general approach to study and control decoherence and leakage. Numerical calculations on an harmonic oscillator coupled to a room temperature harmonic bath show that the leakage can be suppressed by the pulse control technique without requiring ideal pulses.

Angular tunneling ionization probability of fixed-in-space H2 molecules in intense laser pulses

We propose a new approach to obtain molecular frame photoelectron angular distributions from molecules ionized by intense laser pulses. With our method we study the angular tunnel ionization probability of H2 at a wavelength of 800 nm over an intensity range of 2-4.5 x 10(14) W/cm2. We find an anisotropy that is stronger than predicted by any existing model. To explain the observed anisotropy and its strong intensity dependence we develop an analytical model in the framework of the strong-field approximation. It expresses molecular ionization as a product of atomic ionization rate and a Fourier transform of the highest occupied molecular orbital filtered by the strong-field ionization process.

Enhancing high-order above-threshold dissociation of H2+ beams with few-cycle laser pulses

High-order (three-photon or more) above-threshold dissociation (ATD) of H(2)(+) has generally not been observed using 800 nm light. We demonstrate a strong enhancement of its probability using intense 7 fs laser pulses interacting with beams of H(2)(+), HD(+), and D(2)(+) ions. The mechanism invokes a dynamic control of the dissociation pathway. These measurements are supported by theory that additionally reveals, for the first time, an unexpectedly large contribution to ATD from highly excited electronic states.

Numerical simulation of nonadiabatic electron excitation in the strong-field regime. 3. Polyacene neutrals and cations

The electron optical response for a series of linear polyacenes and their molecular ions (mono and dications) in strong laser fields was studied using time-dependent Hartree-Fock theory. The interactions of benzene, naphthalene, anthracene, and tetracene with pulsed fields at a frequency of 1.55 eV and intensities of 8.77 x 10(13), 3.07 x 10(13), 1.23 x 10(13), and 2.75 x 10(12) W/cm2, respectively, were calculated using the 6-31G(d,p) basis set. Nonadiabatic processes, including nonadiabatic time evolution of the dipole moment, Löwden charges, and occupation numbers, were studied. The nonadiabatic response increased with the length of the molecule and was greatest for the molecular monocations. The only exception was tetracene, in which the very strong response of the dication was due to a near resonance with the applied field. The intensity and frequency dependence of the dipole moment response for the monocations of naphthalene and anthracene was also calculated. As the intensity increased, the population of higher-energy excited-states increased, and as the frequency increased, the excitation volume increased in good agreement with the Dykhne approximation.

Nonadiabatic chemical dynamics in an intense laser field: electronic wave packet coupled with classical nuclear motions

Dynamics of molecules in an intense laser field is studied in terms of the quantum electronic wave packet coupled with classical nuclear motions. The equations of motion are derived taking a proper account of molecular interactions with the vector potential of a classical electromagnetic field, along with the nonadiabatic interaction due to the breakdown of the Born-Oppenheimer approximation. With the aid of electronic structure calculations, the present method enables us to track, in an ab initio manner, the dynamics of polyatomic molecules in an intense field. Preliminary calculations are carried out for the vibrational state of LiF and a collision of Li+F under an intense laser pulse, which are limited to the domain of no ionization.

Explosive photodissociation of methane induced by ultrafast intense laser

A new type of molecular fragmentation induced by femtosecond intense laser at the intensity of 2 x 10(14) W/cm2 is reported. For the parent molecule of methane, ethylene, n-butane, and 1-butene, fluorescence from H (n = 3-->2), CH (A 2Delta, B 2Sigma-, and C 2Sigma+-->X 2Pi), or C2 (d 3Pi g-->a 3Pi u) is observed in the spectrum. It shows that the fragmentation is a universal property of neutral molecule in the intense laser field. Unlike breaking only one or two chemical bonds in conventional UV photodissociation, the fragmentation caused by the intense laser undergoes vigorous changes, breaking most of the bonds in the molecule, like an explosion. The fragments are neutral species and cannot be produced through Coulomb explosion of multiply charged ion. The laser power dependence of CH (A-->X) emission of methane on a log-log scale has a slope of 10 +/- 1. The fragmentation is thus explained as multiple channel dissociation of the superexcited state of parent molecule, which is created by multiphoton excitation.