Stark control of electrons across the molecule-semiconductor interface

Controlling matter at the level of electrons using ultrafast laser sources represents an important challenge for science and technology. Recently, we introduced a general laser control scheme (the Stark control of electrons at interfaces or SCELI) based on the Stark effect that uses the subcycle structure of light to manipulate electron dynamics at semiconductor interfaces [A. Garzón-Ramírez and I. Franco, Phys. Rev. B 98, 121305 (2018)]. Here, we demonstrate that SCELI is also of general applicability in molecule-semiconductor interfaces. We do so by following the quantum dynamics induced by non-resonant few-cycle laser pulses of intermediate intensity (non-perturbative but non-ionizing) across model molecule-semiconductor interfaces of varying level alignments. We show that SCELI induces interfacial charge transfer regardless of the energy level alignment of the interface and even in situations where charge exchange is forbidden via resonant photoexcitation. We further show that the SCELI rate of charge transfer is faster than those offered by resonant photoexcitation routes as it is controlled by the subcycle structure of light. The results underscore the general applicability of SCELI to manipulate electron dynamics at interfaces on ultrafast timescales.

Tensor-Train Split-Operator KSL (TT-SOKSL) Method for Quantum Dynamics Simulations

Numerically exact simulations of quantum reaction dynamics, including nonadiabatic effects in excited electronic states, are essential to gain fundamental insights into ultrafast chemical reactivity and rigorous interpretations of molecular spectroscopy. Here, we introduce the tensor-train split-operator KSL (TT-SOKSL) method for quantum simulations in tensor-train (TT)/matrix product state (MPS) representations. TT-SOKSL propagates the quantum state as a tensor train using the Trotter expansion of the time-evolution operator, as in the tensor-train split-operator Fourier transform (TT-SOFT) method. However, the exponential operators of the Trotter expansion are applied using a rank-adaptive TT-KSL scheme instead of using the scaling and squaring approach as in TT-SOFT. We demonstrate the accuracy and efficiency of TT-SOKSL as applied to simulations of the photoisomerization of the retinal chromophore in rhodopsin, including nonadiabatic dynamics at a conical intersection of potential energy surfaces. The quantum evolution is described in full dimensionality by a time-dependent wavepacket evolving according to a two-state 25-dimensional model Hamiltonian. We find that TT-SOKSL converges faster than TT-SOFT with respect to the maximally allowed memory requirement of the tensor-train representation and better preserves the norm of the time-evolving state. When compared to the corresponding simulations based on the TT-KSL method, TT-SOKSL has the advantage of avoiding the need to construct the matrix product state Laplacian by exploiting the linear scaling of multidimensional tensor-train Fourier transforms.

Functional Tensor-Train Chebyshev Method for Multidimensional Quantum Dynamics Simulations

Methods for efficient simulations of multidimensional quantum dynamics are essential for theoretical studies of chemical systems where quantum effects are important, such as those involving rearrangements of protons or electronic configurations. Here, we introduce the functional tensor-train Chebyshev (FTTC) method for rigorous nuclear quantum dynamics simulations. FTTC is essentially the Chebyshev propagation scheme applied to the initial state represented in a continuous analogue tensor-train format. We demonstrate the capabilities of FTTC as applied to simulations of proton quantum dynamics in a 50-dimensional model of hydrogen-bonded DNA base pairs.

Iterative Power Algorithm for Global Optimization with Quantics Tensor Trains

Optimization algorithms play a central role in chemistry since optimization is the computational keystone of most molecular and electronic structure calculations. Herein, we introduce the iterative power algorithm (IPA) for global optimization and a formal proof of convergence for both discrete and continuous global search problems, which is essential for applications in chemistry such as molecular geometry optimization. IPA implements the power iteration method in quantics tensor train (QTT) representations. Analogous to the imaginary time propagation method with infinite mass, IPA starts with an initial probability distribution ρ₀(x) and iteratively applies the recurrence relation ρ_k+1(x) = U(x) ρ_k(x)/∥Uρ_k∥_L¹, where U(x) = e^-V(x) is defined in terms of the potential energy surface (PES) V(x) with global minimum at x = x*. Upon convergence, the probability distribution becomes a delta function δ(x - x*), so the global minimum can be obtained as the position expectation value x* = Tr[x δ(x - x*)]. QTT representations of V(x) and ρ(x) are generated by fast adaptive interpolation of multidimensional arrays to bypass the curse of dimensionality and the need to evaluate V(x) for all possible values of x. We illustrate the capabilities of IPA for global search optimization of two multidimensional PESs, including a differentiable model PES of a DNA chain with D = 50 adenine-thymine base pairs, and a discrete non-differentiable potential energy surface, V(p) = mod(N,p), that resolves the prime factors of an integer N, with p in the space of prime numbers {2, 3,..., p_max} folded as a d-dimensional 2₁ × 2₂ × ··· × 2_d tensor. We find that IPA resolves multiple degenerate global minima even when separated by large energy barriers in the highly rugged landscape of the potentials. Therefore, IPA should be of great interest for a wide range of other optimization problems ubiquitous in molecular and electronic structure calculations.

Vibronic Dynamics of Photodissociating ICN from Simulations of Ultrafast X-Ray Absorption Spectroscopy

Ultrafast UV-pump/soft-X-ray-probe spectroscopy is a subject of great interest since it can provide detailed information about dynamical photochemical processes with ultrafast resolution and atomic specificity. Here, we focus on the photodissociation of ICN in the 1 Π1 excited state, with emphasis on the transient response in the soft-X-ray spectral region as described by the ab initio spectral lineshape averaged over the nuclear wavepacket probability density. We find that the carbon K-edge spectral region reveals a rich transient response that provides direct insights into the dynamics of frontier orbitals during the I-CN bond cleavage process. The simulated UV-pump/soft-X-ray-probe spectra exhibit detailed dynamical information, including a time-domain signature for coherent vibration associated with the photogenerated CN fragment.

Floquet Study of Quantum Control of the Cis-Trans Photoisomerization of Rhodopsin

Understanding how to control reaction dynamics of polyatomic systems by using ultrafast laser technology is a fundamental challenge of great technological interest. Here, we report a Floquet theoretical study of the effect of light-induced potentials on the ultrafast cis-trans photoisomerization dynamics of rhodopsin. The Floquet Hamiltonian involves an empirical 3-state 25-mode model with frequencies and excited-state gradients parametrized to reproduce the rhodopsin electronic vertical excitation energy, the resonance Raman spectrum, and the photoisomerization time and efficiency as probed by ultrafast spectroscopy. We simulate the excited state relaxation dynamics using the time-dependent self-consistent field method, as described by a 3-state 2-mode nuclear wavepacket coupled to a Gaussian ansatz of 23 vibronic modes. We analyze the reaction time and product yield obtained with pulses of various widths and intensity profiles, defining 'dressed states' where the perturbational effect of the pulses is naturally decoupled along the different reaction channels. We find pulses that delay the excited-state photoisomerization for hundreds of femtoseconds, and we gain insights on the underlying control mechanisms. The reported findings provide understanding of quantum control, particularly valuable for the development of ultrafast optical switches based on visual pigments.

Sensing of molecules using quantum dynamics

We design sensors where information is transferred between the sensing event and the actuator via quantum relaxation processes, through distances of a few nanometers. We thus explore the possibility of sensing using intrinsically quantum mechanical phenomena that are also at play in photobiology, bioenergetics, and information processing. Specifically, we analyze schemes for sensing based on charge transfer and polarization (electronic relaxation) processes. These devices can have surprising properties. Their sensitivity can increase with increasing separation between the sites of sensing (the receptor) and the actuator (often a solid-state substrate). This counterintuitive response and other quantum features give these devices favorable characteristics, such as enhanced sensitivity and selectivity. Using coherent phenomena at the core of molecular sensing presents technical challenges but also suggests appealing schemes for molecular sensing and information transfer in supramolecular structures.

Photoinduced Ultrafast Heterogeneous Electron Transfer at Molecule-Semiconductor Interfaces

This Perspective discusses recent developments in ultrafast electron transfer dynamics at interfaces between organic and inorganic materials. Heterogeneous electron transfer (HET) is a key process in important fields like catalysis and solar energy conversion. Furthermore, the solid state nature of the systems gives control over relevant parameters and allows for investigating excited state dynamics and electron transfer processes in unprecedented detail. Progress in synthesis, sample preparation, and instrumentation makes it possible to provide experimental proof of recent prediction from theory concerning the adiabaticity of the reaction and the influence of coherence. A short recapitulation of the field is followed by a discussion of recent experimental efforts that allowed for studying HET, particularly focusing on the influence of energetics and vibrational dynamics.

Electroencephalographic field influence on calcium momentum waves

Macroscopic electroencephalographic (EEG) fields can be an explicit top-down neocortical mechanism that directly drives bottom-up processes that describe memory, attention, and other neuronal processes. The top-down mechanism considered is macrocolumnar EEG firings in neocortex, as described by a statistical mechanics of neocortical interactions (SMNI), developed as a magnetic vector potential A. The bottom-up process considered is Ca(2+) waves prominent in synaptic and extracellular processes that are considered to greatly influence neuronal firings. Here, the complimentary effects are considered, i.e., the influence of A on Ca(2+) momentum, p. The canonical momentum of a charged particle in an electromagnetic field, Π=p+qA (SI units), is calculated, where the charge of Ca(2+) is q=-2e, e is the magnitude of the charge of an electron. Calculations demonstrate that macroscopic EEG A can be quite influential on the momentum p of Ca(2+) ions, in both classical and quantum mechanics. Molecular scales of Ca(2+) wave dynamics are coupled with A fields developed at macroscopic regional scales measured by coherent neuronal firing activity measured by scalp EEG. The project has three main aspects: fitting A models to EEG data as reported here, building tripartite models to develop A models, and studying long coherence times of Ca(2+) waves in the presence of A due to coherent neuronal firings measured by scalp EEG. The SMNI model supports a mechanism wherein the p+qA interaction at tripartite synapses, via a dynamic centering mechanism (DCM) to control background synaptic activity, acts to maintain short-term memory (STM) during states of selective attention.

Quantum dephasing of a two-state system by a nonequilibrium harmonic oscillator

In this paper, we investigate coherent quantum dynamics in a nonequilibrium environment. We focus on a two-state quantum system strongly coupled to a single classical environmental oscillator, and explore the effect of nonstationary statistical properties of the oscillator on the quantum evolution. A simple nonequilibrium model, consisting of an oscillator with a well-defined initial phase which undergoes subsequent diffusion, is introduced and studied. Approximate but accurate analytic expressions for the evolution of the off-diagonal density matrix element of the quantum system are derived in the second-order cumulant approximation. The effect of the initial phase choice on the subsequent quantum evolution is quantified. It is observed that the initial phase can have a significant effect on the preservation of coherence on short time scales, suggesting this variable as a control parameter for optimizing coherence in many-body quantum systems.

Exploring the role of phase modulation on photoluminescence yield

We report an investigation to elucidate the mechanisms of control in phase-sensitive experiments in two molecular systems. A first inspection of optimization procedures yields the same experimental result: increase in the emission efficiency upon excitation by a phase modulated pulse in a two-photon transition. More detailed studies, which include power dependence, spectral response, one and two color pump-probe and pump-pump experiments show that while for one chromophore phase modulation leads to spectral matching between the two-photon cross section and the second order power spectrum for the other it provides a tool to manipulate the wavepacket dynamics in the excited state.

Applying genetic algorithm optimization to a folded geometry acousto-optic modulated spatial pulse shaper

A folded geometry acousto-optic modulator spatial pulse shaper has been designed for shaping individual pulses from a high power amplified laser. The design preserves the capability of computer programmable amplitude and phase modulation of femtosecond laser pulses. An additional application of genetic algorithm optimization approach for compressing a stretched pulse is also demonstrated for such a pulse shaper. Spectrally and temporally resolved optical gating technique is used to characterize the shaped pulses.