Nonadiabatic electron wavepacket dynamics of molecules in an intense optical field: an ab initio electronic state study

A theory of quantum electron wavepacket dynamics that nonadiabatically couples with classical nuclear motions in intense optical fields is studied. The formalism is intended to track the laser-driven electron wavepackets in terms of the linear combination of configuration-state functions generated with ab initio molecular orbitals. Beginning with the total quantum Hamiltonian for electrons and nuclei in the vector potential of classical electromagnetic field, we reduce the Hamiltonian into a mixed quantum-classical representation by replacing the quantum nuclear momentum operators with the classical counterparts. This framework gives equations of motion for electron wavepackets in an intense laser field through the time dependent variational principle. On the other hand, a generalization of the Newtonian equations provides a matrix form of forces acting on the nuclei for nonadiabatic dynamics. A mean-field approximation to the force matrix reduces this higher order formalism to the semiclassical Ehrenfest theory in intense optical fields. To bring these theories into a practical quantum chemical package for general molecules, we have implemented the relevant ab initio algorithms in it. Some numerical results in the level of the semiclassical Ehrenfest-type theory with explicit use of the nuclear kinematic (derivative) coupling and the velocity form for the optical interaction are presented.     

Formation of ultracold SrYb molecules in an optical lattice by photoassociation spectroscopy: theoretical prospects

State-of-the-art ab initio techniques have been applied to compute the potential energy curves for the SrYb molecule in the Born-Oppenheimer approximation for the electronic ground state and the first fifteen excited singlet and triplet states. All the excited state potential energy curves were computed using the equation of motion approach within the coupled-cluster singles and doubles framework and large basis-sets, while the ground state potential was computed using the coupled cluster method with single, double, and noniterative triple excitations. The leading long-range coefficients describing the dispersion interactions at large interatomic distances are also reported. The electric transition dipole moments have been obtained as the first residue of the polarization propagator computed with the linear response coupled-cluster method restricted to single and double excitations. Spin-orbit coupling matrix elements have been evaluated using the multireference configuration interaction method restricted to single and double excitations with a large active space. The electronic structure data were employed to investigate the possibility of forming deeply bound ultracold SrYb molecules in an optical lattice in a photoassociation experiment using continuous-wave lasers. Photoassociation near the intercombination line transition of atomic strontium into the vibrational levels of the strongly spin-orbit mixed b(3)Σ(+), a(3)Π, A(1)Π, and C(1)Π states with subsequent efficient stabilization into the v' = 1 vibrational level of the electronic ground state is proposed. Ground state SrYb molecules can be accumulated by making use of collisional decay from v' = 1 to v' = 0. Alternatively, photoassociation and stabilization to v' = 0 can proceed via stimulated Raman adiabatic passage provided that the trapping frequency of the optical lattice is large enough and phase coherence between the pulses can be maintained over at least tens of microseconds.     

Inelastic collisions of ultracold polar LiCs molecules with caesium atoms in an optical dipole trap

We investigate collisions of ultracold polar LiCs molecules and ultracold caesium atoms. LiCs molecules are formed in an optical dipole trap by photoassociation of caesium and lithium atoms via the B(1)Π excited state followed by spontaneous emission to the X(1)Σ(+) ground state and the lowest triplet state a(3)Σ(+). The molecules are then stored together with caesium atoms in the same optical trap. Rate coefficients for the loss of molecules induced by collisions with surrounding Cs atoms are measured for molecular ensembles produced via different photoassociation resonances. The results are analyzed in terms of the unitarity limit for the inelastic rates and predictions from the universal model of Idziaszek and Julienne (Phys. Rev. Lett., 2010, 104, 113202).     

Nonresonant optical activity of isolated organic molecules

Cavity ring-down polarimetry (CRDP) has been exploited to interrogate the nonresonant optical activity (or circular birefringence) of prototypical organic compounds in the vapor phase, thereby revealing the intrinsic chiro-optical response evoked from isolated (solvent-free) molecules. Specific polarization rotation parameters have been measured at two distinct excitation wavelengths (355 nm and 633 nm) for a variety of gas-phase species drawn from the terpene, epoxide, and alkane/alkene families, with complementary solution-phase polarimetric studies serving to highlight the pronounced influence of solute-solvent interactions. Time-dependent linear response calculations performed at high levels of density functional theory have been enlisted to unravel the structural and electronic origins for observed behavior. Aside from elucidating the complex solvation processes that mediate chiro-optical phenomena taking place in condensed media, this study affords a critical assessment for emerging ab initio predictions of nonresonant optical activity and for their promising ability to assist in the determination of absolute molecular stereochemistry.     

Toward room-temperature optical manipulation of small molecules

Room-temperature optical manipulation of small molecules is a challenging issue in the field of material science. To increase optical force for a single molecule trapping, it has been recognized that resonant excitation of molecules should be controlled under the light illumination. Strongly interacting molecules with solid surfaces at electrified interfaces show the exotic behavior of electronic excitation by localized surface plasmon. In this review, we emphases that surface-enhanced Raman scattering can be used to evaluate the resonant excitation of target molecules at interfaces. Under such excitation, the diffusion of small molecules can be controlled by the optical force generated by the intensity gradient of a highly localized electric field.     

Some features of coherent optical effects in large molecules

In this paper we present a theoretical study of coherent optical effects such as optical free induction decay (OFID) and photon echoes (PE) from an ensemble of collision-free large molecules. We have introduced a generalized effective hamiltonian to account for the temporal characteristics of the macroscopic polarization. Explicit results have been obtained for the OFID and PE under short-time on-resonance optical pulsed excitation elucidating the nature of phase destruction effects originating from intramolecular interstate and intrastate coupling and from inhomogeneous broadening effects.     

Vibrational lifetime measurements of matrix isolated molecules by an optical “double resonance” method

An optical-optical “double-resonance” experiment has been developed and successfully applied to the study of vibrational lifetimes. This pulsed dual laser technique should be applicable to all species possessing resolvable allowed vibronic transitions and it permits the direct measurement of vibrational lifetimes which, at present, may be as short as 20 μs. Application of this technique to the υ″ = 1 and υ″ = 2 states of matrix isolated C₂^? yields half lifes of 0.2 and 1.2 ms, respectively, in an argon matrix at 16K, and 0.3 and 1.3 ms in an N₂ matrix at 14K. No significant temperature variation of these rate constants has been found in the range 14 to 24 K. This is the first direct measurement of the vibrational lifetime of a homonuclear diatomic molecule isolated in a matrix environment.     

Calculation of the nonlinear optical properties of molecules

A finite-field method for the calculation of polarizabilities and hyperpolarizabilities is developed based on both an energy expansion and a dipole moment expansion. This procedure is implemented in the MOPAC semiempirical program. Values and components of the dipole moment (μ), polarizability (α), first hyperpolarizability (β), and second hyperpolarizability (γ) are calculated as an extension of the usual MOPAC run. Applications to benzene and substituted benzenes are shown as test cases utilizing both MNDO and AM1 Hamiltonians.     

Communication: Tailoring the optical gap in light-harvesting molecules

Systematically varying the optical gap that is associated with charge-transfer excitations is an important step in the design of light-harvesting molecules. So far the guidance that time-dependent density functional theory could give in this process was limited by the traditional functionals' inability to describe charge-transfer excitations. We show that a nonempirical range-separated hybrid approach allows to reliably predict charge-transfer excitations for molecules of practically relevant complexity. Calculated absorption energies agree with measured ones. We predict from theory that by varying the number of thiophenes in donor-acceptor-donor molecules, the energy of the lowest optical absorption can be tuned to the lower end of the visible spectrum. Saturation sets in at about five thiophene rings.     

Graphical prediction of quantum interference-induced transmission nodes in functionalized organic molecules

Quantum interference (QI) in molecular transport junctions can lead to dramatic reductions of the electron transmission at certain energies. In a recent work [Markussen et al., Nano Lett., 2010, 10, 4260] we showed how the presence of such transmission nodes near the Fermi energy can be predicted solely from the structure of a conjugated molecule when the energies of the atomic p(z) orbitals do not vary too much. Here we relax the assumption of equal on-site energies and generalize the graphical scheme to molecules containing different atomic species. We use this diagrammatic scheme together with tight-binding and density functional theory calculations to investigate QI in linear molecular chains and aromatic molecules with different side groups. For the molecular chains we find a linear relation between the position of the transmission nodes and the side group π orbital energy. In contrast, the transmission functions of functionalized aromatic molecules generally display a rather complex nodal structure due to the interplay between molecular topology and the energy of the side group orbital.     

Characterization of intensity-dependent optical rotation phenomena in chiral molecules in solution

The rotation of the plane of polarization of linearly polarized light by chiral molecules in solution is due to a forward scattering event. Ordinary optical rotation, a single-photon effect, is independent of intensity. As the light intensity is increased, other effects can appear, such as two-photon scattering or alignment of the molecule by one photon and scattering with a change of polarization by another. Both of these effects result in intensity-dependent (or nonlinear) optical rotation. A polarimeter was used to measure the nonlinear optical rotation of solutions in a heterodyne experiment. No nonlinear optical rotation was found in molecules lacking an absorption band near the laser frequency. In the three pyrimidine nucleosides studied, which do have such an absorption band, a nonlinear optical rotation was identified that was cumulative with each laser pulse. The effect persisted with a time constant that was on the order of seconds and characteristic of the molecule.     

Influence of isomerization on nonlinear optical properties of molecules

The influence of rotational and geometrical isomerism on the nonlinear optical (NLO) properties, specifically the first-order hyperpolarizability beta, of chromophores of current interest has been investigated with density functional theory (DFT). In the first of this two-part study, the rotational isomerism of a linear chromophore was explored. Calculation of the torsion potentials about two of the rotatable and conformation-changing single bonds in a chromophore demonstrated the near equality of the molecular energies at 0 degrees and 180 degrees rotational angles. To explore the consequences of this near conformational energy degeneracy to NLO behavior, the eight low energy rotational isomers of FTC [Robinson, B. H.; et al. Chem. Phys. 1999, 245, 35] were investigated. This study provides the first-reported DFT-based calculation of the statistical mechanical average of beta over the conformational space of a molecule having substantial nonlinear optical behavior. The influence of the solvent reaction field on rotameric populations and on the beta tensor is reported. In the second part, two molecules having two donors and two acceptors bonded respectively in ortho and meta positions on a central benzene ring are shown to have substantially different beta tensors. These two so-called molecular Xs have different highest occupied molecular orbital to lowest unoccupied molecular orbital (HOMO-LUMO) distributions, and consistent with expectations, it is found that the larger beta(zzz) is associated with a large spatial asymmetry between the HOMOs and LUMOs. Large hyperpolarizability correlates with the HOMO concentrated on the donor groups and the LUMO on the acceptor groups.     

Probing electron correlations in molecules by two-dimensional coherent optical spectroscopy

The nonlinear optical signal generated in phenol by three femtosecond pulses with wavevectors k1, k2, and k3 in the phase-matching direction k1 + k2 - k3 is simulated. This two-dimensional coherent spectroscopy (2DCS) signal has a rich pattern containing information on double-excitation states. The signal vanishes for uncorrelated electrons due to interference among quantum pathways and, thus, provides direct signatures of correlated many-electron wavefunctions. This is illustrated by the very different 2DCS signals predicted by two levels of electronic structure calculations: state-averaged complete active space self-consistent field (SA-CASSCF) and multistate multiconfigurational second-order perturbation theory (MS-CASPT2).     

The optical properties of triphenylmethane dye molecules and chromogens

Industrial dye monomers, including malachite green, crystal violet, brilliant green, and methyl violet, were isolated by extraction with the use of heptane. UV light absorption bands characteristic of pure molecules were determined. The molecules of the dyes studied, which were ion pairs (formed by dye cations and oxalate or chlorine anions), did not absorb light in the visible range; that is, they were not chromogens. The conclusion was drawn that chromogen particles responsible for chromaticity were supramolecular dimers of nonchromogenic triphenylmethane series molecules. This conclusion was substantiated by trends in spectral transformations with the participation of immonium hydroxides obtained from dyes and side products of the synthesis of industrial dyes with quinoid molecular structures.     

Raman optical activity of tartaric acid and related molecules

The Raman optical activity spectra of (2R, 3R) (+)- and (2S,3S) (?)-tartaric acid, (2R, 3R) (+)-dimethyl tartrate, (2R,3R) (?)-2,3-butanediol and (2S, 3S) (+)-dibenzoyl tartaric acid are presented. A large couplet at about 500cm^?1 in the first three molecules, which probably originates in deformations of a chiral structural unit, might serve as an indicator of conformation and absolute configuration.