[ Home ] [ Publications ] [ Bio ] [ Research ] [ Jobs ] [ Group ] [ Covers ] [ Cool movies ] [ Photos ] [ Visitor Info ]


Electronic structure and spectroscopy of organic and inorganic nano-materials

Our research is focused on theoretical studies of excited states and optical responses of photoactive organic and inorganic nano-materials such as large organic molecules, biological complexes and semiconductor nano-particles. All electronic processes in these systems are governed by strong electronic correlations and coupling of electrons to molecular structure. We develop quantum chemical approaches for efficient and accurate modeling of electronic structure, dynamics, charge and energy transfer in large molecular systems and apply these methods to different materials such as conjugated polymers, dendrimers, carbon nanotubes, biological light harvesting systems, donor-acceptor complexes, and semiconductor quantum dots. The results of the calculations are then used for modeling of linear (absorption and emission) and nonlinear (including ultrafast frequency and time-resolved spectroscopic probes) optical responses to model experimental data and to understand underlying dynamical photo-physical, photo-chemical and relaxation processes (see our reviews on electronic processes, excited state dynamics, halide perovskites, organic photovoltaics, functionalized carbon nanotubes and machine learning). Our current research projects include:

Non-adiabatic excited state molecular dynamics (NEXMD)

Almost every photochemical, photophysical, charge and energy transport, or spectroscopic process in organic nano-materials involves dynamics of excited vibronic states and ultimately defines molecular electronic properties and functionalities. When excited by photons of light, functional electronic materials show dynamics that are often characterized by large nonadiabatic couplings between multiple excited states through a breakdown of the Born−Oppenheimer approximation. Following photoexcitation, various nonradiative intraband relaxation pathways can lead to a number of complex processes. Therefore, computational simulation of nonadiabatic molecular dynamics is an indispensable tool for understanding complex light-driven processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons. Over the years, we have developed a nonadiabatic excited-state molecular dynamics (NEXMD) framework that efficiently and accurately describes photoinduced phenomena in extended molecular chromophores. We use the fewest-switches surface hopping (FSSH), Erhenfest, and multiconfigurational Ehrenfest ab-initio multiple cloning (MCE-AIMC) algorithms to treat quantum transitions among multiple adiabatic excited state potential energy surfaces. More accurate numerical algorithms for the NEXMD such as MCE-AIMC allows treatment of decoherences and interferences beyond the basic FSSH method. Extended molecular systems often contain hundreds of atoms and involve large densities of excited states that participate in the photoinduced dynamics. The semiempirical NEXMD methodology freely distributed as the NEXMD computational package under the BSD-3 License offers a computationally tractable route for simulating hundreds of atoms on 10 ps time scales where multiple coupled excited states are involved. Furthermore, our aproach has been implemented into NWChem computational package that allows for FSSH, Ehrenfest or MCE-AIMC non-adiabatic molecular dynamics simulations at first-principles time-dependent Density Functional Theory (TD-DFT) level.  The NEXMD technique was successfully applied to study ultrafast excited state dynamics and spectroscopies in conjugated polyfluorenes, phenylene-ethynylenes, dendrimers, cycloparaphenylenes, chlorophylls, explosives and many other molecular materials. Section Movies exemplifies examples of complex evolution of photoexcited wavefunctions in several systems. In addition we are exploring dielectric environment effects on excited state dynamics relying on state-specific solvation within polarizable continuum and quantum/mechanics/molecular mechanics (QM/MM) models.

Design of hybrid perovskite materials for photovoltaics and optoelectronics

The recent discovery of organic-inorganic perovskites offers promising routes for the development of low-cost, solar-based clean global energy solutions for the future Solution-processed organic-inorganic hybrid perovskite planar solar cells, such as CH3NH3PbX3 (X = Cl, Br, I), have achieved high average power conversion efficiency (PCE) values approaching 20%. In 2015 our experimental-theory team at LANL has demonstrated a universally applicable fabrication strategy for the growth of millimeters scale crystals of hybrid perovskites. This allows making hysteresis-free solar cells approaching record 18% light-to-energy conversion efficiency with extraordinary reproducibility and negligible degradation. Our device modeling and DFT calculations characterize charge dynamics at the interface of perovskites in order to aide in materials design and device engineering. Our follow up studies address photodegradation processes in perovskite solar cells. Our theory and experiments elucidate the atomistic origin of light-induced macroscopic degradation tracing this phenomenon to appearance to meta-stable charge trap states (polarons) and then utilize this knowledge to demonstrate photovoltaic devices stable with constant illumination.  Finally, our recent results are on 2D layered hybrid perovskites, which being incorporated into a photovoltaic device exhibited world-record efficiency (13% in the category of 2D materials) and stability for several months under constant light illumination and 65% humidity. The theoretical follow up review aims to provide an interpretive and predictive framework for 3D and 2D layered hybrid perovskite optoelectronic properties. Lately we have been exploring applications of halide perovskite matetrials to X-ray and Gamma-ray detection. Theoretical studies exploring the effects of doping, non-adiabatic dynamics or polarons provide important guidelines to experiment. 

Machine Learning for chemistry and materials properties

Machine Learning (ML) is quickly becoming a premier tool for modeling chemical processes and materials. ML algorithmes for various properties of interest, trained on large datasets of high-quality electronic structure calculations, are particularly attractive due their unique combination of computational efficiency and physical accuracy.  Designing high quality training datasets is crucial to overall model accuracy. One strategy is active learning, in which new data is automatically collected for atomic configurations that produce large ML uncertainties. Another strategy is to use the highest levels of quantum theory possible. Transfer learning allows training to a dataset of mixed fidelity. A model initially trained to a large dataset of density functional theory calculations and can be significantly improved by retraining to a relatively small dataset of expensive coupled cluster theory calculations. The locality approximation underpinning favorable computational scaling of the ML models, is another severe limitation that fails to capture long-range effects that may arise from charge transfer, polarization, electrostatic or dispersion interactions. ML models can overcome nonlocality via introduction of interaction layers, self-consistent cycles, or charge equilibration schemes. These advances are applied to predict such properties as molecular energies, interatomic potentials, charge distributions, charged species, excited states and effective Hamiltonian models,  exemplified by applications to molecules and materials.

Energy and charge transport in macroscopic molecular assemblies

Following absorption (or leading to emission) of light quanta are fundamentals processes of energy (e.g., in a form of excitons) and charge (e.g., in a form of polarons) transport.  Theoretical modeling of these phenomena in large molecular assemblies have evolved from a consideration of simple models to a more or less comprehensive description utilizing first principles calculations of underlining molecular using toward accessing device scale. We have been developing multiscale modeling toolbox for simulation of energy (using Foester theory) and charge carrier mobility (using Marcus theory) in soft electronic materials, to establish a set of structure-property relationships for rational design of such materials and to suggest new materials based on theoretical prediction. Applications were done for small molecule-based organic photovoltaic materials and natural biological systems containing aromatic amino acid residues. This modeling was further extended to non-fullerene acceptor (NFA) materials, where lattice organization defines anisotropic charge transport.

Spectroscopy and electronic processes in dunctionalized carbon nanotube and graphene materials

Single-wall carbon nanotubes (SWCNTs) are one-atom-thick layer of graphene rolled into long rigid cylinders a few nanometers in diameters. The diameter and direction of tube’s rolling (a chiral vector (n,m)) define the main features of SWCNTs, such as a semiconductor or metal-like electronic structure. Carbon nanotubes constitute a class of the most promising technological materials with multifarious applications ranging from electronic chips and optoelectronics to medical technologies and sensing. Over the past two years we pioneered applications of methodologies previously developed for molecular materials to SWCNTs (in contrast to solid-state approaches being standard in the field). These studies made possible to address many aspects of photophysics of an isolated tube. This includes anharmonic coherent phonon dynamics, quantification of exciton-phonon coupling and Huang-Rhys factors,  characterization of delocalized non-excitonic transitions, effects of Peierls distortion and exciton self-trapping, energetics of dark states, prediction of triplet states, study of transverse polarized absorption in nanotubes, and understanding structure of DNA/SWCNT complexes. Many of these results confirm the latest experimental data. It turns out, control over the reorganization energies for polaronic and excitonic states suggest significant potential for tailoring the electronic properties of carbon nanostructures through chemical engineering. In particular, the intentional introduction of fluorescent quantum defects in the sidewalls of semiconducting single-wall carbon nanotubes through chemical reaction is an emerging route to predictably modify nanotube electronic structures and develop advanced photonic functionality such as high photoluminescence quantum yield and room-temperature single-photon emission at telecom wavelengths. Low-level covalent functionalization provides additional synthetic tunability in emission properties and imparts quantum emission function, but also acts as a source of spectral diversity. Quantum-chemical calculations help to establish structure-property relationships  reveling precise binding configurations and dependencies chemical structure of defects, coupling between patterned defects, nanotube chirality, etc. 

 

Electronic interactions in semiconductor quantum dots

Colloidal semiconductor quantum dots (QD) can be viewed as nanometer size objects in which quantum confinement occurs in all three dimensions. In QDs, the bulk band structure collapses into discrete atomic-like levels which makes it possible to view QDs as artificial “tunable” atoms with electronic energies that are controlled by QD dimensions. Following the analogy with atoms, QDs can conceivably be manipulated into complex interacting structures that mimic molecules and atomic or molecular solids. Such complex structures possess huge technological potential associated with the versatile tunability of their properties provided by the combined control of energy levels in the individual nanoparticles and collective interactions of nanoparticles in an assembly. First principle calculations of QD (supermolecular approach) presents a significant scientific challenge since the number of atoms grows rapidly with the system size. We are applying a complementary suit of quantum chemical methods ranging from first principles (TDDFT) to semiempirical and minimal models for quantitative evaluation of electronic structure in QD nanosystems working closely with our experimental colleagues. In particular, we have shown that structurally optimized QD assemblies can feature very efficient transfer rate of excitons from smaller to larger dots via electrostatic coupling approaching picosecond time scales. This leads to interesting energy harvesting applications mimicking natural biological complexes. Using density functional theory (DFT) and time-dependent DFT quantum-chemical methodologies, we recently explored the role of surface ligands on the electronic structure and observe strong surface-ligand interactions in quantum dots leading to formation of hybridized states and polarization effects. This opens new relaxation channels for high energy photoexcitations affecting both population (inelastic scattering) and dephasing (elastic scattering) processes. 

Multiscale modeling of electronic excitations using Exciton Scattering (ES) approach

Modern electronic structure calculations involve solving the Schrodinger equation. Substantial numerical cost [O(N3) - O(N5), N being the number of electrons (orbitals) in the system] makes repetitive computations for large molecules prohibitively expensive. This challenge calls for development of multiscale approaches. Notably, any reduction of computational complexity is highly non-trivial, since one need to take an advantage of quantum-mechanical phenomena, such as delocalized wavefunctions, quantum confinement, and electronic coherences. Using concept of strongly bound excitons, we recently developed multiscale exciton scattering (ES) model, which attributes excited states to standing waves in quasi-one-dimensional structures. This method allows electronic spectra for any structure of arbitrary size within the considered molecular family to be obtained with negligible numerical effort. Complex and non-trivial delocalization patterns of photoexcitations throughout the entire molecular structure can then be universally characterized and understood using the proposed ES method, completely bypassing 'supramolecular' calculations. An agreement is within 1-5 meV for most test cases, when the ES results are compared with the reference TD-DFT calculations has been achieved. Consequently, computational design of molecular structures with desirable electronic and optical properties can be performed in real-time using graph-like representation of molecules. Currently we are extending the ES model to calculate various spectroscopic observables and are applying this approach to study the structure of electronic excitations in macromolecules where it constitutes the only feasible option.

 

Nonlinear optical response of conjugated organics

Optical materials with enhanced nonlinear optical (NLO) response have important technological implications such as optical switching and wave-guiding, compact 3D data storage and micro fabrication, chemical and biological sensing, optical power limiting, up-conversion lasing, bio-imaging, etc. Computational design of nonlinear optical materials constitutes a complex fundamental problem and theoretical modeling of optical polarizabilities of molecular systems requires extensive numerical effort since NLO responses often involve strongly correlated excited electronic states. Recently we developed computational approach for modeling nonlinear optical polarizabilities in the framework of adiabatic time-dependent density functional theory (TDDFT) by deriving the respective nonlinear response functions from the equation of motion for the driven single-electron density matrix. This formalism has been coded and applied to calculate one- and two-photon absorption spectra (related to linear and third-order optical responses, respectively) in a series of large donor-acceptor substituted conjugated molecules. Calculated excitation energies corresponding to one and two-photon-absorption maxima are found to be in excellent agreement with experiment with better than 4% accuracy in average which shows a potential of TDDFT for NLO calculations. Using developed TDDFT extensions, we are studying other nonlinear optical responses and molecules including dipolar, quadrupolar, octupolar, branching and dendrimeric systems. Recently, we have investigated two-photon absorption in conjugated energetic molecules constituting a new class of explosives having potential for controlled direct optical initiation through photochemical pathways.