Engineering Papers⌕ Search

Engineering topics

Li, Xiaosong

Publications and source records attributed to Li, Xiaosong.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Time-dependent equation-of-motion coupled-cluster simulations with a defective Hamiltonian

Here simulations of laser-induced electron dynamics in a molecular system are performed using time-dependent (TD) equation-of-motion (EOM) coupled-cluster (CC) theory. The target system has been chosen to highlight potential shortcomings of truncated TD-EOM-CC methods [represented in this work by TD-EOM-CC with single and double excitations (TD-EOM-CCSD)], where unphysical spectroscopic features can emerge. Specifically, we explore driven resonant electronic excitations in magnesium fluoride in the proximity of an avoided crossing. Near the avoided crossing, the CCSD similarity-transformed Hamiltonian is defective, meaning that it has complex eigenvalues, and oscillator strengths may take on negative values. When an external field is applied to drive transitions to states exhibiting these traits, unphysical dynamics are observed. For example, the stationary states that make up the time-dependent state acquire populations that can be negative, exceed one, or even complex-valued.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Revealing Excited‐State Trajectories on Potential Energy Surfaces with Atomic Resolution in Real Time

Abstract Photoexcited molecular trajectories on potential energy surfaces (PESs) prior to thermalization are intimately connected to the photochemical reaction outcome. The excited‐state trajectories of a diplatinum complex featuring photo‐activated metal–metal σ ‐bond formation and associated Pt−Pt stretching motions were detected in real time using femtosecond wide‐angle X‐ray solution scattering. The observed motions correspond well with coherent vibrational wavepacket motions detected by femtosecond optical transient absorption. Two key coordinates for intersystem crossing have been identified, the Pt−Pt bond length and the orientation of the ligands coordinated with the platinum centers, along which the excited‐state trajectories can be projected onto the calculated PESs of the excited states. This investigation has gleaned novel insight into electronic transitions occurring on the time scales of vibrational motions measured in real time, revealing ultrafast nonadiabatic or non‐equilibrium processes along excited‐state trajectories involving multiple excited‐state PESs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Revealing Excited‐State Trajectories on Potential Energy Surfaces with Atomic Resolution in Real Time

Abstract Photoexcited molecular trajectories on potential energy surfaces (PESs) prior to thermalization are intimately connected to the photochemical reaction outcome. The excited‐state trajectories of a diplatinum complex featuring photo‐activated metal–metal σ ‐bond formation and associated Pt−Pt stretching motions were detected in real time using femtosecond wide‐angle X‐ray solution scattering. The observed motions correspond well with coherent vibrational wavepacket motions detected by femtosecond optical transient absorption. Two key coordinates for intersystem crossing have been identified, the Pt−Pt bond length and the orientation of the ligands coordinated with the platinum centers, along which the excited‐state trajectories can be projected onto the calculated PESs of the excited states. This investigation has gleaned novel insight into electronic transitions occurring on the time scales of vibrational motions measured in real time, revealing ultrafast nonadiabatic or non‐equilibrium processes along excited‐state trajectories involving multiple excited‐state PESs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Scalar Breit interaction for molecular calculations

Variational treatment of the Dirac–Coulomb–Gaunt or Dirac–Coulomb–Breit two-electron interaction at the Dirac–Hartree–Fock level is the starting point of high-accuracy four-component calculations of atomic and molecular systems. In this work, we introduce, for the first time, the scalar Hamiltonians derived from the Dirac–Coulomb–Gaunt and Dirac–Coulomb–Breit operators based on spin separation in the Pauli quaternion basis. While the widely used spin-free Dirac–Coulomb Hamiltonian includes only the direct Coulomb and exchange terms that resemble nonrelativistic two-electron interactions, the scalar Gaunt operator adds a scalar spin–spin term. The spin separation of the gauge operator gives rise to an additional scalar orbit-orbit interaction in the scalar Breit Hamiltonian. Benchmark calculations of Au n (n = 2–8) show that the scalar Dirac–Coulomb–Breit Hamiltonian can capture 99.99% of the total energy with only 10% of the computational cost when real-valued arithmetic is used, compared to the full Dirac–Coulomb–Breit Hamiltonian. In conclusion, the scalar relativistic formulation developed in this work lays the theoretical foundation for the development of high-accuracy, low-cost correlated variational relativistic many-body theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Electron scattering with ethane adsorbed on rare gas multilayers: Hole transfer, coulomb decay, and ion dissociation

Positive ion desorption following electron impact dissociative ionization of ethane adsorbed on Ar, Kr, and Xe multilayers has been studied as a function of incident electron energy from threshold to 100 eV. Based on the dependence of ion yields on the identity of the rare gas, it is likely that the majority of ethane molecules undergo indirect ionization following hole transfer from the ionized underlying rare gas. Further, this has also been corroborated by density of states calculations showing the energetic alignment of the outer valence states of ethane and the condensed rare gas ionization energies. Due to the near-resonant nature of charge transfer for single-hole states, the ethane molecular ion is excited to different final ionic states on different rare gases, which leads to differences in ion desorption yields and branching ratios. The quantitative yields increase with increasing ionization energy gap between the rare gas and ethane, in the order Ar > Kr >Xe. The large increase in yields from 25 eV onwards for all rare gases is likely due to the formation and decay of two-hole states on neighboring rare gas and ethane molecules due to interatomic and intermolecular Coulomb decay (ICD) and not electron transfer mediated decay (ETMD). The ICD and ETMD pathways become accessible when the incoming electron has sufficient energy to excite the inner valence ns level of the rare gas to a Rydberg state or ionize it. The experimental findings are supported by calculations of thresholds, density of states for the final configurations of these processes, and coupling strengths for hole transfer between ethane and rare gases. The fragment ion branching ratios vary with energy from threshold to about 35 eV, showing the fragmentation pattern changes with the mode of hole transfer and availability of excess energy. Sigma C–C bonds are more likely to break than C–H bonds in the mid-20 eV range, and this effect is most pronounced for Xe, followed by Kr, and then Ar.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Solvent Induced Proton Polarization within the Nuclear–Electronic Orbital Framework

To explicitly account for nuclear quantum effects and solvent environments in simulations of chemical processes, the nuclear-electronic orbital approach is coupled with a polarizable continuum model (PCM). This NEO–PCM approach is used to explore the influence of solvation on nuclear polarization through applications to a water dimer and a set of protonated water tetramers. Nuclear polarization in these species is analyzed in terms of changes in proton density and oxygen–hydrogen bond length. Solvation is shown to enhance nuclear polarization with increasing dielectric constant. For the water dimer, the internal, hydrogen-bonded proton is shown to polarize more than the external, free proton. Moreover, proton quantization leads to greater solvent polarization through their mutual polarization. Furthermore, these calculations highlight the complex interplay among electronic, nuclear, and solvent polarization in chemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Bridge-Mediated Metal-to-Metal Electron and Hole Transfer in a Supermolecular Dinuclear Complex: A Computational Study Using Quantum Electron–Nuclear Dynamics

Bimetallic electron donor–acceptor complexes can facilitate electron and energy transfer with excellent structural control through synthetic design. In this work, we investigate the photochemical dynamics in a Ru–Cu bimetallic complex after photoexcitation of the Ru-centered charge transfer state. The physical underpinnings of the metal-to-metal directional charge transfer process are unraveled via analyses of the quantum electronic dynamics and electron–nuclear trajectories. Here, the effects of molecular vibrations in the photoexcited state on the charge transfer processes are also analyzed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field methods

The fully correlated frequency-independent Dirac–Coulomb–Breit Hamiltonian provides the most accurate description of electron–electron interaction before going to a genuine relativistic quantum electrodynamics theory of many-electron systems. In this work, we introduce a correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field method within the frameworks of complete active space and density matrix renormalization group. In this approach, the Dirac–Coulomb–Breit Hamiltonian is included variationally in both the mean-field and correlated electron treatment. Here, we also analyze the importance of the Breit operator in electron correlation and the rotation between the positive- and negative-orbital space in the no-virtual-pair approximation. Atomic fine-structure splittings and lanthanide contraction in diatomic fluorides are used as benchmark studies to understand the contribution from the Breit correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Simultaneous Optimization of Nuclear–Electronic Orbitals

Accurate modeling of important nuclear quantum effects, such as nuclear delocalization, zero-point energy, and tunneling, as well as non-Born-Oppenheimer effects, requires treatment of both nuclei and electrons quantum mechanically. The nuclear–electronic orbital (NEO) method provides an elegant framework to treat specified nuclei, typically protons, on the same level as the electrons. In conventional electronic structure theory, finding a converged ground state can be a computationally demanding task; converging NEO wavefunctions, due to their coupled electronic and nuclear nature, is even more demanding. Herein, we present an efficient simultaneous optimization method that uses the direct inversion in the iterative subspace method to simultaneously converge wavefunctions for both the electrons and quantum nuclei. In conclusion, benchmark studies show that the simultaneous optimization method can significantly reduce the computational cost compared to the conventional stepwise method for optimizing NEO wavefunctions for multicomponent systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Efficient evaluation of the Breit operator in the Pauli spinor basis

The frequency-independent Coulomb–Breit operator gives rise to the most accurate treatment of two-electron interaction in the non-quantum-electrodynamics regime. The Breit interaction in the Coulomb gauge consists of magnetic and gauge contributions. The high computational cost of the gauge term limits the application of the Breit interaction in relativistic molecular calculations. In this work, we apply the Pauli component integral–density matrix contraction scheme for gauge interaction with a maximum spin- and component separation scheme. We also present two different computational algorithms for evaluating gauge integrals. One is the generalized Obara–Saika algorithm, where the Laplace transformation is used to transform the gauge operator into Gaussian functions and the Obara–Saika recursion is used for reducing the angular momentum. The other algorithm is the second derivative of Coulomb interaction evaluated with Rys-quadrature. This work improves the efficiency of performing Dirac–Hartree–Fock with the variational treatment of Breit interaction for molecular systems. In conclusion, we use this formalism to examine relativistic trends in the Periodic Table and analyze the relativistic two-electron interaction contributions in heavy-element complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Ligand-Structure-Dependent Coherent Vibrational Wavepacket Dynamics in Pyrazolate-Bridged Pt(II) Dimers

Bimetallic transition metal complexes have gained increasing attention because of their versatile functions in solar energy conversion and photonics applications arising from inter-metal electronic coupling. In bimetallic platinum (Pt) complexes, electronic communication between the Pt-centered and ligand-centered moieties have been shown to be critical for defining their excited-state dynamic trajectories undergoing either localized ligand centered (LC)/metal-to-ligand-charge-transfer (MLCT) transitions or delocalized metal-metal-to-ligand-charge-transfer (MMLCT) transitions. The branching of the excited-state intersystem crossing (ISC) trajectories are modulated through structural factors that alter the relative energies of the different states. In this study, we investigated the correlation of the structural factors influencing the excited state trajectories. Using femtosecond broadband transient absorption (fs-BBTA) spectroscopy, ultrafast dynamics in the excited state of two select Pt(II) dimers have been mapped out using their coherent vibrational wavepacket signatures in corresponding transient absorption spectra. To examine how the ligand moieties of the Pt(II) dimers influence excited-state dynamics and the coherent vibrational wavepacket behavior, here we carried out comparative studies on two pyrazolate-bridged Pt(II) dimers of the general formula [Pt( t Bu 2 Pz)(N^C)] 2 ( t Bu 2 Pz = 3,5-di-tert-butylpyrazole); N^C = 7,8-benzoquinoline (bzq, 1) or 1-phenylisoquinoline (piq, 2)). We found that photoexcitation into the low energy absorption bands of 1 and 2 respectively induce the formation of 1 MMLCT states from which ultrafast ISC proceeds, resulting in stimulated emission quenching and decoherence of the vibrational wavepacket motions. The results obtained in this study suggest that both energetics and the structural rigidity of the aromatic cyclometalating ligands in 1 and 2 can significantly influence dynamics along the excited state trajectory characterized by dephasing of the coherent oscillations. The collective results provide direct evidence of how ligand structure alters electronic dynamics along excited state trajectories associated with ISC processes, providing insight into using ligand design to steer photochemical processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Proton-Coupled Electron Transfer in a Ruthenium(II) Bipyrimidine Complex in Its Ground and Excited Electronic States

In this work, proton-coupled electron transfer (PCET) was studied for the ground and excited electronic states of a [Ru(terpy)(bpm)(OH 2 )(PF 6 ) 2 ] complex, Ru-bpm. Cyclic voltammetry measurements show that the Ru(II)-aqua moiety undergoes PCET to form a Ru(IV)-oxo moiety in the anodic region, while the bpm ligand undergoes PCET to form bpmH 2 in the cathodic region. The photophysical behavior of Ru-bpm was studied using steady-state and femtosecond transient UV–vis absorption spectroscopy, coupled with density functional theory (DFT) calculations. The lowest-lying excited state of Ru-bpm is described as a (Ru → bpm) metal-to-ligand charge-transfer (MLCT) state, while the metal-centered (MC) excited state was found computationally to be close in energy to the lowest-energy bright MLCT state (MC state was 0.16 eV above the MLCT state). The excited-state kinetics of Ru-bpm were found via transient absorption spectroscopy to be short-lived and were fit well to a biexponential function with lifetimes $τ_1$ = 4 ps and $τ_2$ = 65 ps in aqueous solution. Kinetic isotope effects of 1.75 ($τ_1$) and 1.61 ($τ_2$) were observed for both decay components, indicating that the solvent plays an important role in the excited-state dynamics of Ru-bpm. Based on the pH-dependent studies and the results from prior studies of similar Ru-complexes, we hypothesize that the 3 MLCT state forms an excited-state hydrogen-bond adduct with the solvent molecules and that this process occurs with a 4 ps lifetime. The formation of such a hydrogen-bond complex is consistent with the electronic density accumulation at the peripheral N atoms of the bpm moiety in the 3 MLCT state. The hydrogen-bonded state 3 MLCT decays to the ground state with a 65 ps lifetime. Such a short lifetime is likely associated with the efficient vibrational energy transfer from the 3 MLCT state to the solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Ultrafast branching in intersystem crossing dynamics revealed by coherent vibrational wavepacket motions in a bimetallic Pt($\mathrm{II}$) complex

Ultrafast excited-state processes of transition metal complexes (TMCs) are governed by complicated interplays between electronic and nuclear dynamics, which demand a detailed understanding to achieve optimal functionalities of photoactive TMC-based materials for many applications. Here, in this work, we investigated a cyclometalated platinum(II) dimer known to undergo a Pt-Pt bond contraction in metal-metal-to-ligand-charge-transfer (MMLCT) excited state using femtosecond broadband transient absorption (fs-BBTA) spectroscopy in combination with geometry optimization and normal mode calculations. Using sub-20 fs pump and broadband probe pulses in fs-BBTA spectroscopy, we are able to correlate coherent vibrational wavepacket (CVWP) evolution with stimulated emission (SE) dynamics of the 1 MMLCT state. The results demonstrated that the 145 cm -1 CVWP motions with the damping times of ~ 0.9 ps and ~ 2 ps originate from coherent Pt-Pt stretching vibrations in the singlet and triplet MMLCT states, respectively. On the basis of excited-state potential energy surface calculations in our previous work, we rationalized that the CVWP transfer from Franck-Condon (FC) state to the 3 MMLCT state was mediated by a triplet ligand-centered ( 3 LC) intermediate state through two step intersystem crossing (ISC) on a time scale shorter than a period of the Pt-Pt stretching wavepacket motions. Moreover, it was found that the CVWP motion with 110 cm-1 frequency decays with the damping time of ~ 0.2 ps, well matching the time constant of 0.253 ps corresponding to a redshift in the SE feature at early time. This observation indicates that the Pt-Pt bond contraction changes the stretching frequency from 110 to 145 cm -1 and stabilize the 1 MMLCT state relative to the 3 LC state with the ~ 0.2 ps time scale. Thus, the ultrafast ISC from the 1 MMLCT to the 3 LC states occurs before the Pt-Pt bond shortening. The findings herein provide insight into understanding the impact of Pt-Pt bond contraction on ultrafast branching of the 1 MMLCT population into the direct ( 1 MMLCT → 3 MMLCT) and indirect ISC pathways ( 1 MMLCT → 3 LC → 3 MMLCT) in the Pt(II) dimer. These results revealed intricate excited state electronic and nuclear motions that could steer the reaction pathways with levels of details that have not achieved before

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗