Gauge-invariant expectation values of the energy of a molecule in an electromagnetic field

Origin Independent Current Density Vector Fields Induced by Time-Dependent Magnetic Field. I. The LiH molecule

The electronic current density, induced in a molecule by the optical magnetic field associated to a frequency dependent monochromatic plane wave, assumed to be spatially uniform within the electric quadrupole approximation, has been studied by a theoretical method based on continuous translation of its origin. The induced electronic current density vector field designated by this procedure, invariant of the origin for any point of the molecular domain, is obtained via a computational scheme formally annihilating the diamagnetic contribution of the conventional common-origin approach based on perturbation theory. In a preliminary application of the theoretical methods outlined in the present work, the simple molecule of lithium hydride has been investigated. Particular attention has been paid to the structure of induced electronic current density for several values of the magnetic field frequency, by investigating equilibrium points of four different types, organized in stagnation lines which constitute its stagnation graph, i.e., a topological fingerprint of the vector field conveying complete information in condensed form, to verify the fulfillment of fundamental requirements, e.g., the continuity equation and the Poincaré-Hopf theorem on spherical and toroidal surfaces.

Quantum transition probabilities due to overlapping electromagnetic pulses: Persistent differences between Dirac's form and nonadiabatic perturbation theory

The probability of transition to an excited state of a quantum system in a time-dependent electromagnetic field determines the energy uptake from the field. The standard expression for the transition probability has been given by Dirac. Landau and Lifshitz suggested, instead, that the adiabatic effects of a perturbation should be excluded from the transition probability, leaving an expression in terms of the nonadiabatic response. In our previous work, we have found that these two approaches yield different results while a perturbing field is acting on the system. Here, we prove, for the first time, that differences between the two approaches may persist after the perturbing fields have been completely turned off. We have designed a pair of overlapping pulses in order to establish the possibility of lasting differences, in a case with dephasing. Our work goes beyond the analysis presented by Landau and Lifshitz, since they considered only linear response and required that a constant perturbation must remain as t → ∞. First, a "plateau" pulse populates an excited rotational state and produces coherences between the ground and excited states. Then, an infrared pulse acts while the electric field of the first pulse is constant, but after dephasing has occurred. The nonadiabatic perturbation theory permits dephasing, but dephasing of the perturbed part of the wave function cannot occur within Dirac's method. When the frequencies in both pulses are on resonance, the lasting differences in the calculated transition probabilities may exceed 35%. The predicted differences are larger for off-resonant perturbations.

Quantum transition probabilities during a perturbing pulse: Differences between the nonadiabatic results and Fermi's golden rule forms

For a perturbed quantum system initially in the ground state, the coefficient c_k(t) of excited state k in the time-dependent wave function separates into adiabatic and nonadiabatic terms. The adiabatic term a_k(t) accounts for the adjustment of the original ground state to form the new ground state of the instantaneous Hamiltonian H(t), by incorporating excited states of the unperturbed Hamiltonian H₀ without transitions; a_k(t) follows the adiabatic theorem of Born and Fock. The nonadiabatic term b_k(t) describes excitation into another quantum state k; b_k(t) is obtained as an integral containing the time derivative of the perturbation. The true transition probability is given by b_k(t) ², as first stated by Landau and Lifshitz. In this work, we contrast b_k(t) ² and c_k(t) ². The latter is the norm-square of the entire excited-state coefficient which is used for the transition probability within Fermi's golden rule. Calculations are performed for a perturbing pulse consisting of a cosine or sine wave in a Gaussian envelope. When the transition frequency ω_k0 is on resonance with the frequency ω of the cosine wave, b_k(t) ² and c_k(t) ² rise almost monotonically to the same final value; the two are intertwined, but they are out of phase with each other. Off resonance (when ω_k0 ≠ ω), b_k(t) ² and c_k(t) ² differ significantly during the pulse. They oscillate out of phase and reach different maxima but then fall off to equal final values after the pulse has ended, when a_k(t) ≡ 0. If ω_k0 < ω, b_k(t) ² generally exceeds c_k(t) ², while the opposite is true when ω_k0 > ω. While the transition probability is rising, the midpoints between successive maxima and minima fit Gaussian functions of the form a exp[-b(t - d)²]. To our knowledge, this is the first analysis of nonadiabatic transition probabilities during a perturbing pulse.