Search NASA⌕ Search

Engineering topics

Bhaskaran-Nair, Kiran

Publications and source records attributed to Bhaskaran-Nair, Kiran.

Petascale Simulations of the Morphology and the Molecular Interface of Bulk Heterojunctions

Understanding how additives interact and segregate within bulk heterojunction (BHJ) thin films is critical for exercising control over structure at multiple length scales and delivering improvements in photovoltaic performance. The morphological evolution of poly(3-hexylthiophene) (P3HT) and phenyl-C 61 -butyric acid methyl ester (PCBM) blends that are commensurate with the size of a BHJ thin film is examined using petascale coarse-grained molecular dynamics simulations. When comparing 2 component and 3 component systems containing short P3HT chains as additives undergoing thermal annealing we demonstrate that the short chains alter the morphol- ogy in apparently useful ways: They efficiently migrate to the P3HT/PCBM interface, increasing the P3HT domain size and interfacial area. Simulation results agree with depth profiles determined from neutron reflectometry measurements that reveal PCBM enrichment near substrate and air interfaces, but a decrease in that PCBM enrich- ment when a small amount of short P3HT chains are integrated into the BHJ blend. Atomistic simulations of the P3HT/PCBM blend interfaces show a non-monotonic dependence of the interfacial thickness as a function of number of repeat units in the oligomeric P3HT additive, and the thiophene rings orient parallel to the interfacial plane as they approach the PCBM domain. Using the nanoscale geometries of the P3HT oligomers, LUMO and HOMO energy levels calculated by density functional theory are found to be invariant across the donor/acceptor interface. Finally, these connections between additives, processing, and morphology at all length scales are generally useful for efforts to improve device performance.

36 MATERIALS SCIENCE↗

Coupled cluster Green function: Model involving single and double excitations

In this paper, we report on the development of a parallel implementation of the coupled-cluster (CC) Green function formulation (GFCC) employing single and double excitations in the cluster operator (GFCCSD). A key aspect of this work is the determination of the frequency dependent self-energy, Σ(ω). The detailed description of the underlying algorithm is provided, including approximations used that preserve the pole structure of the full GFCCSD method, thereby reducing the computational costs while maintaining an accurate character of methodology. Furthermore, for systems with strong local correlation, our formulation reveals a diagonally dominate block structure where as the non-local correlation increases, the block size increases proportionally. Here, to demonstrate the accuracy of our approach, several examples including calculations of ionization potentials for benchmark systems are presented and compared against experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sign Learning Kink-based (SiLK) Quantum Monte Carlo for molecular systems

The Sign Learning Kink (SiLK) based Quantum Monte Carlo (QMC) method is used to calculate the ab initio ground state energies for multiple geometries of the H{sub 2}O, N{sub 2}, and F{sub 2} molecules. The method is based on Feynman’s path integral formulation of quantum mechanics and has two stages. The first stage is called the learning stage and reduces the well-known QMC minus sign problem by optimizing the linear combinations of Slater determinants which are used in the second stage, a conventional QMC simulation. The method is tested using different vector spaces and compared to the results of other quantum chemical methods and to exact diagonalization. Our findings demonstrate that the SiLK method is accurate and reduces or eliminates the minus sign problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing microhydration effect on the electronic structure of the GFP chromophore anion: Photoelectron spectroscopy and theoretical investigations

The photophysics of the Green Fluorescent Protein (GFP) chromophore is critically dependent on its local structure and on its environment. Despite extensive experimental and computational studies, there remain many open questions regarding the key fundamental variables that govern this process. One outstanding problem is the role of autoionization as a possible relaxation pathway of the excited state under different environmental conditions. This issue is considered in our work through combined experimental and theoretical studies of microsolvated clusters of the deprotonated p-hydroxybenzylidene-2,3-dimethylimidazolinone anion (HBDI{sup −}), an analog of the GFP chromophore. Through selective generation of microsolvated structures of predetermined size and subsequent analysis of experimental photoelectron spectra by high level ab initio methods, we are able to precisely identify the structure of the system, establish the accuracy of theoretical data, and provide reliable description of auto-ionization process as a function of hydrogen-bonding environment. Our study clearly illustrates the first few water molecules progressively stabilize the excited state of the chromophore anion against the autodetached neutral state, which should be an important trait for crystallographic water molecules in GFPs that has not been fully explored to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Equation of motion coupled cluster methods for electron attachment and ionization potential in polyacenes

Polyacenes have attracted considerable attention due to their use in organic based optoelectronic materials. Polyacenes are polycyclic aromatic hydrocarbons composed of fused benzene rings. Key to understanding and design of new functional materials is an understanding of their excited state properties starting with their electron affinity (EA) and ionization potential (IP). We have developed a highly accurate and com- putationally e*fficient EA/IP equation of motion coupled cluster singles and doubles (EA/IP-EOMCCSD) method that is capable of treating large systems and large basis set. In this study we employ the EA/IP-EOMCCSD method to calculate the electron affinity and ionization potential of naphthalene, anthracene, tetracene, pentacene, hex- acene and heptacene. We have compared our results with other previous theoretical studies and experimental data. Our EA/IP results are in very good agreement with experiment and when compared with the other theoretical investigations our results represent the most accurate calculations as compared to experiment.

benzene↗

Toward enabling large-scale open-shell equation-of-motion coupled cluster calculations: triplet states of β-carotene

In this paper we discuss the application of novel parallel implementation of the coupled cluster (CC) and equation-of-motion coupled cluster methods (EOMCC) in calculations of excitation energies of triplet states in beta-carotene. Calculated excitation energies are compared with experimental data, where available. We also provide a detailed description of the new parallel algorithms for iterative CC and EOMCC models involving single and doubles excitations.

beta-carotene↗

Equation of motion coupled cluster methods for electron attachment and ionization potential in fullerenes C{sub 60} and C{sub 70}

In both molecular and periodic solid-state systems there is a need for the accurate determination of the ionization potential and the electron affinity for systems ranging from light harvesting polymers and photocatalytic compounds to semiconductors. The development of a Green's function approach based on the coupled cluster (CC) formalism would be a valuable tool for addressing many properties involving many-body interactions along with their associated correlation functions. As a first step in this direction, we have developed an accurate and parallel efficient approach based on the equation of motion-CC technique. To demonstrate the high degree of accuracy and numerical efficiency of our approach we calculate the ionization potential and electron affinity for C{sub 60} and C{sub 70}. Accurate predictions for these molecules are well beyond traditional molecular scale studies. We compare our results with experiments and both quantum Monte Carlo and GW calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bridging single and multireference coupled cluster theories with universal state selective formalism

The universal state selective (USS) multireference approach is used to construct new energy functionals which offers a unique possibility of bridging single and multireference coupled cluster theories (SR/MRCC). These functionals, which can be used to develop iterative and non-iterative approaches, utilize a special form of the trial wavefunctions, which assure additive separability (or size-consistency) of the USS energies in the non-interacting subsystem limit. When the USS formalism is combined with approximate SRCC theories, the resulting formalism can be viewed as a size-consistent version of the method of moments of coupled cluster equations (MMCC) employing a MRCC trial wavefunction. Special cases of the USS formulations, which utilize single reference state specific CC (V.V. Ivanov, D.I. Lyakh, L. Adamowicz, Phys. Chem. Chem. Phys. 11, 2355 (2009)) and tailored CC (T. Kinoshita, O. Hino, R.J. Bartlett, J. Chem. Phys. 123, 074106 (2005)) expansions are also discussed.

coupled cluster theories↗

Noniterative Multireference Coupled Cluster Methods on Heterogeneous CPU-GPU Systems

A novel parallel algorithm for non-iterative multireference coupled cluster (MRCC) theories, which merges recently introduced reference-level parallelism (RLP) [K. Bhaskaran-Nair, J.Brabec, E. Aprà, H.J.J. van Dam, J. Pittner, K. Kowalski, J. Chem. Phys. 137, 094112 (2012)] with the possibility of accelerating numerical calculations using graphics processing unit (GPU) is presented. We discuss the performance of this algorithm on the example of the MRCCSD(T) method (iterative singles and doubles and perturbative triples), where the corrections due to triples are added to the diagonal elements of the MRCCSD (iterative singles and doubles) effective Hamiltonian matrix. The performance of the combined RLP/GPU algorithm is illustrated on the example of the Brillouin-Wigner (BW) and Mukherjee (Mk) state-specific MRCCSD(T) formulations.

Brillouin-Wigner↗

Note: Excited State Studies of Ozone using State-Specific Multireference Coupled Cluster Methods

Vertical excitation energies obtained with state-specific multi-reference coupled cluster (MRCC) methods are reported for the ozone molecule. Using state-specific MRCC non-iterative methods with singles, doubles, and non-iterative triples (MRCCSD(T)) we obtain 4.40 eV for the challenging doubly excited 21A1 state when using a reliable model space. This estimate is in good agreement with experiment (4.5 eV). We also compare our MRCC results with the excitation energies obtained with high-order equation-of-motion coupled cluster methods

ab initio↗

Communication: Application of state-specific multireference coupled cluster methods to core-level excitations

The concept of the model space underlying multireference coupled-cluster (MRCC) formulations is a powerful tool to deal with complex correlation effects for various electronic states. Here, we demonstrate that iterative state-specific MRCC methods (SS-MRCC) based on properly defined model spaces can be used to describe core-level excited states even when canonical Hartree-Fock orbitals are utilized. We show that the SS-MRCC models with single and double excitations (SSMRCCSD) are comparable in accuracy to high-level single reference equation-of-motion coupled cluster (EOMCC) formalism. We also demonstrate that the SS-MRCC methods are capable of providing high accuracy results without experiencing numerical problems of the EOMCC methods.

coupled cluster↗

Implementation of the multireference Brillouin-Wigner and Mukherjee’s coupled cluster methods with non-iterative triple excitations utilizing reference-level parallelism

In this paper we discuss the performance of the non-iterative State-Specific Mul- tireference Coupled Cluster (SS-MRCC) methods accounting for the effect of triply excited cluster amplitudes. The corrections to the Brillouin-Wigner and Mukherjee MRCC models based on the manifold of singly and doubly excited cluster amplitudes (BW-MRCCSD and Mk-MRCCSD, respectively) are tested and compared with the exact full configuration interaction results (FCI) for small systems (H2O, N2, and Be3). For larger systems (naphthyne isomers and -carotene), the non-iterative BW-MRCCSD(T) and Mk-MRCCSD(T) methods are compared against the results obtained with the single reference coupled cluster methods. We also report on the parallel performance of the non-iterative implementations based on the use of pro- cessor groups.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Towards large-scale calculations with State-Specific Multireference Coupled Cluster methods: Studies on dodecane, naphthynes, and polycarbenes

Multireference Coupled Cluster (MRCC) methods are essential for accurate descriptions of molecular systems that have nearly degenerate states. A number of these methods have been developed over the last three decades but their use has been hampered by their high computational cost. The development of highly scalableMRCC codes can significantly alleviate this problem. Our recent work on the parallelization of the Brillouin-Wigner MRCC implementation has clearly demonstrated that taking advantage of large computational resources can enable the MRCC calculations on chemically relevant molecules within practical time-to-solution and memory per core constraints. In this letter we report on a parallel implementation of the Mukherjee’s state-specific MRCC model with single and double excitations (MR MkCCSD). In particular the abstraction of a methyl-group from dodecane, the singlet-triplet splitting in naphthyne isomers, and the singlet-high-spin splitting in polycarbenes are considered as demonstrations of systems with significant multireference character and of challenging size that currently can be handled with the MR MkCCSD formalism.

74 ATOMIC AND MOLECULAR PHYSICS↗