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.
An adaptive pseudospectral method for wave packet dynamics
Rent:
Rent this article for
USD
10.1063/1.4737893
/content/aip/journal/jcp/137/4/10.1063/1.4737893
http://aip.metastore.ingenta.com/content/aip/journal/jcp/137/4/10.1063/1.4737893
View: Figures

Figures

Image of FIG. 1.
FIG. 1.

Error as function of the time step for Galerkin (crosses) and collocation (plain lines). Errors compared to a reference solution within the Hagedorn space are shown with dashed lines, and errors compared to a spatially well-resolved reference are shown with solid lines. The TDSE is integrated up to t = 1 with 64 basis functions in (a), and to t = 5 with 8 basis functions in (b).

Image of FIG. 2.
FIG. 2.

Projection of to two different Hagedorn spaces, both with K = 32 basis functions. The figures show u(x) with solid lines and its projection with dashed lines. In both cases, ɛ = 1, q = p = 0, P = i/Q. In (a), Q = 1 and in (b), Q = 0.5. The ℓ-errors are 0.93 and 1.8 × 10−5, respectively. In (b), the two lines are indistinguishable.

Image of FIG. 3.
FIG. 3.

(a) The time evolution of the uncertainty parameter QQ*. The solid line indicates control-based propagation and the dotted semiclassical propagation. The dashed line indicates the control signal. (b) A lin-log plot of the squared modulus of the wave function at time t = 5. The crosses show the location of the collocation points in the control scheme.

Image of FIG. 4.
FIG. 4.

The PES of the IBr molecule with sketched wave packets and grids. The equidistant Fourier grid of the state stays fixed, while the Gauss-Hermite grids in the two excited states follow the wave packets.

Image of FIG. 5.
FIG. 5.

The real part of ψ2 at time t 0, when the laser pulse peaks. The cross marks the coordinate of strongest coupling with the ground electronic state. The circle marks where the classical coordinate q would have been located if no control were used.

Image of FIG. 6.
FIG. 6.

The solution, real parts and moduli, at the three potential energy surfaces at time .

Image of FIG. 7.
FIG. 7.

(a) The front of ψ3 in modulus at time . The dashed line indicates , and the cross x δ. As grows, the dashed line will rise above the top of the hump, and x δ will encounter a discontinuity. The response to this is shown in (b). The dashed line shows the control signal for |Q|, and the solid line its actual value.

Image of FIG. 8.
FIG. 8.

Level curves for the excited state PES of CO2. The location of the saddle point is indicated by a cross, and the centrepoint of the initial Gaussian wave packet by a ring. The level curves are separated by 0.015 hartree.

Image of FIG. 9.
FIG. 9.

Time evolution of the wave function. The modulus of the wave function, |Ψ(s, a)|, is plotted with a separation of 1 between the level curves. The cross indicates the location of the saddle point in the potential. The wave function is plotted at times t = (a) , (b) , (c) , and (d) .

Image of FIG. 10.
FIG. 10.

(a) The modulus of Q 2, which determines the spread of the collocation points along the asymmetric stretch coordinate axis, plotted against time. Semiclassical propagation of Q and P is shown with solid lines, and control-based propagation with dashed lines. (b) The time evolution of the maximum pointwise error for the two propagation schemes.

Image of FIG. 11.
FIG. 11.

The performance of the controller. |Q 1| (below) and |Q 2| (above) are plotted against time. The square root of the control signal is shown with solid lines, and the actual parameter values with dashed lines.

Image of FIG. 12.
FIG. 12.

Level curves for the pointwise error at time for solutions calculated with (a) a control-propagated Hagedorn basis with 322 basis functions, and (b) a Fourier basis with 522 basis functions. The curve separation is 1.0 × 10−5, and the maximum pointwise error is 8.7 × 10−5 and 7.2 × 10−5 for the Hagedorn and Fourier solutions, respectively.

Loading

Article metrics loading...

/content/aip/journal/jcp/137/4/10.1063/1.4737893
2012-07-26
2014-04-16
Loading

Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: An adaptive pseudospectral method for wave packet dynamics
http://aip.metastore.ingenta.com/content/aip/journal/jcp/137/4/10.1063/1.4737893
10.1063/1.4737893
SEARCH_EXPAND_ITEM