A semi-classical model of attosecond electron localization in dissociative ionization of hydrogen

In the development of attosecond molecular science, a series of experiments have recently been performed where ionic fragment asymmetries in the dissociative ionization of H(2) into H(+) + H and that of D(2) into D(+) + D were used to uncover electron localization processes that occur on the attosecond and few-femtosecond timescale. Electron localization was observed both in strong-field dissociative ionization using carrier envelope phase-stable few-cycle laser pulses [Kling et al., Science, 2006, 312, 246] and in a two-color extreme ultra-violet + infrared attosecond pump-probe experiment [Sansone et al., Nature, 2010, 465, 763]. Here we show that the observed electron localization can be well understood using a semi-classical model that describes the dynamics in terms of quasi-static states that take the interaction of the molecule with the laser field instantaneously into account. The electron localization is shown to be determined by the passage of the dissociating molecule through a regime where the laser-molecule interaction is neither diabatic nor adiabatic.     

Molecular photoelectron interference effects by intense circularly polarized attosecond x-ray pulses

We present photoelectron interference effects in molecular photoionization by intense circularly polarized attosecond X-ray laser pulses. Simulations are performed on single electron molecular systems, H\(_{2}^{+}\), HeH²⁺, and H\(_{3}^{2+}\), by numerically solving the corresponding three-dimensional time dependent Schrödinger equations. Photoelectron spectra show that both momentum ring and stripe patterns are obtained, arising respectively from the interference of the direct ionized and scattered electron wave packets. Multi-center electron interference models in molecular ultrafast photoionization are used to describe the ionized electron dynamics. The interference patterns are shown to be dependent on molecular orbital symmetry, thus offering tools for attosecond imaging molecular structure at different molecular geometries.     

Correlated electron dynamics: how aromaticity can be controlled

In this Article, we show that the aromaticity of a molecule can be turned off by controlling the electron dynamics. We present a controlled switching from the aromatic ground state of benzene to two different nonaromatic states, using a laser pulse. The propagation of the molecular wave function is carried out with the time-dependent configuration interaction method. The laser pulse for switching between the ground and excited states is optimized using optimal control theory. Bond orders and Mulliken charges serve as an aromaticity criterion. The nonaromatic target states exhibit localized bonds and partial charges on the carbon atoms; these localized electrons circulate on an attosecond time scale in the ring system.     

Quantum-classical correspondence of a field induced KAM-type transition: A QTM approach

A transition from regular to chaotic behaviour in the dynamics of a classical Henon-Heiles oscillator in the presence of an external field is shown to have a similar quantum signature when studied using the pertaining phase portraits and the associated Kolmogorov-Sinai-Lyapunov entropies obtained through the corresponding Bohmian trajectories.     

Non-adiabatic molecular dynamics with quantum solvent effects

Three novel approaches extending quantum-classical non-adiabatic (NA) molecular dynamics (MD) to include quantum effects of solvent environments are described. In a standard NA-MD the solute subsystem is treated quantum mechanically, while the larger solvent part of a system is treated classically. The three novel approaches presented here are based on the Bohmian formulation of quantum mechanics, the stochastic Schrödinger equation for the evolution of open quantum systems and the quantized Hamilton dynamics generalization of classical mechanics. The approaches extend the standard NA-MD to incorporate the following quantum effects of the solvent. (1) Branching, i.e. the ability of solvent quantum wave packets to split and follow asymptotically diverging trajectories correlated with different quantum states of the solute. (2) Decoherence, i.e. loss of quantum interference within the solute subsystem induced by the diverging solvent trajectories. (3) Zero point energy that contributes to NA coupling and must be preserved during the energy exchange between solvent and solute degrees of freedom. The Bohmian quantum-classical mechanics, stochastic mean-field and quantized mean-field approximations incorporate the quantum solvent effects into the standard quantum-classical NA-MD in a straightforward and efficient way that can be easily applied to quantum dynamics of condensed phase chemical systems.     

Stereocontrol of attosecond time-scale electron dynamics in ABCU using ultrafast laser pulses: a computational study

The attosecond time-scale electronic dynamics induced by an ultrashort laser pulse is computed using a multi configuration time dependent approach in ABCU (C(10)H(19)N), a medium size polyatomic molecule with a rigid cage geometry. The coupling between the electronic states induced by the strong pulse is included in the many electron Hamiltonian used to compute the electron dynamics. We show that it is possible to implement control of the electron density stereodynamics in this medium size molecule by varying the characteristics of the laser pulse, for example by polarizing the electric field either along the N-C axis of the cage, or in the plane perpendicular to it. The excitation produces an oscillatory, non-stationary, electronic state that exhibits localization of the electron density in different parts of the molecule both during and after the pulse. The coherent oscillations of the non-stationary electronic state are also demonstrated through the alternation of the dipole moment of the molecule.     

Quantum dynamics with Bohmian trajectories

We describe the advantages and disadvantages of numerical methods when Bohmian trajectory grids are used for numerical simulations of quantum dynamics. We focus on the crucial noncrossing property of Bohmian trajectories, which, numerically, must be given careful attention. Failure to do so causes instabilities or leads to false simulations.     

强阿秒XUV激走脉冲下H2^＋分子光电子能谱角分布非对称性

利用2D平面模型,求解了描述定向H2^＋分子和阿秒XUV脉冲相互作用的薛定谔方程,并求得光电子的角度分布．在计算模型中,采用基态1sσg和第一激发态2pσu的等比例混合态作为初始态,而激光脉冲的光子能量大于电离势,强度为10^14 W/cm^2． 计算结果表明,光电子角分布的非对称性和脉冲的宽度密切相关．这种非对称性实际上是由于初始态的基态和激发态的相干振荡而导致的．当使用长脉冲时,这种相干振荡的周期效应就会被平均而消失,从而产生的光电子能谱会呈对称角分布．     

Quantum analogue of the Kolmogorov–Arnold–Moser transition in different quantum anharmonic oscillators

The results of numerical computations are presented for the Bohmian trajectories of the family of different one‐ and two‐dimensional anharmonic oscillators, which exhibit regular or chaotic motion in both classical and quantum domains, depending on the values of the parameters appearing in the respective Hamiltonians. Quantum signatures of the Kolmogorov–Arnold–Moser (KAM) transition from the regular to chaotic classical dynamics of these oscillators are studied using a quantum theory of motion (QTM) as developed by de Broglie and Bohm. A phase space distance function between two initially close Bohmian trajectories, the associated Kolmogorov–Sinai–Lyapunov (KSL) entropy, the phase space volume, the autocorrelation function, the associated power spectrum, and the nearest‐neighbor spacing distribution, clearly differentiate the quantum analogues of the corresponding regular and chaotic motions in the classical domain. These quantum anharmonic oscillators are known to be useful in several diverse branches of science. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Ab initio molecular dynamics study on the electron capture processes of protonated methane (CH5+)

Electron capture dynamics of protonated methane (CH5(+)) have been investigated by means of a direct ab initio molecular dynamics (MD) method. First, the ground and two low-lying state structures of CH5 (+) with eclipsed Cs , staggered Cs and C2v symmetries were examined as initial geometries in the dynamics calculation. Next, the initial structures of CH5 (+) in the Franck-Condon (FC) region were generated by inclusion of zero point energy and then trajectories were run from the selected points on the assumption of vertical electron capture. Two competing reaction channels were observed: CH5 (+) + e (-)--> CH4 + H (I) and CH5 (+) + e (-) --> CH3 + H2 (II). Channel II occurred only from structures very close to the s- Cs geometry for which two protons with longer C-H distances are electronically equivalent in CH5 (+). These protons have the highest spin density as hydrogen atoms following vertical electron capture of CH5 (+) and are lost as H2. On the other hand, channel I was formed from a wide structural region of CH5 (+). The mechanism of the electron capture dynamics of CH5 is discussed on the basis of the theoretical results.     

Fragment momentum distributions obtained from coupled electron-nuclear dynamics

We theoretically investigate fragmentation processes induced by femtosecond laser pulses within a model which incorporates electronic and nuclear motion. Single-pulse excitation leads to diffraction patterns in the electron momentum distribution which depend on the nature of the electronic state and also on the nuclear charge distribution. Additional structures appear in the nuclear momentum distribution if two time-delayed pulses produce fragments in the same dissociation channel. It is shown that these functions are modified by the electronic degree-of-freedom. A simultaneous excitation of two different electronic states results in further interferences which are related to electronic wave-packet dynamics on the attosecond time-scale.     

Quantum trajectories in elastic atom-surface scattering: threshold and selective adsorption resonances

The elastic resonant scattering of He atoms off the Cu(117) surface is fully described with the formalism of quantum trajectories provided by Bohmian mechanics. Within this theory of quantum motion, the concept of trapping is widely studied and discussed. Classically, atoms undergo impulsive collisions with the surface, and then the trapped motion takes place covering at least two consecutive unit cells. However, from a Bohmian viewpoint, atom trajectories can smoothly adjust to the equipotential energy surface profile in a sort of sliding motion; thus the trapping process could eventually occur within one single unit cell. In particular, both threshold and selective adsorption resonances are explained by means of this quantum trapping considering different space and time scales. Furthermore, a mapping between each region of the (initial) incoming plane wave and the different parts of the diffraction and resonance patterns can be easily established, an important issue only provided by a quantum trajectory formalism.     

Calculating fifth-order Raman signals for various molecular liquids by equilibrium and nonequilibrium hybrid molecular dynamics simulation algorithms

The fifth-order two-dimensional (2D) Raman signals have been calculated from the equilibrium and nonequilibrium (finite field) molecular dynamics simulations. The equilibrium method evaluates response functions with equilibrium trajectories, while the nonequilibrium method calculates a molecular polarizability from nonequilibrium trajectories for different pulse configurations and sequences. In this paper, we introduce an efficient algorithm which hybridizes the existing two methods to avoid the time-consuming calculations of the stability matrices which are inherent in the equilibrium method. Using nonequilibrium trajectories for a single laser excitation, we are able to dramatically simplify the sampling process. With this approach, the 2D Raman signals for liquid xenon, carbon disulfide, water, acetonitrile, and formamide are calculated and discussed. Intensities of 2D Raman signals are also estimated and the peak strength of formamide is found to be only five times smaller than that of carbon disulfide.     

Bohmian mechanics with complex action: a new trajectory-based formulation of quantum mechanics

In recent years there has been a resurgence of interest in Bohmian mechanics as a numerical tool because of its local dynamics, which suggest the possibility of significant computational advantages for the simulation of large quantum systems. However, closer inspection of the Bohmian formulation reveals that the nonlocality of quantum mechanics has not disappeared-it has simply been swept under the rug into the quantum force. In this paper we present a new formulation of Bohmian mechanics in which the quantum action, S, is taken to be complex. This leads to a single equation for complex S, and ultimately complex x and p but there is a reward for this complexification-a significantly higher degree of localization. The quantum force in the new approach vanishes for Gaussian wave packet dynamics, and its effect on barrier tunneling processes is orders of magnitude lower than that of the classical force. In fact, the current method is shown to be a rigorous extension of generalized Gaussian wave packet dynamics to give exact quantum mechanics. We demonstrate tunneling probabilities that are in virtually perfect agreement with the exact quantum mechanics down to 10(-7) calculated from strictly localized quantum trajectories that do not communicate with their neighbors. The new formulation may have significant implications for fundamental quantum mechanics, ranging from the interpretation of non-locality to measures of quantum complexity.     

Hyperfine coupling tensors of the benzosemiquinone radical anion from Car-Parrinello molecular dynamics.

Based on Car-Parrinello ab initio molecular dynamics simulations of the benzosemiquinone radical anion in both aqueous solution and the gas phase, density functional calculations provide the currently most refined EPR hyperfine coupling (HFC) tensors of semiquinone nuclei and solvent protons. For snapshots taken at regular intervals from the molecular dynamics trajectories, cluster models with different criteria for inclusion of water molecules and an additional continuum solvent model are used to analyse the HFCs. These models provide a detailed picture of the effects of dynamics and of different intermolecular interactions on the spin-density distribution and HFC tensors. Comparison with static calculations allows an assessment of the importance of dynamical effects, and of error compensation in static DFT calculations. Solvent proton HFCs depend characteristically on the position relative to the semiquinone radical anion. A point-dipolar model works well for in-plane hydrogen-bonded protons but deviates from the quantum chemical values for out-of-plane hydrogen bonding.     

Unified derivation of Bohmian methods and the incorporation of interference effects

We present a unified derivation of Bohmian methods that serves as a common starting point for the derivative propagation method (DPM), Bohmian mechanics with complex action (BOMCA), and the zero-velocity complex action method (ZEVCA). The unified derivation begins with the ansatz psi = eiS/Planck's where the action (S) is taken to be complex, and the quantum force is obtained by writing a hierarchy of equations of motion for the phase partial derivatives. We demonstrate how different choices of the trajectory velocity field yield different formulations such as DPM, BOMCA, and ZEVCA. The new derivation is used for two purposes. First, it serves as a common basis for comparing the role of the quantum force in the DPM and BOMCA formulations. Second, we use the new derivation to show that superposing the contributions of real, crossing trajectories yields a nodal pattern essentially identical to that of the exact quantum wavefunction. The latter result suggests a promising new approach to deal with the challenging problem of nodes in Bohmian mechanics.     

Electron trajectories in molecular orbitals

The time-dependent Schrödinger equation can be rewritten so that its interpretation is no longer probabilistic. Two well-known and related reformulations are Bohmian mechanics and quantum hydrodynamics. In these formulations, quantum particles follow real, deterministic trajectories influenced by a quantum force. Generally, trajectory methods are not applied to electronic structure calculations as they predict that the electrons in a ground-state, real, molecular wavefunction are motionless. However, a spin-dependent momentum can be recovered from the nonrelativistic limit of the Dirac equation. Therefore, we developed new, spin-dependent equations of motion for the quantum hydrodynamics of electrons in molecular orbitals. The equations are based on a Lagrange multiplier, which constrains each electron to an isosurface of its molecular orbital, as required by the spin-dependent momentum. Both the momentum and the Lagrange multiplier provide a unique perspective on the properties of electrons in molecules.     

Magnetic relaxation dispersion of lithium ion in solutions of DNA

The magnetic field dependence of the nuclear spin-lattice relaxation rate constant defines the magnetic relaxation dispersion (MRD) and provides a direct characterization of the molecular dynamics that cause fluctuations in the magnetic couplings in the system and may also indicate the dimensional constraints on the motion. The counterion cloud surrounding a linear polyelectrolyte ion, such as DNA in solution, provides an interesting opportunity for ion confinement that helps in understanding the thermodynamics and the dynamics of the interactions between the polyion and other solutes. The MRD profiles of lithium ion and tetramethylammonium ion were recorded in dilute aqueous solutions of native calf thymus DNA, which provides a long, charged rod that reorients slowly. The 7Li ion relaxes through the nuclear electric quadrupole coupling and the proton-lithium dipole-dipole coupling; the protons of the tetramethylammonium ion relax by dipole-dipole coupling. MRD profiles of the 7Li+ ion are dominated by transient interactions with the DNA that yield a linear dependence of the spin-lattice relaxation rate constant on the logarithm of the Larmor frequency. This magnetic field dependence is consistent with diffusive ion motions that modulate two spatial coordinates that characterize the relaxation couplings in the vicinity of the polyion. Motions around the rod and fluctuations in the ion distance from the rod are consistent with these constraints for lithium. The magnetic field dependence of the tetramethylammonium ion proton relaxation rate constant is weak, but also approximately a linear function of the logarithm of the Larmor frequency, which implies that the field dependence is caused in part by local order in the DNA solution.