1887
banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Multimode quantum dynamics using Gaussian wavepackets: The Gaussian-based multiconfiguration time-dependent Hartree (G-MCTDH) method applied to the absorption spectrum of pyrazine
Rent:
Rent this article for
USD
10.1063/1.2996349
/content/aip/journal/jcp/129/17/10.1063/1.2996349
http://aip.metastore.ingenta.com/content/aip/journal/jcp/129/17/10.1063/1.2996349

Figures

Image of FIG. 1.
FIG. 1.

For the four-mode model (Sec. IV B), the figure illustrates the convergence properties of the all-Gaussian (vMCG) method and the comparison to classically evolving GWP particles. (a) vMCG calculations with 40 (dashed trace), 80 (full thin trace), and 200 (full bold trace) 4D GWPs (i.e., 20, 40, and 100 GWPs per state). (b) For comparison, analogous calculations with 200 4D (dashed trace), 3368 1D (full thin trace), and 1D (full bold trace) classically evolving GWPs are shown.

Image of FIG. 2.
FIG. 2.

Autocorrelation functions of the 10-mode model of Sec. IV B (black trace) as compared with the full 24-mode model of Sec. IV C [gray (red) dashed trace] using standard MCTDH calculations. The structure of the first two recurrences is reproduced by the reduced-dimensional ten-mode model.

Image of FIG. 3.
FIG. 3.

Convergence properties of the -mode G-MCTDH hybrid calculations for the 10-mode model of Sec. IV B. (a) The results of two G-MCTDH calculations are compared with the standard MCTDH (exact) result. The first calculation (dashed line), i.e., G-MCTDH I of Table II, employs the same number of particles as compared with the standard MCTDH calculation (thin solid line). The second calculation (bold solid line) has an augmented GWP basis in the secondary subspace (i.e., G-MCTDH II of Table II) and can be considered converged. (b) For comparison, a vMCG calculation (dashed line) with 20 GWPs per particle is shown (see vMCG of Table II). The G-MCTDH II calculation (bold solid line) and reference MCTDH calculation (thin solid line) are also reproduced. The vMCG calculation is significantly harder to converge than the G-MCTDH hybrid calculations.

Image of FIG. 4.
FIG. 4.

(a) Autocorrelation function and (b) spectrum for the nearly converged G-MCTDH calculation (black bold solid line), i.e., G-MCTDH II of Table II, as compared with the standard MCTDH result [gray (red) dashed line]. The spectra are convoluted with a damping function with decay constant .

Image of FIG. 5.
FIG. 5.

Spectrum for a small G-MCTDH calculation [gray (red) dashed trace], i.e., G-MCTDH III of Table II, as compared with the standard MCTDH result (black trace). This illustrates that good results can be obtained even with a cheap calculation.

Image of FIG. 6.
FIG. 6.

Autocorrelation function for the G-MCTDH 24-mode calculation (bold black trace), i.e., G-MCTDH of Table IV, as compared with a standard MCTDH calculation [gray (red) dashed trace], i.e., MCTDH of Table IV Panels (a) and (b) show different time intervals, i.e., 150 fs vs 300 fs. Both calculations are near convergence.

Image of FIG. 7.
FIG. 7.

Spectra obtained by Fourier transformation of the autocorrelation functions of Fig. 6, i.e., G-MCTDH 24 mode of Table IV [gray (red) dashed trace] vs MCTDH 24 mode of Table IV (black trace). The spectra are convoluted with a damping function with a decay constant of .

Image of FIG. 8.
FIG. 8.

Comparison of the G-MCTDH 24-mode calculation (solid line) with the experimental spectrum of Ref. 19 (dotted line). In panel (a), the (high-frequency side) and (low-frequency side) regions are shown, while panel (b) is centered on the region. Note that intensities in the region are not correctly reproduced since the initial condition used in the wavepacket propagation is restricted to one of the diabatic states.

Image of FIG. 9.
FIG. 9.

Diabatic populations obtained for the G-MCTDH 24-mode calculation of Table IV (black trace) as compared with a standard MCTDH calculation of Table IV [gray (red) dashed trace].

Tables

Generic image for table
Table I.

Computational resources required by the four-mode vMCG calculations. is the number of 4D GWPs used to describe the wavepacket in the two electronic states.

Generic image for table
Table II.

MCTDH vs G-MCTDH basis sets for ten-mode calculations. The round brackets denote the combination of vibrational modes and the square brackets the number of SPFs or GWPs used for the representation of the wave function in the and states. The number of modes in one combination defines the dimensionality of the corresponding SPFs. These SPFs are represented on a grid whose size is given by the product of the number of grid points used for each mode of the corresponding combination.

Generic image for table
Table III.

Computational resources required by the ten-mode calculations detailed in Table II. “Num (A-vec)” is the number of expansion coefficients, and “Num (SPFs/GWPs)” is the number of grid points and GWP parameters. Memory and CPU are the computer resources on an Intel Pentium Quad Core Q6600 (2.4 GHz, 8 Gbytes) machine.

Generic image for table
Table IV.

MCTDH vs G-MCTDH basis sets for 24-mode calculation. The conventions are as in Table II.

Generic image for table
Table V.

Computational resources required by the 24-mode calculations. The columns headed “ (MCTDH)” and “ (G-MCTDH)” are for the calculations detailed in Table IV. Conventions are chosen as in Table III.

Loading

Article metrics loading...

/content/aip/journal/jcp/129/17/10.1063/1.2996349
2008-11-05
2014-04-23
Loading

Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Multimode quantum dynamics using Gaussian wavepackets: The Gaussian-based multiconfiguration time-dependent Hartree (G-MCTDH) method applied to the absorption spectrum of pyrazine
http://aip.metastore.ingenta.com/content/aip/journal/jcp/129/17/10.1063/1.2996349
10.1063/1.2996349
SEARCH_EXPAND_ITEM